一种医学图像局部结构的尺度选择方法

0介绍

图像内局部结构的尺度信息是图像处理中一个重要的参数。在医学图像处理领域，为血管、肺结节等局部结构选择适合它们的尺度，对进一步处理、分析这些图像有至关重要的作用。例如，在计算机辅助的肺部诊断中，对肺结节进行分割是一项极为有意义的工作，而合理地估计出肺结节区域的尺度大小将有助于分割这些病灶区域〔1〕。除了图像分割上的应用外，在去除医学图像中的噪声信号时也常常需要对尺度进行估计。例如，利用小波将超声图像分解为细节部分后，通过估计图像内细节部分(噪声)的方差(尺度)，可以达到图像去噪的目的。然而，据我们了解，除了Narayanan等[3]估计待配准图像中未匹配区域的尺度信息，Hahn等在配准多模图像时，提出自适应地选取各向异性核函数的尺度外，大多数的研究者尚未对图像配准中局部结构尺度的计算给予足够关注。这也促使我们在本文中提出一种适用于医学图像局部结构尺度估计的方法，以满足后续医学图像配准中对局部结构尺度信息的需求。总体上说，我们可以通过两大类的方法为图像中的局部结构选择合适的尺度，即插人法(plug-in method)和拟合法(goodness-of-fit method)。插人法在实际信号与估计信号的均方差达到最小时，通过对信号偏差与方差的估计来计算信号的尺度!相对而言，拟合法并不估计信号的偏差和方差，而是在交叉验证(cross-validation、赤池信息量准则0和C>0是用于调节光滑项和残差项对公式(3)作用的经验参数，在本文中它们都被设置为考虑到相邻局部结构间存在尺度相关性，将马尔科夫随机场(Markov random field,MRF)模型也纳人到局部结构尺度的选择中。这里的马尔科夫随机场模型是通过图(graph)G=(v,动定义的.表示图上顶点(vertex)的集合，:为连接这些顶点的边(edge)的集合。在本文中，图像I(x)上的局部结构S被抽象为这里的顶点.

2实验结果及讨论

通过两个医学图像的实例，展示我们尺度估计的结果，以不同灰度等级表示不同的尺度大小(最小尺度为I时，分配最暗的灰度，最大尺度为IS时，分配最亮的灰度)。目前我们仅通过目测不同分辨率下的尺度估计效果来印证尺度估计的合理性。第一组实验是为一幅血管数字减影图像(di gital subtraction an}iography,DSA)中局部结构选择尺度(如图4所示)。在这幅血管的减影图(图4(a))中，不同粗细的血管构成了图像中具有不同尺度的局部结构。图4(b)一(e)展示了不同图像分辨率下计算得到的尺度图。在这些尺度图中，灰度较暗区域被赋予尺度条中较小的尺度值，它们对应于图4(a)中的小尺度结构(如较细的血管)，而尺度图中灰度较亮区域则分配到尺度条中较大的尺度值，它们对应于图4(a)中较为光滑的区域。在低分辨率下，由于超像素局部区域划分得较大，因此图像中一些较细的血管并不可区分，它们只能被赋予与周围区域一致的较大尺度值而呈现较亮的灰度(图4(c)。随着图像分辨的提高，这些较细的血管逐步被分配较小的尺度值而显示更暗的灰度(图4(b)一(e)。相对地，低分辨率下可被识别的粗血管，则逐步被赋予尺度条中较大的尺度值而展现出更亮的灰度。第二组实验对视网膜图像中的局部结构尺度进行计算，如图5所示。图5(a)中视网膜内的血管、视盘等构成了尺度大小不一的局部结构。图5(b)一(e)展示了不同图像分辨率下的尺度图。可以看出，在低分辨率下，局部结构的尺度分配较为单一，因此无法反映出图像内局部结构，特别是视网膜内血管的真正尺度图5(b)。随着图像分辨率的提高，超像素分块逐渐被细化，尺度图中尺度的变换也更丰富(图5(c)一(e))。血管与周围区域逐步可分，并被赋予不同的灰度。到图像最高分辩时，尺度图(图5(e))中血管区域灰度最暗，表明分配给它们的尺度最小，而视盘、背景等相对光滑的均质区域，则呈现出较亮的灰度，说明分配给它们的尺度比较大。3结论与展望本文简要回顾了医学图像中尺度估计的应用，同时提出了一种数据驱动的局部结构尺度选择方法，以满足后续医学图像配准中对局部结构尺度信息的需求。实验结果初步表明了我们方法的合理性。在后续的工作中，将对该方法计算得到的尺度做出更为详细的评价。此外，还会将该方法运用到我们基于联合显著图和自适应核回归的非刚性图像配准，通过将其与描述两幅图像间局部结构的未匹配尺度相结合，自动地调节核回归中核函数的局部尺度，用以控制参与浮动图像形变场重构的位移矢量采样量，实现局部结构匹配更加精确的非刚性医学图像配准。

致谢:感谢R.Achanta等人提供SLIC算法的开源代码对本项目的支持。

参考文献

[1]Okada K,Comaniciu D,Krishnan A.Robust anisotropic Gaussian fitting for volumetric characterization of pulmonary nodules in multislice CT[J].{H}IEEE Transactions on Medical Imaging,2005,(03):409-423.

[2]Gupta S,Chauhan RC,Saxena SC.Locally adaptive wavelet domain Bayesian processor for denoising medical ultrasound images using Speckle modeling based on Rayleigh distribution[J].IEEE Proc Vis Image Signal Proc,2005,(01):129-135.

[3]Narayanan R,Fessler JA,Ma B.Local mismatch location and spatial scale detection in image registration[J].Proc Of SPIE Med Imag,2007.1-8.

[4]Hahn DA,Daum V,Hornegger J.Automatic parameter selection for multimodal image registration[J].IEEE Tran med imag,2010,(05):1140-1155.

[5]Katkovnik V,Foi A,Egiazarian K.From local kernel to nonlocal multiple-model image denoising[J].Int J Compu Vis,2010,(01):1-32.

[6]Gomez G,Marroquin JL,Sucar LE.Probabilistic estimation of local scale[A].2000.790-793.

[7]Sun J,Xu ZB.Scale selection for anisotropic diffusion filter by Markov random field model[J].{H}Pattern Recognition,2010,(08):2630-2645.

[8]Lindeberg.Feature detection with automatic scale selection[J].{H}International Journal of Computer Vision,1998,(02):77-116.

[9]Comaniciu.An algorithm for data-driven bandwidth selection[J].{H}IEEE Transactions on Pattern Analysis and Machine Intelligence,2003,(02):281-288.

[10]Achanta R,Shaji A,Smith K.SLIC superpixels compared to state-of-the-art superpixel methods[J].{H}IEEE Transactions on Pattern Analysis and Machine Intelligence,2012,(11):2274-2281.