- Open Access
An open-source software tool for the generation of relaxation time maps in magnetic resonance imaging
BMC Medical Imaging volume 10, Article number: 16 (2010)
In magnetic resonance (MR) imaging, T1, T2 and T2* relaxation times represent characteristic tissue properties that can be quantified with the help of specific imaging strategies. While there are basic software tools for specific pulse sequences, until now there is no universal software program available to automate pixel-wise mapping of relaxation times from various types of images or MR systems. Such a software program would allow researchers to test and compare new imaging strategies and thus would significantly facilitate research in the area of quantitative tissue characterization.
After defining requirements for a universal MR mapping tool, a software program named MRmap was created using a high-level graphics language. Additional features include a manual registration tool for source images with motion artifacts and a tabular DICOM viewer to examine pulse sequence parameters. MRmap was successfully tested on three different computer platforms with image data from three different MR system manufacturers and five different sorts of pulse sequences: multi-image inversion recovery T1; Look-Locker/TOMROP T1; modified Look-Locker (MOLLI) T1; single-echo T2/T2*; and multi-echo T2/T2*. Computing times varied between 2 and 113 seconds. Estimates of relaxation times compared favorably to those obtained from non-automated curve fitting. Completed maps were exported in DICOM format and could be read in standard software packages used for analysis of clinical and research MR data.
MRmap is a flexible cross-platform research tool that enables accurate mapping of relaxation times from various pulse sequences. The software allows researchers to optimize quantitative MR strategies in a manufacturer-independent fashion. The program and its source code were made available as open-source software on the internet.
Magnetic resonance (MR) imaging is a complex imaging modality, which has gained widespread use in modern medicine. Signal intensity in conventional MR images is influenced by a multitude of physical phenomena. In particular, there are three time constants describing the behavior of the net magnetization vector M in an MR experiment: 1) the longitudinal or spin-lattice relaxation time T1, describing the recovery of the Mz component of M; 2) the transversal or spin-spin relaxation time T2, describing the decay of the Mxy component of M; and 3) T2*, which in contrast to T2 also includes the loss of phase coherence due to field inhomogeneities and susceptibility effects. The degree to which these time constants determine signal intensity in an MR image depend on the technical parameters that are used for image acquisition [1, 2]. In clinical MR imaging, all three time constants represent characteristic magnetic properties of a given tissue , and changes from their normal values can be used to identify pathological states of that tissue. Examples for recent research efforts include investigations into the relaxation behavior of human brain in patients with multiple sclerosis [4, 5], studies of T1 and T2* changes under pharmacological stress in coronary artery disease , quantification of iron overload of the heart and liver in thalassaemia major , and analysis of myocardial fibrosis in aortic regurgitation .
Due to the composite nature of the MR signal, it is not possible to acquire raw images with "pure", quantifiable T1 or T2 properties in a direct fashion. In fact, to obtain pure T1 or T2 information, it is necessary to acquire a set of raw images that use varying acquisition parameters, and to perform multi-parameter curve fitting analysis on this raw data based on the mathematical functions that describe the underlying physical processes . If this is done on a pixel-by-pixel basis, so called parametric "maps" can be created. These maps provide a visualization of the T1 or T2 properties in a quantitative fashion, since the signal intensity of each pixel in such a map directly reflects the relaxation time calculated (typically in milliseconds).
In the past, a number of image acquisition schemes have been developed to enable measurement and mapping of MR relaxation times. A recent offspring of these techniques is called modified Look-Locker inversion recovery (MOLLI) and opens the door for high-resolution T1 mapping in cardiac applications such as the assessment of heart muscle scarring in patients with heart attacks [9, 10]. From a post-processing point of view, most of these techniques have used proprietary software programs for map reconstruction. Where mapping procedures were embedded into standard software packages for image analysis, again only specific image data were processed, and there was rarely any information available about the actual processing algorithms used. So far, both the lack of easily accessible software tools and the uncertainty about the mode of action of "black-box" software packages have posed a significant obstacle for many non-computer-expert researchers in the medical field to study relaxation time changes in diseased tissues.
The aim of our project was to provide a simple, universal software tool that can be used on multiple computer platforms to create relaxation time maps from any image data acquired with multiple mapping schemes including MOLLI .
After a survey among MR scientists at two different MR centers (Franz-Volhard-Klinik, Berlin, Germany; Leeds General Infirmary, Leeds, UK) who were asked to list desirable features of future MR mapping software, the following basic specifications were defined:
– Reading of DICOM (digital imaging and communications in medicine) image data from different MR systems
– Image registration to correct for misregistration of source images
– T1 mapping from standard inversion recovery, conventional Look-Locker or TOMROP , and modified Look-Locker inversion recovery (MOLLI) data sets
– T2 and T2* mapping from single-echo and multi-echo pulse sequences
– Control over basic computing parameters
– Graphical illustration of fitting results
– Export of maps in DICOM format for post-processing with conventional MR software
– Export of maps in standard graphic formats for illustration purposes.
A high-level programming language including routines for handling of complex image data and providing the option to run programs with free runtime licenses on all major computer platforms was selected (IDL 7.0, ITT Visual Information Solutions, Boulder CO, USA) . A software tool named "MRmap" was created according to the specifications requested  (see Additional file 1). Maps are calculated pixel-by-pixel according to the selection of the user, if the type of source images is suitable for the selected technique. Validity checks include the number of source images, common field-of-views, and appropriate timing information within the DICOM headers.
T1 maps can be generated from sets of inversion recovery images (multiple series, each containing one image) with varying inversion time (TI), Look-Locker images with varying effective TI (one series containing multiple images), or MOLLI images with varying effective TI (one series containing multiple images). After sorting the images by their corresponding TI, 3-parameter curve fitting using a Levenberg-Marquardt algorithm is performed for each pixel position where signal intensity of any of the corresponding source pixels is above the user-defined noise level. The following mathematical description of the T1 behavior is used:
with A representing a scaling factor for signal intensity; B reflecting the quality of the inversion (a B value of 2 × A means perfect inversion); and TI = inversion time. If standard source images with magnitude reconstruction are used, the information regarding the initial sign of the MR signal (which might have been negative since inversion recovery-based techniques are used) is lost. Therefore, MRmap performs several cycles of additional curve fittings on the whole data set , where each time the signs of the signal in an increasing number of leading images are switched to negative, starting with the image with the lowest TI and finishing with the last image where TI < 0.67 × limitT1 [ms]. Finally, the combination of signs yielding the best (= lowest) Chi value during the curve fitting process is selected and its curve fitting results are used to reconstruct the map. As a consequence of this procedure, the computing time is depending on the number of cycles necessary for the varying of the signs, which in turn is determined by the T1 limit (limitT1) set by the user. For Look-Locker and MOLLI data sets, where the repetitive image acquisition itself deflects the relaxation curve to be measured, the result of the curve fitting represents T1*. MRmap automatically corrects for this effect  using the remaining two curve fitting parameters A and B according to
For T2 mapping, single-echo (= multiple series with one image each) or multi-echo (= all images within one series) data sets can be used. First, all images are sorted by their echo time (TE). Then, for each pixel set, a Levenberg-Marquardt two-parameter curve fitting  is performed for
The same procedures (but different source images) can be used for T2* mapping. However, there is some evidence in the literature indicating that T2* accuracy might benefit from the use of a constant offset. Thus, MRmap optionally calculates T2* maps by performing three-parameter curve fitting for
Figure 1 shows the main screen of MRmap as implemented according to the predefined criteria. In the upper section of the screen, essential DICOM information of the images is presented in a tabular form to facilitate selection of source images. Preferences include limits of relaxation times and noise, manual registration (Figure 2), anonymization, and color tables can be adjusted on four tabs panels. After selecting data and preferences, mapping procedures or viewing of DICOM headers of the image data can be initiated with the buttons from the bar in the center of the screen. While DICOM data will be listed in tabular form on a separate window, completed maps will be displayed on the lower left of the main screen. Mouse-clicking anywhere on the maps will show the target region with 5 fold magnification in the zoom window next to the viewer, and will produce on the lower right a graph of the raw data, together with the corresponding fitting curve and the fitting results of the pixel that was clicked upon.
The functionality of the software was tested for three different operating systems (Microsoft Windows XP pro, Apple OS X 10.5, and Fedora Linux 10) in different types of image data from multiple MR systems. Accuracy of the automated computation of relaxation times was assessed by comparing mean relaxation times of regions of interest (ROIs) from maps generated by MRmap to results from non-automated curve fitting (Prism 5, GraphPad software, La Jolla/Ca, USA) of corresponding ROIs of raw images from standardized image data sets of gadolinium-doped agarose gel phantoms with multiple T1 and T2 times (8 images per set, matrix 128 × 128, TI for IR, LL, and MOLLI: 100 to 1400 ms) acquired on a clinical 1.5 T MR system (Avanto, Siemens Medical Solutions, Erlangen, Germany). Performance of the sofware was measured in the same data sets using a 2.5 GHz Intel Core 2 Duo processor. The readability of exported maps was assessed for the following standard MR software packages: a) Mass 6.0 (Medis Medical Imaging Systems, Leiden, Netherlands), b) CMR42 (Circle Cardiovascular Imaging Inc., Calgary, Canada), c) Osirix 3.3 (Antoine Rosset, Geneva, Switzerland). MRmap and its source code were made available as open-source software under the GNU General Public License .
Figure 3 contains maps created with MRmap from various sorts of source data. Figure 4 illustrates the agreement of the calculations of T1 and T2(*) times performed by MRmap as compared to non-automated measurements. There were only minimal differences, with a small bias (0.1 ms in MOLLI to 0.74 ms in multi-echo T2*) towards underestimation of relaxation times by MRmap. Differences >1 ms were only detectable with the Look-Locker approach in a phantom with very short T1, where there were multiple pixels within the ROI without converging of the fitting algorithm due to poor quality of the raw images. MRmap automatically assigns a value of zero to such pixels, thus regions with poor fitting quality are visualized as areas of signal drop out, and ROIs of such regions might yield shorter average values. Table 1 lists the computation times of the different mapping procedures for the standardized data sets using different noise levels (Figure 5). Maps exported by MRmap in DICOM format could successfully be read and analyzed by the three software packages tested.
To our knowledge, MRmap is the first open-source software tool that enables parametric T1, T2, and T2* mapping of DICOM source images on a pixel-by-pixel basis from multiple MR systems in a flexible fashion.
Other than conventional solutions that are embedded into vendor-specific application packages, the computing routines used by MRmap are well documented (see Additional file 2). Thus, results can easily be verified and do not come from a "black box", which facilitates their use for research purposes. MRmap provides a whole set of mapping routines that covers the most popular pulse sequence schemes. Our tests show that the resulting maps of the underlying relaxation times achieve the same accuracy as (tedious) non-automated curve fitting does.
By nature, the choice of a high-level graphics language as a software environment causes longer computing times than code that is written in low-level languages such as C++. This is the case in MRmap as well, where some of the mapping procedures - depending on mapping type and resolution, parameters chosen, and computer system used - can exceed one minute. However, research analyses are usually performed off-line, and therefore speed is less critical than in clinical applications. Furthermore, the exclusion of non-relevant pixels by setting appropriate noise levels allows reducing computation times drastically (see Table 1).
Apart from pure mapping procedures, MRmap provides a manual registration tool that helps to optimize mapping results if the source data contains significant motion artifacts. Image registration is a prerequisite for use in cardiac applications, where breathing artifacts are particularly common. As a limitation, MRmap in its current version does not support non-rigid registration, which can be a problem in cases where arrhythmia causes severe mis-triggering of the cardiac cycle. The integrated DICOM viewer facilitates the exploration of the source data, e.g. if images from external sites are to be analyzed. In contrast to most conventional DICOM viewers, DICOM headers from both multiple images and multiple series are listed in a tabular fashion, enabling direct comparison of header data between different images or series.
As a limitation, MRmap does not support less commonly used mapping schemes such as varying flip angle approaches . These might be implemented in future versions of the software.
MRmap is a flexible open-source software tool for the creation of parametric maps of MR relaxation times. Manual registration of source images, visualization of fitting results and data export in multiple image formats are supported. The software might facilitate research activities in the field of quantitative MR tissue analysis.
Availability and requirements
Project name: MRmap
Project home page: http://sourceforge.net/projects/mrmap
Operating systems: platform independent
Program language: IDL 7.0
Other requirements: (free) IDL Virtual Machine 7.0 (or higher) ; X11 (on Mac OS X)
License: GNU General Public License (GPL) Restrictions to use by non-academic: MRmap is intended for research purposes only.
Kaldoudi E, Williams CR: Relaxation time measurements in NMR imgaging. Part I: Longitudinal relaxation time. Concepts in Magnetic Resonance. 1993, 5: 217-242. 10.1002/cmr.1820050303.
Kingsley PB: Methods of measuring spin-lattice (T1) relaxation times: An annotated bibliography. Concepts in Magnetic Resonance. 1999, 11 (4): 243-276. 10.1002/(SICI)1099-0534(1999)11:4<243::AID-CMR5>3.0.CO;2-C.
Bottomley PA, Foster TH, Argersinger RE, Pfeifer LM: A review of normal tissue hydrogen NMR relaxation times and relaxation mechanisms from 1-100 MHz: Dependence on tissue type, NMR frequency, temperature, species, excision, and age. Medical Physics. 1984, 11 (4): 425-448. 10.1118/1.595535.
Vaithianathar L, Tench CR, Morgan PS, Lin X, Blumhardt LD: White matter T(1) relaxation time histograms and cerebral atrophy in multiple sclerosis. J Neurol Sci. 2002, 197 (12): 45-50. 10.1016/S0022-510X(02)00044-8.
Parry A, Clare S, Jenkinson M, Smith S, Palace J, Matthews PM: White matter and lesion T1 relaxation times increase in parallel and correlate with disability in multiple sclerosis. J Neurol. 2002, 249 (9): 1279-1286. 10.1007/s00415-002-0837-7.
Wacker CM, Bock M, Hartlep AW, Beck G, van Kaick G, Ertl G, Bauer WR, Schad LR: Changes in myocardial oxygenation and perfusion under pharmacological stress with dipyridamole: assessment using T2* and T1 measurements. Magn Reson Med. 1999, 41 (4): 686-695. 10.1002/(SICI)1522-2594(199904)41:4<686::AID-MRM6>3.0.CO;2-9.
Anderson LJ, Holden S, Davis B, Prescott E, Charrier CC, Bunce NH, Firmin DN, Wonke B, Porter J, Walker JM, et al: Cardiovascular T2-star (T2*) magnetic resonance for the early diagnosis of myocardial iron overload. Eur Heart J. 2001, 22 (23): 2171-2179. 10.1053/euhj.2001.2822.
Sparrow P, Messroghli DR, Reid S, Ridgway JP, Bainbridge G, Sivananthan MU: Myocardial T1 mapping for detection of left ventricular myocardial fibrosis in chronic aortic regurgitation: pilot study. AJR Am J Roentgenol. 2006, 187 (6): W630-635. 10.2214/AJR.05.1264.
Messroghli DR, Radjenovic A, Kozerke S, Higgins DM, Sivananthan MU, Ridgway JP: Modified Look-Locker inversion recovery (MOLLI) for high-resolution T1 mapping of the heart. Magn Reson Med. 2004, 52 (1): 141-146. 10.1002/mrm.20110.
Messroghli DR, Plein S, Higgins DM, Walters K, Jones TR, Ridgway JP, Sivananthan MU: Human Myocardium: Single-Breath-hold MR T1 Mapping with High Spatial Resolution--Reproducibility Study. Radiology. 2006, 238 (3): 1004-1012. 10.1148/radiol.2382041903.
Graumann R, Barfuß H, Hentschel D, Oppelt A: TOMROP: eine Sequenz zur Bestimmung der Längsrelaxationszeit T1 in der Kernspintomographie. electromedica. 1987, 55 (2): 67-72.
IDL Virtual Machine from ITT Visual Information Solutions. [http://www.ittvis.com/Downloads.aspx]
MRmap project website on sourceforge.net. [http://sourceforge.net/projects/mrmap]
Nekolla S, Gneiting T, Syha J, Deichmann R, Haase A: T1 maps by K-space reduced snapshot-FLASH MRI. J Comput Assist Tomogr. 1992, 16 (2): 327-332. 10.1097/00004728-199203000-00031.
Deichmann R, Haase A: Quantification of T1 values by SNAPSHOT-FLASH NMR imaging. J Magn Reson. 1992, 96: 608-612.
Sperber GO, Ericsson A, Hemmingsson A, Jung B, Thuomas KA: Improved formulae for signal amplitudes in repeated NMR sequences: applications in NMR imaging. Magn Reson Med. 1986, 3 (5): 685-698. 10.1002/mrm.1910030505.
Ghugre NR, Enriquez CM, Coates TD, Nelson MD, Wood JC: Improved R2* measurements in myocardial iron overload. J Magn Reson Imaging. 2006, 23 (1): 9-16. 10.1002/jmri.20467.
Preibisch C, Deichmann R: Influence of RF spoiling on the stability and accuracy of T1 mapping based on spoiled FLASH with varying flip angles. Magn Reson Med. 2009, 61 (1): 125-135. 10.1002/mrm.21776.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2342/10/16/prepub
This work was supported by a Marie Curie Reintegration Grant for D.M. by The European Commission. The DICOM write routines are based on code from "DICOM_WRITER v0.21" by Bhautik Joshi (email@example.com).
The authors declare that they have no competing interests.
DM designed and wrote the software, carried out the phantom tests and drafted the manuscript. AR and HAA acquired the in-vivo data and tested the software. RW participated in the design of the study and tested the software. TK, RD and JSM participated in the study design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Messroghli, D.R., Rudolph, A., Abdel-Aty, H. et al. An open-source software tool for the generation of relaxation time maps in magnetic resonance imaging. BMC Med Imaging 10, 16 (2010). https://doi.org/10.1186/1471-2342-10-16
- Source Image
- Magnetic Resonance System
- Quantitative Magnetic Resonance
- DICOM Viewer
- DICOM Header