We proposed a framework to detect and quantify local tumor morphological changes due to chemo-radiotherapy (CRT) using a Jacobian map and to extract quantitative radiomic features from the Jacobian map to predict the pathologic tumor response in locally advanced esophageal cancer patients. In 20 patients who underwent CRT, a multi-resolution BSpline deformable registration was performed to register the follow-up (post-CRT) CT to the baseline CT image. The Jacobian map ( J ) was computed as the determinant of the gradient of the deformation vector field. The Jacobian map measured the ratio of local tumor volume change where J < 1 indicated tumor shrinkage and J > 1 denoted expansion. The tumor was manually delineated and corresponding anatomical landmarks were generated on the baseline and follow-up images. Intensity, texture and geometry features were then extracted from the Jacobian map of the tumor to quantify tumor morphological changes. The importance of each Jacobian feature in predicting pathologic tumor response was evaluated by both univariate and multivariate analysis. We constructed a multivariate prediction model by using a support vector machine (SVM) classifier coupled with a least absolute shrinkage and selection operator (LASSO) for feature selection. The SVM-LASSO model was evaluated using ten-times repeated 10-fold cross-validation (10 × 10-fold CV). After registration, the average target registration error was 4.30 ± 1.09 mm (LR:1.63 mm AP:1.59 mm SI:3.05 mm) indicating registration error was within two voxels and close to 4 mm slice thickness. Visually, the Jacobian map showed smoothly-varying local shrinkage and expansion regions in a tumor. Quantitatively, the average median Jacobian was 0.80 ± 0.10 and 1.05 ± 0.15 for responder and non-responder tumors, respectively. These indicated that on average responder tumors had 20% median volume shrinkage while non-responder tumors had 5% median volume expansion. In univariate analysis, the minimum Jacobian ( p = 0.009, AUC = 0.98) and median Jacobian ( p = 0.004, AUC = 0.95) were the most significant predictors. The SVM-LASSO model achieved the highest accuracy when these two features were selected (sensitivity = 94.4%, specificity = 91.8%, AUC = 0.94). Novel features extracted from the Jacobian map quantified local tumor morphological changes using only baseline tumor contour without post-treatment tumor segmentation. The SVM-LASSO model using the median Jacobian and minimum Jacobian achieved high accuracy in predicting pathologic tumor response. The Jacobian map showed great potential for longitudinal evaluation of tumor response.
Accurate tumor segmentation in PET is crucial in many oncology applications. We developed an adaptive region-growing (ARG) algorithm with a maximum curvature strategy (ARG_MC) for tumor segmentation in PET. The ARG_MC repeatedly applied a confidence connected region-growing algorithm with increasing relaxing factor f . The optimal relaxing factor (ORF) was then determined at the transition point on the f -volume curve, where the volume just grew from the tumor into the surrounding normal tissues. The ARG_MC along with five widely used algorithms were tested on a phantom with 6 spheres at different signal to background ratios and on two clinic datasets including 20 patients with esophageal cancer and 11 patients with non-Hodgkin lymphoma (NHL). The ARG_MC did not require any phantom calibration or any a priori knowledge of the tumor or PET scanner. The identified ORF varied with tumor types (mean ORF = 9.61, 3.78 and 2.55 respectively for the phantom, esophageal cancer, and NHL datasets), and varied from one tumor to another. For the phantom, the ARG_MC ranked the second in segmentation accuracy with an average Dice similarity index (DSI) of 0.86, only slightly worse than Daisne’s adaptive thresholding method (DSI = 0.87), which required phantom calibration. For both the esophageal cancer dataset and the NHL dataset, the ARG_MC had the highest accuracy with an average DSI of 0.87 and 0.84, respectively. The ARG_MC was robust to parameter settings and region of interest selection, and it did not depend on scanners, imaging protocols, or tumor types. Furthermore, the ARG_MC made no assumption about the tumor size or tumor uptake distribution, making it suitable for segmenting tumors with heterogeneous FDG uptake. In conclusion, the ARG_MC was accurate, robust and easy to use, it provides a highly potential tool for PET tumor segmentation in clinic.
22 August 2022 in Journal of Applied Clinical Medical Physics
This study compared the reproducibility of chestwall and heart position using surface‐guided versus RPM (real‐time position management)‐guided deep inspiration breath hold (DIBH) radiotherapy for left sided breast cancer. Forty DIBH patients under either surface‐guided radiotherapy (SGRT) or RPM guidance were studied. For patients treated with tangential fields, reproducibility was measured as the displacements in central lung distance (CLD) and heart shadow to field edge distance (HFD) between pretreatment MV (megavoltage) images and planning DRRs (digitally reconstructed radiographs). For patients treated with volumetric modulated arc therapy (VMAT), sternum to isocenter (ISO) distance (StID), spine to rib edge distance (SpRD), and heart shadow to central axis (CAX) distance (HCD) between pretreatment kV images and planning DRRs were measured. These displacements were compared between SGRT and RPM‐guided DIBH. In tangential patients, the mean absolute displacements of SGRT versus RPM guidance were 0.19 versus 0.23 cm in CLD, and 0.33 versus 0.62 cm in HFD. With respect to planning DRR, heart appeared closer to the field edge by 0.04 cm with surface imaging versus 0.62 cm with RPM. In VMAT patients, the displacements of surface imaging versus RPM guidance were 0.21 versus 0.15 cm in StID, 0.24 versus 0.19 cm in SpRD, and 0.72 versus 0.41 cm in HCD. Heart appeared 0.41 cm further away from CAX with surface imaging, whereas 0.10 cm closer to field CAX with RPM. None of the differences between surface imaging and RPM guidance was statistically significant. In conclusion, the displacements of chestwall were small and were comparable with SGRT‐ or RPM‐guided DIBH. The position deviations of heart were larger than those of chestwall with SGRT or RPM. Although none of the differences between SGRT and RPM guidance were statistically significant, there was a trend that the position deviations of heart were smaller and more favorable with SGRT than with RPM guidance in tangential patients.
1 January 2023 in Computational and Structural Biotechnology Journal
Lung adenocarcinoma (ADC) is the most common non-small cell lung cancer. Surgical resection is the primary treatment for early-stage lung ADC while lung-sparing surgery is an alternative for non-aggressive cases. Identifying histopathologic subtypes before surgery helps determine the optimal surgical approach. Predominantly solid or micropapillary (MIP) subtypes are aggressive and associated with a higher likelihood of recurrence and metastasis and lower survival rates. This study aims to non-invasively identify these aggressive subtypes using preoperative 18F-FDG PET/CT and diagnostic CT radiomics analysis. We retrospectively studied 119 patients with stage I lung ADC and tumors ≤ 2 cm, where 23 had aggressive subtypes (18 solid and 5 MIPs). Out of 214 radiomic features from the PET/CT and CT scans and 14 clinical parameters, 78 significant features (3 CT and 75 PET features) were identified through univariate analysis and hierarchical clustering with minimized feature collinearity. A combination of Support Vector Machine classifier and Least Absolute Shrinkage and Selection Operator built predictive models. Ten iterations of 10-fold cross-validation (10 ×10-fold CV) evaluated the model. A pair of texture feature (PET GLCM Correlation) and shape feature (CT Sphericity) emerged as the best predictor. The radiomics model significantly outperformed the conventional predictor SUVmax (accuracy: 83.5% vs. 74.7%, p = 9e-9) and identified aggressive subtypes by evaluating FDG uptake in the tumor and tumor shape. It also demonstrated a high negative predictive value of 95.6% compared to SUVmax (88.2%, p = 2e-10). The proposed radiomics approach could reduce unnecessary extensive surgeries for non-aggressive subtype patients, improving surgical decision-making for early-stage lung ADC patients.
16 April 2019 in IEEE Transactions on Radiation and Plasma Medical Sciences
Positron emission tomography (PET) and computed tomography (CT) are widely used imaging modalities in radiation oncology. PET imaging has a high contrast but often leads to blurry tumor edges due to its limited spatial resolution, while CT imaging has a high resolution but a low contrast between a tumor and its surrounding normal soft tissues. Tumor segmentation from either a single PET or CT image is difficult. It is known that co-segmentation methods utilizing the complementary information between PET and CT can improve the segmentation accuracy. This complementary information can be either consistent or inconsistent in the image level. How to correctly localize tumor edges with the inconsistent information is one major challenge for co-segmentation. This paper proposes a novel variational method for tumor co-segmentation in PET/CT, with a fusion strategy specifically designed to handle the information inconsistency between PET and CT in an adaptive way-the method can automatically decide which modality should be more trustful when PET and CT disagree to each other for localizing the tumor boundary. A PET restoration process is further integrated into the co-segmentation process, which further eliminates the uncertainty for tumor segmentation introduced by tumor edge blur in PET. The performance of the proposed method is validated on 50 nonsmall cell lung cancer patients. It achieves an average dice similarity indexes (DSIs) of 0.85 +/- 0.06, volume error (VE) of 0.09 +/- 0.08, and classification error (CE) of 0.31 +/- 0.13. The proposed method has a high accuracy for PET/CT tumor co-segmentation and PET restoration.