Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Current Directions in Biomedical Engineering

Joint Journal of the German Society for Biomedical Engineering in VDE and the Austrian and Swiss Societies for Biomedical Engineering

Editor-in-Chief: Dössel, Olaf

Editorial Board: Augat, Peter / Buzug, Thorsten M. / Haueisen, Jens / Jockenhoevel, Stefan / Knaup-Gregori, Petra / Kraft, Marc / Lenarz, Thomas / Leonhardt, Steffen / Malberg, Hagen / Penzel, Thomas / Plank, Gernot / Radermacher, Klaus M. / Schkommodau, Erik / Stieglitz, Thomas / Urban, Gerald A.

1 Issue per year

Open Access
See all formats and pricing
More options …

Assessing MRI susceptibility artefact through an indicator of image distortion

Alfredo Illanes / Johannes W. Krug / Michael Friebe
Published Online: 2016-09-30 | DOI: https://doi.org/10.1515/cdbme-2016-0095


Susceptibility artefacts in magnetic resonance imaging (MRI) caused by medical devices can result in a severe degradation of the MR image quality. The quantification of susceptibility artefacts is regulated by the ASTM standard which defines a manual method to assess the size of an artefact. This means that the estimated artefact size can be user dependent. To cope with this problem, we propose an algorithm to automatically quantify the size of such susceptibility artefacts. The algorithm is based on the analysis of a 3D surface generated from the 2D MR images. The results obtained by the automatic algorithm were compared to the manual measurements performed by study participants. The results show that the automatic and manual measurements follow the same trend. The clear advantage of the automated algorithm is the absence of the inter- and intra-observer variability. In addition, the algorithm also detects the slice containing the largest artefact which was not the case for the manual measurements.

Keywords: artefact assessment; MRI; susceptibility

1 Introduction

When magnetic resonance imaging (MRI) is used as a modality for guiding a minimally invasive procedures, materials used to fabricate catheters or biopsy needles have to be chosen such that their magnetic susceptibility does not distort the acquired MR images. Such susceptibility artefacts can make MRI unusable or harm its use as a guiding modality in minimally invasive interventions such as liver or breast biopsies. It can affect the image quality if the area of interest is close to the position of the passive instrument [1], [2].

To not significantly affect the quality of the diagnostic information of the MR images, testing of susceptibility artifacts arising from passive instruments is important before introducing them into the MRI environment [3], [4]. Therefore, the precise measurement of the artefact size is important during the development stage of medical tools to be used in MRI. The quantification of artefacts caused by such devices is regulated by the ASTM F2119 standard [5]. The ASTM defines a quantitative but manual method to assess the artefact, where the definition of the artefact is clearly defined but gives only an approximative method for its quantification [6]. Hence, it can be assumed that observed artefact dimensions are user dependent despite a standard exists for reducing subjectivity.

Some authors have proposed different methods to assess artifacts using mainly software measurement tools [7], [8], while other works proposed to automatize the standard measurement objective and reproducible measurement that comply with this standard [6], [9], [10].

This work investigates whether an automatic algorithm can be used as an alternative for measuring susceptibility artefacts of passive instruments used in MRI. The algorithm is based on the computation of a quantitative indicator related to the amount of distortion produced by an instrument on a three-dimensional (3D) surface obtained from the 2D MR image.

For the evaluation of the proposed algorithm, six MRI compatible biopsy needles were used for artefact assessment. The results show that an indicator can be computed automatically from all the slices of a given MRI sequence and that a value computed from the indicator is correlated with the average of the manual measurements.

2 Method

2.1 Data acquisition and manual measurements setup

Six different MR compatible biopsy needles with outer diameters of 1 mm to 1.3 mm (ITP GmbH, Somatex GmbH, Germany) were imaged in a 3T MRI scanner (Magnetom Skyra, Siemens) using a FLASH sequence (see Table 1). A Plexiglas phantom filled with a solution of water and copper sulphate was used for the experiments [5].

Table 1

Parameters of the MRI sequence used to acquire the images of the different biopsy needles.

Artefacts were manually analysed by three MRI or image analysis experts. For each needle configuration, the participants had to select one slice for measuring the artefact size. Each artefact was measured twice by each participant (with a pause of one day between the measurements) following [5] using a DICOM viewer (RadiAnt DICOM). In summary, every participant made six artefact measurements in each of the two measurement series and therefore six measurements are available for each needle.

2.2 Image distortion indicator computation

The standard 2D MR image depicted in Figure 1A shows the artefact of a biopsy needle. The artefact can be represented as a 3D surface version of the MR image by representing the values of each pixel (being located in the x-y plane) along the z-axis (see Figure 1B). In this 3D surface version the susceptibility artefact can be viewed as a distortion of the 3D surface image background. This distortion appears as a downward deflection of the background. It seems that the artefact observed in a 3D surface point of view has a more complicated nature than the traditional 2D image, therefore the main idea is to quantify this deflection of the background using a quantitative indicator of background distortion. For that we propose to calculate an indicator related to the volume resulting from the background deflection produced by the artefact in the 3D surface image.

(A) Standard 2D MR image of a needle, (B) 3D surface representation of the 2D MRI, (C) obtained 3D surface after attenuation of low frequency trend and removal of surface offset.
Figure 1

(A) Standard 2D MR image of a needle, (B) 3D surface representation of the 2D MRI, (C) obtained 3D surface after attenuation of low frequency trend and removal of surface offset.

Before doing this, it is important to observe in Figure 1B that the surface baseline is located above zero and that it has non-homogeneous characteristics with low frequency trends. This last feature can be explained by the location of the Plexiglas phantom with respect to the RF coil position and its inhomogeneous field distribution, which causes a variation of the signal intensity.

Therefore it is required first to attenuate the low frequency trend of the background surface and to remove the surface offset. Thus, the deflection volume, where no artefact is present, will tend to be zero, while the volume where the artefact is present will tend to be negative: the more negative the indicator value is, the more pronounced is the artefact level.

To remove the low frequency trend and the offset, a preprocessing step is performed over each horizontal signal line of the surface image and then the 3D surface is reconstructed using the preprocessed signal lines. In a first step the horizontal signal line passes through a moving average (MA) smoother. Then the deflection waveform (if it is present) is removed from the signal line and an hermite interpolation is used to fill the removed signal segment. A Kalman smoother [11] is then used for the low frequency trend extraction, which is finally subtracted to the MA smoothed signal. In Figure 1C the obtained 3D surface image is shown after preprocessing of every signal line and surface reconstruction.

The resulting preprocessed 3D surface of Figure 1C is only the result for one MRI slice. The main algorithm processes every signal line of every slice for calculating a final distortion indicator. A summary of this algorithm is presented in the following:

  • For each slice the preprocessing algorithm described above is applied.

  • For each line of the preprocessed 3D surface image of each image slice, an indicator of artefact level is computed using a integration operation over the preprocessed signal line.

  • All the values of the indicator of every signal line and slice are stored in the same vector.

  • The elements of the vector are sorted in descending order.

Figure 2 shows an example of the distortion indicator per slice when the algorithm was applied to 13 slices of a needle MRI image, where it is possible to observe how the indicator changes according to the visibility of the artefact. At the right of this figure, the descending order sorted vector is shown, which correspond to the final distortion indicator. In this way the artefact assessment results in a monotonically decreasing curve (MDC) indicator, which goes from zero values representing signal lines, where no artefact is present, passing through middle values and arriving to lower negative values representing the maximal artefact levels assessed in different signal lines. These lower values are the ones which are significant for assessing the artefact created by the tested instrument.

Assesment of artifact caused by a biopsy needle using the distortion indicator in 13 slices and the resulting MDC indicator.
Figure 2

Assesment of artifact caused by a biopsy needle using the distortion indicator in 13 slices and the resulting MDC indicator.

3 Results

First of all in Figure 3 the MDC indicator is shown for every tested biopsy needle (named N1, …, N6 in the sequel). Since the MDC curve is the result of a descending sorted signal, only the last part of this indicator is significant (the maximal value of artefact assessment is in this section of the indicator). That is why only the last 150 values of the MDC signal are shown. Despite that, due to the nature of the indicator the MDC, it has no distance unit but it shows clearly the different artefact sizes produced by each needle.

Last 150 values of the MDC indicator obtained for each one of the six tested needles.
Figure 3

Last 150 values of the MDC indicator obtained for each one of the six tested needles.

For comparison between manual and automatic artefact assessment, the measured manual artefact were evaluated for each needle using the average and standard deviation of all the manual measurements performed per needle. This is shown at the left of Figure 4, where the average and standard deviation of the six measures per needle (2 per each participant) result in different artefacts levels. It is also possible to observe significant variability between observers with sometimes more than 1 mm of measurements differences between them. At the right of the same figure the automatic measurements are shown, where the displayed single values per needle were calculated as the average of the last 100 values of the obtained MDC for each needle. With the exception of the differences in artefact between needles N1 and N6 the curves have an acceptable correlation.

Average and standard deviation of the manual measurements versus the automatic artefact assessment performed from the MDC indicator (absolute values were taken for the right graph).
Figure 4

Average and standard deviation of the manual measurements versus the automatic artefact assessment performed from the MDC indicator (absolute values were taken for the right graph).

The variability between participants shown in the left of Figure 4 is not only caused by different interpretations on how to measure an artefact using the ASTM standard, but also because the observers selected different slices for performing their measurements. This is revealed by Table 2, where the selected slices for each needle by each participant are shown. It is evident that not always the same slice is selected by the participants and that moreover sometimes not the same slice is selected by the same observer in two different measurements of the same artefact. The last row of the table shows the slice where the automatic algorithm finds the largest artefact. It can be noticed that this maximal artefact measurement is in almost all cases located in the slice where there is a manual selection agreement of at least 50% of all the measurements performed in a needle.

Table 2

Selection of slices by the different participants and also by the automatic algorithm for each tested needle.

4 Conclusion

The presented work had the objective of comparing artefact assessment performed by manual measurements with an automatic algorithm where the algorithms measurement principle is based on a different point of view than the proposed standards. Assessing the artefact through the computation of a 3D surface distortion indicator shows a correlation with the manually obtained results.

One question to answer in the future is if the 3D representation of the artefact can reveal more complex structures than the 2D-images can and if this can improve the accuracy in artefact assessment. Also, further work is needed to use the indicator for quantitative information about the artefact in terms of an accepted physical unity.

Author’s Statement

Research funding: This research was funded by the German BMBF (grant 03IPT7100X). Conflict of interest: Authors state no conflict of interest. Material and methods: Informed consent: Informed consent is not applicable. Ethical approval: The conducted research is not related to either human or animal use.


  • [1]

    Kaiser M, Detert M, Rube MA, El-Tahir A, Elle O J, Melzer A, et al. Resonant marker design and fabrication techniques for device visualization during interventional magnetic resonance imaging. Biomed Technik. 2015;60:89–103. Google Scholar

  • [2]

    Patil S, Bieri O, Jhooti P, Scheffler K. Automatic slice positioning (ASP) for passive real-time tracking of interventional devices using projection-reconstruction imaging with echo-dephasing (pride). Magnet Reson Med. 2009;62:935–42. Google Scholar

  • [3]

    Schaefers G. Testing MR safety and compatibility: an overview of the methods and current standards. IEEE Eng Med Biol Mag. 2008;27:23–7. Google Scholar

  • [4]

    Yu N, Gassert R, Riener R. Mutual interferences and design principles for mechatronic devices in magnetic resonance imaging. Int J Comput Assist Radiol Surg. 2011;6:473–88. Google Scholar

  • [5]

    ASTM, Conshohocken PA. ASTM Standard F2119-07, Standard test method for evaluation of MR image artifacts from passive implants. 

  • [6]

    Güttler F, Heinrich A, Teichgräber U. Software development for the determination of susceptibility artefacts in MRI after ASTM F2119. Biomed. Technik. 2012;57:480–480. Google Scholar

  • [7]

    Chen C, Chen W, Goodman S, Hargreaves BA, Koch KM, Lu W, et al. SEMAC and MAVRIC for artifact-corrected MR imaging around metal in the knee. Proceedings of the ISMRM, Stockholm, Sweden, 130; 2010. Google Scholar

  • [8]

    Liao W, Deserno TM, Spitzer K. Evaluation of free non-diagnostic dicom software tools. In Medical Imaging. p. 691903–691903. International Society for Optics and Photonics; 2008. Google Scholar

  • [9]

    Dössel O. Biomechanic models and image guided interventions. Biomed Tech (Berl). 2015;60:519. Google Scholar

  • [10]

    Heinrich A, Guettler F. Software-based automated measurement of susceptibility artifacts on magnetic resonance images. Google Scholar

  • [11]

    Sameni R, Shamsollahi MB, Jutten C, Clifford GD. A nonlinear bayesian filtering framework for ECG denoising. IEEE Trans Biomed Eng. 2007;54:2172–85. Google Scholar

About the article

Published Online: 2016-09-30

Published in Print: 2016-09-01

Citation Information: Current Directions in Biomedical Engineering, Volume 2, Issue 1, Pages 427–431, ISSN (Online) 2364-5504, DOI: https://doi.org/10.1515/cdbme-2016-0095.

Export Citation

©2016 Alfredo Illanes et al., licensee De Gruyter.. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0

Comments (0)

Please log in or register to comment.
Log in