Skin dose investigations on a 0.5 T parallel rotating biplanar linac-MR using Monte Carlo simulations and measurements.

BACKGROUND: The Alberta rotating biplanar linac-MR has a 0.5 T magnetic field parallel to the beamline. When developing a new linac-MR system, interactions of charged particles with the magnetic field necessitate careful consideration of skin dose and tissue interface effects. PURPOSE: To investigate the effect of the magnetic field on skin dose using measurements and Monte Carlo (MC) simulations. METHODS: We develop an MC model of our linac-MR, which we validate by comparison with ion chamber measurements in a water tank. Additionally, MC simulation results are compared with radiochromic film surface dose measurements on solid water. Variations in surface dose as a function of field size are measured using a parallel plate ion chamber in solid water. Using an anthropomorphic computational phantom with a 2 mm-thick skin layer, we investigate dose distributions resulting from three beam arrangements. Magnetic field on and off scenarios are considered for all measurements and simulations. RESULTS: For a 20 × 20 cm2 field size, D 0.2 c c ${D_{0.2cc}}$ (the minimum dose to the hottest contiguous 0.2 cc volume) for the top 2 mm of a simple water phantom is 72% when the magnetic field is on, compared to 34% with magnetic field off (values are normalized to the central axis dose maximum). Parallel plate ion chamber measurements demonstrate that the relative increase in surface dose due to the magnetic field decreases with increasing field size. For the anthropomorphic phantom, D ∼ 0.2 c c ${D_{ \sim 0.2cc}}$ (minimum skin dose in the hottest 1 × 1 × 1 cm3 cube) shows relative increases of 20%-28% when the magnetic field is on compared to when it is off. With magnetic field off, skin D ∼ 0.2 c c ${D_{ \sim 0.2cc}}$ is 71%, 56%, and 21% for medial-lateral tangents, anterior-posterior beams, and a five-field arrangement, respectively. For magnetic field on, the corresponding skin D ∼ 0.2 c c ${D_{ \sim 0.2cc}}$ values are 91%, 67%, and 25%. CONCLUSIONS: Using a validated MC model of our linac-MR, surface doses are calculated in various scenarios. MC-calculated skin dose varies depending on field sizes, obliquity, and the number of beams. In general, the parallel linac-MR arrangement results in skin dose enhancement due to charged particles spiraling along magnetic field lines, which impedes lateral motion away from the central axis. Nonetheless, considering the results presented herein, treatment plans can be designed to minimize skin dose by, for example, avoiding oblique beams and using a larger number of fields.

Clinical reference dosimetry for the 0.5 T inline rotating biplanar Linac-MR.

BACKGROUND: The world's first clinical 0.5 T inline rotating biplanar Linac-MR system is commissioned for clinical use. For reference dosimetry, unique features to device, including an SAD = 120 cm, bore clearance of 60 cm × 110 cm, as well as 0.5 T inline magnetic field, provide some challenges to applying a standard dosimetry protocol (i.e., TG-51). PURPOSE: In this work, we propose a simple and practical clinical reference dosimetry protocol for the 0.5T biplanar Linac-MR and validated its results. METHODS: Our dosimetry protocol for this system is as follows: tissue phantom ratios at 20 and 10 cm are first measured and converted into %dd10x beam quality specifier using equations provided and Kalach and Rogers. The converted %dd10x is used to determine the ion chamber correction factor, using the equations in the TG-51 addendum for the Exradin A12 farmer chamber used, which is cross-calibrated with one calibrated at a standards laboratory. For a 0.5 T parallel field, magnetic field effect on chamber response is assumed to have no effect and is not explicitly corrected for. Once the ion chamber correction factor for a non-standard SAD (kQ,msr) is determined, TG-51 is performed to obtain dose at a depth of 10 cm at SAD = 120 cm. The dosimetry protocol is repeated with the magnetic field ramped down. To validate our dosimetry protocol, Monte Carlo (EGSnrc) simulations are performed to confirm the determined kQ,msr values. MC Simulations and magnetic Field On versus Field Off measurements are performed to confirm that the magnetic field has no effect. To validate our overall dosimetry protocol, external dose audits, based on optical simulated luminescent dosimeters, thermal luminescent dosimeters, and alanine dosimeters are performed on the 0.5 T Linac-MR system. RESULTS: Our EGSnrc results confirm our protocol-determined kQ,msr values, as well as our assumptions about magnetic field effects (kB = 1) within statistical uncertainty for the A-12 chamber. Our external dosimetry procedures also validated our overall dosimetry protocol for the 0.5 T biplanar Linac-MR hybrid. Ramping down the magnetic field has resulted in a dosimetric difference of 0.1%, well within experimental uncertainty. CONCLUSION: With the 0.5 T parallel magnetic field having minimal effect on the ion chamber response, a TPR20,10 approach to determine beam quality provides an accurate method to perform clinical dosimetry for the 0.5 T biplanar Linac-MR.

Time domain principal component analysis for rapid, real-time 2D MRI reconstruction from undersampled data.

PURPOSE: A rapid real-time 2D accelerated method was developed for magnetic resonance imaging (MRI) using principal component analysis (PCA) in the temporal domain. This method employs a moving window of previous dynamic frames to reconstruct the current, real-time frame within this window. This technique could be particularly useful in real-time tracking applications such as in MR-guided radiotherapy, where low latency real-time reconstructions are essential. METHODS: The method was tested retrospectively on 15 fully-sampled data sets of lung patient data acquired on a 3T Philips Achieva system. High frequency data are incoherently undersampled, while the central low-frequency data are always acquired to characterize the temporal fluctuations through PCA. The undersampling pattern is derived in such a way that all of k-space is acquired within a pre-determined number of frames. The missing data in the current frame are then filled in by fitting the temporal characterizations to the acquired undersampled data, using a pre-determined number of PCs. A subset of six patients was used to test the contour ability of the images. Various accelerations between 3x and 8x were tested along with the optimal number of PCs for fitting. A comparison was also performed with previous work from our group proposed by Dietz et al. as well as with a standard low resolution acquisition. In order to determine how the method would perform at lower signal to noise ratio (SNR), noise levels of 2×, 4×, and 6× were added to the 3T data. Metrics such as normalised mean square error and Dice coefficient were used to measure the reconstruction image quality and contour ability. RESULTS: The proposed method demonstrated good temporal robustness as consistent metrics were detected for the duration of the imaging session. It was found that the optimal number of PCs for temporal fitting was dependent on the acceleration rate. For the data tested, five PCs were found to be optimal at the acceleration rates of 3× and 4×. This number decreases to three at accelerations of 5× and 6× and further decreases to two at an acceleration rate of 8×, likely due to greater instability with fewer acquired data points. The use of too many PCs for fitting increased the chances of noisy reconstruction which affected contourability. CONCLUSIONS: The proposed 2D real-time MR acceleration method demonstrated greater robustness in the metrics over time when compared with previous real-time PCA methods using metrics such as normalised mean squared error, peak SNR and structural similarity up to an acceleration of 8x. Improved temporal robustness of image structure contourability and accurate definition was also demonstrated using several metrics including the Dice coefficient. Reconstruction of raw acquired data can be performed at approximately 50 ms per frame using an Intel core i5 CPU. The method has the advantage of being very flexible in terms of hardware requirements as it can operate successfully on a single coil channel and does not require specialized computing power to implement in real-time.

Myelin Water Imaging Demonstrates Lower Brain Myelination in Children and Adolescents With Poor Reading Ability.

Magnetic resonance imaging (MRI) provides a means to non-invasively investigate the neurological links with dyslexia, a learning disability that affects one's ability to read. Most previous brain MRI studies of dyslexia and reading skill have used structural or diffusion imaging to reveal regional brain abnormalities. However, volumetric and diffusion MRI lack specificity in their interpretation at the microstructural level. Myelin is a critical neural component for brain function and plasticity, and as such, deficits in myelin may impact reading ability. MRI can estimate myelin using myelin water fraction (MWF) imaging, which is based on evaluation of the proportion of short T2 myelin-associated water from multi-exponential T2 relaxation analysis, but has not yet been applied to the study of reading or dyslexia. In this study, MWF MRI, intelligence, and reading assessments were acquired in 20 participants aged 10-18 years with a wide range of reading ability to investigate the relationship between reading ability and myelination. Group comparisons showed markedly lower MWF by 16-69% in poor readers relative to good readers in the left and right thalamus, as well as the left posterior limb of the internal capsule, left/right anterior limb of the internal capsule, left/right centrum semiovale, and splenium of the corpus callosum. MWF over the entire group also correlated positively with three different reading scores in the bilateral thalamus as well as white matter, including the splenium of the corpus callosum, left posterior limb of the internal capsule, left anterior limb of the internal capsule, and left centrum semiovale. MWF imaging from T2 relaxation suggests that myelination, particularly in the bilateral thalamus, splenium, and left hemisphere white matter, plays a role in reading abilities. Myelin water imaging thus provides a potentially valuable in vivo imaging tool for the study of dyslexia and its remediation.

Single patient convolutional neural networks for real-time MR reconstruction: coherent low-resolution versus incoherent undersampling.

Accelerated MRI involves undersampling k-space, creating unwanted artifacts when reconstructing the data. While the strategy of incoherent k-space acquisition is proven for techniques such as compressed sensing, it may not be optimal for all techniques. This study compares the use of coherent low-resolution (coherent-LR) and incoherent undersampling phase-encoding for real-time 3D CNN image reconstruction. Data were acquired with our 3 T Philips Achieva system. A retrospective analysis was performed on six non-small cell lung cancer patients who received dynamic acquisitions consisting of 650 free breathing images using a bSSFP sequence. We retrospectively undersampled the data by 5x and 10x acceleration using the two phase-encoding schemes. A quantitative analysis was conducted evaluating the tumor segmentations from the CNN reconstructed data using the Dice coefficient (DC) and centroid displacement. The reconstruction noise was evaluated using the structural similarity index (SSIM). Furthermore, we qualitatively investigated the CNN reconstruction using prospectively undersampled data, where the fully sampled training data set is acquired separately from the accelerated undersampled data. The patient averaged DC, centroid displacement, and SSIM for the tumor segmentation at 5x and 10x was superior using coherent low-resolution undersampling. Furthermore, the patient-specific CNN can be trained in under 6 h and the reconstruction time was 54 ms per image. Both the incoherent and coherent-LR prospective CNN reconstructions yielded qualitatively acceptable images; however, the coherent-LR reconstruction appeared superior to the incoherent reconstruction. We have demonstrated that coherent-LR undersampling for real-time CNN image reconstruction performs quantitatively better for the retrospective case of lung tumor segmentation, and qualitatively better for the prospective case. The tumor segmentation mean DC increased for all six patients at 5x acceleration and the temporal (dynamic) variance of the segmentation was reduced. The reconstruction speed achieved for our current implementation was 54 ms, providing an acceptable frame rate for real-time on-the-fly MR imaging.

Single patient convolutional neural networks for real-time MR reconstruction: a proof of concept application in lung tumor segmentation for adaptive radiotherapy.

Investigate 3D (spatial and temporal) convolutional neural networks (CNNs) for real-time on-the-fly magnetic resonance imaging (MRI) reconstruction. In particular, we investigated the applicability of training CNNs on a patient-by-patient basis for the purpose of lung tumor segmentation. Data were acquired with our 3 T Philips Achieva system. A retrospective analysis was performed on six non-small cell lung cancer patients who received fully sampled dynamic acquisitions consisting of 650 free breathing images using a bSSFP sequence. We retrospectively undersampled the six patient's data by 5× and 10× acceleration. The retrospective data was used to quantitatively compare the CNN reconstruction to gold truth data via the Dice coefficient (DC) and centroid displacement to compare the tumor segmentations. Reconstruction noise was investigated using the normalized mean square error (NMSE). We further validated the technique using prospectively undersampled data from a volunteer and motion phantom. The retrospectively undersampled data at 5× and 10× acceleration was reconstructed using patient specific trained CNNs. The patient average DCs for the tumor segmentation at 5× and 10× acceleration were 0.94 and 0.92, respectively. These DC values are greater than the inter- and intra-observer segmentations acquired by radiation oncologist experts as reported in a previous study of ours. Furthermore, the patient specific CNN can be trained in under 6 h and the reconstruction time was 65 ms per image. The prospectively undersampled CNN reconstruction data yielded qualitatively acceptable images. We have shown that 3D CNNs can be used for real-time on-the-fly dynamic image reconstruction utilizing both spatial and temporal data in this proof of concept study. We evaluated the technique using six retrospectively undersampled lung cancer patient data sets, as well as prospectively undersampled data acquired from a volunteer and motion phantom. The reconstruction speed achieved for our current implementation was 65 ms per image.

A Parametric Computational Model of the Action Potential of Pacemaker Cells.

OBJECTIVE: A flexible, efficient, and verifiable pacemaker cell model is essential to the design of real-time virtual hearts that can be used for closed-loop validation of cardiac devices. A new parametric model of pacemaker action potential is developed to address this need. METHODS: The action potential phases are modeled using hybrid automaton with one piecewise-linear continuous variable. The model can capture rate-dependent dynamics, such as action potential duration restitution, conduction velocity restitution, and overdrive suppression by incorporating nonlinear update functions. Simulated dynamics of the model compared well with previous models and clinical data. CONCLUSION: The results show that the parametric model can reproduce the electrophysiological dynamics of a variety of pacemaker cells, such as sinoatrial node, atrioventricular node, and the His-Purkinje system, under varying cardiac conditions. SIGNIFICANCE: This is an important contribution toward closed-loop validation of cardiac devices using real-time heart models.

Evaluating performance of a user-trained MR lung tumor autocontouring algorithm in the context of intra- and interobserver variations.

PURPOSE: Real-time tracking of lung tumors using magnetic resonance imaging (MRI) has been proposed as a potential strategy to mitigate the ill-effects of breathing motion in radiation therapy. Several autocontouring methods have been evaluated against a "gold standard" of a single human expert user. However, contours drawn by experts have inherent intra- and interobserver variations. In this study, we aim to evaluate our user-trained autocontouring algorithm with manually drawn contours from multiple expert users, and to contextualize the accuracy of these autocontours within intra- and interobserver variations. METHODS: Six nonsmall cell lung cancer patients were recruited, with institutional ethics approval. Patients were imaged with a clinical 3 T Philips MR scanner using a dynamic 2D balanced SSFP sequence under free breathing. Three radiation oncology experts, each in two separate sessions, contoured 130 dynamic images for each patient. For autocontouring, the first 30 images were used for algorithm training, and the remaining 100 images were autocontoured and evaluated. Autocontours were compared against manual contours in terms of Dice's coefficient (DC) and Hausdorff distances (dH ). Intra- and interobserver variations of the manual contours were also evaluated. RESULTS: When compared with the manual contours of the expert user who trained it, the algorithm generates autocontours whose evaluation metrics (same session: DC = 0.90(0.03), dH  = 3.8(1.6) mm; different session DC = 0.88(0.04), dH  = 4.3(1.5) mm) are similar to or better than intraobserver variations (DC = 0.88(0.04), and dH  = 4.3(1.7) mm) between two sessions. The algorithm's autocontours are also compared to the manual contours from different expert users with evaluation metrics (DC = 0.87(0.04), dH  = 4.8(1.7) mm) similar to interobserver variations (DC = 0.87(0.04), dH  = 4.7(1.6) mm). CONCLUSIONS: Our autocontouring algorithm delineates tumor contours (<20 ms per contour), in dynamic MRI of lung, that are comparable to multiple human experts (several seconds per contour), but at a much faster speed. At the same time, the agreement between autocontours and manual contours is comparable to the intra- and interobserver variations. This algorithm may be a key component of the real time tumor tracking workflow for our hybrid Linac-MR device in the future.

Real-time dynamic MR image reconstruction using compressed sensing and principal component analysis (CS-PCA): Demonstration in lung tumor tracking.

PURPOSE: This work presents a real-time dynamic image reconstruction technique, which combines compressed sensing and principal component analysis (CS-PCA), to achieve real-time adaptive radiotherapy with the use of a linac-magnetic resonance imaging system. METHODS: Six retrospective fully sampled dynamic data sets of patients diagnosed with non-small-cell lung cancer were used to investigate the CS-PCA algorithm. Using a database of fully sampled k-space, principal components (PC's) were calculated to aid in the reconstruction of undersampled images. Missing k-space data were calculated by projecting the current undersampled k-space data onto the PC's to generate the corresponding PC weights. The weighted PC's were summed together, and the missing k-space was iteratively updated. To gain insight into how the reconstruction might proceed at lower fields, 6× noise was added to the 3T data to investigate how the algorithm handles noisy data. Acceleration factors ranging from 2 to 10× were investigated using CS-PCA and Split Bregman CS for comparison. Metrics to determine the reconstruction quality included the normalized mean square error (NMSE), as well as the dice coefficients (DC) and centroid displacement of the tumor segmentations. RESULTS: Our results demonstrate that CS-PCA performed superior than CS alone. The CS-PCA patient averaged DC for 3T and 6× noise added data remained above 0.9 for acceleration factors up to 10×. The patient averaged NMSE gradually increased with increasing acceleration; however, it remained below 0.06 up to an acceleration factor of 10× for both 3T and 6× noise added data. The CS-PCA reconstruction speed ranged from 5 to 20 ms (Intel i7-4710HQ CPU @ 2.5 GHz), depending on the chosen parameters. CONCLUSIONS: A real-time reconstruction technique was developed for adaptive radiotherapy using a Linac-MRI system. Our CS-PCA algorithm can achieve tumor contours with DC greater than 0.9 and NMSE less than 0.06 at acceleration factors of up to, and including, 10×. The reconstruction speed for the Split Bregman CS ranged from 200 to 260 ms, whereas the CS-PCA reconstruction speed ranged from 5 to 20 ms implemented using nonoptimized MATLAB code.

Sliding window prior data assisted compressed sensing for MRI tracking of lung tumors.

PURPOSE: Hybrid magnetic resonance imaging and radiation therapy devices are capable of imaging in real-time to track intrafractional lung tumor motion during radiotherapy. Highly accelerated magnetic resonance (MR) imaging methods can potentially reduce system delay time and/or improves imaging spatial resolution, and provide flexibility in imaging parameters. Prior Data Assisted Compressed Sensing (PDACS) has previously been proposed as an acceleration method that combines the advantages of 2D compressed sensing and the KEYHOLE view-sharing technique. However, as PDACS relies on prior data acquired at the beginning of a dynamic imaging sequence, decline in image quality occurs for longer duration scans due to drifts in MR signal. Novel sliding window-based techniques for refreshing prior data are proposed as a solution to this problem. METHODS: MR acceleration is performed by retrospective removal of data from the fully sampled sets. Six patients with lung tumors are scanned with a clinical 3 T MRI using a balanced steady-state free precession (bSSFP) sequence for 3 min at approximately 4 frames per second, for a total of 650 dynamics. A series of distinct pseudo-random patterns of partial k-space acquisition is generated such that, when combined with other dynamics within a sliding window of 100 dynamics, covers the entire k-space. The prior data in the sliding window are continuously refreshed to reduce the impact of MR signal drifts. We intended to demonstrate two different ways to utilize the sliding window data: a simple averaging method and a navigator-based method. These two sliding window methods are quantitatively compared against the original PDACS method using three metrics: artifact power, centroid displacement error, and Dice's coefficient. The study is repeated with pseudo 0.5 T images by adding complex, normally distributed noise with a standard deviation that reduces image SNR, relative to original 3 T images, by a factor of 6. RESULTS: Without sliding window implemented, PDACS-reconstructed dynamic datasets showed progressive increases in image artifact power as the 3 min scan progresses. With sliding windows implemented, this increase in artifact power is eliminated. Near the end of a 3 min scan at 3 T SNR and 5× acceleration, implementation of an averaging (navigator) sliding window method improves our metrics by the following ways: artifact power decreases from 0.065 without sliding window to 0.030 (0.031), centroid error decreases from 2.64 to 1.41 mm (1.28 mm), and Dice coefficient agreement increases from 0.860 to 0.912 (0.915). At pseudo 0.5 T SNR, the improvements in metrics are as follows: artifact power decreases from 0.110 without sliding window to 0.0897 (0.0985), centroid error decreases from 2.92 mm to 1.36 mm (1.32 mm), and Dice coefficient agreements increases from 0.851 to 0.894 (0.896). CONCLUSIONS: In this work we demonstrated the negative impact of slow changes in MR signal for longer duration PDACS dynamic scans, namely increases in image artifact power and reductions of tumor tracking accuracy. We have also demonstrated sliding window implementations (i.e., refreshing of prior data) of PDACS are effective solutions to this problem at both 3 T and simulated 0.5 T bSSFP images.

Neural-network based autocontouring algorithm for intrafractional lung-tumor tracking using Linac-MR.

PURPOSE: To develop a neural-network based autocontouring algorithm for intrafractional lung-tumor tracking using Linac-MR and evaluate its performance with phantom and in-vivo MR images. METHODS: An autocontouring algorithm was developed to determine both the shape and position of a lung tumor from each intrafractional MR image. A pulse-coupled neural network was implemented in the algorithm for contrast improvement of the tumor region. Prior to treatment, to initiate the algorithm, an expert user needs to contour the tumor and its maximum anticipated range of motion in pretreatment MR images. During treatment, however, the algorithm processes each intrafractional MR image and automatically generates a tumor contour without further user input. The algorithm is designed to produce a tumor contour that is the most similar to the expert's manual one. To evaluate the autocontouring algorithm in the author's Linac-MR environment which utilizes a 0.5 T MRI, a motion phantom and four lung cancer patients were imaged with 3 T MRI during normal breathing, and the image noise was degraded to reflect the image noise at 0.5 T. Each of the pseudo-0.5 T images was autocontoured using the author's algorithm. In each test image, the Dice similarity index (DSI) and Hausdorff distance (HD) between the expert's manual contour and the algorithm generated contour were calculated, and their centroid positions were compared (Δd centroid). RESULTS: The algorithm successfully contoured the shape of a moving tumor from dynamic MR images acquired every 275 ms. From the phantom study, mean DSI of 0.95-0.96, mean HD of 2.61-2.82 mm, and mean Δd centroid of 0.68-0.93 mm were achieved. From the in-vivo study, the author's algorithm achieved mean DSI of 0.87-0.92, mean HD of 3.12-4.35 mm, as well as Δd centroid of 1.03-1.35 mm. Autocontouring speed was less than 20 ms for each image. CONCLUSIONS: The authors have developed and evaluated a lung tumor autocontouring algorithm for intrafractional tumor tracking using Linac-MR. The autocontouring performance in the Linac-MR environment was evaluated using phantom and in-vivo MR images. From the in-vivo study, the author's algorithm achieved 87%-92% of contouring agreement and centroid tracking accuracy of 1.03-1.35 mm. These results demonstrate the feasibility of lung tumor autocontouring in the author's laboratory's Linac-MR environment.

48 echo T2 myelin imaging of white matter in first-episode schizophrenia: evidence for aberrant myelination.

Myelin water imaging provides a novel strategy to assess myelin integrity and corresponding clinical relationships in psychosis, of particular relevance in frontal white matter regions. In the current study, T2 myelin water imaging was used to assess the myelin water fraction (MWF) signal from frontal areas in a sample of 58 individuals experiencing first-episode psychosis (FEP) and 44 healthy volunteers. No differences in frontal MWF were observed between FEP subjects and healthy volunteers; however, differences in normal patterns of associations between frontal MWF and age, education and IQ were seen. Significant positive relationships between frontal MWF and age, North American Adult Reading Test (NAART) IQ, and years of completed education were observed in healthy volunteers. In contrast, only the relationship between frontal MWF and NAART IQ was significant after Bonferroni correction in the FEP group. Additionally, significant positive relationships between age and MWF in the anterior and posterior internal capsules, the genu, and the splenium were observed in healthy volunteers. In FEP subjects, only the relationship between age and MWF in the splenium was statistically significant. Frontal MWF was not associated with local white matter volume. Altered patterns of association between age, years of education, and MWF in FEP suggest that subtle disturbances in myelination may be present early in the course of psychosis.

Prior data assisted compressed sensing: a novel MR imaging strategy for real time tracking of lung tumors.

PURPOSE: Hybrid radiotherapy-MRI devices promise real time tracking of moving tumors to focus the radiation portals to the tumor during irradiation. This approach will benefit from the increased temporal resolution of MRI's data acquisition and reconstruction. In this work, the authors propose a novel spatial-temporal compressed sensing (CS) imaging strategy for the real time MRI--prior data assisted compressed sensing (PDACS), which aims to improve the image quality of the conventional CS without significantly increasing reconstruction times. METHODS: Conventional 2D CS requires a random sampling of partial k-space data, as well as an iterative reconstruction that simultaneously enforces the image's sparsity in a transform domain as well as maintains the fidelity to the acquired k-space. PDACS method requires the additional acquisition of the prior data, and for reconstruction, it additionally enforces fidelity to the prior k-space domain similar to viewsharing. In this work, the authors evaluated the proposed PDACS method by comparing its results to those obtained from the 2D CS and viewsharing methods when performed individually. All three methods are used to reconstruct images from lung cancer patients whose tumors move and who are likely to benefit from lung tumor tracking. The patients are scanned, using a 3T MRI, under free breathing using the fully sampled k-space with 2D dynamic bSSFP sequence in a sagittal plane containing lung tumor. These images form a reference set for the evaluation of the partial k-space methods. To create partial k-space, the fully sampled k-space is retrospectively undersampled to obtain a range of acquisition acceleration factors, and reconstructed with 2D-CS, PDACS, and viewshare methods. For evaluation, metrics assessing global image artifacts as well as tumor contour shape fidelity are determined from the reconstructed images. These analyses are performed both for the original 3T images and those at a simulated 0.5T equivalent noise level. RESULTS: In the 3.0T images, the PDACS strategy is shown to give superior results compared to viewshare and conventional 2D CS using all metrics. The 2D-CS tends to perform better than viewshare at the low acceleration factors, while the opposite is true at the high acceleration factors. At simulated 0.5T images, PDACS method performs only marginally better than the viewsharing method, both of which are superior compared to 2D CS. The PDACS image reconstruction time (0.3 s/image) is similar to that of the conventional 2D CS. CONCLUSIONS: The PDACS method can potentially improve the real time tracking of moving tumors by significantly increasing MRI's data acquisition speeds. In 3T images, the PDACS method does provide a benefit over the other two methods in terms of both the overall image quality and the ability to accurately and automatically contour the tumor shape. MRI's data acquisition may be accelerated using the simpler viewsharing strategy at the lower, 0.5T magnetic field, as the marginal benefit of the PDACS method may not justify its additional reconstruction times.

Evaluation of a lung tumor autocontouring algorithm for intrafractional tumor tracking using low-field MRI: a phantom study.

PURPOSE: The first aim of this study is to investigate the feasibility of online autocontouring of tumor in low field MR images (0.2 and 0.5 T) by means of a phantom and simulation study for tumor-tracking in linac-MR systems. The second aim of this study is to develop an MR compatible, lung tumor motion phantom. METHODS: An autocontouring algorithm was developed to determine both the position and shape of a lung tumor from each intra fractional MR image. To initiate the algorithm, an expert user contours the tumor and its maximum anticipated range of motion (herein termed the Background) using pretreatment scan data. During treatment, the algorithm processes each intrafractional MR image and automatically contours the tumor. To evaluate this algorithm, the authors built a phantom that replicates the low field contrast parameters (proton density, T(1), T(2)) of lung tumors and healthy lung parenchyma. This phantom allows simulation of MR images with the expected lung tumor CNR at 0.2 and 0.5 T by using a single 3 T scanner. Dynamic bSSFP images (approximately 4 images per second) are acquired while the phantom undergoes a series of preprogrammed motions based on patient lung tumor motion data. These images are autocontoured off-line using our algorithm. The fidelity of autocontouring is assessed by comparing autocontoured tumor shape and its centroid position to the actual tumor shape and its position. RESULTS: The algorithm successfully contoured the shape of a moving tumor model from dynamic MR images acquired every 275 ms. Dice's coefficients of > 0.96 and > 0.93 are achieved in 0.5 and 0.2 T equivalent images, respectively. Also, the algorithm tracked tumor position during dynamic studies, with root mean squared error (RMSE) values of < 0.55 and < 0.92 mm for 0.5 and 0.2 T equivalent images, respectively. Autocontouring speed is approximately 5 ms for each image. CONCLUSIONS: Dice's coefficients of > 0.96 and > 0.93 are achieved between autocontoured and real tumor shapes, and the position of a tumor can be tracked with RMSE values of < 0.55 and < 0.92 mm in 0.5 and 0.2 T equivalent images, respectively. These results demonstrate the feasibility of lung tumor autocontouring in low field MR images, and, by extension, intrafractional lung tumor tracking with our laboratory's linac-MR system.