Review of Tissue Simulating Phantoms For Optical Spectroscopy, Imaging and Dosimetry
Review of Tissue Simulating Phantoms For Optical Spectroscopy, Imaging and Dosimetry
1 Introduction for minimum performance criteria for new systems and for
1.1 Medical Tissue-Simulating Phantoms routine monitoring of existing systems. The benefit of this
procedure is that system performance can then be made more
The development of all diagnostic imaging systems and most uniform between institutions and over time.
physical therapeutic interventions has required the use of Access to these phantoms is made possible through com-
tissue-simulating objects to mimic the properties of human or mercial distributors who can manufacture them economically.
animal tissues. These so-called “phantoms” are used for a Unfortunately, in the development phase of imaging systems,
number of purposes,1–3 including: the status of tissue phantoms can be inconsistent and change
1. initially testing system designs over time, making comparison of research systems more dif-
2. optimizing signal to noise in existing systems ficult. In addition, considerable wasted effort occurs in rede-
3. performing routine quality control veloping tissue phantoms that have already been well de-
4. comparing performance between systems signed by previous groups. In the case of optical or near-
When systems are established and in routine clinical use with infrared imaging and spectroscopy, the field has developed
regulatory approval, there are generally requirements or rec- considerably over the past several decades, yet routine wide-
ommendations for quality control phantoms that need to be spread clinical use has not been established for many systems.
imaged for validation of system performance and use. Regu- In addition, the spectral range and geometrical range of optics
latory bodies such as the American College of Radiology applications are so diverse that development of systems and
共ACR兲, and medical physics associations such as the Ameri- tissue phantoms has not been a straightforward linear progres-
can Association of Physicists in Medicine 共AAPM兲 and/or sion. In this study, an overview of the various types of tissue
simulating phantoms and their applications is outlined. An
Canadian Organization of Medical Physicists 共COMP兲 make
attempt is made to discuss the strengths and weaknesses of
recommendations on requirements for phantoms to be used
each phantom type, and issues such as system purpose, geom-
etry, and tissue type are included. The tradeoffs between
Address all correspondence to Brian Pogue, Thayer School of Engineering, Dart- structure and biological or chemical function are also in-
mouth College, 8000 Cummings Hall, Dartmouth College Hanover, NH 03755
United States of America; Tel: 603 646-3861; Fax: 603 646-3856; E-mail:
pogu@[Link] 1083-3668/2006/11共4兲/041102/16/$22.00 © 2006 SPIE
cluded, in an effort to provide the most comprehensive listing introduced is there a discernable separation of the effects of
possible at this stage of development. a and s⬘. If the goal is to mimic the tissue transmission,
The history of tissue simulating phantoms for optical or then matching eff can often be sufficient, but in most tissue
near-infrared spectroscopy and imaging of tissue began in the spectroscopy applications where the goal is to separate a共兲
early 1980s with the surge of clinical interest in near-infrared and s⬘共兲 to allow spectral fitting, the tissue must have rep-
transillumination for breast cancer imaging, also termed resentative values for both these parameters. An excellent
diaphanography.4–7 Later interest also arose from applications compendium of tissue optical properties was compiled in the
in photodynamic therapy treatment planning and pulsed laser late 1980’s by Cheong, Prahl, and Welch,41 and updated in
treatment planning,8–14 where knowledge of the optical flu- 1995.42 Since that time, many more spectra have been pro-
ence distribution in tissue was critical to achieving treatment duced for dozens of different tissue types, including
efficacy. In the early 1990s, the introduction of spatially re- breast,43–47 brain,48 skin,49,50 esophagus, and cervix.51–54
solved, time-resolved, and frequency-domain light signals
spurred a larger number of researchers to investigate spectros- 1.3 Molecular/Flow/Structural Complexities of
copy and imaging of tissue, leading to the generation of many Optical Phantoms
different types of tissue phantoms.15–23 In recent years, the
applications of light in medicine have increased dramatically, Most of the early studies in tissue phantoms were focused on
with cosmetic laser surgery being a major commercial driving creating regular-shaped objects that mimicked tissue reduced
force, and fluorescence and reflectance diagnostics emerging scattering s⬘ and absorption a at specific wavelengths. In
as serious contenders for commercial success. Research into the past decade, focus has shifted to providing phantoms that
near-infrared tomography,24–29 photodynamic therapy reproduce tissue properties over broader wavelength ranges,
dosimetry,8,12,30 luminescence imaging,31–33 fluorescence mo- matching the full spectrum of tissue a共兲 and s⬘共兲 values.
lecular imaging,34–36 and optical coherence tomography37 There is also significant interest in developing biochemical
among other applications, keeps the area of tissue phantoms and biologically compatible tissue phantoms, which can uti-
progressing and important. Experimental progress toward mo- lize biologically important molecules such as hemoglobin,
lecular imaging applications requires tissue phantoms that melanin, or endogenous fluorophores such as nicotinamide
have some of the specific molecular features of human tissue. adenine dinucleotide 共NADH兲 and flavin adenine dinucleotide
At the same time, companies are developing tissue imaging 共FAD兲55 or exogenous fluorphores such as porphyrins or cya-
and spectroscopy devices that will require well-calibrated tis- nine dyes.56 Extending the biochemical capacity to measuring
sue phantoms for routine system comparison, evaluation, and transient biochemical species such as radicals or singlet oxy-
quality control. For all of these reasons, the improvement and gen has also been demonstrated, and provides actinometry
standardization of tissue optical phantoms is essential and capabilities for therapy planning and optimization.57,58
likely inevitable, even though this work has low priority in Generation of hybrid phantoms with specific characteris-
most research labs. tics for multimodality imaging, such as elastic properties,59
biochemical properties, water/lipid concentrations,60 electrical
properties,61 magnetic resonance properties, and thermal prop-
1.2 Tissue Optical Properties erties, together with optical properties, is becoming increas-
The key to matching tissue properties in phantoms is a com- ingly useful.62
prehensive understanding of the key physical and biochemical Along the lines of dynamically changeable phantoms,
characteristics of tissue that influence its interaction with there is also a need in some developments to study motion or
light.13 For small scale 共⬍1 mm兲 applications, it is likely im- mass displacement with optical signals. Several methods have
portant to match the absorption coefficient a共兲, the scatter- been developed to image motion in tissue, which ultimately
ing coefficient s⬘共兲, and the anisotropy coefficient g共兲, provides a good measure of mass flow, either by Doppler shift
which is defined as the average cosine of the scattering angle. measurements63,64 or correlation analysis of speckle.65–68
Over larger distances 共more than 3 to 5 scattering lengths, a In recent years, with advances in tissue engineering, a new
scattering length being defined as the reciprocal of the scat- emphasis has been placed on engineered tissue structures as
tering coefficient 1 / s兲 matching the reduced scattering coef- tissue-simulating phantoms for studies that investigate bio-
ficient s⬘ 关also called the transport scattering coefficient, de- logical chemistry or complex biochemical signatures.69 This
fined as s⬘ = 共1-g兲s兴 is all that is required.21 This “reduced” approach and the use of ex-vivo tissue70,71 have become estab-
approximation follows observations in neutral particle scatter- lished areas of investigation, although their use is distinctly
ing that over multiple scattering event lengths, an anisotropic different from the standard concept of a tissue phantom. The
scattering process appears identical to an isotropic scattering ability to better test systems in realistic situations with thin
process with a reduced value for the effective scattering tissue layers, anisotropic properties, and extracellular scaf-
coefficient.12,38–40 In many cases of thick tissue transmission, folding is essential in some applications. Each of these sub-
it is possible to get away with mimicking the effective attenu- jects is addressed in detail in this work.
ation coefficient of a tissue, defined in the wavelength regime
where diffusion theory is accurate as eff = 共3as⬘兲−1/2. This 1.4 Optical Tissue Phantom Composition Choices
is possible because steady-state attenuation in homogeneous In choosing the most useful phantom materials and design, the
media is affected in the same way by the same relative change region of the spectrum to be used is important, as are the
in absorption or scattering. Over long distances, diffusive pro- geometrical design parameters of thickness, heterogeneities,
cesses appear to be attenuated exponentially with this single container, and possible machining constraints. The biological
coefficient, and only when boundaries or temporal signals are compatibility in terms of biochemical action or inclusion of
Particle
Scatterer Biologically Organic Particle Index of distribution
material Permanent compatible chemical size 关nm兴 refraction function Recommended Use References
7,30,60,80,81,148,149
Lipids N Y Y 10 to 500 1.45 Exponentially Intralipid, milk, mixture
nm weighted to Theory/experimental
smaller sizes, tests and multiple phantom
impossible to contrast studies
get a single
size
distribution
Quartz glass Y Y Y 250 nm N/A Single size Used with resin Darmstadt Inc.,
microspheres function, with phantoms Germany152
10% variance
biologically relevant chromorphores and fluorophores is criti- 1. validation of physical models and simulations
cal as well. Since one of the important features of optical and 2. instrument performance testing and optimization
near-infrared 共NIR兲 spectroscopy is the spectral sensitivity to 3. instrument calibration and testing of stability and
molecular features of tissue, it has become increasingly im- reproducibility
portant to develop reliable phantoms that accurately mimic 4. interlaboratory comparison and standardization.
the chemistry of tissue. This requires a shift away from solid The properties of the ideal phantom depend on its intended
nonorganic polymers and silicone phantoms toward biologi- use. For example, validation phantoms need to be precisely
cally compatible structures such as agar, gelatin, or collagen characterized, but stability and reproducibility might not be as
matrixes that allow easy inclusion of cellular constituents important as in phantoms intended for interlaboratory com-
such as blood or fat and fluorescent molecules such as parisons. Thus, as the different phantoms are discussed, these
NADH, FAD, porphyrins,72,73 and other exogenous organic
four uses are kept in mind.
luminescent molecules.74–76.
An “ideal” phantom that could be used for any application
In this survey, the strengths of each approach are put
would have the properties listed as follows. As stated before,
alongside the ease of use, and in Tables 1–8 a summary of
these is included, along with recommendations for use for in real applications only some of these properties are impor-
each type of phantom. Because of the wide variety of phan- tant and the others can be neglected or given a lower priority.
toms and their constituents, it is not possible to have a single 1. Absorption and scattering properties can be varied as in
comprehensive table of constituents without having signifi- different tissues.
cant redundancy and overwhelmingly long tables. In an effort 2. Wavelength dependence of these properties is similar to
to streamline the presentations, the important parameters for tissue.
tissue optical phantoms are separated into scattering particles 3. Molecules of specific interest can be incorporated 共e.g.,
and matrix material. In the sections that follow, more detailed NADH, FAH, collagen, tetrapyrroles, fluorophores, and
discussion of each is provided to include all the pertinent actinometers兲.
details, and to reference the key studies that provide more 4. Properties are stable over time and environmental con-
complete directions of how to make and use these phantoms. ditions 共e.g., temperature, humidity, and photobleaching兲.
5. Index of refraction close to that of tissue 共e.g., index of
1.5 Purposes of Phantoms and the Criteria for tissue ⬇1.4兲.
Determining Their Value 6. Ability to incorporate regions with different optical
In general, the purposes of tissue optical phantoms can be properties 共e.g., inclusions mimicking tumors or layers mim-
roughly divided into the following categories: icking skin兲.
7. Mechanical and surface properties are similar to tissue Table 2 Summary of lipid emulsion-based phantoms.
共e.g., Young’s modulus near 4 to 20 kPa兲.77
8. Ability to incorporate Brownian motion or flow in the Scatterer Function Limitations Stability References
phantom. 6,149,153
9. Ability to include thermal properties similar to tissue. Milk Readily available Not highly Hours
reproducible
10. Ease of manufacturing.
between
11. Inexpensive to produce. samples
12. Easily transported between different sites.
7,60,106
Again, as the different compositions are analyzed, these fea- Oil/fat/lipid Used to custom- Must be Hours
tures are raised and discussed to compare each phantom ma- make scattering emulsified
terial with alternatives. and lipid/water and blended
phantoms reproducibly
80,81,148
2 Scattering Particles in Optical Phantoms Intralipid/ Reproducible Stability better Days
Nultralipid source of lipid than 10%
In most tissue phantoms, the choice of a scattering agent is solution
separate from the choice of matrix composition, as the volume
fraction of the scattering material is typically less than 5% of
the total, and often less than 1%. There have been three main
choices: lipid microparticles, polymer microparticles, and tency in the optical properties between batches. Clearly this is
white metal oxide powders and a brief list is shown in Table an emulsion suspension, and thorough mixing is required for
1. The benefit of lipid microparticles is that they are biologi- homogeneity. The homogeneity lasts for a period of hours,
cally similar to what is thought to cause scattering in tissue, while reuse of the solution has been reported over many days.
namely the bilipid membrane of cells and organelles. The next When used for ultraviolet studies, the lipid content in this
most common choice has been the polymer microsphere, with medium is likely to fluoresce which may interfere with trans-
polystyrene being the most popular. This is an excellent mission or remission studies, and so care must be taken to use
choice from a scientific perspective, because it is produced in nonorganic scattering materials such as are described in the
regular sizes with good quality control over the size and index next two sections.
of refraction. Thus, the repeatibilty and theoretical prediction
of the spectra are excellent. The third choice is common tita- 2.2 Scattering Coefficient Spectrum of Intralipid
nium dioxide or aluminum oxide powder. These are often the An excellent survey of the properties of Intralipid can be
main pigment in white paint and white plastics, due to their found at the website [Link]
high scattering coefficients, and they can be obtained in well- The most cited study by van Staveren et al.81 used measure-
controlled spherical formulations, although the use of these is ments of optical transmission as well as electron microscopy
less well established. Finally, in recent years scattering gold and Mie theory calculations to estimate the scattering spec-
nanoparticles have been developed, and their use in tissue trum. They also proposed a simple power law for the wave-
diagnostics and therapy has considerable promise due to their length dependence of the reduced scattering coefficient, which
high scattering cross section and potential has been utilized by many researchers. For a standard 10%
biocompatibility.78,79 While their use in phantoms is not well stock solution, the formulas for scattering and anisotropy co-
established, their significant Mie scatter cross section makes efficients are:
them a good potential scatterer. Each of these is discussed in
more detail later, and summarized in Table 1. s共兲 = 16−2.4 ,
2.1 Commercially Available Lipid-Based Scatterers units of mm−1, when is in microns, and
The most widely used phantoms for optical imaging and spec- g共兲 = 1.1 − 0.58,
troscopy have been the liquid type, made from milk149 or
emulsified oil suspensions initially, and later being largely re- resulting in the equation for reduced scattering coefficient of:
placed with the well calibrated, commercially available lipid
emulsion with the trade name Intralipid.148,152 These are listed s⬘共兲 = 9.3−1.4 − 1.6−2.4 共units of mm−1兲 .
in Table 2. For a more complete spectrum across the visible range, inclu-
Commerical supplies of calibrated lipid solutions are pos- sion of Rayleigh scattering is likely needed, requiring a third
sible due to their production for intravenous feeding.82 There term having the standard −4 power function, but requiring
is a number of commercial manufacturers, and the trade name fitting for the coefficient. Many analyses retain only the first
of the product varies between manufacturers. Intralipid™
term of this latter formula, and fit s⬘共兲 to the functional
共Kabi-Pharmacia, Erlangen Germany; Pharmacia and Upjohn,
form s⬘共兲 = a−b, with a and b as free parameters. When
Clayton, New Jersey; and Kabi-Vitrum Incorporated, Stock-
fitting is restricted to the near infrared, this can be a reason-
holm, Sweden兲 is the most commonly cited word, with other
able assumption over a limited wavelength range.
versions called Nutralipid™ 共Pharmacia, Quebec, Canada兲
and Liposyn II™ 共Abbott Labs Incorporated, Montreal兲. This
solution is readily available in all hospital pharmacy depart- 2.3 Polymer Microspheres
ments, and the uniformity between batches is thought to be From a scientific viewpoint, polystyrene microspheres pro-
excellent, although there is some contention about the consis- vide the best standard phantom, as they are well controlled in
Table 3 Phantom matrix options to hold the scatterers, absorbers, and fluorophores.
Absorber or
fluorophore Function Limitations Stability References
67,86,156,157
Whole blood Provide realistic Hours to days
tissue spectra and
oxygenation
function
92,97,110,152
Ink Provide nearly Not stable nor Days
flat absorption repeatable unless 共if remixed兲
spectra taken from a
calibrated sample
74,152
Molecular dyes Provide spectra Days to weeks
with wavelength
peaks
158
Fluorophores Compatible with May need to Days to weeks
aqueous avoid aggregation
dissolving effects with
compounds addition of
additional agents
159
Heterogeneities Test Clear enclosures Days
共scattering/absorption tomography and need to be
/fluorescent兲 imaging avoided due to
capabilities light channeling
Used to fill Index of
inclusions in refraction
solid phantoms changes may be
significant for
solid inclusions
2.4 Titanium Dioxide and Aluminum Oxide matrix materials are optimal for different applications, and the
Titanium dioxide 共TiO2兲 powder is perhaps the most common major types are summarized in Table 3. The use of these
choice for scatterering in science and engineering, and this phantom matrix compositions is discussed throughout the re-
stems from its wide availability as the main pigment in com- mainder of the work.
mon white paint. Aluminum oxide or barium oxide powders
are also excellent scatterers, and are commonly used for coat- 4 Aqueous Suspension Phantoms
ing the interior of integrating spheres where exceptionally Water-based phantoms can employ any of the three main scat-
high scatter and low absorption are required. TiO2 powder terers mentioned before lipid, microspheres, or TiO2 power in
comes in several forms and purities, including preformed mi- suspension. The absorption of such phantoms is due mainly to
crospheres, available from Dupont Chemical 共http:// water throughout most of the visible and near-infrared wave-
[Link]/tc/Titanium-Dioxide/兲. lengths. This absorption coefficient is sufficiently low below
The main disadvantage of TiO2 powder is that it resides in 700 nm that it can be ignored 共a ⬍ 0.002 mm−1兲, and ab-
suspension in most media, and so settles when not stirred. sorbers can be added to tailor the absorption coefficient and
This is not a problem in resin or agar phantoms once they are spectrum to that of tissue. The water absorption spectrum can
set, but is an issue for aqueous phantoms. Continuous stirring be reliably assumed to match the measurements of Hale and
of the aqueous suspension produces a homogeneous phantom. Querry,87 and an excellent overview of the water spectra
For resin- or agar-based phantoms, mixing for extended peri- available and their conversion between different units is found
ods is also important to ensure that the particles are uniformly at the website [Link] A brief summary of ab-
distributed. Automated stirring for more than 30 min has been sorbers and fluorophores used is listed in Table 4.
a reliable approach for manufacture of reproducible resin
phantoms. Liquid-based stock supplies of TiO2 are now avail- 4.1 Exogenous Absorbers in Aqueous Phantoms
able from Sigma, and these may be a more reliable additive,
as the scattering properties are better controlled than a powder Addition of absorbers and fluorophores to aqueous phantoms
mixed in suspension. has been demonstrated in hundreds of studies, and this phan-
tom design has proven to be extremely valuable in the initial
validation of an imaging/spectroscopy system. Typically the
3 Bulk Matrix Materials for Optical Phantoms goal has been to mimic tissue, so the addition of
The choice of the bulk material for the phantom has perhaps erythrocytes,58,156–164 whole blood, or hemoglobin have all
the largest impact on how the phantom can be used. Different been reported. When attempting to preserve the oxygen bind-
ing function of hemoglobin, the use of saline rather than dis- used in magnetic-resonance imaging 共MRI兲 and ultrasound
tilled water is important, otherwise the blood cells will lyse imaging for decades,93–95 and they were adopted in optical
and the heme will dissociate from the hemologlobin molecule. tissue phantoms in many laboratories in the mid 1990’s.58,96–99
However, for simplicity, many researchers have chosen to use Agar and gelatin allow inclusion of organic molecules and
ink or molecular absorbing dyes to simulate the absorption of cellular-based constituents, while providing a semisolid object
blood or melanin in tissue. that can have a variety of shapes. Gelatin and agar phantoms
Addition of fluorophores158 has been reported in many have had an equally rich period of development in ultrasound
studies, with hydrophilic molecules used most successfully. imaging, and a large number of papers describe the diversity
Aggregation of certain hydrophobic dyes such as protopor- of phantoms developed here.62,100,101 More recently, the whole
phyrin 9 is possible, but addition of 5% Tween-29 共Fisher area of hydrogels has been studied for biocompatibility and
Scientific, USA兲 as an emulsifying agent has been found to drug delivery applications, and this encompasses most bio-
correct this and result in a monomerized form of the fluoro- logical scaffolds that alter the behavior of water.
phore. The absorption and fluorescence spectra are similar to Polyacrylamide hydrogels have been used as scaffolds for
those observed when the dye is dissolved in a dilute organic collagen and other matrices,102,103 and polyvinyl alcohol hy-
solvent. drogels are reviewed in Sec. 8.1, as they have intrinsic scat-
tering properties as well as being a matrix medium. Polyacry-
4.2 Inclusions and Heterogeneous Lipid Solution lamide hydrogel use has undergone enormous development in
Phantoms biological laboratories for use in electrophoresis and molecu-
An important complication in the use of lipid solution phan- lar separation techniques, yet there are only a few reports of
toms is the choice of container and the possibility of light testing of these matrices as phantom materials.102,103
channeling through the container walls, rather than through
the solution. This is especially problematic over longer dis-
tances and in cases when inclusions or heterogeneities are to 5.1 Scattering Composition
be incorporated in phantoms. Early studies in tomography Since gelatin phantoms are usually used for periods of a day
used containers with thin mylar walls to hold liquid to a week or more and then discarded, they are commonly
inclusions,88,89 but it is apparent that the mylar itself does made with less expensive scattering particles. Construction
perturb the light field. Correction for this effect can be per- with polystyrene microspheres is possible but is quite expen-
formed by filling the inclusion with the same solution as the sive. Use of titanium dioxide 共TiO2兲 power or aluminum ox-
background medium and using this as the “homogeneous” ide 共Al2O3兲 powder is the norm, as they are inexpensive and
reference phantom. However, for smaller inclusions and lower provide a reasonably reliable means to mix a scatterer into the
contrast inclusions, this approach not accurate, and solid liquid gelatin or agar solution while it is cooling. The major
phantoms are preferable. complicating factor in production of these phantoms is the
Light channeling along the top surface of Intralipid over need for careful attention to detail and procedure. For ex-
long distances has also been noted, yet little discussed in pub- ample, the TiO2 scatterers in the phantom readily precipitate
lications. It is important when using this or any optical phan- out, and when ordered in bulk comes in clumpy power form,
tom to shield the surfaces of the phantom properly so that requiring continuous stirring for approximately 20– 30 min to
signals can enter and exit the phantom only at desired loca- ensure homogeneous dispersion in the phantom. Liquid-based
tions. This can be achieved with black masking of the phan- TiO2 is also available and is a reliable method to add con-
tom surfaces, using any opaque acrylic or plastic material. trolled amounts to a solution without the concerns of being
Lipid-based solutions have been used with great success in able to mix and declump the suspension. In addition, TiO2
conjunction with solid phantoms, where holes or channels does settle over time, so the final scattering coefficient of
have been left in the solid phantom to allow dynamic varia- phantoms can vary significantly from one to another. Despite
tion of the heterogeneity optical properties.90–92 This approach careful procedures with TiO2 suspensions, repeated studies in
has been used in many studies to assess detectability of ob- our laboratory show that up to 50% variation can occur. Mak-
jects of differing contrast. This approach allows the use of ing multiple phantoms from a large batch of agar can reduce
contrast-detail analysis as well, to determine the minimum this intersample variation. In addition, the bottom of each
contrast detectable for each size of inclusion in the phantom.90 phantom typically has a large precipitation of TiO2, indicative
There is clearly concern that the transition from solid to liquid of a scattering gradient along the vertical direction. The pre-
matrix involves a change of refractive index, yet experiments cipitated TiO2 is thought to be from larger particles of the
appear to indicate that this is a manageable, if not insignifi- powder, which have higher gravitational force acting on them.
cant, artifact. Increased mixing time reduces the number of these particles.
However, in the end, it is imperative to be able to indepen-
5 Hydrogel-Based Phantoms dently measure the scattering coefficient prior to use in these
types of phantoms, due to inability to exactly predict the scat-
Most substances that encapsulate water as a main component
tering coefficient from a set recipe.
and form a stiff matrix that has limited water mobility are in
the category of hydrogels. Gelatin and agarose are two of the
most common examples, and in biological laboratories there
are hundreds of varieties of these. In this section, agar and 5.2 Additives to Gelatin Phantoms to Improve
gelatin are discussed separately because of their long history Function
as phantom matrix materials. Agar-based phantoms have been Table 5 summarizes some main additives used to improve
the function of gelatin based phantoms. Inclusion of 0.2% agar-based phantoms, but these can become fragile and
formaldehyde in gelatin phantoms increases the melting tem- crumble under applied stress. Gelatin can be ordered from
perature of the gelatin matrix by increasing the crosslinking of different biological origins, and with different bloom levels—
the fibers while preserving the lower Young’s modulus.101,104 increasing the level of bloom results in a stiffer gelatin phan-
This allows the phantom to be used at room temperature with- tom. A pig-skin-based gelatin with a bloom of 175 provides a
out need for refrigeration. This can also be achieved with good stiffness for reliable phantoms.104
Inclusion of biochemically toxic species such as wood pre- scatterers.152,155 The construction of these phantoms requires
servative 共0.01 g / L兲,61 mild acid ethylenediamine tetraacetic mixing a resin and hardener to create a transparent solid resin,
acid 共EDTA兲 at 0.02 g / L,92,104,105 or sodium azide98 provides which typically sets within a few days at room temperature or
a stable phantom that lasts for many days and weeks without within a few hours at elevated temperature. A detailed outline
bacterial growth. The EDTA additive is probably most com- of this procedure can be found at the University College Lon-
mon, because its lower toxicity simplifies handling proce- don website [Link]
dures. Inclusion of penicillin has also been reported for the research/NIRគtopics/phantomគ[Link]. An alternative
same reason.98 While these additives will maintain good bio- recipe can be found at [Link]
logical stability for many days and weeks, they will not keep [Link]. This material can be obtained from a number of
the gelatin from drying out, and the phantoms must be kept manufacturers and in different compositions including: 1.
sealed in airtight enclosures such as plastic bags or containers. Araldite epoxy 共MY753兲 and hardener 共XD716兲, Aeropia
Keeping the phantoms in vegetable oil has also been reported Chemical Supplies, 共Crawley, United Kingdom兲, or 2.
as an excellent way to preserve the water content.104 This Araldite resin 共GY502兲 and hardener 共HY832兲, D. H. Litter
process can provide an intact matrix for years of use of a Incorporated, 共Elmsford, New York兲. Thorough mixing of the
single gelatin phantom, although the other biochemical mol- resin and hardener is critical to obtain a homogeneous volume
ecules included may not last as long as the gelatin matrix that cures in a timely manner. There is significant heat and gas
itself. generated during this process. Degassing of the phantom dur-
Blood has been added to gelatin phantoms55,98 and pro- ing the initial curing process is critical to avoid large numbers
vides an excellent model of tissue spectra in the near infrared, of air bubbles embedded in the phantom. Initial degassing
where the dominant absorbers are hemoglobin and water. In-
during the curing process will cause a massive expansion of
clusion of fat has been reported, but without extensive study
the resin due to the large amount of gas present; however,
of this capability.106
delaying the onset of the degassing, or repeated rapid degas-
For therapeutic study use, these phantoms are ideal, as they
sing and repressurizing cycles, can break the bubbles present
can have the same elastic properties as human tissue and simi-
in the phantom and gradually reduce the phantom volume to
lar thermal properties.107 Inclusion of actinometry agents has
be predominantly resin with little gas.
been demonstrated and used to compare photodynamic dose
A recipe used in our laboratory for a reliable phantom is as
deposition from cw and pulsed laser sources.57,58 Similarly,
follows,90,92 mixing 100 parts Araldite GY502 mixed with 30
measurement of oxygen in phantoms and tissues can be
achieved with fluorescent reporters.108 The potential for bio- parts of the hardener HY 837. Prior to mixing in the hardener,
chemical similarity to tissue, together with the potential for the scatterer and absorber can be mixed into the Araldite thor-
therapy studies, makes these tissue phantoms the best for oughly and degassed to allow a homogeneous mixture of op-
complex tissue geometries and biophysical study. tical properties prior to initiation of the hardening process.
Then when the hardener is added, it can be slowly mixed to
minimize inclusion of air bubbles during the mixing process.
6 Polyester and Polyurethane Resin Phantoms It has been found that 3.5 g of TiO2 powder per liter of resin
Polyester resin phantoms were introduced by Firbank, Delpy, provides a scattering coefficient near 1.0 mm−1 at 800 nm,
and Oda using both TiO223 and polystyrene particle which is proportional to this concentration. Mixing for an
extended period of time with a magnetic stir bar or an electric Table 7 Additives in RTV silicone phantoms.
mixer is strongly recommended.
In previous studies, the bulk absorption coefficient of the Additive Function Limitations Stability References
medium was set by adding 25⫻ 10−6 liters of ink per liter of 115,116,92,117,164
resin, which was found to increase the absorption coefficient TiO2 Scattering Mixing Years
Al2O3 consistency
to a range between 0.006 to 0.009 mm−1, but different ink
is critical
bottles and solutions will vary significantly, so well calibrated Degassing
and mixed samples of ink stock solution must be used. Par- to remove
ticulate ink absorbers, such as India ink, produce a relatively air bubbles
flat absorption spectrum across most of the visible and near is critical
infrared, as they are composed of carbon particles suspended 150
Molecular Possible Years
in an emulsion. It is important to note that particulate inks absorbers
also scatter light,110 so quantification of the absorption coeffi-
cient of ink in standard spectrophotometers is not possible. 118
Heterogeneities Test Years
Instead, it must be measured in a standard “added absorber” tomography
experiment. Use of molecular absorber inks such as 900NP and imaging
has been firmly established.111–113 Many types of these nonor- capabilities
ganic dyes have been successfully added to this matrix and
provide wavelength-dependent absorption across the near in-
frared, and have no significant scattering coefficient as they this material into biologically relevant configurations with the
are smaller molecules. With consistent procedures, it is pos- stiffness of human tissue has been the main argument for its
sible to obtain a process where phantoms produced succes- use. This was demonstrated by Bays et al. for esophageal
sively have absorption and scattering properties within 10% phantoms115 intended for dosimetry for photodynamic treat-
of their target value. A summary of additives typically used is ment planning. Jiang et al. used this material for breast
in Table 6. phantoms92 to help in calibration of an optical tomography
Polyurethane phantoms were more recently described by system. Lualdi et al. have used these phantoms to study im-
Vernon et al.114 and suggested as a superior alternative to aging of skin lesions using melanin and absorbers that mimic
polyester resin, due to their better compatibility with infrared skin lesions.118 The only major drawbacks of this matrix ma-
dyes to better match the absorbing and fluorescent molecules terial are cost and hardening time, but these are not prohibi-
of tissue. It is stated that these resins provide less bleaching of tive, and a pliable tissue phantom can be quite useful for
the dyes, but extensive testing has not been reported. The applications where the mechanical contact to tissue is impor-
transparency and index change are similar to polyester, mak- tant.
ing these phantoms otherwise quite similar.
8 Novel Materials for Optical Phantoms with
7 Room-Temperature-Vulcanizing Silicone Intrinsic Scattering
Phantoms In addition to the materials discussed in the previous two
Room-temperature-volcanizing 共RTV兲 silicone-based soft sections, there are a number of materials that have been used
phantoms were introduced by Bays et al.150 and Beck et al.116 for phantoms that have intrinsic matrix and scattering proper-
The merits of these phantoms are that they are quickly pro- ties that are interlinked. These are less clearly organized than
duced, have a soft rubber texture similar to stiff tissue, and the previous group, but have properties that could make them
can include nonorganic scatterers and absorbers. The RTV- useful options for certain studies. These range from polyvinyl
based compounds can be obtained from a number of manu- alcohol gels, dough, and teflon, to “engineered” or excised
facturers 共RTV Elastosil 604, Wacker, Munich Germany116,117兲 tissues. Each of these are briefly mentioned in Table 8, and
共Rhodorsil RTV 141, Rhone-Poulenc, France115兲, 共RTV-141, summarized in the following subsections.
Medford Silicone, Medford, New Jersey59兲. Preparation of the
material is similar to the resin-based phantoms described in 8.1 Polyvinyl Alcohol Phantoms
the previous section. Mixing the RTV with its hardener ini- Perhaps the most promising and widely used of these options
tiates a chemical process that solidifies the compound, and are the polyvinyl alcohol gels,119,120 sometimes referred to as
heat and gas generation require pumping under vacuum. This cryogels, due to the fact that their scattering coefficient and
degassing removes the bubbles that are generated when it is stiffness increase with repeated freeze/thaw cycles, allowing
curing. them to be tailored for specific applications. These were origi-
A summary of some additives used with RTV phantoms is nally used in ultrasound and MRI121–123 research, and have
used in Table 7. Beck et al.116 examined ways to embed ab- recently been adopted for photoacoustic tomography, where
sorbers and scatterers into the medium, with the conclusion the combination of elastic and optical scattering properties
that certain stable dyes could be added, but organic molecules makes them ideal for this hybrid imaging approach.124
such as porphyrins were not stable in this polymer. Jiang et al. Kharine et al. report reduced scattering coefficients near
examined controlling stiffness by lowering the hardener con- 0.8 mm−1 after seven freeze/thaw cycles. They also demon-
centration. They showed that the elastic modulus of the phan- strated the ability to create pliable phantoms this way, without
tom could be lowered by a factor of 3, 共from 230 to 80 kPa兲, increased scattering, by including dimethyl sufoxide 共DMSO兲,
making it closer to the stiffness of soft human tissues. Shaping thereby reducing the apparent “whiteness” produced by water
Table 8 Phantom materials and tissues with intrinsic scattering within the matrix material.
Solid/ Organic
Phantom liquid/ Biologically chemical Inclusions Adjustable Adjustable Index of Recommended
material Permanent flexible compatible compatible possible? absorption scattering refraction use References
119,120,124
Polyvinyl Y F N N Y Y Y 1.36 Heterogeneity
alcohol studies with
gel deformation
Thermal studies
Acoustic studies
共higher cost OK兲
freezing in the cycles. This phantom can then be considered a used in tomography phantom studies 共unpublished data兲. The
clear matrix, in which microspheres or TiO2 could be embed- absorption coefficient appears to track linearly with the addi-
ded to create well-controlled optical scattering phantoms tion of higher and higher absorber concentration, as would be
while the elastic properties are set independently. These phan- expected.
toms appear highly promising for use in these hybrid applica-
tions where optical and stiffness properties need to be sepa-
rately controlled. 8.3 Engineered Tissues as Phantoms
These gels can be obtained with average molecular weight
Tissue engineering tools have evolved in the past decade to
of 85 to 140 kDa from Sigma-Aldrich 共USA兲 共catalog num-
the point where structures can be created or grown in culture
ber 36 314-6兲, and are dissolved at a concentration of 20% by
that mimic the structural properties of tissues.108,125–132 These
weight in distilled water while being heated to 90 ° C for 2 h
tissues are most important in situations where the subtle com-
with continuous stirring. After cooling for a few hours to al-
plexities of the biochemistry or thin-layered structure of the
low air bubbles to migrate to the surface, it is then poured into
tissue are simply not well understood, and therefore cannot be
a mold and frozen at −20 ° C for 12 h. This gel is then thawed
fully reproduced by inert tissue phantoms. This issue is espe-
at room temperature and refrozen to produce a stiffer gelatin
cially important in optics for anisotropic scattering due to
matrix, and this can be repeated several times to produce
structures such as collagen matrix133 or muscle fibers, or
stiffer and stiffer phantoms. Without the addition of DMSO,
where the layered sequence of tissues affects the light trans-
the scattering coefficient will increase with each cycle as
port into or out of the tissue.
well.124 This matrix is sensitive to humidity and will likely
Study of epithelial squamous tissues has been a primary
require storage and preparation under humidity controlled
area for this approach,126 mainly due to the possibility of
conditions.
growing epithelial cells on a thin collagen matrix, with me-
dium flowing above and below the culture. This is called a
8.2 Dough-Based Phantoms “raft” culture system, because the cells float on a raft of col-
While the concept of using dough or Play-DohTM may appear lagen. The layered structures of squamous epithelium can be
unscientific, these phantoms have considerable promise be- spontaneously developed, allowing in-vitro study of cellular
cause of their ease of construction, ease of use, and long stor- growth, differentiation, and expression of proteins, so this sys-
age time. Composition of these phantoms is based on a recipe tem is “organotypic” in structure and function. Spectroscopy
for the children’s toy, playdough. The standard mixture can be of these layered structures reveals a biochemical spectrum in
obtained from hundreds of websites, but one such recipe is which the influence of the layered features of the tissue is
250-ml flour, 125-ml salt, 15-ml vegetable oil, 30-ml cream unique and not well modeled by a simple phantom.127,128
of tartar, and 250-ml water. After mixing flour, salt, and oil, While this field is arguably still in its infancy, the potential
slowly add the water. Heat slowly and stir until dough be- is reasonably good for these models to become main stream
comes stiff. When a homogeneous dough ball forms, the mix- tools in molecular imaging studies. As engineered tissues be-
ture is then cooled and left to set. Various absorbers can be come more reproducible between laboratories, this becomes a
easily mixed into the dough as well, with India ink being used viable option. Another rationale for the use of these structures
successfully. This composition leads to a pliable phantom is as a replacement for animal studies. Alternatives to animal
with s⬘ = 1.6 mm−1 at 800-nm wavelength. Repeated mix- models are usually welcome in laboratories as long as the
tures had similar scattering coefficients, and were successfully model is a true representation.127,133
9 Conclusions Acknowledgments
This summary of phantoms and phantom materials is an at- This study was supported through grants PO1CA84203,
tempt to identify common themes in a field that has a large PO1CA80139, RO1CA109558, PO1CA43892, and a grant
diversity of applications and methods distributed throughout from the National Cancer Institute of Canada. The authors
hundreds of research laboratories. Major problems exist in wish to gratefully acknowledge assistance and collaboration
tissue phantom work due to the lack of uniformity and the in phantom making and analysis over many years with col-
lack of a “gold standard” for comparison. However, the leagues at the Juravinski Cancer Center 共Hamilton, Ontario,
strengths and weaknesses of phantom technology are best dis- Canada兲, Princess Margaret Hospital 共Toronto, Ontario,
cussed in terms of the application, as mentioned at the begin- Canada兲, Wellman Center for Photomedicine at the Massachu-
ning of this work. setts General Hospital 共Boston, Massachusetts兲, and Thayer
For application of phantoms in validating theoretical or School of Engineering at Dartmouth College 共Hanover, New
experimental systems, optimal choices are based on well cali- Hampshire兲.
brated and known quantities, and so microspheres or In-
tralipid are excellent choices, and allow the use of aqueous
emulsions or gelatin-based solid phantoms. This approach has References
become the de facto standard, although there is really no uni- 1. G. Cohen, “Contrast—detail—dose analysis of six different com-
versally accepted method to measure phantom optical proper- puted tomographic scanners,” J. Comput. Assist. Tomogr. 3共2兲, 197–
203 共1979兲.
ties. Integrating sphere setups for reflectance and transmit- 2. S. E. Seltzer, R. G. Swensson, P. F. Judy, and R. D. Nawfel, “Size
tance measurement are widely considered the best way to discrimination in computed tomographic images. Effects of feature
assess phantom properties, but this approach is still prone to contrast and display window,” Invest. Radiol. 23共6兲, 455–462 共1988兲.
3. J. B. Olsen and E. M. Sager, “Subjective evaluation of image quality 28. B. W. Pogue, J. D. Pitts, M. A. Mycek, R. D. Sloboda, C. M. Wilmot,
based on images obtained with a breast tissue pattern: comparison J. F. Brandsema, and J. A. O’Hara, “In vivo NADH fluorescence
with a conventional image quality phantom,” Br. J. Cancer 68共806兲, monitoring as an assay for cellular damage in photodynamic
160–164 共1995兲. therapy,” Photochem. Photobiol. 74共6兲, 817–824 共2001兲.
4. D. J. Watmough, “Diaphanography: mechanism responsible for the 29. B. Chance, S. Nioka, J. Zhang, E. F. Conant, E. Hwang, S. Briest, S.
images,” Acta Radiol. Oncol. 21共1兲, 11–15 共1982兲. G. Orel, M. Schnall, and B. J. Czerniecki, “Breast cancer detection
5. D. J. Watmough, “Transillumination of breast tissues: factors govern- based on incremental biochemical and physiological properties of
ing optimal imaging of lesions,” Radiology 147共1兲, 89–92 共1983兲. breast cancers: A six-year, two-site study,” Acad. Radiol. 12共8兲, 925–
6. B. Drexler, J. L. Davis, and G. Schofield, “Diaphanography in the 933 共2005兲.
diagnosis of breast cancer,” Radiology 157, 41–44 共1985兲. 30. B. C. Wilson, P. J. Muller, and J. C. Yanch, “Instrumentation and
7. J. Linford, S. Shalev, J. Bews, R. Brown, and H. Schipper, “Devel- light dosimetry for intra-operative photodynamic therapy 共PDT兲 of
opment of a tissue-equivalent phantom for diaphanography,” Med. malignant brain tumours,” Phys. Med. Biol. 31共2兲, 125–133 共1986兲.
Phys. 13共6兲, 869–875 共1986兲. 31. S. Walenta, T. Schroeder, and W. Mueller-Klieser, “Metabolic map-
8. L. I. Grossweiner, “Optical dosimetry in photodynamic therapy,” La- ping with bioluminescence: basic and clinical relevance,” Biomed.
sers Surg. Med. 6共5兲, 462–465 共1986兲. Eng. (NY) 18共6兲, 249–262 共2002兲.
9. S. L. Jacques and S. A. Prahl, “Modeling optical and thermal distri- 32. C. H. Contag and M. H. Bachmann, “Advances in in vivo biolumi-
butions in tissue during laser irradiation,” Lasers Surg. Med. 6共6兲, nescence imaging of gene expression,” Annu. Rev. Biomed. Eng. 4,
494–503 共1987兲. 235–260 共2002兲.
10. L. I. Grossweiner, J. H. Hill, and R. V. Lobraico, “Photodynamic 33. C. H. Contag and B. D. Ross, “It’s not just about anatomy: In vivo
therapy of head and neck squamous cell carcinoma: optical dosimetry bioluminescence imaging as an eyepiece into biology,” J. Magn. Re-
and clinical trial,” Photochem. Photobiol. 46共5兲, 911–917 共1987兲. son Imaging 16, 378–387 共2002兲.
11. S. T. Flock, B. C. Wilson, and M. S. Patterson, “Total attenuation 34. E. M. Sevick-Muraca, J. P. Houston, and M. Gurfinkel,
coefficients and scattering phase functions of tissues and phantom “Fluorescence-enhanced, near infrared diagnostic imaging with con-
materials at 633-nm,” Med. Phys. 14共5兲, 835–841 共1987兲. trast agents,” Curr. Opin. Chem. Biol. 6共5兲, 642–650 共2002兲.
12. W. M. Star, J. P. A. Marijnissen, and M. J. C. van Gemert, “Light 35. V. Ntziachristos, C. H. Tung, C. Bremer, and R. Weissleder, “Fluo-
dosimetry in optical phantoms and in tissues: I. Multiple flux and rescence molecular tomography resolves protease activity in vivo,”
transport theory,” Phys. Med. Biol. 33共4兲, 437–454 共1988兲. Nat. Med. 8共7兲, 757–760 共2002兲.
13. S. L. Jacques, “Laser tissue interactions—photochemical, photother- 36. V. Ntziachristos, C. Bremer, E. E. Graves, J. Ripoll, and R. Weissle-
mal, and photomechanical,” Surg. Clin. North Am. 72共3兲, 531–558 der, “In vivo tomographic imaging of near-infrared fluorescent
共1992兲. probes,” Mol. Imaging 1共2兲, 82–88 共2002兲.
14. W. M. Star, “Light dosimetry in vivo,” Phys. Med. Biol. 42共5兲, 763– 37. J. G. Fujimoto, C. Pitris, S. A. Boppart, and M. E. Brezinski, “Optical
788 共1997兲. coherence tomography: an emerging technology for biomedical im-
15. D. T. Delpy, M. Cope, P. Van der Zee, S. R. Arridge, S. Wray, and J. aging and optical biopsy,” Neoplasia 2共1–2兲, 9–25 共2000兲.
S. Wyatt, “Estimation of optical pathlength though tissue from direct 38. D. R. Wyman, M. S. Patterson, and B. C. Wilson, “Similarity rela-
time of flight measurement,” Phys. Med. Biol. 33, 1433–1442 共1988兲. tions for the interaction parameters in radiation transport,” Appl. Opt.
16. M. S. Patterson, B. Chance, and B. C. Wilson, “Time resolved reflec- 28共24兲, 5243–5249 共1989兲.
tance and transmittance for the non-invasive measurement of tissue 39. M. S. Patterson, B. C. Wilson, and D. R. Wyman, “The propagation
optical properties,” Appl. Opt. 28, 2331–2336 共1989兲. of optical radiation in tissue I. models of radiation transport and their
17. S. L. Jacques, “Time-resolved reflectance spectroscopy in turbid tis- application,” Lasers Med. Sci. 6, 155–168 共1990兲.
sues,” IEEE Trans. Biomed. Eng. 36共12兲, 1155–1161 共1989兲. 40. S. R. Arridge, M. Cope, and D. T. Delpy, “The theoretical basis for
18. S. L. Jacques, “Time resolved propagation of ultrashort laser-pulses the determination of optical pathlengths in tissue: temporal and fre-
within turbid tissues,” Appl. Opt. 28共12兲, 2223–2229 共1989兲. quency analysis,” Phys. Med. Biol. 37共7兲, 1531–1560 共1992兲.
19. B. C. Wilson and S. L. Jacques, “Optical reflectance and transmit- 41. W. F. Cheong, S. A. Prahl, and A. J. Welch, “A review of the optical
tance of tissues—principles and applications,” IEEE J. Quantum properties of biological tissues,” IEEE J. Quantum Electron. 26共12兲,
Electron. 26共12兲, 2186–2199 共1990兲. 2166–2185 共1990兲.
20. B. C. Wilson, T. J. Farrell, and M. S. Patterson, “An optical fiber- 42. W. F. Cheong, “Summary of optical properties,” in Optical-Thermal
based diffuse reflectance spectrometer for non-invasive investigation Response of Laser-Irradiated Tissue, A. J. Welch, Ed., Append. Chap.
of photodynamic sensitizers in vivo,” in Future Directions and Ap- 8, Plenum Press, New York 共1995兲.
plications in PDT, SPIE Inst. Series, pp. 219–232, SPIE Press, Bell- 43. K. Suzuki, Y. Yamashita, K. Ohta, M. Kaneko, M. Yoshida, and B.
ingham, WA 共1990兲. Chance, “Quantitative measurement of optical parameters in normal
21. T. J. Farrell, M. S. Patterson, and B. C. Wilson, “A diffusion theory breasts using time-resolved spectroscopy: In vivo results of 30 Japa-
model of spatially resolved, steady-state diffuse reflectance for the nese women,” J. Biomed. Opt. 1共3兲, 330–334 共1996兲.
noninvasive determination of tissue optical properties,” Med. Phys. 44. R. Cubeddu, C. D’Andrea, A. Pifferi, P. Taroni, A. Torricelli, and G.
19共4兲, 879–888 共1992兲. Valentini, “Effects of the menstrual cycle on the red and near-infrared
22. S. J. Madsen, B. C. Wilson, M. S. Patterson, Y. D. Park, S. L. optical properties of the human breast,” Photochem. Photobiol. 72共3兲,
Jacques, and Y. Hefetz, “Experimental tests of a simple diffusion- 383–391 共2000兲.
model for the estimation of scattering and absorption-coefficients of 45. N. Shah, A. Cerussi, C. Eker, J. Espinoza, J. Butler, J. Fishkin, R.
turbid media from time-resolved diffuse reflectance measurements,” Hornung, and B. Tromberg, “Noninvasive functional optical spectros-
Appl. Opt. 31共18兲, 3509–3517 共1992兲. copy of human breast tissue,” Proc. Natl. Acad. Sci. U.S.A. 98共8兲,
23. M. Firbank and D. T. Delpy, “A design for a stable and reproducible 4420–4425 共2001兲.
phantom for use in near-infrared imaging and spectroscopy,” Phys. 46. S. Srinivasan, B. W. Pogue, S. Jiang, H. Dehghani, C. Kogel, S.
Med. Biol. 38, 847–853 共1993兲. Soho, J. G. Chambers, T. D. Tosteson, S. P. Poplack, and K. D.
24. E. Gratton, S. Fantini, M. A. Franceschini, G. Gratton, and M. Fabi- Paulsen, “Interpreting hemoglobin and water concentration, oxygen
ani, “Measurements of scattering and absorption changes in muscle saturation, and scattering measured by near-infrared tomography of
and brain,” Philos. Trans. R. Soc. London 352共1354兲, 727–735 normal breast in vivo,” Proc. Natl. Acad. Sci. U.S.A. 100共21兲, 12349–
共1997兲. 12354 共2003兲.
25. C. E. Cooper and R. Springett, “Measurement of cytochrome oxidase 47. A. Pifferi, J. Swartling, E. Chikoidze, A. Torricelli, P. Taroni, A.
and mitochondrial energetics by near-infrared spectroscopy,” Philos. Bassi, S. Andersson-Engels, and R. Cubeddu, “Spectroscopic time-
Trans. R. Soc. London 352共1354兲, 669–676 共1997兲. resolved diffuse reflectance and transmittance measurements of the
26. J. C. Hebden, “Advances in optical imaging of the newborn infant female breast at different interfiber distances,” J. Biomed. Opt. 9共6兲,
brain,” Psychophysiology 40共4兲, 501–510 共2003兲. 1143–1151 共2004兲.
27. D. J. Hawrysz and E. M. Sevick-Muraca, “Developments toward di- 48. C. E. Cooper, C. E. ElWell, J. H. Meek, S. J. Matcher, J. S. Wyatt, M.
agnostic breast cancer imaging using near-infrared optical measure- Cope, and D. T. Delpy, “The noninvasive measurement of absolute
ments and fluorescent contrast agents,” Neoplasia 2共5兲, 388–417 cerebral deoxyhemoglobin concentration and mean optical path
共2000兲. length in the neonatal brain by second derivative near infrared spec-
troscopy,” Pediatr. Res. 39共1兲, 32–38 共1996兲. “Near real-time confocal microscopy of amelanotic tissue: detection
49. C. J. Lynn, I. S. Saidi, D. G. Oelberg, and S. L. Jacques, “Gestational- of dysplasia in ex vivo cervical tissue,” 共see comment兲 Acad. Radiol.
age correlates with skin reflectance in newborn-infants of 24–42 9共5兲, 504–512 共2002兲.
weeks gestation,” Biol. Neonat. 64共2–3兲, 69–75 共1993兲. 71. K. Carlson, M. Chidley, K. B. Sung, M. Descour, A. Gillenwater, M.
50. H. Sterenborg, S. Thomsen, S. L. Jacques, M. Duvic, M. Motamedi, Follen, and R. Richards-Kortum, “In vivo fiber-optic confocal reflec-
and R. F. Wagner, “In-vivo fluorescence spectroscopy and imaging of tance microscope with an injection-molded plastic miniature objec-
human skin tumors,” Dermatol. Surg. 21共9兲, 821–822 共1995兲. tive lens,” Appl. Opt. 44共10兲, 1792–1797 共2005兲.
51. U. Utzinger, M. Brewer, E. Silva, D. Gershenson, R. C. Blast, Jr., M. 72. R. Drezek, K. Sokolov, U. Utzinger, I. Boiko, A. Malpica, M. Follen,
Follen, and R. Richards-Kortum, “Reflectance spectroscopy for in and R. Richards-Kortum, “Understanding the contributions of NADH
vivo characterization of ovarian tissue,” Lasers Surg. Med. 28共1兲, and collagen to cervical tissue fluorescence spectra: modeling, mea-
56–66 共2001兲. surements, and implications,” J. Biomed. Opt. 6共4兲, 385–396 共2001兲.
52. R. A. Drezek, R. Richards-Kortum, M. A. Brewer, M. S. Feld, C. 73. T. S. Mang, C. McGinnis, C. Liebow, U. O. Nseyo, D. H. Crean, and
Pitris, A. Ferenczy, M. L. Faupel, and M. Follen, “Optical imaging of T. J. Dougherty, “Fluorescence detection of tumors. Early diagnosis
the cervix,” Cancer 98共9 suppl兲, 2015–2027 共2003兲. of microscopic lesions in preclinical studies,” Cancer 71共1兲, 269–276
53. T. Collier, M. Follen, A. Malpica, and R. Richards-Kortum, “Sources 共1993兲.
of scattering in cervical tissue: determination of the scattering coef- 74. B. Ebert, U. Sukowski, D. Grosenick, H. Wabnitz, K. T. Moesta, K.
ficient by confocal microscopy,” Appl. Opt. 44共11兲, 2072–2081 Licha, A. Becker, W. Semmler, P. M. Schlag, and H. Rinneberg,
共2005兲. “Near-infrared fluorescent dyes for enhanced contrast in optical
54. J. Swartling, J. Svensson, D. Bengtsson, K. Terike, and S. Andersson- mammography: phantom experiments,” J. Biomed. Opt. 6共2兲, 134–
Engels, “Fluorescence spectra provide information on the depth of 140 共2001兲.
fluorescent lesions in tissue,” Appl. Opt. 44共10兲, 1934–1941 共2005兲. 75. C. Bremer, V. Ntziachristos, and R. Weissleder, “Optical-based mo-
55. A. J. Durkin, S. Jaikumar, and R. Richardskortum, “Optically dilute, lecular imaging: contrast agents and potential medical applications,”
absorbing, and turbid phantoms for fluorescence spectroscopy of ho- Eur. Radiol. 13共2兲, 231–243 共2003兲.
mogeneous and inhomogeneous samples,” Appl. Spectrosc. 47共12兲, 76. T. Troy, D. Jekic-McMullen, L. Sambucetti, and B. Rice, “Quantita-
2114–2121 共1993兲. tive comparison of the sensitivity of detection of fluorescent and bi-
56. S. A. Prahl, Oregon Medical Laser Center Website, http:// oluminescent reporters in animal models,” Mol. Imaging 3共1兲, 9–23
[Link] 共2004兲.
57. L. Lilge, T. J. Flotte, I. E. Kochevar, S. L. Jacques, and F. Hillen- 77. T. A. Krouskop, T. M. Wheeler, F. Kallel, B. S. Garra, and D. E. Hall,
kamp, “Photoactivable fluorophores for the measurement of fluence “Elastic modulus of breast and prostate tissue under compression,”
in turbid media,” Photochem. Photobiol. 58共1兲, 37–44 共1993兲. Ultrason. Imaging 20, 260–274 共1998兲.
58. B. W. Pogue, L. Lilge, M. S. Patterson, B. C. Wilson, and T. Hasan, 78. K. Sokolov, J. Aaron, V. Mack, T. Collier, L. Coghlan, A. Gillenwa-
“The absorbed photodynamic dose examined from pulsed and cw ter, M. Follen, and R. Richards-Kortum, “Vital molecular imaging of
light using tissue-simulating dosimeters,” Appl. Opt. 36共28兲, 7257– carcinogenesis with gold bioconjugates,” Med. Phys. 30共6兲, 1539–
7269 共1997兲. 1539 共2003兲.
59. S. Jiang, B. W. Pogue, T. O. McBride, M. M. Doyley, S. P. Poplack, 79. T. S. Tkaczyk, M. Rahman, V. Mack, K. Sokolov, J. D. Rogers, R.
and K. D. Paulsen, “Near-infrared breast tomography calibration with Richards-Kortum, and M. R. Descour, “High resolution, molecular-
optoelastic tissue simulating phantoms,” J. Electron. Imaging 12共4兲, specific, reflectance imaging in optically dense tissue phantoms with
613–620 共2003兲. structured-illumination,” Opt. Express 12共16兲, 3745–3758 共2004兲.
60. S. Merritt, G. Gulsen, G. Chiou, Y. Chu, C. Deng, A. E. Cerussi, A. 80. C. J. M. Moes, M. J. van Gemert, W. M. Star, J. P. A. Marijnissen,
J. Durkin, B. J. Tromberg, and O. Nalcioglu, “Comparison of water and S. A. Prahl, “Measurements and calculations of the energy flu-
and lipid content measurements using diffuse optical spectroscopy ence rate in a scattering and absorbing phantom at 633 nm,” Appl.
and MRI in emulsion phantoms,” Technol. Cancer Res. Treat. 2共6兲, Opt. 28共12兲, 2292–2296 共1989兲.
563–569 共2003兲. 81. H. J. van Staveren, C. J. M. Moes, J. van Marle, S. A. Prahl, and M.
61. D. Li, P. M. Meaney, T. D. Tosteson, S. Jiang, T. E. Kerner, T. O. J. C. van Gemert, “Light scattering in intralipid-10% in the wave-
McBride, B. W. Pogue, A. Hartov, and K. D. Paulsen, “Comparisons length range of 400– 1100 nm,” Appl. Opt. 30共31兲, 4507–4514
of three alternative breast modalities in a common phantom imaging 共1991兲.
experiment,” Med. Phys. 30共8兲, 2194–2205 共2003兲. 82. F. A. Sayeed, H. W. Johnson, K. B. Sukumaran, J. A. Raihle, D. L.
62. W. D. D’Souza, E. L. Madsen, O. Unal, K. K. Vigen, G. R. Frank, Mowles, H. A. Stelmach, and K. R. Majors, “Stability of liposyn II
and B. R. Thomadsen, “Tissue mimicking materials for a multi- fat emulsion in total nutrient admixtures,” Am. J. Hosp. Pharm.
imaging modality prostate phantom,” Med. Phys. 28共4兲, 688–700 43共5兲, 1230–1235 共1986兲.
共2001兲. 83. C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light
63. G. Soelkner, G. Mitic, and R. Lohwasser, “Monte Carlo simulations by Small Particles, Wiley and Sons, New York 共1983兲.
and laser Doppler flow measurements with high penetration depth in 84. K. Vishwanath, B. Pogue, and M. A. Mycek, “Quantitative fluores-
biological tissuelike head phantoms,” Appl. Opt. 36共22兲, 5647–5654 cence lifetime spectroscopy in turbid media: comparison of theoreti-
共1997兲. cal, experimental and computational methods,” Phys. Med. Biol.
64. M. Larsson, W. Steenbergen, and T. Stroemberg, “Influence of optical 47共18兲, 3387–3405 共2002兲.
properties and fiber separation on laser doppler flowmetry,” J. 85. A. A. Oraevsky, R. O. Esenaliev, S. L. Jacques, and F. K. Tittel,
Biomed. Opt. 7共2兲, 236–243 共2002兲. “Laser optic-acoustic tomography for medical diagnostics: prin-
65. S. L. Jacques and S. J. Kirkpatrick, “Acoustically modulatedspeckle ciples,” Proc. SPIE 2676, 22–31 共1996兲.
imaging of biological tissues,” Opt. Lett. 23共11兲, 879–881 共1998兲. 86. E. L. Hull, M. G. Nichols, and T. H. Foster, “Quantitative broadband
66. J. P. Culver, T. Durduran, D. Furuya, C. Cheung, J. H. Greenberg, near-infrared spectroscopy of tissue-simulating phantoms containing
and A. G. Yodh, “Diffuse optical tomography of cerebral blood flow, erythrocytes,” Phys. Med. Biol. 43共11兲, 3381–3404 共1998兲.
oxygenation, and metabolism in rat during focal ischemia,” J. Cereb. 87. G. M. Hale and M. R. Querry, “Optical constants of water in the
Blood Flow Metab. 23共8兲, 911–924 共2003兲. 200-nm to 200-um wavelength region,” Appl. Opt. 12共3兲, 555–563
67. A. Bednov, S. Ulyanov, C. Cheung, and A. G. Yodh, “Correlation 共1973兲.
properties of multiple scattered light: implication to coherent diag- 88. B. W. Pogue, M. S. Patterson, H. Jiang, and K. D. Paulsen, “Initial
nostics of burned skin,” J. Biomed. Opt. 9共2兲, 347–352 共2004兲. assessment of a simple system for frequency domain diffuse optical
68. E. M. C. Hillman, D. A. Boas, A. M. Dale, and A. K. Dunn, “Lami- tomography,” Phys. Med. Biol. 40, 1709–1729 共1995兲.
nar optical tomography:demonstration of millimeter-scale depth- 89. H. Jiang, K. D. Paulsen, U. L. Osterberg, B. W. Pogue, and M. S.
resolved imaging in turbid media,” Opt. Lett. 29共14兲, 1650–1652 Patterson, “Optical image reconstruction using frequency-domain
共2004兲. data: simulations and experiments,” J. Opt. Soc. Am. A 13共2兲, 253–
69. K. Sokolov, M. Follen, and R. Richards-Kortum, “Optical spectros- 266 共1996兲.
copy for detection of neoplasia,” Curr. Opin. Chem. Biol. 6共5兲, 651– 90. B. W. Pogue, C. Willscher, T. O. McBride, U. L. Osterberg, and K. D.
658 共2002兲. Paulsen, “Contrast-detail analysis for detection and characterization
70. T. Collier, A. Lacy, R. Richards-Kortum, A. Malpica, and M. Follen, with near-infrared diffuse tomography,” Med. Phys. 27共12兲, 2693–
132. H. Shin, J. S. Temenoff, G. C. Bowden, K. Zygourakis, M. C. MEDPHOT protocol,” Appl. Opt. 44共11兲, 2104–2114 共2005兲.
Farach-Carson, M. J. Yaszemski, and A. G. Mikos, “Osteogenic 148. S. T. Flock, S. L. Jacques, B. C. Wilson, W. M. Star, and M. J. C.
differentiation of rat bone marrow stromal cells cultured on arg-gly- Vangemert, “Optical-properties of intralipid—a phantom medium
asp modified hydrogels without dexamethasone and -glycerol for light-propagation studies,” Lasers Surg. Med. 12共5兲, 510–519
phosphate,” Biomaterials 26, 3645–3654 共2005兲. 共1992兲.
133. Y. Liu, Y. L. Kim, and V. Backman, “Development of a bioengi- 149. M. D. Waterworth, B. J. Tarte, A. J. Joblin, T. van Doorn, and H. E.
neered tissue model and its application in the investigation of the Niesler, “Optical transmission properties of homogenised milk used
depth selectivity of polarization gating,” Appl. Opt. 44共12兲, 2288– as a phantom material in visible wavelength imaging,” Australas.
2299 共2005兲. Phys. Eng. Sci. Med. 18共1兲, 39–44 共1995兲.
134. A. Kienle, L. Lilge, M. S. Patterson, R. Hibst, R. Steiner, and B. C. 150. R. Bays, G. Wagnieres, D. Robert, J. F. Theumann, A. Vitkin, J. F.
Wilson, “Spatially resolved absolute diffuse reflectance measure- Savary, P. Monnier, and H. vandenBergh, “Three-dimensional opti-
ments for noninvasive determination of the optical scattering and cal phantom and its application in photodynamic therapy,” Lasers
absorption coefficients of biological tissue,” Appl. Opt. 35共13兲, Surg. Med. 21共3兲, 227–234 共1997兲.
2304–2314 共1996兲. 151. J. C. Ramella-Roman, P. R. Bargo, S. A. Prahl, and S. L. Jacques,
135. S. L. Jacques, J. R. Roman, and K. Lee, “Imaging superficial tissues “Evaluation of spherical particle sizes with an asymmetric illumina-
with polarized light,” Lasers Surg. Med. 26共2兲, 119–129 共2000兲. tion microscope,” IEEE J. Sel. Top. Quantum Electron. 9共2兲, 301–
136. S. Nickell, M. Hermann, M. Essenpreis, T. J. Farrell, U. Kramer, 306 共2003兲.
and M. S. Patterson, “Anisotropy of light propagation in human 152. U. Sukowski, F. Schubert, D. Grosenick, and H. Rinneberg, “Prepa-
skin,” Phys. Med. Biol. 45共10兲, 2873–2886 共2000兲. ration of solid phantoms with defined scattering and absorption
137. B. B. Das, K. M. Yoo, and R. R. Alfano, “Ultrafast time-gated properties for optical tomography,” Phys. Med. Biol. 41, 1823–1844
imaging in thick tissues—a step toward optical mammography,” 共1996兲.
Appl. Opt. 18共13兲, 1092–1094 共1993兲.
153. G. Mitic, J. Kolzer, J. Otto, E. Plies, G. Solkner, and W. Zinth,
138. Y. C. Guo, Q. Z. Wang, N. Zhadin, F. Liu, S. Demos, D. Calistru, A.
“Time-gated transillumination of biological tissues and tissue-like
Tirksliunas, A. Katz, Y. Budansky, P. P. Ho, and R. R. Alfano,
phantoms,” Appl. Opt. 33, 6699–6710 共1994兲.
“Two-photon excitation of fluorescence from chicken tissue,” Appl.
154. M. McDonald, S. Lochhead, R. Chopra, and M. J. Bronskill, “Multi-
Opt. 36共4兲, 968–970 共1997兲.
modality tissue-mimicking phantom for thermal therapy,” Phys.
139. A. A. Oraevsky, S. L. Jacques, and F. K. Tittel, “Measurement of
tissue optical properties by time-resolved detection of laser-induced Med. Biol. 49共13兲, 2767–2778 共2004兲.
transient stress,” Appl. Opt. 36共1兲, 402–415 共1997兲. 155. M. Firbank and D. T. Delpy, “A phantom for the testing and cali-
140. W. B. Wang, S. G. Demos, J. Ali, G. Zhang, and R. R. Alfano, bration of near-infrared spectrometers,” Phys. Med. Biol. 39共9兲,
“Visibility enhancement of fluorescent objects hidden in animal tis- 1509–1513 共1994兲.
sues using spectral fluorescence difference method,” Opt. Commun. 156. A. Kienle, M. S. Patterson, L. Ott, and R. Steiner, “Determination of
147共1–3兲, 11–15 共1998兲. the scattering coefficient and the anisotropy factor from laser Dop-
141. S. G. Demos, H. B. Radousky, and R. R. Alfano, “Deep subsurface pler spectra of liquids including blood,” Appl. Opt. 35共19兲, 3404–
imaging in tissues using spectral and polarization filtering,” Appl. 3412 共1996兲.
Opt. 7共1兲, 23–28 共2000兲. 157. S. Srinivasan, B. W. Pogue, S. Jiang, H. Dehghani, and K. D.
142. B. Quistorff and B. Chance, “Simple techniques for freeze clamping Paulsen, “Spectrally constrained chromophore and scattering NIR
and for cutting and milling of frozen tissue at low temperature for tomography improves quantification and robustness of reconstruc-
the purpose of two—or three-dimensional metabolic studies in tion,” Appl. Opt. 44共10兲, 1858–1869 共2004兲.
vivo,” Anal. Biochem. 108共2兲, 237–248 共1980兲. 158. M. S. Patterson and B. W. Pogue, “Mathimatical model for time-
143. B. M. Fenton, E. K. Rofstad, F. L. Degner, and R. M. Sutherland, resolved and frequency-domain fluorescence spectroscopy in bio-
“Cryospectrophotometric determination of tumor intravascular oxy- logical tissues,” Appl. Opt. 33共10兲, 1963–1974 共1994兲.
hemoglobin saturations: dependence on vascular geometry and tu- 159. H. Dehghani, B. W. Pogue, J. Shudong, B. Brooksby, and K. D.
mor growth,” J. Natl. Cancer Inst. 80共20兲, 1612–1619 共1988兲. Paulsen, “Three-dimensional optical tomography: resolution in
144. Z. Zhang, D. Blessington, H. Li, T. Busch, J. Glickson, Q. Luo, B. small-object imaging,” Appl. Opt. 42共16兲, 3117–3128 共2003兲.
Chance, and G. Zheng, “Redox ratio of mitochondria as an indicator 160. B. Brooksby, S. Srinivasan, S. Jiang, H. Dehghani, B. W. Pogue, K.
for the response of photodynamic therapy,” J. Biomed. Opt. 9共4兲, D. Paulsen, J. B. Weaver, C. Kogel, and S. P. Poplack, “Spectral
772–778 共2004兲. priors improve near-infrared diffuse tomography more than spatial
145. H. E. Savage, R. K. Halder, U. Kartazayeu, R. B. Rosen, T. Gayen, priors,” Opt. Lett. 30共15兲, 1968–1970 共2005兲.
S. A. McCormick, N. S. Patel, A. Katz, H. D. Perry, M. Paul, and R. 161. M. Kohl, R. Watson, and M. Cope, “Optical properties of highly
R. Alfano, “NIR laser tissue welding of in vitro porcine cornea and scattering media determined from changes in attenuation, phase, and
sclera tissue,” Lasers Surg. Med. 35共4兲, 293–303 共2004兲. modulation depth,” Appl. Opt. 36共1兲, 105–115 共1997兲.
146. J. Swartling, J. S. Dam, and S. Andersson-Engels, “Comparison of 162. C. D. Kurth, H. Liu, W. S. Thayer, and B. Chance, “A dynamic
spatially and temporally resolved diffuse-reflectance measurement phantom brain model for near-infrared spectroscopy,” Phys. Med.
systems for determination of biomedical optical properties,” Appl. Biol. 40共12兲, 2079–2092 共1995兲.
Opt. 42共22兲, 4612–4620 共2003兲. 163. R. Lohwasser and G. Soelkner, “Experimental and theoretical laser-
147. A. Pifferi, A. Torricelli, A. Bassi, P. Taroni, R. Cubeddu, H. Wab- doppler frequency spectra of a tissuelike model of a human head
nitz, D. Grosenick, M. Moller, R. Macdonald, J. Swartling, T. with capillaries,” Appl. Opt. 38共10兲, 2128–2137 共1999兲.
Svensson, S. Andersson-Engels, R. L. P. van Veen, H. Sterenborg, J. 164. S. Jiang, B. W. Pogue, K. D. Paulsen, C. Kogel, and S. P. Poplack,
M. Tualle, H. L. Nghiem, S. Avrillier, M. Whelan, and H. Stamm, “In vivo near-infrared spectral detection of pressure-induced
“Performance assessment of photon migration instruments: the changes in breast tissue,” Opt. Lett. 28共14兲, 1212–1214 共2003兲.