Partial volume correction for improved PET quantification in 18F-NaF imaging of atherosclerotic plaques

Background Accurate quantification of plaque imaging using 18F-NaF PET requires partial volume correction (PVC). Methods PVC of PET data was implemented by the use of a local projection (LP) method. LP-based PVC was evaluated with an image quality (NEMA) and with a thorax phantom with “plaque-type” lesions of 18-36 mL. The validated PVC method was then applied to a cohort of 17 patients, each with at least one plaque in the carotid or ascending aortic arteries. In total, 51 calcified (HU > 110) and 16 non-calcified plaque lesions (HU < 110) were analyzed. The lesion-to-background ratio (LBR) and the relative change of LBR (ΔLBR) were measured on PET. Results Following PVC, LBR of the spheres (NEMA phantom) was within 10% of the original values. LBR of the thoracic lesions increased by 155% to 440% when the LP-PVC method was applied to the PET images. In patients, PVC increased the LBR in both calcified [mean = 78% (−8% to 227%)] and non-calcified plaques [mean = 41%, (−9%-104%)]. Conclusions PVC helps to improve LBR of plaque-type lesions in both phantom studies and clinical patients. Better results were obtained when the PVC method was applied to images reconstructed with point spread function modeling. Electronic supplementary material The online version of this article (doi:10.1007/s12350-017-0778-2) contains supplementary material, which is available to authorized users.


INTRODUCTION
Cardiovascular diseases (CVDs) are frequently caused by the rupture of a vulnerable atherosclerotic plaque, thus, resulting in thrombotic occlusion or distal embolization. 1,2 Early calcium deposits are common constituents of atherosclerotic plaques, and, therefore, serve as a surrogate marker for atherosclerosis. 3 CT imaging has been used to assess the severity of vascular calcification and the atherosclerotic plaque burden. 4 However, CT calcium scores do not represent the pathophysiologic behavior of plaques and, therefore, do not correlate well with cardiovascular risk factors. 5 In view of the similarities between calcification extent and osteogenesis, 3 18 F-NaF uptake in PET has been described as a marker of ongoing calcium deposition in vulnerable plaque in the carotid arteries. [6][7][8][9][10] These studies, however, showed significant discrepancies between PET and CT findings. While the areas with positive PET uptake correlated with the CT-based calcification, only a fraction of the calcifications identified on CT images presented a significant 18 F-NaF uptake. Recently, Fiz et al. demonstrated an inverse correlation between calcification and 18 F-NaF uptake of plaques in the carotid arteries in a cohort of 64 patients. 11 These studies share one potential limitation. Due to the small size of the vulnerable plaques, which is comparable to the typical spatial resolution of PET, partial volume effects (PVE) may arise and compromise the quantification accuracy of 18 F-NaF PET uptake in plaques.
Numerous approaches have been proposed for partial volume correction (PVC) in nuclear medicine applications. [12][13][14][15][16] These techniques can be divided into two main categories: within-reconstruction and postreconstruction methods, 12 which include projectionbased methods, such as the local projection (LP) method proposed for SPECT 17,18 and for pre-clinical PET. 19 While numerous studies have demonstrated the importance of PVC in neurology, psychiatry, oncology, or cardiology, 16 few studies have addressed this question regarding PET quantification of vulnerable atherosclerotic plaques, while merely focusing on [ 18 F]fluorodeoxyglucose (FDG). 2,20 More recently, Huet et al. stated the importance of implementing advanced PVC methods for a proper assessment of plaque inflammation or calcification with PET. 21 The aim of this work was to evaluate a PVC method for clinical PET imaging of atherosclerotic plaques. More specifically, this was done by validating the LP method 17,18 on phantom scans and subsequently applying it to clinical 18 F-NaF PET images of patients with plaques in the carotid or aortic arteries.

Partial Volume Correction Method
Consider J different tissue compartments, including the lesion of interest, within a small VOI delineated in the object as shown in Figure 1a (for J = 2). The measured emission projection counts for each sinogram bin, k i , can be modeled as the sum of the projection counts from each of the j-segmented tissue (with unit activity concentration), scaled by their respective activity, plus the counts originating from the global background outside the VOI, as shown in Eq. (1): where k i are the expected counts per sinogram bin, i, A j the activity for each segmented tissue j within the VOI, P ij is the resolution-blurred tissue shape function for tissue j and sinogram bin i, and g out;i represents the background counts coming from the region outside the VOI in each sinogram bin i. Note that the indices i and j increase through the sinogram bins in the projected space and the segmented tissues in the image space, respectively. The tissue activities A j can be determined by fitting the measured projection data to the model in Eq. (1). Further details can be found in Appendix 1. These tissue activities (A j ) can be used as a prior in a new reconstruction procedure, which is in this study performed using the software for tomographic image reconstruction (STIR). 22 This step leads to a locally ''partial volume corrected (PVC)'' PET image, with improved visual appearance and quantification of the lesions of interest. A schematic view of the procedure is presented in Figure 1b.
• First, the system matrix values (P ij ) for each tissue and the global background g out,i that affects the VOI activity estimates are computed.
• Later, the LP activities for each segmented tissue are obtained, and the activity in each voxel within the VOI is substituted with the tissue activities obtained with the LP method.
• The resulting image is forward-projected using the STIR projector. A realistic physical modeling of the system is included within the projector, with attenuation and scatter estimates from the anatomical CT image and realistic Poisson noise (taking into account the activity concentration and the acquisition time) included into the forward projection. The scatter was estimated using a version of the Single Scatter Simulation (SSS) algorithm available in STIR. 22 This step results in realistic ''simulated'' projection data for the given acquisition.
• The resulting sinogram is reconstructed using a conventional ''Ordered Subsets Expectation Maximization'' (OSEM) algorithm, 23 available within the STIR library, with 5 iterations and 21 subsets. No further changes of the LP activity values were made in the new reconstruction procedure. Scatter and attenuation corrections were included into the iterative reconstruction algorithm as additive and multiplicative terms to the estimated data, respectively. This additional reconstruction procedure yields the PVC image.
The above-described procedure for PVC was employed in all the phantom and patient studies. In all cases, the standard PET images were obtained using the vendor software. Two different reconstruction algorithms were used for the reconstruction of the standard PET images: ordered subsets expectation maximization (OSEM) algorithm 23 and an OSEM reconstruction with point spread function modeling (PSF). 24 All the relevant physical corrections (attenuation, scatter, normalization, decay, dead time) were included in the vendor OSEM and PSF algorithms. After reconstruction, the LP method for PVC was applied, and the LP activity values were used in an additional STIR reconstruction in order to obtain the PVC image.

Phantom Evaluations
The clinical implementation of the LP method was evaluated using acquisitions of a NEMA NU2-2012 IQ phantom 25 and a human-sized thorax phantom. 26 The NEMA phantom contains six fillable spheres with internal diameters of 10, 13, 17, 22, 28, and 37 mm (Figure 2a) with an experimental lesion-to-background ratio (LBR) of 4.95, which is in accordance with the LBR values recommended by the NEMA NU-2 2012 protocol for the measurement of image quality (two acquisitions with LBR 4:1 and 8:1, respectively 33 ). The thorax phantom has three spherical ''plaque-type'' lesions of 36, 31, and 18 mm 3 inserted ( Figure 2b) with a LBR of 70:1, following the values suggested in Delso et al. 27 NEMA IQ phantom data and the thorax phantom were acquired using a Biograph true-point true-view (TPTV) PET/ CT 28 and Biograph mCT PET/CT system, respectively. 29 The parameters for the acquisition and the reconstruction of the data are summarized in Table 1. OSEM and PSF algorithms, both available from the vendor software, were used for the reconstruction of the acquired images.

F-NaF PET/CT Patient Studies
The LP-based PVC method ( Figure 1) was applied retrospectively to a cohort of patients with multiple myeloma, who underwent 18 F-NaF PET/CT whole-body imaging to characterize bone lesions. In this work, we evaluated the 18 F-NaF PET uptake in the carotid or ascending aortic arteries. In Figure 1. Illustration of the procedure employed to improve the quantification of the reconstructed image within the VOI using the activities computed with the LP method. First, the LP method is applied to obtain the PVC tissue activities of each segmented tissue j and voxel v within the VOI. Then, the original activities in each voxel within the VOI, IMG (v), are substituted with the LPbased tissue activities, A j (v). The resulting image is forward-projected to obtain a simulated sinogram, which is then reconstructed with STIR, yielding a simulated PVC image.
total, 17 patients (12 male, 5 female, mean age: (64 ± 9) years, range: (47-77) years) with at least one positive plaque were analyzed. Plaques were classified as calcified (HU [ 110 within the plaque) and non-calcified (HU \ 110), and calcification was defined as the area with a minimum density of 110 HU on CT. Given the fact that a low-dose attenuationcorrected CT (AC_CT) image was used for anatomical reference, a comparatively low HU value, 10,11,30 was used as a threshold level for the definition of calcified plaque. In total, 51 calcified (HU [ 110 within the plaque) and 16 noncalcified plaque lesions (HU \ 110) were analyzed. This retrospective study was approved by the Institutional Ethics Committee and was in accordance with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
All scans were performed on a Biograph TPTV system. 27 The patients were injected with (4.3 ± 1.0) MBq/kg (range: 3.1-6.0 MBq/kg) of 18 F-NaF. The post-injection delay interval was (56 ± 12) minute (range: 34-70 minute), and the PET acquisition time was 2 minute per bed position. 3D PET data were reconstructed using a PSF reconstruction with resolution modeling available from the vendor (4 iterations, 21 subsets). One-bed position image centered in the head-neck region, with 336 9 336 9 109 matrix size and 2.03 9 2.03 9 2.03 mm 3 voxel size, was obtained for each patient.

Segmentation Methods
In the cases of calcified plaques and NEMA spheres, the segmentation of the local VOI was performed using a low-dose CT image (512 9 512 9 109 voxels and 1.37 9 1.37 9 2.03 mm 3 per voxel) co-registered to the PET image. A threshold-based segmentation of the PET image was performed in non-calcified plaque lesions and the thorax  The voxel size used in the reconstructions was 2.03 9 2.03 9 2.03 mm 3 in all cases phantom data. This segmentation was made using the 3D isocontour half-way between the maximum voxel activity and the mean background activity, as defined by Boellaard et al. 31 The dependence of the PVC results on the method used for segmentation was evaluated with NEMA IQ phantom acquisitions.

Data Analysis
Phantom data The quantification accuracy was evaluated by measuring the LBR and the relative change of LBR (DLBR) after applying the PVC method, for each hot sphere and for each plaque-type lesion. The LBR for each hot lesion was measured using the maximum voxel activity within the sphere (LBR max ) and the mean activity within a 3D isocontour at 50% of the maximum voxel activity adapted to the mean background activity (LBR A50 ). The DLBR was calculated as where LBR(PVC) is the lesion-to-background ratio after applying the PVC (measured from the tissue activities obtained with the LP method or from the PVC image) and LBR(noPVC) is the lesion-to-background ratio measured in the image reconstructed with the vendor software (OSEM or PSF). The activity of the background region, on other hand, was determined by drawing several VOIs in uniform regions. Patient data The maximum HU value within the plaque was evaluated from the CT images. All plaque lesions were . OSEM (top) and PSF (center) reconstructions of the NEMA IQ phantom. From left to right: standard images reconstructed with the vendor software, PVC image reconstructed with STIR software using the LP tissue activities corresponding to the 17-mm-diameter sphere, the 13-mm sphere, and the smallest 10-mm sphere. Bottom: Activity profiles through the 17-, 13-, and 10-mmdiameter spheres, as depicted in the PVC images. These profiles were obtained from images reconstructed with an OSEM algorithm (blue), with (dashed line) and without (solid line) local PVC, and from PSF images (black), with (dashed) and without (solid) PVC. The measured activity for the spheres was 23.8 kBq/mL, and is marked as a EXP in the plots. The PVC images (dashed lines) showed activity values for the spheres close to the measured ones.  For calcified plaques, the volume was determined from the CT image by defining all the voxels within the plaque. For noncalcified plaques, we measured the volume of the segmented plaque lesion from the PET threshold-based segmentation. From the PET images, the chosen figures of merit were the LBR using the maximum (LBR max ) and mean (LBR mean ) pixel values within the segmented plaque lesion, and the DLBR after applying the partial volume correction [Eq. (2)]. The background region was depicted in an arterial region where neither calcium deposition nor increased 18 F-NaF uptake was detected. This region was placed between 10 mm and 15 mm below the location of the plaque lesion. The LBR and DLBR for each group of patients were reported as mean ± SD. The dependence of the LBR max and DLBR max on the density of the plaque (in HU) and the segmented volume was evaluated by fitting the data to the following analytical expression: where the dependent variable ''y'' is LBR or DLBR and the independent variable ''x'' is the density or the segmented volume of the plaque. A and B are the fitting parameters. The first term of the equation represents the non-linear dependence close to x = 0 (small volume or low HU value of the plaque), while the second term represents the uniform LBR or DLBR values observed for large, heavy plaques. Pearson and Spearman coefficients were evaluated to test correlations between the measured variables, and a one-sided paired t test was used to evaluate statistical significant changes in the LBR max and DLBR max values obtained with the PSF, PSF ? PVC, and LP methods. In 6 of 17 patients, a small misalignment was observed between the PET and the CT image volumes in the carotid region. In these cases, we performed an additional manual finetuning of spatial alignment of the images following a rigid affine translation. Figure 3 shows the comparison of the standard OSEM and PSF reconstructions with their corresponding PVC images, for the acquisitions of the NEMA IQ phantom. The local PVC performance is illustrated for the 10-, 13-, and 17-mm spheres. The activity profiles across these spheres are shown in Figure 3(bottom). The quantification of the spheres was significantly improved when the PVC method was applied. Table 2 summarizes the performance of the PVC method by means of the dependence of the LBR max and LBR A50 values on the sphere size (for the NEMA IQ phantom). The LBR values were obtained from the standard images (OSEM and PSF columns) and from the PVC images obtained using the CT-based (OSEM ? PVC -CT, PSF ? PVC -CT) and the PET-based (OSEM ? PVC -PET, PSF ? PVC -PET) segmentation approaches. The three columns to the right provide the sphere-to-background ratio obtained using the segmented tissue activities obtained with the LP method (LP -CT and LP -PET, see Figure 1B step 2) and the experimental values measured in a well counter (EXP). The quantification accuracy for the spheres was significantly improved when PVC was applied, yielding LBR max and LBR A50 values close to the reference value (EXP). The performance of the PVC method for the plaquetype lesions in the thorax phantom is summarized in Table 3. Figure 4 shows coronal views of OSEM and OSEM ? PVC reconstructions of the thorax phantom with the three plaque-type lesions. Similar results were obtained for the PSF reconstructions. A significant increase in the LBR max and LBR A50 values was observed when PVC was applied to the images (DLBR max values between 115% and 328% for the three plaque lesions). Furthermore, even larger DLBR max values than the ones reported above were observed when comparing the LBR values in the uncorrected images with the tissue activities obtained from the LP method (values between 155% and 475%). Figure 5 shows transaxial images of patients with calcified and non-calcified plaque in the carotids. After PVC, both plaque uptake and delineation of the calcified plaques improve. More specifically, a significant increase of the LBR max was observed in both calcified [mean = 78%, (-8% to 227%)] and non-calcified plaques [mean = 41%, (-9% to 104%)], when the LP method was applied. The relation between LBR max and DLBR max with the plaque segmented volume is presented in Figure 6a, b. As expected, the DLBR max increases when the volume of the plaque decreases. Figure 6c, d shows the dependence of the LBR max and DLBR max on the density of the plaque for calcified plaque lesions. In that case, the LBR max does not demonstrate a significant dependence on plaque density in the absence of PVC. However, we observe a higher LBR for lighter plaques for the case when PVC is applied. This observation is reinforced when the identified plaques are classified into the four groups mentioned above (see materials and methods, data analysis), as can be seen in Tables 4 and 5. Note that the empirical fits in Figure 6 are presented solely to guide the eye of the readers; they do not imply a theoretical dependence, following the fitted function, of the LBR and DLBR values with the volume or the HU of the plaque. Table 6 shows the fitting parameters for all the fits in Figure 6 and the Pearson and Spearman correlation coefficients between the measured variables. As expected, the fitting parameter A for the LBR max vs HU curves was close to zero when no PVC was applied, demonstrating that the LBR does not depend significantly on the HU value in the PSF images. In contrast, for the PSF ? PVC and LP methods, a significant dependence was observed. A similar behavior was found in the LBR max vs volume curves. Positive correlations were found between the LBR/DLBR and the segmented volume or HU of the plaque in the PSF images. However, for the PSF ? PVC images and the LP method, the correlation was found to be negative and significantly higher for the LP method than for the PSF ? PVC images. The higher values for the Spearman correlation factor in the LP curves confirm the non-linear behavior of these curves. Significant differences between the PSF and the PSF ? PVC mean LBR were found (t-statistic = -12.3, P = 5 9 10 -19 ) and between the PSF ? PVC and the LP mean LBR values (t-statistic = -6.1, P = 3 9 10 -8 ). Figure 6. a Lesion-to-background ratio (LBR) obtained using the maximum pixel value within the plaque versus the volume of the plaque. b Activity recovered after applying the PVC method versus volume of the plaque. c Lesion-to-background ratio (LBR) obtained using the maximum pixel value within the plaque versus the HU value of the plaque (only calcified plaques). d DLBR after applying PVC versus the HU of the plaque (only calcified plaques). The empirical fits of the data to the function y ¼ A=x þ B, for each of the evaluated methods, are also shown. Note, the segmented volume in panels A and B is represented using a logarithmic scale.

DISCUSSION
In this work, we assess the ability and usefulness of PVC for PET imaging of plaque-type lesions. Based on phantom and patient data, we are able to demonstrate that the LBR of plaque-type lesions increases by up to 475% and 227% in phantoms and patients, respectively, when adopting a PVC method that is based on a previously proposed methodology. 19 The demonstrated improvements in LBR should be seen in the light of recent studies by Derlin et al. and Fiz et al., who advocate the use of 18F-NaF PET imaging for the detection and characterization of vulnerable plaques. 7,10,11 Of note, the clinical implementation of the LP method differed from the pre-clinical version evaluated in previous work. 19 Here, the LP method was implemented only as a post-processing step which, together with the additional STIR reconstruction, results in a PVC image. In the pre-clinical implementation, the LP method was also implemented within the reconstruction process (PVC reconstruction). 19 In this work, we decided to evaluate a version of the LP method that would be easier to implement in clinical practice since it is based on a single post-processing step. This approach should make the algorithm more useful in clinical practice, where many retrospective studies do not have the projection data available and not all institutions have access to their own reconstruction algorithm. As we showed in our pre-clinical implementation, PVC reconstruction approaches can be implemented within the STIR reconstruction framework by calling the LP algorithm after each iteration.
The validation of our PVC method was performed by means of acquisition of a NEMA IQ phantom and a human-sized thorax phantom with three plaque-type lesions. Of note, the PSF matrix of the system must be known in order to ensure the best possible performance of the PVC method. Here, we made the assumption that the PSF can be described by a uniform Gaussian function, which is reasonable when the primary structures of interest (e.g., carotids) are located near the center of the PET transaxial field of view (FOV). The Full Width Half Maximum (FWHM) of the Gaussian blurring was obtained by fitting the LP tissue activities for the 37-mm sphere of the NEMA IQ phantom. As expected, the resulting FWHM values were different for the OSEM (FWHM * 8 mm) and PSF (FWHM *4 mm) reconstructions. The spatially invariant PSF approximation used in this study may not work properly for lesions located far away from the center of the transaxial FOV, given the spatial variance of the PSF. 32 In consequence, for lesions located far from the center of the FOV, we expect a reduced accuracy of the LP method. This limitation can be solved with a more accurate description of the PSF.
Both NEMA and thorax phantom experiments showed a significant improvement in quantification accuracy of the lesions when the PVC was applied in OSEM or PSF images (Figures 3, 4; Tables 2, 3). The LBR values for each lesion were closer to the reference in the PVC images. For lesions below 10 mm diameter, total recovery of the PVE was not achieved (Tables 2,  3), thus giving LBR values from the PVC images or from the LP tissue activities well below the experimental values. This is mainly due to two reasons: First, the PET-based segmentation will be significantly bigger than the real size of the lesion due to the spread of activity of these very small sources. Second, for these lesions, the Nyquist sampling condition (lesion sizes bigger than 3 voxels in each spatial direction, voxel size 2.03 mm) is not satisfied, and, therefore, a full recovery of the lesion activities using the LP method is not possible. 19 Nonetheless, large activity recovery values were obtained when applying the PVC method. These values could be further enhanced by reducing the voxel size in the PET image and by using a more accurate segmentation from a high-resolution anatomical image.
In general, the measured LP tissue activity ratios were more accurate than those obtained after the additional reconstruction with STIR. While the lesion contrast always appeared to be quantitatively improved in the PVC image (compared to OSEM or PSF alone), the best quantitative results were nevertheless generally obtained using the directly computed LP activity values with no further image reconstruction process.
The generally higher LBR values observed in the PSF reconstructions are more likely due to Gibbs artifacts (overshoots around the boundaries of hot lesions or organs), which are usually inherent to reconstruction algorithms with resolution recovery. 32,33 We observed that the proposed PVC method corrected for the overestimation of activities, giving LBR closer to the experimental values in the PSF ? PVC images and in the LP tissue activities obtained from the PSF images ( Table 2).
The performance of the LP method for different segmentation algorithms was evaluated with the NEMA IQ phantom acquisition ( Table 2). In general, slightly higher LBR values were observed when the PET threshold-based segmentation was used during PVC. This behavior was more pronounced for smaller spheres. This is due to the fact that the segmented volume of the spheres, using the 50% threshold adapted to background, was always smaller than the actual volume of the sphere (about 70-80% of the actual volume). This is in concordance with similar studies in the literature. 34 However, for the plaque-type lesions in the thorax phantom and the smallest non-calcified plaques in the patients, as their size is similar to the voxel size, the LP method using a PET-based segmentation is expected to under-perform the LP method using a CT-based segmentation, due to the spread of activity of very small hot sources.
In patient images, no significant correlations were observed between 18 F-NaF plaque uptake in the uncorrected images and CT-based calcifications. However, in the new reconstructed images with PVC, a significant correlation of 18 F-NaF uptake and calcification density of the atherosclerotic plaque was observed, due to the better quantitative evaluation of 18 F-NaF uptake in the calcified plaques (Table 6). Furthermore, the quantification and inverse correlation further improved when considering the tissue activity values obtained with the LP method ( Table 6). As expected, the effect of the PVC was more significant in small calcifications, thus resulting in higher DLBR for the smaller light plaque lesions ( Figure 6). The non-calcified plaques evaluated in this work were relatively large in size (Tables 4, 5). This result in relatively low DLBR values, compared with the ones obtained for light plaque lesions.
For the calcified plaques, a CT-based segmentation of the local VOI was used for the PVC. A limitation of this approach, as for the majority of PVC techniques, is that it can be affected by miss-registration errors. 35 We attempted to limit such errors by carefully verifying the spatial alignment of the PET and CT images prior to applying the PVC. In non-calcified lesions, a PET threshold-based segmentation of the local VOI was used for the PVC. This method works well provided the lesion is clearly visible in the non-corrected image and, therefore, can be segmented accurately. This was the case for the non-calcified lesions analyzed in this work. Table 6. Fitting parameters for the fits performed in Figure 6, and Pearson (R) and Spearman (S) correlation parameters for the data evaluated The main limitation of the method is that noncalcified lesions must be present and detectable on either PET or CT in order for the proposed PVC to be applicable. Another potential limitation of this study is the possibility of non-compensated patient motion during the acquisitions, which could also affect the coregistration of the PET and CT images, and the quantitative properties of the PET image. However, we consider involuntary motion effect to be of minor importance for this evaluation, considering the short PET acquisition time (2 minute per bed position).
Although pathological features in carotid or aortic arteries are likely representative of general atherosclerotic disease, it is known that most of the deaths from cardiovascular events result from complications of inflammatory plaques in the smaller coronary plaques. 36 The application of the PVC method for the evaluation of 18F-NaF uptake in vulnerable plaque in the coronaries, however, represents a more challenging problem, due to the presence of respiratory and cardiac motion, which have to be compensated before applying the PVC. 27 The evaluation of an image reconstruction framework, which incorporates motion-compensation (MoCo) and partial volume correction (PVC) for 18F-NaF PET imaging of vulnerable plaque in the coronary arteries, is currently work in progress, with first promising results in simulated data. 37

NEW KNOWLEGDE GAINED
We have demonstrated the need for partial volume correction when imaging small structures of interest, such as atherosclerotic plaques in the carotid and aortic arteries by means of PET. PVC can lead to increased LBR max values up to a factor of 3 in small plaque lesions.

CONCLUSION
The good performance of the LP-based PVC method was confirmed in both OSEM and PSF reconstructions of phantom data and clinical patients. LBR of plaque-type lesions increased by up to 475% (phantom) and 227% (patients) when PVC was applied. Results from this study following PVC further support reports of an inverse correlation of 18 F-NaF uptake in calcified plaques with plaque density.

Open Access
This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

APPENDIX 1: DESCRIPTION OF THE LP METHOD
Consider a reconstructed image with a detected lesion (for instance, an atherosclerotic plaque that shows increased uptake of the administered radiotracer). Our goal is to improve the quantification of the activity concentration in that lesion and removing spill-over activity from neighboring tissues. If we have a precise segmentation of organs and tissues, we can take advantage of this information and improve PET quantification based on the following: • It is reasonable to assume that voxels belonging to the same tissue have a more similar activity than voxels belonging to different tissues. In consequence, voxels from different tissues would be represented by different average activity values.
• We can also estimate the PET counts assigned to the inside of the VOI which actually originated from outside the VOI.
With these assumptions, the measured data, k i , can be modeled as where k i are the expected counts per sinogram bin, i, A j the activity for each segmented tissue j within the VOI, P ij is the resolution-blurred tissue shape function for tissue j and sinogram bin i, and g out;i represents the background counts coming from the region outside the VOI. The joint likelihood of measuring a given projection data set, n, is given by the product of the Poisson probability density function for each measured projection ray, i, over all rays traversing the VOI: where n i are the measured counts per LOR. Taking the