Mayerhoefer Marius E, Szomolanyi Pavol, Jirak Daniel, Berg Andreas, Materka Andrzej, Dirisamer Albert, Trattnig Siegfried
Department of Radiology, MR Center, Medical University of Vienna, Vienna, Austria.
Invest Radiol. 2009 Jul;44(7):405-11. doi: 10.1097/RLI.0b013e3181a50a66.
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.
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 x 30 mm(2)), but with variations of matrix size (MTX): 16 x 16; 32 x 32; 64 x 64; 128 x 128; and 256 x 256. Original images were interpolated to higher matrix sizes (up to 256 x 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.
For images based on a very low original resolution (MTX = 16 x 16), misclassification rates remained high, despite the use of interpolation. For higher resolution images (MTX = 32 x 32 and 64 x 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 x 16), zero-fill interpolation outperformed the two pixel interpolation techniques, for images based on higher original resolutions (MTX = 32 x 32 and 64 x 64). We observed the most considerable improvements after interpolation by a factor of 2 or 4.
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.
(1)确定在像素或k空间层面进行磁共振(MR)图像插值是否能改善基于纹理的模式分类结果;(2)比较图像插值对不同类别的纹理特征的影响,以及它们区分不同模式的能力。
我们使用配备专用显微成像梯度插入件的3.0特斯拉扫描仪,获取了两组各3个具有不同结节模式的聚苯乙烯球体和琼脂凝胶(PSAG)体模的T2加权多切片多回波MR图像(球体直径:PSAG - 1,0.8 - 1.25毫米;PSAG - 2,1.25 - 2.0毫米;PSAG - 3,2.0 - 3.15毫米)。图像数据集由每组20个连续的轴向切片组成,视野恒定(30×30毫米²),但矩阵大小(MTX)有所变化:16×16;32×32;64×64;128×128;以及256×256。原始图像通过线性和三次B样条(像素层面)以及零填充(k空间层面)插值法插值到更高的矩阵大小(最高256×256)。对于原始图像数据集和插值后的图像数据集,分别计算从共生矩阵(COC)、游程长度矩阵(RUN)、绝对梯度(GRA)、自回归模型和小波变换(WAV)导出的纹理特征。基于使用第一组PSAG体模(训练数据集)的图像计算的Fisher系数确定的每个类别的3个最佳纹理特征,进行k均值聚类,以区分属于第二组体模(测试数据集)的PSAG - 1、PSAG - 2和PSAG - 3图像。对所有原始图像数据集和插值后的图像数据集独立进行此操作。错误分类数据向量的比率用作主要结果指标。
对于基于非常低的原始分辨率(MTX = 16×16)的图像,尽管使用了插值,错误分类率仍然很高。对于更高分辨率的图像(MTX = 32×32和64×64),插值增强了除WAV之外的所有类别中纹理特征区分3种体模的能力。这种积极效果在COC和RUN特征中尤为明显,在较小程度上,GRA特征也是如此。插值后,未观察到WAV特征有一致的改善,甚至有一些负面影响。虽然在非常低的分辨率(MTX = 16×16)下没有任何一种插值技术具有明显优势,但对于基于更高原始分辨率(MTX = 32×32和64×64)的图像,零填充插值优于两种像素插值技术。我们观察到插值2倍或4倍后改善最为显著。
基于COC、RUN和GRA特征,MR图像插值有潜力改善模式分类结果。除非空间分辨率非常差,否则零填充是首选的插值技术,建议最大插值因子为4。