Lung Phantom Author: Marc Kachelrieß |
Figure 1: Lung phantom up to the 10^{th} generation, i.e. 2047 branch objects and thus 8188 geometric objects. Bounding box: [-212mm, 205mm] ´ [-132mm, 132mm] ´ [-23mm, 297mm]. |
Aims
To provide semi–antropomorphic cone–beam rawdata for CT–angiography and –bronchoscopy.
Phantom Description
The Lung Phantom is a fractal model of the bronchial tree which has been described first in [1]. Each branch of the Lung Phantom consists of a hollow cylinder with a hollow half–sphere at each end. Out of one end of the branch grow two sprouts, a small and a big one. The plane in which lie the two sprouts is the so–called branch plane. The branch plane is rotated from generation to generation by a given angle to extent the structure from two to three dimensions and thus to homogeneously fill a given volume (the lung).
Phantom Definition
The phantom structure can mathematically be modelled as a L–System (Lindenmayer System). It consists of the starting axiom
T[L[s][b]][R[s][b]]
and the simultaneous substitutions
b
B[b][s]
s
S[b][s].
Substituting times yiels a bronchial
tree of up to the –th generation
since already generations 0, 1 and 2 are included in the axiom. The result
of the substitutions is the so–called product. Capitalized symbols are
the terminal symbols. They will not be substituted anymore. For example the
product for is
T[L[S[B][S]][B[B][S]]][R[S[B][S]][B[B][S]]]
where the remaining active symbols have been capitalized after the last substitutions to yield terminal symbols only. Each terminal symbol represents one of the following branch objects used to build the lung:
Trachea:
--------
BranchAngle: 0°
RotationAngle: 0°
Length: 100mm
OuterRadius: 9mm
InnerRadius: 7.9mm
Scalefactor: 0.86
RightMainBronchus:
------------------
BranchAngle: 55°
RotationAngle: 180°
Length: 50mm
OuterRadius: 7.5mm
InnerRadius: 6.4mm
Scalefactor: 0.86
LeftMainBronchus:
-----------------
BranchAngle: -55°
RotationAngle: 0°
Length: 40mm
OuterRadius: 6.5mm
InnerRadius: 5.4mm
Scalefactor: 0.86
BigSprout:
----------
BranchAngle: 25°
RotationAngle: 90°
Length: 50mm
OuterRadius: 8.5mm
InnerRadius: 7.4mm
Scalefactor: 0.86
SmallSprout:
------------
BranchAngle: -45°
RotationAngle: 10°
Length: 40mm
OuterRadius: 7mm
InnerRadius: 5.9mm
Scalefactor: 0.86
The BranchAngle determines the angle between the current branch and the new sprout. The RotationAngle describes the rotation of the branch plane after the new sprout has been added. Length obviously gives the length of the branch and OuterRadius as well as InnerRadius give the radii of the branch. The most important feature is the Scalefactor which determines how to scale the dimensions of the branch from generation to generation: If is the number of the generation considered () then its dimensions will be multiplied by Scalefactor.
The number of branch objects used to build the lung is given by . Nevertheless, as each branch object will be modeled by two cylinders and two spheres at the lower end (the spheres at the upper end can be left away since they completely lie within the mother branch of the last generation before) with radius OuterRadius and InnerRadius the complete lung will be made of basic objects.
Materials and Densities
Currently the phantom is designed to be monochromatic. The wall of each branch has the density of water, the airways have a density of 10% water.
Drawbacks
The program does not check for a self–penetration of the lung yet.
[1] Markus Blank. Untersuchungen zur Morphologie der Lunge anhand Spiral–CT–Daten. Dissertation, Friedrich–Alexander–Universität Erlangen–Nürnberg, May 1998.
View Phantom Definition (3-10 seperate generations - download as zip, 702k) |
phantom syntax explanation |
Other Phantoms | |||
Abdomen | Head | Hip | Jaw |
Lung | Thorax | High Contrast | Low Contrast |