Lung Phantom
Author: Marc Kachelrieß

marc@imp.uni-erlangen.de

Figure 1: Lung phantom up to the 10th 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 pfeil.gif (855 Byte) B[b][s]
s pfeil.gif (855 Byte) 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