Effects of Magnetic Resonance Image Interpolation on the Results of Texture-Based Pattern Classification: A Phantom Study

Objectives:To (1) determine whether magnetic resonance (MR) image interpolation at the pixel or k-space level can improve the results of texture-based pattern classification, and (2) compare the effects of image interpolation on texture features of different categories, with regard to their ability to distinguish between different patterns. Materials and Methods:We obtained T2-weighted, multislice multiecho MR images of 2 sets of each 3 polystyrene spheres and agar gel (PSAG) phantoms with different nodular patterns (sphere diameter: PSAG-1, 0.8-1.25 mm; PSAG-2, 1.25-2.0 mm; PSAG-3, 2.0-3.15 mm), using a 3.0 Tesla scanner equipped with a dedicated microimaging gradient insert. Image datasets, which consisted of 20 consecutive axial slices each, were obtained with a constant field of view (30 × 30 mm2), but with variations of matrix size (MTX): 16 × 16; 32 × 32; 64 × 64; 128 × 128; and 256 × 256. Original images were interpolated to higher matrix sizes (up to 256 × 256) by means of linear and cubic B-spline (pixel level) as well as zero-fill (k-space level) interpolation. For both original and interpolated image datasets, texture features derived from the co-occurrence (COC) and run-length matrix (RUN), absolute gradient (GRA), autoregressive model, and wavelet transform (WAV) were calculated independently. Based on the 3 best texture features of each category, as determined by calculation of Fisher coefficients using images from the first set of PSAG phantoms (training dataset), k-means clustering was performed to separate PSAG-1, PSAG-2, and PSAG-3 images belonging to the second set of phantoms (test dataset). This was done independently for all original and interpolated image datasets. Rates of misclassified data vectors were used as primary outcome measures. Results:For images based on a very low original resolution (MTX = 16 × 16), misclassification rates remained high, despite the use of interpolation. For higher resolution images (MTX = 32 × 32 and 64 × 64), interpolation enhanced the ability of texture features, in all categories except WAV, to discriminate between the 3 phantoms. This positive effect was particularly pronounced for COC and RUN features, and to a lesser degree, also GRA features. No consistent improvements, and even some negative effects, were observed for WAV features, after interpolation. Although there was no clear superiority of any single interpolation techniques at very low resolution (MTX = 16 × 16), zero-fill interpolation outperformed the two pixel interpolation techniques, for images based on higher original resolutions (MTX = 32 × 32 and 64 × 64). We observed the most considerable improvements after interpolation by a factor of 2 or 4. Conclusions:MR image interpolation has the potential to improve the results of pattern classification, based on COC, RUN, and GRA features. Unless spatial resolution is very poor, zero-filling is the interpolation technique of choice, with a recommended maximum interpolation factor of 4.