Assessment of a digital and an analog PET/CT system for accurate myocardial perfusion imaging with a flow phantom.

In Myocardial Perfusion Imaging (MPI) with Positron Emission Tomography/Computed Tomography (PET/CT) systems, accurate quantification is essential. We assessed flow quantification accuracy over various injected activities using a flow phantom.
METHODS: The study was performed on the digital 4-ring Discovery MI (DMI-20) and analog Discovery 690 (D690) PET/CT systems, using 325-1257 MBq of [15O]H2O. PET performance and flow quantification accuracy were assessed in terms of count-rates, dead-time factors (DTF), scatter fractions (SF), time-activity curves (TACs), areas-under-the-curves (AUCs) and flow values.
RESULTS: On DMI-20, prompts of 12.8 Mcps, DTF of 2.06 and SF of 46.1% were measured with 1257 MBq of activity. On the D690, prompts of 6.85 Mcps, DTF of 1.57 and SF of 32.5% were measured with 1230 MBq of activity. AUC values were linear over all activities. Mean wash-in flow error was - 9% for both systems whereas wash-out flow error was - 5% and - 6% for DMI-20 and D690. With the highest activity, wash-out flow error was - 12% and - 7% for the DMI-20 and D690.
CONCLUSION: DMI-20 and D690 preserved accurate flow quantification over all injected activities, with maximum error of - 12%. In the future, flow quantification accuracy over the activities and count-rates evaluated in this study should be assessed.

Introduction

Myocardial Perfusion Imaging (MPI) with dynamic Positron Emission Tomography (PET) allows quantitative measurements of absolute Myocardial Blood Flow (MBF).1 MPI assessed using PET and short-lived tracers is highly reproducible.2-5 However, physical performance of PET systems has been a key issue in achieving accurate flow quantification over a wide range of activities.6-9 PET systems are exposed to high count-rates and dead-time when activity is present in the injected vein and heart chamber, which may compromise image accuracy and quality.10 Due to dead-time effects, time-activity curves (TACs) from left ventricle blood pool and myocardial tissue may be underestimated, resulting in MBF quantification inaccuracy.2,11-13

With the introduction of digital PET/Computed Tomography (PET/CT) systems, improvements in count-rate performance and image quality have been shown with National Electrical Manufacturers Association (NEMA) tests.14,15 Therefore, there might be potential advantages in applying digital PET/CT systems in MPI quantification. Thus, the count-rate capabilities of digital systems need to be assessed with phantom studies at high count-rates representative to MPI measurements. Evaluation of flow quantification at various injected activities is of paramount importance to understand the effects of count-rate performance, which differs among PET systems.6,7,9,16

Most importantly, flow quantification accuracy of digital PET/CT systems over various count-rates should be evaluated not only on static phantoms but with phantoms simulating dynamic imaging conditions. In 16 the authors investigated perfusion quantification accuracy at various activities by using a PET/MR compatible flow phantom. A similar phantom modelling realistic cardiac flow was presented and validated for MPI in PET in Reference 17, allowing flow modelling and quantification against a ground truth value. Thus, PET system assessment using a flow phantom allows to evaluate both count-rate performance and flow quantification accuracy.

In this study, we aimed at assessing flow quantification accuracy of a novel digital PET/CT system in comparison to an analog PET/CT system. A state-of-the-art flow phantom was selected to evaluate the accuracy of flow quantification for various activities up to 1260 MBq of [15O]H2O at scan start time.

Materials and Methods

PET/CT Systems

The digital Discovery MI with four detector rings (DMI-20) 15 and analog Discovery 690 (D690) (GE Healthcare, US) 18 PET/CT systems were assessed. Performance characteristics are summarized in Supplemental Table 1.

Table 1

Activities at scan start time with flow meter readings. Qpump was 200 mL/min and Qcyl was 60% of Qpump

Measurement	Activity at scan start time		Qcyl		Qtube
	MBq		mL/min		ML/min
	DMI-20	D690	DMI-20	D690	DMI-20	D690
1	325	359	123	120	79.0	91.0
2	488	400	121	117	82.0	91.0
3	546	532	121	117	83.0	93.0
4	621	607	120	121	83.0	92.0
5	655	729	122	120	82.0	92.0
6	691	833	119	117	84.0	93.0
7	906	995	118	119	85.0	93.0
8	1060	1130	117	118	86.0	93.0
9	1257	1230	117	115	87.0	94.0
Mean ± SD	-	-	120 ± 2.2	118 ± 2.0	83.4 ± 2.4	92.4 ± 1.0

DMI-20 has 4 rings of silicon photomultiplier (SiPM) detectors. One detector block comprises 4 × 9 array of 3.95 × 5.3 × 25 mm LBS-crystals coupled to 3 × 6 array of SiPM with 2 × 3 pixels. Transaxial and axial field-of-views (FOV) are 70 cm and 20 cm.15

D690 has 24 detector rings with photomultiplier tube (PMT) blocks. One block comprises 9 × 6 array of 4.2 × 6.3 × 25 mm LBS-crystals coupled to a single square PMT with 4 anodes. Transaxial and axial FOVs are 70 cm and 15.7 cm.18

Both systems employ similar randoms, scatter and dead-time corrections in image reconstruction: randoms from singles,19 model-based scatter correction with extension to time-of-flight (TOF) 20,21 and dead-time measured from the detector pile-up losses.22,23

Flow Phantom System

To investigate flow quantification accuracy, a flow phantom (DCE Dynamic Flow Phantom, Shelley Medical Imaging Technologies, Canada) was used. Phantom validation for MPI in PET is presented in Reference 17.

The phantom set-up contains water container, peristaltic pump, injection port, phantom shell, flow constrictor valves and flow meters. The phantom schematic is shown in Supplemental Figure. 1. The shell has similar dimensions as the NEMA image quality phantom (diameter of 31 cm). An input chamber (volume of 15.7 mL) and an exchange cylinder (volume of 161 mL) containing a perforated tube (volume of 35 mL) are located inside the shell. The input chamber models the left ventricle blood pool whereas the exchange cylinder and perforated tube model tracer exchange in myocardial tissue. These allow to measure input and tissue TACs for kinetic modelling.

Flow inside the system is controlled with the peristaltic pump . Flow inside the exchange cylinder () and perforated tube () are controlled with flow constrictor valves and recorded with microturbine flow meters (Omega Engineering Inc., US).24 Ideally, and should be equal to . Reference flow () values are related to and derived from flow meter calibration using a lookup table. Further information from the flow values is given in Supplemental Figure. 1.

Kinetic Modeling of Flow Values

The input and tissue TACs can be measured by a volume-of-interest (VOI) analysis, defined for the input chamber and exchange cylinder, namely and . VOIs were selected as in Reference 5. The TACs are used for image-based flow () modelling based on a two-compartment kinetic (one-tissue compartment) model as implemented in the phantom flow quantification software (QuantifyDCE 1.1, Shelley Medical Imaging Technologies, Canada).17

The software models tracer concentration in the exchange cylinder as:where represents discrete convolution operation, the time-dependent tracer concentration in the exchange cylinder, the tracer wash-in rate to the exchange cylinder (min−1), and the tracer wash-out rate from the exchange cylinder (min−1). represents tracer concentration in the perforated tube, which can be estimated from the input chamber TAC by transport delay of when Moreover, as contains contributions from both and , a signal mixing correction function is implemented as:

Thus, complete kinetic model consists of following free parameters: , and . Model fitting is performed with a weighted (by frame length) least squares algorithm.

The modeled parameters and can be considered analogous to rate parameters and in kinetic modelling of blood-tissue exchange. Final image-derived flow values and can be computed based on and , multiplied with the cylinder volume (), representing flow from the perforated tube to the exchange cylinder () and flow out from the exchange cylinder (). Ideally, these flow values should be equal to the reference flow value, i.e. (mL/min). 17

PET and CT Data Acquisition

At the start of the measurement, the flow meters, peristaltic pump (MasterFlex L/S, 07522-20, US) and [15O]H2O dispenser (Hidex Oy, Finland) were calibrated. Flow meter calibration and flow rates were fixed for all studies. was set to 200 mL/min, while was adjusted to 60% of , following previous protocols.5,17

A CT-based attenuation correction (CTAC) using a tube voltage of 120 kV and current of 67mAs was acquired. The flow phantom was imaged in the center of the PET FOV. Dynamic list-mode PET acquisition was started 50 s after the dose injection, with a duration of 4 min and 40 s, following the clinical protocol.25 Measurements were repeated 9 times with individual doses. Injected doses were measured automatically by the [15O]H2O dispenser (Hidex Oy, Finland), which is cross-calibrated to a Veenstra VDC-404 (Veenstra Instruments, Netherlands) dose calibrator. Activities at scan start times and flow rates are reported in Table 1.

Data Reconstruction

List-mode data was binned into 24 time-frames of 14 × 5 s, 3 × 10 s, 3 × 20 s and 4 × 30 s. On both systems, reconstructions were performed using three-dimensional ordered-subset-expectation-maximization (3D-OSEM) algorithm using point-spread function modelling (PSF) and TOF (vendor name: VUE Point FX-S), with 5 mm Gaussian post-filter and 35 cm FOV. On DMI-20, 3 iterations, 16 subsets and matrix size of 192 × 192 were used. On D690, 2 iterations, 24 subsets and matrix size of 128 × 128 were used. All data corrections including decay, attenuation, scatter, randoms and dead-time were performed.

Flow Quantification Accuracy

Count-rate performance was assessed for each measurement by extracting count-rate data from DICOM headers. Specifically, peak prompts (P), randoms (R), dead-time factors (DTF) and scatter fractions (SF) corresponding the peak prompts were extracted.

Scatter rates (S) and trues rates (T) were derived based on the extracted P, R and SF. The scatter rate was derived as:and the trues rate was derived as:

We report the peak count-rates (P, T, R, S), DTF and SF as function of activity at scan start time.

Flow quantification accuracy was assessed in terms of shape of the input and tissue TACs and area-under-the-curve (AUC) of the TACs. Linearity of AUCs was determined with regression analysis of the AUC values versus activities at scan start time. We report AUCs and goodness-of-fits (R2) of the regression slopes.

Finally, accuracy of the modeled flow values and with respect to the reference flow was calculated as relative error:where represents either or for each measurement. We report the accuracy of flow values with activity at scan start time. All flow values are reported as mean ± standard deviation (SD).

Results

Figure 1 shows peak prompts, scatter, randoms and trues as a function of activity at scan start time on DMI-20 and D690. On both systems, count-rates increased in relation to activities, as expected. For DMI-20, peak prompts were 12.8 Mcps with 1257 MBq. On D690, peak prompts were 6.85 Mcps with 1230 MBq. Additional data of the DMI prompts and randoms for the measurement with the highest injected activity (1257 MBq) is given in Supplementary Data 1.

Figure 1

Count-rates for (a) prompts, scatter, randoms and trues and (b) zoomed-in plot for scatter and trues for DMI-20 and D690 systems for all flow phantom measurements as a function of activity at scan start time (Table 1)

DTFs as a function of activity at scan start time are presented in Figure 2A. DTFs increased nearly linearly in relation to activities, reflecting the count-rates (Figure 1). DTFs were higher with DMI-20 compared to D690. The highest DTFs 2.06 and 1.57 were measured with 1257 MBq and 1230 MBq on DMI-20 and D690.

Figure 2

(a) Dead-time factors (DTFs) and (b) scatter fraction factors (SFs) at peak prompts as a function of activity at scan start time for all flow phantom measurements

Figure 2B shows SFs as a function of activity at scan start time. Due to larger axial FOV, SFs were higher on DMI-20 compared to D690. SFs increased up to the three highest activities (906 MBq, 1060 MBq, 1257 MBq at scan start time), after which SFs were nearly constant on DMI-20. On D690, SFs increased in relation to all activities.

Figure 3 shows TACs derived from the input chamber and exchange cylinder of the flow phantom. No visible distortions in TAC shapes were seen. Two-compartmental model fitting R2 values were 0.998 and 0.999 on both DMI-20 and D690.

Figure 3

Time-activity curves (TACs) measured from the input chamber (input curve) and exchange cylinder (tissue curve) on (a) Discovery MI (DMI-20) and (b) Discovery 690 (D690) PET/CT system for all flow phantom measurements

AUCs as a function of activity at scan start time with regression analysis are presented in Figure 4. R2 values for the input and tissue AUCs were 0.995 and 0.997 for DMI-20 and 0.998 and 0.998 for D690. AUCs increased linearly in relation to all activities on both systems. Tissue AUCs followed the line-of-identity closely whereas the input AUCs were above the identity line on both systems.

Figure 4

Areas-under-the-curves (AUCs) of the time-activity curves (TACs) on the (a) Discovery MI (DMI-20) and (b) Discovery 690 (D690) PET/CT system as a function of activity at scan start time with regression slope lines and line-of-identities for all flow phantom measurements

Prompts, DTFs and SFs with the accuracy of the modeled flow values are presented in Supplemental Table 2. Errors of the modeled flow values with respect to are presented as a function of activity at scan start time in Figure 5. A slight increase in error was seen above activities of 906 MBq and 995 MBq at scan start time on DMI-20 and D690. In general, the error in flow quantification accuracy was small over all count-rates and activities. Also, the error magnitude of and can be considered small (below 5%) between DMI-20 and D690.

Figure 5

Errors of the modeled and flow values with respect to the reference flow measured on the (a) Discovery MI (DMI-20) and (b) Discovery 690 (D690) as a function of activity at scan start time for all flow phantom measurements

Discussion

We assessed the accuracy of flow quantification on digital (DMI-20) and analog (D690) PET/CT systems over various injected activities with a flow phantom, for the first time.

DMI-20 showed improved performance with higher prompts compared to D690 over similar activities (Figure 1), in line with 9 and NEMA performance studies.15 DMI-20 limits the amount of data passed to the acquisition sorter after reaching a throttle limit, due to its sensitivity. The throttle is specific for different ring configurations on DMI systems, limiting the amount of data transmitted for histogramming. The number of prompts in the DICOM header will be also capped at the scanner-specific limit due to the throttle. The DMI-20 specific 12.8 Mcps limit was reached at activity of 1257 MBq. Therefore, a sudden notch in scatter and trues rate in Figure 1A and B is seen, similarly as in Reference 16.

However, during reconstruction, a scaling factor based on the prompts seen by the coincidence processor and prompts transmitted to histogramming will be applied to preserve quantification despite the throttle. Moreover, the actual RFS estimate will be slightly lower than the value for randoms in the DICOM header above throttle, as no count losses are taken into account for the randoms value saved in the DICOM header. The interested reader is referred to Supplementary Data 1 for more information.

Higher DTFs and SFs (Figure 2) were seen on DMI-20 over all activities. SFs increased in relation to activities, stabilizing at the highest activities (Figure 2B). This might be caused by biased scatter scaling at high activities by counts outside the body contour or by increased detector scatter, which is usually seen in NEMA testing.15 However, we detected no CTAC-PET misalignment that could affect scatter scaling. In general, no artifacts were seen resulting from scatter overcorrection on either system.

For flow quantification, TACs showed no visible artifacts (Figure 3). Similarly, the input and tissue AUCs showed high linearity (R2 close to 0.99) and followed the line-of-identity well (Figure 4). Saturation of TACs and non-linearity of AUCs over activities of 594 MBq of [18F]F− 16 and saturation of input TAC with activities of 1006-1491 MBq of [82Rb] 13 have been shown previously. Saturation or non-linearity of AUCs were not seen on either system, contrary to 16 and.13 However, the input AUCs were slightly higher compared to tissue AUCs, which can contribute to the difference seen between and values (Figure 5), as we have previously discussed in.5

In line with the AUC linearity, errors in the modeled flow values were similar over all activities and systems (Figure 5). Mean error was -9% and mean error was -5% and -6% on DMI-20 and D690. The small difference between systems is contributed by different reconstruction parameters, which were selected based on clinical protocols on each system. error increased when activity increased, up to − 12% on DMI-20 at the highest activity, and up to 7% on D690. Our finding is similar to 13 measured with [82Rb], although the maximum error in (− 12%) was smaller compared to the reported error of 22%. Similarly to Reference 5, there was a systematic offset between in vs (Supplementary Table 1) and offset of 10 mL/min in Qcyl and Qtube (Table 1), which results in systematic negative bias irrespective of count-rate performance or the PET systems.

As a limitation, flow quantification was assessed only up to 1257 MBq at scan start time due to reaching the maximum limit of our [15O]H2O cyclotron production rate. Most importantly, image quantification was preserved up to highest activities. Therefore, we recommend limiting the injected activities to within the investigated range pending further investigation into values outside of this range on DMI-20 and D690. In clinical MPI studies, the injected doses for e.g. [82Rb] might exceed the limit of 1260 MBq, as both European and North American guidelines recommend injected doses from 1100 to 1500 MBq of [82Rb].26-28 At our institute, the injected activity is approximately 500 MBq, which decays to approximately 400 MBq at scan start time on both systems in clinical [15O]H2O MPI studies.

As an extension, other DMI systems should be evaluated. However, we expect that different DMI systems should behave similarly in terms of count-rate linearity despite the system-specific throttling limit. Additionally, we encourage future assessments with other perfusion tracers, such as [82Rb] or [13N] and an alternative MPI software estimation, as we only assessed [15O]H2O and the flow quantification software designed specifically for the phantom (QuantifyDCE).

Finally, our study indicates accurate MPI with [15O]H2O can be achieved if PET is acquired under the activities and count-rates specified above, as in References 6, 7, 9, and 16. The advantage of using a flow phantom is two-fold: the reference flow is always known, and CT-PET data can be collected free of any motion artefacts or misalignment. Thus, the phantom protocol used in this study could be applied for MPI harmonization studies for several PET/CT systems, according to their count-rate performance and reconstruction methods.

New Knowledge Gained

Flow quantification accuracy in MPI with PET was preserved over various injected activities on DMI-20 and D690 PET/CT systems. Mean flow quantification error was -9% for wash-in flow and − 5% and − 6% for wash-out flow on DMI-20 and D690, indicating similar flow quantification errors on both systems. As flow quantification error was − 12% at 12.8 Mcps prompts on the DMI and − 7% at 6.85 Mcps prompts on the D690, we assume DMI-20 system might offer potential benefits in MPI in PET.

Conclusions

Flow quantification accuracy was preserved over activities of 325-1257 MBq and 359-1230 MBq at scan start time on the digital 4-ring Discovery MI PET/CT system and analog Discovery 690 PET/CT system using [15O]H2O. The measurement with the highest prompts rate of 12.8 Mcps with corresponding DTF of 2.06 showed flow quantification error of − 12% on DMI-20. On D690 the highest prompts rate of 6.85 Mcps with corresponding DTF of 1.57 showed an error of − 7%. Generally, Discovery MI 4-ring system might offer potential benefits in MPI PET, although future studies should be conducted to assess the system performance at higher injected activities and prompt rates. We recommend performing MPI under 1257 MBq of activity at scan start time and a prompts rate of 12.8 Mcps on DMI-20.

Below is the link to the electronic supplementary material.

Supplementary material 1 (DOCX 264 kb)

Supplementary material 2 (PPTX 7419 kb)