Prognostic Value of Early Fluorodeoxyglucose-Positron Emission Tomography Response Imaging and Peripheral Immunologic Biomarkers: Substudy of a Phase II Trial of Risk-Adaptive Chemoradiation for Unresectable Non-Small Cell Lung Cancer

Purpose

We sought to examine the prognostic value of fluorodeoxyglucose-positron emission tomography (PET) imaging during chemoradiation for unresectable non-small cell lung cancer for survival and hypothesized that tumor PET response is correlated with peripheral T-cell function.

Methods and Materials

Forty-five patients with American Joint Committee on Cancer version 7 stage IIB-IIIB non-small cell lung cancer enrolled in a phase II trial and received platinum-doublet chemotherapy concurrent with 6 weeks of radiation (NCT02773238). Fluorodeoxyglucose-PET was performed before treatment start and after 24 Gy of radiation (week 3). PET response status was prospectively defined by multifactorial radiologic interpretation. PET responders received 60 Gy in 30 fractions, while nonresponders received concomitant boosts to 74 Gy in 30 fractions. Peripheral blood was drawn synchronously with PET imaging, from which germline DNA sequencing, T-cell receptor sequencing, and plasma cytokine analysis were performed.

Results

Median follow-up was 18.8 months, 1-year overall survival (OS) 82%, 1-year progression-free survival 53%, and 1-year locoregional control 88%. Higher midtreatment PET total lesion glycolysis was detrimental to OS (1 year 87% vs 63%, P < .001), progression-free survival (1 year 60% vs 26%, P = .044), and locoregional control (1 year 94% vs 65%, P = .012), even after adjustment for clinical/treatment factors. Twenty-nine of 45 patients (64%) were classified as PET responders based on a priori definition. Higher tumor programmed death-ligand 1 expression was correlated with response on PET (P = .017). Higher T-cell receptor richness and clone distribution slope were associated with improved OS (P = .018-0.035); clone distribution slope was correlated with PET response (P = .031).

Conclusions

Midchemoradiation PET imaging is prognostic for survival; PET response may be linked to tumor and peripheral T-cell biomarkers.

Sensitivity analysis of FDG PET tumor voxel cluster radiomics and dosimetry for predicting mid-chemoradiation regional response of locally advanced lung cancer

We investigated the sensitivity of regional tumor response prediction to variability in voxel clustering techniques, imaging features, and machine learning algorithms in 25 patients with locally advanced non-small cell lung cancer (LA-NSCLC) enrolled on the FLARE-RT clinical trial. Metabolic tumor volumes (MTV) from pre-chemoradiation (PETpre) and mid-chemoradiation fluorodeoxyglucose-positron emission tomography (FDG PET) images (PETmid) were subdivided into K-means or hierarchical voxel clusters by standardized uptake values (SUV) and 3D-positions. MTV cluster separability was evaluated by CH index, and morphologic changes were captured by Dice similarity and centroid Euclidean distance. PETpre conventional features included SUVmean, MTV/MTV cluster size, and mean radiation dose. PETpre radiomics consisted of 41 intensity histogram and 3D texture features (PET Oncology Radiomics Test Suite) extracted from MTV or MTV clusters. Machine learning models (multiple linear regression, support vector regression, logistic regression, support vector machines) of conventional features or radiomic features were constructed to predict PETmid response. Leave-one-out-cross-validated root-mean-squared-error (RMSE) for continuous response regression (ΔSUVmean) and area-under-receiver-operating-characteristic-curve (AUC) for binary response classification were calculated. K-means MTV 2-clusters (MTVhi, MTVlo) achieved maximum CH index separability (Friedman p < 0.001). Between PETpre and PETmid, MTV cluster pairs overlapped (Dice 0.70-0.87) and migrated 0.6-1.1 cm. PETmid ΔSUVmean response prediction was superior in MTV and MTVlo (RMSE = 0.17-0.21) compared to MTVhi (RMSE = 0.42-0.52, Friedman p < 0.001). PETmid ΔSUVmean response class prediction performance trended higher in MTVlo (AUC = 0.83-0.88) compared to MTVhi (AUC = 0.44-0.58, Friedman p = 0.052). Models were more sensitive to MTV/MTV cluster regions (Friedman p = 0.026) than feature sets/algorithms (Wilcoxon signed-rank p = 0.36). Top-ranked radiomic features included GLZSM-LZHGE (large-zone-high-SUV), GTSDM-CP (cluster-prominence), GTSDM-CS (cluster-shade) and NGTDM-CNT (contrast). Top-ranked features were consistent between MTVhi and MTVlo cluster pairs but varied between MTVhi-MTVlo clusters, reflecting distinct regional radiomic phenotypes. Variability in tumor voxel cluster response prediction can inform robust radiomic target definition for risk-adaptive chemoradiation in patients with LA-NSCLC. FLARE-RT trial: NCT02773238.

Impact of contouring variability on oncological PET radiomics features in the lung

Radiomics features extracted from oncological PET images are currently under intense scrutiny within the context of risk stratification for a variety of cancers. However, the lack of robustness assessment poses problems for their application across institutions and for broader patient populations. The objective of the current study was to examine the extent to which radiomics parameters from oncological PET vary in response to manual contouring variability in lung cancer. Imaging data employed in the study consisted of 26 PET scans with lesions in the lung being created through the use of an anthropomorphic phantom in conjunction with Monte Carlo simulations. From each of the simulated lesions, 25 radiomics features related to the gray-level co-occurrence matrices (GLCOM), gray-level size zone matrices (GLSZM), and gray-level neighborhood difference matrices (GLNDM) were extracted from ground truth contour and from manual contours provided by 10 raters in regard to four intensity discretization schemes with number of gray levels of 32, 64, 128, and 256, respectively. The impact of interrater variability in tumor delineation upon the agreement between raters on radiomics features was examined via interclass correlation and leave-p-out assessment. Only weak and moderate correlations were found between segmentation accuracy as measured by the Dice coefficient and percent feature error from ground truth for the vast majority of the features being examined. GLNDM-based texture parameters emerged as the top performing category of radiomcs features in terms of robustness against contouring variability for discretization schemes engaging number of gray levels of 32, 64, and 128 while GLCOM-based parameters stood out for discretization scheme engaging 256 gray levels. How and to what extent interrater reliability of radiomics features vary in response to the number of raters were largely feature-dependent. It was concluded that impact of contouring variability on PET-based radiomics features is present to varying degrees and could be experienced as a barrier to convey PET-based radiomics research to clinical relevance.

Quantitative imaging: Erring patterns in manual delineation of PET-imaged lung lesions

BACKGROUND AND PURPOSE

Uncertainty and variability in manual contouring of PET-imaged tumor targets are well recognized; however, the error patterns associated with it were little known and rarely investigated. The present study is aimed to quantitatively assess the erring patterns inherent to manual delineation of PET-imaged lung lesions in a setting with complete ground truth.

MATERIALS AND METHODS

Images being used for assessment consisted of 26 synthetic PET datasets created by using the anthropomorphic Zubal phantom in conjunction with the Monte Carlo based SimSET computational package. Each dataset included one PET-positive lesion differing in shape, dimension, uptake heterogeneity, and anatomical location inside the lung. Target contours were provided by 10 raters and the contour accuracy was evaluated using 12 metrics from five categories - spatial overlap, pair counting, information theory, distance, and volume.

RESULTS

In terms of spatial overlap, manual contouring results intersect substantially with the ground truth whereas tend to oversegment the lesions. Shapes of the segmented tumor volumes are in general geometrically consistent with the ground truth but lack sensitivity in characterizing topographical details. No complete consensus could be achieved between manual contours and the ground truth for any of the given lesions being examined when assessing using pair counting- and informatics-based metrics thus indicating an intrinsic stochastic component of manual contouring. Evaluation based on metrics related to distance and volume demonstrated that it is at the borderline areas between the lesions and the normal tissues where the majority part of manual delineation errors occurred and the extent of volume being identified false positively as cancerous by the raters is appreciable.

CONCLUSION

Quantification of segmentation errors associated with expert manual contouring of PET positive lesion in the lung reveals general patterns in what otherwise might be thought of as randomness. Findings from the current study may allow for the formation of new hypotheses towards improving the accuracy and precision of manual delineation of PET positive tumor targets in the lung.

Quantitative radiomics: Validating image textural features for oncological PET in lung cancer

BACKGROUND AND PURPOSE

Radiomics textural features derived from PET imaging are of broad and current interest due to recent evidence of their prognostic value during cancer management. An inherent assumption is the link between these imaging features and the underlying tumoral phenotypic spatial heterogeneity. The purpose of this work was to validate this assumption for tumors within the lung through a comparison of image based textural features and the ground truth activity distribution from which the images were created. A second purpose was to assess the level at which PET imaging introduces spatial texture not present in the associated ground truth activity distribution.

MATERIALS AND METHODS

25 lung lesions were created using an anthropomorphic phantom. Ten of the lesions had a spherical shape with a uniform activity distribution. The remaining 15 had an irregular shape with a heterogeneous activity distribution. PET images were created for each lesion using Monte Carlo simulation. 79 textural features related to the gray-level intensity histograms, co-occurrence matrices, neighborhood difference matrices, run length, and size zone matrices were derived from both the simulated PET images and ground truth activity maps. A comparison was made between the two datasets using statistical analysis.

RESULTS

For homogenous lesions, features extracted from the PET images were largely irrelevant to the underlying uniform activity distribution. Additionally, the majority of these features assumed substantial values implying that an extensive amount of spatial texture had been introduced into the final imaging data. For heterogeneous lesions, complex trends were observed in the deviation between features extracted from PET images and those extracted from the ground truth activity maps. Moreover, the extent of both the deviation and the associated dynamic range was seen to be greatly feature-dependent.

CONCLUSION

The use of image based textural features as a surrogate for tumoral phenotypic spatial heterogeneity could not be clearly validated. The association between the two is complex and a significant amount of uncertainty exist due to the introduction of incidental texture during image acquisition and reconstruction.

The first MICCAI challenge on PET tumor segmentation

INTRODUCTION

Automatic functional volume segmentation in PET images is a challenge that has been addressed using a large array of methods. A major limitation for the field has been the lack of a benchmark dataset that would allow direct comparison of the results in the various publications. In the present work, we describe a comparison of recent methods on a large dataset following recommendations by the American Association of Physicists in Medicine (AAPM) task group (TG) 211, which was carried out within a MICCAI (Medical Image Computing and Computer Assisted Intervention) challenge.

MATERIALS AND METHODS

Organization and funding was provided by France Life Imaging (FLI). A dataset of 176 images combining simulated, phantom and clinical images was assembled. A website allowed the participants to register and download training data (n = 19). Challengers then submitted encapsulated pipelines on an online platform that autonomously ran the algorithms on the testing data (n = 157) and evaluated the results. The methods were ranked according to the arithmetic mean of sensitivity and positive predictive value.

RESULTS

Sixteen teams registered but only four provided manuscripts and pipeline(s) for a total of 10 methods. In addition, results using two thresholds and the Fuzzy Locally Adaptive Bayesian (FLAB) were generated. All competing methods except one performed with median accuracy above 0.8. The method with the highest score was the convolutional neural network-based segmentation, which significantly outperformed 9 out of 12 of the other methods, but not the improved K-Means, Gaussian Model Mixture and Fuzzy C-Means methods.

CONCLUSION

The most rigorous comparative study of PET segmentation algorithms to date was carried out using a dataset that is the largest used in such studies so far. The hierarchy amongst the methods in terms of accuracy did not depend strongly on the subset of datasets or the metrics (or combination of metrics). All the methods submitted by the challengers except one demonstrated good performance with median accuracy scores above 0.8.

Quantitative imaging: Correlating image features with the segmentation accuracy of PET based tumor contours in the lung

The purpose of this study was to investigate the correlation between image features extracted from PET images and the accuracy of manually drawn lesion contours in the lung. Such correlations are interesting in that they could potentially be used in predictive models to help guide physician contouring. In this work, 26 synthetic PET datasets were created using an anthropomorphic phantom and Monte Carlo simulation. Manual contours of simulated lesions were provided by 10 physicians. Contour accuracy was quantified using five commonly used similarity metrics which were then correlated with several features extracted from the images. Features were sub-divided into three groups using intensity, geometry, and texture as categorical descriptors. When averaged among the participants, the results showed relatively strong correlations with complexity and contrast_I (r≥0.65, p<0.001), and moderate correlations with several other image features (r≥0.5, p<0.01). The predictive nature of these correlations was improved through stepwise regression and the creation of multi-feature models. Imaging features were also correlated with the standard deviation of contouring error in order to investigate inter-observer variability. Several features were consistently identified as influential including integral of mean curvature and complexity. These relationships further the understanding as to what causes variation in the contouring of PET positive lesions.

A Dirichlet process mixture model for automatic (18)F-FDG PET image segmentation: Validation study on phantoms and on lung and esophageal lesions

PURPOSE

The aim of this study was to implement a Dirichlet process mixture (DPM) model for automatic tumor edge identification on (18)F-fluorodeoxyglucose positron emission tomography ((18)F-FDG PET) images by optimizing the parameters on which the algorithm depends, to validate it experimentally, and to test its robustness.

METHODS

The DPM model belongs to the class of the Bayesian nonparametric models and uses the Dirichlet process prior for flexible nonparametric mixture modeling, without any preliminary choice of the number of mixture components. The DPM algorithm implemented in the statistical software package R was used in this work. The contouring accuracy was evaluated on several image data sets: on an IEC phantom (spherical inserts with diameter in the range 10-37 mm) acquired by a Philips Gemini Big Bore PET-CT scanner, using 9 different target-to-background ratios (TBRs) from 2.5 to 70; on a digital phantom simulating spherical/uniform lesions and tumors, irregular in shape and activity; and on 20 clinical cases (10 lung and 10 esophageal cancer patients). The influence of the DPM parameters on contour generation was studied in two steps. In the first one, only the IEC spheres having diameters of 22 and 37 mm and a sphere of the digital phantom (41.6 mm diameter) were studied by varying the main parameters until the diameter of the spheres was obtained within 0.2% of the true value. In the second step, the results obtained for this training set were applied to the entire data set to determine DPM based volumes of all available lesions. These volumes were compared to those obtained by applying already known algorithms (Gaussian mixture model and gradient-based) and to true values, when available.

RESULTS

Only one parameter was found able to significantly influence segmentation accuracy (ANOVA test). This parameter was linearly connected to the uptake variance of the tested region of interest (ROI). In the first step of the study, a calibration curve was determined to automatically generate the optimal parameter from the variance of the ROI. This "calibration curve" was then applied to contour the whole data set. The accuracy (mean discrepancy between DPM model-based contours and reference contours) of volume estimation was below (1 ± 7)% on the whole data set (1 SD). The overlap between true and automatically segmented contours, measured by the Dice similarity coefficient, was 0.93 with a SD of 0.03.

CONCLUSIONS

The proposed DPM model was able to accurately reproduce known volumes of FDG concentration, with high overlap between segmented and true volumes. For all the analyzed inserts of the IEC phantom, the algorithm proved to be robust to variations in radius and in TBR. The main advantage of this algorithm was that no setting of DPM parameters was required in advance, since the proper setting of the only parameter that could significantly influence the segmentation results was automatically related to the uptake variance of the chosen ROI. Furthermore, the algorithm did not need any preliminary choice of the optimum number of classes to describe the ROIs within PET images and no assumption about the shape of the lesion and the uptake heterogeneity of the tracer was required.

(90)Y-PET/CT Imaging Quantification for Dosimetry in Peptide Receptor Radionuclide Therapy: Analysis and Corrections of the Impairing Factors

PURPOSE

We evaluated the possibility to assess (90)Y-PET/CT imaging quantification for dosimetry in (90)Y-peptide receptor radionuclide therapy.

METHODS

Tests were performed by Discovery 710 Elite (GE) PET/CT equipment. A body-phantom containing radioactive-coplanar-spheres was filled with (90)Y water solution to reproduce different signal-to-background-activity-ratios (S/N). We studied minimum detectable activity (MDA) concentration, contrast-to-noise ratio (CNR), and full-width-at-half-maximum (FWHM). Subsequently, three recovery coefficients (RC)-based correction approaches were evaluated: maximum-RC, resolution-RC, and isovolume-RC. The analysis of the volume segmentation thresholding method was also assessed to derive a relationship between the true volume of the targets and the threshold to be applied to the PET images. (90)Y-PET/CT imaging quantification was then achieved on some patients and related with preclinical tests. Moreover, the dosimetric evaluation was obtained on the target regions.

RESULTS

CNR value was greater than 5 if the MDA was greater than 0.2 MBq/mL with no background activity and 0.5-0.7 MBq/mL with S/N ranging from 3 to 6. FWHM was equal to 7 mm. An exponential fitting of isovolume RCs-based correction technique was adopted for activity quantification. Adaptive segmentation thresholding exponential curves were obtained and applied for target volume identification in three signal-to-background-activity-ratios. The imaging quantification study and dosimetric evaluations in clinical cases was feasible and the results were coherent with those obtained in preclinical tests.

CONCLUSIONS

(90)Y-PET/CT imaging quantification is possible both in phantoms and in patients. Absorbed dose evaluations in clinical applications are strongly related to targets activity concentration.

A gel tumour phantom for assessment of the accuracy of manual and automatic delineation of gross tumour volume from FDG-PET/CT

INTRODUCTION

Our primary aim was to make a phantom for PET that could mimic a highly irregular tumour and provide true tumour contours. The secondary aim was to use the phantom to assess the accuracy of different methods for delineation of tumour volume from the PET images.

MATERIAL AND METHODS

An empty mould was produced on the basis of a contrast enhanced computed tomography (CT) study of a patient with a squamous cell carcinoma in the head and neck region. The mould was filled with a homogeneous fast-settling gel that contained both (18)F for positron emission tomography (PET) and an iodine contrast agent. This phantom (mould and gel) was scanned on a PET/CT scanner. A series of reference tumour contours were obtained from the CT images in the PET/CT. Tumour delineation based on the PET images was achieved manually, by isoSUV thresholding, and by a recently developed three-dimensional (3D) Difference of Gaussians algorithm (DoG). Average distances between the PET-derived and reference contours were assessed by a 3D distance transform.

RESULTS

The manual, thresholding and DoG delineation methods resulted in volumes that were 146%, 86% and 100% of the reference volume, respectively, and average distance deviations from the reference surface were 1.57 mm, 1.48 mm and 1.0, mm, respectively.

DISCUSSION

Manual drawing as well as isoSUV determination of tumour contours in geometrically irregular tumours may be unreliable. The DoG method may contribute to more correct delineation of the tumour. Although the present phantom had a homogeneous distribution of activity, it may also provide useful knowledge in the case of inhomogeneous activity distributions.

CONCLUSION

The geometric irregular tumour phantom with its inherent reference contours was an important tool for testing of different delineation methods and for teaching delineation.

MRI Determination of Knee Effusion Volume: A Cadaveric Study

Background

There is currently limited literature on quantitative determination of knee effusion volume using magnetic resonance imaging (MRI).

Purpose

To describe a method of knee effusion volume determination using MRI generated models and to demonstrate accuracy of this technique.

Materials and methods

Using axial T2-weighted turbo spin echo and sagittal SPACE sequences, MRIs of three cadaver knees with multiple saline loads were obtained. Effusions models were created and effusion volumes were estimated using the Rhinoceros software. Estimated and known effusion volumes were compared using a bivariate correlation analysis.

Results

The SPACE sequence and T2WTSE estimates were highly correlated with the known volumes (R = 0.996 and 0.993 respectively, p < 0.001).

Conclusion

MRI-generated models of knee effusions provide accurate estimates of knee effusion volumes.

Clinical relevance

MRI determination of knee effusion volume may provide a useful clinical outcomes tool.

Wittstein J, Spritzer C, Garrett WE. MRI Determination of Knee Effusion Volume: A Cadaveric Study. The Duke Orthop J 2013;3(1):67-70.

Modeling and Simulation of 4D PET-CT and PET-MR Images

The driving force in the research and development of new hybrid PET-CT/MR imaging scanners is the production of images with optimum quality, accuracy, and resolution. However, the acquisition process is limited by several factors. Key issues are the respiratory and cardiac motion artifacts that occur during an imaging session. In this article the necessary tools for modeling and simulation of realistic high-resolution four-dimensional PET-CT and PET-MR imaging data are described. Beyond the need for four-dimensional simulations, accurate modeling of the acquisition process can be included within the reconstruction algorithms assisting in the improvement of image quality and accuracy of estimation of physiologic parameters from four-dimensional hybrid PET imaging.

The use of zeolites to generate PET phantoms for the validation of quantification strategies in oncology

PURPOSE

In recent years, segmentation algorithms and activity quantification methods have been proposed for oncological (18)F-fluorodeoxyglucose (FDG) PET. A full assessment of these algorithms, necessary for a clinical transfer, requires a validation on data sets provided with a reliable ground truth as to the imaged activity distribution, which must be as realistic as possible. The aim of this work is to propose a strategy to simulate lesions of uniform uptake and irregular shape in an anthropomorphic phantom, with the possibility to easily obtain a ground truth as to lesion activity and borders.

METHODS

Lesions were simulated with samples of clinoptilolite, a family of natural zeolites of irregular shape, able to absorb aqueous solutions of (18)F-FDG, available in a wide size range, and nontoxic. Zeolites were soaked in solutions of (18)F-FDG for increasing times up to 120 min and their absorptive properties were characterized as function of soaking duration, solution concentration, and zeolite dry weight. Saturated zeolites were wrapped in Parafilm, positioned inside an Alderson thorax-abdomen phantom and imaged with a PET-CT scanner. The ground truth for the activity distribution of each zeolite was obtained by segmenting high-resolution finely aligned CT images, on the basis of independently obtained volume measurements. The fine alignment between CT and PET was validated by comparing the CT-derived ground truth to a set of zeolites' PET threshold segmentations in terms of Dice index and volume error.

RESULTS

The soaking time necessary to achieve saturation increases with zeolite dry weight, with a maximum of about 90 min for the largest sample. At saturation, a linear dependence of the uptake normalized to the solution concentration on zeolite dry weight (R(2) = 0.988), as well as a uniform distribution of the activity over the entire zeolite volume from PET imaging were demonstrated. These findings indicate that the (18)F-FDG solution is able to saturate the zeolite pores and that the concentration does not influence the distribution uniformity of both solution and solute, at least at the trace concentrations used for zeolite activation. An additional proof of uniformity of zeolite saturation was obtained observing a correspondence between uptake and adsorbed volume of solution, corresponding to about 27.8% of zeolite volume. As to the ground truth for zeolites positioned inside the phantom, the segmentation of finely aligned CT images provided reliable borders, as demonstrated by a mean absolute volume error of 2.8% with respect to the PET threshold segmentation corresponding to the maximum Dice.

CONCLUSIONS

The proposed methodology allowed obtaining an experimental phantom data set that can be used as a feasible tool to test and validate quantification and segmentation algorithms for PET in oncology. The phantom is currently under consideration for being included in a benchmark designed by AAPM TG211, which will be available to the community to evaluate PET automatic segmentation methods.

A segmentation framework towards automatic generation of boost subvolumes for FDG-PET tumors: a digital phantom study

Potential benefits of administering nonuniform radiation dose to heterogeneous tumors imaged with FDG-PET have been widely demonstrated; whereas the number of discrete dose levels to be utilized and corresponding locations for prescription inside tumors vary significantly with current existing methods. In this paper, an automated and unsupervised segmentation framework constituted mainly by an image restoration mechanism based on variational decomposition and a voxel clustering scheme based on spectral clustering was presented towards partitioning FDG-PET imaged tumors into subvolumes characterized with the total intra-subvolume activity similarity and the total inter-subvolume activity dissimilarity being simultaneously maximized. Experiments to evaluate the proposed system were carried out with using FDG-PET data generated from a digital phantom that employed SimSET (Simulation System for Emission Tomography) to simulate PET acquisition of tumors. The obtained results show the feasibility of the proposed system in dividing FDG-PET imaged tumor volumes into subvolumes with intratumoral heterogeneity being properly characterized, irrespective of variation in tumor morphology as well as diversity in intratumoral heterogeneity pattern.

Delineation of FDG-PET tumors from heterogeneous background using spectral clustering

This paper explored the feasibility of using spectral clustering to segment FDG-PET tumor in the presence of heterogeneous background. Spectral clustering refers to a class of clustering methods which employ the eigenstructure of a similarity matrix to partition image voxels into disjoint clusters. The similarity between two voxels was measured with the intensity distance scaled by voxel-varying factors capturing local statistics and the number of clusters was inferred based on rotating the eigenvector matrix for the maximally sparse representation. Metrics used to evaluate the segmentation accuracy included: Dice coefficient, Jaccard coefficient, false positive dice, false negative dice, symmetric mean absolute surface distance, and absolute volumetric difference. Comparison of segmentation results between the presented method and the adaptive thresholding method on the simulated PET data shows the former attains an overall better detection accuracy. Applying the presented method on patient data gave segmentation results in fairly good agreement with physician manual annotations. These results indicate that the presented method have the potential to accurately delineate complex shaped FDG-PET tumors containing inhomogeneous activities in the presence of heterogeneous background.

What is the best way to contour lung tumors on PET scans? Multiobserver validation of a gradient-based method using a NSCLC digital PET phantom

PURPOSE

To evaluate the accuracy and consistency of a gradient-based positron emission tomography (PET) segmentation method, GRADIENT, compared with manual (MANUAL) and constant threshold (THRESHOLD) methods.

METHODS AND MATERIALS

Contouring accuracy was evaluated with sphere phantoms and clinically realistic Monte Carlo PET phantoms of the thorax. The sphere phantoms were 10-37 mm in diameter and were acquired at five institutions emulating clinical conditions. One institution also acquired a sphere phantom with multiple source-to-background ratios of 2:1, 5:1, 10:1, 20:1, and 70:1. One observer segmented (contoured) each sphere with GRADIENT and THRESHOLD from 25% to 50% at 5% increments. Subsequently, seven physicians segmented 31 lesions (7-264 mL) from 25 digital thorax phantoms using GRADIENT, THRESHOLD, and MANUAL.

RESULTS

For spheres <20 mm in diameter, GRADIENT was the most accurate with a mean absolute % error in diameter of 8.15% (10.2% SD) compared with 49.2% (51.1% SD) for 45% THRESHOLD (p < 0.005). For larger spheres, the methods were statistically equivalent. For varying source-to-background ratios, GRADIENT was the most accurate for spheres >20 mm (p < 0.065) and <20 mm (p < 0.015). For digital thorax phantoms, GRADIENT was the most accurate (p < 0.01), with a mean absolute % error in volume of 10.99% (11.9% SD), followed by 25% THRESHOLD at 17.5% (29.4% SD), and MANUAL at 19.5% (17.2% SD). GRADIENT had the least systematic bias, with a mean % error in volume of -0.05% (16.2% SD) compared with 25% THRESHOLD at -2.1% (34.2% SD) and MANUAL at -16.3% (20.2% SD; p value <0.01). Interobserver variability was reduced using GRADIENT compared with both 25% THRESHOLD and MANUAL (p value <0.01, Levene's test).

CONCLUSION

GRADIENT was the most accurate and consistent technique for target volume contouring. GRADIENT was also the most robust for varying imaging conditions. GRADIENT has the potential to play an important role for tumor delineation in radiation therapy planning and response assessment.

Artificial Neural Network-Based System for PET Volume Segmentation

Tumour detection, classification, and quantification in positron emission tomography (PET) imaging at early stage of disease are important issues for clinical diagnosis, assessment of response to treatment, and radiotherapy planning. Many techniques have been proposed for segmenting medical imaging data; however, some of the approaches have poor performance, large inaccuracy, and require substantial computation time for analysing large medical volumes. Artificial intelligence (AI) approaches can provide improved accuracy and save decent amount of time. Artificial neural networks (ANNs), as one of the best AI techniques, have the capability to classify and quantify precisely lesions and model the clinical evaluation for a specific problem. This paper presents a novel application of ANNs in the wavelet domain for PET volume segmentation. ANN performance evaluation using different training algorithms in both spatial and wavelet domains with a different number of neurons in the hidden layer is also presented. The best number of neurons in the hidden layer is determined according to the experimental results, which is also stated Levenberg-Marquardt backpropagation training algorithm as the best training approach for the proposed application. The proposed intelligent system results are compared with those obtained using conventional techniques including thresholding and clustering based approaches. Experimental and Monte Carlo simulated PET phantom data sets and clinical PET volumes of nonsmall cell lung cancer patients were utilised to validate the proposed algorithm which has demonstrated promising results.

A novel fuzzy C-means algorithm for unsupervised heterogeneous tumor quantification in PET

PURPOSE

Accurate and robust image segmentation was identified as one of the most challenging issues facing PET quantification in oncological imaging. This difficulty is compounded by the low spatial resolution and high noise characteristics of PET images. The fuzzy C-means (FCM) clustering algorithm was largely used in various medical image segmentation approaches. However, the algorithm is sensitive to both noise and intensity heterogeneity since it does not take into account spatial contextual information.

METHODS

To overcome this limitation, a new fuzzy segmentation technique adapted to typical noisy and low resolution oncological PET data is proposed. PET images smoothed using a nonlinear anisotropic diffusion filter are added as a second input to the proposed FCM algorithm to incorporate spatial information (FCM-S). In addition, a methodology was developed to integrate the a trous wavelet transform in the standard FCM algorithm (FCM-SW) to allow handling of heterogeneous lesions' uptake. The algorithm was applied to the simulated data of the NCAT phantom, incorporating heterogeneous lesions in the lung and clinical PET/CT images of 21 patients presenting with histologically proven nonsmall-cell lung cancer (NSCLC) and 7 patients presenting with laryngeal squamous cell carcinoma (LSCC) to assess its performance for segmenting tumors with arbitrary size, shape, and tracer uptake. For NSCLC patients, the maximal tumor diameters measured from the macroscopic examination of the surgical specimen served as the ground truth for comparison with the maximum diameter estimated by the segmentation technique, whereas for LSCC patients, the 3D macroscopic tumor volume was considered as the ground truth for comparison with the corresponding PET-based volume. The proposed algorithm was also compared to the classical FCM segmentation technique.

RESULTS

There is a good correlation (R2 = 0.942) between the actual maximal diameter of primary NSCLC tumors estimated using the proposed PET segmentation procedure and those measured from the macroscopic examination, and the regression line agreed well with the line of identity (slope = 1.08) for the group analysis of the clinical data. The standard FCM algorithm seems to underestimate actual maximal diameters of the clinical data, resulting in a mean error of -4.6 mm (relative error of -10.8 +/- 23.1%) for all data sets. The mean error of maximal diameter estimation was reduced to 0.1 mm (0.9 +/- 14.4%) using the proposed FCM-SW algorithm. Likewise, the mean relative error on the estimated volume for LSCC patients was reduced from 21.7 +/- 22.0% for FCM to 8.6 +/- 28.3% using the proposed FCM-SW technique.

CONCLUSIONS

A novel unsupervised PET image segmentation technique that allows the quantification of lesions in the presence of heterogeneity of tracer uptake was developed and evaluated. The technique is being further refined and assessed in clinical setting to delineate treatment volumes for the purpose of PET-guided radiation therapy treatment planning but could find other applications in clinical oncology such as the assessment of response to treatment.

An automatic method for PET target segmentation using a lookup table based on volume and concentration ratio

Accurate evaluation of functionally significant target volumes in combination with anatomic imaging is of primary importance for effective radiation therapy treatment planning. In this study, a method for rapid and accurate PET image segmentation and volumetrics based on phantom measurements and independent of scanner calibration was developed. A series of spheres ranging in volume from 0.5 mL to 95 mL were imaged in an anthropomorphic phantom of human thorax using two commercial PET and CT/PET scanners. The target to background radioactivity concentration ratio ranged from 3:1 to 12:1 in 11 separate phantom scanning experiments. The results confirmed that optimal segmentation thresholding depends on target volume and radioactivity concentration ratio. This information can be derived from a generalized pre-determined "lookup table" of volume and contrast dependent threshold values instead of using fitted curves derived from machine specific information. A three-step method based on the PET image intensity information alone was used to delineate volumes of interest. First, a mean intensity segmentation method was used to generate an initial estimate of target volume, and the radioactivity concentration ratio was computed by a family of recovery coefficient curves to compensate for the partial volume effect. Next, the appropriate threshold value was obtained from a phantom-generated threshold lookup table. Lastly, a threshold level set method was performed on the threshold value to further refine the target contour by reducing the limitation of global thresholding. The segmentation results were consistent for spheres greater than 2.5 mL which yielded volume average uncertainty of 11.2% in phantom studies. The results of segmented volumes were comparable to those determined by contrast-oriented method and iterative threshold method (ITM). In addition, the new volume segmentation method was applied clinically to ten patients undergoing PET/CT volume analysis for radiation therapy treatment planning of solitary lung metastases. For these patients, the average PET segmented volumes were within 8.0% of the CT volumes and were highly dependent on the extension of functionally inactive tumor volume. In summary, the current method does not require fitted threshold curves or a priori knowledge of the CT/MRI target volume. This threshold method can be universally applied to radiation therapy treatment planning with comparable accuracy, and may be useful in the rapid identification and assessment of plans containing multiple targets.

PET-guided delineation of radiation therapy treatment volumes: a survey of image segmentation techniques

Historically, anatomical CT and MR images were used to delineate the gross tumour volumes (GTVs) for radiotherapy treatment planning. The capabilities offered by modern radiation therapy units and the widespread availability of combined PET/CT scanners stimulated the development of biological PET imaging-guided radiation therapy treatment planning with the aim to produce highly conformal radiation dose distribution to the tumour. One of the most difficult issues facing PET-based treatment planning is the accurate delineation of target regions from typical blurred and noisy functional images. The major problems encountered are image segmentation and imperfect system response function. Image segmentation is defined as the process of classifying the voxels of an image into a set of distinct classes. The difficulty in PET image segmentation is compounded by the low spatial resolution and high noise characteristics of PET images. Despite the difficulties and known limitations, several image segmentation approaches have been proposed and used in the clinical setting including thresholding, edge detection, region growing, clustering, stochastic models, deformable models, classifiers and several other approaches. A detailed description of the various approaches proposed in the literature is reviewed. Moreover, we also briefly discuss some important considerations and limitations of the widely used techniques to guide practitioners in the field of radiation oncology. The strategies followed for validation and comparative assessment of various PET segmentation approaches are described. Future opportunities and the current challenges facing the adoption of PET-guided delineation of target volumes and its role in basic and clinical research are also addressed.

Molecular PET/CT imaging-guided radiation therapy treatment planning

The role of positron emission tomography (PET) during the past decade has evolved rapidly from that of a pure research tool to a methodology of enormous clinical potential. (18)F-fluorodeoxyglucose (FDG)-PET is currently the most widely used probe in the diagnosis, staging, assessment of tumor response to treatment, and radiation therapy planning because metabolic changes generally precede the more conventionally measured parameter of change in tumor size. Data accumulated rapidly during the last decade, thus validating the efficacy of FDG imaging and many other tracers in a wide variety of malignant tumors with sensitivities and specificities often in the high 90 percentile range. As a result, PET/computed tomography (CT) had a significant impact on the management of patients because it obviated the need for further evaluation, guided further diagnostic procedures, and assisted in planning therapy for a considerable number of patients. On the other hand, the progress in radiation therapy technology has been enormous during the last two decades, now offering the possibility to plan highly conformal radiation dose distributions through the use of sophisticated beam targeting techniques such as intensity-modulated radiation therapy (IMRT) using tomotherapy, volumetric modulated arc therapy, and many other promising technologies for sculpted three-dimensional (3D) dose distribution. The foundation of molecular imaging-guided radiation therapy lies in the use of advanced imaging technology for improved definition of tumor target volumes, thus relating the absorbed dose information to image-based patient representations. This review documents technological advancements in the field concentrating on the conceptual role of molecular PET/CT imaging in radiation therapy treatment planning and related image processing issues with special emphasis on segmentation of medical images for the purpose of defining target volumes. There is still much more work to be done and many of the techniques reviewed are themselves not yet widely implemented in clinical settings.

Does running cause metatarsophalangeal joint effusions? A comparison of synovial fluid volumes on MRI in athletes before and after running

OBJECTIVE

The metatarsophalangeal joints (MTPJ) are the only joints that bear weight directly through synovium. The purpose of this study was to determine whether there is an association between synovial stresses during running and increases in volume of joint fluid.

MATERIALS AND METHODS

This was a prospective case controlled study (nine healthy athlete volunteers acting as own controls). High-resolution coronal 3D T2W magnetic resonance imaging of the MTPJs were obtained following 24 h rest and after a 30-min run. The volume of joint fluid in each MTPJ (n = 90) was measured by two independent observers using an automated propagating segmentation tool.

RESULTS

The median volume of synovial fluid in the MTPJs at rest was 0.018 ml (inter-quartile range (IQ) range 0.005-0.04) and after running 0.019 ml (IQ range 0.005-0.04, p = 0.34, 99% confidence interval (CI), 0.330.35). The volume of fluid in the MTPJs of the great toes was substantially larger than other toes (0.152 ml at rest, 0.154 ml after exercise, p = 0.903). Median volumes decrease from second to fifth MTPJs (0.032-0.007 ml at rest and 0.035-0.004 ml after exercise). Subset analysis for each toe revealed no significant difference in volumes before and after running (p = 0.39 to p = 0.9). The inter-rater reliability for observer measurements was good with an intra-class correlation of 0.70 (95% CI, 0.60 to 0.78).

CONCLUSION

It appears to be normal to find synovial fluid, particularly in the MTPJs of the great toes, of athletes at rest and after running. There does not appear to be an association between moderate distance running and an increase in the volume of synovial fluid.