Interactive Segmentation of MR Images from Brain Tumor Patients

Medical doctors often do not trust the result of fully automatic segmentations because they have no possibility to make corrections if necessary. On the other hand, manual corrections can introduce a user bias. In this work, we propose to integrate the possibility for quick manual corrections into a fully automatic segmentation method for brain tumor images. This allows for necessary corrections while maintaining a high objectiveness. The underlying idea is similar to the well-known Grab-Cut algorithm, but here we combine decision forest classification with conditional random field regularization for interactive segmentation of 3D medical images. The approach has been evaluated by two different users on the BraTS2012 dataset. Accuracy and robustness improved compared to a fully automatic method and our interactive approach was ranked among the top performing methods. Time for computation including manual interaction was less than 10 minutes per patient, which makes it attractive for clinical use.

RADIX: a minimal-resources rapid-access drilling system

Determining the expected age at a potential ice-core drilling site on a polar ice sheet generally depends on a combination of information from remote-sensing methods, estimates of current accumulation and modelling. This poses irreducible uncertainties in retrieving an undisturbed ice core of the desired age. Although recently perfected radar techniques will improve the picture of the ice sheet below future drilling sites, rapid prospective drillings could further increase the success of deep drilling projects. Here we design and explore a drilling system for a minimum-size rapid-access hole. The advantages of a small hole are the low demand for drilling fluid, low overall weight of the equipment, fast installing and de-installing and low costs. We show that, in theory, drilling of a 20mm hole to a depth of 3000m is possible in ~4 days. First concepts have been realized and verified in the field. Both the drill cuttings and the hole itself can be used to characterize the properties of the ice sheet and its potential to provide a trustworthy palaeo-record. A candidate drilling site could be explored in ~2 weeks, which would enable the characterization of several sites in one summer season.

Postmortem MR quantification of the heart for characterization and differentiation of ischaemic myocardial lesions.

OBJECTIVES Recently, an MRI quantification sequence has been developed which can be used to acquire T1- and T2-relaxation times as well as proton density (PD) values. Those three quantitative values can be used to describe soft tissue in an objective manner. The purpose of this study was to investigate the applicability of quantitative cardiac MRI for characterization and differentiation of ischaemic myocardial lesions of different age. MATERIALS AND METHODS Fifty post-mortem short axis cardiac 3 T MR examinations have been quantified using a quantification sequence. Myocardial lesions were identified according to histology and appearance in MRI images. Ischaemic lesions were assessed for mean T1-, T2- and proton density values. Quantitative values were plotted in a 3D-coordinate system to investigate the clustering of ischaemic myocardial lesions. RESULTS A total of 16 myocardial lesions detected in MRI images were histologically characterized as acute lesions (n = 8) with perifocal oedema (n = 8), subacute lesions (n = 6) and chronic lesions (n = 2). In a 3D plot comprising the combined quantitative values of T1, T2 and PD, the clusters of all investigated lesions could be well differentiated from each other. CONCLUSION Post-mortem quantitative cardiac MRI is feasible for characterization and discrimination of different age stages of myocardial infarction. KEY POINTS • MR quantification is feasible for characterization of different stages of myocardial infarction. • The results provide the base for computer-aided MRI cardiac infarction diagnosis. • Diagnostic criteria may also be applied for living patients.

Temperature dependence of postmortem MR quantification for soft tissue discrimination.

OBJECTIVES To investigate and correct the temperature dependence of postmortem MR quantification used for soft tissue characterization and differentiation in thoraco-abdominal organs. MATERIAL AND METHODS Thirty-five postmortem short axis cardiac 3-T MR examinations were quantified using a quantification sequence. Liver, spleen, left ventricular myocardium, pectoralis muscle and subcutaneous fat were analysed in cardiac short axis images to obtain mean T1, T2 and PD tissue values. The core body temperature was measured using a rectally inserted thermometer. The tissue-specific quantitative values were related to the body core temperature. Equations to correct for temperature differences were generated. RESULTS In a 3D plot comprising the combined data of T1, T2 and PD, different organs/tissues could be well differentiated from each other. The quantitative values were influenced by the temperature. T1 in particular exhibited strong temperature dependence. The correction of quantitative values to a temperature of 37 °C resulted in better tissue discrimination. CONCLUSION Postmortem MR quantification is feasible for soft tissue discrimination and characterization of thoraco-abdominal organs. This provides a base for computer-aided diagnosis and detection of tissue lesions. The temperature dependence of the T1 values challenges postmortem MR quantification. Equations to correct for the temperature dependence are provided. KEY POINTS • Postmortem MR quantification is feasible for soft tissue discrimination and characterization • Temperature dependence of the T1 values challenges the MR quantification approach • The results provide the basis for computer-aided postmortem MRI diagnosis • Diagnostic criteria may also be applied for living patients.

Metastases in normal-sized pelvic lymph nodes: detection with diffusion-weighted MR imaging

PURPOSE To prospectively assess the diagnostic performance of diffusion-weighted (DW) magnetic resonance (MR) imaging in the detection of pelvic lymph node metastases in patients with prostate and/or bladder cancer staged as N0 with preoperative cross-sectional imaging. MATERIALS AND METHODS This study was approved by an independent ethics committee. Written informed consent was obtained from all patients. Patients with no enlarged lymph nodes on preoperative cross-sectional images who were scheduled for radical resection of the primary tumor and extended pelvic lymph node dissection were enrolled. All patients were examined with a 3-T MR unit, and examinations included conventional and DW MR imaging of the entire pelvis. Image analysis was performed by three independent readers blinded to any clinical information. Metastases were diagnosed on the basis of high signal intensity on high b value DW MR images and morphologic features (shape, border). Histopathologic examination served as the standard of reference. Sensitivity and specificity were calculated, and bias-corrected 95% confidence intervals (CIs) were obtained with the bootstrap method. The Fleiss and Cohen κ and median test were applied for statistical analyses. RESULTS A total of 4846 lymph nodes were resected in 120 patients. Eighty-eight lymph node metastases were found in 33 of 120 patients (27.5%). Short-axis diameter of these metastases was less than or equal to 3 mm in 68, more than 3 mm to 5 mm in 13, more than 5 mm to 8 mm in five; and more than 8 mm in two. On a per-patient level, the three readers correctly detected metastases in 26 (79%; 95% CI: 64%, 91%), 21 (64%; 95% CI: 45%, 79%), and 25 (76%; 95% CI: 60%, 90%) of the 33 patients with metastases, with respective specificities of 85% (95% CI: 78%, 92%), 79% (95% CI: 70%, 88%), and 84% (95% CI: 76%, 92%). Analyzed according to hemipelvis, lymph node metastases were detected with histopathologic examination in 44 of 240 pelvic sides (18%); the three readers correctly detected these on DW MR images in 26 (59%; 95% CI: 45%, 73%), 19 (43%; 95% CI: 27%, 57%), and 28 (64%; 95% CI: 47%, 78%) of the 44 cases. CONCLUSION DW MR imaging enables noninvasive detection of small lymph node metastases in normal-sized nodes in a substantial percentage of patients with prostate and bladder cancer diagnosed as N0 with conventional cross-sectional imaging techniques.

Endovascular treatment for acute ischemic stroke patients: implications and interpretation of IMS III, MR RESCUE, and SYNTHESIS EXPANSION trials: A report from the Working Group of International Congress of Interventional Neurology

OBJECTIVE The results of Interventional Management of Stroke (IMS) III, Magnetic Resonance and REcanalization of Stroke Clots Using Embolectomy (MR RESCUE), and SYNTHESIS EXPANSION trials are expected to affect the practice of endovascular treatment for acute ischemic stroke. The purpose of this report is to review the components of the designs and methods of these trials and to describe the influence of those components on the interpretation of trial results. METHODS A critical review of trial design and conduct of IMS III, MR RESCUE, and SYNTHESIS EXPANSION is performed with emphasis on patient selection, shortcomings in procedural aspects, and methodology of data ascertainment and analysis. The influence of each component is estimated based on published literature including multicenter clinical trials reporting on endovascular treatment for acute ischemic stroke and myocardial infarction. RESULTS We critically examined the time interval between symptom onset and treatment and rates of angiographic recanalization to differentiate between "endovascular treatment" and "parameter optimized endovascular treatment" as it relates to the IMS III, MR RESCUE, and SYNTHESIS EXPANSION trials. All the three trials failed to effectively test "parameter optimized endovascular treatment" due to the delay between symptom onset and treatment and less than optimal rates of recanalization. In all the three trials, the magnitude of benefit with endovascular treatment required to reject the null hypothesis was larger than could be expected based on previous studies. The IMS III and SYNTHESIS EXPANSION trials demonstrated that rates of symptomatic intracerebral hemorrhages subsequent to treatment are similar between IV thrombolytics and endovascular treatment in matched acute ischemic stroke patients. The trials also indirectly validated the superiority/equivalence of IV thrombolytics (compared with endovascular treatment) in patients with minor neurological deficits and those without large vessel occlusion on computed tomographic/magnetic resonance angiography. CONCLUSIONS The results do not support a large magnitude benefit of endovascular treatment in subjects randomized in all the three trials. The possibility that benefits of a smaller magnitude exist in certain patient populations cannot be excluded. Large magnitude benefits can be expected with implementation of "parameter optimized endovascular treatment" in patients with ischemic stroke who are candidates for IV thrombolytics.

Outcome study of real-time MR-guided cervical periradicular injection therapy in an open 1.0 Tesla MRI system

PURPOSE To evaluate the accuracy, safety, and efficacy of cervical nerve root injection therapy using magnetic resonance guidance in an open 1.0 T MRI system. METHODS Between September 2009 and April 2012, a total of 21 patients (9 men, 12 women; mean age 47.1 ± 11.1 years) underwent MR-guided cervical periradicular injection for cervical radicular pain in an open 1.0 T system. An interactive proton density-weighted turbo spin echo (PDw TSE) sequence was used for real-time guidance of the MR-compatible 20-gauge injection needle. Clinical outcome was evaluated on a verbal numeric rating scale (VNRS) before injection therapy (baseline) and at 1 week and 1, 3, and 6 months during follow-up. RESULTS All procedures were technically successful and there were no major complications. The mean preinterventional VNRS score was 7.42 and exhibited a statistically significant decrease (P < 0.001) at all follow-up time points: 3.86 ± 1.53 at 1 week, 3.21 ± 2.19 at 1 month, 2.58 ± 2.54 at 3 months, and 2.76 ± 2.63 at 6 months. At 6 months, 14.3 % of the patients reported complete resolution of radicular pain and 38.1 % each had either significant (4-8 VNRS score points) or mild (1-3 VNRS score points) relief of pain; 9.5 % experienced no pain relief. CONCLUSION Magnetic resonance fluoroscopy-guided periradicular cervical spine injection is an accurate, safe, and efficacious treatment option for patients with cervical radicular pain. The technique may be a promising alternative to fluoroscopy- or CT-guided injections of the cervical spine, especially in young patients and in patients requiring repeat injections.

Diffusion-weighted MR imaging of upper abdominal organs: field strength and intervendor variability of apparent diffusion coefficients

PURPOSE To determine the variability of apparent diffusion coefficient (ADC) values in various anatomic regions in the upper abdomen measured with magnetic resonance (MR) systems from different vendors and with different field strengths. MATERIALS AND METHODS Ten healthy men (mean age, 36.6 years ± 7.7 [standard deviation]) gave written informed consent to participate in this prospective ethics committee-approved study. Diffusion-weighted (DW) MR imaging was performed in each subject with 1.5- and 3.0-T MR systems from each of three vendors at two institutions. Two readers independently measured ADC values in seven upper abdominal regions (left and right liver lobe, gallbladder, pancreas, spleen, and renal cortex and medulla). ADC values were tested for interobserver differences, as well as for differences related to field strength and vendor, with repeated-measures analysis of variance; coefficients of variation (CVs) and variance components were calculated. RESULTS Interreader agreement was excellent (intraclass coefficient, 0.876). ADC values were (77.5-88.8) ×10(-5) mm(2)/sec in the spleen and (250.6-278.5) ×10(-5) mm(2)/sec in the gallbladder. There were no significant differences between ADC values measured at 1.5 T and those measured at 3.0 T in any anatomic region (P >.10 for all). In two of seven regions at 1.5 T (left and right liver lobes, P < .023) and in four of seven regions at 3.0 T (left liver lobe, pancreas, and renal cortex and medulla, P < .008), intervendor differences were significant. CVs ranged from 7.0% to 27.1% depending on the anatomic location. CONCLUSION Despite significant intervendor differences in ADC values of various anatomic regions of the upper abdomen, ADC values of the gallbladder, pancreas, spleen, and kidney may be comparable between MR systems from different vendors and between different field strengths.

Comparison of (31)P saturation and inversion magnetization transfer in human liver and skeletal muscle using a clinical MR system and surface coils.

(31)P MRS magnetization transfer ((31)P-MT) experiments allow the estimation of exchange rates of biochemical reactions, such as the creatine kinase equilibrium and adenosine triphosphate (ATP) synthesis. Although various (31)P-MT methods have been successfully used on isolated organs or animals, their application on humans in clinical scanners poses specific challenges. This study compared two major (31)P-MT methods on a clinical MR system using heteronuclear surface coils. Although saturation transfer (ST) is the most commonly used (31)P-MT method, sequences such as inversion transfer (IT) with short pulses might be better suited for the specific hardware and software limitations of a clinical scanner. In addition, small NMR-undetectable metabolite pools can transfer MT to NMR-visible pools during long saturation pulses, which is prevented with short pulses. (31)P-MT sequences were adapted for limited pulse length, for heteronuclear transmit-receive surface coils with inhomogeneous B1 , for the need for volume selection and for the inherently low signal-to-noise ratio (SNR) on a clinical 3-T MR system. The ST and IT sequences were applied to skeletal muscle and liver in 10 healthy volunteers. Monte-Carlo simulations were used to evaluate the behavior of the IT measurements with increasing imperfections. In skeletal muscle of the thigh, ATP synthesis resulted in forward reaction constants (k) of 0.074 ± 0.022 s(-1) (ST) and 0.137 ± 0.042 s(-1) (IT), whereas the creatine kinase reaction yielded 0.459 ± 0.089 s(-1) (IT). In the liver, ATP synthesis resulted in k = 0.267 ± 0.106 s(-1) (ST), whereas the IT experiment yielded no consistent results. ST results were close to literature values; however, the IT results were either much larger than the corresponding ST values and/or were widely scattered. To summarize, ST and IT experiments can both be implemented on a clinical body scanner with heteronuclear transmit-receive surface coils; however, ST results are much more robust against experimental imperfections than the current implementation of IT.

3D morphometry using automated aortic segmentation in native MR angiography: an alternative to contrast enhanced MRA?

INTRODUCTION Native-MR angiography (N-MRA) is considered an imaging alternative to contrast enhanced MR angiography (CE-MRA) for patients with renal insufficiency. Lower intraluminal contrast in N-MRA often leads to failure of the segmentation process in commercial algorithms. This study introduces an in-house 3D model-based segmentation approach used to compare both sequences by automatic 3D lumen segmentation, allowing for evaluation of differences of aortic lumen diameters as well as differences in length comparing both acquisition techniques at every possible location. METHODS AND MATERIALS Sixteen healthy volunteers underwent 1.5-T-MR Angiography (MRA). For each volunteer, two different MR sequences were performed, CE-MRA: gradient echo Turbo FLASH sequence and N-MRA: respiratory-and-cardiac-gated, T2-weighted 3D SSFP. Datasets were segmented using a 3D model-based ellipse-fitting approach with a single seed point placed manually above the celiac trunk. The segmented volumes were manually cropped from left subclavian artery to celiac trunk to avoid error due to side branches. Diameters, volumes and centerline length were computed for intraindividual comparison. For statistical analysis the Wilcoxon-Signed-Ranked-Test was used. RESULTS Average centerline length obtained based on N-MRA was 239.0±23.4 mm compared to 238.6±23.5 mm for CE-MRA without significant difference (P=0.877). Average maximum diameter obtained based on N-MRA was 25.7±3.3 mm compared to 24.1±3.2 mm for CE-MRA (P<0.001). In agreement with the difference in diameters, volumes obtained based on N-MRA (100.1±35.4 cm(3)) were consistently and significantly larger compared to CE-MRA (89.2±30.0 cm(3)) (P<0.001). CONCLUSIONS 3D morphometry shows highly similar centerline lengths for N-MRA and CE-MRA, but systematically higher diameters and volumes for N-MRA.

MR Angiography at 3 Tesla to Assess Proximal Internal Carotid Artery Stenoses: Contrast-Enhanced or 3D Time-of-Flight MR Angiography?

PURPOSE The aim of this study was to compare the diagnostic accuracy of 3D time-of-flight (TOF-MRA) and contrast-enhanced (CE-MRA) magnetic resonance angiography at 3 T for detection and quantification of proximal high-grade stenosis using multidetector computed tomography angiography (MDCTA) as reference standard. METHODS The institutional ethics committee approved this prospective study. A total of 41 patients suspected of having internal carotid artery (ICA) stenosis underwent both MDCTA and MRA. CE-MRA and TOF-MRA were performed using a 3.0-T imager with a dedicated eight-element cervical coil. ICA stenoses were measured according to the North American Symptomatic Carotid Endarterectomy Trial criteria and categorized as 0-25 % (minimal), 25-50 % (mild), 50-69 % (moderate), 70-99 % (high grade), and 100 % (occlusion). Sensitivity and specificity for the detection of high-grade ICA stenoses (70-99 %) and ICA occlusions were determined. In addition, intermodality agreement was assessed with κ-statistics for detection of high-grade ICA stenoses (70-99 %) and ICA occlusions. RESULTS A total of 80 carotid arteries of 41 patients were reviewed. Two previously stented ICAs were excluded from analysis. On MDCTA, 7 ICAs were occluded, 12 ICAs presented with and 63 without a high-grade ICA stenosis (70-99 %). For detecting 70-99 % stenosis, both 3D TOF-MRA and CE-MRA were 91.7 % sensitive and 98.5 % specific, respectively. Both MRA techniques were highly sensitive (100 %), and specific (CE-MRA, 100 %; TOF-MRA, 98.7 %) for the detection of ICA occlusion. However, TOF-MRA misclassified one high-grade stenosis as occlusion. Intermodality agreement for detection of 70-99 % ICA stenoses was excellent between TOF-MRA and CE-MRA [κ = 0.902, 95 % confidence interval (CI) = 0.769-1.000], TOF-MRA and MDCTA (κ = 0.902, 95 % CI = 0.769-1.000), and CE-MRA and MDCTA (κ = 0.902, 95 % CI = 0.769-1.000). CONCLUSION Both 3D TOF-MRA and CE-MRA at 3 T are reliable tools for detecting high-grade proximal ICA stenoses (70-99 %). 3D TOF-MRA might misclassify pseudo-occlusions as complete occlusions. If there are no contraindications for CE-MRA, CE-MRA is recommended as primary MR imaging modality.

3D Intervertebral Disc Localization and Segmentation from MR Images by Data-Driven Regression and Classification

In this paper we propose a new fully-automatic method for localizing and segmenting 3D intervertebral discs from MR images, where the two problems are solved in a unified data-driven regression and classification framework. We estimate the output (image displacements for localization, or fg/bg labels for segmentation) of image points by exploiting both training data and geometric constraints simultaneously. The problem is formulated in a unified objective function which is then solved globally and efficiently. We validate our method on MR images of 25 patients. Taking manually labeled data as the ground truth, our method achieves a mean localization error of 1.3 mm, a mean Dice metric of 87%, and a mean surface distance of 1.3 mm. Our method can be applied to other localization and segmentation tasks.

Localization and Segmentation of 3D Intervertebral Discs in MR Images by Data Driven Estimation

This paper addresses the problem of fully-automatic localization and segmentation of 3D intervertebral discs (IVDs) from MR images. Our method contains two steps, where we first localize the center of each IVD, and then segment IVDs by classifying image pixels around each disc center as foreground (disc) or background. The disc localization is done by estimating the image displacements from a set of randomly sampled 3D image patches to the disc center. The image displacements are estimated by jointly optimizing the training and test displacement values in a data-driven way, where we take into consideration both the training data and the geometric constraint on the test image. After the disc centers are localized, we segment the discs by classifying image pixels around disc centers as background or foreground. The classification is done in a similar data-driven approach as we used for localization, but in this segmentation case we are aiming to estimate the foreground/background probability of each pixel instead of the image displacements. In addition, an extra neighborhood smooth constraint is introduced to enforce the local smoothness of the label field. Our method is validated on 3D T2-weighted turbo spin echo MR images of 35 patients from two different studies. Experiments show that compared to state of the art, our method achieves better or comparable results. Specifically, we achieve for localization a mean error of 1.6-2.0 mm, and for segmentation a mean Dice metric of 85%-88% and a mean surface distance of 1.3-1.4 mm.