Learn more about Image Processing Toolbox

phantom - Create head phantom image

Syntax

P = phantom(def, n) P = phantom(E, n) [P, E] = phantom(...)

Description

P = phantom(def, n) generates an image of a head phantom that can be used to test the numerical accuracy of radon and iradon or other two-dimensional reconstruction algorithms. P is a grayscale intensity image that consists of one large ellipse (representing the brain) containing several smaller ellipses (representing features in the brain).

def is a string that specifies the type of head phantom to generate. Valid values are

'Shepp-Logan' — Test image used widely by researchers in tomography
'Modified Shepp-Logan' (default) — Variant of the Shepp-Logan phantom in which the contrast is improved for better visual perception

n is a scalar that specifies the number of rows and columns in P. If you omit the argument, n defaults to 256.

P = phantom(E, n) generates a user-defined phantom, where each row of the matrix E specifies an ellipse in the image. E has six columns, with each column containing a different parameter for the ellipses. This table describes the columns of the matrix.

Column	Parameter	Meaning
Column 1	`A`	Additive intensity value of the ellipse
Column 2	`a`	Length of the horizontal semiaxis of the ellipse
Column 3	`b`	Length of the vertical semiaxis of the ellipse
Column 4	`x0`	x-coordinate of the center of the ellipse
Column 5	`y0`	y-coordinate of the center of the ellipse
Column 6	`phi`	Angle (in degrees) between the horizontal semiaxis of the ellipse and the x-axis of the image

For purposes of generating the phantom, the domains for the x- and y-axes span [-1,1]. Columns 2 through 5 must be specified in terms of this range.

[P, E] = phantom(...) returns the matrix E used to generate the phantom.

Class Support

All inputs and all outputs must be of class double.

Remarks

For any given pixel in the output image, the pixel's value is equal to the sum of the additive intensity values of all ellipses that the pixel is a part of. If a pixel is not part of any ellipse, its value is 0.

The additive intensity value A for an ellipse can be positive or negative; if it is negative, the ellipse will be darker than the surrounding pixels. Note that, depending on the values of A, some pixels can have values outside the range [0,1].

Examples

P = phantom('Modified Shepp-Logan',200);
imshow(P)

Reference

[1] Jain, Anil K., Fundamentals of Digital Image Processing, Englewood Cliffs, NJ, Prentice Hall, 1989, p. 439.