Patrik Brynolfsson

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.

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.

Motion corrected DWI with integrated T2-mapping for simultaneous estimation of ADC, T2-relaxation and perfusion in prostate cancer

OBJECTIVE Multiparametric magnetic resonance imaging (MRI) and PI-RADS (Prostate Imaging - Reporting and Data System) has become the standard to determine a probability score for a lesion being a clinically significant prostate cancer. T2-weighted and diffusion-weighted imaging (DWI) are essential in PI-RADS, depending partly on visual assessment of signal intensity, while dynamic-contrast enhanced imaging is less important. To decrease inter-rater variability and further standardize image evaluation, complementary objective measures are in need. METHODS We here demonstrate a sequence enabling simultaneous quantification of apparent diffusion coefficient (ADC) and T2-relaxation, as well as calculation of the perfusion fraction f from low b-value intravoxel incoherent motion data. Expandable wait pulses were added to a FOCUS DW SE-EPI sequence, allowing the effective echo time to change at run time. To calculate both ADC and f, b-values 200s/mm2 and 600s/mm2 were chosen, and for T2-estimation 6 echo times between 64.9ms and 114.9ms were used. RESULTS Three patients with prostate cancer were examined and all had significantly decreased ADC and T2-values, while f was significantly increased in 2 of 3 tumors. T2 maps obtained in phantom measurements and in a healthy volunteer were compared to T2 maps from a SE sequence with consecutive scans, showing good agreement. In addition, a motion correction procedure was implemented to reduce the effects of prostate motion, which improved T2-estimation. CONCLUSIONS This sequence could potentially enable more objective tumor grading, and decrease the inter-rater variability in the PI-RADS classification.

Density estimation of grey-level co-occurrence matrices for image texture analysis

The Haralick texture features are common in the image analysis literature, partly because of their simplicity and because their values can be interpreted. It was recently observed that the Haralick texture features are very sensitive to the size of the GLCM that was used to compute them, which led to a new formulation that is invariant to the GLCM size. However, these new features still depend on the sample size used to compute the GLCM, i.e. the size of the input image region-of-interest (ROI). The purpose of this work was to investigate the performance of density estimation methods for approximating the GLCM and subsequently the corresponding invariant features. Three density estimation methods were evaluated, namely a piece-wise constant distribution, the Parzen-windows method, and the Gaussian mixture model. The methods were evaluated on 29 different image textures and 20 invariant Haralick texture features as well as a wide range of different ROI sizes. The results indicate that there are two types of features: those that have a clear minimum error for a particular GLCM size for each ROI size, and those whose error decreases monotonically with increased GLCM size. For the first type of features, the Gaussian mixture model gave the smallest errors, and in particular for small ROI sizes (less than about [Formula: see text]). In conclusion, the Gaussian mixture model is the preferred method for the first type of features (in particular for small ROIs). For the second type of features, simply using a large GLCM size is preferred.

Technical Note: Adapting a GE SIGNA PET/MR scanner for radiotherapy

PURPOSE Simultaneous collection of PET and MR data for radiotherapy purposes are useful for, for example, target definition and dose escalations. However, a prerequisite for using PET/MR in the radiotherapy workflow is the ability to image the patient in treatment position. The aim of this work was to adapt a GE SIGNA PET/MR scanner to image patients for radiotherapy treatment planning and evaluate the impact on signal-to-noise (SNR) of the MR images, and the accuracy of the PET attenuation correction. METHOD A flat tabletop and a coil holder were developed to image patients in the treatment position, avoid patient contour deformation, and facilitate attenuation correction of flex coils. Attenuation corrections for the developed hardware and an anterior array flex coil were also measured and implemented to the PET/MR system to minimize PET quantitation errors. The reduction of SNR in the MR images due to the added distance between the coils and the patient was evaluated using a large homogenous saline-doped water phantom, and the activity quantitation errors in PET imaging were evaluated with and without the developed attenuation corrections. RESULT We showed that the activity quantitation errors in PET imaging were within ±5% when correcting for attenuation of the flat tabletop, coil holder, and flex coil. The SNR of the MRI images were reduced to 74% using the tabletop, and 66% using the tabletop and coil holders. CONCLUSION We present a tabletop and coil holder for an anterior array coil to be used with a GE SIGNA PET/MR scanner, for scanning patients in the radiotherapy work flow. Implementing attenuation correction of the added hardware from the radiotherapy setup leads to acceptable PET image quantitation. The drop in SNR in MR images may require adjustment of the imaging protocols.

Evaluation of advanced MR techniques for development of early biomarkers for treatment efficacy in malignant brain tumors

Evaluation of advanced MR techniques for development of early biomarkers for treatment efficacy in malignant brain tumors