Enhancement of Hounsfield unit distribution in cone-beam CT images for adaptive radiation therapy: Evaluation of a hybrid correction approach

PurposeThis study aims to evaluate the accuracy of a hybrid approach combining the histogram matching (HM) and the multilevel threshold (MLT) to correct the Hounsfield Unit (HU) distribution in cone-beam CT (CBCT) images.Methods and MaterialsCBCT images acquired for ten prostate cancer patients were processed by matching their histograms to those of deformed planning CT (pCT) images obtained after applying a deformable registration (DR) process. Then, HU values corresponding to five tissue types in the pCT were assigned to the obtained CBCT images (CBCT_HM-MLT). Finally, the CBCT_HM-MLT images were compared to the deformed pCT visually and using different statistical metrics.ResultsThe visual assessment and the profiles comparison showed that the high discrepancies in the CBCT images were significantly reduced when using the proposed approach. Furthermore, the correlation values indicated that the CBCT_HM-MLT were in good agreement with the deformed pCT with correlation values ranging from 0.9893 to 0.9962. In addition, the root mean squared error (RMSE) over the entire volume was reduced from 64.15 ± 9.50 to 51.20 ± 6.76 HU. Similarly, the mean absolute error in specific tissue classes was significantly reduced especially in the soft tissue-air interfaces. These results confirmed that applying MLT after HM worked better than using only HM for which the correlation values were ranging from 0.9878 to 0.9955 and the RMSE was 55.95 ± 10.43 HU.ConclusionEvaluation of the proposed approach showed that the HM + MLT correction can improve the HU distribution in the CBCT images and generate corrected images in good agreement with the pCT.     

Accuracy of target delineation by positron emission tomography-based auto-segmentation methods after deformable image registration: A phantom study

PurposeDiagnostic positron emission tomography and computed tomography (PET/CT) images can be fused to the planning CT images by a deformable image registration (DIR). The aim of this study was to evaluate the standardized uptake value (SUV) and target delineation on deformed PET images.MethodsWe used a cylindrical phantom and removable inserts of four spheres (16–38 mm in diameter) and three ellipsoids with a volume equal to the 38-mm-diameter sphere (S38) in each. S38 was filled with 18F-fluorodeoxyglucose activity, and then PET/CT images were acquired. The contours of S38 were generated using original PET images by PET auto-segmentation (PET-AS) methods of (1) SUV2.5, (2) 40% of maximum SUV (SUV40%max), and (3) gradient-based (GB), and were deformed to the other inserts by DIR. We compared the volumes and the SUV_max with the generated contours using the deformed PET images.ResultsThe SUV_max was slightly decreased by DIR; the mean absolute difference was −0.10 ± 0.04. For SUV2.5 and SUV40%max, the differences in S38 volumes between the original and deformed PET images were less than 5%, regardless of deformation type. For the GB, the contoured volumes obtained from deformed PET images were larger than those of the original PET images for the deformation type of ellipsoids. When the S38 was deformed to the 16-mm-diameter sphere, the maximum volume difference was −22.8%.ConclusionsAlthough SUV fluctuations by DIR were negligible, the target delineation on deformed PET images by the GB should be carefully considered owing to the distortion of intensity profiles.     

Evaluation of a cycle-generative adversarial network-based cone-beam CT to synthetic CT conversion algorithm for adaptive radiation therapy

PurposeImage-guided radiation therapy could benefit from implementing adaptive radiation therapy (ART) techniques. A cycle-generative adversarial network (cycle-GAN)-based cone-beam computed tomography (CBCT)-to-synthetic CT (sCT) conversion algorithm was evaluated regarding image quality, image segmentation and dosimetric accuracy for head and neck (H&N), thoracic and pelvic body regions.MethodsUsing a cycle-GAN, three body site-specific models were priorly trained with independent paired CT and CBCT datasets of a kV imaging system (XVI, Elekta). sCT were generated based on first-fraction CBCT for 15 patients of each body region. Mean errors (ME) and mean absolute errors (MAE) were analyzed for the sCT. On the sCT, manually delineated structures were compared to deformed structures from the planning CT (pCT) and evaluated with standard segmentation metrics. Treatment plans were recalculated on sCT. A comparison of clinically relevant dose-volume parameters (D₉₈, D₅₀ and D₂ of the target volume) and 3D-gamma (3%/3mm) analysis were performed.ResultsThe mean ME and MAE were 1.4, 29.6, 5.4 Hounsfield units (HU) and 77.2, 94.2, 41.8 HU for H&N, thoracic and pelvic region, respectively. Dice similarity coefficients varied between 66.7 ± 8.3% (seminal vesicles) and 94.9 ± 2.0% (lungs). Maximum mean surface distances were 6.3 mm (heart), followed by 3.5 mm (brainstem). The mean dosimetric differences of the target volumes did not exceed 1.7%. Mean 3D gamma pass rates greater than 97.8% were achieved in all cases.ConclusionsThe presented method generates sCT images with a quality close to pCT and yielded clinically acceptable dosimetric deviations. Thus, an important prerequisite towards clinical implementation of CBCT-based ART is fulfilled.     

The effect of tube focal spot size and acquisition mode on task-based image quality performance of a GE revolution HD dual energy CT scanner

PurposeTo assess the task-based performance of images obtained under different focal spot size and acquisition mode on a dual-energy CT scanner.MethodsAxial CT image series of the Catphan phantom were obtained using a tube focus at different sizes. Acquisitions were performed in standard single-energy, high resolution (HR) and dual-energy modes. Images were reconstructed using conventional and high definition (HD) kernels. Task-based transfer function at the 50% level (TTF_50%) for teflon, delrin, low density polyethylene (LDPE) and acrylic, as well as image noise and noise texture, were assessed across all focal spots and acquisition modes using Noise Power Spectrum (NPS) analysis. A non-prewhitening mathematical observer model was used to calculate detectability index ($d_{\mathit{NPW}}^{\prime }$).ResultsTTF_50% degraded with increasing focal spot size. TTF_50% ranged from 0.67 mm⁻¹ for teflon to 0.25 mm⁻¹ for acrylic. For standard kernel, image noise and NPS-determined average spatial frequency were 8.3 HU and 0.29 mm⁻¹, respectively in single-energy, 12.0 HU and 0.37 mm⁻¹ in HR, and 7.9 HU and 0.26 mm⁻¹ in dual-energy mode. For standard kernel, $d_{\mathit{NPW}}^{\prime }$ was 61 in single-energy and HR mode and reduced to 56 in dual-energy mode.ConclusionsThe task-based image quality assessment metrics have shown that spatial resolution is higher for higher image contrast materials and detectability is higher in the standard single-energy mode compared to HR and dual-energy mode. The results of the current study provide CT operators the required knowledge to characterize their CT system towards the optimization of its clinical performance.     

Evaluation of the Block Matching deformable registration algorithm in the field of head-and-neck adaptive radiotherapy

Background and purposeTo compare the accuracy of the Block Matching deformable registration (DIR) against rigid image registration (RIR) for head-and-neck multi-modal images CT to cone-beam CT (CBCT) registration.Material and methodsPlanning-CT and weekly CBCT of 10 patients were used for this study. Several volumes, including medullary canal (MC), thyroid cartilage (TC), hyoid bone (HB) and submandibular gland (SMG) were transposed from CT to CBCT images using either DIR or RIR. Transposed volumes were compared with the manual delineation of these volumes on every CBCT. The parameters of similarity used for analysis were: Dice Similarity Index (DSI), 95%-Hausdorff Distance (95%-HD) and difference of volumes (cc).ResultsWith DIR, the major mean difference of volumes was −1.4 cc for MC, revealing limited under-segmentation. DIR limited variability of DSI and 95%-HD. It significantly improved DSI for TC and HB and 95%-HD for all structures but SMG. With DIR, mean 95%-HD (mm) was 3.01 ± 0.80, 5.33 ± 2.51, 4.99 ± 1.69, 3.07 ± 1.31 for MC, TC, HB and SMG, respectively. With RIR, it was 3.92 ± 1.86, 6.94 ± 3.98, 6.44 ± 3.37 and 3.41 ± 2.25, respectively.ConclusionBlock Matching is a valid algorithm for deformable multi-modal CT to CBCT registration. Values of 95%-HD are useful for ongoing development of its application to the cumulative dose calculation.     

Simultaneous achievement of accurate CT number and image quality improvement for myocardial perfusion CT at 320-MDCT volume scanning

PurposeTo investigate differences in image-to-image variations between full- and half-scan reconstruction on myocardial CT perfusion (CTP) study.MethodsUsing a cardiac phantom we performed ECG-gated myocardial CTP on a second-generation 320-multidetector CT volume scanner. The heart rate was set at 60 bpm; once per second for a total of 24 s were performed. CT images were acquired at 80- and 120 kVp and subjected to full- and half-scan reconstruction. On images acquired at the same slice level we then measured image-to-image variations, coefficients of variance (CV), and image noise.ResultsThe image-to-image variations with full- and half-scan reconstruction were 1.3 HU vs. 27.2 HU at 80 kVp (p < 0.001) and 0.70 HU vs. 9.3 HU at 120 kVp (p < 0.001) even though the mean HU value was almost the same for both reconstruction methods. The CV of 80- and 120-kVp images of the left ventricular cavity decreased by 0.16% and 0.17%, respectively, with full-scan reconstruction; with half-scan reconstruction it decreased by 3.34% and 2.30%, respectively. Compared with half-scan reconstruction, the image noise was reduced by 27.2% at 80 kVp and by 28.0% at 120 kVp with full-scan reconstruction.ConclusionMyocardial CTP with full-scan reconstruction substantially decreased image-to-image variations and provided accurate CT attenuation.     

Dose warping performance in deformable image registration in lung

PurposeIt is unclear that spatial accuracy can reflect the impact of deformed dose distribution. In this study, we used dosimetric parameters to compare an in-house deformable image registration (DIR) system using NiftyReg, with two commercially available systems, MIM Maestro (MIM) and Velocity AI (Velocity).MethodsFor 19 non-small-cell lung cancer patients, the peak inspiration (0%)-4DCT images were deformed to the peak expiration (50%)-4DCT images using each of the three DIR systems, which included computation of the deformation vector fields (DVF). The 0%-gross tumor volume (GTV) and the 0%-dose distribution were also then deformed using the DVFs. The agreement in the dose distributions for the GTVs was evaluated using generalized equivalent uniform dose (gEUD), mean dose (D_mean), and three-dimensional (3D) gamma index (criteria: 3 mm/3%). Additionally, a Dice similarity coefficient (DSC) was used to measure the similarity of the GTV volumes.ResultsD_mean and gEUD demonstrated good agreement between the original and deformed dose distributions (differences were generally less than 3%) in 17 of the patients. In two other patients, the Velocity system resulted in differences in gEUD of 50.1% and 29.7% and in D_mean of 11.8% and 4.78%. The gamma index comparison showed statistically significant differences for the in-house DIR vs. MIM, and MIM vs. Velocity.ConclusionsThe finely tuned in-house DIR system could achieve similar spatial and dose accuracy to the commercial systems. Care must be taken, as we found errors of more than 5% for D_mean and 30% for gEUD, even with a commercially available DIR tool.     

Effects of lung tissue characterization in radiotherapy of breast cancer under deep inspiration breath hold when using Monte Carlo dosimetry

PurposeTo investigate the sensitivity of Monte Carlo (MC) calculated lung dose distributions to lung tissue characterization in external beam radiotherapy of breast cancer under Deep Inspiration Breath Hold (DIBH).MethodsEGSnrc based MC software was employed. Mean lung densities for one hundred patients were analysed. CT number frequency and clinical dose distributions were calculated for 15 patients with mean lung density below 0.14 g/cm³. Lung volume with a pre-defined CT numbers was also considered. Lung tissue was characterized by applying different CT calibrations in the low-density region and air-lung tissue thresholds. Dose impact was estimated by Dose Volume Histogram (DVH) parameters.ResultsMean lung densities below 0.14 g/cm³ were found in 10% of the patients. CT numbers below −960 HU dominated the CT frequency distributions with a high rate of CT numbers at −990 HU. Mass density conversion approach influenced the DVH shape. V_4Gy and V_8Gy varied by 7% and 5% for the selected patients and by 9% and 3.5% for the pre-defined lung volume. V_16Gy and V_20Gy, were within 2.5%. Regions above 20 Gy were affected. Variations in air- lung tissue differentiation resulted in DVH parameters within 1%. Threshold at −990 HU was confirmed by the CT number frequency distributions.ConclusionsLung dose distributions were more sensitive to variations in the CT calibration curve below lung (inhale) density than to air-lung tissue differentiation. Low dose regions were mostly affected. The dosimetry effects were found to be potentially important to 10% of the patients treated under DIBH.     

Evaluation of four-dimensional cone beam computed tomography ventilation images acquired with two different linear accelerators at various gantry speeds using a deformable lung phantom

We evaluated four-dimensional cone beam computed tomography (4D-CBCT) ventilation images (VI_CBCT) acquired with two different linear accelerator systems at various gantry speeds using a deformable lung phantom.The 4D-CT and 4D-CBCT scans were performed using a computed tomography (CT) scanner, an X-ray volume imaging system (Elekta XVI) mounted in Versa HD, and an On-Board Imager (OBI) system mounted in TrueBeam. Intensity-based deformable image registration (DIR) was performed between peak-exhale and peak-inhale images. VI_CBCT- and 4D-CT-based ventilation images (VI_CT) were derived by DIR using two metrics: one based on the Jacobian determinant and one on changes in the Hounsfield unit (HU). Three different DIR regularization values (λ) were used for VI_CBCT. Correlations between the VI_CBCT and VI_CT values were evaluated using voxel-wise Spearman’s rank correlation coefficient (r).In case of both metrics, the Jacobian-based VI_CBCT with a gantry speed of 0.6 deg/sec in Versa HD showed the highest correlation for all the gantry speeds (e.g., λ = 0.05 and r = 0.68). Thus, the r value of the Jacobian-based VI_CBCT was greater or equal to that of the HU-based VI_CBCT. In addition, the ventilation accuracy of VI_CBCT increased at low gantry speeds.Thus, the image quality of VI_CBCT was affected by the change in gantry speed in both the imaging systems. Additionally, DIR regularization considerably influenced VI_CBCT in both the imaging systems. Our results have the potential to assist in designing CBCT protocols, incorporating VI_CBCT imaging into the functional avoidance planning process.     

Dose calculation accuracy of different image value to density tables for cone-beam CT planning in head & neck and pelvic localizations

PurposeWe aimed to identify the most accurate combination of phantom and protocol for image value to density table (IVDT) on volume-modulated arc therapy (VMAT) dose calculation based on kV-Cone-beam CT imaging, for head and neck (H&N) and pelvic localizations.MethodsThree phantoms (Catphan^®600, CIRS^®062M (inner phantom for head and outer phantom for body), and TomoTherapy^® “Cheese” phantom) were used to create IVDT curves of CBCT systems with two different CBCT protocols (Standard-dose Head and Standard Pelvis). Hounsfield Unit (HU) time stability and repeatability for a single On-Board-Imager (OBI) and compatibility of two distinct devices were assessed with Catphan^®600. Images from the anthropomorphic phantom CIRS ATOM^® for both CT and CBCT modalities were used for VMAT dose calculation from different IVDT curves. Dosimetric indices from CT and CBCT imaging were compared.ResultsIVDT curves from CBCT images were highly different depending on phantom used (up to 1000 HU for high densities) and protocol applied (up to 200 HU for high densities). HU time stability was verified over seven weeks. A maximum difference of 3% on the dose calculation indices studied was found between CT and CBCT VMAT dose calculation across the two localizations using appropriate IVDT curves. One IVDT curve per localization can be established with a bi-monthly verification of IVDT-CBCT.ConclusionsThe IVDT-CBCT_CIRS-Head phantom with the Standard-dose Head protocol was the most accurate combination for dose calculation on H&N CBCT images. For pelvic localizations, the IVDT-CBCT_Cheese established with the Standard Pelvis protocol provided the best accuracy.     

Can different Catphan phantoms be used in a multi-centre audit of radiotherapy CT image quality?

PurposeTo determine the variation between Catphan image quality CT phantoms, specifically for use in a future multi-centre image quality audit.Method14 Catphan phantoms (models 503, 504 and 604) were scanned on a Canon Aquilion Prime CT scanner using a single scan protocol. Measurements were made of noise in the uniformity section, visibility of low contrast targets and contrast, x-ray attenuation and CT number for 5 materials in the sensitometry section. Scans were also acquired using one phantom and varying reconstruction field of view, image slice thickness, effective tube-current-time product and iterative reconstruction settings to determine how the degree of inter-phantom variability compared with the magnitude of changes from scan parameter alteration.ResultsAcross all phantoms the mean CT value in the uniformity section was 7.0 (SD 0.9) range: 4.9–8.1 HU. For the different materials the CT numbers were air: −1004 ± 5, Polymethylpentene: −190 ± 2, Polystyrene: −42 ± 2, Delrin: 321 ± 5 and Teflon: 898 ± 8 HU. Consistency of low contrast targets through visual scoring was good. Measured contrast was lower (p < 0.001) with more variability for 504 versus 604 models. All phantoms produced identical tube current settings with x-ray tube current modulation, indicating no x-ray attenuation differences. The degree of change in image quality metrics between phantoms was small compared with results when scan parameters were varied.ConclusionCatphan phantoms model 604 showed minimal differences and will be used for multi-centre inter-comparison work, with the consistency between phantoms appropriate for measuring possible variations in image quality.     

The effectiveness of gadolinium MRI to improve target delineation for radiotherapy in hepatocellular carcinoma: A comparative study of rigid image registration techniques

To achieve consistent target delineation in radiotherapy for hepatocellular carcinoma (HCC), image registration between simulation CT and diagnostic MRI was explored.Twenty patients with advanced HCC were included. The median interval between MRI and CT was 11 days. CT was obtained with shallow free breathing and MRI at exhale phase. On each CT and MRI, the liver and the gross target volume (GTV) were drawn. A rigid image registration was taken according to point information of vascular bifurcation (Method[A]) and pixel information of volume of interest only including the periphery of the liver (Method[B]) and manually drawn liver (Method[C]). In nine cases with an indefinite GTV on CT, a virtual sphere was generated at the epicenter of the GTV. The GTV from CT (V_GTV[CT]) and MRI (V_GTV[MR]) and the expanded GTV from MRI (V_+GTV[MR]) considering geometrical registration error were defined. The underestimation (uncovered V[CT] by V[MR]) and the overestimation (excessive V[MR] by V[CT]) were calculated. Through a paired T-test, the difference between image registration techniques was analyzed.For method[A], the underestimation rates of V_GTV[MR] and V_+GTV[MR] were 16.4 ± 8.9% and 3.2 ± 3.7%, and the overestimation rates were 16.6 ± 8.7% and 28.4 ± 10.3%, respectively. For V_GTV[MR] and V_+GTV[MR], the underestimation rates and overestimation rates of method[A] were better than method[C]. The underestimation rates and overestimation rates of the V_GTV[MR] were better in method[B] than method[C]. By image registration and additional margin, about 97% of HCC could be covered. Method[A] or method[B] could be recommended according to physician preference.     

A new strategy for craniospinal axis localization and adaptive dosimetric evaluation using cone beam CT

Background and AimComputational complexities encountered in craniospinal irradiation (CSI) have been widely investigated with different planning strategies. However, localization of the entire craniospinal axis (CSA) and evaluation of adaptive treatment plans have traditionally been ignored in CSI treatment. In this study, a new strategy for CSI with comprehensive CSA localization and adaptive plan evaluation has been demonstrated using cone beam CT with extended longitudinal field-of-view (CBCT_eLFOV).Materials and MethodsMulti-scan CBCT images were acquired with fixed longitudinal table translations (with 1 cm cone-beam overlap) and then fused into a single DICOM-set using the custom software coded in MatLab™. A novel approach for validation of CBCT_eLFOV was demonstrated by combined geometry of Catphan-504 and Catphan-604 phantoms. To simulate actual treatment scenarios, at first, the end-to-end workflow of CSI with VMAT was investigated using an anthropomorphic phantom and then applied for two patients (based on random selection).ResultsThe fused CBCT_eLFOV images were in excellent agreement with planning CT (pCT). The custom developed software effectively manages spatial misalignments arising out of the uncertainties in treatment/setup geometry. Although the structures mapped from pCT to CBCT_eLFOV showed minimal variations, a maximum spatial displacement of up to 1.2 cm (and the mean of 0.8 ± 0.3 cm) was recorded in phantom study. Adaptive plan evaluation of patient paradigms showed the likelihood of under-dosing the craniospinal target.ConclusionOur protocol serves as a guide for precise localization of entire CSA and to ensure adequate dose to the large and complex targets. It can also be adapted for other complex treatment techniques such as total-marrow-irradiation and total-lymphoid-irradiation.     

Machine log file-based dose verification using novel iterative CBCT reconstruction algorithm in commercial software during volumetric modulated arc therapy for prostate cancer patients

PurposeTo evaluate the utility of the use of iterative cone-beam computed tomography (CBCT) for machine log file-based dose verification during volumetric modulated arc therapy (VMAT) for prostate cancer patients.MethodsAll CBCT acquisition data were used to reconstruct images with the Feldkamp-Davis-Kress algorithm (FDK-CBCT) and the novel iterative algorithm (iCBCT). The Hounsfield unit (HU)-electron density curves for CBCT images were created using the Advanced Electron Density Phantom. The I’mRT and anthropomorphic phantoms were irradiated with VMAT after CBCT registration. Subsequently, fourteen prostate cancer patients received VMAT after CBCT registration. Machine log files and both CBCT images were exported to the PerFRACTION software, and a 3D patient dose was reconstructed. Mean dose for planning target volume (PTV), the bladder, and rectum and the 3D gamma analysis were evaluated.ResultsFor the phantom studies, the variation of HU values was observed at the central position surrounding the bones in FDK-CBCT. There were almost no changes in the difference of doses at the isocenter between measurement and reconstructed dose for planning CT (pCT), FDK-CBCT, and iCBCT. Mean dose differences of PTV, rectum, and bladder between iCBCT and pCT were approximately 2% lower than those between FDK-CBCT and pCT. For the clinical study, average gamma analysis for 2%/2 mm was 98.22% ± 1.07 and 98.81% ± 1.25% in FDK-CBCT and iCBCT, respectively.ConclusionsA similar machine log file-based dose verification accuracy is obtained for FDK-CBCT and iCBCT during VMAT for prostate cancer patients.     

A steady state model of sequential irreversible enzyme reactions

Summary One of the questions which arises in the study of certain inborn errors of metabolism as well as in the field of enzyme kinetics is: what are the quantitative relationships between parameters of enzyme activity and substrate pool sizes in a metabolic pathway? A steady state model has been devised to answer this question for a homogeneous system of non-branched sequential irreversible enzyme reactions which follow Michaelis-Menten kinetics. The concentration of a substrate in such a pathway, [S_i], is a function of 5 variables: (a) the K_M of the enzyme which forms the substrate (K_{M (i–1)}), (b) the K_M of the enzyme which utilizes the substrate (K_{M i}), (c) the V_max of the enzyme which forms the substrate (V_{m (i–1)}), (d) the V_max of the enzyme which utilizes the substrate (V_{m i}) and (e) the immediate precursor concentration [S_(i–1)] where [S_i] = K_{M i} V_{m (i–1)} [S_(i–1)]/[S_(i–1)] (V_mi -V_{m (i–1)}) + K_{M (i–1)}) V_mi The model introduces and defines the concept of and conditions for amplification. An input in the form of a steady state concentration of precursor [S_(i–1)] may be amplified as an output in the form of an increased steady state concentration of product [S_i]. The model also defines the values of the above 5 parameters which do not allow attainment of a steady state for the type of pathway considered.From the Metabolism Branch, National Cancer Institute, National Institutes of Health, Bethesda, Maryland 20014.     

Quantitative evaluation of internal marks made using MRgFUS as seen on MRI,CT, US,and digital color images – A pilot study

This pilot study compared the detectability of internal thermal marks produced with MRI-guided focused ultrasound (MRgFUS) on MRI, computed tomography (CT), ultrasonography (US), and color images from digital scanning. Internal marks made using MRgFUS could potentially guide surgical, biopsy or radiotherapy procedures. New Zealand White rabbits (n = 6) thigh muscle were marked using a Philips MRgFUS system. Before and after sonications, rabbits were imaged using T₁- and T₂-weighted MRI. Then rabbits were sacrificed and imaging was performed using CT and US. After surgical excision specimens were scanned for color conspicuity analysis. Images were read by a radiologist and quantitative analysis of signal intensity was calculated for marks and normal muscle. Of a total of 19 excised marks, approximately 79%, 63%, and 62% were visible on MRI, CT, and US, respectively. The average maximum temperature elevation in the marks during MRgFUS was 39.7 ± 10.1 °C, and average dose diameter (i.e., the diameter of the area that achieved a thermal dose greater than 240 cumulative equivalent minutes at 43 °C) of the mark at the focal plane was 7.3 ± 2.1 mm. On MRI the average normalized signal intensities were significantly higher in marks compared to normal muscle (p < 0.05). On CT, the marked regions were approximately 10 HU lower than normal muscle (p < 0.05). The results demonstrate that MRgFUS can be used to create internal marks that are visible on MRI, CT and US.     

Commissioning and validation of commercial deformable image registration software for adaptive contouring

PurposeTo report the commissioning and validation of deformable image registration(DIR) software for adaptive contouring.MethodsDIR (SmartAdapt®v13.6) was validated using two methods namely contour propagation accuracy and landmark tracking, using physical phantoms and clinical images of various disease sites. Five in-house made phantoms with various known deformations and a set of 10 virtual phantoms were used. Displacement in lateral, anterio-posterior (AP) and superior-inferior (SI) direction were evaluated for various organs and compared with the ground truth. Four clinical sites namely, brain (n = 5), HN (n = 9), cervix (n = 18) and prostate (n = 23) were used. Organs were manually delineated by a radiation oncologist, compared with the deformable image registration (DIR) generated contours. 3D slicer v4.5.0.1 was used to analyze Dice Similarity Co-efficient (DSC), shift in centre of mass (COM) and Hausdorff distances Hf_95%/avg.ResultsMean (SD) DSC, Hf_95% (mm), Hf_avg (mm) and COM of all the phantoms 1–5 were 0.84 (0.2) mm, 5.1 (7.4) mm, 1.6 (2.2) mm, and 1.6 (0.2) mm respectively. Phantom-5 had the largest deformation as compared to phantoms 1–4, and hence had suboptimal indices. The virtual phantom resulted in consistent results for all the ROIs investigated. Contours propagated for brain patients were better with a high DSC score (0.91 (0.04)) as compared to other sites (HN: 0.84, prostate: 0.81 and cervix 0.77). A similar trend was seen in other indices too. The accuracy of propagated contours is limited for complex deformations that include large volume and shape change of bladder and rectum respectively. Visual validation of the propagated contours is recommended for clinical implementation.ConclusionThe DIR algorithm was commissioned and validated for adaptive contouring.     

Electrophysiologic Characteristics of the Crista Terminalis and Implications on Atrial Tachycardia in Rabbits

The objective of this study was to evaluate the electrophysiologic characteristics of Crista terminalis (CT) and their implication in the pathogenesis of atrial tachycardia in rabbits. For this purpose, 27 New Zealand rabbits were used. Using standard glass microelectrode technique, cellular action potentials (APs) of CT and pectinate muscle (PM) were recorded in normal Tyrode’s perfusion and Tyrode’s perfusion with 4 μM isoproterenol. Longitudinal conduction velocity (V _L) and transverse conduction velocity (V _T) of CT were measured. As our data show, CT tissue had a trend of spontaneous phase IV depolarization. Conduction anisotropy (V _L/V _T) of CT was 4.53 ± 0.91. The duration of the AP of CT was longer than that of PM cells. APD₂₀ and APD₉₀ for CT were 28.1 ± 3.5 and 145.3 ± 7.1 ms; and for PM cells were 21.8 ± 4.1 and 125.3 ± 6.3 ms, respectively (all P values < 0.01). The early and delayed action depolarizations were recorded after isoproterenol perfusion. A fast paroxysmal irregular rhythm was recorded which could be arrested by 0.1 mmol/l Isoptin. It was, therefore, concluded that the latent autorhythmicity, trigger activity, and conduction properties of CT might provide the electrophysiologic basis for the occurrence and sustenance of atrial arrhythmia.     

Optimization of image quality and accuracy of low iodine concentration quantification as function of kVp pairs for abdominal imaging using dual-source CT: A phantom study

PurposeTo determine the suitable kVp pair for optimal image quality of the virtual monochromatic images (VMIs) and iodine quantification accuracy at low concentration, using a third generation dual-source CT (DSCT).Materials and methodsMulti-energy CT phantoms with and without body rings were scanned with a DSCT using four kVp pairs (tube “A”/“B” voltage): 100/Sn150, 90/Sn150, 80/Sn150 and 70/Sn150 kVp. The reference mAs was adjusted to obtain a CTDI_vol close to 11 mGy. HU values accuracy (RMSD_HU), noise (SD) and contrast-to-noise ratio (CNR) of iodine inserts of 0.5, 1, 2 and 5 mg/mL concentrations were assessed on VMIs at 40/50/60/70 keV. Iodine quantification accuracy was assessed using the RMSD_iodine and iodine bias (IB_iodine).ResultsThe RMSD_HU decreased when the tube “A” voltage increased. The mean noise value increased significantly with tube “A” voltage (p < 0.001) but decreased between 80/Sn150 and 90/Sn150 kVp for the small phantom (1.1 ± 0.1%; p = 0.047). The CNR significantly decreased with tube “A” voltage (p < 0.001), except between 80/Sn150 and 90/Sn150 kVp for all inserts and between 90/Sn150 kVp and 100/Sn150 kVp for the 1.0 and 0.5 mg/mL inserts in the large phantom. In the small phantom, no significant difference was found between 80/Sn150 kVp and 90/Sn150 kVp for all inserts and between 80/Sn150, 90/Sn150 and 100/Sn150 kVp for the 1 and 0.5 mg/mL inserts. The RMSD_iodine and IB_iodine decreased as the tube “A” voltage of the kVp pair increased.ConclusionThe kVp pair of 70/Sn150 led to better image quality in VMIs and sufficient iodine accuracy.     

Zero echo time MRI-only treatment planning for radiation therapy of brain tumors after resection

Using magnetic resonance imaging (MRI) as the sole imaging modality for patient modeling in radiation therapy (RT) is a challenging task due to the need to derive electron density information from MRI and construct a so-called pseudo-computed tomography (pCT) image. We have previously published a new method to derive pCT images from head T1-weighted (T1-w) MR images using a single-atlas propagation scheme followed by a post hoc correction of the mapped CT numbers using local intensity information. The purpose of this study was to investigate the performance of our method with head zero echo time (ZTE) MR images. To evaluate results, the mean absolute error in bins of 20 HU was calculated with respect to the true planning CT scan of the patient. We demonstrated that applying our method using ZTE MR images instead of T1-w improved the correctness of the pCT in case of bone resection surgery prior to RT (that is, an example of large anatomical difference between the atlas and the patient).