T. J. Arai

A continuous surface reconstruction method on point cloud captured from a 3D surface photogrammetry system.

PURPOSE To accurately and efficiently reconstruct a continuous surface from noisy point clouds captured by a surface photogrammetry system (VisionRT). METHODS The authors have developed a level-set based surface reconstruction method on point clouds captured by a surface photogrammetry system (VisionRT). The proposed method reconstructs an implicit and continuous representation of the underlying patient surface by optimizing a regularized fitting energy, offering extra robustness to noise and missing measurements. By contrast to explicit/discrete meshing-type schemes, their continuous representation is particularly advantageous for subsequent surface registration and motion tracking by eliminating the need for maintaining explicit point correspondences as in discrete models. The authors solve the proposed method with an efficient narrowband evolving scheme. The authors evaluated the proposed method on both phantom and human subject data with two sets of complementary experiments. In the first set of experiment, the authors generated a series of surfaces each with different black patches placed on one chest phantom. The resulting VisionRT measurements from the patched area had different degree of noise and missing levels, since VisionRT has difficulties in detecting dark surfaces. The authors applied the proposed method to point clouds acquired under these different configurations, and quantitatively evaluated reconstructed surfaces by comparing against a high-quality reference surface with respect to root mean squared error (RMSE). In the second set of experiment, the authors applied their method to 100 clinical point clouds acquired from one human subject. In the absence of ground-truth, the authors qualitatively validated reconstructed surfaces by comparing the local geometry, specifically mean curvature distributions, against that of the surface extracted from a high-quality CT obtained from the same patient. RESULTS On phantom point clouds, their method achieved submillimeter reconstruction RMSE under different configurations, demonstrating quantitatively the faith of the proposed method in preserving local structural properties of the underlying surface in the presence of noise and missing measurements, and its robustness toward variations of such characteristics. On point clouds from the human subject, the proposed method successfully reconstructed all patient surfaces, filling regions where raw point coordinate readings were missing. Within two comparable regions of interest in the chest area, similar mean curvature distributions were acquired from both their reconstructed surface and CT surface, with mean and standard deviation of (μrecon=-2.7×10(-3) mm(-1), σrecon=7.0×10(-3) mm(-1)) and (μCT=-2.5×10(-3) mm(-1), σCT=5.3×10(-3) mm(-1)), respectively. The agreement of local geometry properties between the reconstructed surfaces and the CT surface demonstrated the ability of the proposed method in faithfully representing the underlying patient surface. CONCLUSIONS The authors have integrated and developed an accurate level-set based continuous surface reconstruction method on point clouds acquired by a 3D surface photogrammetry system. The proposed method has generated a continuous representation of the underlying phantom and patient surfaces with good robustness against noise and missing measurements. It serves as an important first step for further development of motion tracking methods during radiotherapy.

MR-CBCT image-guided system for radiotherapy of orthotopic rat prostate tumors

Multi-modality image-guided radiotherapy is the standard of care in contemporary cancer management; however, it is not common in preclinical settings due to both hardware and software limitations. Soft tissue lesions, such as orthotopic prostate tumors, are difficult to identify using cone beam computed tomography (CBCT) imaging alone. In this study, we characterized a research magnetic resonance (MR) scanner for preclinical studies and created a protocol for combined MR-CBCT image-guided small animal radiotherapy. Two in-house dual-modality, MR and CBCT compatible, phantoms were designed and manufactured using 3D printing technology. The phantoms were used for quality assurance tests and to facilitate end-to-end testing for combined preclinical MR and CBCT based treatment planning. MR and CBCT images of the phantoms were acquired utilizing a Varian 4.7 T scanner and XRad-225Cx irradiator, respectively. The geometry distortion was assessed by comparing MR images to phantom blueprints and CBCT. The corrected MR scans were co-registered with CBCT and subsequently used for treatment planning. The fidelity of 3D printed phantoms compared to the blueprint design yielded favorable agreement as verified with the CBCT measurements. The geometric distortion, which varied between -5% and 11% throughout the scanning volume, was substantially reduced to within 0.4% after correction. The distortion free MR images were co-registered with the corresponding CBCT images and imported into a commercial treatment planning software SmART Plan. The planning target volume (PTV) was on average 19% smaller when contoured on the corrected MR-CBCT images relative to raw images without distortion correction. An MR-CBCT based preclinical workflow was successfully designed and implemented for small animal radiotherapy. Combined MR-CBCT image-guided radiotherapy for preclinical research potentially delivers enhanced relevance to human radiotherapy for various disease sites. This novel protocol is wide-ranging and not limited to the orthotopic prostate tumor study presented in the study.

Comparison of Quantitative Multiple Breath Specific Ventilation Imaging Using Co-localized 2D Oxygen Enhanced MRI and Hyperpolarized 3He MRI.

Two magnetic resonance specific ventilation imaging (SVI) techniques, namely, oxygen-enhanced proton (OE-1H) and hyperpolarized 3He (HP-3He), were compared in eight healthy supine subjects [age 32 (6) yr]. An in-house radio frequency coil array for 1H configured with the 3He transmit-receive coil in situ enabled acquisition of SVI data from two nuclei from the same slice without repositioning the subjects. After 3 × 3 voxel downsampling to account for spatial registration errors between the two SV images, the voxel-by-voxel correlation coefficient of two SV maps ranged from 0.11 to 0.63 [0.46 mean (0.17 SD); P < 0.05]. Several indexes were analyzed and compared from the tidal volume-matched SV maps: the mean of SV log-normal distribution (SVmean), the standard deviation of the distribution as a measure of SV heterogeneity (SVwidth), and the gravitational gradient (SVslope). There were no significant differences in SVmean [OE-1H: 0.28 (0.08) and HP-3He: 0.32 (0.14)], SVwidths [OE-1H: 0.28 (0.08) and HP-3He: 0.27 (0.10)], and SVslopes [OE-1H: -0.016 (0.006) cm-1 and HP-3He: -0.013 (0.007) cm-1]. Despite the statistical similarities of the population averages, Bland-Altman analysis demonstrated large individual intertechnique variability. SDs of differences in these indexes were 42% (SVmean), 46% (SVwidths), and 62% (SVslopes) of their corresponding overall mean values. The present study showed that two independent, spatially coregistered, SVI techniques presented a moderate positive voxel-by-voxel correlation. Population averages of SVmean, SVwidth, and SVslope were in close agreement. However, the lack of agreement when the data sets were analyzed individually might indicate some fundamental mechanistic differences between the techniques. NEW & NOTEWORTHY To the best of our knowledge, this is the first cross-comparison of two different specific ventilation (SV) MRI techniques in the human lung (i.e., oxygen-enhanced proton and hyperpolarized 3He). The present study showed that two types of spatially coregistered SV images presented a modest positive correlation. The two techniques also yielded similar population averages of SV indexes such as log-normal mean, SV heterogeneity, and the gravitational slope, albeit with some intersubject variability.

SU-D-207A-05: Investigating Sparse-Sampled MRI for Motion Management in Thoracic Radiotherapy

PURPOSE Sparse sampling and reconstruction-based MRI techniques represent an attractive strategy to achieve sufficiently high image acquisition speed while maintaining image quality for the task of radiotherapy guidance. In this study, we examine rapid dynamic MRI using a sparse sampling sequence k-t BLAST in capturing motion-induced, cycle-to-cycle variations in tumor position. We investigate the utility of long-term MRI-based motion monitoring as a means of better characterizing respiration-induced tumor motion compared to a single-cycle 4DCT. METHODS An MRI-compatible, programmable, deformable lung motion phantom with eleven 1.5 ml water marker tubes was placed inside a 3.0 T whole-body MR scanner (Philips Ingenia). The phantom was programmed with 10 lung tumor motion traces previously recorded using the Synchrony system. 2D+t image sequences of a coronal slice were acquired using a balanced-SSFP sequence combined with k-t BLAST (accn=3, resolution=0.66×0.66×5 mm3; acquisition time = 110 ms/slice). kV fluoroscopic (ground truth) and 4DCT imaging was performed with the same phantom setup and motion trajectories. Marker positions in all three modalities were segmented and tracked using an opensource deformable image registration package, NiftyReg. RESULTS Marker trajectories obtained from rapid MRI exhibited <1 mm error compared to kv Fluoro trajectories in the presence of complex motion including baseline shifts and changes in respiratory amplitude, indicating the ability of MRI to monitor motion with adequate geometric fidelity for the purpose of radiotherapy guidance. In contrast, the trajectory derived from 4DCT exhibited significant errors up to 6 mm due to cycle-to-cycle variations and baseline shifts. Consequently, 4DCT was found to underestimate the range of marker motion by as much as 50%. CONCLUSION Dynamic MRI is a promising tool for radiotherapy motion management as it permits for longterm, dose-free, soft-tissue-based monitoring of motion, yielding richer and more accurate information about tumor position and motion range compared to the current state-of-the-art, 4DCT. This work was partially supported through research funding from National Institutes of Health (R01CA169102).

WE-G-BRD-02: Characterizing Information Loss in a Sparse-Sampling-Based Dynamic MRI Sequence (k-T BLAST) for Lung Motion Monitoring

Purpose: Rapid MRI is an attractive, non-ionizing tool for soft-tissue-based monitoring of respiratory motion in thoracic and abdominal radiotherapy. One big challenge is to achieve high temporal resolution while maintaining adequate spatial resolution. K-t BLAST, sparse-sampling and reconstruction sequence based on a-priori information represents a potential solution. In this work, we investigated how much “true” motion information is lost as a-priori information is progressively added for faster imaging. Methods: Lung tumor motions in superior-inferior direction obtained from ten individuals were replayed into an in-house, MRI-compatible, programmable motion platform (50Hz refresh and 100microns precision). Six water-filled 1.5ml tubes were placed on it as fiducial markers. Dynamic marker motion within a coronal slice (FOV: 32×32cm2, resolution: 0.67×0.67mm2, slice-thickness: 5mm) was collected on 3.0T body scanner (Ingenia, Philips). Balanced-FFE (TE/TR: 1.3ms/2.5ms, flip-angle: 40degrees) was used in conjunction with k-t BLAST. Each motion was repeated four times as four k-t acceleration factors 1, 2, 5, and 16 (corresponding frame rates were 2.5, 4.7, 9.8, and 19.1Hz, respectively) were compared. For each image set, one average motion trajectory was computed from six marker displacements. Root mean square error (RMS) was used as a metric of spatial accuracy where measured trajectories were compared to original data. Results: Tumor motion was approximately 10mm. The mean(standard deviation) of respiratory rates over ten patients was 0.28(0.06)Hz. Cumulative distributions of tumor motion frequency spectra (0–25Hz) obtained from the patients showed that 90% of motion fell on 3.88Hz or less. Therefore, the frame rate must be a double or higher for accurate monitoring. The RMS errors over patients for k-t factors of 1, 2, 5, and 16 were.10(.04),.17(.04), .21(.06) and.26(.06)mm, respectively. Conclusions: K-t factor of 5 or higher can cover the high frequency component of tumor respiratory motion, while the estimated error of spatial accuracy was approximately.2mm.