Anders Garpebring

Improved quality of computed tomography substitute derived from magnetic resonance (MR) data by incorporation of spatial information – potential application for MR-only radiotherapy and attenuation correction in positron emission tomography

Abstract Background. Estimation of computed tomography (CT) equivalent data, i.e. a substitute CT (s-CT), from magnetic resonance (MR) images is a prerequisite both for attenuation correction of positron emission tomography (PET) data acquired with a PET/MR scanner and for dose calculations in an MR-only radiotherapy workflow. It has previously been shown that it is possible to estimate Hounsfield numbers based on MR image intensities, using ultra short echo-time imaging and Gaussian mixture regression (GMR). In the present pilot study we investigate the possibility to also include spatial information in the GMR, with the aim to improve the quality of the s-CT. Material and methods. MR and CT data for nine patients were used in the present study. For each patient, GMR models were created from the other eight patients, including either both UTE image intensities and spatial information on a voxel by voxel level, or only UTE image intensities. The models were used to create s-CT images for each respective patient. Results. The inclusion of spatial information in the GMR model improved the accuracy of the estimated s-CT. The improvement was most pronounced in smaller, complicated anatomical regions as the inner ear and post-nasal cavities. Conclusions. This pilot study shows that inclusion of spatial information in GMR models to convert MR data to CT equivalent images is feasible. The accuracy of the s-CT is improved and the spatial information could make it possible to create a general model for the conversion applicable to the whole body.

ADC texture—An imaging biomarker for high‐grade glioma?

PURPOSE Survival for high-grade gliomas is poor, at least partly explained by intratumoral heterogeneity contributing to treatment resistance. Radiological evaluation of treatment response is in most cases limited to assessment of tumor size months after the initiation of therapy. Diffusion-weighted magnetic resonance imaging (MRI) and its estimate of the apparent diffusion coefficient (ADC) has been widely investigated, as it reflects tumor cellularity and proliferation. The aim of this study was to investigate texture analysis of ADC images in conjunction with multivariate image analysis as a means for identification of pretreatment imaging biomarkers. METHODS Twenty-three consecutive high-grade glioma patients were treated with radiotherapy (2 Gy/60 Gy) with concomitant and adjuvant temozolomide. ADC maps and T1-weighted anatomical images with and without contrast enhancement were collected prior to treatment, and (residual) tumor contrast enhancement was delineated. A gray-level co-occurrence matrix analysis was performed on the ADC maps in a cuboid encapsulating the tumor in coronal, sagittal, and transversal planes, giving a total of 60 textural descriptors for each tumor. In addition, similar examinations and analyses were performed at day 1, week 2, and week 6 into treatment. Principal component analysis (PCA) was applied to reduce dimensionality of the data, and the five largest components (scores) were used in subsequent analyses. MRI assessment three months after completion of radiochemotherapy was used for classifying tumor progression or regression. RESULTS The score scatter plots revealed that the first, third, and fifth components of the pretreatment examinations exhibited a pattern that strongly correlated to survival. Two groups could be identified: one with a median survival after diagnosis of 1099 days and one with 345 days, p = 0.0001. CONCLUSIONS By combining PCA and texture analysis, ADC texture characteristics were identified, which seems to hold pretreatment prognostic information, independent of known prognostic factors such as age, stage, and surgical procedure. These findings encourage further studies with a larger patient cohort.

Effects of inflow and radiofrequency spoiling on the arterial input function in dynamic contrast-enhanced MRI: A combined phantom and simulation study

The arterial input function is crucial in pharmacokinetic analysis of dynamic contrast‐enhanced MRI data. Among other artifacts in arterial input function quantification, the blood inflow effect and nonideal radiofrequency spoiling can induce large measurement errors with subsequent reduction of accuracy in the pharmacokinetic parameters. These errors were investigated for a 3D spoiled gradient‐echo sequence using a pulsatile flow phantom and a total of 144 typical imaging settings. In the presence of large inflow effects, results showed poor average accuracy and large spread between imaging settings, when the standard spoiled gradient‐echo signal equation was used in the analysis. For example, one of the investigated inflow conditions resulted in a mean error of about 40% and a spread, given by the coefficient of variation, of 20% for Ktrans. Minimizing inflow effects by appropriate slice placement, combined with compensation for nonideal radiofrequency spoiling, significantly improved the results, but they remained poorer than without flow (e.g., 3–4 times larger coefficient of variation for Ktrans). It was concluded that the 3D spoiled gradient‐echo sequence is not optimal for accurate arterial input function quantification and that correction for nonideal radiofrequency spoiling in combination with inflow minimizing slice placement should be used to reduce the errors. Magn Reson Med, 2011.

A Novel Estimation Method for Physiological Parameters in Dynamic Contrast-Enhanced MRI: Application of a Distributed Parameter Model Using Fourier-Domain Calculations

Dynamic contrast-enhanced magnetic resonance imaging (MRI) is a promising tool in the evaluation of tumor physiology. From rapidly acquired images and a model for contrast agent pharmacokinetics, physiological parameters are derived. One pharmacokinetic model, the tissue homogeneity model, enables estimation of both blood flow and vessel permeability together with parameters that describe blood volume and extracellular extravascular volume fraction. However, studies have shown that parameter estimation with this model is unstable. Therefore, several initial guesses are needed for accurate estimates, which makes the estimation slow. In this study a new estimation algorithm for the tissue homogeneity model, based on Fourier domain calculations, was derived and implemented as a Matlab program. The algorithm was tested with Monte-Carlo simulations and the results were compared to an existing method that uses the adiabatic approximation. The algorithm was also tested on data from a metastasis in the brain. The comparison showed that the new algorithm gave more accurate results on the 2.5th and 97.5th percentile levels, for instance the error in blood volume was reduced by 21%. In addition, the time needed for the computations was reduced with a factor 25. It was concluded that the new algorithm can be used to speed up parameter estimation while accuracy can be gained at the same time.

INTERNAL FIDUCIAL MARKERS AND SUSCEPTIBILITY EFFECTS IN MRI—SIMULATION AND MEASUREMENT OF SPATIAL ACCURACY

BACKGROUND It is well-known that magnetic resonance imaging (MRI) is preferable to computed tomography (CT) in radiotherapy target delineation. To benefit from this, there are two options available: transferring the MRI delineated target volume to the planning CT or performing the treatment planning directly on the MRI study. A precondition for excluding the CT study is the possibility to define internal structures visible on both the planning MRI and on the images used to position the patient at treatment. In prostate cancer radiotherapy, internal gold markers are commonly used, and they are visible on CT, MRI, x-ray, and portal images. The depiction of the markers in MRI are, however, dependent on their shape and orientation relative the main magnetic field because of susceptibility effects. In the present work, these effects are investigated and quantified using both simulations and phantom measurements. METHODS AND MATERIALS Software that simulated the magnetic field distortions around user defined geometries of variable susceptibilities was constructed. These magnetic field perturbation maps were then reconstructed to images that were evaluated. The simulation software was validated through phantom measurements of four commercially available gold markers of different shapes and one in-house gold marker. RESULTS Both simulations and phantom measurements revealed small position deviations of the imaged marker positions relative the actual marker positions (<1 mm). CONCLUSION Cylindrical gold markers can be used as internal fiducial markers in MRI.

CT substitutes derived from MR images reconstructed with parallel imaging.

PURPOSE Computed tomography (CT) substitute images can be generated from ultrashort echo time (UTE) MRI sequences with radial k-space sampling. These CT substitutes can be used as ordinary CT images for PET attenuation correction and radiotherapy dose calculations. Parallel imaging allows faster acquisition of magnetic resonance (MR) images by exploiting differences in receiver coil element sensitivities. This study investigates whether non-Cartesian parallel imaging reconstruction can be used to improve CT substitutes generated from shorter examination times. METHODS The authors used gridding as well as two non-Cartesian parallel imaging reconstruction methods, SPIRiT and CG-SENSE, to reconstruct radial UTE and gradient echo (GE) data into images of the head for 23 patients. For each patient, images were reconstructed from the full dataset and from a number of subsampled datasets. The subsampled datasets simulated shorter acquisition times by containing fewer radial k-space spokes (1000, 2000, 3000, 5000, and 10,000 spokes) than the full dataset (30,000 spokes). For each combination of patient, reconstruction method, and number of spokes, the reconstructed UTE and GE images were used to generate a CT substitute. Each CT substitute image was compared to a real CT image of the same patient. RESULTS The mean absolute deviation between the CT number in CT substitute and CT decreased when using SPIRiT as compared to gridding reconstruction. However, the reduction was small and the CT substitute algorithm was insensitive to moderate subsampling (≥ 5000 spokes) regardless of reconstruction method. For more severe subsampling (≤ 3000 spokes), corresponding to acquisition times less than a minute long, the CT substitute quality was deteriorated for all reconstruction methods but SPIRiT gave a reduction in the mean absolute deviation of down to 25 Hounsfield units compared to gridding. CONCLUSIONS SPIRiT marginally improved the CT substitute quality for a given number of radial spokes as compared to gridding. However, the increased reconstruction time of non-Cartesian parallel imaging reconstruction is difficult to motivate from this improvement. Because the CT substitute algorithm was insensitive to moderate subsampling, data for a CT substitute could be collected in as little as minute and reconstructed with gridding without deteriorating the CT substitute quality.

Uncertainty estimation in dynamic contrast-enhanced MRI

Using dynamic contrast‐enhanced MRI (DCE‐MRI), it is possible to estimate pharmacokinetic (PK) parameters that convey information about physiological properties, e.g., in tumors. In DCE‐MRI, errors propagate in a nontrivial way to the PK parameters. We propose a method based on multivariate linear error propagation to calculate uncertainty maps for the PK parameters. Uncertainties in the PK parameters were investigated for the modified Kety model. The method was evaluated with Monte Carlo simulations and exemplified with in vivo brain tumor data. PK parameter uncertainties due to noise in dynamic data were accurately estimated. Noise with standard deviation up to 15% in the baseline signal and the baseline T1 map gave estimated uncertainties in good agreement with the Monte Carlo simulations. Good agreement was also found for up to 15% errors in the arterial input function amplitude. The method was less accurate for errors in the bolus arrival time with disagreements of 23%, 32%, and 29% for Ktrans, ve, and vp, respectively, when the standard deviation of the bolus arrival time error was 5.3 s. In conclusion, the proposed method provides efficient means for calculation of uncertainty maps, and it was applicable to a wide range of sources of uncertainty. Magn Reson Med 69:992–1002, 2013.

Registration accuracy for MR images of the prostate using a subvolume based registration protocol

BackgroundIn recent years, there has been a considerable research effort concerning the integration of magnetic resonance imaging (MRI) into the external radiotherapy workflow motivated by the superior soft tissue contrast as compared to computed tomography. Image registration is a necessary step in many applications, e.g. in patient positioning and therapy response assessment with repeated imaging. In this study, we investigate the dependence between the registration accuracy and the size of the registration volume for a subvolume based rigid registration protocol for MR images of the prostate.MethodsTen patients were imaged four times each over the course of radiotherapy treatment using a T2 weighted sequence. The images were registered to each other using a mean square distance metric and a step gradient optimizer for registration volumes of different sizes. The precision of the registrations was evaluated using the center of mass distance between the manually defined prostates in the registered images. The optimal size of the registration volume was determined by minimizing the standard deviation of these distances.ResultsWe found that prostate position was most uncertain in the anterior-posterior (AP) direction using traditional full volume registration. The improvement in standard deviation of the mean center of mass distance between the prostate volumes using a registration volume optimized to the prostate was 3.9 mm (p < 0.001) in the AP direction. The optimum registration volume size was 0 mm margin added to the prostate gland as outlined in the first image series.ConclusionsRepeated MR imaging of the prostate for therapy set-up or therapy assessment will both require high precision tissue registration. With a subvolume based registration the prostate registration uncertainty can be reduced down to the order of 1 mm (1 SD) compared to several millimeters for registration based on the whole pelvis.

Haralick texture features from apparent diffusion coefficient (ADC) MRI images depend on imaging and pre-processing parameters

In recent years, texture analysis of medical images has become increasingly popular in studies investigating diagnosis, classification and treatment response assessment of cancerous disease. Despite numerous applications in oncology and medical imaging in general, there is no consensus regarding texture analysis workflow, or reporting of parameter settings crucial for replication of results. The aim of this study was to assess how sensitive Haralick texture features of apparent diffusion coefficient (ADC) MR images are to changes in five parameters related to image acquisition and pre-processing: noise, resolution, how the ADC map is constructed, the choice of quantization method, and the number of gray levels in the quantized image. We found that noise, resolution, choice of quantization method and the number of gray levels in the quantized images had a significant influence on most texture features, and that the effect size varied between different features. Different methods for constructing the ADC maps did not have an impact on any texture feature. Based on our results, we recommend using images with similar resolutions and noise levels, using one quantization method, and the same number of gray levels in all quantized images, to make meaningful comparisons of texture feature results between different subjects.

Combining phase and magnitude information for contrast agent quantification in dynamic contrast-enhanced MRI using statistical modeling

The purpose of this study was to investigate, using simulations, a method for improved contrast agent (CA) quantification in DCE‐MRI.