Crohn’s disease is an inflammatory bowel disease mostly affecting motility in terminal ileum of small bowel. In this study, cine magnetic resonance enterography scans were used to assess the terminal ileum motility. Motility was quantified using optical flow based and gradient based analysis. ROC statistical analysis showed that immotility and motility were separable with 87% accuracy when analyzed with optical flow based algorithm and 89% accuracy with gradient based algorithm. The best classification accuracy of 90.5% was obtained when both optical flow and gradient based analysis results were used as features to train a kNN algorithm with 15-fold cross validation.
2D coronal plane balanced-SSFP sequence (True-FISP/Siemens), standard body matrix coil, and 3-tesla MR scanner (Siemens Skyra) were used for dynamic cine MRI acquisitions. An experienced radiologist defined terminal ileum as the region of interest (ROI) in each dynamic series. For motility assessment, 55 datasets of 31 patients (14 female, 17 male, aged 15-62) were used. Terminal ileum was marked as immotile in 16 datasets (9 patients), having reduced motility in 16 datasets (9 patients), and normal to high motility in 23 datasets (13 patients) by the radiologist. The motility was quantified by using two different methods, which were optical flow based analysis and gradient based analysis. All findings of entire MRE data and clinical work-ups were accepted as gold standard results for comparison.
Optical Flow Analysis: Dynamic images were analyzed with optical flow algorithm, which calculates velocity vector based on brightness of each pixel inplane to quantify motility 5,6. Magnitude of velocity was taken as a measure of motility for each pixel. Average motility was calculated within the defined ROI of terminal ileum marked by the radiologist.
Gradient Based Analysis: In each dynamic series, signal intensity–time curves were generated for each pixel. Generated curves were smoothed with a smoothing filter and gradients were calculated to generate gradient–time curves to quantify the change in each pixel. The average of gradients at each time point was taken as a measure of motility for each pixel. The mean motility within the defined ROI was calculated. Motility maps were also generated for the slice showing the small bowel by using each method (Figure 1).
For classification, weighted and medium k nearest neighbor (kNN) and medium Gaussian support vector machine (SVM) algorithms were trained with optical flow based and gradient based motility scores separately and together to classify patients as immotile and motile (normal-to-high or reduced motility). 15-folds cross validation was used. Mean scores were also analyzed with receiver operating characteristic curve (ROC) to find a cut-off value between each groups.
1. Algin, O., et al., Magnetic resonance enterography findings of chronic radiation enteritis. Cancer Imaging, 2011. 11: p. 189-94.
2. Evrimler, S. and Algin O., MR enterography with oral contrast agent composed of methylcellulose, low-dose barium sulfate, sorbitol, and lactulose: assessment of diagnostic performance, reliability, image quality, and patient tolerance. Clin Imaging, 2016. 40(3): p. 523-30. 3. Odille, F., et al., Quantitative assessment of small bowel motility by nonrigid registration of dynamic MR images. Magn Reson Med, 2012. 68(3): p. 783-93.
4. Hahnemann, M.L., et al., Quantitative assessment of small bowel motility in patients with Crohn's disease using dynamic MRI. Neurogastroenterol Motil, 2015. 27(6): p. 841-8.
5. Hahnemann, M.L., et al., Motility mapping as evaluation tool for bowel motility: initial results on the development of an automated color-coding algorithm in cine MRI. J Magn Reson Imaging, 2015. 41(2): p. 354-60.
6. MATLABCENTRAL. Horn-Schunck Optical Flow Method. MathWorks, Natick, MA. Available at: http://www.mathworks.de/matlabcentral/fileexchange/22756-horn-schunck-optical-flow-method