Improved accuracy of quantitative parameter estimates in dynamic contrast-enhanced CT study with low temporal resolution.

PURPOSE: A previously proposed method to reduce radiation dose to patient in dynamic contrast-enhanced (DCE) CT is enhanced by principal component analysis (PCA) filtering which improves the signal-to-noise ratio (SNR) of time-concentration curves in the DCE-CT study. The efficacy of the combined method to maintain the accuracy of kinetic parameter estimates at low temporal resolution is investigated with pixel-by-pixel kinetic analysis of DCE-CT data. METHODS: The method is based on DCE-CT scanning performed with low temporal resolution to reduce the radiation dose to the patient. The arterial input function (AIF) with high temporal resolution can be generated with a coarsely sampled AIF through a previously published method of AIF estimation. To increase the SNR of time-concentration curves (tissue curves), first, a region-of-interest is segmented into squares composed of 3 × 3 pixels in size. Subsequently, the PCA filtering combined with a fraction of residual information criterion is applied to all the segmented squares for further improvement of their SNRs. The proposed method was applied to each DCE-CT data set of a cohort of 14 patients at varying levels of down-sampling. The kinetic analyses using the modified Tofts' model and singular value decomposition method, then, were carried out for each of the down-sampling schemes between the intervals from 2 to 15 s. The results were compared with analyses done with the measured data in high temporal resolution (i.e., original scanning frequency) as the reference. RESULTS: The patients' AIFs were estimated to high accuracy based on the 11 orthonormal bases of arterial impulse responses established in the previous paper. In addition, noise in the images was effectively reduced by using five principal components of the tissue curves for filtering. Kinetic analyses using the proposed method showed superior results compared to those with down-sampling alone; they were able to maintain the accuracy in the quantitative histogram parameters of volume transfer constant [standard deviation (SD), 98th percentile, and range], rate constant (SD), blood volume fraction (mean, SD, 98th percentile, and range), and blood flow (mean, SD, median, 98th percentile, and range) for sampling intervals between 10 and 15 s. CONCLUSIONS: The proposed method of PCA filtering combined with the AIF estimation technique allows low frequency scanning for DCE-CT study to reduce patient radiation dose. The results indicate that the method is useful in pixel-by-pixel kinetic analysis of DCE-CT data for patients with cervical cancer.

Quality assurance of asymmetric jaw alignment using 2D diode array.

PURPOSE: A method using a 2D diode array is proposed to measure the junction gap (or overlap) and dose with high precision for routine quality assurance of the asymmetric jaw alignment. METHODS: The central axis (CAX) of the radiation field was determined with a 15 × 15 cm(2) photon field at four cardinal collimator angles so that the junction gap (or overlap) can be measured with respect to the CAX. Two abutting fields having a field size of 15 cm (length along the axis parallel to the junction) × 7.5 cm (width along the axis perpendicular to the junction) were used to irradiate the 2D diode array (MapCHECK2) with 100 MU delivered at the photon energy of 6 MV. The collimator was slightly rotated at 15° with respect to the beam central axis to increase the number of diodes effective on the measurement of junction gap. The junction gap and dose measured in high spatial resolution were compared to the conventional methods using an electronic portal imaging device (EPID) and radiochromic film, respectively. In addition, the reproducibility and sensitivity of the proposed method to the measurements of junction gap and dose were investigated. RESULTS: The junction gap (or overlap) and dose measured by MapCHECK2 agreed well to those measured by the conventional methods of EPID and film (the differences ranged from -0.01 to 0 cm and from -1.34% to 0.6% for the gap and dose, respectively). No variation in the repeat measurements of the junction gap was found whereas the measurements of junction dose were found to vary in quite a small range over the days of measurement (0.21%-0.35%). While the sensitivity of the measured junction gap to the actual junction gap applied was the ideal value of 1 cm∕cm as expected, the sensitivity of the junction dose to the actual junction gap increased as the junction gap (or overlap) decreased (maximum sensitivity: 201.7%∕cm). CONCLUSIONS: The initial results suggest that the method is applicable for a comprehensive quality assurance of the asymmetric jaw alignment.

Volumetric-modulated arc therapy for the treatment of a large planning target volume in thoracic esophageal cancer.

Recently, volumetric-modulated arc therapy (VMAT) has demonstrated the ability to deliver radiation dose precisely and accurately with a shorter delivery time compared to conventional intensity-modulated fixed-field treatment (IMRT). We applied the hypothesis of VMAT technique for the treatment of thoracic esophageal carcinoma to determine superior or equivalent conformal dose coverage for a large thoracic esophageal planning target volume (PTV) with superior or equivalent sparing of organs-at-risk (OARs) doses, and reduce delivery time and monitor units (MUs), in comparison with conventional fixed-field IMRT plans. We also analyzed and compared some other important metrics of treatment planning and treatment delivery for both IMRT and VMAT techniques. These metrics include: 1) the integral dose and the volume receiving intermediate dose levels between IMRT and VMATI plans; 2) the use of 4D CT to determine the internal motion margin; and 3) evaluating the dosimetry of every plan through patient-specific QA. These factors may impact the overall treatment plan quality and outcomes from the individual planning technique used. In this study, we also examined the significance of using two arcs vs. a single-arc VMAT technique for PTV coverage, OARs doses, monitor units and delivery time. Thirteen patients, stage T2-T3 N0-N1 (TNM AJCC 7th edn.), PTV volume median 395 cc (range 281-601 cc), median age 69 years (range 53 to 85), were treated from July 2010 to June 2011 with a four-field (n = 4) or five-field (n = 9) step-and-shoot IMRT technique using a 6 MV beam to a prescribed dose of 50 Gy in 20 to 25 F. These patients were retrospectively replanned using single arc (VMATI, 91 control points) and two arcs (VMATII, 182 control points). All treatment plans of the 13 study cases were evaluated using various dose-volume metrics. These included PTV D99, PTV D95, PTV V9547.5Gy(95%), PTV mean dose, Dmax, PTV dose conformity (Van't Riet conformation number (CN)), mean lung dose, lung V20 and V5, liver V30, and Dmax to the spinal canal prv3mm. Also examined were the total plan monitor units (MUs) and the beam delivery time. Equivalent target coverage was observed with both VMAT single and two-arc plans. The comparison of VMATI with fixed-field IMRT demonstrated equivalent target coverage; statistically no significant difference were found in PTV D99 (p = 0.47), PTV mean (p = 0.12), PTV D95 and PTV V9547.5Gy (95%) (p = 0.38). However, Dmax in VMATI plans was significantly lower compared to IMRT (p = 0.02). The Van't Riet dose conformation number (CN) was also statistically in favor of VMATI plans (p = 0.04). VMATI achieved lower lung V20 (p = 0.05), whereas lung V5 (p = 0.35) and mean lung dose (p = 0.62) were not significantly different. The other OARs, including spinal canal, liver, heart, and kidneys showed no statistically significant differences between the two techniques. Treatment time delivery for VMATI plans was reduced by up to 55% (p = 5.8E-10) and MUs reduced by up to 16% (p = 0.001). Integral dose was not statistically different between the two planning techniques (p = 0.99). There were no statistically significant differences found in dose distribution of the two VMAT techniques (VMATI vs. VMATII) Dose statistics for both VMAT techniques were: PTV D99 (p = 0.76), PTV D95 (p = 0.95), mean PTV dose (p = 0.78), conformation number (CN) (p = 0.26), and MUs (p = 0.1). However, the treatment delivery time for VMATII increased significantly by two-fold (p = 3.0E-11) compared to VMATI. VMAT-based treatment planning is safe and deliverable for patients with thoracic esophageal cancer with similar planning goals, when compared to standard IMRT. The key benefit for VMATI was the reduction in treatment delivery time and MUs, and improvement in dose conformality. In our study, we found no significant difference in VMATII over single-arc VMATI for PTV coverage or OARs doses. However, we observed significant increase in delivery time for VMATII compared to VMATI.

A comparison of dynamic contrast-enhanced CT and MR imaging-derived measurements in patients with cervical cancer.

This work is to compare the kinetic parameters derived from the DCE-CT and -MR data of a group of 37 patients with cervical cancer. The modified Tofts model and the reference tissue method were applied to estimate kinetic parameters. In the MR kinetic analyses using the modified Tofts model for each patient data set, both the arterial input function (AIF) measured from DCE-MR images and a population-averaged AIF from the literature were applied to the analyses, while the measured AIF was used for the CT kinetic analysis. The kinetic parameters obtained from both modalities were compared. Significant moderate correlations were found in modified Tofts parameters [volume transfer constant(K(trans) ) and rate constant (k(ep) )] between CT and MR analysis for MR with the measured AIFs (R = 0·45, P<0·01 and R = 0·40, P<0·01 in high-K(trans) region; R = 0·38, P<0·01 and R = 0·80, P<0·01 in low-K(trans) region) as well as with the population-averaged AIF (R = 0·59, P<0·01 and R = 0·62, P<0·01 in high-K(trans) region; R = 0·50, P<0·01 and R = 0·63, P<0·01 in low-K(trans) region), respectively. In addition, from the Bland-Altman plot analysis, it was found that the systematic biases (the mean difference) between the modalities were drastically reduced in magnitude by adopting the population-averaged AIF for the MR analysis instead of the measured ones (from 51·5% to 18·9% for K(trans) and from 21·7% to 4·1% for k(ep) in high-K(trans) region; from 73·0% to 29·4% for K(trans) and from 63·4% to 24·5% for k(ep) in low-K(trans) region). The preliminary results showed the feasibility in the interchangeable use of the two imaging modalities in assessing cervical cancers.

Fast transit portal dosimetry using density-scaled layer modeling of aSi-based electronic portal imaging device and Monte Carlo method.

PURPOSE: Fast and accurate transit portal dosimetry was investigated by developing a density-scaled layer model of electronic portal imaging device (EPID) and applying it to a clinical environment. METHODS: The model was developed for fast Monte Carlo dose calculation. The model was validated through comparison with measurements of dose on EPID using first open beams of varying field sizes under a 20-cm-thick flat phantom. After this basic validation, the model was further tested by applying it to transit dosimetry and dose reconstruction that employed our predetermined dose-response-based algorithm developed earlier. The application employed clinical intensity-modulated beams irradiated on a Rando phantom. The clinical beams were obtained through planning on pelvic regions of the Rando phantom simulating prostate and large pelvis intensity modulated radiation therapy. To enhance agreement between calculations and measurements of dose near penumbral regions, convolution conversion of acquired EPID images was alternatively used. In addition, thickness-dependent image-to-dose calibration factors were generated through measurements of image and calculations of dose in EPID through flat phantoms of various thicknesses. The factors were used to convert acquired images in EPID into dose. RESULTS: For open beam measurements, the model showed agreement with measurements in dose difference better than 2% across open fields. For tests with a Rando phantom, the transit dosimetry measurements were compared with forwardly calculated doses in EPID showing gamma pass rates between 90.8% and 98.8% given 4.5 mm distance-to-agreement (DTA) and 3% dose difference (DD) for all individual beams tried in this study. The reconstructed dose in the phantom was compared with forwardly calculated doses showing pass rates between 93.3% and 100% in isocentric perpendicular planes to the beam direction given 3 mm DTA and 3% DD for all beams. On isocentric axial planes, the pass rates varied between 95.8% and 99.9% for all individual beams and they were 98.2% and 99.9% for the composite beams of the small and large pelvis cases, respectively. Three-dimensional gamma pass rates were 99.0% and 96.4% for the small and large pelvis cases, respectively. CONCLUSIONS: The layer model of EPID built for Monte Carlo calculations offered fast (less than 1 min) and accurate calculation for transit dosimety and dose reconstruction.

A method for patient dose reduction in dynamic contrast enhanced CT study.

PURPOSE: In dynamic contrast enhanced CT (DCE-CT) study, prolonged CT scanning with high temporal resolution is required to give accurate and precise estimates of kinetic parameters. However, such scanning protocol could lead to substantial radiation dose to the patient. A novel method is proposed to reduce radiation dose to patient, while maintaining high accuracy for kinetic parameter estimates in DCE-CT study. METHODS: The method is based on a previous investigation that the arterial impulse response (AIR) in DCE-CT study can be predicted using a population-based scheme. In the proposed method, DCE-CT scanning is performed with relatively low temporal resolution, hence, giving rise to reduction in patient dose. A novel method is proposed to estimate the arterial input function (AIF) based on the coarsely sampled AIF. By using the estimated AIF in the tracer kinetic analysis of the coarsely sampled DCE-CT study, the calculated kinetic parameters are able to achieve a high degree of accuracy. The method was tested on a DCE-CT data set of 48 patients with cervical cancer scanned at high temporal resolution. A random cohort of 34 patients was chosen to construct the orthonormal bases of the AIRs via singular value decomposition method. The determined set of orthonormal bases was used to fit the AIFs in the second cohort (14 patients) at varying levels of down sampling. For each dataset in the second cohort, the estimated AIF was used for kinetic analyses of the modified Tofts and adiabatic tissue homogeneity models for each of the down-sampling schemes between intervals from 2 to 15 s. The results were compared with analyses done with the "raw" down-sampled AIF. RESULTS: In the first group of 34 patients, there were 11 orthonormal bases identified to describe the AIRs. The AIFs in the second group were estimated in high accuracy based on the 11 orthonormal bases established in the first group along with down-sampled AIFs. Using the 11 orthonormal bases, the estimated AIFs for the second group were found to have an averaged maximal percentage error of 3.4% ± 7.5% in all sampling schemes up to 15 s. The results of kinetic analysis with the proposed method compared with down sampling alone showed that the proposed method is superior in maintaining the accuracy in volume transfer constant (K(trans) ) after 9 s down-sampling interval, blood volume (v(b) ) for almost all down-sampling intervals, and blood flow (F) after 11 s down-sampling interval. The preliminary results suggested that the proposed method is able to support scanning intervals of 10-15 s at a cost of 6.2%-10.0% loss in accuracy of K(trans) and 10.9%-19.4% in v(b), and the scanning intervals of 12-15 s at a cost of 9.7%-14.6% for F in DEC-CT studies for patients with cervix cancer. CONCLUSIONS: The proposed method of AIF estimation allows low scanning frequency in DCE-CT study to reduce radiation dose to patient, while maintaining relatively high accuracy in the kinetic parameter estimates. The initial results suggested that the method is applicable for DCE-CT studies for patients with cervical cancer.