Hybrid Imaging Instrumentation and Data Processing

Download as pdf or txt
Download as pdf or txt
You are on page 1of 30

REVIEW

published: 18 May 2018


doi: 10.3389/fphy.2018.00047

Hybrid Imaging: Instrumentation and


Data Processing
Jacobo Cal-Gonzalez 1*, Ivo Rausch 1 , Lalith K. Shiyam Sundar 1 , Martin L. Lassen 1 ,
Otto Muzik 2 , Ewald Moser 1,3 , Laszlo Papp 1 and Thomas Beyer 1
1
QIMP Team, Center for Medical Physics and Biomedical Engineering, Medical University Vienna, Vienna, Austria,
2
Department of Radiology, Children’s Hospital of Michigan, The Detroit Medical Center, Wayne State University School of
Medicine, Detroit, MI, United States, 3 MR Center of Excellence, Medical University of Vienna, Vienna, Austria

State-of-the-art patient management frequently requires the use of non-invasive imaging


methods to assess the anatomy, function or molecular-biological conditions of patients
or study subjects. Such imaging methods can be singular, providing either anatomical
or molecular information, or they can be combined, thus, providing “anato-metabolic”
information. Hybrid imaging denotes image acquisitions on systems that physically
combine complementary imaging modalities for an improved diagnostic accuracy and
confidence as well as for increased patient comfort. The physical combination of formerly
independent imaging modalities was driven by leading innovators in the field of clinical
Edited by:
research and benefited from technological advances that permitted the operation of
Claudia Kuntner, PET and MR in close physical proximity, for example. This review covers milestones
Austrian Institute of Technology, of the development of various hybrid imaging systems for use in clinical practice and
Austria
small-animal research. Special attention is given to technological advances that helped
Reviewed by:
Jose M. Perez, the adoption of hybrid imaging, as well as to introducing methodological concepts that
Centro de Investigaciones benefit from the availability of complementary anatomical and biological information, such
Energéticas, Medioambientales y
Tecnológicas, Spain
as new types of image reconstruction and data correction schemes. The ultimate goal of
Marc Huisman, hybrid imaging is to provide useful, complementary and quantitative information during
Medical Center, VU University patient work-up. Hybrid imaging also opens the door to multi-parametric assessment
Amsterdam, Netherlands
of diseases, which will help us better understand the causes of various diseases that
*Correspondence:
Jacobo Cal-Gonzalez currently contribute to a large fraction of healthcare costs.
[email protected]
Keywords: hybrid imaging, combined imaging, instrumentation, nuclear medicine, data processing, data
corrections
Specialty section:
This article was submitted to
Biomedical Physics,
a section of the journal
INTRODUCTION
Frontiers in Physics
Since the discovery of X-rays by Wilhelm Conrad Roentgen in 1895 [1], non-invasive medical
Received: 28 December 2017 imaging has become a standard tool for the diagnosis and staging of numerous diseases. X-ray
Accepted: 30 April 2018
Computed Tomography (CT) and Magnetic Resonance (MR), first introduced in the early [1]
Published: 18 May 2018
and late 1970s [2, 3], respectively, are the most widely used tomographic imaging techniques
Citation: for depicting morphological changes of the human anatomy [4–6]. Metabolic or functional
Cal-Gonzalez J, Rausch I, Shiyam
changes, which may occur without a corresponding change of anatomy, can be depicted by
Sundar LK, Lassen ML, Muzik O,
Moser E, Papp L and Beyer T (2018)
functional imaging, which has proven to provide essential information for the diagnosis and staging
Hybrid Imaging: Instrumentation and of many diseases. The first tomographic functional imaging modality was Single Photon Emission
Data Processing. Front. Phys. 6:47. Tomography (SPECT), introduced in the early 1960s [7], followed be the first Positron Emission
doi: 10.3389/fphy.2018.00047 Tomography (PET) system in 1972 [8] and the first MR system in 1977 [3].

Frontiers in Physics | www.frontiersin.org 1 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

By combining anatomical and functional imaging within Clinical SPECT/CT, PET/CT and PET/MR
a single, hybrid imaging system, complementary diagnostic Systems
information can be obtained in order to gather a comprehensive Figure 1 shows representative designs for SPECT/CT, PET/CT,
picture of the disease. First attempts for obtaining such “anato- and PET/MR systems available on the market today, and Table 1
metabolic images” [9] were based on sophisticated software summarizes their most important technical specifications.
techniques to co-register structural and functional information
[10, 11]. In the late 1990s, imaging systems that combine two SPECT/CT
complementary imaging techniques within the same gantry (e.g., The first combined SPECT/CT design was proposed by
PET/CT, SPECT/CT) became available [12, 13]. This approach is Mirshanov in 1987 [31]; however, it took a decade until
known as “hardware fusion”, in contrast to the software fusion SPECT/CT became commercially available following some key
approaches mentioned above. contributions by Blankespoor et al. [12]. Since then, SPECT/CT
This review briefly describes the developmental paths of has advanced rapidly and several commercial system designs are
hardware fusion systems and discusses their future in clinical available today.
routine and research. In section Basic Concepts of Hybrid Two facets of the design of the integrated CT components
Imaging we describe the basic concepts of clinical and preclinical can be identified: first, SPECT/CT systems include fully-
hybrid imaging systems. Section Hybrid Imaging Technology diagnostic CT systems with fast-rotation detectors that permit
highlights technological advances, such as novel detectors for the simultaneous acquisition of 16 or 64 detector rows, while
PET, SPECT, and CT, time-of-flight PET and organ-specific the X-ray tubes provide sufficient tube voltages together with
or total-body PET systems. In section Data handling we high tube current and automatic exposure control (General
discuss the data handling in hybrid imaging systems, including Electric Discovery NM/CT 670 and the Siemens Symbia [32]).
data acquisition, data storage, image reconstruction and data These systems also support advanced CT capabilities such as
correction techniques. Section Joint Data Exploration focuses cardiac gating, calcium scoring and iterative reconstruction [18,
on the potential of joint data exploration, and section Multi- 33].
Center Standardization summarizes ongoing efforts regarding On the other hand, SPECT/CT systems include rotating
the standardization of hybrid imaging systems. Finally, a SPECT components that come with adapted rotation times
summary of the state-of-art and an outlook on future trends and dose-optimized acquisition modes, while employing CT-
in hybrid imaging is presented in section Outlook and Future type components with limited acquisition flexibility, such as
Trends. lower tube voltages and currents (General Electric Hawkeye,
Philips Brightview-XCT [32]). The Hungarian company Mediso,
for example, offers a triple-modality SPECT/CT/PET system
BASIC CONCEPTS OF HYBRID IMAGING (Mediso AnyScan), which combines all three modalities within
a single device (Figure 1).
Among the range of existing functional imaging techniques Most clinical SPECT systems are based on planar detectors
(functional magnetic resonance imaging—fMRI, perfusion consisting of two-dimensional (2D) array of photomultiplier
MR imaging, magnetic particle imaging—MPI, Near-infrared tubes (PMTs) attached to the back of the scintillation crystal
spectroscopy—NIRS, etc.), PET and SPECT present as very with en-face collimators. The location of a photon interaction
sensitive methods for the non-invasive and quantitative site is computed as the center-of-gravity of the position-
investigation of physiological processes at a molecular level. dependent energy signals from the PMTs according to the so-
Nevertheless, PET and SPECT provide mostly functional called Anger logic [34]. In recent years, alternatives to the Anger
information that may not always be directly associated with logic have been proposed (see section Detector Technology
well-defined anatomical structures. The lack of high-contrast for CT, SPECT, and PET). However, the comparably high
anatomical information in either SPECT or PET image data costs involved, mainly from the use of solid state detectors,
(independent of the radiotracer) is a major limitation of these restricts adoption in smaller systems designed for special
imaging techniques. In order to maximize the potential of applications, such as cardiac [35], brain [36] or pre-clinical
PET or SPECT, it has been recognized that either nuclear imaging [37].
medicine modality could be combined with a high-resolution The benefit of SPECT/CT has been proven for a wide
anatomical imaging modality. In that regard, PET/CT and range of clinical applications [38–41]. The main advantages
PET/MR were developed with the aim of aligning functional of SPECT/CT include improved attenuation correction and
and anatomical information to improve the clinical outcome accurate anatomical allocation of the SPECT/CT findings, both
of these studies, while SPECT/CT was conceived primarily resulting in better diagnostic performance. Moreover, combined
for the purpose of providing routinely CT-based attenuation SPECT/CT imaging has demonstrated its value particularly in the
and scatter correction of SPECT data [14–16]. Nonetheless, clinical management of patients with cardiovascular disease [35].
the adoption of CT-based attenuation and scatter correction
has proven to yield similar benefits for SPECT and PET PET/CT
alike; these include shorter transmission times as well as The very first PET/CT prototype was proposed in 1984 at Gunma
higher-quality data for post-injection transmission imaging University in Japan [5]. The two tomographs were situated next
[17]. to each other with the patient table moving sideways between the

Frontiers in Physics | www.frontiersin.org 2 May 2018 | Volume 6 | Article 47


TABLE 1 | PET, SPECT, CT, and MR specifications of selected dual and triple modality clinical systems commercially available.

SPECT/CT PET/CT PET/SPECT/ PET/MR


CT
Cal-Gonzalez et al.

Company Siemens GE Philips Siemens Philips GE Toshiba Kindsway United imaging Mediso Siemens Philips GE
biotech

System name Intevo Discovery BrightView mCT and Ingenuity Vereos Discovery Discovery Celestion PoleStar uMI 510 uMI 780 Anyscan mMR Ingenuity Signa
NM/CT 670 XCT mCT flow PET/CT 690 IQ m660 PET/MR

References [18] [19, 20] [21] [22] [23] [24] [25] [26] [27] Vendor Vendor [28] [29] [30]
webpage

Frontiers in Physics | www.frontiersin.org


PET Scintillator – – – LSO LYSO LYSO LYSO BGO LYSO LYSO LYSO LYSO LYSO LSO LYSO LBS
specifications
Crystal size (mm3 ) – – – 4 × 4 × 20 4 × 4 × 22 4 × 4 × 19 4.2 × 6.3 × 6.3 × 6.3 × 4 × 4 × 12 3.63 × 3.63 2.35 × 2.35 × 15 2.76 × 2.59 × 18 3.9 × 3.9 × 20 4 × 4 × 20 4 × 4 × 22 4.2 × 5.3 × 25
25 30 × 20
Total detector elements – – – 32 448 28 336 23 040 13 824 11 520 30 720 37 632 110 592 101 920 39 672 28 672 28 336 20 160
Photo-detector – – – PMT PMT dSiPM PMT PMT PMT PMT PMT SiPM PMT APD PMT SiPM
ToF capability – – – Yes Yes Yes Yes No Yes Yes Yes Yes No No Yes Yes
Patient bore (cm) 78.0 71.7 70.0 70.0 74.0 88.0 N.A. 70.0 70.0 70.0 60.0 70.7 60.0
Transaxial FOV (cm) – – – 81.5 67.6 N.A. 70.0 70.0 70.0 N.A. 70.0 70.0 55.0 59.4 N.A. 60.0
Axial FOV (cm) – – – 21.8 18.0 16.4 15.7 26.0 19.6 N.A. 23.6 30.0 23.0 25.8 18.0 25.0
Energy window (keV) – – – 435–650 440–665 N.A. 425–650 425–650 425–650 425–650 N.A. LLD, 430 N.A. 430–610 460–665 425–650
Energy resolution (%) – – – 11.5 11.1 11.1 12.4 N.A. 12.4 N.A. N.A. N.A. N.A. 14.5 11.6 11.0
Time coincidence – – – 4.1 4.5 4.0 4.9 9.5 1.6–3.2 4.1 N.A. N.A. 5.0 5.9 6.0 4.6
window (ns)
Time resolution (ns) – – – 0.5 0.5 0.3 0.5 N.A. 0.4 0.43 0.485 N.A. N.A. 2.9 0.5 0.4
Transaxial resolution 1 – – – 4.4/4.9 4.8/5.1 4.1/4.5 4.7/5.0 4.5/5.4 5.1/5.1 3.76/4.56 2.85/3.07 Max resol. 2.9 4.1/4.9 4.3/5.0 4.7/5.1 4.2/5.2
cm/10 cm (mm)

3
Axial resolution 1 – – – 4.4/5.9 4.7/5.2 4.0/4.3 4.7/5.6 4.8/4.8 5.0 /5.4 3.64/5.29 3.01/2.97 4.2/5.1 4.3/6.6 4.6/5.0 5.8/7.1
cm/10 cm (mm)
Sensitivity (kcps/kBq) – – – 9.7 7.3 5.7 7.4 22.8 3.8 10.9 8.3 16 8.1 15.0 7.0 21.0
Scatter Fraction (%) – – – 33.2 36.7 30.0 37.0 36.2 37.3 N.A. 38.4 40.0 N.A. 37.9 at peak 26.0 43.6 at peak
NECR NECR
Peak NEC (kcps @ – – – 180 @ 28 124 @ 20.3 171 @ 50 139 @ 29 124 @ 9.1 70 @ 29.6 224.6 @ 109 @ 21.5 170 @ 16.0 150 @ N.A. 184 @ 23.1 89 @ 13.7 210 @ 17.5
kBq/mL) 29.0

SPECT Detector type 3.8 in. NaI 3.8 in. NaI 3.8 in. NaI – – – – – – – – – NaI – – –
specifications
Photo-detector PMT PMT PMT – – – – – – – – – 60(48) PMT – – –
Detector size (cm) 38.7 × 53.3 40.6 × 54.0 40.0 × 54.0 – – – – – – – – – 58.5(55.8) × – – –
47.0(41.8)

CT specifications Max CT slices 2,6,16 16 Cone beam 128 128 128 64 16 16 64 16 128 16 – – –
CT tube max voltage 130 140 120 140 140 140 140 140 135 140 140 140 140 – – –
(kVp)
CT tube max current 345 440 80 800 665 665 800 440 600 667 420 833 500 – – –
(mA)
Max CT rotation (s) 0.5 0.5 12 0.3 0.4 0.4 0.35 0.5 0.5 0.39 0.5 0.3 0.4 – – –

MR specifications PET/MR integration – – – – – – – – – – – – – Fully integrated Sequential Fully integrated


Magnet – – – – – – – – – – – – – Superconductor Superconductor Superconductor
Magnetic field – – – – – – – – – – – – – 3T 3T 3T
Magnet length (cm) – – – – – – – – – – – – – 163 157 N.A.
Magnet bore (cm) – – – – – – – – – – – – – 60 60 60
Maximum FOV (cm3 ) – – – – – – – – – – – – – 60 50 × 50 × 45 50 × 50× 50

N.A., data not available.


Hybrid Imaging: Instrumentation and Data Processing

May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

FIGURE 1 | Images of commercially available hybrid imaging systems for clinical use. The figure shows three SPECT/CT, five PET/CT, three PET/MR and one triple
modality SPECT/PET/CT systems. Systems include the first CTM-PET(/CT) system (Siemens mCT Flow, see section Organ-Specific System Design and Total-Body
Systems); two SiPM-based PET systems (GE Signa, uMI 780), the first digital SiPM based PET system (Philips Vereos); a BGO-based system (GE Discovery IQ) and
the first commercially available fully-integrated PET/MR (Siemens mMR). Images taken from vendor’s web pages and published with the permission of the copyright
holders (the respective vendors).

two units. The first whole-body PET/CT prototype was proposed vendors have been presented following the advances in CT and
by Townsend and colleagues in the late 1990’s [13, 42]. Over PET instrumentation. To date, several vendors worldwide offer a
the years, PET/CT designs from various medical imaging system broad range of PET/CT designs.

Frontiers in Physics | www.frontiersin.org 4 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

The success of PET/CT imaging is based on several factors: Quantification in PET and SPECT
first, an anatomical and functional whole-body survey in a One of the most important advantages of tomographic nuclear
single session is logistically efficient for both the patient imaging techniques is the ability to accurately quantify the
and the healthcare provider. Second, due to the interplay of amount of radio-labeled biomolecules (radiotracers) in vivo. In
complementing data streams, the diagnostic information of a PET, this ability is based on the properties of the positron
PET/CT scan is superior to that of PET or CT alone [43, 44]. A emission coincidence detection, which allows for correction of
third advantage is the possibility to use CT information to correct photon attenuation in the emission signal by using information
PET data for photon attenuation and scatter or utilize it to correct gained from a separate transmission measurement [66]. As a
for partial volume effects [16]. result, emission images represent tracer concentration in units
of Bq/mL. However, for the assessment of functional processes,
PET/MR the sole knowledge of tracer concentration is not sufficient. The
An interesting alternative to PET/CT is to combine PET with tracer concentration in, for example, an organ depends on the
Magnetic Resonance Imaging (MR), since the range of MR amount of available radiotracer, which has led to the introduction
examinations is complementary to that of PET. MR imaging of the standardized uptake value (SUV) in the clinical reports.
reveals structural, functional and metabolic information through The SUV is a semi-quantitative measure that normalizes the
the interaction of three different magnetic fields (static, gradient, measured tracer concentration in tissue to a surrogate of the
dynamic) with the protons present in the tissues [2, 3]. The available tracer concentration in arterial blood / plasma. The
wide variety of imaging sequences, along with better soft-tissue most commonly used surrogate is the ratio of injected radiotracer
contrast compared to CT, renders MR an efficient diagnostic tool. to patient weight. The basic expression for SUV is [67]:
The potential of reducing patient dose by employing MR in lieu of
CT can be also considered an advantage compared to CT, mainly AMeasured
SUV = (1)
for pediatric examinations and cardiac imaging. AInjected /BW
However, the physical combination of PET and MR represents
a major technological challenge [45–48]. Conventional PET where AMeasured is the measured activity decay corrected to the
systems use PMTs to detect the scintillation light. These PMTs time of injection (kBq/mL). AInjected is the activity injected into
are highly sensitive to magnetic fields and, therefore, cannot the patient (kBq/mL) and BW is the body weight of the patient
be operated inside an MR magnet. One possible approach to (kg).
overcome this issue is to spatially separate the PET and the MR As of today, SUV remains a standard metric for PET-
system in combination with an active shielding of the PMTs based diagnosis and therapy response assessment in oncology
against the magnetic field from the MR. This design was proposed [68, 69]. However, the SUV represents only a simplified
in 2010 for the first whole-body PET/MR, the Philips Ingenuity measure to describe a physiological process. For absolute
TF PET/MR system [29]. Alternative approaches make use of quantification, i.e., the determination of physiological parameters
the potential of solid state photo-detectors, such as APDs or (such as metabolic rate of glucose consumption in mmol/g/min),
SiPMs [49–52], see section Detector Technology for CT, SPECT, pharmaco-kinetic modeling is required [68]; this entails the
and PET for details. This technology enables the design of knowledge of the time course of tracer concentration in arterial
fully integrated PET/MR systems that permit the simultaneous blood as well as in the tissue. From that, the pharmaco-kinetic
acquisition of PET and MR data within the same axial field-of- parameters can be derived from a temporal relationship between
view (aka co-planar FOV). Today, two integrated systems are these functions [70]. However, both the requirement of invasive
available commercially: the Siemens Biograph mMR [28] and the arterial blood sampling and the complexity of kinetic modeling
GE Signa PET/MR [30]. are significant limitations for the routine adoption of absolute
quantification in clinical routine so far [71].
Preclinical Hybrid Imaging In contrast to PET imaging, and given the more challenging
Animal models of human disease are the basis of many approaches toward attenuation and scatter correction, SPECT
research efforts to understand disease processes and the has not been considered a quantitative modality until recently.
development of new pharmaceuticals [53]. Similar to human However, methodological advances have led to an increase in
imaging, the combination of molecular imaging techniques the number of papers reporting high quantitative accuracy with
and CT or MR has been proven beneficial for small-animal SPECT [18, 72].
imaging [54]. Therefore, a large variety of dedicated small-
animal hybrid systems, has been developed since the 1990s, HYBRID IMAGING TECHNOLOGY
with all hybrid combinations (PET/CT, PET/MR, SPECT/CT,
SPECT/MR, and PET/SPECT/MR) commercially available. This section reviews the basic constituent technologies for
Figure 2 shows a number of dual- and triple-modality systems hybrid imaging systems. We start with a discussion on radiation
for small animal imaging that are available from several vendors. detectors used in CT, SPECT, and PET systems (section
Table 2 summarizes the main system specifications. A complete Detector Technology for CT, SPECT, and PET). Time-of-flight
description of all the available preclinical hybrid systems is and continuous table motion PET are discussed in sections
outside the focus of this review paper, please see [46, 61–65] for Time-of-Flight PET and Spiral PET, Continuous Table Motion,
details. respectively, while in section Organ-Specific System Design and

Frontiers in Physics | www.frontiersin.org 5 May 2018 | Volume 6 | Article 47


TABLE 2 | PET, SPECT, CT and MR specifications of selected dual and triple modality preclinical systems.

PET/CT SPECT/PET/CT SPECT /MR PET/MR

Company Sofie Sedecal MI Labs Mediso Siemens Bruker FLEXTM Mediso Mediso Bruker MR solutions
Cal-Gonzalez et al.

biosciences

System name GNEXT SuperArgus VECTor5 NanoScan Inveon Albira Si Triumph NanoScan NanoScan Bruker PET MRS-PETTM
PET/CT 4R SPECT/PET/CT SPECT/PET/CT SPECT/PET/CT SPECT/PET/CT SPECT/MR PET/MR /MR

References Vendor [55] [56–58] Vendor [59] Vendor [60] Vendor webpage Vendor Vendor
webpage webpage webpage webpage webpage

Frontiers in Physics | www.frontiersin.org


PET Detector type Pixelated Pixelated Uses the same Pixelated LYSO Pixelated LSO Continuous Pixelated - Pixelated LYSO Continuous Dual-layer
specifications LYSO/BGO LYSO/GSO detector (122S/82S) block detector LYSO/LGSO (122S/82S) block detector LYSO matrix
system than coupled to a 12 coupled to a 12 with 1/2 pixel
SPECT, with a × 12 SiPM × 12 SiPM offset
clustered- array array
pinhole
collimator
designed for
high-energy
photons, which
makes the
coincidence
detection not
needed.
Crystal size 1.01 × 1.01 × 1.45 × 1.45 × 1.12 × 1.12 × 1.5 × 1.5 × 10 N.A. 2.0 × 2.0 × – 1.12 × 1.12 × N.A. N.A.
(mm3 ) 6.1 mm3 LYSO (7+8) 13/1.51 × 1.51 (12+14) 13/1.51 × 1.51

6
1.55 × 1.55 × × 10 × 10
8.9 mm3 BGO
Total #detector 23 040/10 240 32 448 36 504/13 456 25 600 24 6 144 – 36 504/13 456 24 N.A.
elements
Photo-detector PMT PMT N.A. PMT SiPM PMT – N.A. SiPM N.A.
Depth-of- 2 layers 2 layers Yes (Monte No Yes 2 layers – Yes (Monte Yes Yes
interaction LYSO/BGO LYSO/GSO Carlo - based) LYSO/LGSO Carlo - based)
(DOI)
Transaxial FOV 13.0 12.0 1.2 central 12.0/8.0 10.0 8.0 11.0 – 12.0/8.0 8.0 7.0
(cm) (hourglass-
shaped)
Axial FOV (cm) 10.4 10.0 0.9 central 10.0 12.7 14.8 (3 rings 10.0 – 10.0 14.8 (3 rings 13.5 (3 rings
system) system) system)
Energy window N.A. 250–700 409–613 N.A. 350–625 350–650 N.A. – N.A. 350–650 N.A.
(keV)
Max spatial ? 1.0 0.8 (OSEM3D) 0.6 0.9 (OSEM3D) 1.4 (FBP) 0.7 (OSEM3D) 0.9 (MLEM) – 0.9 (OSEM3D) 0.7 (OSEM3D) ?0.9 (OSEM3D)
resolution (mm)
Sensitivity (%) N.A. 8.9 0.3 8.0/7.0 3.2 12 (3 rings N.A. – 8.0/7.0 12 (3 rings N.A.
system) system)
Scatter fraction N.A. N. A. N.A. N.A. 13.6/19.2 N.A. N.A. – N.A. N.A. N.A.
(%): mouse/rat
Peak NEC: N.A. N.A. N.A. 850 @ N.A. and 1394 @ 146 560 @ 35 and N.A. – 850 @ N.A. and 560 @ 35 and N.A.
mouse and rat 230 @N.A./460 and 560 @ 97 330 @ 43 230 @N.A./460 330 @ 43
(kcps @ MBq) @ N.A. and 130 @ N.A. and 130
@ N.A. @ N.A.

(Continued)
Hybrid Imaging: Instrumentation and Data Processing

May 2018 | Volume 6 | Article 47


TABLE 2 | Continued

PET/CT SPECT/PET/CT SPECT /MR PET/MR

Company Sofie Sedecal MI Labs Mediso Siemens Bruker FLEXTM Mediso Mediso Bruker MR solutions
Cal-Gonzalez et al.

biosciences

System name GNEXT SuperArgus VECTor5 NanoScan Inveon Albira Si Triumph NanoScan NanoScan Bruker PET MRS-PETTM
PET/CT 4R SPECT/PET/CT SPECT/PET/CT SPECT/PET/CT SPECT/PET/CT SPECT/MR PET/MR /MR

References Vendor [55] [56–58] Vendor [59] Vendor [60] Vendor webpage Vendor Vendor
webpage webpage webpage webpage webpage

Frontiers in Physics | www.frontiersin.org


SPECT Detector type – – 3 head NaI(Tl) 4-head 2-head – 5-by-5 CZT Continuous NaI – – –
specifications continuous NaI pixelated NaI modules
Photo-detector – – 55 PMT per PMT PMT – – PMT – – –
detector
Detector size – – 50.8 × 38.1 28 × 28 15.0 × 15.0 – 12.7 × 12.7 28 × 28 – – –
(cm)

CT CT voxel size 20 20 80 10 20 90 50 – – – –
specifications (µm)
CT tube max N.A. 110 65 70 80 50 N.A. – – – –
voltage (kVp)
CT tube max N.A. N.A. 615 N.A. 500 1,000 N.A. – – – –
current (µA)
Max CT ? 60 15 40 36 N.A. N.A. N.A. – – – –

7
rotation (s)

MR PET/MR – – – – – – – Sequential Integrated and Sequential


specifications integration PET insert
Magnet – – – – – – – Superconductor Superconductor Superconductor
Magnetic field – – – – – – – 1.5/3T 3T 3 or 7 T
Magnet bore – – – – – – – N.A. 18 17 /24
(cm)
Maximum – – – – – – – 100 N.A. N.A.
resolution (µm)
Maximum FOV – – – – – – – N.A. N.A. Eliptical 10 ×
(cm3 ) 7/13.5 × 9.8
(3T) 11 ×
7/15.4 × 9.8
(7T)

Of note, the MI Labs, Mediso and Bruker companies also offer SPECT/CT and PET/CT systems with the same sub-system characteristics than the ones shown in the SPECT/PET/CT configurations. The Inveon SPECT/PET/CT system
was discontinued and it is not commercially available anymore. MR solutions also offer and SPECT/MR system with multi-pinhole SPECT technology and similar MR capabilities than the ones shown in their PET/MR. N.A., data not
available.
Hybrid Imaging: Instrumentation and Data Processing

May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

FIGURE 2 | Selected commercially available, dual- and triple-modality preclinical systems. Images taken from vendor’s web pages and published with the permission
of the copyright holders (the respective vendors).

Total-Body Systems we discuss organ-specific and total body also characterized by the highest light output (i.e., number of
hybrid imaging systems. Finally, in section MR Technology we photons emitted per unit of deposited energy), which is directly
provide a brief overview of the technology used in MR systems. related to the energy resolution. Most PET and SPECT systems
are composed of pixelated scintillators to assign the interaction
Detector Technology for CT, SPECT, and position of the gamma photon [74], although blocks made of
PET continuous crystal have been proposed as well [75, 76]. Table 3
The most commonly used radiation detectors for CT, SPECT, and lists some of the scintillators commonly used in CT, SPECT and
PET imaging systems are based on scintillator crystals, because PET detector systems together with their key properties [77–
they are fast, they provide a high stopping power for photons 81]. New scintillator materials are being developed continuously.
across a range of emission energies and their cost-per-volume Some examples include scintillators based on cerium doping
rate is far superior to other detector materials. A scintillation of lanthanide and transition metal elements, such as LuAP:Ce,
detector consist of a crystal that produces scintillation light after CeBr3, LuBO3:Ce, and others based on lead (Pb), tungsten (W),
the interaction with radiation and a photo-detector that converts and gadolinium (Gd) [82].
the scintillation light into an electrical signal [73], which is then Photomultiplier tubes (PMTs) represent the most common
processed by the electronic system. mean to measure and detect the scintillation light emitted by the
There are several scintillating materials that are currently scintillator crystal. They consist of a vacuum enclosure within
being used. They can be classified into organic or inorganic a thin photocathode layer at the entrance window and several
scintillators that come in solid, liquid or gaseous state. For electrodes, called dynodes, which amplify the electrical signal
medical imaging applications, solid inorganic scintillators are created in the photocathode by means of secondary electron
preferred, due to their higher density (i.e., stopping power), emissions. A PMT produces an electric pulse with amplitude
which is required for high sensitivity imaging systems. They are proportional to the number of scintillation photons that reach

Frontiers in Physics | www.frontiersin.org 8 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

TABLE 3 | Physical properties of some scintillators used in medical imaging applications (adapted from Lewellen [77]).

NaI BGO GSO LSO LYSO LaBr3 BaF2 LuAP

Composition NaI:Tl Bi4 Ge3 O12 Gd2 SiO5 :Ce Lu2 SiO5 :Ce (LuY)2 SiO5 :Ce LaBr3 :Ce BaF2 LuAP:Ce
Density (g/cm3 ) 6.7 7.1 6.7 7.4 7.1 5.3 4.9 8.3
Effective atomic number 51 74 59 66 60 47 54 65
Attenuation coefficient (cm−1 ) 0.34 0.92 0.62 0.87 0.86 0.47 0.44 0.90
Refractive index 1.85 2.15 1.85 1.82 1.81 1.88 – 1.95
Light yield (%NaI) 100 15 41 75 80 160 5 16
Wavelenght for max. Emission (nm) 410 480 430 420 420 370 220 365
Decay constant (ns) 230 300 56 40 41 25 0.8 18
Hygroscopic Yes No No No No No Slight No

TABLE 4 | Characteristics of photo-detectors used in medical imaging matrix of small APD cells biased to be operated above avalanche
applications (adapted from Lecomte [79]). breakdown in the so-called Geiger mode. Since the Geiger-mode
PMT APD SiPM
operation yields a high gain (105 –106 ), a multi-cell structure
can provide a proportional output for moderate photon flux by
Active area (mm2 ) 1–2,000 cm2 1–2,000 cm2 1–10 mm2 summing the signal of all cells that have been activated. In 2009,
Gain 105 –106 102 105 -106 Philips Healthcare introduced the digital SiPMs, also known as
Dynamic range 106 104 103 Digital Photon Counter (DPC), which combines a conventional
Excess noise factor 0.1–0.2 >2 1.1–1.2 array of Geiger-mode photodiodes with a fully-digital electronic
Rise time (ns) <1 2–3 1 readout system, thus, allowing timing resolutions of up to 60 ps
Time jitter (ns FWHM) 0.3 >1 0.1 Full Width Half Maximum (FWHM) [84–87].
Dark current <0.1 nA/cm2 1–10 nA/mm2 0.1–1 MHz/mm2 Semiconductor detectors represent the main alternative
Photon detection efficiency 25 60–80 <40 to scintillator-based detectors in CT and SPECT imaging
@ 420 nm systems. In comparison to scintillators the main advantage of
Bias-voltage (V) 1,000–2,000 100–1,500 >100 semiconductor detectors is the direct conversion from radiation
Power consumption 100 mW/ch 10 µW/mm2 <50 µW/mm2 to an electric pulse, which avoids the degrading effects associated
Gain dependence with <1%/◦ C 2–3%/◦ C 3–5%/◦ C with scintillation light production, propagation and conversion
temperature to an electrical signal in the photodetector. Their use in PET
Gain dependence with <1%/V 10%/V 100%/V systems is much less frequent, due to their lower stopping power
voltage for the high-energy 511 keV annihilation photons and their
Magnetic susceptibility Very high (mT) No (up to 9.4 T) No (up to 15 T) higher costs.
To date, the most widely investigated semiconductors for
nuclear medicine and CT imaging are CdTe and CdZnTe (CZT).
Both offer a relatively high stopping power (similar to the one
the photocathode, which, in turn, is proportional to the deposited offered by NaI:Tl scintillator crystals) and they can be operated
energy. at room temperature. Furthermore, with the use of highly
The main advantage of PMTs as scintillation light detector granulated detectors high spatial, energy and time resolution
is its high amplification capability (in the order to 106 –107 ). can be achieved [88]. Table 5 summarizes important properties
On the other hand, their most significant drawback of PMTs of some semiconductor materials used for medical imaging
is their sensitivity to magnetic fields, which makes them not applications [80].
suitable for use in combined PET/MR or SPECT/MR systems
[83]. A valid alternative to PMTs are solid-state detectors such Time-of-Flight PET
as Avalanche Photodiodes (APDs) or Silicon Photomultipliers Another important technological advance for PET systems,
(SiPMs). Solid-state detectors have several inherent advantages referred to as Time-of-Flight (ToF) PET, was made possible
over PMTs (Table 4), such as high quantum efficiency, compact with the introduction of very fast detectors. The concept of ToF
and flexible shape, ruggedness, demonstrated insensitivity to assumes that the annihilation point of two photons originating
magnetic fields up to 9.4 T and potentially inexpensive mass from a single positron annihilation can be calculated from their
production [49–51]. Solid-state detectors are semiconductor travel time differences [89–91].
devices with a low-field depleted region where visible or near- State-of-the-art PET systems achieve time resolutions of
UV photons can create electron-hole pairs by photoelectric effect. about 500 ps FWHM (using photo multiplier tubes) [19,
APDs exist as small discrete devices or as monolithic arrays, 29], down to about 300 ps FWHM (using SiPMs) [22,
which can be used for individual or multiplexed crystal readouts. 30], which corresponds to an uncertainty of determining
Silicon Photomultipliers (SiPMs) consists of a densely packed the origin of the annihilation along the Line Of Response

Frontiers in Physics | www.frontiersin.org 9 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

TABLE 5 | Physical properties of some semiconductor detectors used in medical technologies were made in recent years. For example, Siemens
imaging applications (adapted from Peterson and Furenlid [80]). Healthcare did propose a “Brain PET” system that was based on
Si Ge CdTe CdZnTe
a PET ring insert for a 3 T MR system [104]. More recently, a
SiPM-based PET insert has been developed for integration with
Density (g cm−3 ) 2.33 5.32 5.85 5.82 any MR scanner [105]. Alternative efforts aim at developing MR-
Attenuation @ 140 keV (cm−1 ) 0.02 0.72 3.22 3.07 compatible PET inserts for breast imaging. Initial prototypes
Energy e-h pair (eV) 3.61 2.98 4.43 5 have been proposed [106] with promising results, and a breast
Electron mobility lifetime (cm2 V−1 ) 0.42 0.72 0.003 0.003 hybrid PET/RF insert is being developed based on digital SiPMs
Hole mobility lifetime (cm2 V−1 ) 0.22 0.84 0.0005 0.00005 (dSiPMs) for enhanced diagnosis of breast cancer (HYPMED
2016: http://www.hypmed.eu/). Similar approaches have been
proposed for the diagnosis of prostate and pancreatic cancer,
whereby PET and Ultrasound (US) components are combined in
(LOR) of about 15 and 9 cm, respectively. With this additional an endoscope device [107].
information, the SNR in the reconstructed PET images can be The concept of a large axial FOV (aFOV) whole-body PET
improved significantly [92]. Thanks to recent improvements in was introduced by Crosetto in the 1990’s [108]. Following the
detector technology, new scintillator materials and fast digital introduction of this concept and after several simulation studies,
SiPMs photodetectors [93], ToF resolution can be significantly the EXPLORER (Extreme Performance Long Research scanner)
improved; the latest developments in PET detector technology consortium was set up with the task of building the world’s first
aim at timing resolutions of 100 ps FWHM, or less [92, total-body PET/CT system—the EXPLORER PET/CT [109, 110].
94]. The EXPLORER has a 200 cm aFOV based on mainstream PET
detector technology. It will consist of a total of 400,000 crystals,
Spiral PET, Continuous Table Motion thus, yielding about 100 times more LORs than a state-of-the-art
Commercially available whole-body PET systems offer fixed whole-body PET(/CT) system.
axial FOVs between 15 cm and 25 cm. To cover larger axial The elongated axial coverage of the EXPLORER type system
areas, multiple bed positions are acquired in a step-and- is expected to yield a 30–40 times increase in sensitivity
shoot protocol [95]. An alternative to the traditional step- over current generation PET systems. This translates to
and-shoot approach is a continuous table motion (CTM) detecting over 40% of the counts emitted from a point
during the acquisition. This technique was first proposed source in air located at the center of the FOV [111]. This
in 1992 by Dahlbom et al. [95] aiming to increase the impressive increase in sensitivity provides a number of key
uniformity of the sensitivity across the examined axial area. advantages: imaging at very low radioactivity levels (∼10 MBq)
First attempts to implement CTM in a clinical PET only [112], increased throughput from shorter examination times
system were made as early as 2001 [96] and followed by and whole-body imaging that can be performed within a
several implementations of CTM techniques in PET/CT systems single breath hold. With the total-body coverage provided
a few years later [97, 98]. However, it took until 2013 by the EXPLORER, it is also possible to perform whole-
to implement CTM into a commercially available system body dynamic PET imaging, thus, avoiding temporal gaps,
[19]. which are present in current multi-bed, multi-pass imaging
CTM offers two main advantages in comparison to step- protocols [113, 114]. However, efficient image reconstruction
and-shoot acquisitions: a uniform sensitivity across the entire and data handling schemes are required for the EXPLORER
axial FOV and the possibility to customize the axial scan range. system. Zhang et al. [112] proposed a quantitative image
While in step-and shoot acquisitions uniform sensitivity can reconstruction method that demonstrated a 7-fold increase
also be achieved by an appropriate selection of bed overlap, in the signal-to-noise ratio in comparison to current PET
the axial scan range is restricted to a discrete number of systems.
bed positions. Another advantage of CTM acquisitions is the Independent to the EXPLORER, a high-resolution 100
possibility of new acquisition protocols that support whole- cm-aFOV PET system, called “PET20.0,” is currently being
body parametric imaging [99]. Given the above benefits, CTM configured in cooperation between Ghent University and
can be regarded a particular asset in combined PET/CT Vrije Universiteit Brussels [115]. In comparison to the PET
imaging. EXPLORER, PET20.0 has a shorter aFOV of 100 cm. However,
with monolithic detector crystals and improved positioning
Organ-Specific System Design and methods, the spatial resolution of PET20.0 (2.0 mm) is expected
Total-Body Systems to be better than the spatial resolution of the EXPLORER
The development of state-of-the-art detectors allows the design (3.5 mm).
of high-resolution organ-specific hybrid systems, as well as high- In contrast to the aforementioned PET systems, a lower-cost
sensitivity total body PET/CT with extended axial field-of-view alternative based on resistive plate chamber (RPC) detectors has
(aFOV). been proposed also. In addition to being economic and reliable,
Although most attempts toward organ-specific imaging were RPC-PET detectors offer good spatial-temporal resolution, high
based on stand-alone PET or SPECT systems [100–103], time-of-flight resolution and high energy sensitivity. Further,
significant efforts in fully-integrated, organ-specific imaging the RPC detectors allow for the accurate measurement of

Frontiers in Physics | www.frontiersin.org 10 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

Depth-of-interaction (DOI), which renders RPC-PET parallax- There are two main categories for Image Reconstruction (IR)
error free. These properties make RPC-detectors an excellent algorithms in CT, SPECT and PET: analytical reconstructions and
option for an extended aFOV PET system [116]. iterative reconstructions [124–126]. Analytical approaches are
based on a closed-form equation that directly yields one solution
MR Technology based on the input (Figure 3A), while iterative techniques are
Currently over 60% of clinical MR systems operate at 1.5 Tesla based on a more accurate description of the imaging process,
for cost and safety reasons [3, 117]. To date, the conventional thereby resulting in a more complicated mathematical solution
upper limit for the field strength of clinical MR is 3 Tesla. The that requires multiple steps to achieve an image. Iterative
growing adoption of MR systems with 7 Tesla (over 80 units methods model the data collection process in a tomographic
installed worldwide) could be seen to represent the new frontier system and search for the image that is most consistent with the
of MR imaging1 . Of note, singular testing of MR systems with measured data (Figure 3B).
10.5 T (University of Minnesota, USA) and 11.7 T (CEA, Saclay,
France) is work-in-progress [3], but it is unlikely that these Analytical IR
ultra-high-field MR systems will soon be employed in hybrid The filtered-backprojection (FBP) reconstruction [127] method
imaging. is the standard method for image reconstruction of CT data.
For hybrid imaging of human subjects and patients, i.e., The FBP method combines a back-projection operation, which
combined PET and MR, actively shielded 3 T magnets have been describes the propagation of the measured projection data into
chosen by all manufacturers. In addition to the wide bore magnet, the image domain, with a filter component that compensates
fast gradient coils (slew rate 200 mT/m/ms) and dedicated RF for the low-pass blurring inherent to the back-projection
coils (proton frequency at 128 MHz) are required [3, 118]. approach. With the advances in CT technology, different
Recently, maximum gradient strength of up to 300 mT/m were adaptations of the FBP algorithm (interpolation methods, 3D-
realized at 3 T [119], although state-of-the-art clinical systems reconstruction methods: Feldkamp algorithm, etc. [128–132])
operate with a gradient strength of 70 mT/m, which is sufficient have been proposed to compensate for fan-beam and cone-beam
for most routine pulse sequences prescribed. Dedicated RF- geometries. Likewise, exact [128, 133] and approximate [134]
coil arrays are being developed to further improve sensitivity 3D reconstruction methods have been proposed for helical cone
and speed and to avoid costly ultra-high field magnets, thus, beam reconstruction, which are being used in most commercial
rendering a suitable option for combined PET/MR systems [120] CT systems due to their flexibility and computational efficiency
as well as for various parallel imaging techniques [121–123]. [125].
Parallel imaging enables more efficient data sampling to either
increase acquisition speed or the amount of information collected Iterative IR
per time interval. Iterative methods offer improvements over the analytical
approach because they allow the noise structure of the data to be
DATA HANDLING accounted for, thus, incorporating a more realistic model of the
system. All iterative reconstruction methods consist of three steps
This section provides an overview of standard data acquisition that are iterated until a convergence condition is met (Figure 3B).
and image reconstruction techniques for functional (SPECT, First, an estimated raw data set is created from an estimated initial
PET) and anatomical (MR, CT) imaging modalities (section image using a realistic model of the tomographic system. Second,
Image Acquisition and Reconstruction in CT, PET, and SPECT the artificial raw data is compared with the real data acquired by
to Quantitative Data Corrections in PET and SPECT). Advanced the system in order to obtain a correction term for each bin in the
data processing methods are described in sections Anatomically- projected data. Finally, these correction terms are back-projected
Driven PET/SPECT Image Reconstruction to Outlook and and used to update the image estimate.
Future Trends, which make use of the anatomical information Iterative reconstruction methods can be classified in
provided by the CT or MR images to improve the PET-SPECT deterministic and statistical approaches [135]. The most
quantification). frequently used deterministic approaches are the Algebraic
Reconstruction Technique (ART) [136], the simultaneous ART
Image Acquisition and Reconstruction in (SART) [137], and the ordered-subsets Simultaneous Iterative
CT, PET, and SPECT Reconstruction technique (OS-SIRT) [138]. Statistical methods
A CT image represents the tissue-dependent attenuation of incorporate counting statistics of the detected photons into
X-rays of the investigated object, thereby making use of an the reconstruction process. The most well-known methods are
X-ray source of known source strength and an opposite detector the Maximum likelihood expectation-maximization (ML-EM)
array, both of which rotate at a fixed speed around the center [139], its accelerated version by means of the ordered-subsets
of the FOV. Alternatively, SPECT and PET images represent expectation-maximization (OSEM) [140], and the model-
the distribution of activity concentrations (per voxel) of a single based iterative reconstruction (MBIR) [141], which is used
photon (SPECT) or a positron (PET) emitter. predominantly in multi-slice helical CT.
Image reconstruction methods in SPECT and PET are similar
1 Siemens Healthineers (Erlangen, Germany) obtained both, CE labeling and FDA to the methods employed in CT although the physical nature of
approval for their 7 T Terra system in 2017. the acquired data is different. These differences relate mainly to

Frontiers in Physics | www.frontiersin.org 11 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

FIGURE 3 | (A) Simplified view of the direct FBP image reconstruction technique: the projected data is filtered using different filter kernels and back-projected into the
image domain. Multiple projections are required to obtain the final CT image. (B) Schematic view of the iterative reconstruction process: first, a forward projection of
the initial estimated image is used to create the estimated projected data; then the estimated data is compared to the acquired raw data and a set of correction
factors is derived. These correction factors are back-projected and used to update the initial estimated image. All three steps describe an iterative loop, which is
repeated until a predefined condition is satisfied and the final image is generated. Similar reconstruction processes apply for SPECT and PET data.

Frontiers in Physics | www.frontiersin.org 12 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

the geometrical model (e.g., CT—fan-beam, cone-beam or helical to 16-bit digital resolution) and visualized on a computer screen
vs. PET—cylindrical rings of detectors vs. SPECT—rotating flat as gray scale values/images of tissue specific image contrast.
detector panels plus collimators), and, thus, to the level of In a mathematical framework, a Fourier transformation of the
mandatory a priori data corrections. detected signal is stored in k-space (2D-space consisting of all
Statistical IRs have been used extensively in SPECT and frequencies and phases detected in the measured RF signal) and
PET due to their benefits with regards to noise reduction and then transformed into image space (x,y) via the inverse Fourier
improved accuracy in the reconstructed images [124]. However, transformation. The result is a spatially resolved magnetization
these approaches are not yet widely adopted in clinical CT. This is vector M(r).
due to the numerous pre-processing and calibration steps in CT,
which change the statistical properties of the measured data; the Quantitative Data Corrections in PET and
high spatial resolution in CT, which requires the application of SPECT
edge-preserving regularization techniques; and the large amount In order to obtain quantitative and artifact-free PET or SPECT
of data that requires long reconstruction times. Nonetheless, images, several corrections must be applied to the acquired data,
the exponential growth of computer technology and the recent which include, for example: decay correction, photon attenuation
introduction of IR methods that can be implemented on a correction (AC), scatter correction, normalization, dead time
Graphic Processing Unit (GPU, [142]) render the use of iterative correction and randoms subtraction (only in PET). Here we
reconstruction in CT more and more a clinical commodity. will discuss the most widely methods used for quantitative data
corrections in PET and SPECT images, and their application in
PET/CT and SPECT/CT systems. In PET/MR, as the MR image
information is not related to the attenuation properties of the
Image Acquisition and Reconstruction in material, novel methodologies for AC need to be developed, and
MR are discussed in detail in section Novel MR-Based Attenuation
MR imaging is a non-invasive tomographic imaging technique Correction (AC) Methods for PET/MR. For a general overview
that acquires data from the inside of an object using on the quantitative correction methods for PET we refer the
radio-frequency (RF) excitation of protons that subsequently reader to Zaidi [144], while for SPECT we suggest [72, 145].
retransmit the absorbed energy, depending on the biological
environment. Nuclear induction of RFs is performed using coils Attenuation and Scatter Corrections in PET/CT
or coil arrays, with a wavelength larger or similar to the object The emitted photons in PET and SPECT are subject to
dimensions (at 3 T/128 MHz the proton wavelength in biological attenuation as they travel through the patient. As consequence,
tissue is about 30–50 cm). In order to allow spatial encoding of the number of detected photons in each LOR will be reduced.
RF emission from bulk tissue, gradient coils need to be used that The most widely used method for AC in PET/CT was proposed
limit the resonance frequency to a small portion of image space. by Kinahan et al. [17, 146] and it is known as “bilinear”
This is commonly performed using three independent linear segmentation-scaling method. Here, the PET attenuation image
gradient fields, although other approaches, using, e.g., radial or is estimated by first using a threshold to separate the bone
spiral gradients, exist. Following a combination of excitations and component from the soft tissues of the CT image, and then using
measurements, imaging and contrast modules can be selected separate scaling factors for the bone and non-bone component.
independently, thus, making MR a highly versatile and flexible Alternatively, the emitted photons may be scattered in the patient
imaging technique with excellent soft tissue contrast [143]. The body or in the detector itself, suffering a deviation of its trajectory.
imaging module defines the details (i.e., matrix size, 2D/3D) of As a consequence, the LOR of a scattered event will not be
the data collection from an object of a given size in a given longer aligned with the emission point. The most extensively used
time (i.e., measurement time). Anatomical MR data acquisition method for scatter correction in PET/CT (and PET/MR) is the
is performed for the purpose of high-resolution imaging (i.e., Single Scatter Simulation (SSS) method [147, 148], where only a
matrix size 256 × 256 up to 512 × 512) in a matter of a single scatter event is considered and multiple or out of the FoV
few minutes. For functional data, reduced spatial resolution is scatter contributions are included as scaling factors.
acceptable (i.e., 64 × 64 up to 128 × 128) in view of acquiring
3D frames within 10–50 ms and with short repetition times, TR, Attenuation and Scatter Corrections in SPECT/CT
of 100–2,000 ms over the course of several minutes. Similar to in PET/CT, the bilinear method proposed for CT-based
Given the variety of pulse sequences and sequence parameters AC can be used in combined SPECT/CT systems [18, 149].
(e.g., TR, TE, MA, FOV, etc.) data processing and analysis On the other hand, the most common approach for estimating
pipelines, particularly in research, can be rather complex. In scatter events in SPECT/CT is the measure of counts in
addition, differences in pulse sequences and processing software additional energy windows adjacent to the photopeak window.
can be observed between manufacturers [118], often depending The most common examples are the Dual Energy Window
on hardware performance or Intellectual Property (IP) issues. By (DEW) approach that neglects upper scatter or the Triple Energy
combining short (ms) magnetic pulses produced by an RF-coil Window (TEW) [150–152]. More complex solutions, based on
with orthogonal magnetic gradients in x, y, z; voxel localization Monte Carlo estimations of the scattered events, have been also
in 2D or 3D is enabled and the voltage induced in the receive coil proposed [153, 154] and are expected to be widely used in the
can be detected. This voltage is converted into numbers (12-bit near future.

Frontiers in Physics | www.frontiersin.org 13 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

Randoms Subtraction in PET correction [168] and the geometric transfer-matrix (GTM)
Random events occur when two uncorrelated photons, from two method [169]. Other post-reconstruction corrections consist of
different decays, are detected in coincidence within the timing voxel-based methods, such as the popular Mueller-Gaertner
window. The randoms events can be estimated by using a delayed method (MGM) [170], image deconvolution [171], and the
coincidence window [155]. “region-based voxel-wise correction” (RBV) [172]. Alternative
approaches for PVC operate in sinogram space, where statistical
Normalization in PET and SPECT noise is spatially uncorrelated and easier to incorporate in
Because of variations in the gain of PMTs, inaccuracies/tolerances the PVC procedure [173–176]. Reconstruction-based methods
in detector block building, physical variation of scintillator enhance the spatial resolution by using the implicit PVC
efficiency, etc., the detection sensitivity varies from detector compensation through system-response modeling [177]. More
to detector. Information on these variations, often known as advanced PVC methods make use of anatomical priors during the
normalization [145, 156], is required for the reconstruction of reconstruction process [159, 161], as described in the previous
quantitative and artifact-free images. section.
Despite the large number of PVC methodologies proposed
Anatomically-Driven PET/SPECT Image in the literature, explicit corrections for PVE are still far from
Reconstruction being fully utilized in clinical routine [178]. This is due to several
Iterative IR methods for PET and SPECT are challenged in challenges associated with their automated implementation.
two ways. First, PET, and SPECT images have an inherently First, the accurate delineation of the tissues of interest from the
limited spatial resolution due to the physical process that morphological images is a challenging task, especially in cases
leads to the image formation. And, second, noisy images are where the metabolic tissue boundaries do not correspond well to
obtained from over-iterated reconstructions, given the lack of the morphologic boundaries. Second, some of the PVC methods
an objective stopping rule in the iterative reconstruction. The require segmentation of tissues into different VOIs that span
latter problem is addressed in clinical routine by using a Gaussian the entire cross-sectional extent of the subject, which can be
post-reconstruction filter, which further degrades the perceived computationally intensive.
resolution and contrast of the images in return for lower noise
levels. Motion Compensation
Bayesian methods attempt to improve the quality of the The hybrid imaging protocols for SPECT/CT, PET/CT, and
reconstructed image by taking advantage of prior knowledge PET/MR systems often span acquisition times of several minutes
[157], that can be obtained, for example, from a co-registered and, thus, cover multiple breath-holds or cardiac contractions.
anatomical image (CT or MR). This information is known a Therefore, these hybrid images are subject to involuntary patient
priori and is often incorporated into a maximum a posteriori motion during the acquisition [179]. The induced respiratory
(MAP) objective function [158]. These reconstruction methods motion as well as cardiac contraction introduce motion blur in
lead to improvements of contrast and noise properties of the the acquired images (Figure 5), which may affect the clinical
reconstructed SPECT and PET images [159–161]. This a-priori interpretation of the acquired images [182].
known anatomical information can be introduced into the Patient motion has traditionally been detected through the use
reconstruction algorithm by using probabilistic image models of external markers, such as ECG-electrodes for cardiac motion
[162] or similarity measures between the anatomical and the or respiratory belts [183] as well as infrared/optical systems [184]
functional images [163]. for respiration-induced motion. Recently, the use of data-driven
methods has gained substantial interest for both stand-alone
Partial Volume Correction in PET and and multi-modality imaging systems [185–188]. Continuous
SPECT motion tracking can be used for motion compensation of PET
The partial volume effect (PVE) is related to the limited and MR images [189, 190]. In thoracic PET imaging, several
spatial resolution of the nuclear imaging system and affects the motion detection techniques have been proposed, through built-
quantitative accuracy of the images, particularly for small lesions in readout of the respiratory position and special radial-self
and brain structures (Figure 4). PVEs are caused by two different gating sequences, respectively [186, 187, 191]. Cardiac motion is
phenomena: the limited spatial resolution of the imaging system most commonly estimated by using cine-MR imaging and tagged
(Figure 4A) and the discrete image sampling of a continuous 3D MR imaging [192].
activity distribution (Figure 4B). PVE results in smearing of the Approaches toward respiratory and/or cardiac gating are
reconstructed activity levels and, thus, has a significant impact on commonly used to reduce motion-induced blurring in PET
their qualitative and quantitative assessment. and SPECT images, whereby the emission data is divided
A wide variety of partial volume correction (PVC) methods into gates that correspond to different respiratory or cardiac
have been developed in the past, most of them employing high- phases [185]. However, a more preferable solution is to
resolution anatomical information from a co-registered CT or perform motion compensation (MoCo) using the full data-
MR image as reference [164–167]. PVC methods can be grouped acquisition, both to reduce acquisition times and to improve
into two main categories: post-reconstruction and during- the signal-to-noise ratios in the PET or SPECT images.
reconstruction methods. Post-reconstruction methods include Therefore, different MoCo techniques have been proposed:
the region-based methods, such as the recovery-coefficient MoCo before image reconstruction (projection-based methods)

Frontiers in Physics | www.frontiersin.org 14 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

FIGURE 4 | Illustration of the partial volume effect (PVE) in medical imaging. (A) A circular source of uniform activity in a warm background region yields a measured
image in which part of signal emanating from the source is observed outside (spill-out) while part of the background signal is observed within the lesion (spill-in).
(B) The tracer distribution is sampled on a voxel grid and the contours of the voxels do not match the actual contours of the tracer distribution (Tissue-fraction).

FIGURE 5 | Effect of respiratory and cardiac motion on PET image quality. (A) Blurring effect due to cardiac motion in a patient scanned with 13N-NH3.
(B) Misalignment artifact between the PET-emission data and the corresponding MR-AC map. Misalignment artifacts like these are common in cardiac PET-imaging,
where they are known to cause false-positive findings in up to 40% of all studies [180]. Reprinted from Rausch et al. [181] with permission from Elsevier.

[193], during PET-image reconstruction—motion compensation (RTA). Recent studies have been performed to evaluate the
during image reconstruction (MCIR) [194], and after PET- performance of these three methodologies [196, 197]. More
image reconstruction [195]—Reconstruction Transform Average advanced methods, that use the synergistic information derived

Frontiers in Physics | www.frontiersin.org 15 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

from the simultaneous PET and MR data acquisition, have been [225, 226], such as body truncation, metal and respiratory
proposed to estimate [198] and to correct for motion in both artifacts, MR coils attenuation, etc. Some promising solutions
PET and MR images [189]. Other attempts focused on the for whole-body AC are currently under development, for further
implementation of joint MoCo and PVC approaches to further details we refer the reader to Leynes et al. [223], Mehranian et al.
improve the PET quantification of small lesions in regions that [225], Heußer et a. [227].
are susceptible to motion [199, 200].
In addition to respiratory and cardiac motion, non-periodic
and unpredictable motion can also occur during the scan. JOINT DATA EXPLORATION
Organ specific motion detection techniques have been developed
This section describes main efforts toward a fully-synergistic
for head and neck examinations, through the use of frequent
use of the anatomical and functional information available from
3D-MR navigators, which can be employed to apply motion
hybrid imaging examinations. In section Kinetic Modeling and
compensation (MoCo) to the PET images [201].
Image Derived Input Function (IDIF) we discuss the use of
While most of the aforementioned references were related to
image-derived input functions in the context of kinetic modeling,
PET/CT and PET/MR imaging, it is important to note that most
in section Multi-Parametric Imaging (MPI) we provide an
of these techniques are applicable also to SPECT [202], while
overview of the multi-parametric imaging studies performed
standard data-driven motion detection techniques for PET do
using SPECT/CT, PET/CT, and PET/MR systems. Finally, in
not work well for SPECT [203]. Nevertheless, the more advanced
sections In Vivo Disease Characterization and Image-Derived
Principal Component Analysis (PCA) and Laplacian Eigenmaps
Prediction Models we explore recent efforts in obtaining
(LE) methods for data-driven respiratory motion detection have
in vivo disease characterization and image-derived prediction
demonstrated good results for both PET and SPECT imaging
models.
[203, 204].

Kinetic Modeling and Image Derived Input


Novel MR-Based Attenuation Correction Function (IDIF)
(AC) Methods for PET/MR As discussed in section Quantification in PET and SPECT,
In PET/MR examinations the anatomical information used for absolute PET quantification requires the measurement
attenuation correction (AC) originates from the complementary of an input function, which is typically done by means
MR image information, which is not related to the attenuation of an arterial blood sampling. The measurement of an
properties of the material [205]. Moreover, the imaging of cortical arterial input function (AIF) is invasive, laborious and
bone is a challenge in MR given the fast relaxation times of stressful for the patient. Obtaining an image-derived input
solid materials [206]. As a consequence, it is not possible to function (IDIF) is a non-invasive alternative, whereby
distinguish between cortical bone and air in most MR sequences the input function (IF) can be directly obtained by
[207]. To address this issue, new AC methods had to be defining a volume-of-interest (VOI) in the PET images
developed, that can be categorized into three basic concepts [181] and using anatomical information from the MR or CT
(Figure 6): (i) segmentation-based methods, which are based on images to segment the tissue of interest to obtain the IDIF
a segmentation of MR images into different tissue classes; (ii) (Figure 7A).
atlas- or model-based approaches, which incorporate a-priory The concept of an IDIF has been successfully implemented
knowledge of attenuation properties of the investigated subject in clinical routine examinations for cardiovascular studies, due
from data bases; and (iii) reconstruction-based methods, which to the availability of large blood pools (i.e., left ventricle)
are based on a direct reconstruction of the attenuation and in the PET FOV [229, 230]. However in PET studies of
activity values from the emission data. the brain, calculation of the IDIF is challenged by two factors:
At present, standard AC approaches for whole-body imaging PVE [231, 232] due to the small size of the blood pools,
in clinical systems rely on segmentation-based methods [218, and subject motion [232, 233]. Various approaches have been
219], recently also in combination with atlas-based approaches proposed to extract an accurate IDIF that can be classified
incorporating the major bone structures into the AC [220, 221]. into PET-only based methods [232, 234–244], stand-alone PET
In general, these approaches perform reasonably well, although and MR-based methods [228, 245, 246], and fully-integrated
for specific application their accuracy may not be sufficient PET/MR-based methods [247–249]. Combined PET/MR can
[222]. Thus, several advanced AC methods have been explored potentially allow researchers to address the aforementioned
in conjunction with various types of PET/MR examinations challenges (PVE and subject motion) in an automatic way [228,
[217, 223, 224]. Most of the developed AC methods are tailored 248, 249], in addition to bearing a logistic advantage. It should
to brain examinations and provide acceptable accuracies (with be noted that various corrections (delay, dispersion, metabolite
deviations below 5% respect to CT-based AC) [217]. A selection correction) may be required to convert the measured blood IDIF
of available MRI based AC methods for the head is given to plasma IDIF.
in Figure 6. A detailed comparison of these methods can be An IF can be employed together with a physiological model
found in Ladefoged et al. [217]. However, for AC of body parts in order to calculate physiological parameters, such as metabolic
other than the brain as well as for AC of non-rigid hardware rate of glucose (umol/100 g/min) or blood flow (ml/100 g/min)
components (e.g., flexible surface coils), challenges still persist [250]. The creation of such parametric images can be divided into

Frontiers in Physics | www.frontiersin.org 16 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

FIGURE 6 | Illustration of popular MR-based AC methods using the same axial slice of AC maps of a patient processed with different AC methods: (A) Reference
CT-based AC; Standard methods: (B) standard DIXON-based MR-AC implemented in the Siemens mMR system (SW v. VB20), (C) UTE-based standard MR-AC for
brain examinations implemented in the Siemens mMR system (SW v. VB20); Template-based methods: (D) Koesters et al. [208], (E) Anazodo et al. [209],
(F) Izquierdo-Garcia et al. [210], (G) Burgos et al. [211], (H) Merida et al. [212]; MLAA-based methods: (I) Benoit et al. [213]; Segmentation-based methods: (J)
Cabello et al. [214], (K) Juttukonda et al [215], (L) Ladefoged et al. [216]. Figure adapted from Ladefoged et al. [217] (Courtesy of Claes N. Ladefoged, Rigshospitalet
Copenhagen, Denmark; Modified from the original image published under the creative commons license http://creativecommons.org/licenses/by-nc-nd/4.0/). Figure
published with the permission of the copyright holders (Claes N. Ladefoged).

three different methods: (i) graphical analysis, (ii) compartment where the cerebellum or the pons is often used as reference
model, and (iii) reference tissue model. The first two methods tissues.
require information on the non-bound tracer-activity available in
the blood.
Graphical analyses methods include the Gjedde-Patlak Multi-Parametric Imaging (MPI)
equation [251, 252], which describes irreversibly bound tracers, Multi-parametric imaging (MPI) has been of interest since
such as [18 F]FDG and the Logan-plot, which can be used to the adoption of standalone PET, CT, and MR imaging, but
describe the reversibly bound tracers [253]. Tissue-compartment gained further attention following the availability of combined
models have been proposed by Kety and Smith in 1948 for the imaging systems [257]. Combined, or dual-modality imaging
calculation of the blood-flow in the brain using nitric oxide systems facilitate multi-faceted and complex evaluations of
[254, 255]. Figure 7B represents a 2-tissue compartment model, tumor phenotypes, thus, promising an improved characterization
which is the most appropriate to evaluate the tracer dynamics of lesions and pathophysiology [258].
of FDG. Here, K1 is the influx-constant (amount of tracer In addition to the anatomical information provided by CT-
moving from the blood to the tissue), k2 the efflux-rate (amount images, quantitative assessment of physiological parameters
of tracer that returns to blood), CP the blood activity and extracted from CT-images in combination with functional PET-
CT the tissue concentration over time t. Compartment 1 (C1 ) images increase the diagnostic value [259]. For example, texture
represents the un-metabolized tracer and compartment 2 (C2 ) analyses (TA) of the tumor heterogeneity obtained from CT-
the metabolized tracer (where k3 represents phosphorylation images is frequently used in combination with functional PET
by hexokinase and k4 dephosphorylation) [256]. Reference data [259]. However, the TA varies significantly with the quality
tissue models are primarily employed for neurological studies of the underlying image data, thus, limiting the reproducibility

Frontiers in Physics | www.frontiersin.org 17 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

FIGURE 7 | Schematic view of the IDIF concept. (A) 1- PET coronal view showing the carotid arteries, 2- T1-MPRAGE image, 3- TOF MR angiography image of neck
region used to segment the carotid arteries, 4- PET coronal view co-registered with the segmented carotid arteries. (B) Illustration of the 2-tissue compartment model
commonly used for FDG tracer kinetics. Here, CP represents the blood activity and C1 , C2 the activity concentration of each tissue over time t. K1, k2, etc., are the rate
constants that define the rate of tracer movement between compartments. This figure was adapted from Sari et al. [228], with permission of SAGE Publications, Ltd.

of the procedure, particularly in the context of PET imaging PET/MR systems can help to shed new light on the mechanisms
[260]. that underlie structural changes in the brain through the
Fully-integrated PET/MR systems permit the quasi- use of targeted radio-labeled tracers in combination with
simultaneous acquisition and evaluation of functional and gadolinium-based contrast enhancement of MR [262]. Likewise,
structural information of patients, and, thus, support the the combination of PET and MR images has proven to add
multi-parametric assessment of these subjects, as shown for prognostic value for assessments of myocardial examinations
the patients with neurological disorders [261]. Specifically, [189, 263, 264], while the value of the multi-parametric images
the use of parametric imaging [section Kinetic Modeling obtained for oncological studies is still under evaluation
and Image Derived Input Function (IDIF)] of data from [265].

Frontiers in Physics | www.frontiersin.org 18 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

In Vivo Disease Characterization examination. Quantification in nuclear medicine examinations


Hybrid imaging plays an important role in the characterization of is not only a matter of proper imaging system performance.
several diseases [4]. Thanks to hybrid imaging lesion detection, The outcome of a study is influenced also by various factors
tissue delineation and therapy monitoring have become more related to imaging technology, biology and imaging physics
accurate [266]. Hybrid imaging has been key to cancer patient [290]. In PET/CT operations, for example, a significant variability
management [267], thanks also to its ability to depict tumor of system performance, examination protocols, and quantitative
heterogeneity in a morphological and functional context [260]. reporting across different centers has been described [291]. A
Given that conventional indices, such as standard uptake value similar variability in workflow has been observed in SPECT/CT
(SUV), maximum or metabolic tumor volume (MTV) are operations [292]. Although in SPECT/CT most readings are
insufficient to characterize malignancies [268], several groups based on a qualitative assessment of the tracer distribution,
have started to investigate the feasibility of in vivo tumor significant variation in quantitative evaluations of SPECT/CT
characterization in the context of textural evaluation [269]. examination can be expected from the use of CT-AC.
As a result of this initiative, first reports have appeared that This variability in quantification between systems and across
point at promising results of a quantitative assessment of sites is a major barrier to the utilization of quantitative nuclear
tumor heterogeneity in light of therapy response prediction medicine and hybrid imaging in larger cohort and multi-
[270, 271], disease-specific survival [272] as well as prognostic center studies [293]. To gain reproducible and statistically
stratification [273]. Meanwhile, challenges using textural features significant results, e.g., for the evaluation of a new drug
remain, since such calculations are highly sensitive to acquisition, or the usability of an imaging modality to assess treatment
reconstruction and sample size variations [274–276]. Overall, the response, a sufficient number of examinations has to be collected.
need of reproducibility evaluation as well as standardization of Since this is usually not possible within a single institution,
textural features is being acknowledged in the field [274]. data from several centers need to be pooled following the
The approach that considers medical images as data mining harmonization of imaging protocols [293]. Several guidelines
sources for large-scale in vivo feature evaluations is called have been put forward [294–296] and specific programs, namely
radiomics [277]. It is expected that with the introduction the accreditation programs of the EANM or the ACR, have been
of hybrid imaging, radiomics holds the potential to lead to launched to sensitize the imaging community and to help address
automated tumor characterization [278–280] and personalized this challenge toward pooled data evaluation.
treatment regiments in the future [281]. For more details, we In general, most of the factors influencing simple quantitative
refer to the companion review by Papp et al.2 in this journal issue. readings, as SUV, are understood and—in many cases—
corrections and workarounds exist [290]. However, with the
Image-Derived Prediction Models availability of more advanced image evaluation techniques,
Machine learning (ML) algorithms have the ability to identify such as ML-based approaches employing textural parameters,
key patterns and characteristics of large scale datasets [282]. new challenges for a standardized evaluation arise. Textural
While the concept of ML is not new, it has only recently parameters are strongly influenced by the choice of image
become popular thanks to recent advances in computational reconstruction, voxel size and post processing steps (e.g.,
power and capacities that made it feasible to handle large Gaussian filtering). Furthermore, the choice and definitions of
datasets. ML holds also a great potential for dealing with “big the textural analysis metrics vary widely and the evaluation
data” generated by hybrid imaging methods [280]. For example, techniques differ substantially, thus, rendering a comparison
ML has been successfully applied in hybrid imaging to predict between studies almost impossible [260, 297]. The need for
survival [283], treatment outcome [284], and tumor grading standardization efforts in the field of textural analyses has been
[285]. Furthermore, unsupervised ML allows the identification highlighted previously [260, 298, 299], and efforts are currently
of breast cancer subtypes [286]. For more details, we refer to the under way to address this challenge [260, 300–302].
companion review by Papp et al.4 in this journal issue.
OUTLOOK AND FUTURE TRENDS
MULTI-CENTRE STANDARDIZATION
Open Research Data, Sharing Knowledge
The performance of a SPECT or PET system is related to factors Clinical research is an essential building block for the concept
as detector type and geometry, the properties of the built in of efficient patient management. Research studies are generally
electronics and the implemented data processing. To characterize complex and the resulting data are valuable, not only to the
the system performance, standardized procedures are published principal investigator but to society as a whole [303]. Although
by the National Electrical Manufacturers Association (NEMA) today’s ubiquitous multi-modality imaging protocols produce an
[287–289]. While these standards are a good example how to abundance of qualitatively diverse data sets, there is a lack of
set up a basic set of measurements for the purpose of system genuinely integrative analysis schemes that could provide added
characterization and comparison, their results are generally value to a pure additive analysis.
not suitable to predict the quantitative accuracy of a clinical Moreover, many researchers remain reluctant to share their
data with an expert audience [304, 305] beyond describing
2 Papp et al. Personalizing medicine through hybrid imaging and medical big data them as part of peer-reviewed publications. In contrast, sharing
analysis. research data in a structured and tangible way has been shown

Frontiers in Physics | www.frontiersin.org 19 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

to yield benefits for both, the principal investigators of the study disentangle the coincidences originating from each tracer [317,
who generated the data as well as other experts [306] in the field 318].
who may re-use the data with alternative evaluation algorithms
to extract new information that may subsequently benefit patient Dual-Modality Tracers/Contrasts
management [303]. Obviously, a reliable anonymization or even Contrast-enhanced CT data as part of a combined PET/CT
pseudo-anonymization of the relevant data must precede any examination provides additional, helpful information in
type of data accrual and accessibility as part of public data comparison to non-enhanced PET/CT studies [319–321]. Here,
repositories. Frequently, the challenges associated with prior data the main benefit relates to a more precise anatomical localization
handling alone render such repositories a mighty wish rather of pathologies by differentiating lesions from surrounding
than a useful reality. structures. On the other hand, MR provides highly-detailed
Still, after data are properly anonymized and made available non-invasive anatomical image information and excellent soft
under an open data policy, the quality of such public tissue contrast [322], even in comparison to contrast-enhanced
Supplementary Material collections of published studies, at least, CT.
is frequently variable, and many times the re-use of these data Efforts are under way to develop dual-modality contrast in
is not possible [307]. The same holds true for the quality of PET/MR based on nanoparticle-based probes. Currently, super
alternative public data archives that were shown to contain paramagnetic iron oxide nanoparticles (SPIONs) are used that
incomplete data and data archived only partially in over 56% incorporate positron emitters (e.g., 52-Mn, 59-Fe, and 124-I) and
cases that prevented re-use [305]. contain small surface molecules (peptides, organic molecules)
Nevertheless, future trends in multi-modality imaging point that bind to tumor tissue [323, 324]. Although the physiological
toward the emergence of advanced database structures that functionality of the PET probe is limited by its attachment
include interactive 3D visualization tools and powerful data to a large nanoparticle, the combination provides high spatial
mining applications. The leading candidate for generating resolution (MR) with high sensitivity (PET).
such flexible database structures is the XML (eXtensible With the advent of hybrid microPET/SPECT/CT systems,
Markup Language) model [308] due to its flexibility and efforts to develop multi-tracer PET-SPECT studies have been
scalability. The flexibility of XML database structures allows made. Although the scattered photons from PET preclude to
the combination of patient data with similar disease within a obtain useful information in SPECT with simultaneously injected
common reference frame, allowing meta-analysis of data patterns PET/SPECT tracers [325, 326], sequential studies where the
distributed over many modalities. An example of such an SPECT radiotracer is injected and imaged first and the PET
overarching effort is the National Database for Autism Research radiotracer is injected and imaged in the presence of the SPECT
(NDAR, http://data-archive.nimh.nih.gov/#NDAR-anchor), an tracer have shown promising results [326].
NIH-funded data repository that aims to accelerate progress in
autism spectrum disorder (ASD) research through data sharing, Hybrid Imaging and Therapy (Planning,
data harmonization, and the reporting of research results. Follow-Up)
One of the major developments in radiotherapy in the coming
Multi-Tracer PET and SPECT Imaging years is the use of hybrid imaging for individualized biology-
Another interesting application is the multi-tracer SPECT or guided radiation therapy (RT) [327]. Given the information
PET imaging, which aims to image two or more tracers in a on the tumors provided by PET/CT and PET/MR imaging,
single scan, simultaneously characterizing multiple aspects of dual-modality PET-based imaging represents an optimal
physiology and function without the need for repeat imaging basis for RT individualization, especially in the context of
visits [309]. Tracer separation in multi-tracer SPECT takes intensity-modulated radiotherapy (IMRT). The combination of
advantage of the differences in the spectro-temporal properties IMRT and hybrid imaging enables the modulation of radiation
(differences in the gamma emission energies and dynamic dose distribution according to local phenotypic or micro-
behavior) of the radionuclides to identify the tracer distributions environmental variations in an individual tumor (aka “dose
[310, 311]. painting”) [328, 329]. The goal of the dose painting process is to
Multi-tracer PET must rely on differences in the kinetics select and delineate target volumes (and organs at risk) on the
and/or spatial distributions of the tracers, because the energy basis of complementary diagnostic information, thus, aiding in
of annihilation photons of all tracers is 511 keV. The different heterogeneous delivery of radiation within the tumor volume
techniques to separate multiple-tracer PET may be based by targeting radio-resistant areas. This approach is particularly
on the analysis of the different radioactive half-lives of the promising in treatment monitoring, by taking into account
radionuclides [312, 313], multi-tracer compartment models [314] the individual patient response and dynamically adjusting the
or generalized methods that make use of PCA to separate the treatment plan based on current outcome. Given the broad
tracers [315]. A recent development in multi-tracer PET is the selection of established PET tracers, hybrid imaging provides
combination of a pure and a non-pure positron emitter [316]. RT planning with patho-physiological information pertaining to
This technique relies upon detecting the auxiliary prompt gamma various molecular pathways of the tumor, including metabolism,
in coincidence with an annihilation event in order to measure proliferation, oxygen delivery and consumption as well as
triple-coincidence events originating from only one of the tracers, receptor or gene expression. Examples are [18 F]FDG [330] and
and use the spatial distribution of the triple coincidences to [18 F]/[11 C]choline [331–333] as surrogates for tumor burden,

Frontiers in Physics | www.frontiersin.org 20 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

[18 F]FLT as a surrogate for proliferation (or cellular growth Hybrid imaging is supported by continuous cross-specialty
fraction) [334], and hypoxia sensitive tracers such as [18 F]FMISO efforts on integrating hardware components and leveraging
[335] and [64 Cu]ATSM [336] as surrogates of cellular complementary signals for improved image quality and
hypoxia. quantitative accuracy of non-invasive assessments of the patients
On other hand, PET monitoring is a popular method for and subjects.
range verification of particle therapy: proton or hadron therapy To date, multiple evidence exists for the benefits of hybrid
[337]. Positron emitters are generated inside the patient during imaging that extend from increased diagnostic accuracy, faster
treatment, and therefore, by using PET imaging the delivered examinations to cross-fertilization of know-how from formerly
dose can be indirectly measured [338, 339]. This can be done distinct professional groups.
during the irradiation (on-line or in-beam PET), which is Hybrid imaging is able to provide unique information on
beneficial when imaging short-lived emitters [340], or after the diseases that—when pooled in larger-scale studies—can be
treatment in an off-site system (off-line), which can provide utilized in a “big data” approach to design automated, computer-
better counting statistics for relatively long-lived radionuclides based models for disease and therapy response prediction,
[341]. In on-line systems, the integration of PET imaging both of which benefit also from ongoing technological and
into the treatment environment poses geometric constraints, methodological improvements of hybrid imaging systems.
conditioning the imaging performance. In off-line imaging, full- Hybrid imaging invites stakeholders from a variety of
ring scanners are used, offering higher detection efficiency. scientific disciplines in an effort to provide diagnostic means that
However, transferring the patient to a separate system may cause help us understand the complex biological processes involved in
alignment errors or isotope washout [337]. a disease better and, thus, help patients today and in the years to
come.
Future Trends
In our opinion, future trends in the field of hybrid imaging relate
AUTHOR CONTRIBUTIONS
to:
1. The development of new advanced imaging technologies that JC-G collected the data and wrote sections Introduction,
allow for more precise and accurate information, such as the Basic Concepts of Hybrid Imaging, Detector Technology for
use of faster ToF detectors, the design of total-body systems CT, SPECT and PET, Image Acquisition and Reconstruction
with very high sensitivity, the introduction of organ-specific in CT, PET, and SPECT; Quantitative Data Corrections in
systems with very high spatial resolution and new probes for PET and SPECT; Anatomically-Driven PET/SPECT Image
contrast-enhanced PET/CT and PET/MR. Reconstruction; Partial Volume Correction in PET and SPECT;
2. The development of advanced quantitative corrections for and Multi-Tracer PET and SPECT Imaging. He also merged
hybrid imaging data. This encloses accurate techniques for and revised all the individual contributions from the other
attenuation and scatter correction of PET and SPECT data in authors and prepared the full version of the manuscript
hybrid systems, more efficient multi-modality techniques for with all figures and tables. IR wrote sections Time-of-Flight
patient motion correction, as well as advanced and efficient PET; Spiral PET, Continuous Table Motion; Novel MR-Based
methodologies for PVC of the lower resolution functional Attenuation Correction (AC) Methods for PET/MR; and Multi-
images. Center Standardization. LS wrote sections Organ-Specific System
3. The development of efficient and clinically viable multi- Design and Total-Body Systems; Kinetic Modeling and Image
modality protocols and computational frameworks that Derived Input Function (IDIF); and prepared all the figures of
can handle various large data cohorts. An essential feature the manuscript. ML wrote sections Motion Compensation and
of such integration efforts is the quantitative description Multi- Parametric Imaging (MPI). OM wrote section Outlook
of data properties, so that relationships in anatomical and and Future Trends. EM wrote sections MR Technology, Image
functional domains between complementing modalities Acquisition and Reconstruction in MR, and contributed to
can be expressed in a mathematically concise way. Thus, others. LP wrote sections In vivo Disease Characterization
taking advantage of high-level integration schemes, and Image-Derived Prediction Models. TB conceived the topic,
quantitative results are combined into data structures prepared the abstract and conclusions from the full version of
that provide a consistent framework for the application the manuscript and revised all the individual contributions. All
of both machine learning and advanced data mining authors participated in manuscript review.
techniques.
ACKNOWLEDGMENTS
CONCLUSION
The authors would like to thank Miklos Kovacs for the useful
Hybrid imaging is the physical combination of complementary discussions about SPECT systems and configurations, Claes
anato-metabolic imaging modalities, and, to date, relates mainly N. Ladefoged for kindly providing us with the image data
to combined SPECT/CT, PET/CT, and PET/MR imaging, for Figure 6 and Jun Bao for providing us with the technical
which are available for clinical use and small-animal research. specifications of the uMI PET/CT systems.

Frontiers in Physics | www.frontiersin.org 21 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

REFERENCES 23. Bettinardi V, Presotto L, Rapisarda E, Picchio M, Gianolli L, Gilardi MC.


Physical performance of the new hybrid PET/CT Discovery-690. Med Phys.
1. Friedland GW, Thurber BD. The birth of CT. AJR Am J Roentgenol. (1996) (2011) 38:5394. doi: 10.1118/1.3635220
167:1365–1370. doi: 10.2214/ajr.167.6.8956560 24. Reynés-Llompart G, Gámez-Cenzano C, Romero-Zayas I, Rodríguez-Bel L,
2. Rinck PA. Magnetic Resonance in Medicine: The Basic Textbook of the Vercher-Conejero JL, Martí-Climent JM. Performance characteristics of the
European Magnetic Resonance Forum. 12th revised and enlarged edition. whole-body discovery IQ PET/CT system. J Nucl Med. (2017) 58:1155–61.
Hoboken, NJ: Blackwell Scientific Publications (2018). Available online at: doi: 10.2967/jnumed.116.185561
www.magnetic-resonance.org 25. Kaneta T, Ogawa M, Motomura N, Iizuka H, Arisawa T, Hino-Shishikura A.
3. Moser E, Laistler E, Schmitt F, Kontaxis G. Ultra-high field NMR and MRI— Initial evaluation of the Celesteion large-bore PET/CT scanner in accordance
the role of magnet technology to increase sensitivity and specificity. Front with the NEMA NU2-2012 standard and the Japanese guideline for oncology
Phys. (2017) 5:33. doi: 10.3389/fphy.2017.00033 FDG PET/CT data acquisition protocol version 2.0. EJNMMI Res. (2017)
4. Beyer T, Freudenberg LS, Townsend DW, Czernin J. The future of hybrid 7:83. doi: 10.1186/s13550-017-0331-y
imaging—part 1: hybrid imaging technologies and SPECT/CT. Insights 26. Montgomery DW, Amira A, Zaidi H. Performance evaluation of a new
Imaging (2011) 2:161–9. doi: 10.1007/s13244-010-0063-2 high-sensitivity TOF clinical PET/CT system. J Nucl Med. (2015) 56:433–3.
5. Beyer T, Townsend DW, Czernin J, Freudenberg LS. The future of 27. Xu B, Liu C, Dong Y, Tang R, Liu Y, Yang H, et al. Performance evaluation of
hybrid imaging—part 2: PET/CT. Insights Imaging. (2011) 2:225–34. a high-resolution TOF clinical PET/CT. J Nucl Med. (2016) 57:202.
doi: 10.1007/s13244-011-0069-4 28. Delso G, Fürst S, Jakoby B, Ladebeck R, Ganter C, Nekolla SG.
6. Beyer T, Freundenberg LS, Czernin J, Townsend DW. The future of hybrid Performance measurements of the siemens mMR integrated whole-body
imaging-part 3: Pet/mr, small-animal imaging and beyond. Insights Imaging PET/MR scanner. J Nucl Med. (2011) 52:1914–22. doi: 10.2967/jnumed.111.0
(2011) 3:189. doi: 10.1007/s13244-011-0136-x 92726
7. Jaszczak RJ. The early years of single photon emission computed tomography 29. Zaidi H, Ojha N, Morich M, Griesmer J, Hu Z, Maniawski P. Design
(SPECT): an anthology of selected reminiscences. Phys Med Biol. (2006) and performance evaluation of a whole-body Ingenuity TF PET-MRI
51:R99–115. doi: 10.1088/0031-9155/51/13/R07 system. Phys Med Biol. (2011) 56:3091–106. doi: 10.1088/0031-9155/56/
8. Nutt R. The history of positron emission tomography. Mol Imaging Biol. 10/013
(2002) 4:11–26. doi: 10.1016/S1095-0397(00)00051-0 30. Grant AM, Deller TW, Khalighi MM, Maramraju SH, Delso G, Levin
9. Wahl RL, Quint LE, Cieslak RD, Aisen AM, Koeppe RA, Meyer CR. CS. NEMA NU 2-2012 performance studies for the SiPM-based ToF-
“Anatometabolic” tumor imaging: fusion of FDG PET with CT or MRI to PET component of the GE SIGNA PET/MR system. Med Phys. (2016)
localize foci of increased activity. J Nucl Med. (1993) 34:1190–7. 43:2334–43. doi: 10.1118/1.4945416
10. Hutton BF, Braun M. Software for image registration: algorithms, accuracy, 31. Mirshanov DM. Transmission-Emission Computer Tomograph (1987). USSR
efficacy. Semin Nucl Med. (2003) 33:180–92. doi: 10.1053/snuc.2003.127309 Patent No. 621.386:616-07320.01.87-SU-181935.
11. Slomka PJ, Baum RP. Multimodality image registration with software: 32. Seret A, Nguyen D, Bernard C. Quantitative capabilities of four
state-of-the-art. Eur J Nucl Med Mol Imaging (2009) 36:44–55. state-of-the-art SPECT-CT cameras. EJNMMI Res. (2012) 2:45.
doi: 10.1007/s00259-008-0941-8. doi: 10.1186/2191-219X-2-45
12. Blankespoor SC, Wu X, Kalki K, Brown JK, Cann CE, Hasegawa BH. 33. Buck AK, Nekolla S, Ziegler S, Beer A, Krause BJ, Herrmann K, et al.
Attenuation correction of SPECT using X-ray CT on an emission- SPECT/CT. J Nucl Med. (2008) 49:1305–19. doi: 10.2967/jnumed.107.050195
transmission CT system: myocardial perfusion assessment. In: IEEE Nuclear 34. Anger HO. Scintillation camera. Rev Sci Instrum. (1958) 29:27–33.
Science Symposium and Medical Imaging Conference Record. Vol 2. IEEE doi: 10.1063/1.1715998
(1995). p. 1126–30. 35. Slomka PJ, Berman DS, Germano G. New cardiac cameras: single-
13. Beyer T, Townsend DW, Brun T, Kinahan PE, Charron M, Roddy R. photon emission CT and PET. Semin Nucl Med. (2014) 44:232–51.
A combined PET/CT scanner for clinical oncology. J Nucl Med. (2000) doi: 10.1053/j.semnuclmed.2014.04.003
41:1369–79. 36. Suzuki A, Takeuchi W, Ishitsu T, Tsuchiya K, Morimoto Y, Ueno
14. Lang TF, Hasegawa BH, Liew SC, Brown JK, Blankespoor SC, Reilly SM. Y. High-sensitivity brain SPECT system using cadmium telluride
Description of a prototype emission-transmission computed tomography (CdTe) semiconductor detector and 4-pixel matched collimator.
imaging system. J Nucl Med. (1992) 33:1881–7. Phys Med Biol. (2013) 58:7715–31. doi: 10.1088/0031-9155/58/2
15. LaCroix KJ, Tsui BMW, Hasegawa BH, Brown JK. Investigation of the use of 1/7715
X-ray CT images for attenuation compensation in SPECT. IEEE Trans Nucl 37. Tan J-W, Cai L, Meng L-J. A prototype of the MRI-compatible ultra-high
Sci. (1994) 41:2793–9. doi: 10.1109/23.340649 resolution SPECT for in vivo mice brain imaging. In: 2009 IEEE Nuclear
16. Townsend DW. Multimodality imaging of structure and function. Phys Med Science Symposium Conference Record (NSS/MIC). IEEE (2009). p. 2800–5.
Biol. (2008) 53:R1–39. doi: 10.1088/0031-9155/53/4/R01 38. Mariani G, Bruselli L, Kuwert T, Kim EE, Flotats A, Israel O. A review on the
17. Kinahan PE, Hasegawa BH, Beyer T. X-ray-based attenuation correction for clinical uses of SPECT/CT. Eur J Nucl Med Mol Imaging (2010) 37:1959–85.
positron emission tomography/computed tomography scanners. Semin Nucl doi: 10.1007/s00259-010-1390-8
Med. (2003) 33:166–79. doi: 10.1053/snuc.2003.127307 39. Bhargava P, He G, Samarghandi A, Delpassand ES. Pictorial review of
18. Ritt P, Sanders J, Kuwert T. SPECT/CT technology. Clin Transl Imaging SPECT/CT imaging applications in clinical nuclear medicine. Am J Nucl Med
(2014) 2:445–57. doi: 10.1007/s40336-014-0086-7 Mol Imaging (2012) 2:221–31.
19. Rausch I, Cal-González J, Dapra D, Gallowitsch HJ, Lind P, Beyer T. 40. Ahmadzadehfar H, Biersack H, (eds.). Clinical Applications of SPECT-CT.
Performance evaluation of the Biograph mCT Flow PET/CT system Berlin; Heidelberg: Springer-Verlag (2014).
according to the NEMA NU2-2012 standard. EJNMMI Phys. (2015) 2:26. 41. Ljungberg M, Pretorius PH. SPECT/CT: an update on technological
doi: 10.1186/s40658-015-0132-1 developments and clinical applications. Br J Radiol. (2018) 91:20160402.
20. Jakoby BW, Bercier Y, Conti M, Casey ME, Bendriem B, Townsend DW. doi: 10.1259/bjr.20160402
Physical and clinical performance of the mCT time-of-flight PET/CT 42. Townsend DW, Beyer T, Blodgett TM. PET/CT scanners: a hardware
scanner. Phys Med Biol. (2011) 56:2375–89. doi: 10.1088/0031-9155/56/8/004 approach to image fusion. Semin Nucl Med. (2003) 33:193–204.
21. Kolthammer JA, Su K-H, Grover A, Narayanan M, Jordan DW, Muzic doi: 10.1053/snuc.2003.127314
RF. Performance evaluation of the Ingenuity TF PET/CT scanner with a 43. Czernin J, Allen-Auerbach M, Schelbert HR. Improvements in cancer staging
focus on high count-rate conditions. Phys Med Biol. (2014) 59:3843–59. with PET/CT: literature-based evidence as of September 2006. J Nucl Med.
doi: 10.1088/0031-9155/59/14/3843 (2007) 48(Suppl 1):78S−88S.
22. Miller M, Zhang J, Binzel K, Griesmer J, Laurence, T, Narayanan M, et al. 44. Walrand S, Hesse M, Jamar F. Update on novel trends in PET/CT
Characterization of the vereos digital photon counting PET system. J Nucl technology and its clinical applications. Br J Radiol. (2018) 91:20160534.
Med. (2015) 56(Suppl. 3):434. doi: 10.1259/bjr.20160534

Frontiers in Physics | www.frontiersin.org 22 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

45. Herzog H, Lerche C. Advances in clinical PET/MRI instrumentation. PET 66. Huang SC, Hoffman EJ, Phelps ME, Kuhl DE. Quantitation in
Clin. (2016) 11:95–103. doi: 10.1016/j.cpet.2015.09.001 positron emission computed tomography: 2. Effects of inaccurate
46. Vandenberghe S, Marsden PK. PET-MRI: a review of challenges and attenuation correction. J Comput Assist Tomogr. (1979) 3:804–14.
solutions in the development of integrated multimodality imaging. Phys Med doi: 10.1097/00004728-197903060-00018
Biol. (2015) 60:R115–54. doi: 10.1088/0031-9155/60/4/R115 67. Fletcher JW, Kinahan PE. PET/CT Standardized uptake values (SUVs) in
47. Zaidi H, Del Guerra A. An outlook on future design of hybrid PET/MRI clinical practice and assessing response to therapy. NIH Public Access. (2010)
systems. Med Phys. (2011) 38:5667–89. doi: 10.1118/1.3633909 31:496–505. doi: 10.1053/j.sult.2010.10.001.PET/CT
48. Cabello J, Ziegler SI. Advances in PET/MR instrumentation and image 68. Wahl RL, Jacene H, Kasamon Y, Lodge MA. From RECIST to PERCIST:
reconstruction. Br J Radiol. (2018) 91:20160363. doi: 10.1259/bjr.20160363 evolving considerations for PET response criteria in solid tumors. J Nucl
49. Pichler BJ, Judenhofer MS, Catana C, Walton JH, Kneilling M, Nutt RE. Med. (2009) 50(Suppl. 1):122S−50S. doi: 10.2967/jnumed.108.057307
Performance test of an LSO-APD detector in a 7-T MRI scanner for 69. Yankeelov TE, Abramson RG, Quarles CC. Quantitative multimodality
simultaneous PET/MRI. J Nucl Med. (2006) 47:639–47. imaging in cancer research and therapy. Nat Rev Clin Oncol. (2014)
50. Spanoudaki VC, Mann AB, Otte AN, Konorov I, Torres-Espallardo 11:670–80. doi: 10.1038/nrclinonc.2014.134
I, Paul S, et al. Use of single photon counting detector arrays in 70. Kotasidis FA, Tsoumpas C, Rahmim A. Advanced kinetic modelling
combined PET/MR: characterization of LYSO-SiPM detector modules strategies: towards adoption in clinical PET imaging. Clin Transl Imaging
and comparison with a LSO-APD detector. J Instrum. (2007) 2:P12002. (2014) 2:219–37. doi: 10.1007/s40336-014-0069-8
doi: 10.1088/1748-0221/2/12/P12002 71. Lammertsma AA. Forward to the past: the case for quantitative PET imaging.
51. Judenhofer MS, Wehrl HF, Newport DF, Catana C, Siegel SB, Becker M. J Nucl Med. (2017) 58:1019–24. doi: 10.2967/jnumed.116.188029
Simultaneous PET-MRI: a new approach for functional and morphological 72. Bailey DL, Willowson KP. An evidence-based review of quantitative SPECT
imaging. Nat Med. (2008) 14:459–65. doi: 10.1038/nm1700 imaging and potential clinical applications. J Nucl Med. (2013) 54:83–9.
52. Spanoudaki VC, Levin CS. Photo-detectors for time of flight positron doi: 10.2967/jnumed.112.111476
emission tomography (ToF-PET). Sensors (2010) 10:10484–505. 73. Wernick MN, Aarsvold JN. Emission Tomography : The
doi: 10.3390/s101110484 Fundamentals of PET and SPECT. Elsevier Academic Press (2004).
53. de Jong M, Essers J, van Weerden WM. Imaging preclinical tumour doi: 10.1016/B978-0-12-744482-6.50030-2
models: improving translational power. Nat Rev Cancer. (2014) 14:481–93. 74. Casey ME, Nutt R. A Multicrystal two dimensional BGO detector system
doi: 10.1038/nrc3751 for positron emission tomography. IEEE Trans Nucl Sci. (1986) 33:460–3.
54. Meikle SR, Kench P, Kassiou M, Banati RB. Small animal SPECT and its doi: 10.1109/TNS.1986.4337143
place in the matrix of molecular imaging technologies. Phys Med Biol. (2005) 75. Joung J, Miyaoka RS, Lewellen TK. cMiCE:a high resolution animal PET
50:R45–61. doi: 10.1088/0031-9155/50/22/R01 using continuous LSO with a statistics based positioning scheme. In: IEEE
55. Lee T-S, Rittenbach A, Fernandez CG, Lopez-Longas J, Arco JM, Tsui Nuclear Science Symposium Conference Record (Cat. No.01CH37310). IEEE
BMW. Initial evaluation of a state-of-the-art commercial preclinical PET/CT (2001). p. 1137–41.
scanner. In: 2016 IEEE Nuclear Science Symposium, Medical Imaging 76. Carles M, Lerche CW, Sánchez F, Mora F, Benlloch JM. Position correction
Conference and Room-Temperature Semiconductor Detector Workshop with depth of interaction information for a small animal PET system. Nucl
(NSS/MIC/RTSD). Strasbourg: IEEE (2016). p. 1–4. Instruments Methods Phys Res Sect A Accel Spectrometers, Detect Assoc Equip.
56. van der Have F, Vastenhouw B, Ramakers RM, Branderhorst W, (2011) 648:S176–80. doi: 10.1016/j.nima.2010.11.192
Krah JO, Ji C. U-SPECT-II: an ultra-high-resolution device for 77. Lewellen TK. Recent developments in PET detector technology. Phys Med
molecular small-animal imaging. J Nucl Med. (2009) 50:599–605. Biol. (2008) 53:R287–317. doi: 10.1088/0031-9155/53/17/R01
doi: 10.2967/jnumed.108.056606 78. Madsen MT. Recent advances in SPECT imaging. J Nucl Med. (2007)
57. Goorden MC, van der Have F, Kreuger R, Ramakers RM, Vastenhouw 48:661–73. doi: 10.2967/jnumed.106.032680
B, Burbach JP. VECTor: a preclinical imaging system for simultaneous 79. Lecomte R. Novel detector technology for clinical PET. Eur J Nucl Med Mol
submillimeter SPECT and PET. J Nucl Med. (2013) 54:306–12. Imaging (2009) 36:69–85. doi: 10.1007/s00259-008-1054-0
doi: 10.2967/jnumed.112.109538 80. Peterson TE, Furenlid LR. SPECT detectors: the Anger Camera and beyond.
58. Walker MD, Goorden MC, Dinelle K, Ramakers RM, Blinder S, Phys Med Biol. (2011) 56:R145–82. doi: 10.1088/0031-9155/56/17/R01
Shirmohammad M. Performance assessment of a preclinical PET 81. Shefer E, Altman A, Behling R, Goshen R, Gregorian L, Roterman Y, et al.
scanner with pinhole collimation by comparison to a coincidence- State of the art of CT detectors and sources: a literature review. Curr Radiol
based small-animal PET scanner. J Nucl Med. (2014) 55:1368–74. Rep. (2013) 1:76–91. doi: 10.1007/s40134-012-0006-4
doi: 10.2967/jnumed.113.136663 82. Lecoq P. Development of new scintillators for medical applications. Nucl
59. Magota K, Kubo N, Kuge Y, Nishijima K-I, Zhao S, Tamaki N. Performance Instruments Methods Phys Res Sect A Accel Spectrometers Detect Assoc Equip.
characterization of the Inveon preclinical small-animal PET/SPECT/CT (2016) 809:130–9. doi: 10.1016/j.nima.2015.08.041
system for multimodality imaging. Eur J Nucl Med Mol Imaging. (2011) 83. Vaquero JJ, Sánchez JJ, Udías JM, Cal-González J, Desco M. MRI
38:742–752. doi: 10.1007/s00259-010-1683-y compatibility of position-sensitive photomultiplier depth-of-interaction
60. Larsson Åkerman L. A Technical Validation of The PET/SPECT/CT PET detectors modules for in-line multimodality preclinical studies. Nucl
(Triumph) Scanner. (2011). Available online at: https://uu.diva-portal.org/ Instruments Methods Phys Res Sect A Accel Spectrometers Detect Assoc Equip.
smash/get/diva2:407708/FULLTEXT01.pdf. (2013) 702:83–7. doi: 10.1016/j.nima.2012.08.046
61. Kuntner C, Stout D. Quantitative preclinical PET imaging: opportunities and 84. Degenhardt C, Prescher G, Frach T, Thon A, de Gruyter R, Schmitz A,
challenges. Front Phys. (2014) 2:12. doi: 10.3389/fphy.2014.00012 et al. The digital silicon photomultiplier - a novel sensor for the detection
62. Binderup T, El-Ali H, Ambrosini V, Skovgaard D, Jensen MM, Li F, et al. of scintillation light. In: IEEE Nuclear Science Symposium Conference Record.
Molecular imaging with small animal PET/CT. Curr Med Imaging Rev. Orlando, FL (2009). p. 2383–6.
(2011) 7:234–7. doi: 10.2174/157340511796411221 85. Frach T, Prescher G, Degenhardt C, De Gruyter R, Schmitz A, Ballizany
63. Goertzen AL, Bao Q, Bergeron M, Blankemeyer E, Blinder S, Cañadas M. R. The digital silicon photomultiplier - principle of operation and intrinsic
NEMA NU 4-2008 comparison of preclinical pet imaging systems. J Nucl detector performance. IEEE Nucl Sci Symp Conf Rec. (2009). p.1959–65.
Med. (2012) 53:1300–9. doi: 10.2967/jnumed.111.099382 86. Degenhardt C, Zwaans B, Frach T, de Gruyter R. Arrays of digital Silicon
64. Khalil MM, Tremoleda JL, Bayomy TB, Gsell W. Molecular SPECT imaging: Photomultipliers—intrinsic performance and application to scintillator
an overview. Int J Mol Imaging (2011) 2011:1–15. doi: 10.1155/2011/796025 readout. In: IEEE Nuclear Science Symposuim and Medical Imaging
65. Wehrl HF, Judenhofer MS, Wiehr S, Pichler BJ. Pre-clinical PET/MR: Conference. Knowville, TN: IEEE (2010). p. 1954–6.
technological advances and new perspectives in biomedical research. Eur J 87. Düppenbecker PM, Weissler B, Gebhardt P, Schug D, Wehner J, Marsden
Nucl Med Mol Imaging (2009) 36:56–68. doi: 10.1007/s00259-009-1078-0 PK. Development of an MRI-compatible digital SiPM detector stack

Frontiers in Physics | www.frontiersin.org 23 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

for simultaneous PET/MRI. Biomed Phys Eng Express. (2016) 2:15010. 107. Varela J. EndoTOFPET-US: multi-modal endoscope for ultrasound and time
doi: 10.1088/2057-1976/2/1/015010 of flight PET. In: 2014 IEEE Nuclear Science Symposium and Medical Imaging
88. Chmeissani M, Kolstein M, Gabriel Macias-Montero J, Puigdengoles Conference, NSS/MIC 2014. Seattle, WA (2016).
C, García J, Prats X. First results of a highly granulated 3D 108. Crosetto DB. A modular VME or IBM PC based data acquisition
CdTe detector module for PET. Phys Med Biol. (2018) 63:25032. system for multi-modality PET/CT scanners of different sizes and
doi: 10.1088/1361-6560/aaa44c detector types. Nucl Sci Symp Conf Rec 2000 IEEE. (2000) 2:12/78–12/97.
89. Conti M. Focus on time-of-flight PET: the benefits of improved doi: 10.1109/NSSMIC.2000.949946
time resolution. Eur J Nucl Med Mol Imaging (2011) 38:1147–57. 109. Cherry SR, Badawi RD, Karp JS, Moses WW, Price P, Jones T. Total-body
doi: 10.1007/s00259-010-1711-y imaging: transforming the role of positron emission tomography. Sci Transl
90. Surti S, Karp JS. Advances in time-of-flight PET. Phys Med. (2016) 32:12–22. Med. (2017) 9:eaaf6169. doi: 10.1126/scitranslmed.aaf6169
doi: 10.1016/j.ejmp.2015.12.007 110. Cherry SR, Jones T, Karp JS, Qi J, Moses W, Badawi R. Total-body PET:
91. Vandenberghe S, Mikhaylova E, D’Hoe E, Mollet P, Karp JS. Recent maximizing sensitivity to create new opportunities for clinical research
developments in time-of-flight PET. EJNMMI Phys. (2016) 3:3. and patient care. J Nucl Med. (2018) 59:3–12. doi: 10.2967/jnumed.116.
doi: 10.1186/s40658-016-0138-3 184028
92. Lecoq P. Pushing the limits in time-of-flight PET imaging. IEEE Trans Radiat 111. Eriksson L, Townsend D, Conti M, Eriksson M, Rothfuss H, Schmand M,
Plasma Med Sci. (2017) 1:473–85. doi: 10.1109/TRPMS.2017.2756674 et al. An investigation of sensitivity limits in PET scanners. Nucl Instrum
93. Gundacker S, Acerbi F, Auffray E, Ferri A, Gola A, Nemallapudi MV, et Methods Phys Res Sect A Accel Spectrometers Detect Assoc Equip. (2007)
al. State of the art timing in TOF-PET detectors with LuAG, GAGG and 580:836–42. doi: 10.1016/j.nima.2007.06.112
L(Y)SO scintillators of various sizes coupled to FBK-SiPMs. J Instrum. (2016) 112. Zhang X, Zhou J, Cherry SR, Badawi RD, Qi J. Quantitative
11:P08008. doi: 10.1088/1748-0221/11/08/P08008 image reconstruction for total-body PET imaging using the 2-
94. Nemallapudi MV, Gundacker S, Lecoq P, Auffray E, Ferri A, Gola A. meter long EXPLORER scanner. Phys Med Biol. (2017) 62:2465–85.
Sub-100 ps coincidence time resolution for positron emission tomography doi: 10.1088/1361-6560/aa5e46
with LSO:Ce codoped with Ca. Phys Med Biol. (2015) 60:4635–49. 113. Karakatsanis NA, Garibotto V, Rager O, Zaidi H. Continuous bed motion Vs.
doi: 10.1088/0031-9155/60/12/4635 step-and-shoot acquisition on clinical whole-body dynamic and parametric
95. Dahlbom M, Yu D-C, Cherry SR, Chatziioannou A, Hoffman EJ. Methods PET imaging. In: 2015 IEEE Nuclear Science Symposium and Medical Imaging
for improving image quality in whole body PET scanning. IEEE Trans Nucl Conference, NSS/MIC 2015. San Diego, CA (2016).
Sci. (1992) 39:1079–83. doi: 10.1109/23.159763 114. Karakatsanis NA, Lodge M A, Tahari AK, Zhou Y, Wahl RL, Rahmim
96. Dahlbom M, Reed J, Young J. Implementation of true continuous bed A. Dynamic whole-body PET parametric imaging: I. Concept, acquisition
motion in 2-D and 3-D whole-body PET scanning. IEEE Trans Nucl Sci. protocol optimization and clinical application. Phys Med Biol. (2013)
(2001) 48:1465–9. doi: 10.1109/23.958381 58:7391–418. doi: 10.1088/0031-9155/58/20/7391
97. Brasse D, Newport D, Carney JP, Yap J, Reynolds C, Reed J, et al. Continuous 115. Randall I. Total Body PET Scanner Targets Paediatric Patients.
bed motion acquisition on a whole body combined PET/CT system. In: EEE MedicalPhysicsWeb. Available online at http://medicalphysicsweb.org/
Nuclear Science Symposium Conference Record. Vol 2. Norfolk, VA: IEEE cws/article/research/69657 (Accessed November 14, 2017).
(2002). p. 951–5. 116. Martins PM, Crespo P, Couceiro M, Ferreira NC, Marques RF, Seco J, et
98. Townsend DW, Reed J, Newport DF, Carney JPJ, Tolbert SH, Newby D et al. al. Fast full-body reconstruction for a functional human RPC-PET imaging
Continuous bed motion acquisition for an LSO PET/CT scanner. In: IEEE system using list-mode simulated data and its applicability to radiation
Symposium Conference Record Nuclear Science. Vol 4. Rome: IEEE (2004). oncology and radiology. arXiv. 2017:1706.07075.
p. 2383–7. 117. Moser E, Stahlberg F, Ladd ME, Trattnig S. 7-T MR-from research to clinical
99. Karakatsanis NA, Zhou Y, Lodge MA, Casey ME, Wahl RL, Zaidi applications? NMR Biomed. (2012) 25:695–716. doi: 10.1002/nbm.1794
H. Generalized whole-body Patlak parametric imaging for enhanced 118. Moser E, Stadlbauer A, Windischberger C, Quick HH, Ladd ME. Magnetic
quantification in clinical PET. Phys Med Biol. (2015) 60:8643–73. resonance imaging methodology. Eur J Nucl Med Mol Imaging (2009)
doi: 10.1088/0031-9155/60/22/8643 36:30–41. doi: 10.1007/s00259-008-0938-3
100. Wienhard K, Schmand M, Casey ME, Baker K, Bao J, Eriksson L, et al. 119. Le Bihan D, Mangin JF, Poupon C, Clark CA, Pappata S, Molko N. Diffusion
The ECAT HRRT: performance and first clinical application of the new tensor imaging: concepts and applications. J Magn Reson Imaging (2001)
high resolution research tomograph. IEEE Trans Nucl Sci. (2002) 49:104. 13:534–46. doi: 10.1002/jmri.1076
doi: 10.1109/TNS.2002.998689 120. Sander CY, Keil B, Chonde DB, Rosen BR, Catana C, Wald LL. A 31-channel
101. Majewski S, Proffitt J, Brefczynski-Lewis J, Stolin A, Weisenberger AG, Xi MR brain array coil compatible with positron emission tomography. Magn
W, et al. HelmetPET: a silicon photomultiplier based wearable brain imager. Reson Med. (2015) 73:2363–75. doi: 10.1002/mrm.25335
In: IEEE Nuclear Science Symposium Conference Record. Valencia (2011). p. 121. Pruessmann KP, Weiger M, Scheidegger MB, Boesiger P. SENSE: sensitivity
4030–4. encoding for fast MRI. Magn Reson Med. (1999) 42:952–62.
102. Yamamoto S, Honda M, Oohashi T, Shimizu K, Senda M. Development of a 122. Griswold MA, Jakob PM, Heidemann RM, Nittka M, Jellus V, Wang J.
brain PET system, PET-Hat: a wearable PET system for brain research. IEEE Generalized autocalibrating partially parallel acquisitions (GRAPPA). Magn
Trans Nucl Sci. (2011) 58(3 Pt 1):668–73. doi: 10.1109/TNS.2011.2105502 Reson Med. (2002) 47:1202–10. doi: 10.1002/mrm.10171
103. Tashima H, Yamaya T. Proposed helmet PET geometries with add- 123. Larkman DJ, Nunes RG. Parallel magnetic resonance imaging. Phys Med Biol.
on detectors for high sensitivity brain imaging. Phys Med Biol. (2016) (2007) 52:R15–55. doi: 10.1088/0031-9155/52/7/R01
61:7205–20. doi: 10.1088/0031-9155/61/19/7205 124. Alessio A, Kinahan P. PET image reconstruction. Nucl Med. (2006) 1:1–22.
104. Kolb A, Wehrl HF, Hofmann M, Judenhofer MS, Eriksson L, Ladebeck R. doi:10.1088/0031-9155/54/12/007.Iterative
Technical performance evaluation of a human brain PET/MRI system. Eur 125. Hsieh J, Nett B, Yu Z, Sauer K, Thibault J-B, Bouman CA. Recent
Radiol. (2012) 22:1776–88. doi: 10.1007/s00330-012-2415-4 advances in CT image reconstruction. Curr Radiol Rep. (2013) 1:39–51.
105. González AJ, Majewski S, Sánchez F, Aussenhofer S, Aguilar A, Conde P, et al. doi: 10.1007/s40134-012-0003-7
The MINDView brain PET detector, feasibility study based on SiPM arrays. 126. Geyer LL, Schoepf UJ, Meinel FG, Nance JW, Bastarrika G, Leipsic JA.
Nucl Instruments Methods Phys Res Sect A Accel Spectrometers Detect Assoc State of the art: iterative CT reconstruction techniques. Radiology. (2015)
Equip. (2016) 818:82–90. doi: 10.1016/j.nima.2016.02.046 276:339–57. doi: 10.1148/radiol.2015132766
106. Ravindranath B, Junnarkar SS, Purschke ML, Maramraju S, Hong X, Tomasi 127. Herman GT. Fundamentals of Computerized Tomography. London: Springer
D, et al. Results from prototype II of the BNL simultaneous PET-MRI London (2009).
dedicated breast scanner. In: IEEE Nuclear Science Symposium Conference 128. Feldkamp LA, Davis LC, Kress JW. Practical cone-beam algorithm. J Opt Soc
Record. Orlando, FL (2009). p. 3315–7. doi: 10.1109/NSSMIC.2009.5401742 Am A. (1984) 1:612–9. doi: 10.1364/JOSAA.1.000612

Frontiers in Physics | www.frontiersin.org 24 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

129. Katsevich A. A general scheme for constructing inversion algorithms 151. Ogawa K, Harata Y, Ichihara T, Kubo A, Hashimoto S. A practical method
for cone beam CT. Int J Math Math Sci. (2003) 2003:1305–21. for position-dependent Compton-scatter correction in single photon
doi: 10.1155/S0161171203209315 emission CT. IEEE Trans Med Imaging. (1991) 10:408–12. doi: 10.1109/42.
130. Zhuang T, Leng S, Nett BE, Chen G-H. Fan-beam and cone-beam 97591
image reconstruction via filtering the backprojection image of 152. Zaidi H, Koral KF. Scatter modelling and compensation in emission
differentiated projection data. Phys Med Biol. (2004) 49:5489–503. tomography. Eur J Nucl Med Mol Imaging (2004) 31:761–82. doi: 10.1007/
doi: 10.1088/0031-9155/49/24/007 s00259-004-1495-z
131. Pack JD, Noo F, Clackdoyle R. Cone-beam reconstruction using the 153. Peterson M, Gustafsson J, Ljungberg M. Monte Carlo-based quantitative
backprojection of locally filtered projections. IEEE Trans Med Imaging pinhole SPECT reconstruction using a ray-tracing back-projector. EJNMMI
(2005) 24:70–85. doi: 10.1109/TMI.2004.837794 Phys. (2017) 4:32. doi: 10.1186/s40658-017-0198-z
132. Grimmer R, Oelhafen M, Elstrøm U, Kachelrieß M. Cone-beam CT image 154. Elschot M, Lam MGEH, van den Bosch MAAJ, Viergever MA, de Jong
reconstruction with extended z range. Med Phys. (2009) 36:3363–70. HWAM. Quantitative Monte Carlo-based 90Y SPECT reconstruction. J Nucl
doi: 10.1118/1.3148560 Med. (2013) 54:1557–63. doi: 10.2967/jnumed.112.119131
133. Katsevich A. Theoretically exact filtered backprojection-type inversion 155. Bailey DL, Townsend DW. Positron Emission Tomography: Basic Sciences.
algorithm for spiral CT. SIAM J Appl Math. (2002) 62:2012–26. Basel: Springer (2005).
doi: 10.1137/S0036139901387186 156. Badawi RD, Marsden PK. Developments in component-based
134. Tang X, Hsieh J, Nilsen RA, Dutta S, Samsonov D, Hagiwara A. A normalization for 3D PET. Phys Med Biol. (1999) 44:571–94.
three-dimensional-weighted cone beam filtered backprojection (CB-FBP) doi: 10.1088/0031-9155/44/2/020
algorithm for image reconstruction in volumetric CT—helical scanning. Phys 157. Guobao W, Jinyi Q. Penalized likelihood PET image reconstruction using
Med Biol. (2006) 51:855–74. doi: 10.1088/0031-9155/51/4/007 patch-based edge-preserving regularization. IEEE Trans Med Imaging.
135. Beister M, Kolditz D, Kalender WA. Iterative reconstruction methods (2012) 31:2194–204. doi: 10.1109/TMI.2012.2211378
in X-ray CT. Phys Medica. (2012) 28:94–108. doi: 10.1016/j.ejmp.2012. 158. Green PJ. Bayesian reconstructions from emission tomography data using
01.003 a modified EM algorithm. IEEE Trans Med Imaging (1990) 9:84–93.
136. Gordon R, Bender R, Herman GT. Algebraic reconstruction techniques doi: 10.1109/42.52985
(ART) for three-dimensional electron microscopy and x-ray photography. 159. Alessio AM, Kinahan PE. Improved quantitation for PET/CT image
J Theor Biol. (1970) 29:471–81. doi: 10.1016/0022-5193(70)90109-8 reconstruction with system modeling and anatomical priors. Med Phys.
137. Andersen AH, Kak AC. Simultaneous algebraic reconstruction technique (2006) 33:4095–103. doi: 10.1118/1.2358198
(SART): a superior implementation of the art algorithm. Ultrason Imaging 160. Chun SY. The use of anatomical information for molecular image
(1984) 6:81–94. doi: 10.1177/016173468400600107 reconstruction algorithms: attenuation/scatter correction, motion
138. Gilbert P. Iterative methods for the three-dimensional reconstruction compensation, and noise reduction. Nucl Med Mol Imaging (2010).
of an object from projections. J Theor Biol. (1972) 36:105–17. (2016) 50:13–23. doi: 10.1007/s13139-016-0399-8
doi: 10.1016/0022-5193(72)90180-4 161. Bai B, Li Q, Leahy RM. Magnetic resonance-guided positron emission
139. Shepp LA, Vardi Y. Maximum likelihood reconstruction for tomography image reconstruction. Semin Nucl Med. (2013) 43:30–44.
emission tomography. IEEE Trans Med Imaging (1982) 1:113–22. doi: 10.1053/j.semnuclmed.2012.08.006
doi: 10.1109/TMI.1982.4307558 162. Hebert T, Leahy R. A generalized EM algorithm for 3-D Bayesian
140. Hudson HM, Larkin RS. Accelerated image reconstruction using ordered reconstruction from Poisson data using Gibbs priors. IEEE Trans Med
subsets of projection data. IEEE Trans Med Imaging (1994) 13:601–9. Imaging. (1989) 8:194–202. doi: 10.1109/42.24868
doi: 10.1109/42.363108 163. Ardekani BA, Braun M, Hutton BF, Kanno I, Iida H. Minimum cross-entropy
141. Thibault J-B, Sauer KD, Bouman CA, Hsieh J. A three-dimensional statistical reconstruction of PET images using prior anatomical information. Phys Med
approach to improved image quality for multislice helical CT. Med Phys. Biol. (1996) 41:2497–517. doi: 10.1088/0031-9155/41/11/018
(2007) 34:4526–44. doi: 10.1118/1.2789499 164. Soret M, Bacharach SL, Buvat I. Partial-volume effect in PET tumor imaging.
142. Xu F, Xu W, Jones M, Keszthelyi B, Sedat J, Agard D. On the J Nucl Med. (2007) 48:932–45. doi: 10.2967/jnumed.106.035774
efficiency of iterative ordered subset reconstruction algorithms for 165. Rousset O, Rahmim A, Alavi A, Zaidi H. Partial volume correction strategies
acceleration on GPUs. Comput Methods Programs Biomed. (2010) 98:261–70. in PET. PET Clin. (2007) 2:235–49. doi: 10.1016/j.cpet.2007.10.005
doi: 10.1016/j.cmpb.2009.09.003 166. Erlandsson K, Buvat I, Pretorius PH, Thomas BA, Hutton BF. A review of
143. Liang Z-P, Lauterbur PC. Principles of Magnetic Resonance Imaging : A partial volume correction techniques for emission tomography and their
Signal Processing Perspective. New York, NY: SPIE Optical Engineering Press applications in neurology, cardiology and oncology. Phys Med Biol. (2012)
(2000). 57:R119–59. doi: 10.1088/0031-9155/57/21/R119
144. Zaidi H. Quantitative Analysis in Nuclear Medicine Imaging. Basel: Springer 167. Erlandsson K, Dickson J, Arridge S, Atkinson D, Ourselin S, Hutton BF. MR
(2006). imaging–guided partial volume correction of PET data in PET/MR imaging.
145. Willowson K, Bailey DL, Baldock C. Quantitative SPECT reconstruction PET Clin. (2015) 11:161–77. doi: 10.1016/j.cpet.2015.09.002
using CT-derived corrections. Phys Med Biol. (2008) 53:3099–12. 168. Hoffman EJ, Huang SC, Phelps ME. Quantitation in positron emission
doi: 10.1088/0031-9155/53/12/002 computed tomography: 1. Effect of object size. J Comput Assist Tomogr.
146. Kinahan PE, Townsend DW, Beyer T, Sashin D. Attenuation correction (1979) 3:299–308. doi: 10.1097/00004728-197906000-00001
for a combined 3D PET/CT scanner. Med Phys. (1998) 25:2046–53. 169. Rousset OG, Ma Y, Evans AC. Correction for partial volume effects in PET:
doi: 10.1118/1.598392 principle and validation. J Nucl Med. (1998) 39:904–11.
147. Watson CC. New, faster, image-based scatter correction for 3D PET. IEEE 170. Müller-Gärtner HW, Links JM, Prince JL, Bryan RN, McVeigh E,
Trans Nucl Sci. (2000) 47:1587–94. doi: 10.1109/23.873020 Leal JP. Measurement of radiotracer concentration in brain gray
148. Watson CC, Casey ME, Michel C, Bendriem B. Advances in scatter matter using positron emission tomography: MRI-based correction for
correction for 3D PET/CT. IEEE Symp Conf Rec Nucl Sci 2004. (2004) 5:3–7. partial volume effects. J Cereb Blood Flow Metab. (1992) 12:571–83.
doi: 10.1109/NSSMIC.2004.1466317 doi: 10.1038/jcbfm.1992.81
149. Patton JA, Turkington TG. SPECT/CT physical principles and 171. Teo BK, Seo Y, Bacharach SL, Carrasquillo JA, Libutti SK, Shukla H. Partial-
attenuation correction. J Nucl Med Technol. (2008) 36:1–10. volume correction in PET: validation of an iterative postreconstruction
doi: 10.2967/jnmt.107.046839 method with phantom and patient data. J Nucl Med. (2007) 48:802–10.
150. Jaszczak RJ, Greer KL, Floyd CE, Harris CC, Coleman RE. Improved SPECT doi: 10.2967/jnumed.106.035576
quantification using compensation for scattered photons. J Nucl Med. (1984) 172. Thomas BA, Erlandsson K, Modat M, Thurfjell L, Vandenberghe R, Ourselin
25:893–900. S. The importance of appropriate partial volume correction for PET

Frontiers in Physics | www.frontiersin.org 25 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

quantification in Alzheimer’s disease. Eur J Nucl Med Mol Imaging. (2011) 192. Huang C, Petibon Y, Ouyang J, Reese TG, Ahlman MA, Bluemke DA.
38:1104–19. Accelerated acquisition of tagged MRI for cardiac motion correction in
173. Moore SC, Southekal S, Park MA, McQuaid SJ, Kijewski MF, Muller SP. simultaneous PET-MR: phantom and patient studies. Med Phys. (2015)
Improved regional activity quantitation in nuclear medicine using a new 42:1087–97. doi: 10.1118/1.4906247
approach to correct for tissue partial volume and spillover effects. IEEE Trans 193. Livieratos L, Stegger L, Bloomfield PM, Schafers K, Bailey DL, Camici
Med Imaging. (2012) 31:405–16. doi: 10.1109/TMI.2011.2169981 PG. Rigid-body transformation of list-mode projection data for respiratory
174. Southekal S, McQuaid SJ, Kijewski MF, Moore SC. Evaluation of a method for motion correction in cardiac PET. Phys Med Biol. (2005) 50:3313–22.
projection-based tissue-activity estimation within small volumes of interest. doi: 10.1088/0031-9155/50/14/008
Phys Med Biol. (2012) 57:685–701. doi: 10.1088/0031-9155/57/3/685 194. Rahmim A, Tang J, Zaidi H. Four-Dimensional image reconstruction
175. Cal-González J, Moore SC, Park MA, Herraiz JL, Vaquero JJ, strategies in cardiac-gated and respiratory-gated PET imaging. PET Clin.
Desco M. Improved quantification for local regions of interest (2013) 8:51–67. doi: 10.1016/j.cpet.2012.10.005
in preclinical PET imaging. Phys Med Biol. (2015) 60:7127–49. 195. Picard Y, Thompson CJ. Motion correction of PET images using
doi: 10.1088/0031-9155/60/18/7127 multiple acquisition frames. IEEE Trans Med Imaging (1997) 16:137–44.
176. Cal-Gonzalez J, Li X, Heber D, Rausch I, Moore SC, Schäfers K. doi: 10.1109/42.563659
Partial volume correction for improved PET quantification in 18F- 196. Feng T, Wang J, Fung G, Tsui B. Non-rigid dual respiratory and
NaF imaging of atherosclerotic plaques. J Nucl Cardiol. (2017). cardiac motion correction methods after, during, and before image
doi: 10.1007/s12350-017-0778-2. [Epub ahead of print]. reconstruction for 4D cardiac PET. Phys Med Biol. (2016) 61:151–68.
177. Panin VY, Kehren F, Michel C, Casey M. Fully 3-D PET reconstruction with doi: 10.1088/0031-9155/61/1/151
system matrix derived from point source measurements. IEEE Trans Med 197. Polycarpou I, Tsoumpas C, Marsden PK. Analysis and comparison of two
Imaging (2006) 25:907–21. doi: 10.1109/TMI.2006.876171 methods for motion correction in PET imaging. Med Phys. (2012) 39:6474.
178. Cysouw MCF, Kramer GM, Schoonmade LJ, Boellaard R, de Vet HCW, doi: 10.1118/1.4754586
Hoekstra OS. Impact of partial-volume correction in oncological PET 198. Fieseler M, Gigengack F, Jiang X, Schäfers KP. Motion correction of whole-
studies: a systematic review and meta-analysis. Eur J Nucl Med Mol Imaging body PET data with a joint PET-MRI registration functional. Biomed Eng
(2017) 44:2105–16. doi: 10.1007/s00259-017-3775-4 Online. (2014) 13(Suppl. 1):S2. doi: 10.1186/1475-925X-13-S1-S2
179. Nehmeh SA, Erdi YE. Respiratory motion in positron emission 199. Petibon Y, Guehl NJ, Reese TG, Ebrahimi B, Normandin MD, Shoup TM.
tomography/computed tomography: a review. Semin Nucl Med. (2008) Impact of motion and partial volume effects correction on PET myocardial
38:167–76. doi: 10.1053/j.semnuclmed.2008.01.002 perfusion imaging using simultaneous PET-MR. Phys Med Biol. (2017)
180. Gould KL, Pan T, Loghin C, Johnson NP, Guha A, Sdringola S. Frequent 62:326–43. doi: 10.1088/1361-6560/aa5087
diagnostic errors in cardiac PET/CT due to misregistration of CT attenuation 200. Cal-González J, Tsoumpas C, Lassen ML, Rasul S, Koller L, Hacker M.
and emission PET images: a definitive analysis of causes, consequences, and Impact of motion compensation and partial volume correction for 18 F-
corrections. J Nucl Med. (2007) 48:1112–21. doi: 10.2967/jnumed.107.039792 NaF PET/CT imaging of coronary plaque. Phys Med Biol. (2017) 63:015005.
181. Rausch I, Quick HH, Cal-Gonzalez J, Sattler B, Boellaard R, Beyer T. doi: 10.1088/1361-6560/aa97c8
Technical and instrumentational foundations of PET/MRI. Eur J Radiol. 201. Keller SH, Hansen C, Hansen C, Andersen FL, Ladefoged C, Svarer C.
(2017) 94:A3–13. doi: 10.1016/j.ejrad.2017.04.004 Sparsely sampled MR navigators as a practical tool for quality control and
182. Kesner AL, Chung JH, Lind KE, Kwak JJ, Lynch D, Burckhardt D. Validation correction of head motion in simultaneous PET/MR. EJNMMI Phys. (2015)
of software gating: a practical technology for respiratory motion correction 1(Suppl. 1):A36. doi: 10.1186/2197-7364-1-S1-A36
in PET. Radiology. (2016) 281:239–48. doi: 10.1148/radiol.2016152105 202. Rahmim A, Zaidi H. PET versus SPECT: strengths, limitations
183. Lassen ML, Rasmussen T, Christensen TE, Kjær A, Hasbak P. Respiratory and challenges. Nucl Med Commun. (2008) 29:193–207.
gating in cardiac PET: effects of adenosine and dipyridamole. J Nucl Cardiol. doi: 10.1097/MNM.0b013e3282f3a515
(2016) 24:1941–9. doi: 10.1007/s12350-016-0631-z 203. Sanders JC, Ritt P, Kuwert T, Vija AH, Maier AK. Fully automated
184. Olesen OV, Paulsen RR, Højgaard L, Roed B, Larsen R. Motion tracking for data-driven respiratory signal extraction from SPECT images using
medical imaging: a nonvisible structured light tracking approach. IEEE Trans laplacian eigenmaps. IEEE Trans Med Imaging. (2016) 35:2425–35.
Med Imaging. (2012) 31:79–87. doi: 10.1109/TMI.2011.2165157 doi: 10.1109/TMI.2016.2576899
185. Kesner AL, Schleyer PJ, Büther F, Walter MA, Schäfers KP, Koo PJ. On 204. Thielemans K, Schleyer P, Marsden PK, Manjeshwar RM, Wollenweber S,
transcending the impasse of respiratory motion correction applications Ganin A. Comparison of different methods for data-driven respiratory gating
in routine clinical imaging – a consideration of a fully automated of PET data. In: 2013 IEEE Nuclear Science Symposium and Medical Imaging
data driven motion control framework. EJNMMI Phys. (2014) 1:8. Conference (2013 NSS/MIC). Seoul: IEEE (2013), p. 1–4.
doi: 10.1186/2197-7364-1-8 205. Beyer T, Moser E. MR/PET or PET/MRI: does it matter? Magn Reson Mater
186. Grimm R, Fürst S, Souvatzoglou M, Forman C, Hutter J, Dregely I. Phys Biol Med. (2013) 26:1–4. doi:10.1007/s10334-012-0365-0
Self-gated MRI motion modeling for respiratory motion compensation 206. Träber F, Block W, Layer G, Bräucker G, Gieseke J, Kretzer S, et al.
in integrated PET/MRI. Med Image Anal. (2014) 19:110–20. Determination of H relaxation times of water in human bone marrow
doi: 10.1016/j.media.2014.08.003 by fat-suppressed turbo spin echo in comparison to MR spectroscopic
187. Munoz C, Kolbitsch C, Reader AJ, Marsden P, Schaeffter T, Prieto C. MR- methods. J Magn Reson Imaging (1996) 6:541–8. doi: 10.1002/jmri.18800
based cardiac and respiratory motion-compensation techniques for PET-MR 60318
imaging. PET Clin. (2016) 11:179–91. doi: 10.1016/j.cpet.2015.09.004 207. Samarin A, Burger C, Wollenweber SD, Crook DW, Burger IA, Schmid
188. Kesner AL. The relevance of data driven motion correction DT. PET/MR imaging of bone lesions - Implications for PET quantification
in diagnostic PET. Eur J Nucl Med Mol Imaging 44:2326–7. from imperfect attenuation correction. Eur J Nucl Med Mol Imaging. (2012)
doi: 10.1007/s00259-017-3794-1 39:1154–60. doi: 10.1007/s00259-012-2113-0
189. Kolbitsch C, Ahlman MA, Davies-Venn C, Evers R, Hansen M, Peressutti D. 208. Koesters T, Friedman KP, Fenchel M, Zhan Y, Hermosillo G, Babb J. Dixon
Cardiac and respiratory motion correction for simultaneous cardiac PET - sequence with superimposed model-based bone compartment provides
MR. J Nucl Med. (2017) 58:846–52. doi: 10.2967/jnumed.115.171728 highly accurate PET/MR attenuation correction of the brain. J Nucl Med.
190. Munoz C, Neji R, Cruz G, Mallia A, Jeljeli S, Reader AJ. Motion-corrected (2016) 57:918–24. doi: 10.2967/jnumed.115.166967
simultaneous cardiac positron emission tomography and coronary MR 209. Anazodo UC, Thiessen JD, Ssali T, Mandel J, Günther M, Butler J. Feasibility
angiography with high acquisition efficiency. Magn Reson Med. (2017) of simultaneous whole-brain imaging on an integrated PET-MRI system
79:339–50. doi: 10.1002/mrm.26690 using an enhanced 2-point Dixon attenuation correction method. Front
191. Gillman A, Rose S, Smith J, Thomas P, Dowson N. PET motion correction Neurosci. (2015) 8:434. doi: 10.3389/fnins.2014.00434
in context of integrated PET/MR: current techniques, limitations, and future 210. Izquierdo-Garcia D, Hansen AE, Förster S, Benoit D, Schachoff
projections. Med Phys. (2017). 44:e430–45. doi: 10.1002/mp.12577. S, Fürst S. An SPM8-based approach for attenuation correction

Frontiers in Physics | www.frontiersin.org 26 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

combining segmentation and nonrigid template formation: application 227. Heußer T, Rank CM, Berker Y, Freitag MT, Kachelrieß M. MLAA-based
to simultaneous PET/MR brain imaging. J Nucl Med. (2014) 55:1825–30. attenuation correction of flexible hardware components in hybrid PET/MR
doi: 10.2967/jnumed.113.136341 imaging. EJNMMI Phys. (2017) 4:12. doi: 10.1186/s40658-017-0177-4
211. Burgos N, Cardoso MJ, Thielemans K, Modat M, Pedemonte S, 228. Sari H, Erlandsson K, Law I, Larsson HB, Ourselin S, Arridge S. Estimation
Dickson J. Attenuation correction synthesis for hybrid PET-MR scanners: of an image derived input function with MR-defined carotid arteries in FDG-
application to brain studies. IEEE Trans Med Imaging (2014) 33:2332–41. PET human studies using a novel partial volume correction method. J Cereb
doi: 10.1109/TMI.2014.2340135 Blood Flow Metab. (2017) 37:1398–409. doi: 10.1177/0271678X16656197
212. Merida I, Costes N, Heckemann R, Hammers A. Pseudo-CT 229. Iida H, Rhodes CG, de Silva R, Yamamoto Y, Araujo LI, Maseri A. Myocardial
generation in brain MR-PET attenuation correction: comparison of tissue fraction–correction for partial volume effects and measure of tissue
several multi-atlas methods. EJNMMI Phys. (2015) 2(Suppl. 1):A29. viability. J Nucl Med. (1991) 32:2169–2175.
doi: 10.1186/2197-7364-2-S1-A29 230. Iida H, Rhodes CG, de Silva R, Araujo LI, Bloomfield PM, Lammertsma AA.
213. Benoit D, Ladefoged CN, Rezaei A, Keller SH, Andersen FL, Højgaard L. Use of the left ventricular time-activity curve as a noninvasive input function
Optimized MLAA for quantitative non-TOF PET/MR of the brain. Phys Med in dynamic oxygen-15-water positron emission tomography. J Nucl Med.
Biol. (2016) 61:8854–74. doi: 10.1088/1361-6560/61/24/8854 (1992) 33:1669–77.
214. Cabello J, Lukas M, Forster S, Pyka T, Nekolla SG, Ziegler SI. 231. Zanotti-Fregonara P, Chen K, Liow J-S, Fujita M, Innis RB. Image-derived
MR-Based attenuation correction using ultrashort-echo-time pulse input function for brain PET studies: many challenges and few opportunities.
sequences in dementia patients. J Nucl Med. (2015) 56:423–9. J Cereb Blood Flow Metab. (2011) 31:1986–98. doi: 10.1038/jcbfm.2011.107
doi: 10.2967/jnumed.114.146308 232. Zanotti-Fregonara P, Liow JS, Comtat C, Zoghbi SS, Zhang Y, Pike VW.
215. Juttukonda MR, Mersereau BG, Chen Y, Su Y, Rubin BG, Benzinger TLS. Image-derived input function in PET brain studies. Nucl Med Commun.
MR-based attenuation correction for PET/MRI neurological studies (2012) 33:982–9. doi: 10.1097/MNM.0b013e328356185c
with continuous-valued attenuation coefficients for bone through 233. Mourik JEM, Lubberink M, Lammertsma AA, Boellaard R. Image derived
a conversion from R2∗ to CT-Hounsfield units. Neuroimage. (2015) input functions: effects of motion on tracer kinetic analyses. Mol Imaging
112:160–8. doi: 10.1016/j.neuroimage.2015.03.009 Biol. (2011) 13:25–31. doi: 10.1007/s11307-010-0301-5
216. Ladefoged CN, Benoit D, Law I, Holm S, Kjær A, Højgaard L. Region specific 234. Mourik JE, Lubberink M, Schuitemaker A, Tolboom N, van Berckel BN,
optimization of continuous linear attenuation coefficients based on UTE Lammertsma AA. Image-derived input functions for PET brain studies. Eur
(RESOLUTE): application to PET/MR brain imaging. Phys Med Biol. (2015) J Nucl Med Mol Imaging (2009) 36:463–71. doi: 10.1007/s00259-008-0986-8
60:8047–65. doi: 10.1088/0031-9155/60/20/8047 235. Mourik JEM, Lubberink M, Klumpers UMH, Comans EF, Lammertsma
217. Ladefoged CN, Law I, Anazodo U, St Lawrence K, Izquierdo-Garcia D, AA, Boellaard R. Partial volume corrected image derived input
Catana C. A multi-centre evaluation of eleven clinically feasible brain functions for dynamic PET brain studies: methodology and
PET/MRI attenuation correction techniques using a large cohort of patients. validation for [11C]flumazenil. Neuroimage. (2008) 39:1041–50.
Neuroimage (2017) 147:346–59. doi: 10.1016/j.neuroimage.2016.12.010 doi: 10.1016/j.neuroimage.2007.10.022
218. Martinez-Möller A, Souvatzoglou M, Delso G, Bundschuh RA, Chefd’hotel 236. Zhou S, Chen K, Reiman EM, Li D, Shan B. A method for generating image-
C, Ziegler SI. Tissue classification as a potential approach for attenuation derived input function in quantitative 18F-FDG PET study based on the
correction in whole-body PET/MRI: evaluation with PET/CT data. J Nucl monotonicity of the input and output function curve. Nucl Med Commun.
Med. (2009) 50:520–6. doi: 10.2967/jnumed.108.054726 (2012) 33:362–70. doi: 10.1097/MNM.0b013e32834f262e
219. Rausch I, Rust P, DiFranco MD, Lassen M, Stadlbauer A, Mayerhoefer ME. 237. Wahl LM, Asselin MC, Nahmias C. Regions of interest in the venous sinuses
Reproducibility of MRI dixon-based attenuation correction in combined as input functions for quantitative PET. J Nucl Med. (1999) 40:1666–75.
PET/MR with applications for lean body mass estimation. J Nucl Med. (2016) 238. Guo H, Renaut R, Chen K, Reiman E. Clustering huge data
57:1096–101. doi: 10.2967/jnumed.115.168294 sets for parametric PET imaging. BioSystems. (2003) 71:81–92.
220. Sekine T, Buck A, Delso G, Ter Voert EE, Huellner M, Veit-Haibach P. doi: 10.1016/S0303-2647(03)00112-6
Evaluation of atlas-based attenuation correction for integrated PET/MR 239. Fang Y-H, Kao T, Liu R-S, Wu L-C. Estimating the input function non-
in human brain: application of a head atlas and comparison to invasively for FDG-PET quantification with multiple linear regression
true CT-based attenuation correction. J Nucl Med. (2016) 57:215–220. analysis: simulation and verification with in vivo data. Eur J Nucl Med Mol
doi: 10.2967/jnumed.115.159228 Imaging (2004) 31:692–702. doi: 10.1007/s00259-003-1412-x
221. Paulus DH, Quick HH, Geppert C, Fenchel M, Zhan Y, Hermosillo G. 240. Liptrot M, Adams KH, Martiny L, Pinborg LH, Lonsdale MN, Olsen
Whole-body PET/MR imaging: quantitative evaluation of a novel model- NV. Cluster analysis in kinetic modelling of the brain: a noninvasive
based MR attenuation correction method including bone. J Nucl Med. (2015) alternative to arterial sampling. Neuroimage (2004) 21:483–93.
56:1061–6. doi: 10.2967/jnumed.115.156000 doi: 10.1016/j.neuroimage.2003.09.058
222. Rausch I, Rischka L, Ladefoged CN, Furtner J, Fenchel M, Hahn 241. Naganawa M, Kimura Y, Ishii K, Oda K, Ishiwata K, Matani A. Extraction
A. PET/MRI for oncologic brain imaging: a comparison of standard of a plasma time-activity curve from dynamic brain PET images based on
MR-based attenuation corrections with a model-based approach for independent component analysis. IEEE Trans Biomed Eng. (2005) 52:201–10.
the siemens mMR PET/MR system. J Nucl Med. (2017) 58:1519–25. doi: 10.1109/TBME.2004.840193
doi: 10.2967/jnumed.116.186148 242. Su Y, Shoghi KI. Single-input–dual-output modeling of image-based
223. Leynes AP, Yang J, Wiesinger F, Kaushik S, Shanbhag DD, Seo Y, et al. Direct input function estimation. Mol Imaging Biol. (2010) 12:286–94.
PseudoCT generation for pelvis PET/MRI attenuation correction using deep doi: 10.1007/s11307-009-0273-5
convolutional neural networks with multi-parametric MRI: zero echo-time 243. Lyoo CH, Zanotti-Fregonara P, Zoghbi SS, Liow J-S, Xu R, Pike VW,
and dixon deep pseudoCT (ZeDD-CT). J Nucl Med. (2018). 59:852–8. et al. Image-derived input function derived from a supervised clustering
doi: 10.2967/jnumed.117.198051 algorithm: methodology and validation in a clinical protocol using [11C](R)-
224. Arabi H, Zaidi H. Whole-body bone segmentation from MRI for PET/MRI rolipram. PLoS ONE (2014) 9:e89101. doi: 10.1371/journal.pone.0089101
attenuation correction using shape-based averaging. Med Phys. (2016) 244. Simončič U, Zanotti-Fregonara P. Image-derived input function with factor
43:5848–61. doi: 10.1118/1.4963809 analysis and a-priori information. Nucl Med Commun. (2015) 36:187–93.
225. Mehranian A, Arabi H, Zaidi H. Vision 20/20: magnetic resonance imaging- doi: 10.1097/MNM.0000000000000231
guided attenuation correction in PET/MRI: challenges, solutions, and 245. Litton JE. Input function in PET brain studies using MR-
opportunities. Med Phys. (2016) 43:1130–55. doi: 10.1118/1.4941014 defined arteries. J Comput Assist Tomogr. (1997) 21:907–9.
226. Olin A, Ladefoged CN, Langer NH, Keller SH, Löfgren JO, Hansen doi: 10.1097/00004728-199711000-00012
AE. Reproducibility of MR-based attenuation maps in PET/MRI and 246. Fung EK, Carson RE. Cerebral blood flow with [15 O]water PET studies using
the impact on PET quantification in lung cancer. J Nucl Med. (2017). an image-derived input function and MR-defined carotid centerlines. Phys
doi: 10.2967/jnumed.117.198853. [Epub ahead of print]. Med Biol. (2013) 58:1903–23. doi: 10.1088/0031-9155/58/6/1903

Frontiers in Physics | www.frontiersin.org 27 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

247. Da Silva NA, Herzog H, Weirich C, Tellmann L, Rota Kops E, Hautzel H, 267. Kjaer A. Hybrid imaging with PET/CT and PET/MR. Cancer Imaging (2014)
et al. Image-derived input function obtained in a 3TMR-brainPET. Nucl 14(Suppl. 1):O32. doi: 10.1186/1470-7330-14-S1-O32
Instruments Methods Phys Res Sect A Accel Spectrometers Detect Assoc Equip. 268. Visvikis D, Hatt M, Tixier F, Le Rest CC. The age of reason for FDG
(2013) 702:22–5. doi: 10.1016/j.nima.2012.08.030 PET image-derived indices. Eur J Nucl Med Mol Imaging (2012) 39:1670–2.
248. Jochimsen TH, Zeisig V, Schulz J, Werner P, Patt M, Patt J. Fully automated doi: 10.1007/s00259-012-2239-0
calculation of image-derived input function in simultaneous PET/MRI in a 269. Haralick RM, Shanmugam K, Dinstein I. Textural features
sheep model. EJNMMI Phys. (2016) 3:2. doi: 10.1186/s40658-016-0139-2 for image classification. Syst Man Cybern. (1973) 3:610–21.
249. Khalighi MM, Deller TW, Fan AP, Gulaka PK, Shen B, Singh P. Image- doi: 10.1109/TSMC.1973.4309314
derived input function estimation on a TOF-enabled PET/MR for cerebral 270. Tixier F, Le Rest CC, Hatt M, Albarghach N, Pradier O, Metges
blood flow mapping. J Cereb Blood Flow Metab. (2018) 38:126–35. JP. Intratumor heterogeneity characterized by textural features on
doi: 10.1177/0271678X17691784 baseline 18F-FDG PET images predicts response to concomitant
250. van Assema DM, Lubberink M, Bauer M, van der Flier WM, Schuit RC, radiochemotherapy in esophageal cancer. J Nucl Med. (2011) 52:369–78.
Windhorst AD. Blood-brain barrier P-glycoprotein function in Alzheimer’s doi: 10.2967/jnumed.110.082404
disease. Brain. (2012) 135:181–9. doi: 10.1093/brain/awr298 271. Bundschuh RA, Dinges J, Neumann L, Seyfried M, Zsótér N, Papp L, et al.
251. Gjedde A. Calculation of cerebral glucose phosphorylation from brain uptake Textural parameters of tumor heterogeneity in 18F-FDG PET/CT for therapy
of glucose analogs in vivo: a re-examination. Brain Res Rev. (1982) 4:237–74. response assessment and prognosis in patients with locally advanced rectal
doi: 10.1016/0165-0173(82)90018-2 cancer. J Nucl Med. (2014) 55:891–7. doi: 10.2967/jnumed.113.127340
252. Patlak CS, Blasberg RG, Fenstermacher JD. Graphical evaluation of blood-to- 272. Pyka T, Bundschuh RA, Andratschke N, Mayer B, Specht HM, Papp L.
brain transfer constants from multiple-time uptake data. J Cereb Blood Flow Textural features in pre-treatment [F18]-FDG-PET/CT are correlated with
Metab. (1983) 3:1–7. doi: 10.1038/jcbfm.1983.1 risk of local recurrence and disease-specific survival in early stage NSCLC
253. Logan J, Fowler JS, Volkow ND, Wolf AP, Dewey SL, Schlyer DJ, patients receiving primary stereotactic radiation therapy. Radiat Oncol.
et al. Graphical analysis of reversible radioligand binding from time- (2015) 10:100. doi: 10.1186/s13014-015-0407-7
activity measurements applied to [N-11C-methyl]-(-)-cocaine PET studies 273. Desseroit MC, Visvikis D, Tixier F, Majdoub M, Perdrisot R, Guillevin R.
in human subjects. J Cereb Blood Flow Metab. (1990) 10:740–7. Development of a nomogram combining clinical staging with 18F-FDG
doi: 10.1038/jcbfm.1990.127 PET/CT image features in non-small-cell lung cancer stage I-III. Eur J Nucl
254. Kety SS, Schmidt CF. The determination of cerebral blood flow in man by Med Mol Imaging (2016) 43:1477–85. doi: 10.1007/s00259-016-3325-5
the use of nitrous oxide in low concentrations. Am J Physiol Content. (1945) 274. Brooks FJ, Grigsby PW. The effect of small tumor volumes on studies of
143:53–66. doi: 10.1152/ajplegacy.1945.143.1.53 intratumoral heterogeneity of tracer uptake. J Nucl Med. (2014) 55:37–42.
255. Kety SS, Schmidt CF. The nitrous oxide method for the quantitative doi: 10.2967/jnumed.112.116715
determination of cerebral blood flow in man: theory, procedure and normal 275. Hatt M, Majdoub M, Vallières M, Tixier F, Le Rest CC, Groheux D. 18F-
values. J Clin Invest. (1948) 27:476–83. doi: 10.1172/JCI101994 FDG PET uptake characterization through texture analysis: investigating
256. Mintun MA, Raichle ME, Kilbourn MR, Wooten GF, Welch MJ. A the complementary nature of heterogeneity and functional tumor volume
quantitative model for the in vivo assessment of drug binding sites in a multi-cancer site patient cohort. J Nucl Med. (2015) 56:38–44.
with positron emission tomography. Ann Neurol. (1984) 15:217–27. doi: 10.2967/jnumed.114.144055
doi: 10.1002/ana.410150302 276. Zhao B, Tan Y, Tsai WY, Qi J, Xie C, Lu L. Reproducibility of radiomics
257. Bailey DL, Pichler BJ, Gückel B, Barthel H, Beer AJ, Bremerich J. Combined for deciphering tumor phenotype with imaging. Sci Rep. (2016) 6:23428.
PET/MRI: multi-modality multi-parametric imaging is here. Mol Imaging doi: 10.1038/srep23428
Biol. (2015) 17:595–608. doi: 10.1007/s11307-015-0886-9 277. Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures,
258. Padhani AR, Miles KA. Multiparametric imaging of tumor response to they are data. Radiology. (2015) 278:151169. doi: 10.1148/radiol.2015151169
therapy 1. Radiology (2010) 256:348–64. doi: 10.1148/radiol.10091760 278. Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho
259. Singh D, Miles K. Multiparametric PET/CT in oncology. Cancer Imaging. S. Decoding tumour phenotype by noninvasive imaging using a quantitative
(2012) 12:336–44. doi: 10.1102/1470-7330.2012.9007 radiomics approach. Nat Commun. (2014) 5:4006. doi: 10.1038/ncomms5006
260. Hatt M, Tixier F, Pierce L, Kinahan PE, Le Rest CC, Visvikis D. 279. Parmar C, Leijenaar RT, Grossmann P, Rios Velazquez E, Bussink J, Rietveld
Characterization of PET/CT images using texture analysis: the past, the D. Radiomic feature clusters and Prognostic Signatures specific for Lung and
present. . . any future? Eur J Nucl Med Mol Imaging (2017) 44:151–65. Head and Neck cancer. Sci Rep. (2015) 5:1–10. doi: 10.1038/srep11044
doi: 10.1007/s00259-016-3427-0 280. Gillies RJ, Beyer T. PET and MRI: is the whole greater than the sum of its
261. Heiss W-D. Hybrid PET/MR imaging in neurology: present parts? Cancer Res. (2016) 76:6163–6. doi: 10.1158/0008-5472.CAN-16-2121
applications and prospects for the future. J Nucl Med. (2016) 57:993–5. 281. Court LE, Fave X, Mackin D, Lee J, Yang J, Zhang L. Computational
doi: 10.2967/jnumed.116.175208 resources for radiomics. Transl Cancer Res. (2016) 5:340–8.
262. Puttick S, Bell C, Dowson N, Rose S, Fay M. PET, MRI, and doi: 10.21037/tcr.2016.06.17
simultaneous PET/MRI in the development of diagnostic and 282. Kourou K, Exarchos TP, Exarchos KP, Karamouzis M V., Fotiadis DI.
therapeutic strategies for glioma. Drug Discov Today (2015) 20:306–17. Machine learning applications in cancer prognosis and prediction. Comput
doi: 10.1016/j.drudis.2014.10.016 Struct Biotechnol J. (2015) 13:8–17. doi: 10.1016/j.csbj.2014.11.005
263. Nensa F, Poeppel TD, Tezgah E, Heusch P, Nassenstein K, Forsting M, et al. 283. Kickingereder P, Burth S, Wick A, Götz M, Eidel O, Schlemmer HP.
Integrated assessment of cardiac PET/MRI: co-registered PET and MRI polar Radiomic profiling of glioblastoma: identifying an imaging predictor
plots by mutual mr-based segmentation of the left ventricular myocardium. of patient survival with improved performance over established
World J Cardiovasc Dis (2017) 7:91–104. doi: 10.4236/wjcd.2017.74010 clinical and radiologic risk models. Radiology (2016) 280:880–9.
264. Rischpler C, Nekolla SG, Dregely I, Schwaiger M. Hybrid PET/MR imaging doi: 10.1148/radiol.2016160845
of the heart: potential, initial experiences, and future prospects. J Nucl Med. 284. Ypsilantis PP, Siddique M, Sohn HM, Davies A, Cook G, Goh V.
(2013) 54:402–15. doi: 10.2967/jnumed.112.105353 Predicting response to neoadjuvant chemotherapy with PET imaging
265. Spick C, Herrmann K, Czernin J. 18F-FDG PET/CT and PET/MRI perform using convolutional neural networks. PLoS ONE (2015) 10:e0137036.
equally well in cancer: evidence from studies on more than 2,300 patients. J doi: 10.1371/journal.pone.0137036
Nucl Med. (2016) 57:420–30. doi: 10.2967/jnumed.115.158808 285. Wu W, Parmar C, Grossmann P, Quackenbush J, Lambin P, Bussink J.
266. Even-Sapir E, Keidar Z, Bar-Shalom R. Hybrid imaging (SPECT/CT and Exploratory study to identify radiomics classifiers for lung cancer histology.
PET/CT)-improving the diagnostic accuracy of functional/metabolic Front Oncol. (2016) 6:71. doi: 10.3389/fonc.2016.00071
and anatomic imaging. Semin Nucl Med. (2009) 39:264–75. 286. Wu J, Cui Y, Sun X, Cao G, Li B, Ikeda DM. Unsupervised clustering of
doi: 10.1053/j.semnuclmed.2009.03.004 quantitative image phenotypes reveals breast cancer subtypes with distinct

Frontiers in Physics | www.frontiersin.org 28 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

prognoses and molecular pathways. Clin Cancer Res. (2017) 23:3334–2. 305. Roche DG, Kruuk LEB, Lanfear R, Binning SA. Public data archiving
doi: 10.1158/1078-0432.CCR-16-2415 in ecology and evolution: how well are we doing? PLoS Biol. (2015)
287. National Electrical Manufacturers Association. NEMA Standards Publication 13:e1002295. doi: 10.1371/journal.pbio.1002295
NU 4 - 2008 Performance Measurements of Small Animal Positron Emission 306. Tenopir C, Allard S, Douglass K, Aydinoglu AU, Wu L, Read E. Data
Tomographs. Rosslyn, VA (2008). sharing by scientists: practices and perceptions. PLoS ONE (2011) 6:e21101.
288. National Electrical Manufacturers Association. NEMA Standards Publication doi:10.1371/journal.pone.0021101
NU-2 2012 Performance Measurements of Positron Emission Tomographs. 307. Santos C, Blake J, States DJ. Supplementary data need to be kept in public
Rosslyn, VA (2012). repositories. Nature (2005) 438:738. doi: 10.1038/438738a
289. National Electrical Manufacturers Association. NEMA Standards Publication 308. Bray T, Paoli J, Sperberg-McQueen CM, Maler E, Yergeau F. Extensible
NU-1 2012 Performance Measurements of Gamma Cameras. Rosslyn, VA markup language (XML). World Wide Web J. (1997) 2:27–66.
(2012). 309. Kadrmas DJ, Hoffman JM. Methodology for quantitative rapid multi-
290. Boellaard R. Standards for PET image acquisition and quantitative tracer PET tumor characterizations. Theranostics (2013) 3:757–73.
data analysis. J Nucl Med. (2008) 50(Suppl. 1):11S−20S. doi: 10.7150/thno.5201
doi: 10.2967/jnumed.108.057182 310. El Fakhri G, Moore SC, Maksud P, Aurengo A, Kijewski MF. Absolute activity
291. Beyer T, Czernin J, Freudenberg LS. Variations in clinical PET/CT quantitation in simultaneous 123I/99mTc brain SPECT. J Nucl Med. (2001)
operations: results of an international survey of active PET/CT users. J Nucl 42:300–8.
Med. (2011) 52:303–10. doi: 10.2967/jnumed.110.079624 311. El Fakhri G, Sitek A, Zimmerman RE, Ouyang J. Generalized five-
292. Wieder H, Freudenberg LS, Czernin J, Navar BN, Israel O, Beyer dimensional dynamic and spectral factor analysis. Med Phys. (2006)
T. Variations of clinical SPECT/CT operations. Nuklearmedizin (2012) 33:1016–24. doi: 10.1118/1.2179168
51:154–60. doi: 10.3413/Nukmed-0467-12-01 312. Verhaeghe J, D’Asseler Y, Staelens S, Lemahieu I. Noise properties of
293. Chauvie S, Bergesio F. The strategies to homogenize PET/CT metrics: simultaneous dual tracer PET imaging. In: IEEE Nuclear Science Symposium
the case of onco-haematological clinical trials. Biomedicines. (2016) 4:26. Conference Record, Vol 5. Fajardo: IEEE (2005). p. 2611–4.
doi: 10.3390/biomedicines4040026 313. Figueiras FP, Jiménez X, Pareto D, Gómez V, Llop J, Herance R.
294. Boellaard R, Delgado-Bolton R, Oyen WJ, Giammarile F, Tatsch K, Simultaneous dual-tracer PET imaging of the rat brain and its application
Eschner W. FDG PET/CT: EANM procedure guidelines for tumour in the study of cerebral ischemia. Mol Imaging Biol. (2011) 13:500–10.
imaging: version 2.0. Eur J Nucl Med Mol Imaging (2014) 42:328–54. doi: 10.1007/s11307-010-0370-5
doi: 10.1007/s00259-014-2961-x 314. Koeppe RA, Ficaro EP, Raffel DM, Minoshima S, Kilbourn MR. Temporally
295. Varrone A, Asenbaum S, Vander Borght T, Booij J, Nobili F, Någren overlapping dual-tracer PET studies. In: Carson RE, Daube-Witherspoon
K, et al. EANM procedure guidelines for PET brain imaging using ME, Herscovitch P, editors. Quantitative Functional Brain Imaging with
[18F]FDG, version 2. Eur J Nucl Med Mol Imaging. (2009) 36:2103–10. Positron Emission Tomography. San Diego, CA: Elsevier (1998). p. 359–66.
doi: 10.1007/s00259-009-1264-0 315. Kadrmas DJ, Rust TC. Feasibility of rapid multitracer PET tumor imaging.
296. Flotats A, Knuuti J, Gutberlet M, Marcassa C, Bengel FM, Kaufmann IEEE Trans Nucl Sci. (2005) 52:1341–7. doi: 10.1109/TNS.2005.858230
PA. Hybrid cardiac imaging: SPECT/CT and PET/CT. A joint position 316. Andreyev A, Celler A. Dual-isotope PET using positron-gamma emitters.
statement by the European Association of Nuclear Medicine (EANM), the Phys Med Biol. (2011) 56:4539–56. doi: 10.1088/0031-9155/56/14/020
European Society of Cardiac Radiology (ESCR) and the European Council of 317. Lage E, Herraiz JL, Parot V. Multiplexable Emission Tomography.
Nuclear Cardiology (ECNC). Eur J Nucl Med Mol Imaging. (2011) 38:201–12. (2013). Available online at: https://patents.google.com/patent/
doi: 10.1007/s00259-010-1586-y WO2013188011A1/en.
297. Chalkidou A, O’Doherty MJ, Marsden PK. False discovery rates in PET and 318. Cal-González J, Lage E, Herranz E, Vicente E, Udias JM, Moore SC.
CT studies with texture features: a systematic review. PLoS ONE (2015) Simulation of triple coincidences in PET. Phys Med Biol. (2015) 60:117–36.
10:e0124165. doi: 10.1371/journal.pone.0124165 doi: 10.1088/0031-9155/60/1/117
298. Molina D, Pérez-Beteta J, Martínez-González A, Martino J, Velasquez 319. Strobel K, Heinrich S, Bhure U, Soyka J, Veit-Haibach P, Pestalozzi BC, et al.
C, Arana E. Lack of robustness of textural measures obtained from 3D Contrast-enhanced 18 F-FDG PET/CT: 1-stop- shop imaging for assessing
brain tumor MRIs impose a need for standardization. PLoS ONE (2017) the resectability of pancreatic cancer. J Nucl Med. (2008) 49:1408–13.
12:e0178843. doi: 10.1371/journal.pone.0178843 doi: 10.2967/jnumed.108.051466
299. Nyflot MJ, Yang F, Byrd D, Bowen SR, Sandison GA, Kinahan PE. 320. Badiee S, Franc BL, Webb EM, Chu B, Hawkins RA, Coakley F. Role of IV
Quantitative radiomics: impact of stochastic effects on textural feature iodinated contrast material in 18 F-FDG PET/CT of liver metastases. Am J
analysis implies the need for standards. J Med Imaging. (2015) 2:41002. Roentgenol. (2008) 191:1436–9. doi: 10.2214/AJR.07.3750
doi: 10.1117/1.JMI.2.4.041002 321. Asagi A, Ohta K, Nasu J, Tanada M, Nadano S, Nishimura R. Utility of
300. van Velden FH, Kramer GM, Frings V, Nissen IA, Mulder ER, contrast-enhanced FDG-PET/CT in the clinical management of pancreatic
de Langen AJ. Repeatability of radiomic features in non-small-cell cancer. Pancreas (2013) 42:11–19. doi: 10.1097/MPA.0b013e3182550d77
lung cancer [18F]FDG-PET/CT studies: impact of reconstruction and 322. Kuhn FP, Hullner M, Mader CE, Kastrinidis N, Huber GF, von Schulthess
delineation. Mol Imaging Biol. (2016) 18:788–95. doi: 10.1007/s11307-016- GK, et al. Contrast-enhanced PET/MR imaging versus contrast-enhanced
0940-2 PET/CT in head and neck cancer: how much mr information is needed? J
301. Tixier F, Hatt M, Le Rest CC, Le Pogam A, Corcos L, Visvikis D. Nucl Med. (2014) 55:551–8. doi: 10.2967/jnumed.113.125443
Reproducibility of tumor uptake heterogeneity characterization through 323. Thakor AS, Gambhir SS. Nanooncology: the future of cancer diagnosis and
textural feature analysis in 18F-FDG PET. J Nucl Med. (2012) 53:693–700. therapy. CA Cancer J Clin. (2013) 63:395–418. doi: 10.3322/caac.21199
doi: 10.2967/jnumed.111.099127 324. Yen SK, Padmanabhan P, Selvan ST. Multifunctional iron oxide
302. Leijenaar RT, Carvalho S, Velazquez ER, van Elmpt WJ, Parmar C, Hoekstra nanoparticles for diagnostics, therapy and macromolecule delivery.
OS. Stability of FDG-PET radiomics features: AN integrated analysis of test- Theranostics (2013) 3:986–1003. doi: 10.7150/thno.4827
retest and inter-observer variability. Acta Oncol (Madr) (2013) 52:1391–7. 325. Fakhri G El. Ready for prime time? Dual tracer PET and SPECT imaging. Am
doi: 10.3109/0284186X.2013.812798 J Nucl Med Mol Imaging (2012) 2:415–7.
303. Costello MJ. Motivating online publication of data. Bioscience (2009) 326. Chapman SE, Diener JM, Sasser TA, Correcher C, González AJ, Avermaete
59:418–27. doi: 10.1525/bio.2009.59.5.9 TV. Dual tracer imaging of SPECT and PET probes in living mice using a
304. Wicherts JM, Bakker M, Molenaar D. Willingness to share research sequential protocol. Am J Nucl Med Mol Imaging. (2012) 2:405–14.
data is related to the strength of the evidence and the quality 327. Shi X, Meng X, Sun X, Xing L, Yu J. PET/CT imaging-guided
of reporting of statistical results. PLoS ONE (2011) 6:e26828. dose painting in radiation therapy. Cancer Lett. (2014) 355:169–75.
doi: 10.1371/journal.pone.0026828 doi: 10.1016/j.canlet.2014.07.042

Frontiers in Physics | www.frontiersin.org 29 May 2018 | Volume 6 | Article 47


Cal-Gonzalez et al. Hybrid Imaging: Instrumentation and Data Processing

328. Thomas HM, Kinahan PE, Samuel JJE, Bowen SR. Impact of tumour motion 336. Lapi SE, Lewis JS, Dehdashti F. Evaluation of hypoxia with copper-
compensation and delineation methods on FDG PET-based dose painting labeled diacetyl-bis(N-methylthiosemicarbazone). Semin Nucl Med. (2015)
plan quality for NSCLC radiation therapy. J Med Imaging Radiat Oncol. 45:177–85. doi: 10.1053/j.semnuclmed.2014.10.003
(2018) 62:81–90. doi: 10.1111/1754-9485.12693 337. Knopf A-C, Lomax A. In vivo proton range verification: a review. Phys Med
329. Daniel M, Andrzejewski P, Sturdza A, Majercakova K, Baltzer P, Pinker Biol. (2013) 58:R131–60. doi: 10.1088/0031-9155/58/15/R131
K. Impact of hybrid PET/MR technology on multiparametric imaging and 338. Sportelli G, Belcari N, Camarlinghi N, Cirrone GA, Cuttone G,
treatment response assessment of cervix cancer. Radiother Oncol. (2017) Ferretti S. First full-beam PET acquisitions in proton therapy with a
125:420–5. doi: 10.1016/j.radonc.2017.10.036 modular dual-head dedicated system. Phys Med Biol. (2014) 59:43–60.
330. Bradley J, Thorstad WL, Mutic S, Miller TR, Dehdashti F, Siegel BA. Impact doi: 10.1088/0031-9155/59/1/43
of FDG-PET on radiation therapy volume delineation in non–small-cell 339. Shao Y, Sun X, Lou K, Zhu XR, Mirkovic D, Poenisch F. In-beam PET
lung cancer. Int J Radiat Oncol. (2004) 59:78–86. doi: 10.1016/j.ijrobp.2003. imaging for on-line adaptive proton therapy: an initial phantom study. Phys
10.044 Med Biol. (2014) 59:3373–88. doi: 10.1088/0031-9155/59/13/3373
331. Souvatzoglou M, Krause BJ, Pürschel A, Thamm R, Schuster T, Buck AK, 340. Ferrero V, Fiorina E, Morrocchi M, Pennazio F, Baroni G, Battistoni G.
et al. Influence of 11C-choline PET/CT on the treatment planning for Online proton therapy monitoring: clinical test of a Silicon-photodetector-
salvage radiation therapy in patients with biochemical recurrence of prostate based in-beam PET. Sci Rep. (2018) 8:4100. doi: 10.1038/s41598-018-
cancer. Radiother Oncol. (2011) 99:193–200. doi: 10.1016/j.radonc.2011. 22325-6
05.005 341. Parodi K, Bortfeld T, Haberer T. Comparison between in-beam and
332. Pasqualetti F, Panichi M, Sainato A, Matteucci F, Galli L, Cocuzza P. offline positron emission tomography imaging of proton and carbon ion
[(18)F]Choline PET/CT and stereotactic body radiotherapy on treatment therapeutic irradiation at synchrotron- and cyclotron-based facilities.
decision making of oligometastatic prostate cancer patients: preliminary Int J Radiat Oncol. (2008) 71:945–56. doi: 10.1016/j.ijrobp.2008.
results. Radiat Oncol. (2016) 11:9. doi: 10.1186/s13014-016-0586-x 02.033
333. Alongi F, Fersino S, Giaj Levra N, Mazzola R, Ricchetti F, Fiorentino
A. Impact of 18F-choline PET/CT in the decision-making strategy Conflict of Interest Statement: The authors declare that the research was
of treatment volumes in definitive prostate cancer volumetric conducted in the absence of any commercial or financial relationships that could
modulated radiation therapy. Clin Nucl Med. (2015) 40:e496–500. be construed as a potential conflict of interest.
doi: 10.1097/RLU.0000000000000841
334. McKinley ET, Ayers GD, Smith RA, Saleh SA, Zhao P, Washington MK, et al. Copyright © 2018 Cal-Gonzalez, Rausch, Shiyam Sundar, Lassen, Muzik, Moser,
Limits of [18F]-FLT PET as a biomarker of proliferation in oncology. PLoS Papp and Beyer. This is an open-access article distributed under the terms of
ONE (2013) 8:e58938. doi: 10.1371/journal.pone.0058938 the Creative Commons Attribution License (CC BY). The use, distribution or
335. Segard T, Robins PD, Yusoff IF, Ee H, Morandeau L, Campbell EM, et al. reproduction in other forums is permitted, provided the original author(s) and
Detection of hypoxia with 18F-fluoromisonidazole (18F-FMISO) PET/CT the copyright owner are credited and that the original publication in this journal
in suspected or proven pancreatic cancer. Clin Nucl Med. (2013) 38:1–6. is cited, in accordance with accepted academic practice. No use, distribution or
doi: 10.1097/RLU.0b013e3182708777 reproduction is permitted which does not comply with these terms.

Frontiers in Physics | www.frontiersin.org 30 May 2018 | Volume 6 | Article 47

You might also like