# Parametric Images in Assessing Bone Grafts Using  Dynamic18F-Fluoride PET

Authors: Lingfeng Wen; Stefan Eberl; (David) Dagan Feng; Paul Stalley; Gang Huang; Michael J. Fulham
Journal: International Journal of Molecular Imaging
 (2011)

---

## Abstract

The early identification of graft failure would improve patient management.18F-fluoride is a suitable tracer for quantifying bone metabolism. Performance of parametric images constructed by Patlak graphical analysis (PGA) with various time periods was evaluated in the analysis of dynamic 18F-fluoride PET studies of eight patients with fibula bone grafts after limb salvage surgery. The PGA parametric image approach tended to underestimate influx rate. The linear portion of PGA analysis was found to be from 10 to 50 min. It shows promise in providing a quantitative assessment of the viability of bone grafts.


---

## Body

## 1. Introduction

Free vascularised fibula grafting is well described in limb salvage surgery after resection of bone tumours. Nonunion at the osteotomy site and fracture of the graft are recognised complications, which may be related to poor blood flow within the graft. Currently there is no diagnostic modality that can reliably assess graft viability. Graft viability is usually judged by evidence of bony bridging and graft hypertrophy on plain radiographs, but this is problematic and delayed union and infection cannot be reliably predicted. The early identification of graft failure would improve patient management, and a reliable assessment of graft viability would also help identify characteristics associated with successful grafts.Functional imaging techniques, such as positron emission tomography (PET), can visualise subtle metabolic changes and have the potential for assessing graft viability. The PET tracer18F-fluoride has been used to evaluate regional bone metabolism and skeletal kinetics    [1–5] as well as monitoring therapeutic response [6]. The dynamic behaviour of complex biological systems can be described by suitable kinetic models in functional imaging. Rate constants and macroparameters of kinetic models are related to biological processes; hence, these parameters could be used to aid clinical evaluation.Kinetic modelling of PET imaging usually involves obtaining a series of blood samples for the input function (IF), and a tissue time activity curve (TTAC), derived from dynamic imaging data as the output function. Rate constants are determined by parameter estimation methods, which fit the TTAC according to an underlying kinetic model. Traditionally, manually defined regions of interest (ROI) are used to derive average, but reduced noise, TTACs for the selected regions by deriving the mean activities for the defined ROI. Alternatively, parameter estimation can be conducted voxel by voxel to form a three-dimensional parametric image volume whose voxel values represent quantitative functional parameters. The parametric image approach reduces operator-dependant bias in the manual delineation of ROI and eliminates the need to know the spatial distribution of newly developed tracers [7]. Furthermore, the parametric image approach readily visualizes the spatial distribution of quantitative functional parameters.The nonlinear least square (NLS) analysis, a standard parameter estimation method, provides estimates with optimum statistical accuracy in kinetic analysis by iteratively adjusting the kinetic parameters of the nonlinear model equations [8]. There is intrinsic noise in functional imaging data, thus an appropriate weight function is usually chosen to derive the objective function in NLS fitting. Thus, NLS is often referred to as the weighted nonlinear least square (WNLS). However, NLS/WNLS are not suitable for the construction of parametric images due to their high computational cost, proneness to be trapped in local minima, and a requirement for appropriate initial parameters. Graphical analysis (GA) methods, such as the Patlak (PGA) [9] and Logan approaches    [10], transform kinetic model equations into a simple linear plot of selected variables. The slope and intercept of GA methods are usually related to parameters of interest. Computational efficiency and reliable parameter estimates make GA methods suitable for deriving parametric images.A number of investigators reported on the kinetic analysis of18F-fluoride in the evaluation of vascularisation of allogenic bone grafts [11–16]. Different kinetic quantitative approaches are also reported [5, 12, 15]. However, these comparisons were only limited to the TTACs derived by manually placing ROI on the images, which is time consuming and prone to subjective bias. Thus, the performance of these techniques for generating three-dimensional parametric images has not been fully evaluated. In addition, there is no consensus on the optimum time period in defining the linear portion of the PGA plot, which may also introduce bias in PGA plots.Our aim was to systematically investigate the performance of parametric images derived by the PGA method and evaluate the optimum linear portion of PGA using dynamic18F-fluoride PET imaging. Eight patients with limb salvage surgery and bone grafting were included with volume of interest (VOI) defined on bone grafts and normal bone tissues. Quantitative rate constants and net influx rates were also derived for VOI-derived TTACs with NLS and PGA analysis for comparison.


## 2. Materials and Methods

### 2.1. Kinetic Model for Fluoride Bone Metabolism

A three-compartment and four-parameter model has been used in the kinetic analysis for fluoride bone metabolism (Figure1) [1]. It consists of plasma, extracellular fluid, and bone mineral compartments. Rate constants describe the transport of 18F-fluoride between the compartments with K1 representing the unidirectional clearance of fluoride from plasma to bone tissue, k2 the reverse transport of fluoride from bone tissue to plasma, k3 reflecting the incorporation, and k4 the release of fluoride from the bone mineral compartment. The unit of K1 is mL·min−1·mL−1, while k2,k3, and k4 have units of min−1.Figure 1
Three-compartment and four-parameter kinetic model for fluoride bone metabolism.Cp(t) represents the plasma concentration of 18F-fluoride, while Ce(t) and Cm(t) are 18F-fluoride concentrations in extracellular fluid and bone. Ct(t) is the sum of Ce(t) and Cm(t) and represents total concentration of 18F-fluoride in bone tissue. Through kinetic modelling technique, the output function, Ct(t), can be represented by the input function, Cp(t), in the following:(1)Ct(t)=K1α2-α1[(k4-α1)⋅e-α1t+(α2-k4)⋅e-α2t]⊗Cp(t)+K1k3α2-α1(e-α1t-e-α2t)⊗Cp(t),
where α1,2=(k2+k3+k4∓(k2+k3+k4)2-4k2k4)/2and⊗ denotes the mathematical convolution operator.For the ROI-based approach, a fifth parameter of fractional blood volume (fBV) is often included in parameter estimation to address spillover from capillary blood activity within the tissue regions as shown in(2)Ct(t)=(1-fBV)⋅K1α2-α1×[(k4-α1)⋅e-α1t+(α2-k4)⋅e-α2t]⊗Cp(t)+K1k3α2-α1(e-α1t-e-α2t)⊗Cp(t)+fBV⋅Cp(t).The influx rate macroparameter,Ki, represents net influx or uptake of plasma fluoride activity by the bone matrix and has units of mL·min−1·mL−1 and is given by (3)Ki=K1×k3k2+k3.Ki is an important macroparameter for describing the level of osteoblastic activity and is a measure of bone metabolism.


### 2.2. Image Acquisition

The studies and protocols were approved by the Ethics Committee of the Central Sydney Area Health Service. Eight patients with previous limb salvage surgery and bone grafting were studied; there were 3 men and 5 women with age range of 20–53 years. The bone grafts were taken from the fibula. The time between the bone graft surgery and PET scans was ranging from 1.3 to 4.9 years. All imaging studies were carried out on an ECAT 951R whole body PET scanner (Siemens/CTI, Knoxville, Tenn, USA) in Royal Prince Alfred Hospital. The scanner acquired 31 planes simultaneously with a slice separation of 3.375 mm and axial field of view (FOV) of 10.8 cm.206.2±18.1 MBq 18F-sodium fluoride was infused at a constant rate over a 3-minute period. Simultaneously with the beginning of tracer infusion, a 60-minute dynamic sequence of 29 emission scans was started with twelve 10-second, six 30-second, and eleven 5-minute frames. A postinjection transmission method was used to correct for photon attenuation [17].12 arterialised-venous blood samples    [8, 18] were obtained every 30 s for the first 6 min after injection, followed by 4 samples at 1-minute, 2 samples at 5-minute, and 4 samples at 10-minute intervals. Plasma samples were counted in a γ-well-counter cross-calibrated with the PET scanner.Transmission images were reconstructed by an ordered-subset median-root-prior (OS-MRP) iterative reconstruction algorithm with 2 iterations and 16 subsets, which were then segmented into lung, bone, and soft tissue. The dynamic emission images were reconstructed into128×128 matrices by an ordered-subset expectation-maximization (OS-EM) iterative reconstruction technique with 1 iteration and 16 subsets. Attenuation correction was implemented in the OS-EM reconstruction with the attenuation coefficients derived from the segmented reconstructed transmission images.


### 2.3. VOI Delineation

The images from the last six frames, ranging from 30 to 60 min, were averaged to aid definition of VOI, which were defined using the PMOD package (version 3.1, PMOD Technologies Ltd., Zurich, Switzerland) by manually delineating ROIs on voxels with similar values across a sequence of image planes shown in Figure2. VOIs were defined on the bone graft and also on normal bone in the ilium and lumbar vertebra. The intervertebral discs were not included in the VOIs for lumbar vertebra.Figure 2
Example of one transaxial plane from the average images of the last six frames. The arrows indicate manually drawn ROI.Mean value and standard deviation (SD) value were then derived for each VOI on each frame of reconstructed images to form TTAC and SD curve. Figure3 plots typical TTACs and plasma time activity curve (PTAC) of one patient in Figure 3(a) and SD curves over time in Figure 3(b).(a) PTAC and TTACs for graft and ilium. (b) Standard deviation curves.
(a)(b)


### 2.4. Constructing Parametric Images

The Patlak graphical analysis (PGA) assumes a three-compartment and three-parameter model, where the release of fluoride from the bony matrix is considered negligible, that is,k4=0 [9]. This leads to (4) to describe the relations between the input and output function in Figure 1:(4)Ct(t)Cp(t)=K1⋅k3k2+k3⋅∫0tCp(τ)dτCp(t)+k2k2+k3Ce(t)Cp(t).When equilibrium has been reached after a sufficient time(t>t*), (4) can be simplified to (5), where Const represents a constant value:(5)Ct(t)Cp(t)≈Ki⋅∫0tCp(τ)dτCp(t)+Const,t>t*,Ki can then be readily derived from the slope of the linear plot of ∫0tCp(τ)dτ/Cp(t) versus Ct(t)/Cp(t).PGA was used to derive the parametric images ofKi according to (5) by an in-house package implemented in the IDL language (version 6.0, ITT Visual Information Solutions, Boulder, Colo, USA). Four linear periods of the PGA plot reported in the literature were used to derive Ki voxel by voxel for 4 to 60 min [6], 10 to 50 min [12, 15], 20 to 50 min [11], and 20 to 60 min [13, 14]. The value of a voxel was set to zero if the derived Ki was negative.For comparison, the derived VOIs in Section2.3 were used to derive mean and SD values of Ki from the parametric images, respectively, for the same patient. PGA was also applied to derive Ki for VOI-derived TTACs with the same linear periods investigated as in the parametric images.


### 2.5. NLS Analysis

The weights used in NLS analysis were defined by the inverse of the VOI-derived SD curve. The range of parameter variation was set from 0 to 1 for fBV,K1,k2, and k3 and 0 to 0.1 for k4, while the initial parameters were 0.1 mL·mL−1, 0.1 mL·min−1·mL−1, 0.15 min−1, 0.1 min−1, and 0.01 min−1 for fBV, K1,k2,k3, and k4 which were specified    [12]. The generated TTACs were fitted by using the Marquardt-Levenberg algorithm to adjust the rate constant according to (2) using the PMOD package. The maximum number of iterations was set at 200.For the equivalent comparison with the PGA method, NLS analysis was also applied to fit the VOI-derived TTAC with fBV andk4 predefined to be 0.


## 2.1. Kinetic Model for Fluoride Bone Metabolism

A three-compartment and four-parameter model has been used in the kinetic analysis for fluoride bone metabolism (Figure1) [1]. It consists of plasma, extracellular fluid, and bone mineral compartments. Rate constants describe the transport of 18F-fluoride between the compartments with K1 representing the unidirectional clearance of fluoride from plasma to bone tissue, k2 the reverse transport of fluoride from bone tissue to plasma, k3 reflecting the incorporation, and k4 the release of fluoride from the bone mineral compartment. The unit of K1 is mL·min−1·mL−1, while k2,k3, and k4 have units of min−1.Figure 1
Three-compartment and four-parameter kinetic model for fluoride bone metabolism.Cp(t) represents the plasma concentration of 18F-fluoride, while Ce(t) and Cm(t) are 18F-fluoride concentrations in extracellular fluid and bone. Ct(t) is the sum of Ce(t) and Cm(t) and represents total concentration of 18F-fluoride in bone tissue. Through kinetic modelling technique, the output function, Ct(t), can be represented by the input function, Cp(t), in the following:(1)Ct(t)=K1α2-α1[(k4-α1)⋅e-α1t+(α2-k4)⋅e-α2t]⊗Cp(t)+K1k3α2-α1(e-α1t-e-α2t)⊗Cp(t),
where α1,2=(k2+k3+k4∓(k2+k3+k4)2-4k2k4)/2and⊗ denotes the mathematical convolution operator.For the ROI-based approach, a fifth parameter of fractional blood volume (fBV) is often included in parameter estimation to address spillover from capillary blood activity within the tissue regions as shown in(2)Ct(t)=(1-fBV)⋅K1α2-α1×[(k4-α1)⋅e-α1t+(α2-k4)⋅e-α2t]⊗Cp(t)+K1k3α2-α1(e-α1t-e-α2t)⊗Cp(t)+fBV⋅Cp(t).The influx rate macroparameter,Ki, represents net influx or uptake of plasma fluoride activity by the bone matrix and has units of mL·min−1·mL−1 and is given by (3)Ki=K1×k3k2+k3.Ki is an important macroparameter for describing the level of osteoblastic activity and is a measure of bone metabolism.


## 2.2. Image Acquisition

The studies and protocols were approved by the Ethics Committee of the Central Sydney Area Health Service. Eight patients with previous limb salvage surgery and bone grafting were studied; there were 3 men and 5 women with age range of 20–53 years. The bone grafts were taken from the fibula. The time between the bone graft surgery and PET scans was ranging from 1.3 to 4.9 years. All imaging studies were carried out on an ECAT 951R whole body PET scanner (Siemens/CTI, Knoxville, Tenn, USA) in Royal Prince Alfred Hospital. The scanner acquired 31 planes simultaneously with a slice separation of 3.375 mm and axial field of view (FOV) of 10.8 cm.206.2±18.1 MBq 18F-sodium fluoride was infused at a constant rate over a 3-minute period. Simultaneously with the beginning of tracer infusion, a 60-minute dynamic sequence of 29 emission scans was started with twelve 10-second, six 30-second, and eleven 5-minute frames. A postinjection transmission method was used to correct for photon attenuation [17].12 arterialised-venous blood samples    [8, 18] were obtained every 30 s for the first 6 min after injection, followed by 4 samples at 1-minute, 2 samples at 5-minute, and 4 samples at 10-minute intervals. Plasma samples were counted in a γ-well-counter cross-calibrated with the PET scanner.Transmission images were reconstructed by an ordered-subset median-root-prior (OS-MRP) iterative reconstruction algorithm with 2 iterations and 16 subsets, which were then segmented into lung, bone, and soft tissue. The dynamic emission images were reconstructed into128×128 matrices by an ordered-subset expectation-maximization (OS-EM) iterative reconstruction technique with 1 iteration and 16 subsets. Attenuation correction was implemented in the OS-EM reconstruction with the attenuation coefficients derived from the segmented reconstructed transmission images.


## 2.3. VOI Delineation

The images from the last six frames, ranging from 30 to 60 min, were averaged to aid definition of VOI, which were defined using the PMOD package (version 3.1, PMOD Technologies Ltd., Zurich, Switzerland) by manually delineating ROIs on voxels with similar values across a sequence of image planes shown in Figure2. VOIs were defined on the bone graft and also on normal bone in the ilium and lumbar vertebra. The intervertebral discs were not included in the VOIs for lumbar vertebra.Figure 2
Example of one transaxial plane from the average images of the last six frames. The arrows indicate manually drawn ROI.Mean value and standard deviation (SD) value were then derived for each VOI on each frame of reconstructed images to form TTAC and SD curve. Figure3 plots typical TTACs and plasma time activity curve (PTAC) of one patient in Figure 3(a) and SD curves over time in Figure 3(b).(a) PTAC and TTACs for graft and ilium. (b) Standard deviation curves.
(a)(b)


## 2.4. Constructing Parametric Images

The Patlak graphical analysis (PGA) assumes a three-compartment and three-parameter model, where the release of fluoride from the bony matrix is considered negligible, that is,k4=0 [9]. This leads to (4) to describe the relations between the input and output function in Figure 1:(4)Ct(t)Cp(t)=K1⋅k3k2+k3⋅∫0tCp(τ)dτCp(t)+k2k2+k3Ce(t)Cp(t).When equilibrium has been reached after a sufficient time(t>t*), (4) can be simplified to (5), where Const represents a constant value:(5)Ct(t)Cp(t)≈Ki⋅∫0tCp(τ)dτCp(t)+Const,t>t*,Ki can then be readily derived from the slope of the linear plot of ∫0tCp(τ)dτ/Cp(t) versus Ct(t)/Cp(t).PGA was used to derive the parametric images ofKi according to (5) by an in-house package implemented in the IDL language (version 6.0, ITT Visual Information Solutions, Boulder, Colo, USA). Four linear periods of the PGA plot reported in the literature were used to derive Ki voxel by voxel for 4 to 60 min [6], 10 to 50 min [12, 15], 20 to 50 min [11], and 20 to 60 min [13, 14]. The value of a voxel was set to zero if the derived Ki was negative.For comparison, the derived VOIs in Section2.3 were used to derive mean and SD values of Ki from the parametric images, respectively, for the same patient. PGA was also applied to derive Ki for VOI-derived TTACs with the same linear periods investigated as in the parametric images.


## 2.5. NLS Analysis

The weights used in NLS analysis were defined by the inverse of the VOI-derived SD curve. The range of parameter variation was set from 0 to 1 for fBV,K1,k2, and k3 and 0 to 0.1 for k4, while the initial parameters were 0.1 mL·mL−1, 0.1 mL·min−1·mL−1, 0.15 min−1, 0.1 min−1, and 0.01 min−1 for fBV, K1,k2,k3, and k4 which were specified    [12]. The generated TTACs were fitted by using the Marquardt-Levenberg algorithm to adjust the rate constant according to (2) using the PMOD package. The maximum number of iterations was set at 200.For the equivalent comparison with the PGA method, NLS analysis was also applied to fit the VOI-derived TTAC with fBV andk4 predefined to be 0.


## 3. Results

### 3.1. Effect of Parametric Images

Figure4 plots the net influx rates derived from VOI-based TTACs versus those from the parametric images by PGA with 10 to 50 min selected as the linear fitting period. A significant, high correlation (r>0.999) was observed between the values of parametric images and values of VOI-based TTAC for the four PGA approaches.Figure 4
Plots of values ofKi derived by PGA from VOI-based TTAC fitting and values derived from VOI placed on the parametric images. The linear fitting period was 10 to 50 min.Comparably, low values of SD were observed in all the regions for parametric image ofKi (Table 1). This indicated that PGA provided reliable and consistent estimates of Ki for region averaged TTACs as well as voxel-by-voxel TTACs.Table 1
Summary of standard deviation (mean ± SD) in all the regions for parametric images ofKi.
Linear Period4–60 min10–50 min20–50 min20–60 minMean ± SD0.004±0.0020.005±0.0030.007±0.0040.005±0.003Figure5 shows the parametric images derived by PGA approach for the investigated linear fitting periods. More detail is achieved by the PGA with the linear period of 20–50 min (Figure 5(c)) and 20–60 min (Figure 5(d)). However, this comes at the expense of increased noise and slightly less reliable parameter estimation reliability.Figure 5
Parametric images ofKi derived by PGA with the linear portion of (a) 4–60 min, (b) 10–50 min, (c) 20–50 min, and (d) 20–60 min, for the transaxial plane of the average image shown in Figure 2.


### 3.2. Rate Constants of NLS Analysis

Table2 lists the rate constants derived by NLS for the three-compartment and five-parameter model as well as the relevant influx rate (Ki-5p). The mean values of parametric images of Ki derived by PGA are also listed for the comparison. The influx rate for three-compartment and three-parameter model with fBV = 0 and k4=0 is also given in the table and is referred to as Ki-3p.Table 2
Kinetic parameters for bone metabolism.
No.fBVK1k2k3k4Ki-5pKi-3pKi4–60Ki10–15Ki20–50Ki20–60Graft10.0000.0370.1180.0910.0160.0160.002#0.0140.0160.0160.01220.0000.0510.3890.1200.0000.0120.0120.0130.0130.0130.01230.0000.0180.0600.0980.0000.0110.0120.0120.0110.0110.01240.0840.0130.000#0.4530.0060.0130.0120.0120.0120.0120.01150.0420.0160.000#1.000#0.1000.0160.0160.0160.0160.0170.01860.0000.0440.0640.0480.0000.0190.0200.0220.0210.0220.02170.0000.0120.000#0.000#0.000#0.000#0.0110.0110.0120.0120.01080.0000.0570.0300.0010.0000.002#0.0270.0290.0260.0260.02790.0000.0070.000#0.9100.0900.007#0.007#0.0220.0230.0220.015100.0000.0140.000#0.000#0.0000.006#0.013#0.0230.0220.0200.018110.0000.0700.8440.1360.0060.0100.0080.0080.0090.0100.008120.0090.1200.4440.0440.0080.0110.0090.0080.0080.0090.007130.0000.0320.0540.0370.0040.0130.0120.0150.0150.0150.014Ilium140.0000.1720.4860.1350.0070.0370.0330.0330.0320.0300.030150.0000.1250.2290.0780.0040.0320.0300.0330.0330.0340.031160.0000.2480.3020.0620.0060.0420.0380.0450.0440.0410.037170.0000.4080.8950.1220.0040.0490.0440.0450.0480.0470.041180.0000.1800.2200.0740.0070.0460.0400.0440.0420.0400.039190.0000.2340.5920.1700.0140.0520.0390.0400.0410.0340.029200.0480.0670.0790.1260.0150.0420.0350.0350.0350.0340.032210.0000.1450.2120.1310.0060.0550.0500.0520.0490.0480.049220.0020.1360.4840.1650.0030.0340.0320.0320.0320.0320.032230.0630.1100.0690.0740.0080.0570.0520.0520.0500.0480.047240.0000.0120.000#0.9990.0000.012#0.011#0.0200.0180.0180.020250.0240.0671.000#0.6090.0090.0250.0220.0200.0200.0180.018260.0000.0750.1970.2160.0130.0390.0290.0330.0350.0350.030270.0000.0240.0160.0230.0000.0140.0140.0170.0170.0170.016280.0000.0360.3120.0980.0040.0090.0080.0080.0080.0080.008290.0000.0390.3370.1700.0130.0130.0100.0110.0100.0090.010300.0000.0580.8010.1670.0080.0100.0080.0090.0080.0080.008310.0000.0350.3260.1780.0110.0120.0100.0110.0100.0100.010320.0000.1080.1710.0700.0000.0310.0320.0340.0320.0320.032330.0000.0240.000#0.0900.0900.024#0.014#0.0350.0380.0380.032Lumbar vertebra340.0650.0900.1060.0480.0000.0280.0270.0290.0230.0210.026350.0030.1630.2260.1210.0110.0570.0460.0480.0460.0430.043360.0000.1460.1750.0900.0050.0490.0420.0500.0470.0450.047370.0000.1030.2830.1630.0100.0380.0310.0310.0320.0330.028380.0000.0370.5130.2010.0140.0100.0080.0080.0080.0070.007390.0000.0400.4910.2010.0130.0120.0090.0090.0090.0080.008#Unsuccessful fit.Some unsuccessful fits were observed by NLS analysis such as whenk2 was close to zero or one. This is not unexpected because the success of physiological meaningful NLS analysis is dependent on not only appropriate initial parameters but also proper noise model and underlying kinetic model. Insufficient data quality in TTAC may give rise to such unsuccessful fit for the three-compartment and five-parameter kinetic model. Interestingly, most of estimated fBV is observed to be zero, which implies that the spillover contribution from surrounding vascular system can be ignored.Table3 lists the mean values and SD for all the VOIs with unsuccessful fits excluded. The lowest value of K1 was observed for the bone grafts, which implies reduced blood flow. Similar values of k2,k3, and k4 were observed in the graft, ilium, and lumbar vertebra. The lower K1 of bone graft in contrast to normal bone suggests that reduced osteoblastic activity, indicated by reduced Ki, of the bone graft was mainly caused by poor vascularisation and reduced blood flow. This is consistent with the findings in    [13].Table 3
Mean and standard deviation (mean ± SD) of derived rate constants for successful NLS analysis.
fBVK1k2k3k4Flux-KiGraft0.001±0.0030.053±0.0340.282±0.2970.082±0.0390.005±0.0060.013±0.003Ilium0.007±0.0030.129±0.0760.337±0.2760.121±0.0420.007±0.0030.034±0.014Lumbar vertebra0.011±0.0260.097±0.0520.299±0.1680.137±0.0620.009±0.0050.032±0.019


### 3.3. Comparison of NLS and PGA for Parametric Image

Linear regression analysis was applied to compare the accuracy of the net influx rate derived by PGA compared with NLS for the three-compartment and five-parameter model with unsuccessful fits excluded. Table4 lists the results of linear regression for the studied four PGA linear periods. One example of linear regression analysis is plotted in Figure 6 for the PGA with the linear period of 20 to 60 min.Table 4
Linear regression between NLS and PGA for parametric images.
Linear portion4–60 min10–50 min20–50 min20–60 minr0.9790.9830.9700.961y = ax + b0.890×Ki-5p+0.0010.876×Ki-5p+0.0010.818×Ki-5p+0.0020.787×Ki-5p+0.002Figure 6
Linear regression analysis between NLS for VOI-based TTAC and PGA for parametric image with linear period of 20 to 60 min. Black line: the plot of obtained linear regression; grey line: theoretical ideal linear regression.A high linear correlation was observed forKi derived from PGA parametric images compared with NLS with the linear correlation r≥0.961. However, PGA tends to underestimate Ki with the slope of linear analysis ranging from 0.787 to 0.890 as compared with Ki-5p derived by NLS. The values of Ki-3p in Table 2 are almost identical to the values of original Ki-5p when fitted k4 is 0, while Ki-3p is lower than Ki-5p when fitted k4 is >0. The same underlying model is assumed in deriving Ki-3p for NLS and Ki for PGA. Ignoring the contribution of k4 is likely to be the main reason for the underestimation of the influx rate with PGA.The potential for underestimatingKi with PGA should be recognised when interpreting the kinetic analysis of bone metabolism. However, the simplicity of PGA and the benefit of parametric images make PGA a more attractive approach for the quantitative analysis of bone metabolism. Overall, taking reliability, accuracy, and linear regression analysis into account, the PGA analysis with the linear period of 10 to 50 min should provide best parametric images of Ki in the quantitative analysis of bone metabolism.


## 3.1. Effect of Parametric Images

Figure4 plots the net influx rates derived from VOI-based TTACs versus those from the parametric images by PGA with 10 to 50 min selected as the linear fitting period. A significant, high correlation (r>0.999) was observed between the values of parametric images and values of VOI-based TTAC for the four PGA approaches.Figure 4
Plots of values ofKi derived by PGA from VOI-based TTAC fitting and values derived from VOI placed on the parametric images. The linear fitting period was 10 to 50 min.Comparably, low values of SD were observed in all the regions for parametric image ofKi (Table 1). This indicated that PGA provided reliable and consistent estimates of Ki for region averaged TTACs as well as voxel-by-voxel TTACs.Table 1
Summary of standard deviation (mean ± SD) in all the regions for parametric images ofKi.
Linear Period4–60 min10–50 min20–50 min20–60 minMean ± SD0.004±0.0020.005±0.0030.007±0.0040.005±0.003Figure5 shows the parametric images derived by PGA approach for the investigated linear fitting periods. More detail is achieved by the PGA with the linear period of 20–50 min (Figure 5(c)) and 20–60 min (Figure 5(d)). However, this comes at the expense of increased noise and slightly less reliable parameter estimation reliability.Figure 5
Parametric images ofKi derived by PGA with the linear portion of (a) 4–60 min, (b) 10–50 min, (c) 20–50 min, and (d) 20–60 min, for the transaxial plane of the average image shown in Figure 2.


## 3.2. Rate Constants of NLS Analysis

Table2 lists the rate constants derived by NLS for the three-compartment and five-parameter model as well as the relevant influx rate (Ki-5p). The mean values of parametric images of Ki derived by PGA are also listed for the comparison. The influx rate for three-compartment and three-parameter model with fBV = 0 and k4=0 is also given in the table and is referred to as Ki-3p.Table 2
Kinetic parameters for bone metabolism.
No.fBVK1k2k3k4Ki-5pKi-3pKi4–60Ki10–15Ki20–50Ki20–60Graft10.0000.0370.1180.0910.0160.0160.002#0.0140.0160.0160.01220.0000.0510.3890.1200.0000.0120.0120.0130.0130.0130.01230.0000.0180.0600.0980.0000.0110.0120.0120.0110.0110.01240.0840.0130.000#0.4530.0060.0130.0120.0120.0120.0120.01150.0420.0160.000#1.000#0.1000.0160.0160.0160.0160.0170.01860.0000.0440.0640.0480.0000.0190.0200.0220.0210.0220.02170.0000.0120.000#0.000#0.000#0.000#0.0110.0110.0120.0120.01080.0000.0570.0300.0010.0000.002#0.0270.0290.0260.0260.02790.0000.0070.000#0.9100.0900.007#0.007#0.0220.0230.0220.015100.0000.0140.000#0.000#0.0000.006#0.013#0.0230.0220.0200.018110.0000.0700.8440.1360.0060.0100.0080.0080.0090.0100.008120.0090.1200.4440.0440.0080.0110.0090.0080.0080.0090.007130.0000.0320.0540.0370.0040.0130.0120.0150.0150.0150.014Ilium140.0000.1720.4860.1350.0070.0370.0330.0330.0320.0300.030150.0000.1250.2290.0780.0040.0320.0300.0330.0330.0340.031160.0000.2480.3020.0620.0060.0420.0380.0450.0440.0410.037170.0000.4080.8950.1220.0040.0490.0440.0450.0480.0470.041180.0000.1800.2200.0740.0070.0460.0400.0440.0420.0400.039190.0000.2340.5920.1700.0140.0520.0390.0400.0410.0340.029200.0480.0670.0790.1260.0150.0420.0350.0350.0350.0340.032210.0000.1450.2120.1310.0060.0550.0500.0520.0490.0480.049220.0020.1360.4840.1650.0030.0340.0320.0320.0320.0320.032230.0630.1100.0690.0740.0080.0570.0520.0520.0500.0480.047240.0000.0120.000#0.9990.0000.012#0.011#0.0200.0180.0180.020250.0240.0671.000#0.6090.0090.0250.0220.0200.0200.0180.018260.0000.0750.1970.2160.0130.0390.0290.0330.0350.0350.030270.0000.0240.0160.0230.0000.0140.0140.0170.0170.0170.016280.0000.0360.3120.0980.0040.0090.0080.0080.0080.0080.008290.0000.0390.3370.1700.0130.0130.0100.0110.0100.0090.010300.0000.0580.8010.1670.0080.0100.0080.0090.0080.0080.008310.0000.0350.3260.1780.0110.0120.0100.0110.0100.0100.010320.0000.1080.1710.0700.0000.0310.0320.0340.0320.0320.032330.0000.0240.000#0.0900.0900.024#0.014#0.0350.0380.0380.032Lumbar vertebra340.0650.0900.1060.0480.0000.0280.0270.0290.0230.0210.026350.0030.1630.2260.1210.0110.0570.0460.0480.0460.0430.043360.0000.1460.1750.0900.0050.0490.0420.0500.0470.0450.047370.0000.1030.2830.1630.0100.0380.0310.0310.0320.0330.028380.0000.0370.5130.2010.0140.0100.0080.0080.0080.0070.007390.0000.0400.4910.2010.0130.0120.0090.0090.0090.0080.008#Unsuccessful fit.Some unsuccessful fits were observed by NLS analysis such as whenk2 was close to zero or one. This is not unexpected because the success of physiological meaningful NLS analysis is dependent on not only appropriate initial parameters but also proper noise model and underlying kinetic model. Insufficient data quality in TTAC may give rise to such unsuccessful fit for the three-compartment and five-parameter kinetic model. Interestingly, most of estimated fBV is observed to be zero, which implies that the spillover contribution from surrounding vascular system can be ignored.Table3 lists the mean values and SD for all the VOIs with unsuccessful fits excluded. The lowest value of K1 was observed for the bone grafts, which implies reduced blood flow. Similar values of k2,k3, and k4 were observed in the graft, ilium, and lumbar vertebra. The lower K1 of bone graft in contrast to normal bone suggests that reduced osteoblastic activity, indicated by reduced Ki, of the bone graft was mainly caused by poor vascularisation and reduced blood flow. This is consistent with the findings in    [13].Table 3
Mean and standard deviation (mean ± SD) of derived rate constants for successful NLS analysis.
fBVK1k2k3k4Flux-KiGraft0.001±0.0030.053±0.0340.282±0.2970.082±0.0390.005±0.0060.013±0.003Ilium0.007±0.0030.129±0.0760.337±0.2760.121±0.0420.007±0.0030.034±0.014Lumbar vertebra0.011±0.0260.097±0.0520.299±0.1680.137±0.0620.009±0.0050.032±0.019


## 3.3. Comparison of NLS and PGA for Parametric Image

Linear regression analysis was applied to compare the accuracy of the net influx rate derived by PGA compared with NLS for the three-compartment and five-parameter model with unsuccessful fits excluded. Table4 lists the results of linear regression for the studied four PGA linear periods. One example of linear regression analysis is plotted in Figure 6 for the PGA with the linear period of 20 to 60 min.Table 4
Linear regression between NLS and PGA for parametric images.
Linear portion4–60 min10–50 min20–50 min20–60 minr0.9790.9830.9700.961y = ax + b0.890×Ki-5p+0.0010.876×Ki-5p+0.0010.818×Ki-5p+0.0020.787×Ki-5p+0.002Figure 6
Linear regression analysis between NLS for VOI-based TTAC and PGA for parametric image with linear period of 20 to 60 min. Black line: the plot of obtained linear regression; grey line: theoretical ideal linear regression.A high linear correlation was observed forKi derived from PGA parametric images compared with NLS with the linear correlation r≥0.961. However, PGA tends to underestimate Ki with the slope of linear analysis ranging from 0.787 to 0.890 as compared with Ki-5p derived by NLS. The values of Ki-3p in Table 2 are almost identical to the values of original Ki-5p when fitted k4 is 0, while Ki-3p is lower than Ki-5p when fitted k4 is >0. The same underlying model is assumed in deriving Ki-3p for NLS and Ki for PGA. Ignoring the contribution of k4 is likely to be the main reason for the underestimation of the influx rate with PGA.The potential for underestimatingKi with PGA should be recognised when interpreting the kinetic analysis of bone metabolism. However, the simplicity of PGA and the benefit of parametric images make PGA a more attractive approach for the quantitative analysis of bone metabolism. Overall, taking reliability, accuracy, and linear regression analysis into account, the PGA analysis with the linear period of 10 to 50 min should provide best parametric images of Ki in the quantitative analysis of bone metabolism.


## 4. Discussion

18F-fluoride PET with quantitative kinetic analysis has allowed the quantitative assessment of regional lesions and treatment response in metabolic bone disease and also enabled early identification of bone viability and discriminated uneventful and impaired healing processes of fractures, bone grafts, and osteonecrosis    [19]. Despite its desirable characteristics of high and rapid bone uptake with very rapid blood clearance, 18F-fluoride PET scanning has yet to become part of routine clinical practice    [20].The Gaussian noise model is usually assumed in the NLS fitting with the weights defined according to the activity of TTAC and the duration of imaging frame. In this study, the weight is defined by the inverse of the VOI-derived standard deviation instead, which provides an approximate estimation of noise distribution for the VOI-derived TTAC. However, due to the low signal-to-noise ratio in early, short imaging frames, the iterative reconstruction caused the values of voxels with low activity to be estimated to be zero and explains why the values of SD are shown to be zero for the early frames in Figure3(b). Thus, the extraordinary high weights assigned for the early zero points would give rise to unsuccessful NLS fits, especially for the lower 18F-fluoride uptake in the grafts. With the first sixteen data points (first 4 min) excluded, NLS analysis was found to achieve successful fitting for VOI-based TTACs as summarised in Table 2. For example, the rate constants (region no. 7, graft region in Figure 3) are 0.000 for fBV, 0.025 for K1, 0.179, 0.173, and 0.006, respectively, for k2,k3, and k4 while Ki-5p is 0.012, which is identical to the mean value of influx rate derived from the corresponding parametric images VOI.The estimation ofK1 in NLS analysis was sensitive to data points in early imaging time. The obtained zero points in the first 4-minute period may lead to biased estimation of K1. However, the parameter of interest, Ki, was insensitive to the early imaging time. The NLS analysis was expected to provide reasonable estimate of Ki as the reference in the evaluation of PGA method.The late time period inherently used for the PGA fit avoids the problem of the low count early frames for the PGA method. We observed high reliability for the studied PGA with four linear periods in the generated parametric images ofKi, which are highly correlated with NLS-derived Ki values. Interestingly, the underestimation of PGA-derived Ki tends to be higher with later starting times of the linear period. Thus, the inclusion of early staring time in PGA analysis is necessary to reduce the underestimation of Ki due to the simplifying assumptions of the underlying kinetic model.Hybrid PET/CT scanner has now become widespread with improved image quality and shorter imaging time with coregistered CT images. Combining with the detailed anatomical information readily available from CT images, the state-of-the-art PET/CT imaging may lead to more accurate assessment of graft viability. This requires in-depth investigation on improving the accuracy of CT-based attenuation correction for PET scans, especially for bone and graft tissues, to take advantage of the benefits offered by the hybrid scanners.


## 5. Conclusions

Dynamic18F-fluoride PET imaging studies in patients with limb salvage surgery and fibula bone grafts were analysed using NLS and PGA methods to derive TTACs and parametric images. 39 regions were analysed respectively for the bone graft, ilium, and lumbar vertebrae. The results show that parametric images derived by PGA are consistent with VOI-based parameter estimation by PGA with high reliability. The PGA approach tended to underestimate influx rate because the assumption that k4=0 was not valid in all cases. The linear portion for PGA analysis was suggested to be from 10 to 50 min. This quantitative approach using dynamic 18F-fluoride PET imaging allows the assessment of bone metabolism, and additional clinical studies will clarify its role in identifying bone graft viability.



---

*Source: 189830-2011-05-22.xml*