Problems and Solutions in Medical Physics - Nuclear Medicine
Problems and Solutions in Medical Physics - Nuclear Medicine
in Medical Physics
Series in Medical Physics and Biomedical Engineering
Series Editors: John G. Webster, E. Russell Ritenour, Slavik Tabakov, and
Kwan Hoong Ng
Recent books in the series:
Clinical Radiotherapy Physics with MATLAB®: A Problem-Solving Approach
Pavel Dvorak
Advances in Particle Therapy: A Multidisciplinary Approach
Manjit Dosanjh and Jacques Bernier (Eds)
Radiotherapy and Clinical Radiobiology of Head and Neck Cancer
Loredana G. Marcu, Iuliana Toma-Dasu, Alexandru Dasu, and Claes Mercke
Problems and Solutions in Medical Physics: Diagnostic Imaging Physics
Kwan Hoong Ng, Jeannie Hsiu Ding Wong, and Geoffrey D. Clarke
Advanced and Emerging Technologies in Radiation Oncology Physics
Siyong Kim and John W. Wong (Eds)
A Guide to Outcome Modeling in Radiotherapy and Oncology: Listening
to the Data
Issam El Naqa (Ed)
Advanced MR Neuroimaging: From Theory to Clinical Practice
Ioannis Tsougos
Quantitative MRI of the Brain: Principles of Physical Measurement,
Second Edition
Mara Cercignani, Nicholas G. Dowell, and Paul S. Tofts (Eds)
A Brief Survey of Quantitative EEG
Kaushik Majumdar
Handbook of X-ray Imaging: Physics and Technology
Paolo Russo (Ed)
Graphics Processing Unit-Based High-Performance Computing in
Radiation Therapy
Xun Jia and Steve B. Jiang (Eds)
Targeted Muscle Reinnervation: A Neural Interface for Artificial Limbs
Todd A. Kuiken, Aimee E. Schultz Feuser, and Ann K. Barlow (Eds)
Emerging Technologies in Brachytherapy
William Y. Song, Kari Tanderup, and Bradley Pieters (Eds)
Environmental Radioactivity and Emergency Preparedness
Mats Isaksson and Christopher L. Rääf
The Practice of Internal Dosimetry in Nuclear Medicine
Michael G. Stabin
Problems and Solutions
in Medical Physics
Nuclear Medicine Physics
Kwan Hoong Ng
Chai Hong Yeong
Alan Christopher Perkins
CRC Press
Taylor & Francis Group
6000 Broken Sound Parkway NW, Suite 300
Boca Raton, FL 33487-2742
This book contains information obtained from authentic and highly regarded sources. Reasonable
efforts have been made to publish reliable data and information, but the author and publisher
cannot assume responsibility for the validity of all materials or the consequences of their use. The
authors and publishers have attempted to trace the copyright holders of all material reproduced in
this publication and apologize to copyright holders if permission to publish in this form has not
been obtained. If any copyright material has not been acknowledged please write and let us know
so we may rectify in any future reprint.
Except as permitted under U.S. Copyright Law, no part of this book may be reprinted, reproduced,
transmitted, or utilized in any form by any electronic, mechanical, or other means, now known or
hereafter invented, including photocopying, microfilming, and recording, or in any information
storage or retrieval system, without written permission from the publishers.
For permission to photocopy or use material electronically from this work, please access www.
copyright.com (http://www.copyright.com/) or contact the Copyright Clearance Center,
Inc. (CCC), 222 Rosewood Drive, Danvers, MA 01923, 978-750-8400. CCC is a not-for-profit
organization that provides licenses and registration for a variety of users. For organizations
that have been granted a photocopy license by the CCC, a separate system of payment has been
arranged.
v
vi Contents
7 Radionuclide Therapy 93
7.1 Sealed and Unsealed Source Therapy 93
7.2 Therapeutic Procedures in Nuclear Medicine 94
7.3 Hyperthyroidism Absorbed Dose Calculation 94
7.4 I-131 Treatment Guidelines 95
7.5 Calculation of Administered Activity for I-131 Treatment 96
7.6 Radioiodine Therapy and Pregnancy 97
7.7 Safe Administration of I-131 98
7.8 Guidance Level for Hospitalisation of I-131 Patients 99
7.9 Radioembolisation 99
7.10 Radioimmunotherapy (RIT) 100
Bibliography 139
About the series
• Artificial organs
• Assistive technology
• Bioinformatics
• Bioinstrumentation
• Biomaterials
• Biomechanics
• Biomedical engineering
• Clinical engineering
• Imaging
• Implants
• Medical computing and mathematics
• Medical/surgical devices
• Patient monitoring
• Physiological measurement
• Prosthetics
• Radiation protection, health physics and dosimetry
• Regulatory issues
• Rehabilitation engineering
• Sports medicine
• Systems physiology
• Telemedicine
• Tissue engineering
• Treatment
xi
xii About the series
xv
Authors
Kwan Hoong Ng, PhD, FinstP, DABMP, received his MSc (medical phys-
ics) from the University of Aberdeen and PhD (medical physics) from the
University of Malaya, Malaysia. He is certified by the American Board of
Medical Physicists. Professor Ng was honoured as one of the top 50 medical
physicists in the world by the International Organization of Medical Physics
(IOMP) in 2013. He also received the International Day of Medical Physics
Award in 2016. He has authored/co-authored over 230 papers in peer-reviewed
journals and 25 book chapters and has co-edited 5 books. He has presented
over 500 scientific papers and has more than 300 invited lectures. He has also
organised and directed several workshops on radiology quality assurance,
digital imaging and scientific writing. He has directed research initiatives in
breast imaging, intervention radiology, radiological safety and radiation dosim-
etry. Professor Ng serves as a consultant for the International Atomic Energy
Agency (IAEA) and is a member of the International Advisory Committee of
the World Health Organization (WHO), in addition to previously serving as a
consulting expert for the International Commission on Non-Ionizing Radiation
Protection (ICNIRP). He is the founding and emeritus president of the South
East Asian Federation of Organizations for Medical Physics (SEAFOMP) and
is a past president of the Asia-Oceania Federation of Organizations for Medical
Physics (AFOMP).
xvii
xviii Authors
xix
List of abbreviations
A mass number
ACD annihilation coincidence detection
C-11 carbon-11
CFOV central field-of-view
COR centre of rotation
Cr-51 chromium-51
CT computed tomography
CZT cadmium zinc telluride
ECG electrocardiogram
ELISA enzyme-linked immunosorbent assay
F-18 fluorine-18
FDA Food and Drug Administration
FDG fluorodeoxyglucose
FOV field-of-view
FWHM full width at half maximum
He helium
HVL half value layer
I-125 iodine-125
I-131 iodine-131
IAEA International Atomic Energy Agency
ICRP International Commission on Radiological Protection
In-111 indium-111
LEHR low-energy high-resolution
LOR line-of-response
LSF line spread function
mAbs monoclonal antibodies
MIRD medical internal radiation dosimetry
Mo-99 molebdenum-99
MRI magnetic resonance imaging
MTF modulation transfer function
MUGA multiple-gated acquisition
NEMA National Electrical Manufacturers Association
xxi
xxii List of abbreviations
N-13 nitrogen-13
NaI sodium iodide
NaI(Tl) sodium iodide doped with thallium
O-15 oxygen-15
O-18 oxygen-18
P-32 phosphorus-32
Pb lead
PET positron emission tomography
PHA pulse height analyser
PMT photomultiplier tube
Po-210 polonium-210
PSF point spread function
PSPMT position sensitive photomultiplier tube
QC quality control
Ra-223 radium-223
Re-188 rhenium-188
RIA radioimmunoassay
RIT radioimmunotherapy
ROI region-of-interest
RPO radiation protection officer
Sm-153 samarium-153
SPECT single photon emission computed tomography
SUV standardised uptake value
t1/2 half-life
Te effective half-life
Tb biological half-life
Tp physical half-life
Tc-99m technetium-99m
TOF time-of-flight
TVL tenth value layer
U-235 uranium-235
UFOV useful field-of-view
Y-90 yittrium-90
Z atomic number
List of physical
constants
xxiii
Radioactivity
and Nuclear
Transformation
1
1.1 NUCLEAR STABILITY CURVE
PROBLEM
Figure 1.1 shows the nuclear stability curve. Name the types of radioactive
decays that tend to happen at regions A, B and C.
Solution
A: Beta decay
B: Positron emission or electron capture
C: Alpha decay
1
2 Problems and Solutions in Medical Physics
Solution
Alpha decay is the spontaneous emission of an alpha particle (24 He or α)
from an unstable nucleus (Figure 1.2). It typically occurs for heavy nuclides
(A > 150 or Z > 83) when the nuclear binding energy can no longer hold the
nucleons together. Alpha decay is often followed by gamma or characteris-
tic X-ray emission resulting from the competing processes, such as internal
conversion and Auger electron emission. An alpha particle has two protons
and two neutrons and carries an electronic charge of 2+. Alpha decay can be
expressed by the following equation:
A
Z X → AZ−−42Y + 24 He + transition energy
Solution
When a nucleus has an excessive number of neutrons compared to protons
(i.e. high N/Z ratio), the excess neutron will be converted into a proton, an
electron and an antineutrino (ν e ). The proton remains in the nucleus, but the
electron is emitted as a beta (β−) particle (Figure 1.3). Beta decay increases the
number of protons by 1 and hence transforms the atom into a different element
with atomic number Z + 1. The mass number remains unchanged because of
simultaneous decrease of a neutron only. Beta decay can be expressed by the
following equation:
A
Z X → Z +A1Y + β − +ν e + energy
Solution
When a nucleus has excess of protons compared to neutrons (i.e. low N/Z
ratio), the excessive proton will be converted into a neutron, a positron (β+),
and a neutrino (ν) (Figure 1.4). The neutron will remain in the nucleus, but a
positron will be emitted. Positron decay reduces the number of protons by 1
and hence transforms the atom into a different element with atomic number
Z–1. The mass number remains unchanged because of the simultaneous incre-
ment of a neutron only.
4 Problems and Solutions in Medical Physics
A
Z X → Z −A1Y + β + + ν e + energy
Solution
Electron capture is an alternative process to positron decay when the nucleus
has excessive number of protons compared to neutrons (i.e. low N/Z ratio).
During electron capture, the nucleus captures an inner orbital (i.e. K- or
L-shell) electron, thereby changing a proton to a neutron with the simultane-
ous emission of a neutrino (νe). The capture of an electron from the inner shell
creates a vacancy in the orbit, which is then filled by an outer shell electron,
and the energy difference between the electron shells results in the emission of
a characteristic X-ray (Figure 1.5).
Electron capture can be expressed by the following equation:
A
Z X + e − → Z −1AY + νe + energy
1 • Radioactivity and Nuclear Transformation 5
Solution
During radioactive decay, a daughter is often formed in an excited (metastable
or isomeric) state. Gamma rays are emitted when the daughter nucleus under-
goes an internal rearrangement and transitions from the excited state to a lower
energy state (Figure 1.6). This process is called isomeric transition. Therefore,
isomeric transition is a decay process that yields gamma radiation without the
emission or capture of a particle by the nucleus. There is no change in atomic
number (Z), mass number (A), or neutron number (N). Thus, this decay mode
is isobaric and isotonic, and it occurs between two nuclear energy states with
no change in the neutron-to-proton (N/Z) ratio. Isomeric transition can be
expressed by the following equation:
A
Z X * → ZA X + Energy
Solution
An alpha particle has the least penetrability, followed by beta and gamma radi-
ation. Alpha can be stopped by a sheet of paper or skin, beta by aluminium,
and gamma radiation by very thick concrete or lead. The penetrability of a
neutron is dependent on its energy. Neutrons may penetrate through concrete
and lead but can be slowed down by water or low atomic number materials,
such as boron.
Solution
ln 2
λ= = 5.795 × 10 −8 s−1
138.4 d × 24 h × 60 min × 60 s
A = 370 MBq = 3.7 × 108 Bq = 3.7 × 108 dps
A = λN
Solution
10000 counts
=
• Sample plus background count rate, A = 1000 cpm
10 min
1296 counts
=
• Background count rate, B = 216 cpm
6 min
• Net sample count rate = 1000 – 216 cpm = 784 cpm
1000 counts
=
• Standard deviation for A = 10 cpm
10 min
8 Problems and Solutions in Medical Physics
216 counts
=
• Standard deviation for B = 6 cpm
6 min
Solution
• Net standard count = 85500 − 952 = 84548
• Standard deviation for the standard counts, σ s = 85500 + 952 = 294
• Net thyroid count = 42100 – 2200 = 39900
• Standard deviation for the thyroid counts, σ t = 42100 + 2200 = 210
39900
• % thyroid uptake = 84548 × 100% = 47.2%
( 84548 ) + ( 39900 )
2 2
• Standard deviation of the uptake = 39900 × 294 210
84548
= 0.472 × 0.00001209 + 0.00002770
= 0.472 × 0.006308
= 0.003
• % standard deviation of the uptake = 00..472
003 × 100% = 0.6%
Solution
Decay constant, λ = 0.25 h−1
0.693
λ =
t1 2
0.693 0.693
Therefore, the physical half-life, t1 2 = = = 2.77 h
λ 0.25 h −1
Solution
1
activity equals to four half-lives.
16
24 h
= 6h
4 half-lives
The physical half-life of the source is six hours. The source is most prob-
ably Tc-99m, which is the most common radionuclide used in nuclear
medicine.
10 Problems and Solutions in Medical Physics
Solution
a. The effective half-life is defined as the time interval required to
reduce the radioactivity level of an internal organ or the whole body
to one-half of its original activity due to radioactive decay (physical
half-life) and biological elimination (biological half-life). Effective
half-life, Te, is expressed as the following:
1 1 1
= +
Te Tp Tb
or
T p × Tb
Te =
T p + Tb
where:
Tp = physical half-life, which is defined as the time interval
required for a radioactive substance to decay to half of its
original activity.
Tb = biological half-life, which is defined as the time interval
required for an organ or the whole body to eliminate 50% of
its original activity by normal routes of elimination: metabolic
turnover and excretion.
T p × Tb
b. Te =
T p + Tb
6h × 4h
Therefore, Te =
6h + 4h
Te = 2.4 h
1 • Radioactivity and Nuclear Transformation 11
Note: Effective half-life is always smaller than the biological half-life and
physical half-life. For a permanently implanted radionuclide (which no
biological half-life is applied), the effective half-life is the same as the physical
half-life.
Solution
• Total biological elimination = 10 + 25 + 3 = 38%
• Given that 38% elimination in six hours.
• To achieve 50% elimination (biological half-life, Tb), it takes:
Tb = 50
38 × 6 h = 7.89 h
Tb × T p 7.89 h × 6 h
Te = = = 3.4 h
Tb + T p 7.89 + 6 h
Solution
Given:
dN
A = − λN =
dt
dN
− = λ dt
N
Nt t
dN
−
∫
N0
N ∫
= λ dt
0
N
ln t = −λt
N0
Nt
= e( )
− λt
N0
N t = N 0e(
− λt )
Solution
Given:
A0 = 1110 MBq
Tp1/ 2 = 6.0 2 h
50%
Tb1/ 2 = × 2h = 5 h
20%
T × Tb 6.02 × 5
Te = p = = 2.73 h
Tp + Tb 6.02 + 5
ln 2 ln 2
λ= = = 0.254 h −1
Te 2.73 h
A = A 0e − λ t
−0.254 h −1 ( 2 h )
= 1110 MBq e
= 667.88 MBq
Therefore, the remaining activity at the time of imaging is approximately
668 MBq.
Solution
T p × Tb 8 d × 2 d
Te = = = 1.6 d
T p + Tb 8 d + 2 d
ln 2 ln 2
λ= = = 0.433 d −1
Te 1.6 d
A = A 0e − λt
14 Problems and Solutions in Medical Physics
−0.433 d −1 ( 3 d )
= 3700 MBq e
= 1009 MBq
Therefore, the remaining activity in patient’s body on the third day of the treat-
ment is approximately 1009 MBq.
Solution
ln 2
λ=
t1/ 2
ln 2
λ1 = = 0.0525 h −1
13.2 h
ln 2
λ2 = = 0.1155 h −1
6h
e(
−0.0525t )−( −0.1155t )
=5
e 0.063t = 5
0.063t = 1.609
t = 25.5 h
1 • Radioactivity and Nuclear Transformation 15
Therefore, 370 MBq I-123 and 1850 MBq Tc-99m will achieve the same activ-
ity at 25.5 hours later.
1.19 ATTENUATION
PROBLEM
The mass attenuation coefficient of bone with a density of 1.8 g · cm−3 is
0.2 cm2g−1 for an 80 keV gamma ray. Calculate the percentage of a photon
beam attenuated by 5 cm thickness of bone.
Solution
µ
Mass attenuation coefficient =
ρ
µ
= 0.2 cm 2g −1
1.8 g ⋅ cm −3
μ = 0.36 cm−1
Fraction transmitted:
Nt
= e−µ x
N0
Nt −0.36 cm −1 ( 5 cm )
=e = 0.165
N0
Solution
a. The specific gamma ray constant (Γ) is the exposure rate at a specific
distance from a given amount of a photon-emitting radionuclide.
These constants are used frequently for dosimetry and radiation
protection purposes. The SI unit for Γ is Ckg−1s−1Bq−1 at 1 m, but it
is often expressed in its traditional units, R h−1mCi−1 at 1 cm.
b. A high Γ means there are large numbers of gamma rays emitted
and available to form the image. Some radionuclides decay by more
than one mechanism, so the radionuclides that produce more useful
gamma rays are preferable for nuclear medicine imaging.
Solution
a. To calculate alpha particle range in air for 4 < E < 8 MeV, the
following equation is used:
R (cm) = 0.325 E3/2 (MeV)
R = 0.325(6.3)3/2
R = 5.14 cm
0.001293 11.81
Rplastic = (5.14 cm)
1.13 14.8
0.00444
Rplastic = (5.14 cm )
4.3472
Rplastic = 0.0053 cm
0.001293 12.88
Rmylar = (5.14 cm )
1.4 14.8
0.00464
Rmylar = (5.14 cm )
5.3859
Rmylar = 0.0044 cm
0.001293 13.6
Rhuman tissue = (5.14 cm )
1.04 14.8
0.00477
Rhuman tissue = (5.14 cm )
4.0009
Solution
The criteria of an ideal radiopharmaceutical for diagnostic nuclear medicine
include:
19
20 Problems and Solutions in Medical Physics
Solution
The criteria of an ideal radiopharmaceutical for therapeutic nuclear medicine
include:
Solution
a. Tc-99m is a colourless metallic chemical element. It has an atomic
number, Z, of 43, and mass number, A, of 56. It is metastable (as
indicated by the symbol ‘m’) and thereby decays into Tc-99 through
gamma emission (88%) and isomeric transition (12%). It has a physi-
cal half-life of 6.02 hours and emits 140 keV gamma rays, which
2 • Radionuclide Production and Radiopharmaceuticals 21
2.4 CYCLOTRON
PROBLEM
With the aid of a diagram, explain the basic principles of radionuclide produc-
tion in a cyclotron. Use In-111 as an example.
Solution
Figure 2.2 shows the basic structure of a cyclotron. A cyclotron consists of
two large dipole magnets (known as Dees) that produce a uniform magnetic
field. An oscillating voltage is applied to produce an electric field across the
gap between the Dees. Charged particles, such as protons and electrons, are
injected into the magnetic field region of the Dees. The electric field in the
gap then accelerates the particles as they pass through it. The particles now
have higher energy so they follow a semi-circular path in the next Dee with
a larger radius and then reach the gap again. The field in the gap accelerates
them, and they enter the first Dee again. Thus, the particles gain energy as
they spiral outwards until they gain sufficient velocity and are deflected into
a target.
22 Problems and Solutions in Medical Physics
In the production of In-111, the accelerated particles are protons. The tar-
get atoms are Cd-111. When a proton enters the nucleus of a Cd-111 atom, the
Cd-111 is transformed into In-111 by discharging a neutron. This interaction in
the form of nuclear reaction is as follows:
Target atom (bombarding particle, emitted particle) product nuclide
For example:
2.5 CYCLOTRON-PRODUCED
RADIONUCLIDES
PROBLEM
Name three radionuclides other than In-111 that are produced by a cyclotron.
State their nuclear reaction equations and physical half-lives.
Solution
Solution
Nuclear fission is the splitting of a heavy nuclide into daughter products
(smaller atoms) with the release of energy (Figure 2.3).
2.7 REACTOR-PRODUCED
RADIONUCLIDES
PROBLEM
• Name three radionuclides other than Mo-99 that are produced by a
nuclear reactor.
• State their nuclear reaction equations and physical half-lives.
Solution
2.8 Mo-99/Tc-99m
RADIONUCLIDE GENERATOR (I)
PROBLEM
Sketch a labelled diagram of a Mo-99/Tc-99m radionuclide generator and
describe the process of Tc-99m elution.
Solution
Figure 2.4 shows the basic structure of a Mo99/Tc-99m generator. The
generator is heavily shielded in a lead container known as the ‘lead pig’.
Chemically purified Mo-99 is loaded in the form of ammonium molybdenate
(NH4+ MoO4−) onto a porous column containing 5–10 g of alumina (Al2O3)
resin. The ammonium molybdenate is attached to the surface of the alumina
2.9 Mo-99/Tc-99m
RADIONUCLIDE GENERATOR (II)
PROBLEM
Differentiate between the wet and dry systems of the Mo-99/Tc-99m
generator.
Solution
The wet system has a saline reservoir connected by tubing to one end of the
column, and an evacuated vial draws the saline through the column, hence
keeping the column wet. During eluation, a vacuum vial is placed at the exit
port to collect the Tc-99m eluent.
The dry system does not have a saline reservoir in the generator; therefore,
it requires a small vial containing normal saline to be attached at the entry port
and a vacuum vial at the collection port during eluation.
Note: Evacuated (vacuum) vials are used since this creates a condition of nega-
tive pressure during generator elution. This prevents the egress of solution and
the possibility of area contamination should any parts of the generator or tub-
ing develop a puncture or leak.
2 • Radionuclide Production and Radiopharmaceuticals 27
Solution
a. Mo-99 is the parent radionuclide of Tc-99m and is bound to the
alumina column in the generator. It is essential that no Mo-99 sepa-
rates from the column during elution, as this would both degrade the
gamma camera image and contribute unnecessary additional radia-
tion doses to the patient. Mo-99 has a physical half-life of 66 hours
and emits high energy of 740 keV gamma rays.
b. Al3+ is a trace metal that should not be administered into the patients.
It will interfere with the preparation of Tc-99m radiopharmaceuti-
cals, such as Tc-99m-labelled colloid and Tc-99m bone-scanning
agents by forming larger particles, which will be trapped in the
small capillaries of the lungs following intravenous injection. It will
also agglutinate the red blood cells during labelling.
Note: According to the FDA rules and regulations, the concentration of
aluminium ion should not be more than 10 μg/mL of the generator eluate.
dN d (2.1)
= λ p N p − λd N d
dt
28 Problems and Solutions in Medical Physics
λd ( Ap )0 − λ pt −λd t
( Ad )t = λd N d = (e −e ) (2.2)
(λd − λ p )
where:
Np = number of particles of the parent radionuclide
Nd = number of particles of the daughter radionuclide
λp = decay constant for the parent radionuclide
λd = decay constant for the daughter radionuclide
Ap = activity of the parent radionuclide
Ad = activity of the daughter radionuclide
t = time
From Equations 2.1 and 2.2, derive an equation that shows transient equilib-
rium. What is the condition for transient equilibrium to happen?
Solution
If λd > λp, that is (t1/2)d < (t1/2)p, then e−λdt is negligible compared to e−λpt when
t is sufficiently long.
The equation becomes:
λd ( Ap )0
( λd − λ p ) ( )
( Ad )t = e −λ pt
λd
= ( Ap )t
( λd − λ p )
0.693
(t1/ 2 ) p
= (A )
0.693 0.693 p t
−
(t1/ 2 ) p (t1/ 2 )d
(t1/ 2 ) p
= ( Ap ) t
(t1/ 2 ) p − (t1/ 2 )d
Solution
i. Curve X indicates the radioactive decay of the parent radionuclide,
that is, Mo-99.
ii. Curve Y indicates the radioactivity of the daughter (i.e. Tc-99m) at
subsequent time intervals.
iii. The plot demonstrates the phenomenon of transient equilibrium.
The activity of Mo-99 (curve X) and Tc-99m (curve Y) versus time
were plotted on a logarithm graph. The activity of the daughter,
Tc-99m, initially builds up as a result of the decay of the parent,
30 Problems and Solutions in Medical Physics
PROBLEM
Mo-99 (t1 /2 = 66 hours) and Tc-99m (t1/2 = 6 hours) are in transient equilibrium
in a Mo generator. If 1000 MBq of Mo-99 is present in the generator, what
would be the activity of Tc-99m after 82 hours, assuming that 87% of Mo-99
decays to Tc-99m?
Solution
0.693
− ( 82 h )
( Ap ) t = ( Ap ) 0 e 66 h
= 1000 MBq e(
−0.861)
= 422.74 MBq
(t1 2 ) p
( Ad )t = ( ( Ap )t ( 0.87 )
(t1 2 ) p − (t1 2 )d
66
= ( 422.74 MBq)(0.87)
60
= 404.56 MBq
2 • Radionuclide Production and Radiopharmaceuticals 31
Solution
When λd is much greater than λp, that (t1/2)d is much shorter than (t1/2)p, the λd
can be neglected compared to λp. Equation 2.2 now becomes:
λd ( Ap )0
( Ad )t =
λd
(e )−λ pt
( Ad )t = ( Ap )0 ( e −λ t )
p
= ( Ap ) t
The condition for secular equilibrium to happen is when the physical half-life
of the parent is much longer than the daughter, by a factor of 100 or greater.
Solution
Figure 2.6 illustrates the method for the Mo-99 breakthrough test. The vial
containing the Tc-99m eluate is placed in a thick (~6 mm) lead container
and then placed into a radionuclide calibrator. The high-energy photons of
Mo-99 can be detected, whereas virtually all of the Tc-99m 140 keV photons
are attenuated by the lead container. The acceptance limit for Mo-99 break-
through is 5.55 kBq Mo-99 per 37 MBq of Tc-99m eluate.
32 Problems and Solutions in Medical Physics
Note: This procedure should be undertaken on the first elution from each new
generator and at any time when it is suspected that the generator has not been
used according to accepted normal standard procedures. It is also good prac-
tice to check for breakthrough at the end of use of the generator to confirm the
integrity of the column.
2.16 PREPARATION OF
RADIOPHARMACEUTICAL
PROBLEM
A vial containing 5 mL of Tc-99m-MDP is prepared in the nuclear medicine
hot-lab for radionuclide bone imaging. At 8.00 a.m., the calibrated activity
of the radiopharmaceutical is 3700 MBq. What volume of the preparation is
needed to give a 740 MBq injection at 12.00 p.m., given the physical half-life
of Tc-99m is six hours?
2 • Radionuclide Production and Radiopharmaceuticals 33
Solution
Decay constant λ for Tc-99m:
0.693 0.693
λ = = = 0.115 h −1
t1/ 2 6h
Using the decay equation,
N = Noe−λt
N = 3700 e
( )
− 0.115 h −1 ( 4 h )
= 2335 MBq
5 mL → 2335 MBq
? mL → 740 MBq
740 MBq
Volume of preparation = × 5 mL = 1.58 mL
2335 MBq
2.17 ADMINISTRATION OF
RADIOPHARMACEUTICALS
PROBLEM
A patient was given a radioactive meal labelled with 1110 MBq In-111 for a
colonic transit study. The patient was imaged 24 hours later. Assuming that
none of the activity was excreted, what is the remaining total activity at the
time of imaging, given the physical half-life for In-111 is 67.9 hours?
Solution
A = Aoe −λt
Ao = 1110 MBq
0.693
λ= = 0.01 h −1
67.9 h
A = 1110 MBq e
( −0.01 h−1 )( 24 h ) = 873 MBq
34 Problems and Solutions in Medical Physics
Note: Non-absorbable tracers are commonly used for gastric emptying and
gastrointestinal transit.
Solution
2.19 F-18-Fluorodeoxyglucose
PROBLEM
What is F-18-FDG? Explain the metabolism of F-18-FDG in human tissues
that makes it a useful radiopharmaceutical in detecting and staging of malig-
nant tumours.
Solution
F-18-FDG stands for fluorine-18-flurodeoxyglucose. It is a glucose ana-
logue, with the positron emitting radionuclide, F-18 substituted for the
normal hydroxyl group (-OH) at the C-2 position in the glucose molecule.
Figure 2.7 shows the difference between the F-18-FDG and glucose mol-
ecule structures.
2 • Radionuclide Production and Radiopharmaceuticals 35
Solution
Radiochemical impurities may be due to:
Solution
‘Dead time’ of a counting system refers to the period where the counter cannot
process a second event after the first event has been detected. Typical values of
dead time are between 250 and 1000 ns.
37
38 Problems and Solutions in Medical Physics
Solution
Figure 3.1 shows the differences between paralysable and non-paralysable
counting systems. The black bars indicate processed signals, whereas the grey
bars indicate non-processed signals. In non-paralysable systems, a fixed dead
time, τ, follows each event that occurs during the live period of the detector.
Signals that arrive during the dead time are not recorded and have no effect on
the system (i.e. signals 3 and 5). However, in paralysable systems, in addition
to the usual dead time effect (signals 3 and 5), the arrival of a secondary signal
before complete processing of the previous signals will extend the dead time,
causing the secondary signal to be rejected (i.e. signal 6).
Note: In extreme cases, where the count rate is high, paralysable systems can
be completely unresponsive (switched off) as the signals and dead times over-
lap and no signal is registered, hence the term ‘paralysed’.
FIGURE 3.2 Observed count rate versus input count rate of different types of
detection systems.
Solution
B is the non-paralysable system. The dead time increases with increasing input
count rate and becomes saturated in the plateau region.
C is the paralysable system. As the input count rate increases beyond the
saturation region, the system becomes paralysed and loses its detection effi-
ciency rapidly beyond the saturation region.
Solution
Figure 3.3 illustrates the basic operational principle of a gas-filled detector.
A gas-filled detector contains a volume of compressed gas (normally a noble
gas such as helium or argon) between two electrodes. A potential difference
(voltage) is applied between the electrodes. When ionising radiation interacts
with the atoms of the gas, it will form ion pairs due to the ionisation process.
The positive ions (cations) are attracted to the negative electrode, whereas the
40 Problems and Solutions in Medical Physics
negative ions (anions) are attracted to the positive electrode. In most detectors,
the cathode is the wall of the container, and the anode is an insulated wire
placed at the centre of the container. When the electrons reach the anode, they
will travel through the circuit to the cathode, where they recombine with the
cations. The electrical current that is formed by the flow of the electrons can
then be measured with a sensitive ammeter. The response is proportional to
the ionisation rate. Hence, the activity, exposure or dose rate can be calculated.
Solution
Noble gases, such as argon and xenon, are commonly used in gas-filled detec-
tors because of their inert nature. Chemical reactions will not occur within a
noble gas following ionisation events. Such reactions could change the charac-
teristics and performance of the detector.
3 • Non-imaging Detectors and Counters 41
Note: The noble gas is normally placed under pressure within the chamber.
Loss of pressure can affect the performance of the detector.
FIGURE 3.4 Plot of number of ion pairs formed versus applied voltage.
Solution
a. Region A: Recombination region
The applied voltage is relatively low in this region. When a
small voltage is applied, ion pairs form where electrons are
attracted to the anode, and positive ions are attracted to the
cathode before they are recombined. As the voltage increases,
more ions are collected and fewer ions are recombined. The cur-
rent increases as the voltage increases until a saturation point.
42 Problems and Solutions in Medical Physics
Solution
a. The activity calibrator is a well-typed ionisation chamber used to
measure the activity of administered radiopharmaceuticals. The
activity calibrator does not measure the radiation dose contributed
by the radionuclide, but it measures the amount of radioactivity con-
tained within the ‘dose’ of a radiopharmaceutical drawn up to be
administered to a patient.
Note: An activity calibrator is commonly known as the ‘dose’ calibrator.
Strictly speaking, the term ‘dose’ should not be used since dose is measured in
Gy; however, it is a common term still used by many people.
b. Figure 3.5 shows the internal structure of an activity calibrator.
An activity calibrator comprises an ionisation chamber, high-voltage
power supply, electrometer, a radionuclide selector, and a display
unit. The ionisation chamber contains pressurised gas (such as argon
gas), and the hermetically sealed chamber contains two co-axial
cylindrical electrodes, which are connected to a high-voltage supply.
The wall of the chamber is connected to the cathode, and the collector
electrodes are connected to the anode. When a sample is placed into
the chamber, the gas is ionised. The ion pairs then migrate towards
the anode and cathode and generate an electrical current. The current
is proportional to the activity of the sample. An electrometer is used
to measure the very small current (about a few μA). A calibration
factor is applied according to the selected radionuclide to convert
the current to activity. The output is usually displayed in MBq or
mCi. The activity calibrator is operated in current mode; therefore,
dead-time effect is minimal. It can accurately assay activity up to
74 GBq (or 2 Ci); however, it is relatively insensitive and cannot
measure activity less than 3.7 × 105 GBq (or 1 μCi).
44 Problems and Solutions in Medical Physics
Solution
Factors that may affect the measurement accuracy of an activity calibrator
include:
Solution
The common quality control tests include:
Solution
Scintillators are materials that emit visible or ultraviolet light after absorb-
ing ionising radiation. When ionising radiation interacts with a scintillator,
the electrons of the scintillator atoms are raised to an excited energy level.
Eventually, these electrons will fall back to its original energy state, with the
emission of visible or ultraviolet light. A scintillation detector is obtained
when a scintillator is coupled to an electronic light sensor, such as a photomul-
tiplier tube (PMT), photodiode or silicon photomultiplier. The PMT absorbs
the light emitted by the scintillator and reemits it in the form of electrons via
a photoelectric effect. The subsequent multiplication of these electrons results
in an electrical pulse, which can then be analysed and correlate to the inten-
sity of the radiation that struck the scintillator. In all scintillators, the amount
of light emitted increases in direct proportion to the energy deposited in the
scintillator.
Note: Most scintillators have more than one mode of light emission, and each
mode has its characteristic decay constant. For example, luminescence is the
emission of light after excitation, fluorescence is the prompt emission of light,
and phosphorescence (or commonly known as afterglow) is the delayed emis-
sion of light and varies from a few nanoseconds to hours depending on the
material.
Solution
• Gamma well counter (sample counter) for measuring activity within
standard preparations and biological samples, such as blood and
urine.
• Probe detector, such as a thyroid probe detector, for measuring
organ uptake.
• Gamma camera (or scintillation camera) for diagnostic imaging.
3 • Non-imaging Detectors and Counters 47
Solution
Figure 3.6 shows the schematic diagram of a gamma well counter. A gamma well
counter consists of a well-shaped scintillation detector (usually NaI[TI] crystal)
with a hole in the centre. A sample is placed inside the hole for increased geo-
metric and counting efficiency. This design gives extremely high efficiency and
allows detection of very low activity, that is, less than 37 Bq (or 1 nCi). The scin-
tillation detector is coupled to a photomultiplier tube (PMT) and connected to a
pre-amplifier and amplifier. When ionising radiation interacts with the scintillation
crystal, lights are emitted and converted to electrons by a thin photocathode layer in
front of the PMT. The electrons undergo a multiplication process in the PMT and
are then sent to a pre-amplifier and amplifier to amplify the electronic signals. The
signals are then sent to a pulse height analyser to discriminate the energy informa-
tion. Only the signals in the selected energy window are used for counting.
Solution
i. Glomerular function rate (GFR) test
The test uses a radiopharmaceutical that is rapidly excreted by the
kidneys, such as Cr-51-EDTA or Tc-99m-DTPA. The total renal
function is measured from the clearance of the activity in a series of
blood samples taken following intravenous injection of the radio-
pharmaceutical. The activity of the blood samples is counted using
a gamma well counter. This is a valuable test in patients undergoing
chemotherapy with cytotoxic drugs that may result in renal toxicity.
ii. Plasma and red blood cells volume determination
Plasma and red blood cells volume can be determined using iso-
tope dilution method. By diluting a known activity and volume
of radiotracer in an unknown volume (plasma or red blood cells),
and then measuring the activity in a volume withdrawn after
adequate mixing, the unknown volume can be calculated. For
plasma volume determination, I-125-Human Serum Albumin is
used; whereas for red blood cell volume determination, Cr-51 is
commonly used. The most common clinical application of this
test is when evaluating patients for polycythemia vera.
Note: Several precautions need to be taken during the test. For example, samples
should be taken with the patient in the same position each time, the radiotracer
must be stably bound and the volume must be small enough as to not affect the
volume being measured. In addition, equilibrium must be reached before samples
are drawn; usually a 20-minute sample is used. However, it should be noted that
the rate of mixing varies from patient to patient, especially in conditions such as
shock or elevated haematocrits in which mixing may be delayed considerably.
iii. Radioimmunoassay (RIA)
RIA is an in-vitro assay technique used to measure the concen-
trations of an antigen (e.g. hormone levels in the blood) with
very high sensitivity. In this technique, a known quantity of
3 • Non-imaging Detectors and Counters 49
Solution
Figure 3.7 shows the schematic diagram of a thyroid probe. A thyroid probe
consists of a cylindrical NaI(Tl) crystal coupled to a photomultiplier tube
(PMT) which in turn is connected to a pre-amplifier and other associated
electronics, such as amplifier, pulse height analyser (PHA), timer, scaler, rate
meter and a recorder. The operation principle is similar to a gamma well coun-
ter except that a long cylindrical or conical collimator is applied on the thy-
roid probe to limit the field-of-view (FOV) of the thyroid. This reduces the
background radiation from areas outside the thyroid to reach the detector. The
probe is connected to a high-voltage power supply.
50 Problems and Solutions in Medical Physics
When ionising radiation interacts with the NaI(Tl) crystal, electrons are raised
to an excited energy state. Eventually the electrons will return to their ground
energy state by emitting the excessive energy in terms of visible or ultravio-
let light. The light photons will interact with the photocathode located at the
entrance of the PMT and are converted to electrons. The electrons are then
accelerated and focused towards a series of dynodes in the PMT. The increas-
ing voltage electrons will then impinge the anode. The resultant output signal
will then be amplified and analysed by the PHA. Only the energies in the
selected energy range will be used for counting. The timer, scaler and rate
meter are used for counting of the radioactivity signals.
Solution
Thyroid uptake measurement can be performed using one or two capsules
of I-123 or I-131 sodium iodide. A neck phantom, consisting of a Lucite or
Perspex cylinder of diameter similar to the neck and containing a hole parallel
to its axis for radioiodine capsule placement is required.
3 • Non-imaging Detectors and Counters 51
Two-capsule method:
Two radioiodine capsules with almost identical activities are used. Firstly,
each capsule is placed in the neck phantom and counted. Secondly, one cap-
sule is given to the patient, and another one is used as the ‘standard’. The
gamma emissions from the patient’s neck are counted at 4–6 hours after
administration and repeated at 24 hours post-administration. During each
measurement, the patient’s distal thigh and background activities are also
obtained with the same duration of counting. The measurements are repeated
with the ‘standard’ capsule placed in the neck phantom for counting. Finally,
the uptake is calculated for each neck measurement as follows:
One-capsule method:
A single radioiodine capsule is used. The capsule is first counted in the
neck phantom, then ingested by the patient. As in the two-capsule method,
the patient’s neck, distal thigh and background activities are counted for the
same period of time at 4–6 hours and again at 24 hours post-administration.
The times of the capsule administration and neck counts are recorded.
Finally, the uptake is calculated as follows:
where:
t1/2 = physical half-life of the radionuclide
t = time elapsed between the measurement of capsule in the phantom
and the patient’s neck
Note: The single-capsule method reduces the cost of the examination and
requires fewer measurements, but it is more susceptible to instability of the
equipment, technologist errors and dead-time effects.
Instrumentation
for Gamma
Imaging
4
4.1 X-RAY VERSUS GAMMA-
RAY IMAGING
PROBLEM
Compare and contrast the principles of image formation between projection
X-ray and gamma scintigraphic imaging.
Solution
Figure 4.1 illustrates the differences between projection X-ray and gamma
scintigraphic imaging. In X-ray imaging, the radiation beams are emitted from
a single source (X-ray tube) in a fan or cone shape. The X-ray beams are pro-
jected in one direction only, that is, from the X-ray tube to the detector. Each
radiation beam will be attenuated by the tissues as it passes through the body
and reaches the detector. The final image formed on the detector is therefore
representing the attenuation properties of the tissues/organs.
53
54 Problems and Solutions in Medical Physics
Solution
a. A: Scintillating detector or sodium iodide doped with thallium
(NaI[Tl]) detector.
B: X-, Y-positioning circuit.
C: Energy discriminator circuit or pulse height analyser (PHA).
b. Energy discriminator circuit (known as pulse height analyser) is
used to filter the scattered photons that contain energies outside of
the pre-selected energy window(s). Firstly, the four output signals
(±X and ±Y) from the X-, Y-positioning circuit are weighted and
summed to form a voltage pulse (known as the ‘Z’ pulse), which
represents the intensity of a scintillation. Secondly, the Z-pulse is
then sent to the energy discriminator circuit and analysed if it falls
within the pre-set energy window(s) (with upper and lower energy
levels). In theory, scattered or coincidence photons can be rejected
this way because they fall outside the energy window. However, the
filtration is not perfect since the energy loss due to Compton scat-
tering is not very large, and the energy window cannot be too strict
as other factors, such as time and sensitivity need to be taken into
account.
Solution
When 140 keV gamma rays strike on a NaI(Tl) crystal, the gamma rays will
interact with the crystal via coherent scattering, photoelectric absorption or
Compton scattering, whereby the NaI molecules are raised to higher energy
states through ionisation and excitation. The ionised or excited atoms then
return to the ground state by emitting light photons. The light output is propor-
tional to the energy of the incident radiation. Approximately 38 light photons
are produced per 1 keV of deposited energy. The light photons then enter the
photomultiplier tubes (PMTs) which are attached to the NaI(Tl) crystal.
Solution
NaI(Tl) has relatively high density (3.67 g∙cm−3) and high stopping power
coefficient due to the high atomic number of iodine (Z = 53). Pure NaI does
not produce any scintillation after interacting with gamma ray at room tem-
perature. A trace amount (0.1%–0.4%) of thallium is therefore added as an
activator. The activating material also influences the wavelength (colour) of
the light produced by the scintillator. In NaI(Tl), the scintillation light is blue
(λmax = 415 nm). NaI(Tl) crystals are hygroscopic and fragile. Once the crystal
4 • Instrumentation for Gamma Imaging 57
absorbs water, it would cause a colour change that reduces light transmission
to the photomultiplier tubes (PMTs). Hence, the crystal is sealed in an alu-
minium container in a gamma camera system with a Pyrex glass backing to
allow the light to be transmitted to the PMTs.
Solution
In a scintillation-based gamma camera, incident gamma rays deposit their
energy in the scintillator where it is converted into visible (or near-UV) light
photons, which are detected by the photomultiplier tubes. In a solid-state
gamma camera (e.g. cadmium zinc telluride [CZT] detector), radiation depos-
its energy in the crystal lattice where it results in the generation of pairs of
charge carriers (direct conversion). The application of an electric field causes
the charge carriers to be swept to the cathode and anode of the device where
they induce a current pulse that can be detected and positional information
collected. The collimators in these systems are registered to the detector array
to maximise efficiency.
Solution
The gamma emissions from the patient’s body are isotropic; hence, a collima-
tor is necessary to form a useful image by permitting gamma rays approaching
the camera from certain directions to reach the detector while absorbing most
of the scattered photons. Without a collimator, the acquisition of all the emit-
ted radiation would result in an almost incoherent flood image.
58 Problems and Solutions in Medical Physics
Note: Collimators consist of a series of holes in a lead plate that allows the
gamma rays to pass through to the crystal. In this way, it behaves a little like
a lens on a camera.
Solution
• Sensitivity—Longer, narrower holes reduce the photon transmis-
sion, hence reduce sensitivity of the system.
• Resolution—Longer, narrower holes reduce the angle of accep-
tance, hence increase resolution of the image.
• Energy of the radioisotope—A higher energy source needs thicker
septa to block the scattered radiation again resulting in loss of
sensitivity.
Note: Collimators are designed for specific gamma energies. Usually below
180 keV are considered low energies (e.g. Tc-99m), between 180 and 370 keV
are considered medium energies (e.g. In-111), and 370–550 keV are considered
high energies (e.g. I-131). In each case, a high-resolution or high-sensitivity
collimator can be used depending on the clinical application required.
Solution
• Parallel-hole collimator: It consists of a lead plate with all of the
holes parallel to each other. The function is to allow only photons
4 • Instrumentation for Gamma Imaging 59
When the distance between the object and the pinhole is the same
as the distance between the pinhole and the crystal, there is no mag-
nification. If the object is moved yet farther from the pinhole, the
object will be minimised. The pinhole collimator is most commonly
used in thyroid imaging and paediatric nuclear medicine.
Solution
a. The collimator of a gamma camera is typically designed for less
than 5% photon transmission. The thickness of the septa is deter-
mined according to the purpose of the study; for example thicker or
longer septa are used for higher spatial resolution image. However,
a minimum thickness of the septa is required to form a useful image
which can be calculated using the following equation:
6d
µ
t≥
3
L−
µ
where:
t = thickness of the septa
d = diameter of the hole
L = length of the septa
4 • Instrumentation for Gamma Imaging 61
Solution
Therefore,
I (50 cm − 5 cm )
= = 0.75
O (50 cm +10 cm )
I
= 0.75
30 cm
I = 22.5 cm
Therefore, the image size is 22.5 cm.
Note: A diverging collimator contains holes that diverge from the detector
face. The holes diverge from a point (divergence point) typically 40–50 cm
behind the collimator, projecting a minified, non-inverted image of the source
onto the detector.
Solution
Magnification factor,
I
=
( f +t)
O ( f + t − b)
where:
I and O = image and object size
f = distance from the front of the collimator to the convergence point
t = thickness of the collimator
b = front of the collimator to the source
Therefore,
I ( 40 cm + 5 cm )
= = 1.5
O ( 40 cm + 5 cm − 15 cm )
4 • Instrumentation for Gamma Imaging 63
I
= 1.5
10 cm
Therefore, I = 15 cm
Note: A converging collimator contains holes that converge to a point at about
40–50 cm in front of the collimator. A converging collimator projects a magni-
fied, non-inverted image of the source onto the detector.
Solution
• A photomultiplier tube (PMT) consists of an evacuated glass tube
containing a photocathode, typically 10–12 electrodes (known as
dynodes) and an anode (see Figure 4.4).
• The PMT is fixed on to the NaI(TI) crystal with the photocathode
facing the crystal.
• The photocathode is usually an alloy of cesium and antimony that
releases electrons after absorption of light photons.
• A high-voltage power supply of approximately 1000 V is applied
between the photocathode and anode, with a series of resistors
dividing the voltage into equal increments.
• When a light photon strikes the photocathode, photoelectrons are
emitted via the photoemission process. The electrons are attracted
to the first dynode and are accelerated to the kinetic energies equal
to the potential difference between the photocathode and the first
dynode. For example, if the potential difference is 100 V, the kinetic
energy of each electron is 100 eV.
• Approximately five electrons are ejected for each electron that
strikes on the dynode. These electrons are then attracted to the next
dynode, reaching kinetic energy equal to the potential difference
between the first and second dynodes, and causing five electrons to
be ejected for each electron that strikes on it.
64 Problems and Solutions in Medical Physics
Solution
Scintillation events produced in the NaI(Tl) crystal are detected by a large
number of photomultiplier tubes (PMTs) which are arranged in a 2D array.
There are typically 30–90 PMTs in a modern gamma camera. The output
4 • Instrumentation for Gamma Imaging 65
voltages generated by the PMTs are fed to a position circuit which produces
four output signals, namely X+, X−, Y+ and Y−. These position signals con-
tain information about where the scintillations were produced within the
crystal, as well as the intensity of each scintillation. The intensity infor-
mation can be derived from the position signals by summing up the four
position signals (X+, X−, Y+ and Y−) to generate a voltage pulse (known as
the Z-pulse) which represents the intensity of a scintillation. The Z-pulse
is then sent to the pulse height analyser (PHA) to filter the scattered and
coincident photons.
Solution
Although a collimator can block >99% of the scattered radiation, there may be
conditions that some of the scattered photons escape the septa and reach the
detector. These ‘false’ signals are illustrated in Figure 4.5.
The gamma photons that reach the detector after passing through the col-
limator may undergo three possible interactions: (1) some may be absorbed
in the NaI(Tl) crystal, (2) some scatter from the crystal and (3) some pass
through the crystal without any interaction. The energy discrimination cir-
cuit filters the scattered or coincidence photons by summing the signals from
all the photomultiplier tubes (PMTs) (Z-pulse) and check if it falls within
the pre-set energy window. In principle, the Compton scattered photons will
FIGURE 4.5 Illustration of ‘false signals’ that reach the NaI(Tl) detector.
66 Problems and Solutions in Medical Physics
lose energy, and the coincidence interactions will contribute extra energy.
A pulse height analyser (PHA) is used to filter the signals that fall outside
of the preselected energy window or ‘channels’; hence, only the signals that
fall within the selected energy will be registered as a ‘true’ and count in the
image.
Solution
a. In frame mode, a matrix with size approximate to the area of the
detector is defined. A position (X, Y) in the detector corresponds to
a pixel position in the matrix. Before acquisition, all pixels within
the matrix are set to zero. During acquisition, every time when a
new signal is detected, it is added to the corresponding (X, Y) pixel
in the matrix. The data acquisition continues until a preselected
time or total count is reached. Therefore, in frame mode, the size
and depth of the matrix, number of frames per study and time per
frame or total counts must be specified.
In list mode, the (X, Y) signals are stored in a list instead of
being immediately formed into an image. Each signal is coded with
‘time mark’ according to the time it is detected and stored as indi-
vidual event. After acquisition is completed, the data can be sorted
and displayed to form the desired images. Data can be manipulated
by changing matrix size, time per frame, physiological markers, etc.
after the acquisition. The ‘bad’ signals can also be discarded manu-
ally by the user.
b. The advantages and disadvantages of frame and list modes are sum-
marised below:
4 • Instrumentation for Gamma Imaging 67
ADVANTAGES DISADVANTAGES
Frame The data are immediately The matrix (image) format
Mode displayed in a designated must be specified before
image format; hence, acquisition, and no
limited post-processing manipulation of the data is
work is required. It requires possible after the image is
lesser memory and disk formed.
space for acquisition and
storage compared to the
list mode.
List It provides flexibility to It generates a large amount of
Mode manipulating data data which requires larger
according to the way the memory for acquisition,
users want it to be greater processing power and
displayed. disk space for storage.
Solution
a. A profile is obtained from two small radioactive sources placed at
a known separation, or from a single source moved to known posi-
tions on the X or Y axis of the camera. The profile will give two
spatial peaks. The distance between the maxima of the two peaks
is measured in terms of pixels and used to calculate the distance
between the two source positions. Knowing the distance and the
number of pixels allows calculation of the individual pixel size.
40 cm
=
b. Pixelsize = 0.3 cm per pixel
128 pixel
68 Problems and Solutions in Medical Physics
Solution
For 128 × 128 matrix size:
100%
Statisical =
noise = 18%
30.5
100%
=
Statisical noise = 36%
7.6
Hence, the matrix size of 256 × 256 has 50% more noise than the matrix size
of 128 × 128.
To achieve the same noise level in 256 × 256 and 128 × 128 matrix sizes:
Total counts
= 30.5 counts per pixel
256 × 256 pixels
Solution
In static study, a single image is acquired for either a pre-set time interval or until
the total numbers of counts reaches a pre-set number. The data is acquired in frame
mode, and matrix size is predefined before acquisition. Static planar imaging is
used for studies in which the distribution of the radiopharmaceutical is effectively
static throughout the acquisition. It provides information, such as organ morphol-
ogy (size, shape and position), and regions of increased or decreased uptake.
Clinical imaging examples of static imaging studies include Tc-99m-
HDP bone scan, I-123-iodide thyroid scintigraphy and renal imaging using
Tc-99m-DMSA.
Solution
In a dynamic study, a series of images are acquired in sequence one after
another, for a pre-set time per image. This is usually done when the biokinetics
of the tracer is relatively rapid. The frame rate (frame per second) and matrix
size can be varied between images; however, the patient’s position cannot be
changed during the whole course of study.
The clinical imaging examples of dynamic studies include oesophageal
transit, hepatobiliary imaging and renogram.
70 Problems and Solutions in Medical Physics
Solution
A gated study is performed when the dynamic process of an organ occurs too
rapidly to be effectively captured by a dynamic image acquisition, and when
the dynamic process is repetitive, such as the cardiac or respiratory cycle.
Gated acquisitions require a physiological monitor that provides a trigger pulse
to mark the beginning of each cycle of the process being studied. For example,
in a gated cardiac study, an electrocardiogram (ECG) monitor provides a trig-
ger pulse to the computer whenever the monitor detects a QRS complex. Using
this technique, the image data from each desired phase of the cycle can be
stored in specific location (known as ‘bins’). When all the data in a bin are
added together, the image represents the specific phase of the cycle.
Examples of gated studies include cardiac blood pool study for calculation
of left ventricular ejection fraction, gated SPECT myocardial perfusion study
for regional wall uptake and motion and a respiratory gated study of the lungs.
SPECT and
PET Imaging 5
5.1 PHYSICAL PRINCIPLES OF SPECT
PROBLEM
Explain briefly the basic principles of single-photon emission computed
tomography (SPECT).
Solution
SPECT is performed using a gamma camera that consists of one or more
detectors. Each detector is comprised of a collimator, scintillation crystal,
array of photomultiplier tubes (PMTs) and digital positional network. During
SPECT imaging, the detectors are rotated around the patient at small-angle
increments, typically 3°–10°. The rotation can be in continuous motion or
using the ‘step-and-shoot’ method.
In most cases, a full 360° rotation is used to obtain an optimal reconstruc-
tion; however, in cardiac imaging, a rotation of 180° is sufficient. Individual
planar images are acquired at each angle for a set period of time, for example
30 seconds. In each detector, gamma rays pass through the collimator and
are absorbed by the crystal. The gamma ray energies are converted into light
photons and detected by the PMTs. The images are displayed in the form of a
digital matrix, for example, a 128 × 128 or 256 × 256 matrix. Each pixel con-
tains position and count rate information. The two-dimensional (2D) planar
images are pre-filtered and reconstructed using either an iterative process or
filtered back projection to form a three-dimensional (3D) image. The data may
then be manipulated to show thin slices along any chosen axis (e.g. transverse,
sagittal or coronal) of the body. The count rates of each voxel are commonly
displayed using a colour scale.
71
72 Problems and Solutions in Medical Physics
Solution
In planar imaging, the radioactivity in the tissues in front and behind of an
organ reduces contrast. If the activity in these overlapping structures is not
uniform, the pattern of activity distribution from different structures will be
superimposed. This causes a source of structural noise that hinders the ability
to distinguish the activity distribution in the tissues. SPECT imaging can
greatly improve the subject contrast and reduce structural noise by eliminating
the activity in overlapping structures. In addition, SPECT imaging also allows
partial correction of tissue attenuation and photon scattering in the patient.
The possible disadvantage of SPECT includes the slight decrease of spatial
resolution compared to planar imaging. This is because the detector is usually
closer to the patient in planar imaging than in SPECT; however, the increase
in image information at depth easily compensates for any loss of resolution.
In addition, using a shorter time per angle in SPECT may necessitate the use of
a lower resolution/higher sensitivity collimator to obtain an adequate number
of counts within the time frame.
Solution
a. A collimator is required to project the gamma photons onto the
detector crystal. Collimators are essentially constructed as lead
5 • SPECT and PET Imaging 73
plates with holes or lead foils to channel the photons. Any photon
not running parallel to the hole will be absorbed by the septa. For
diagnostic imaging, collimators are classified according to three
energy ranges:
EXAMPLES OF
ENERGY LEVEL ENERGY RANGE (KeV) RADIOISOTOPES
Low energy 100–200 Tc-99m (140 keV)
I-123 (159 keV)
Medium energy 200–300 In-111 (171 and
245 keV)
High energy 300–400 I-131 (364 keV)
The higher the energy, the thicker the lead septa are required to
attenuate the photon. Hence, the overall sensitivity is reduced.
In addition, the spatial resolution is also reduced at higher energy.
b. For SPECT imaging, the camera rotates around the patient. The
peripheral margins of the patient are closer to the detector, and hence
the photons originating from that area will be subjected to less atten-
uation than those arising from the centre of the patient. By defining
the size of the patient from the outline of the image (an X-ray com-
puted tomography [CT] image may be used), attenuation correction
could be applied to correct for photon attenuation with depth.
N = N oe − µ t
where:
N = number of transmitted photons
No = number of incident photons
μ = linear attenuation coefficient of the tissue
t = thickness of the tissue
The attenuation effect will be more pronounced for lower energy
gamma rays, that is, Tc-99m (140 keV). The resulting image data will
therefore have a lower count density in the centre of the patient in
comparison to the periphery. This can be corrected by modifying the
image data using the appropriate attenuation coefficient for the energy
of the gamma rays used (i.e. Tc-99m = 0.12). In this way, counts are
added to the central regions to give a corrected uniform image.
74 Problems and Solutions in Medical Physics
Solution
The Fourier noise-filtering method is done in the spatial domain. It is used
to eliminate the blurring function (star artefact) from simple back projection.
In this method, the amplitude of different frequencies is modulated. The princi-
ple is to preserve the broad structures (the image) represented by low frequency
and remove the fine structure (noise) represented by high frequency. There are a
number of filters classified according to their functions. For example, a ‘low pass
filter’ (e.g. Butterworth, Parzen, and Weiner filters) rejects the high-frequency
data, and a ‘high pass filter’ (e.g. ramp filter) rejects the low to medium fre-
quency data. Often the low pass filters are applied as pre-filters to remove high-
frequency noise, and the ramp filter is then applied during back projection to
remove star artefact. This combination of filter is called ‘high pass-low pass
filter’, for example, Ramp-Parzen filter. All of the Fourier filters are character-
ised by a maximum frequency, known as the Nyquist frequency. The Nyquist
frequency is calculated as 0.5 cycle per pixel, which is the highest fundamental
frequency useful for imaging. The cut-off frequency is the maximum frequency
the filter will pass. If the cut-off frequency is higher than the Nyquist frequency,
the data will be filtered at the Nyquist frequency.
Solution
The convolution method is carried out in the spatial domain. It uses a smooth-
ing ‘kernel’, which is basically a mathematical function to remove the 1/r
5 • SPECT and PET Imaging 75
In the nine-point smoothing method, the individual pixel value in the acquisi-
tion matrix is multiplied with a 3 × 3 smoothing matrix. The sum of the product
is then divided by the sum of all pixel values of the smoothing matrix. In this
example, the pixel value ‘7’ is converted to ‘5’ after smoothing. Similarly, all
pixel values in the acquisition matrix are smoothed, and a smoothed image is
produced.
Solution
• In the iterative method, an initial estimate activity is made of each ele-
ment in the N × N matrix, and values are compared to the measured
values.
• Corrections are applied according to the estimated values.
76 Problems and Solutions in Medical Physics
Solution
Figure 5.1 illustrates the basic principle of positron emission tomography.
In PET imaging, a positron emitting radiopharmaceutical, such as F-18-FDG,
is administered into the patient’s body. The positrons lose most of their energy
in matter by causing ionisation and excitation. When a positron has lost most
5.8 ANNIHILATION
COINCIDENCE DETECTION
PROBLEM
With the aid of a diagram, discuss briefly the annihilation coincidence detec-
tion (ACD) method used in PET imaging.
Solution
Figure 5.2 illustrates the principles of coincident annihilation detection method
used in PET imaging. When a positron is emitted by nuclear transformation,
it transverses through matter and loses energy. After it loses most of its energy,
it interacts with an electron, resulting in two 511 keV photons that are emitted
in nearly opposite directions. This phenomenon is known as annihilation.
When the two photons are simultaneously detected within a ring of detec-
tors surrounding the patient, it is assumed that the annihilation event occurred
on the line connecting the interactions (i.e. the line-of-response [LOR]).
78 Problems and Solutions in Medical Physics
Each detector will generate a timed pulse when it detects an incident photon.
These pulses are then combined in coincidence circuitry, and if the pulses fall
within a short time window (within approximately 5ns), they are deemed to be
coincident and the signal will be registered in the image.
Solution
a. True coincidence: the two 511 keV photons that are detected from
opposite detectors are from the actual source location (Figure 5.3).
b. Scatter coincidence: the photons are scattered in the patient’s body
before being detected by the PET detectors. The annihilation coinci-
dence detection (ACD) method does not present the actual location
of the source.
5 • SPECT and PET Imaging 79
Solution
A PET scanner measures time-of-flight under the assumption that the location
of the annihilation can be localised along the line of flight of the coincident
photons by measuring time of arrival of each photon at the opposing crystal.
Unless the event occurs in the exact centre of the ring, one of the photons will
arrive before the other. Typical coincidence time is 12 ns. The time difference is
proportional to the difference in distance and can be used to calculate the posi-
tion of the event.
The disadvantages of this method include lower accuracy in determining
the location of the source due to possibilities of scatter events and random
coincidence, as well as reducing spatial resolution of the image.
80 Problems and Solutions in Medical Physics
Solution
The resolution of the PET image is determined mainly by the detector size.
Smaller detectors give better resolution. The ring diameter and the detector
material also affect the intrinsic resolution of the PET image. The energy of
the positron plays a minor role as higher energy positrons travel farther away
from the site of origin before annihilating.
Solution
In 2D acquisition, coincidence events are collected by detectors in a single ring
around the source, and each detector is in coincidence with any other detector
in the ring and defines a 2D plane. In 3D acquisition, coincidence events are
collected by detectors in multiple rings or in 2D arrays of detectors around
the source. Each detector is in coincidence with any other detectors in all the
rings or in the entire opposing 2D array and helps to define a 3D image. 2D
acquisition is less prone to scatter coincidences and requires less electronics
for volumetric acquisition. 3D acquisition is more prone to false coincidence
effects and requires fast processing electronics.
5 • SPECT and PET Imaging 81
Solution
Comparison of physical properties of Tc-99m and F-18:
Tc-99m F-18
Production method Mo-99 generator Cyclotron
Mode of decay Isomeric transition (88%) Positron emission (97%)
Internal conversion (12%) Electron capture (3%)
Decay product Tc-99 O-18 (stable)
Physical half-life 6.02 hours 110 min
Detection photon 140 keV 511 keV
energy
Detection principle Single photon detection Dual photon (511 keV each)
travelling in opposite
directions
c. Radionuclide
d. Spatial resolution
e. Attenuation
Solution
SPECT PET
a. Principle of Using physical collimator Annihilation coincidence
projection data detection
collection
b. Transverse image Filtered back projection or Filtered back projection or
reconstruction iterative reconstruction iterative reconstruction
c. Radionuclide Any radionuclides emitting Positron emitters only
X-rays, gamma rays or
annihilation photons.
d. Spatial resolution Depends on collimator and Relatively constant across
camera orbit. Within a transaxial image, best at
transaxial image, the the centre. Typically
resolution in the radial 4.5–5 mm full width at
direction is relatively half maximum (FWHM)
uniform, but the at the centre.
tangential resolution is
degraded towards the
centre. Typically, 10 mm
full width at half
maximum (FWHM) at the
centre for a 30-cm
diameter orbit using
Tc-99m. Larger camera
orbits produce lower
spatial resolution.
e. Attenuation Attenuation is less Attenuation is more
significant. Attenuation significant. Accurate
correction methods are attenuation correction is
available. possible with a
transmission source.
5 • SPECT and PET Imaging 83
Solution
In SPECT, the spatial resolution and the detection efficiency are primarily
determined by the collimator. Both are ultimately limited by the compromise
between collimator efficiency and collimator spatial resolution that is a con-
sequence of collimated image formation. Larger camera orbits result in lower
resolution. This causes the spatial resolution to deteriorate from the edge to the
centre in transverse SPECT images.
In comparison, PET has better spatial resolution and detection efficiency.
ACD method is used as an electronic collimator for the PET scanner by deter-
mining the path of the detected photons. ACD allows detection of two annihila-
tion photons given off by a single positron annihilation event. The two photons
(511 keV) are considered to be from the same event if they are counted within a
defined coincidence timing window and recorded as occurring along the line-
of-response (LOR). The location of the event can then be calculated using the
principle of time-of-flight (TOF). This combination approach effectively filters
the scattered photons and determines the location of the source, hence improv-
ing the spatial resolution and detection efficiency of PET imaging.
Imaging
Techniques
in Nuclear
6
Medicine
Solution
Principle
• Whole body bone SPECT imaging is a sensitive method for detect-
ing a variety of anatomical and physiological abnormalities of the
musculoskeletal system. A Tc-99m-labelled radiopharmaceutical
(such as methylene diphosphonate [MDP]), which localises in the
bone, is injected intravenously to the patient. The activity accu-
mulated in the bone is then evaluated through the SPECT images.
It helps to diagnose a number of bone conditions, including bone
cancer or metastasis, inflammation, fractures and bone infection.
Radiopharmaceutical
• The most commonly used radiopharmaceutical for bone SPECT
imaging is Tc-99m labelled with methylene diphosphonate (MDP).
85
86 Problems and Solutions in Medical Physics
Image Acquisition
• SPECT imaging is performed at two to three hours following injec-
tion of the radiopharmaceutical to allow optimal concentration of
the radiotracer uptake by active osteoblasts. The delay also allows
time for background activity to be excreted from the body. Patients
are required to empty their bladders before the scan.
• A low-energy high-resolution (LEHR) collimator is used. The
acquisition energy window is set to 140% ± 10% keV and a matrix
size of 256 × 1024 is typically used. The patient is positioned in the
supine, arms by the side and feet first position. Anterior and pos-
terior views are acquired simultaneously if a dual-head camera is
used. The total imaging time is approximately 30 minutes. Patients
are asked to stay still during the scanning.
Solution
Multiple-gated acquisition (MUGA) or gated cardiac blood pool imaging using
Tc-99m-labelled red blood cells (Figure 6.1).
FIGURE 6.1 Basic principles of multiple gated acquisitions (MUGA) for cardiac
ejection fraction study. The numbers 1 to 8 represent the frames allocated over
the cardiac cycle/heart beat.
88 Problems and Solutions in Medical Physics
6.3 RENOGRAM
PROBLEM
An adult male patient is referred for renal scintigraphy or renogram to assess
his kidney function.
a. Describe briefly the procedure of this examination including patient
preparation, radiopharmaceutical, administered activity and imag-
ing technique.
b. Explain briefly the method of quantitative analysis of the renogram
to assess the kidney function.
c. Name two examples of clinical information that can be derived
from a renogram.
Solution
a. A renogram refers to serial imaging after intravenous administra-
tion of Tc-99m-DTPA or Tc-99m-MAG3. It is used for qualitative
6 • Imaging Techniques in Nuclear Medicine 89
Solution
Radioimmunoassay (RIA) is a sensitive in vitro assay technique used to mea-
sure the concentration of an antigen (e.g. hormone level in the blood) using a
specific antibody. During the test, a known amount of antigen is labelled with a
radioactive source (e.g. I-125), and the radiolabelled antigen is then mixed with
a known amount of antibody. The radiolabelled antigen will bind together with
the antibody. Subsequently, a sample of serum from the patient containing
an unknown amount of the same antigen is added. This causes a competitive
binding of two antigens with the antibody. The antigen, which has a higher
concentration, will bind extensively with the antibody, displacing the other.
In this case, if the antigen in the patient’s serum has a higher concentra-
tion, it would bind with the antibody and displace the radiolabelled antigen.
The bound antigens are then separated from the unbound ones, and the
radioactivity of the unbound antigen remaining in the supernatant is measured
using a gamma counter. Using known standards and formulas, the amount of
antigen in the patient’s serum can be derived.
Note: A number of non-radioactive methods can now be used to replace RIA,
such as the enzyme-linked immunosorbent assay (ELISA).
6.5 STANDARDISED
UPTAKE VALUE (SUV)
PROBLEM
Describe the use of standardised uptake value (SUV) in positron emission
tomography (PET) and its limitations.
Solution
The SUV is used as a relative measure of radiopharmaceutical uptake, such
as FDG, standardised for the injected activity and patient weight. An ROI is
6 • Imaging Techniques in Nuclear Medicine 91
defined on the area of the PET image, and the SUV within the ROI is calcu-
lated using the following formula:
The main purposes of SUVs are to provide a good correlation between the his-
tological grades of the tumour and to facilitate comparisons between patients.
However, there are several limitations of the use of an SUV. First, it has a
large degree of variability due to physical and biological sources of errors.
Furthermore, the suggestion of a threshold value of an SUV to characterise
benign or malignancy of a lesion remains a clinical challenge. The sensitiv-
ity and specificity of using the SUV threshold usually decreases with lesions
smaller than 7 mm. Another limitation of the SUV is in identifying the lesion
boundaries and determining the counts within the ROI. The SUV may be
influenced by image noise, low image resolution and variable user-biased ROI
selection.
Solution
In oncology, F-18-FDG PET imaging is very useful for the following
applications:
• To assess regional tumour metabolic activity by evaluating the dis-
tribution and uptake as standard uptake values (tumour staging).
• To monitor the success of therapy by comparing the pre- and post-
therapy images.
• To detect early recurrent tumours.
• To provide a whole-body survey for cancer that may have metasta-
sised, as the FDG radiopharmaceutical is distributed throughout the
body and whole-body imaging is commonly done without adding
additional activity or dose to the patient.
• To identify benign and malignant growths.
• To guide radiotherapy dose planning.
92 Problems and Solutions in Medical Physics
Solution
Unsealed source therapy relates the use of soluble forms of radioactive sources,
which are administered to the patient body by either injection or ingestion.
Unsealed source therapy is commonly known as radionuclide therapy or
molecular targeted radiotherapy. Sealed source therapy uses radionuclides,
which are sealed in a capsule, rod or metal wire during therapy. Sealed source
therapy is more commonly known as ‘brachytherapy’.
The main factors to be considered when choosing a suitable radionuclide
for unsealed source therapy include its effective half-life, types of radiation,
energy range and its biochemical characteristics. Effective half-life includes
both the physical and biological half-life. A suitable range of physical half-life
is between six hours to seven days, whereas the biological half-life depends on
a number of factors, such as radiotracer delivery, uptake kinetics, metabolism,
clearance and excretion. For therapeutic purposes, radiations with high linear
energy transfer (LET), such as alpha and beta particles, are preferred because
they allow very high ionisation per length of travel and reduce the need for
additional radiological protection. For biochemical characteristics, a suitable
radionuclide should provide selective concentrations and prolonged retention
in the tumour, while maintaining minimum uptake in the normal tissue.
93
94 Problems and Solutions in Medical Physics
Solution
i. Hyperthyroidism treatment using I-131-sodium iodide radiopharma-
ceutical.
ii. Thyroid carcinoma treatment using high-dose I-131-sodium iodide.
iii. Treatment of polycythemia using P-32-phosphate.
iv. Bone metastases palliative treatment using Sm-153-EDTMP
(methylenephosphonic acid) or Ra-223-dichloride.
v. Treatment for hepatocellular carcinoma using Y-90-labelled micro-
spheres or Re-188-lipiodol.
Solution
T p × Tb
• Effective half-life, Te = = 2.18 days
T p + Tb
• Cumulated activity, A = 1.44 A0·Te·fh = 1.44 × 925 MBq ×
2.18 days × 0.85
= 2468.196 MBq·days × 24 h·day−1
= 59236.70 MBq·h
• Absorbed dose = A·S = 59236.70 MBq·h × 2.2 × 10 −3 mGy(MBq·h)−1
= 130.32 mGy
7 • Radionuclide Therapy 95
Solution
Before the Treatment
• The patient will be given a list of prohibited food, dietary informa-
tion and medication that could potentially affect the thyroid uptake
of the administered radiopharmaceutical.
• Female patients who are of childbearing age must be tested for preg-
nancy prior to administration of the I-131, unless pregnancy is clini-
cally proven to be impossible in the patient.
• Patients who are breastfeeding should be excluded or feeding should
be carried out by someone else using bottled formula milk.
• The treating nuclear medicine physician should confirm that the
patient is continent of urine or an arrangement should be made to
prevent any contamination caused by urinary incontinence.
• Written informed consent must be obtained from the patient before
therapy.
Post-Treatment
• Before discharge from the hospital, the patient needs to be briefed
on the methods to reduce unnecessary radiation exposure to the
family members and members of the public. Written instructions
should be provided to the patient.
• The patient should be advised to avoid pregnancy for a period of
time stated in the local radiation protection guidelines.
• Patients must be provided with a written document stating that they
have been given a radionuclide treatment and details of the radionu-
clide since the radioactivity may be detected by monitoring devices
at borders such as airports during travel.
Solution
Using the Marinelli-Quimby formula,
21.4 AW
Administered activity (MBq) =
UTe
where:
A = prescribed dose (Gy)
W = weight of the thyroid glands (g)
U = percentage uptake of I-131 at 24 h (%)
Te = effective half-life (d)
7 • Radionuclide Therapy 97
Therefore,
π abc
Volume ( mL) =
6
where:
a = length (cm), b = width (cm) and c = thickness (cm).
The total thyroid volume is obtained by adding the volume of both lobes.
Solution
Radioiodine can easily cross the maternal placenta, and foetal thyroid uptake
following therapeutic doses can pose significant problems for the foetus, par-
ticularly with respect to causing permanent hypothyroidism. As a general rule,
pregnant women should not be treated with any radioactive substance unless
the therapy is required to save her life. Thyroid cancers are relatively unag-
gressive compared to most other cancers. As a result, radioiodine treatment
should be delayed until after pregnancy. However, if any therapy is to be per-
formed during pregnancy, this should be restricted to after thyroid surgery
during the second or third trimester of pregnancy.
98 Problems and Solutions in Medical Physics
Note: According to the ICRP Publication 84, a major problem occurs when
a female, who is not thought to be pregnant, is treated for thyroid carcinoma
and is found to be pregnant after the administration of radioiodine. If a patient
is discovered to be pregnant shortly after a therapeutic radioiodine admin-
istration, maternal hydration and frequent voiding should be encouraged to
help eliminate maternal radioactivity and to reduce radioiodine residence time
in the bladder. If the conceptus is more than 8 weeks’ post-conception and
the pregnancy is discovered within 12 hours of radioiodine administration,
giving the mother 60–130 mg of stable potassium iodide will partially block
the foetal thyroid and reduce thyroid dose. After 12 hours post-radioiodine
administration, this intervention is not very effective.
Solution
The safe oral administration of I-131 should be practiced as following:
Solution
The guidance level as recommended by the IAEA Safety Report Series No. 63
is illustrated in Table 7.1.
TABLE 7.1 IAEA recommended guidance levels for I-131 therapy patients
7.9 RADIOEMBOLISATION
PROBLEM
Explain the principle, techniques and advantages of radioembolisation in
cancer treatment.
Solution
Radioembolisation is a cancer treatment in which radioactive particles of specific
sizes (about 15–60 μm) are delivered to the tumour through intra-arterial injection.
The radioactive particles will be trapped in the microvessels and reside within the
tumour volume to emit radiation that kill the cancer cells. Radioembolisation is
most often used to treat unresectable or metastasised liver cancer.
100 Problems and Solutions in Medical Physics
Solution
RIT is a type of cancer-cell-targeted therapy that uses monoclonal antibodies
(mAbs) directed against tumour-associated antigens labelled with a therapeu-
tic radionuclide. It is a combination of radiation therapy and immunotherapy.
When injected into the blood stream, the radiolabelled mAbs will travel to and
bind to the cancer cells, allowing a high radiation dose to be delivered directly
to the tumour. The major benefit of RIT is the ability for the antibody to spe-
cifically bind to a tumour-associated antigen, hence increasing the radiation
dose delivered to the tumour cells while decreasing the radiation dose to the
surrounding normal tissues.
RIT can be used to treat non-Hodgkin’s B-cell lymphoma, chronic
lymphocytic leukaemia and immune diseases.
Internal
Radiation
Dosimetry
8
8.1 INTERNAL RADIATION DOSIMETRY
PROBLEM
Explain why the estimation of absorbed dose in internal radiation dosim-
etry presents greater challenges than the measurement of external radiation
exposure.
Solution
The estimation of absorbed dose in internal radiation dosimetry is more com-
plicated due to the following factors:
101
102 Problems and Solutions in Medical Physics
Solution
Factors that are likely to influence the absorbed dose to an organ include:
• Type of radionuclide.
• Activity administered.
• Activity in the source organ.
• Activity in the other organs.
• The size and shape of the organs.
• The kinetics of the radiopharmaceutical.
• The quality of the radiopharmaceutical in terms of the radiochemi-
cal purity.
Solution
Source organ(s) is/are the organ(s) of interest with significant uptake of
radioisotope.
Target organ(s) is/are the organ(s) that receives radiation from the other
source organs, for which the absorbed dose is to be calculated. The source and
target organs may be the same when the radiation dose, due to the radioactivity
in the target itself, is calculated.
8 • Internal Radiation Dosimetry 103
Solution
Monte Carlo modelling is a ray-tracing method in which the fate of an indi-
vidual photon or particle is determined. The method is based on randomly
sampling a probability distribution for each successive interaction. For exam-
ple, Monte Carlo modelling can be used to simulate the random trajectories of
individual photons or particles in a matter by using the knowledge of probabil-
ity distribution governing the individual interactions of photons or particles in
materials.
One keeps track of the physical quantities of interest for a large number
of histories to provide the required information about the average quantities.
Monte Carlo modelling in radiation dosimetry requires detailed knowledge of
the absorption and scattering coefficients for the specific radiation energies
and for the properties of various types of tissues.
D rk = ∑ A S(r ← r )
h
h k h
A h = A0 ⋅ fh ⋅1.44Te
Name all the quantities used in the formulas above and state their SI units.
104 Problems and Solutions in Medical Physics
Solution
SYMBOLS QUANTITIES SI UNITS
Solution
Absorbed fraction is the fraction of energy emitted by the source organ that is
absorbed in the target organ.
The value is determined by:
8.7 S-VALUE
PROBLEM
Define S-value used in the MIRD formula. What factors determine the
S-value?
Solution
S-value is defined as the dose to the target organ per unit of cumulated activity
in a specified source organ. It is determined by:
• The mass of the target organ.
• The type and amount of ionising radiation emitted per disintegration.
• The fraction of the emitted radiation energy that reaches and is
absorbed by the target organ.
• Each S factor is specific to a particular source organ/target organ
combination.
Solution
• Effective half-life, Te = (T p ×Tb ) / (T p + Tb )
= 2.67 days
• Cumulated activity, A = 1.44A0 f hT
= 1.44 × 1110 MBq × 2.67 × 0.6
= 2560.64 MBq days × 24 h(day)−1
= 61455 MBq ⋅ h
S = 61455 MBq ⋅ h × 2.2 × 10 −3 mGy(MBq ⋅ h) −1
• Absorbed dose = A⋅
= 135.2 mGy
106 Problems and Solutions in Medical Physics
Solution
A = 1.44 A0 f h Te
D = A S
= 4.2 × 10 −3 Gy
= 1.4 × 10 −4 Gy
8 • Internal Radiation Dosimetry 107
Solution
Physical half-life of Tc-99m = 6 h. The effective half-life of two biological
clearances:
Te1 = (3 × 6/3 + 6) = 2 h
Te2 = (7 × 6/7 + 6) = 3.2 h
à = 1.44 A0 f hTe
= 1.44 × 148 × 106 × 0.99 × (0.45 × 2 + 0.55 × 3.2)
= 0.56 GBq·h
D =÷S
= 0.56 × 109 × 1.40 × 10 −11 = 7.89 mGy
Solution
Assume there is no biological excretion of the radioisotope, hence the
Te = Tp = 6 h
Hence, the total absorbed dose to the liver is the sum of all D above
Solution
i. The radioactivity is assumed to be uniformly distributed in each
source organ and absorbing target organ.
ii. The organ size and geometries are simplified for mathematical
computational.
iii. The tissue density and composition are assumed to be homogenous
in each organ.
iv. The phantoms for the ‘reference’ adult and child are only approxi-
mations of the physical dimensions of any given individual.
v. The energy deposition is averaged over the entire mass of the target
organs.
vi. Dose contributions from Bremsstrahlung and other secondary
sources are ignored.
vii. With a few exceptions, low-energy photons and all particulate radia-
tions are assumed to be absorbed locally (i.e. non-penetrating).
Quality
Control in
Nuclear
9
Medicine
Solution
• Constancy test—daily.
• Accuracy test—at installation, annually, and after repairs.
• Linearity test—at installation, quarterly, and after repairs.
• Geometry test—at installation and after repairs.
111
112 Problems and Solutions in Medical Physics
Solution
Extrinsic measurement is the measurement of a scintillation camera perfor-
mance with collimator attached; while intrinsic measurement is the measure-
ment of a scintillation camera performance when the collimator is removed.
a. Uniformity
b. Spatial resolution
c. Energy resolution
d. Spatial linearity
e. System efficiency
Solution
a. Uniformity is a measure of a camera’s response to uniform irradia-
tion of the detector surface. The intrinsic uniformity is measured
by placing a point radioactive source (typically 3.7–18.5 MBq
Tc-99m) in front of an uncollimated camera. The source should
be placed at a distance more than four times the largest dimension
of the crystal, and at least five times away if the uniformity image
is to be analysed quantitatively. Extrinsic uniformity is assessed
by placing a uniform planar radioactive source in front of a col-
limated camera. The planar source should be large enough to cover
the crystal of the camera.
9 • Quality Control in Nuclear Medicine 113
Es = Ec × Ei × f
where:
Ec = fraction of photons emitted by a source that penetrate the
collimator holes.
Ei = fraction of photons penetrating the collimator that interact
with the detector.
f = fraction of interacting photons accepted by the energy
discrimination circuits.
9.4 UNIFORMITY
PROBLEM
Discuss the common causes of non-uniformity of a gamma camera image and
ways to overcome them.
Solution
The common causes of non-uniformity of a gamma camera image include:
Solution
i. The collimator hole’s size and length will affect the range of
accepted photon angles. Thicker or longer septa will block more
scattered photons, hence increasing the spatial resolution.
ii. Scintillation crystal thickness can affect the probability of photons
scattering. Although thicker crystal increases the sensitivity, it also
decreases the spatial resolution of the system.
iii. Gamma ray energy may affect the probability of scatter in crystal
and probability of septal penetration. Higher energies increase the
probability of Compton scattering and septal penetration, hence
reducing the spatial resolution.
116 Problems and Solutions in Medical Physics
iv. The size and number of PMTs may affect positional accuracy. Better
spatial resolution can be achieved by using more and smaller-size
PMTs.
v. The distance of the source-to-collimator can affect the spatial
resolution. More scattered radiation can pass through the collimator
holes when the source is farther away from the collimator, hence
reducing the spatial resolution.
9.6 SENSITIVITY
PROBLEM
a. Define the sensitivity of a gamma camera.
b. What are the factors affecting the sensitivity of a gamma camera?
Solution
a. The sensitivity of a gamma camera is defined as the number of
counts per unit time detected by the camera for each unit of activity
from a radioactive source (cps∙MBq−1).
b. The factors affecting the sensitivity of a gamma camera include the
geometric efficiency of the collimator, detection efficiency of the
detector, pulse height analyser (PHA) discrimination settings and
the dead time of the system.
Solution
a. Collimator efficiency is also known as geometric efficiency. It is
defined as the number of gamma ray photons passing through the
collimator holes per unit activity presents in a radioactive source.
9 • Quality Control in Nuclear Medicine 117
d4
Ec = K ⋅
t ( d + a )2
2
e
where:
K = constant, which is a function of the shape and arrangement
of holes in the collimator. It varies between 0.24 for round
holes in a hexagonal array to 0.28 for square holes in a
square array.
d = hole diameter
te = effective length of the collimator hole
a = septal thickness
b. The collimator efficiency, Ec, for parallel-hole collimator increases
with increasing diameter of the collimator holes (d) and decreases
with increasing collimator thickness (t) and septa thickness (a). For
parallel-hole collimators, Ec is not affected by the source-to-detector
distance for an extended planar source. However, for other types of
collimators, Ec varies with the types of collimators as a function of
source-to-collimator distance.
Solution
A—Low-resolution collimator
B—Medium-resolution collimator
C—High-resolution collimator
Solution
a. MTF of an imaging system is a plot of the imaging system’s modu-
lation versus spatial frequency, where modulation is essentially the
output contrast normalised by the input contrast. The principle of
MTF is illustrated in Figure 9.3.
9 • Quality Control in Nuclear Medicine 119
Amax − Amin
Ms =
Amax + Amin
the valley, which are smaller in magnitude than the Amax and Amin.
Therefore, the modulation in the image (MI ) is given as:
Cmax − Cmin
MI =
Cmax + Cmin
MI
MTF(v) =
MS
Solution
a. Multienergy spatial registration is a measure of the camera’s ability
to maintain the same image magnification, regardless of the ener-
gies deposited in the crystal by the incident photons.
b. The multienergy spatial registration can be assessed by imaging
several point sources of Ga-67. Ga-67 has three useful photopeak
energies for imaging at 93, 184 and 300 keV. The source is placed
offset from the centre of the camera. Only one photopeak energy
is tested at a time. The centroid of the count distribution of each
source should be at the same position in the image for all three
energies.
9 • Quality Control in Nuclear Medicine 121
Solution
i. Off-centred window setting on pulse height analyser
ii. Defective photomultiplier tube
iii. Metal in the field-of-view (metal-like jewellery)
iv. NaI crystal defect (such as a crack)
Solution
a. Ideally, the COR of a SPECT system must be aligned exactly to the
centre, in the X-direction of all projection images (shown in the left
image in Figure 9.4). If the COR is misaligned and not corrected,
it may cause a loss of spatial resolution in the resultant transverse
images. If the misalignment is large, it can cause a point source to
appear like a ‘doughnut’ shape (as illustrated in Figure 9.4).
b. The causes of COR misalignments include improper shifting in camera
tuning (electronics), mechanics of the rotating gantry (e.g. the camera
head may not be exactly centred in the gantry, refer to Figure 9.5) and
misaligned attachment of the collimator to the detector.
c. COR can be assessed by placing a point or line source within the field-
of-view (FOV) of a camera with a collimator attached. The position
122 Problems and Solutions in Medical Physics
of the source is off-centre (at some distance away from the central
axis of the detector). If a line source is used, it is placed parallel to the
axis of rotation (AOR). A set of projection images are then acquired
from different projection angles for 360° and stored in a 64 × 64
matrix. The position of the source along the X- and Y-axis from the
centre is then computed. The X-axis plot should give a symmetrical
bell-shaped curve, and the Y-axis plot should be a straight line pass-
ing through the pixel on which the source is positioned.
9 • Quality Control in Nuclear Medicine 123
Solution
a. Partial volume effect is the loss of apparent activity in small objects
or regions when the object partially occupies the sensitive volume
of the imaging device (in space or time). In gamma imaging, when
a hot spot relative to a ‘cold’ background is smaller than twice the
spatial resolution of the camera, the activity around the hot object
is smeared over a larger area (pixel in two dimensions [2D], and
voxel in three dimensions [3D]) than it occupies in the reconstructed
image. While the total counts are preserved, the object appears to be
larger and has a lower activity concentration than its actual value.
Analogously, a cold spot relative to a hot background would appear
smaller with higher activity around the object.
b. Partial volume effect can be corrected by applying a correction factor,
known as the recovery coefficient, on the image. The recovery coef-
ficient is the ratio of the reconstructed count density to the true count
density of the region of interest smaller than the spatial resolution of
the system. The recovery coefficient can be determined by measuring
the count densities of different objects containing the same activity but
with sizes larger and smaller than the spatial resolution of the system.
124 Problems and Solutions in Medical Physics
Solution
One of the most common used phantoms in SPECT QC is the ‘Jaszczak’
phantom. It can be used for semi-quantitative assessment of uniformity, spatial
resolution and image contrast of the system. This phantom may be filled with a
solution of the imaging radionuclide and is usually used for acceptance testing
and periodic QC testing of SPECT systems.
Radiation
Protection
in Nuclear
10
Medicine
Solution
a. Equivalent dose (HT ) is a measure of the radiation dose to tissue tak-
ing into consideration the different relative biological effects of dif-
ferent types of ionising radiation. The equivalent dose in tissue, T, is
given by the expression:
125
126 Problems and Solutions in Medical Physics
HT = ∑W ⋅ D
R T, R
where:
DT,R = absorbed dose averaged over the tissue or organ T, due to
radiation R
WR = radiation weighting factor due to radiation R
HE = ∑W ⋅ H
T T
where:
HT = equivalent dose in tissue or organ T
W T = weighting factor for tissue T
a. Radiation workers
b. Apprentices or students
c. Members of the public
10 • Radiation Protection in Nuclear Medicine 127
Solution
a. 20 mSv per year, averaged over defined periods of five years, with
no single year exceeding 50 mSv.
b. 6 mSv per year.
c. 1 mSv per year.
10.3 CLASSIFICATION OF
RADIATION WORK AREAS
PROBLEM
Define the following classification of radiation work areas:
a. Clean area
b. Supervised area
c. Controlled area
Solution
a. Clean area is a work area where the annual dose received by a
worker is not likely to exceed the dose limit for a member of the
public, that is, 1 mSv·y−1.
b. Supervised area is a work area for which the occupational expo-
sure conditions are kept under review, even though specific protec-
tive measures and safety provision are not normally needed. The
area must be demarcated with radiation warning signs, and legible
notices must be clearly posted.
c. Controlled area is a work area where specific protection measures
and safety provisions could be required for controlling normal
exposures or preventing the spread of contamination during normal
working conditions and preventing or limiting the extent of poten-
tial exposures. Annual dose received by a worker in this area is
likely to exceed 3/10 of the annual occupational dose limit. The
area must be demarcated with radiation warning signs, and legible
notices must be clearly posted.
128 Problems and Solutions in Medical Physics
Solution
ΓA
Exposure rate =
d2
where:
Γ = gamma constant (R·cm2 [mCi·h]−1)
A = activity of the source (mCi)
t = time (h)
d = distance from the source (cm)
Therefore,
cm 2 mCi
0.45 R ⋅ 1110 MBq ×
mCi ⋅ h 37 MBq
Exposure rate = 2
(100 cm )
= 0.00135 R·h−1
= 1.35 mR·h−1
Solution
R ⋅ cm 2 27 mCi
2.2 370 GBq ×
mCi ⋅ h GBq
Exposure rate at 1 m =
(100 cm )2
= 2.198 R/h = 2198 mR·h−1
2198
mR⋅h −1 = 219.8 mR·h−1
10
A factor of 220 is required to reduce the exposure rate from 2198 to 10 mR·h−1.
In terms of HVL,
2n > 220
While 27 = 128 and 28 = 256, the number of HVLs required to reduce the
exposure rate to 10 mR·h−1 is 8.
Given that 1 HVL = 3 mm, therefore the thickness of Pb required to
reduce the exposure rate to 10 mR·h−1 is 8 × 3 mm = 24 mm.
Solution
An approximate dose from a small source can be calculated as following:
ΓAt
Dose =
d2
130 Problems and Solutions in Medical Physics
where:
Γ = gamma constant (mSv·cm2 [MBq·h]−1)
A = activity of the source (MBq)
t = time (h)
d = distance from the source (cm)
Therefore,
h
1.3
mSv
MBq ⋅ h
(
)
3.7 ×106 MBq 15 min ×
60 min
Dose = = 120 mSv
(100 cm )2
Note: Although this is an approximation, it is adequate given the usual uncer-
tainty with the individual’s exact distance from the source and the time spent near
it. While the radiation dose has exceeded the annual dose limit for a radiation
worker (20 mSv per year averaged for a period of five consecutive years), this
incident would result in no significant acute medical consequence to the person-
nel, as the acute radiation syndrome is only noticeable for a dose of approximately
2 Sv. Considering the stochastic effects, a dose of approximately 1 Sv would
increase the chance of cancer by 5% compared to the non-exposed population.
a. Two weeks after the radioiodine therapy, the patient discovered that
she was pregnant. How would you estimate the gestation age?
b. Assuming that the foetus was four weeks old during the administra-
tion of I-131, what is the estimated dose to the foetus? Given that the
S-value (uterus ← thyroid) at 12th week is
Solution
a. The gestation age can be estimated from the first day of the last
menstrual period (LMP) of the mother. If the LMP is unknown or
not regular, then the gestation age can be estimated using ultrasound
measurements of the foetus combined with the dates of first foetal
heart tones and other developmental milestones.
b. Effective half-life, Te = (Tp·Tb) (Tp + Tb)−1
= (8 d × 80 d) (8 d + 80 d)−1
= 7.3 d
= 63720 s
INITIAL
AT TIME ACTIVITY, FRACTION Ã = A0 ⋅ fh ⋅1.44Te
(h) A0 (Bq) UPTAKE, fh Te (s) (Bq·s)
2 1.85 × 108 0.05 1.44 630,720 8.40 × 1012
24 1.85 × 108 0.24 1.44 630,720 4.03 × 1013
24 2.04 × 109 0.25 1.44 630,720 4.62 × 1014
48 2.04 × 109 0.15 1.44 630,720 2.77 × 1014
Total à 7.88 × 1014
D = Ã.S
= (7.88 × 1014 Bq·s) (3.70 × 10−17 Gy[Bq·s]−1)
= 29.2 × 10−3 Gy
~30 mGy
Solution
For most nuclear medicine diagnostic procedures, there is no necessity for
pregnant staff to take any additional precautions other than limiting their
direct contact to the radioactive source. As the radiation exposure from the
patients who have been administered with radiopharmaceuticals is quite low,
there is no radiological reason for the pregnant worker to refrain from under-
taking the imaging procedures. According to the Basic Safety Standards,
notification of pregnancy shall not be considered a reason to exclude a female
worker from work. However, the pregnant worker should refrain from han-
dling the high-dose therapeutic radioactivities, such as I-131 therapy for thy-
roid cancer.
Solution
Effective dose = ∑wT HT
where:
wT = tissue weighting factor
HT = equivalent dose to organ tissue T
134 Problems and Solutions in Medical Physics
The ICRP 103 recommendation dose limit for radiation workers is 20 mSv
per year.
Thus, the worker could still receive a whole-body equivalent dose of 4.8 mSv
in that particular year.
Solution
i. Decay storage: This method is usually used for radionuclides with
physical half-lives less than 120 days. The radionuclides are allowed
to decay in storage for a minimum period of 10 half-lives. After
10 half-lives, if the radioactivity cannot be distinguished from the
background activity, it can be disposed of in the normal waste after
removing all the radiation labels.
ii. Release into sewage drain system: A small amount of water-soluble
radioactive material may be disposed of in a designated sink, and
the total amount of activity does not exceed the regulatory limits.
Detailed records of all radioactive material disposals must be kept
for inspection by regulatory bodies. Radioactive excreta from
patients may be exempted from these limits and may be disposed
into the sanitary sewer. Items that are contaminated with radioactive
excreta (e.g. linen and diapers) have to follow the same limitations.
10 • Radiation Protection in Nuclear Medicine 135
Solution
i. Choose wet decontamination over dry decontamination, because
the dry method could create airborne dust hazards.
ii. Use mild decontamination with non-abrasive agents first, as other
methods could damage the surface involved.
iii. Take precautions to prevent further spread of contamination during
decontamination operations.
iv. Isolate and separate contamination with short-lived activity to allow
it to decay naturally.
Solution
The following advice would be given to the patient:
Solution
Advice to the nursing staff would include:
Solution
• According to the ICRP-94 recommendation, the activity limit
for burial and cremation of patients following I-131 treatment is
400 MBq.
• The preparation for burial or cremation should be supervised by a
competent authority or radiation protection officer (RPO).
• The funeral directors will need to be advised of any necessary pre-
cautions, and notification of the relevant competent authorities is
required.
• Relatives and friends should not be allowed to come into close con-
tact with the corpse.
• Clear communication is needed between the authorities, hospital
staff, funeral directors and the family members to ensure that ade-
quate controls are implemented without compromising dignity.
• All personnel involved in handling the corpse should be instructed
by the RPO and monitored.
• All objects, clothes, documents, etc. that might have been in contact
with the deceased must be monitored for any contamination.
• It may be expedient to wrap the corpse in waterproof material
immediately after death to prevent the spread of contaminated body
fluids.
• Embalming of the corpse should be avoided.
• Autopsy of a highly radioactive corpse should be avoided.
Bibliography
Bailey, D.L., Humm, J.L., Todd-Pokropek, A., van Aswegen, A., 2014. Nuclear
Medicine Physics: A Handbook for Teachers and Students. Vienna, Austria:
IAEA.
Cherry, S.R., Sorenson, J.A. and Phelps, M.E., 2012. Physics in Nuclear Medicine
e-Book. Philadelphia, PA: Elsevier Health Sciences.
De Lima, J.J. Ed., 2016. Nuclear Medicine Physics. Boca Raton, FL: CRC Press.
Flower, M.A. Ed., 2012. Webb’s Physics of Medical Imaging. Boca Raton, FL: CRC
Press.
Griffiths, H.J., 1988. “Radiation dose to patients from radiopharmaceuticals.” ICRP
Publication 53. Radiology, 169(3), pp. 652–652.
International Atomic Energy Agency. Nuclear Medicine Physics: A Handbook for
Teachers and Students. Vienna, Austria: IAEA, 2014.
International Commission on Radiological Protection, 1998. Radiation Dose to
Patients from Radiopharmaceuticals. ICRP Publication 80, Addendum to ICRP
53. Oxford, UK: Pergamon Press.
International Commission on Radiological Protection, 2000. “Pregnancy and medical
radiation.” ICRP Publication 84. Ann ICRP, 30(1), pp. 1–43.
Lin, G.S., Hines, H.H., Grant, G., Taylor, K. and Ryals, C., 2006. “Automated quantifi-
cation of myocardial ischemia and wall motion defects by use of cardiac SPECT
polar mapping and 4-dimensional surface rendering.” J Nucl Med Technol 34(1),
pp. 3–17.
Powsner, R.A., Palmer, M.R. and Powsner, E.R., 2013. Essentials of Nuclear Medicine
Physics and Instrumentation. Chichester, UK: John Wiley & Sons.
Russell, J.R., Stabin, M.G., Sparks, R.B. and Watson, E., 1997. “Radiation absorbed
dose to the embryo/fetus from radiopharmaceuticals.” Health Phys. 73(5),
pp. 756–769.
Saha, G.B., 2012. Physics and Radiobiology of Nuclear Medicine. New York: Springer
Science & Business Media.
Theobald, T. Ed., 2011. Sampson’s Textbook of Radiopharmacy. London, UK:
Pharmaceutical Press.
139