Automated slice positioning for 2D MRA in bolus tracking of DCE-MRI
Takao Goto1 and Mirai Araki1

1MR Engineering, GE Healthcare, Hino-shi, Japan

Synopsis

Accurate placement of a 2D plane across the aorta while examining scout images is a complex task and makes the operator's workflow difficult in bolus tracking of DCE-MRI. We present a new method for automated slice positioning for 2D MRA used to monitor bolus arrival. The 2D plane was planned by aorta detection using both Hough Forests and AdaBoost classifiers following the classification of axial images. A dataset with 40 patients was tested, and 35 cases depicted the cross section of the aorta clearly. This automation will help the operator and decrease the total study time.

PURPOSE

In bolus tracking techniques monitoring the arrival of the bolus, the operator prescribes the 2D plane of MRA crossing the aorta and manually starts the dynamic contrast-enhanced (DCE) scan. The prescription of the 2D plane is a complex process, and its accuracy is highly dependent on the skill of the operator. Our goal is to position the 2D plane automatically, thus resulting in the improvement of operator workflow.

METHODS

To prescribe a 2D plane cross-section of the aorta, accurate detection of the aorta location is necessary. Ordinary 2D SSFSE scout images are analyzed to detect the aorta. Fig. 1 shows the flow of our algorithm. Firstly, axial images (> 6 slices) are classified by anatomy, including the lung, the boundary of the lung, the liver, and the liver with two kidneys. We used Laplacian eigenmaps1 for this classification. The axial image is cropped, inscribed outside of the thorax, and then resized to a rectangular shape (44 × 22 pixels). Following classification, one of four atlases that correspond to the four classifications is applied to the axial image, which can narrow down the existing area of the aorta. The axial images are also used for the detection of cerebral spinal fluid (CSF) that is searched by the AdaBoost classifier and becomes an anchor point in searching for the aorta2. To detect the aorta, we used the Hough Forests classifier3 that utilizes the feature value in the neighboring patches of objects, and creates a map of votes from each feature patch. The patch (16 × 16 pixels) is extracted from a sub-window (48 × 48 pixels) in the axial image, and the sub-window is rotated around the CSF to detect the aorta in the next axial image (Fig. 2). The Hough map is weighted by the atlas. The location of the maximum value in the map shows the aorta location. There are two aorta classes in SSFSE axial images. One appears bright and the other appears dark depending on whether flow spins are re-phased or de-phased due to the blood velocity. The Hough Forests classifier detects both the bright and dark aortas, and the AdaBoost classifier determines whether it is dark or bright2. The esophagus, which is near the aorta, has a similar contrast and shape to the aorta. This can cause misdetections and false positives. We built a model expressing the location distribution of the aorta in relation to the esophagus (Fig. 3). With this model, another peak apart from the incorrectly detected peak in the Hough map is searched, and another peak is identified as the aorta peak if the peak value is more than 50% of the esophagus. The oblique plane is positioned by the line crossing the center of the detected aortas (7 aortas is needed at least) using the regression method (Fig. 1d). The oblique plane was reformatted from 3D datasets of DCE liver images acquired in the aortic phase to simulate an oblique MRA plane. We tested 40 patients (mean age = 64 years; age range = 24 – 85 years) after informed consent was obtained under institutional review and approval.

RESULTS

The Laplacian eigenmaps consisted of 218 cropped images. The accuracy rate was approximately 80%. For aorta detection, 10,000 patches extracted from 200 axial images of 43 patient datasets were used as the training step for the Hough Forests classifier. Aside from these learning datasets, 358 axial images were used for testing from 40 patients. Table 1 shows the error of the detected center position of the aorta, which was calculated and compared to the manual results (ground truth). In the resulting reformatted oblique images, 5 of 40 datasets did not show 50% of the cross-section of the aorta. Fig. 4 shows an example of the 2D oblique planes for success and failure cases.

DISCUSSION

In aorta detection, most of the error was caused in the detection of the bright aorta because of the difficult identification of the center of the aorta. In particular, the occurrence frequency of the bright aorta in patients tended to be greater than in volunteers that we tested previously2 due to slow vessel flow. Therefore, it is critical to decrease the occurrence of the bright aorta.

CONCLUSION

We proposed a new method for automated positioning of MRA for bolus tracking, and demonstrated that our algorithm was able to plan the 2D scan plane for MRA. This automation will help the operator and decrease the total study time.

Acknowledgements

No acknowledgement found.

References

[1]. Belkin M. Niyogi P. Neural Compt, 2003;15:1373-1396

[2]. Goto T. Kabasawa H. MRI, 2015;33:63–71

[3]. Gall J. Yao Angela. et al. PAMI, 2011;33:2188-2201

Figures

Flow of algorithm. (a) Body classifier using Laplacian eigenmaps chooses atlas from four that restrict the search region of aorta. (b,c) Hough Forests classifier detects aortas following the CSF detection. (d) Oblique plane of MRA determined by detected aorta locations using linear regression method.

Pipeline detection of the aorta in every axial image from the superior to the inferior direction (a – c). Hough map area is rotated counterclockwise around the CSF.

Method for removal of miss-detection of esophagus. (a) Detected esophagus and the corresponding Hough map. (b) Distribution model of aorta locations for the detected esophagus peak point. (c) Detected peak of aorta by suppressing the esophagus peak.

Reformatted images of 3D DCE in aortic phase using linear regression. The DCE image was acquired in the same examination as scout images. (a, c) are computed from manually determined aortas. (b, d) are computed from automatically detected aortas. (b) depicted the cross section of the aorta better than (d).

Euclidian distance error



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
0248