The Stockwell Transform has the potential to perform multi-resolution texture analysis in magnetic resonance imaging (MRI). However, it is computationally intensive and memory demanding. The polar Stockwell Transform (PST) is rotation-invariant and relatively memory efficient, but still computationally demanding. The new Discrete Orthogonal Stockwell Transform (DOST) appears to have addressed both the computation and storage challenges; however, its utility in localized texture analysis remains unclear. Our goal was to investigate the theory and texture analysis ability of the DOST versus PST using both synthetic and MR images, and explore the relative importance of the associated texture features using a simple classification example based on clinical brain MRI of six multiple sclerosis patients. MRI texture analysis focused on FLAIR images, and the classification used a machine learning algorithm, random forest, that differentiated regions of interest (ROIs) into 2 classes: white matter lesions, and the contralateral normal-appearing white matter (control). Our results showed that the PST features had a greater ability in detecting subtle changes in image structure than the DOST and polar-index DOST (PDOST). Quantitatively, based on 187 lesion and 187 control ROIs, both the PST and the rotation-invariant radial PST performed better in the classification than the DOST and PDOST, where the latter were no better than guessing (p = 0.65 and 0.98). Further analysis using a hierarchical random forest showed that combining MRI signal intensity with the PST or DOST predictions increased the classification performance, with the accuracy, sensitivity, and specificity all improved to >85% in the tests. Collectively, the DOST is less competitive than the PST in localized image texture analysis. The PST features may help with texture-based lesion classification in MS based on clinical brain MRI scans following further verification.
Read full abstract