Medicine:HARP (algorithm)

From HandWiki
Short description: Medical image analysis algorithm
HARP
Developer(s)Image Analysis and Communications Laboratory
Operating systemLinux, Mac OS X, Windows
TypeCardiac Motion Tracking
WebsiteHARP Overview (Software Download)

Harmonic phase (HARP) algorithm[1] is a medical image analysis technique capable of extracting and processing motion information from tagged magnetic resonance image (MRI) sequences. It was initially developed by N. F. Osman and J. L. Prince at the Image Analysis and Communications Laboratory at Johns Hopkins University. The method uses spectral peaks in the Fourier domain of tagged MRI, calculating the phase images of their inverse Fourier transforms, which are called harmonic phase (HARP) images. The motion of material points through time is then tracked, under the assumption that the HARP value of a fixed material point is time-invariant. The method is fast and accurate, and has been accepted as one of the most popular tagged MRI analysis methods in medical image processing.

Background

In cardiac magnetic resonance imaging, tagging techniques[2][3][4][5] make it possible to capture and store the motion information of myocardium in vivo. MR tagging uses a special pulse sequence to create temporary features – tags in the myocardium. Tags deform together with the myocardium as the heart beats and are captured by MR imaging. Analysis of the motion of the tag features in many images taken from different orientations and at different times can be used to track material points in the myocardium.[6][7] Tagged MRI is widely used to develop and refine models of normal and abnormal myocardial motion[8][9][10][11] to better understand the correlation of coronary artery disease with myocardial motion abnormalities and the effects of treatment after myocardial infarction. However, suffered from long imaging and post-processing times,[12] tagged MRI was slow in entering into routine clinical use until the HARP algorithm was developed and published in 1999.[13]

Description

HARP processing

Demonstration of the HARP processing of a tagged cardiac MRI slice. (a) An MR image with vertical SPAMM tags. (b) shows the magnitude of its Fourier transform. By extracting the spectral peak inside the circle, a complex image is produced with a magnitude (c) and a phase (d).

A tagged MRI showing motion of a human heart is shown in the image (a). The effect of tagging can be described as a multiplication of the underlying image by a sinusoid tag pattern having a certain fundamental frequency, causing an amplitude modulation of the underlying image and replicating its Fourier transform into the pattern shown in (b).

HARP processing uses a bandpass filter to isolate one of the spectral peaks. For example, the circle drawn in (b) is the -3 dB isocontour of the bandpass filter used to process this data. Selection of the filters for optimal performance is discussed in this paper.[14] The inverse Fourier transform of the filtered image yields a complex harmonic image [math]\displaystyle{ I_k(\mathbf{y}, t) }[/math] at image coordinates [math]\displaystyle{ \mathbf{y} = [y_1, y_2]^T }[/math] and time [math]\displaystyle{ t }[/math]:

[math]\displaystyle{ I_k(\mathbf{y}, t)=D_k(\mathbf{y}, t)e^{j\phi_k(\mathbf{y}, t)} }[/math]

where [math]\displaystyle{ D_k }[/math] is called the harmonic magnitude image and [math]\displaystyle{ \phi_k }[/math] is called the harmonic phase image. The harmonic magnitude image in (c) extracted from a using the filter in (b) shows the geometry of the heart. And the harmonic phase image in (d) contains the motion of the myocardium in horizontal direction. In practice, tagged images from two directions (both horizontal and vertical, i.e., [math]\displaystyle{ k }[/math] is 1 and 2) are processed to provide a 2D motion map in the image plane. Notice that the harmonic phase images are computed by taking the inverse tangent of the imaginary part divided by the real part of [math]\displaystyle{ I_k(\mathbf{y}, t) }[/math], such that the range of this computation is only in [math]\displaystyle{ [-\pi, +\pi) }[/math]. In other words, d is only the wrapped value of the actual phase. We denote this principle value by [math]\displaystyle{ a_k(\mathbf{y}, t) }[/math]; it is mathematically related to the true phase by:

[math]\displaystyle{ a_k(\mathbf{y}, t)=mod(\phi_k(\mathbf{y}, t)+\pi,2\pi)-\pi }[/math]

Either [math]\displaystyle{ \phi_k }[/math] or [math]\displaystyle{ a_k }[/math] might be called a harmonic phase (HARP) image, but only [math]\displaystyle{ a_k }[/math] can be directly calculated and visualized. It is the basis for HARP tracking.

HARP tracking

Result of HARP tracking of a tagged cardiac MRI slice

For a fixed material point with a HARP value, only one of the points sharing the same HARP value in a later time frame is the correct match. If the apparent motion is small from one image to the next, it is likely that the nearest of these points is the correct point. The tracking result is very accurate in this case.[citation needed]

Consider a material point located at [math]\displaystyle{ \mathbf{y}_m }[/math] at time [math]\displaystyle{ t_m }[/math]. If [math]\displaystyle{ \mathbf{y}_{m+1} }[/math] is the apparent position of this point at time [math]\displaystyle{ t_{m+1} }[/math], we have:

[math]\displaystyle{ \phi_k(\mathbf{y}_{m+1},t_{m+1})=\phi_k(\mathbf{y}_m,t_m) }[/math]

The Newton–Raphson interactive method is used to find a solution, which is:

[math]\displaystyle{ y^{(n+1)}=y^{(n)}-[\nabla\phi_k(\mathbf{y}^{(n)},t_{m+1})]^{-1}[\phi_k(\mathbf{y}^{(n)},t_{m+1})-\phi_k(\mathbf{y}_m,t_m)] }[/math]

In practice, since [math]\displaystyle{ \phi_k }[/math] is not available, [math]\displaystyle{ a_k }[/math] is used in its place. This equation can be rewritten after a few derivations considering the "wrapping" relation between [math]\displaystyle{ \phi_k }[/math] and [math]\displaystyle{ a_k }[/math].

The result of HARP tracking of one frame of cardiac MRI is shown in the figure. It is obtained by calculating both motions from horizontal direction and vertical direction, resulting in a 2D vector field showing the motion of every material point on the myocardium at this time frame.

The entire HARP algorithm takes only a few minutes to perform on a normal computer and the motion tracking result is accurate (with a typical error range of [math]\displaystyle{ \pm1 }[/math] pixel). As a result, it is now widely adopted by the medical image analysis community as a standard processing technique for tagged MRI.

See also

References

  1. ↑ "Imaging Heart Motion Using Harmonic Phase MRI". IEEE Trans Med Imaging 19 (3): 186–202. 2000. doi:10.1109/42.845177. PMID 10875703. 
  2. ↑ Zerhouni, E.A.; Parish, D.M.; Rogers, W.J.; Yang, A.; Shapiro, E.P. (1988). "Human heart: tagging with MR imaging—a method for noninvasive assessment of myocardial motion". Radiology 169 (1): 59–63. doi:10.1148/radiology.169.1.3420283. PMID 3420283. 
  3. ↑ Axel, L.; Dougherty, L. (1989). "MR imaging of motion with spatial modulation of magnetization". Radiology 171 (3): 841–845. doi:10.1148/radiology.171.3.2717762. PMID 2717762. 
  4. ↑ McVeigh, E.R.; Atalar, E. (1992). "Cardiac tagging with breath-hold cine MRI". Magn Reson Med 28 (2): 318–327. doi:10.1002/mrm.1910280214. PMID 1461130. 
  5. ↑ Fischer, S.E.; McKinnon, G.C.; Maier, S.E.; Boesiger, P. (1993). "Improved myocardial tagging contrast". Magn Reson Med 30 (2): 191–200. doi:10.1002/mrm.1910300207. PMID 8366800. 
  6. ↑ McVeigh, E.R. (1998). "Regional myocardial function". Cardiol Clin 16 (2): 189–206. doi:10.1016/s0733-8651(05)70008-4. PMID 9627756. 
  7. ↑ McVeigh, E.R. (1996). "MRI of myocardial function: motion tracking techniques". Magn Reson Imag 14 (2): 137–150. doi:10.1016/0730-725x(95)02009-i. PMID 8847969. 
  8. ↑ Young, A.A.; Axel, L. (1992). "Three-dimensional motion and deformation of the heart wall: estimation with spatial modulation of magnetization—a model-based approach". Radiology 185 (2): 241–247. doi:10.1148/radiology.185.1.1523316. PMID 1523316. 
  9. ↑ Moore, C.; O'Dell, W.; McVeigh, E.R.; Zerhouni, E. (1992). "Calculation of three-dimensional left ventricular strains from biplanar tagged MR images". Magn Reson Imag 2 (2): 165–175. doi:10.1002/jmri.1880020209. PMID 1562767. 
  10. ↑ Clark, N.R.; Reichek, N.; Bergey, P.; Hoffman, E.A.; Brownson, D.; Palmon, L.; Axel, L. (1991). "Circumferential myocardial shortening in the normal human left ventricle". Circulation 84 (1): 67–74. doi:10.1161/01.cir.84.1.67. PMID 2060124. 
  11. ↑ McVeigh, E.R.; Zerhouni, E.A. (1991). "Noninvasive measurements of transmural gradients in myocardial strain with MR imaging". Radiology 180 (3): 677–683. doi:10.1148/radiology.180.3.1871278. PMID 1871278. 
  12. ↑ Budinger, T.F.; Berson, A.; McVeigh, E.R.; Pettigrew, R.I.; Pohost, G.M.; Watson, J.T.; Wickline, S.A. (1998). "Cardiac MR imaging: report of a working group sponsored by the National Heart, Lung, and Blood Institute". Radiology 208 (3): 573–576. doi:10.1148/radiology.208.3.9722831. PMID 9722831. 
  13. ↑ Osman, N.F.; Kerwin, W.S.; McVeigh, E.R.; Prince, J.L. (1999). "Cardiac Motion Tracking Using CINE Harmonic Phase (HARP) Magnetic Resonance Imaging". Magn Reson Med 42 (6): 1048–1060. doi:10.1002/(sici)1522-2594(199912)42:6<1048::aid-mrm9>3.3.co;2-d. PMID 10571926. 
  14. ↑ Osman, N.F.; Prince, J.L. (1998). "Motion estimation from tagged MR images using angle images". Proc Int Conf Imag Proc: 704–708. 

External links



Retrieved from "https://handwiki.org/wiki/index.php?title=Medicine:HARP_(algorithm)&oldid=3270025"