Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “guide”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,495 records · Page 83

SU-E-J-151: Dosimetric Evaluation of DIR Mapped Contours for Image Guided Adaptive Radiotherapy with 4D Cone-Beam CT

Purpose: To estimate dosimetric errors resulting from using contours deformably mapped from planning CT to 4D cone beam CT (CBCT) images for image-guided adaptive radiotherapy of locally advanced non-small cell lung cancer (NSCLC). Methods: Ten locally advanced non-small cell lung cancer (NSCLC) patients underwent one planning 4D fan-beam CT (4DFBCT) and weekly 4DCBCT scans. Multiple physicians delineated the gross tumor volume (GTV) and normal structures in planning CT images and only GTV in CBCT images. Manual contours were mapped from planning CT to CBCTs using small deformation, inverse consistent linear elastic (SICLE) algorithm for two scans in each patient. Two physicians reviewed and rated the DIR-mapped (auto) and manual GTV contours as clinically acceptable (CA), clinically acceptable after minor modification (CAMM) and unacceptable (CU). Mapped normal structures were visually inspected and corrected if necessary, and used to override tissue density for dose calculation. CTV (6mm expansion of GTV) and PTV (5mm expansion of CTV) were created. VMAT plans were generated using the DIR-mapped contours to deliver 66 Gy in 33 fractions with 95% and 100% coverage (V66) to PTV and CTV, respectively. Plan evaluation for V66 was based on manual PTV and CTV contours. Results: Mean PTV V66 was 84% (range 75% – 95%) and mean CTV V66 was 97% (range 93% – 100%) for CAMM scored plans (12 plans); and was 90% (range 80% – 95%) and 99% (range 95% – 100%) for CA scored plans (7 plans). The difference in V66 between CAMM and CA was significant for PTV (p = 0.03) and approached significance for CTV (p = 0.07). Conclusion: The quality of DIR-mapped contours directly impacted the plan quality for 4DCBCT-based adaptation. Larger safety margins may be needed when planning with auto contours for IGART with 4DCBCT images. Reseach was supported by NIH P01CA116602.

60 APPLIED LIFE SCIENCES↗

SU-E-J-162: Quality Assurance Procedures for MR Guided Focused Ultrasound Treatment of Bone Metastasis

Purpose: The purpose of this work is to develop and verify our quality assurance (QA) procedures to ensure the safety and efficacy of MR-guided focused ultrasound (MRgFUS) treatment of bone metastases. Methods: A practical QA program was developed. Monthly and daily QA (DQA) procedures were performed. The major QA items included the checks of the machine hardware, software and patient safety features. Briefly, these checks/tests include: 1) the cooling system reservoir and treatment table; 2) power to the treatment table; 3) the MR coil; 4) the transducer position with MRI; 5) image display on the treatment work station; 6) the effective focal spot in 3 directions using MR thermometry; and 7) all the safety devices including a sonication lamp, and the emergency stop-sonication switches. In order to avoid patient skin burn, it is important to remove gas bubbles in the interfaces between the treatment table and the gel pad, and the gel pad and patients skin during the patient setup. Our QA procedures have been verified and evaluated through patient treatments. Seven patients with scapula, humeral head, sacrum, ilium, pubic ramus and acetabular bone metastases were treated using MRgFUS. Results: Our study showed that all seven patients tolerated the MRgFUS treatment well. No skin toxicity or other complications were observed. The pain score (0–10) using the visual analog scale (VAS) was significantly reduced from 8.0 ± 1.1 before treatment to 4.7 ± 3.0, 3.0 ± 1.5, 3.2 ± 2.8 and 3.4 ± 1.5 at one day, one month, two months and three months after the MRgFUS treatment, respectively. Conclusion: We demonstrated that with the appropriate QA procedures, MRgFUS is a safe, effective and noninvasive treatment modality for palliation of bone metastases.

60 APPLIED LIFE SCIENCES↗

SU-E-T-154: Establishment and Implement of 3D Image Guided Brachytherapy Planning System

Purpose: Cannot observe the dose intuitionally is a limitation of the existing 2D pre-implantation dose planning. Meanwhile, a navigation module is essential to improve the accuracy and efficiency of the implantation. Hence a 3D Image Guided Brachytherapy Planning System conducting dose planning and intra-operative navigation based on 3D multi-organs reconstruction is developed. Methods: Multi-organs including the tumor are reconstructed in one sweep of all the segmented images using the multiorgans reconstruction method. The reconstructed organs group establishs a three-dimensional visualized operative environment. The 3D dose maps of the three-dimentional conformal localized dose planning are calculated with Monte Carlo method while the corresponding isodose lines and isodose surfaces are displayed in a stereo view. The real-time intra-operative navigation is based on an electromagnetic tracking system (ETS) and the fusion between MRI and ultrasound images. Applying Least Square Method, the coordinate registration between 3D models and patient is realized by the ETS which is calibrated by a laser tracker. The system is validated by working on eight patients with prostate cancer. The navigation has passed the precision measurement in the laboratory. Results: The traditional marching cubes (MC) method reconstructs one organ at one time and assembles them together. Compared to MC, presented multi-organs reconstruction method has superiorities in reserving the integrality and connectivity of reconstructed organs. The 3D conformal localized dose planning, realizing the 'exfoliation display' of different isodose surfaces, helps make sure the dose distribution has encompassed the nidus and avoid the injury of healthy tissues. During the navigation, surgeons could observe the coordinate of instruments real-timely employing the ETS. After the calibration, accuracy error of the needle position is less than 2.5mm according to the experiments. Conclusion: The speed and quality of 3D reconstruction, the efficiency in dose planning and accuracy in navigation all can be improved simultaneously.

60 APPLIED LIFE SCIENCES↗

SU-E-T-157: Evaluation and Comparison of Doses to Pelvic Lymph Nodes and to Point B with 3D Image Guided Treatment Planning for High Dose Brachytherapy for Treatment of Cervical Cancer

Purpose: To estimate and compare the doses received by the obturator, external and internal iliac lymph nodes and point Methods: CT-MR fused image sets of 15 patients obtained for each of 5 fractions of HDR brachytherapy using tandem and ring applicator, were used to generate treatment plans optimized to deliver a prescription dose to HRCTV-D90 and to minimize the doses to organs at risk (OARs). For each set of image, target volume (GTV, HRCTV) OARs (Bladder, Rectum, Sigmoid), and both left and right pelvic lymph nodes (obturator, external and internal iliac lymph nodes) were delineated. Dose-volume histograms (DVH) were generated for pelvic nodal groups (left and right obturator group, internal and external iliac chains) Per fraction DVH parameters used for dose comparison included dose to 100% volume (D100), and dose received by 2cc (D2cc), 1cc (D1cc) and 0.1 cc (D0.1cc) of nodal volume. Dose to point B was compared with each DVH parameter using 2 sided t-test. Pearson correlation were determined to examine relationship of point B dose with nodal DVH parameters. Results: FIGO clinical stage varied from 1B1 to IIIB. The median pretreatment tumor diameter measured on MRI was 4.5 cm (2.7– 6.4cm).The median dose to bilateral point B was 1.20 Gy ± 0.12 or 20% of the prescription dose. The correlation coefficients were all <0.60 for all nodal DVH parameters indicating low degree of correlation. Only 2 cc of obturator nodes was not significantly different from point B dose on t-test. Conclusion: Dose to point B does not adequately represent the dose to any specific pelvic nodal group. When using image guided 3D dose-volume optimized treatment nodal groups should be individually identified and delineated to obtain the doses received by pelvic nodes.

60 APPLIED LIFE SCIENCES↗

SU-E-T-245: MR Guided Focused Ultrasound Increased PARP Related Apoptosis On Prostate Cancer in Vivo

Purpose: Our previous study demonstrated that significant tumor growth delay was observed in the mice treated with pulsed high intensity focused ultrasound (pHIFU). The purpose of this study is to understand the cell killing mechanisms of pHIFU. Methods: Prostate cancer cells (LNCaP), were grown orthotopically in 17 nude mice. Tumor-bearing mice were treated using pHIFU with an acoustic power of 25W, pulse width 100msec and 300 pulses in one sonication under MR guidance. Mutiple sonications were used to cover the whole tumor volume. Temperature (less than 40 degree centigrade in the focal spot) was monitored using MR thermometry. Animals were euthanized at pre-determined time points (n=2) after treatment: 0 hours; 6 hrs; 24 hrs; 48 hrs; 4 days and 7 days. Two tumorbearing mice were used as control. Three tumor-bearing mice were treated with radiation (RT, 2 Gy) using 6 MV photon beams. RT treated mice were euthanized at 0 hr, 6 hrs and 24 hrs. The tumors were processed for immunohistochemical (IHC) staining for PARP (a surrogate of apoptosis). A multispectral imaging analysis system was used to quantify the expression of PARP staining. Cell apoptosis was calculated based on the PARP expression level, which is the intensity of the DAB reaction. Results: Our data showed that PARP related apoptosis peaked at 48 hrs and 7 days in pHIFU treated mice, which is comparable to that for the RT group at 24 hrs. The preliminary results from this study were consistent with our previous study on tumor growth delay using pHIFU. Conclusion: Our results demonstrated that non-thermal pHIFU increased apoptotic tumor cell death through the PARP related pathway. MR guided pHIFU may have a great potential as a safe, noninvasive treatment modality for cancer therapy. This treatment modality might be able to synergize with PARP inhibitors to achieve better result.

60 APPLIED LIFE SCIENCES↗

SU-E-T-255: Optimized Supine Craniospinal Irradiation with Image-Guided and Field Matched Beams

Purpose: Conventional craniospinal irradiation (CSI) challenges include dose inhomogeneity at field junctions and position uncertainty due to the field divergence, particular for the two spinal fields. Here we outline a new supine CSI technique to address these difficulties. Methods: Patient was simulated in supine position. The cranial fields had isocenter at C2/C3 vertebral and were matched with 1st spinal field. Their inferior border was chosen to avoid the shoulder, as well as chin from the 1st spine field. Their collimator angles were dependent on asymmetry jaw setting of the 1st spinal field. With couch rotation, the spinal field gantry angles were adjusted to ensure, the inferior border of 1st and superior border of 2nd spinal fields were perpendicular to the table top. The radio-opaque wire position for the spinal junction was located initially by the light field from an anterior setup beam, and was finalized by the portal imaging of the 1st spinal field. With reference to the spinal junction wire, the fields were matched by positioning the isocenter of the 2nd spinal field. A formula was derived to optimize supine CSI treatment planning, by utilizing the relationship among the Yjaw setting, the spinal field gantry angles, cranial field collimator angles, and the spinal field isocenters location. The plan was delivered with portal imaging alignment for the both cranial and spinal junctions. Results: Utilizing this technique with matching beams, and conventional technique such as feathering and forwarding planning, a homogenous dose distribution was achieved throughout the entire CSI treatment volume including the spinal junction. Placing the spinal junction wire visualized in both spinal portals, allows for precise determination and verification of the appropriate match line of the spine fields. Conclusion: This technique of optimization supine CSI achieved a homogenous dose distributions and patient localization accuracy with image-guided and matched beams.

60 APPLIED LIFE SCIENCES↗

SU-E-T-306: Study of the Reduction Technique for the Secondary Cancer Risk Due to Cone Beam CT in Image Guided Radiotherapy

Purpose: This study evaluated the effectiveness of a thin lead sheet based simple shielding method for imaging doses from cone beam computed tomography (CBCT) in image-guided radiotherapy (IGRT). Methods: The entire body, except for the region scanned by CBCT, was shielded by wrapping in a 2 mm lead sheet. Reduction of secondary doses from CBCT was measured using a radio-photoluminescence glass dosimeter (RPLGD) placed inside an anthropomorphic phantom and changes in secondary cancer risk due to the shielding effect were estimated using BEIR VII model. Results: Doses to out-of-field organs for head-and-neck, chest, and pelvis scans were decreased 15∼100 %, 23∼90 %, and 23∼98 %, respectively, and the average reductions in lifetime secondary cancer risk due to the 2 mm lead shielding were 1.61, 10.4, and 12.8 persons per 100,000, respectively. Conclusion: This study suggests that a simple thin lead sheet based shielding method results in a non-negligible reduction of secondary doses to out-of-field regions for CBCT.

60 APPLIED LIFE SCIENCES↗

SU-E-T-376: 3-D Commissioning for An Image-Guided Small Animal Micro- Irradiation Platform

Purpose: A 3-D radiochromic plastic dosimeter has been used to cross-test the isocentricity of a high resolution image-guided small animal microirradiation platform. In this platform, the mouse stage rotating for cone beam CT imaging is perpendicular to the gantry rotation for sub-millimeter radiation delivery. A 3-D dosimeter can be used to verify both imaging and irradiation coordinates. Methods: A 3-D dosimeter and optical CT scanner were used in this study. In the platform, both mouse stage and gantry can rotate 360° with rotation axis perpendicular to each other. Isocentricity and coincidence of mouse stage and gantry rotations were evaluated using star patterns. A 3-D dosimeter was placed on mouse stage with center at platform isocenter approximately. For CBCT isocentricity, with gantry moved to 90°, the mouse stage rotated horizontally while the x-ray was delivered to the dosimeter at certain angles. For irradiation isocentricity, the gantry rotated 360° to deliver beams to the dosimeter at certain angles for star patterns. The uncertainties and agreement of both CBCT and irradiation isocenters can be determined from the star patterns. Both procedures were repeated 3 times using 3 dosimeters to determine short-term reproducibility. Finally, dosimeters were scanned using optical CT scanner to obtain the results. Results: The gantry isocentricity is 0.9 ± 0.1 mm and mouse stage rotation isocentricity is about 0.91 ± 0.11 mm. Agreement between the measured isocenters of irradiation and imaging coordinates was determined. The short-term reproducibility test yielded 0.5 ± 0.1 mm between the imaging isocenter and the irradiation isocenter, with a maximum displacement of 0.7 ± 0.1 mm. Conclusion: The 3-D dosimeter can be very useful in precise verification of targeting for a small animal irradiation research. In addition, a single 3-D dosimeter can provide information in both geometric and dosimetric uncertainty, which is crucial for translational studies.

60 APPLIED LIFE SCIENCES↗

SU-E-T-453: A Novel Daily QA System for Robotic Image Guided Radiosurgery with Variable Aperture Collimator

Purpose: A novel end-to-end system using a CCD camera and a scintillator based phantom that is capable of measuring the beam-by-beam delivery accuracy of Robotic Radiosurgery has been developed and reported in our previous work. This work investigates its application to end-to-end type daily QA for Robotic Radiosurgery (Cyberknife) with Variable Aperture Collimator (Iris). Methods: The phantom was first scanned with a CT scanner at 0.625 slice thickness and exported to the Cyberknife Muliplan (v4.6) treatment planning system. An isocentric treatment plan was created consisting of ten beams of 25 Monitor Units each using Iris apertures of 7.5, 10, 15, 20, and 25 mm. The plan was delivered six times in two days on the Cyberknife G4 system with fiducial tracking on the four metal fiducials embedded in phantom with re-positioning between the measurements. The beam vectors (X, Y, Z) are measured and compared with the plan from the machine delivery file (XML file). The Iris apertures (FWHM) were measured from the beam flux map and compared with the commissioning data. Results: The average beam positioning accuracies of the six deliveries are 0.71 ± 0.40 mm, 0.72 ± 0.44 mm, 0.74 ± 0.42 mm, 0.70 ± 0.40 mm, 0.79 ± 0.44 mm and 0.69 ± 0.41 mm respectively. Radiation beam width (FWHM) variations are within ±0.05 mm, and they agree with the commissioning data within 0.22 mm. The delivery time for the plan is about 7 minutes and the results are given instantly. Conclusion: The experimental results agree with stated sub-millimeter delivery accuracy of Cyberknife system. Beam FWHM variations comply with the 0.2 mm accuracy of the Iris collimator at SAD. The XRV-100 system has proven to be a powerful tool in performing end-to-end type tests for Robotic Image Guided Radiosurgery Daily QA.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-E-T-500: Dose Escalation Strategy for Lung Cancer Patients Using a Biologically- Guided Target Definition

Purpose: Dose escalation strategy for lung cancer patients can lead to late symptoms such as pneumonitis and cardiac injury. We propose a strategy to increase radiation dose for improving local tumor control while simultaneously striving to minimize the injury of organs at risk (OAR). Our strategy is based on defining a small, biologically-guided target volume for receiving additional radiation dose. Methods: 106 patients with lung cancer treated with radiotherapy were selected for patients diagnosed with stage II and III disease. Previous research has shown that 50% of the maximum SUV threshold in FDG-PET imaging is appropriate for delineation of the most aggressive part of a tumor. After PET- and CT-derived targets were contoured, an IMRT treatment plan was designed to deliver 60 Gy to the GTV as delineated on a 4D CT (Plan 1). A second plan was designed with additional dose of 18 Gy to the PET-derived volume (Plan 2). A composite plan was generated by the addition of Plan 1 and Plan 2. Results: Plan 1 was compared to the composite plan and increases in OAR dose were assessed. For seven patients on average, lung V5 was increased by 1.4% and V20 by 4.2% for ipsilateral lung and by 13.5% and 7% for contralateral lung. For total lung, V5 and V20 were increased by 4.5% and 4.8% respectively. Mean lung dose was increased by 9.7% for the total lung. The maximum dose to the spinal cord increased by 16% on average. For the heart, V20 increased by 4.2% and V40 by 5.2%. Conclusion: It seems feasible that an additional 18 Gy of radiation dose can be delivered to FDG PET-derived subvolume of the CT-based GTV of the primary tumor without significant increase in total dose to the critical organs such as lungs, spinal cord and heart.

60 APPLIED LIFE SCIENCES↗

TU-F-BRF-09: Feasibility Study of Spatial and Temporal Fractionation Using a Table-Top Image-Guided MRT System

Purpose: Microbeam radiation therapy is a promising experimental radiotherapy method for cancer treatment. Synchrotron studies have shown that MRT can preferentially ablate tumors while mostly preserving the surrounding normal tissues. Our purpose is to develop a compact microbeam irradiator that can be accessible for laboratory research on MRT's therapeutic mechanism, with the ultimate goal of translating this technique for clinical applications. Utilizing a carbon nanotube field emission X-ray source array, our lab has developed a first of its kind table-top microbeam irradiator. The purpose of this study is to investigate the possibility of increasing the total dose to the tumor region using spatial and temporal fractionation. Methods: MRI and CT scans were acquired and registered to locate both the target region and the landmark in the treatment coordinate space. A crossbeam treatment configuration was planned to achieve the desired dosage and tumor coverage. One array of 300 μm-thick microbeams was delivered to the target on the first day. Another array, perpendicular to the first one, was delivered the next day targeting the same region. The feasibility of this image-guided MRT procedure was evaluated using a phantom. Gafchromic EBT2 film was employed to characterize the targeting accuracy and the delivered tumor dose. Results: The preliminary phantom study showed that a cross-pattern of microbeam arrays was successfully delivered with sub-millimeter targeting accuracy. The delivered dose in the target region was doubled. The fraction of tumor volume that received the peak dose was increased. Conclusion: Temporal fractionation of crosspatterned microbeams can be delivered to the same target region. The feasibility of this procedure was investigated with a phantom study. Once implemented, this procedure has the potential of increasing the tumor control effect of MRT and of achieving better normal tissue recovery.

60 APPLIED LIFE SCIENCES↗

WE-A-17A-05: Differences in Applicator Configuration and Dwell Loading Between Standard and Image-Guided Tandem and Ring (T and R) HDR Brachytherapy

Purpose: To investigate differences in: (i) relative location of the tandem and the ring compared to a rigid standard applicator model; and (ii) relative loading and changes in loading pattern between standard and image-guided planning. Methods: All T and R insertions performed in 2013 in our institution under CT- or MR-guidance were analyzed. Standard plans were generated using library applicator models with a fixed relationship between ring and tandem, standardized uniform dwell loading and normalization to point A. The graphic plans and the associated standard-plan dwell configurations were compared: the rings were rigidly registered, and the residual tandem shift, rotation and maximum distance between plan tandem dwell and corresponding model tandem dwell were calculated. The normalization ratio (NR = the ratio of graphic versus standard-plan total reference air kerma [TRAK]), the general loading difference (GLD = the difference between graphic and standard ratios of the tandem versus the ring TRAK), and the percent standard deviation (SD% = SD/mean) of the tandem and the ring TRAK for the graphic plan (all standard-plans SD% = 0) were calculated. Results: 71 T and R were analyzed. Residual tandem shift, rotation and maximum corresponding dwell distance were 1.2±0.8mm (0.4±0.4mm lateral, 0.9±0.8mm craniocaudal, 0.4±0.3mm anterior-posterior), 2.3±1.9deg and 3.4±2.3mm. NR was 0.86±0.11 indicating a lower overall loading of the graphic compared to the standard plans. GLD was -0.12±0.16 indicating a modest increased ring loading relative to the tandem in the graphic plans. SD% was 2.1±1.6% for tandem and 2.8±1.9% for ring, indicating small deviations from uniform loading. Conclusion: Variability in the relative locations of the tandem and the ring necessitates the independent registration of each component model for accurate digitization. Our clinical experience suggests that graphically planned T and R results on average in a lower total dose to the patient, and in an increased relative loading of the ring component.

60 APPLIED LIFE SCIENCES↗

WE-D-9A-02: Automated Landmark-Guided CT to Cone-Beam CT Deformable Image Registration

Purpose: The anatomical changes that occur between the simulation CT and daily cone-beam CT (CBCT) are investigated using an automated landmark-guided deformable image registration (LDIR) algorithm with simultaneous intensity correction. LDIR was designed to be accurate in the presence of tissue intensity mismatch and heavy noise contamination. Method: An auto-landmark generation algorithm was used in conjunction with a local small volume (LSV) gradient matching search engine to map corresponding landmarks between the CBCT and planning CT. The LSVs offsets were used to perform an initial deformation, generate landmarks, and correct local intensity mismatch. The landmarks act as stabilizing controlpoints in the Demons objective function. The accuracy of the LDIR algorithm was evaluated on one synthetic case with ground truth and data of ten head and neck cancer patients. The deformation vector field (DVF) accuracy was accessed using a synthetic case. The Root mean square error of the 3D canny edge (RMSECE), mutual information (MI), and feature similarity index metric (FSIM) were used to access the accuracy of LDIR on the patient data. The quality of the corresponding deformed contours was verified by an attending physician. Results: The resulting 90 percentile DVF error for the synthetic case was within 5.63mm for the original demons algorithm, 2.84mm for intensity correction alone, 2.45mm using controlpoints without intensity correction, and 1.48 mm for the LDIR algorithm. For the five patients the mean RMSECE of the original CT, Demons deformed CT, intensity corrected Demons CT, control-point stabilized deformed CT, and LDIR CT was 0.24, 0.26, 0.20, 0.20, and 0.16 respectively. Conclusion: LDIR is accurate in the presence of multimodal intensity mismatch and CBCT noise contamination. Since LDIR is GPU based it can be implemented with minimal additional strain on clinical resources. This project has been supported by a CPRIT individual investigator award RP11032.

60 APPLIED LIFE SCIENCES↗

WE-G-17A-01: Improving Tracking Image Spatial Resolution for Onboard MR Image Guided Radiation Therapy Using the WHISKEE Technique

Purpose: To evaluate the feasibility of using the weighted hybrid iterative spiral k-space encoded estimation (WHISKEE) technique to improve spatial resolution of tracking images for onboard MR image guided radiation therapy (MR-IGRT). Methods: MR tracking images of abdomen and pelvis had been acquired from healthy volunteers using the ViewRay onboard MRIGRT system (ViewRay Inc. Oakwood Village, OH) at a spatial resolution of 2.0mm*2.0mm*5.0mm. The tracking MR images were acquired using the TrueFISP sequence. The temporal resolution had to be traded off to 2 frames per second (FPS) to achieve the 2.0mm in-plane spatial resolution. All MR images were imported into the MATLAB software. K-space data were synthesized through the Fourier Transform of the MR images. A mask was created to selected k-space points that corresponded to the under-sampled spiral k-space trajectory with an acceleration (or undersampling) factor of 3. The mask was applied to the fully sampled k-space data to synthesize the undersampled k-space data. The WHISKEE method was applied to the synthesized undersampled k-space data to reconstructed tracking MR images at 6 FPS. As a comparison, the undersampled k-space data were also reconstructed using the zero-padding technique. The reconstructed images were compared to the original image. The relatively reconstruction error was evaluated using the percentage of the norm of the differential image over the norm of the original image. Results: Compared to the zero-padding technique, the WHISKEE method was able to reconstruct MR images with better image quality. It significantly reduced the relative reconstruction error from 39.5% to 3.1% for the pelvis image and from 41.5% to 4.6% for the abdomen image at an acceleration factor of 3. Conclusion: We demonstrated that it was possible to use the WHISKEE method to expedite MR image acquisition for onboard MR-IGRT systems to achieve good spatial and temporal resolutions simultaneously. Y. Hu and O. green receive travel reimbursement from ViewRay. S. Mutic has consulting and research agreements with ViewRay. Q. Zeng, R. Nana, J.L. Patrick, S. Shvartsman and J.F. Dempsey are ViewRay employees.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

TH-E-17A-11: Tracking Tumor’s Boundary in MV Image Sequences for Image-Guided Radiation Therapy

Purpose: To develop a level set method (LSM)-based algorithm to track the tumor’s boundary in MV image sequences for image-guided radiation therapy (IGRT). Method: Four MV image sequences, each of which consists of 100 frames at frame rate of 7.5 Hz, are acquired by using the electronic portal imaging device (EPID) (Varian Medical Systems, Pal Alto, USA) during the treatment of lung cancer. In the first frame of each sequence, we roughly delineate an initial contour of the target tumor by hand. Using a LSM-based algorithm, the initial contour can automatically shape itself to fit the tumor, and eventually detect the tumor’s boundary. We then employ the tumor’s boundary obtained from the previous frame as the initial contour in the subsequent frame, so that the LSM-based method can drive this initial contour to the tumor’s boundary quickly and fulfill a tracking task. Results: The proposed method has been evaluated on four MV image sequences. The mean tracking errors were 0.23, 0.29, 0.37, and 1.18 mm, and their corresponding standard deviations were 0.97, 0.75, 1.2 and 1.48 mm, respectively. Conclusion: Compared with conventional tumor tracking techniques, the proposed system is capable not only of tracking the tumor’s position, but also of detecting the tumor’s boundary varying with the respiration during the treatment. Considering current radiation therapy technique, for example, dynamical multi-leaf collimator (DMLC) has been widely applied in clinical treatment, this study indicates the potential for significant accuracy improvement in radiation therapy. This work was partially supported by a research grant from Varian Medical Systems (Palo Alto, California)

60 APPLIED LIFE SCIENCES↗

Sampling strategies for subsampled segmented EPI PRF thermometry in MR guided high intensity focused ultrasound

Purpose: To investigate k-space subsampling strategies to achieve fast, large field-of-view (FOV) temperature monitoring using segmented echo planar imaging (EPI) proton resonance frequency shift thermometry for MR guided high intensity focused ultrasound (MRgHIFU) applications. Methods: Five different k-space sampling approaches were investigated, varying sample spacing (equally vs nonequally spaced within the echo train), sampling density (variable sampling density in zero, one, and two dimensions), and utilizing sequential or centric sampling. Three of the schemes utilized sequential sampling with the sampling density varied in zero, one, and two dimensions, to investigate sampling the k-space center more frequently. Two of the schemes utilized centric sampling to acquire the k-space center with a longer echo time for improved phase measurements, and vary the sampling density in zero and two dimensions, respectively. Phantom experiments and a theoretical point spread function analysis were performed to investigate their performance. Variable density sampling in zero and two dimensions was also implemented in a non-EPI GRE pulse sequence for comparison. All subsampled data were reconstructed with a previously described temporally constrained reconstruction (TCR) algorithm. Results: The accuracy of each sampling strategy in measuring the temperature rise in the HIFU focal spot was measured in terms of the root-mean-square-error (RMSE) compared to fully sampled “truth.” For the schemes utilizing sequential sampling, the accuracy was found to improve with the dimensionality of the variable density sampling, giving values of 0.65 °C, 0.49 °C, and 0.35 °C for density variation in zero, one, and two dimensions, respectively. The schemes utilizing centric sampling were found to underestimate the temperature rise, with RMSE values of 1.05 °C and 1.31 °C, for variable density sampling in zero and two dimensions, respectively. Similar subsampling schemes with variable density sampling implemented in zero and two dimensions in a non-EPI GRE pulse sequence both resulted in accurate temperature measurements (RMSE of 0.70 °C and 0.63 °C, respectively). With sequential sampling in the described EPI implementation, temperature monitoring over a 192 × 144 × 135 mm{sup 3} FOV with a temporal resolution of 3.6 s was achieved, while keeping the RMSE compared to fully sampled “truth” below 0.35 °C. Conclusions: When segmented EPI readouts are used in conjunction with k-space subsampling for MR thermometry applications, sampling schemes with sequential sampling, with or without variable density sampling, obtain accurate phase and temperature measurements when using a TCR reconstruction algorithm. Improved temperature measurement accuracy can be achieved with variable density sampling. Centric sampling leads to phase bias, resulting in temperature underestimations.

42 ENGINEERING↗

Towards the clinical implementation of iterative low-dose cone-beam CT reconstruction in image-guided radiation therapy: Cone/ring artifact correction and multiple GPU implementation

Purpose: Compressed sensing (CS)-based iterative reconstruction (IR) techniques are able to reconstruct cone-beam CT (CBCT) images from undersampled noisy data, allowing for imaging dose reduction. However, there are a few practical concerns preventing the clinical implementation of these techniques. On the image quality side, data truncation along the superior–inferior direction under the cone-beam geometry produces severe cone artifacts in the reconstructed images. Ring artifacts are also seen in the half-fan scan mode. On the reconstruction efficiency side, the long computation time hinders clinical use in image-guided radiation therapy (IGRT). Methods: Image quality improvement methods are proposed to mitigate the cone and ring image artifacts in IR. The basic idea is to use weighting factors in the IR data fidelity term to improve projection data consistency with the reconstructed volume. In order to improve the computational efficiency, a multiple graphics processing units (GPUs)-based CS-IR system was developed. The parallelization scheme, detailed analyses of computation time at each step, their relationship with image resolution, and the acceleration factors were studied. The whole system was evaluated in various phantom and patient cases. Results: Ring artifacts can be mitigated by properly designing a weighting factor as a function of the spatial location on the detector. As for the cone artifact, without applying a correction method, it contaminated 13 out of 80 slices in a head-neck case (full-fan). Contamination was even more severe in a pelvis case under half-fan mode, where 36 out of 80 slices were affected, leading to poorer soft tissue delineation and reduced superior–inferior coverage. The proposed method effectively corrects those contaminated slices with mean intensity differences compared to FDK results decreasing from ∼497 and ∼293 HU to ∼39 and ∼27 HU for the full-fan and half-fan cases, respectively. In terms of efficiency boost, an overall 3.1 × speedup factor has been achieved with four GPU cards compared to a single GPU-based reconstruction. The total computation time is ∼30 s for typical clinical cases. Conclusions: The authors have developed a low-dose CBCT IR system for IGRT. By incorporating data consistency-based weighting factors in the IR model, cone/ring artifacts can be mitigated. A boost in computational efficiency is achieved by multi-GPU implementation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A simulation technique for 3D MR-guided acoustic radiation force imaging

Purpose: In magnetic resonance-guided focused ultrasound (MRgFUS) therapies, the in situ characterization of the focal spot location and quality is critical. MR acoustic radiation force imaging (MR-ARFI) is a technique that measures the tissue displacement caused by the radiation force exerted by the ultrasound beam. This work presents a new technique to model the displacements caused by the radiation force of an ultrasound beam in a homogeneous tissue model. Methods: When a steady-state point-source force acts internally in an infinite homogeneous medium, the displacement of the material in all directions is given by the Somigliana elastostatic tensor. The radiation force field, which is caused by absorption and reflection of the incident ultrasound intensity pattern, will be spatially distributed, and the tensor formulation takes the form of a convolution of a 3D Green’s function with the force field. The dynamic accumulation of MR phase during the ultrasound pulse can be theoretically accounted for through a time-of-arrival weighting of the Green’s function. This theoretical model was evaluated experimentally in gelatin phantoms of varied stiffness (125-, 175-, and 250-bloom). The acoustic and mechanical properties of the phantoms used as parameters of the model were measured using independent techniques. Displacements at focal depths of 30- and 45-mm in the phantoms were measured by a 3D spin echo MR-ARFI segmented-EPI sequence. Results: The simulated displacements agreed with the MR-ARFI measured displacements for all bloom values and focal depths with a normalized RMS difference of 0.055 (range 0.028–0.12). The displacement magnitude decreased and the displacement pattern broadened with increased bloom value for both focal depths, as predicted by the theory. Conclusions: A new technique that models the displacements caused by the radiation force of an ultrasound beam in a homogeneous tissue model theory has been rigorously validated through comparison with experimentally obtained 3D displacement data in homogeneous gelatin phantoms using a 3D MR-ARFI sequence. The agreement of the experimentally measured and simulated results demonstrates the potential to use MR-ARFI displacement data in MRgFUS therapies.

60 APPLIED LIFE SCIENCES↗