基于改进mrf凝胶图像分割算法

《基于改进mrf凝胶图像分割算法》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、基于改进MRF凝胶图像分割算法摘要：为了得到更好的凝胶图像分割效果，文章提出了一种基于马尔科夫随机场(MRF)模型分割算法的改进算法，首先用非局部均值(NL-means)对图像进行滤波，然后通过模糊c均值算法(FCM)实现了图像的初始聚类分割,进一步通过二阶逻辑模型(MLL)模糊聚类作为先验知识获取其先验概率并因此得到后验概率，在此过程中引入了灰度点密度权值来实现聚类和均值方差的更新，最终实现凝胶图像的最优化分割。关键词：马尔科夫随机场；分割；灰度点密度；凝胶图像引言图像分割就是把一幅图像分成若干个特定的区域并提出感兴趣目标的技术和过程，也是从图像处理进入图像分

2、析的关键步骤，属于基本计算机视觉技术范畴。只有在图像分割成功的基础上才有可能实现对目标特征的提取和参数的测量，才能使更高层图像的分析和图像理解成为可能。20世纪70年代中期，二维凝胶电泳这个概念由O'Farrell与Klose等人建立[1],但那时并未有蛋白质组学这个概念，随着二维凝胶技术的逐渐成熟发展，蛋白质组学的发展也突飞猛进，直至今天的辉煌成就。二维凝胶图像［2］可以被看作是一幅数字图像，在这幅图像上，蛋白质点被扫描成不同形状、不同大小和不同灰度的点，分割是其分析的关键步骤之一，它的目的是有效的将目标蛋白质点从凝胶图像中分割出来。基于马尔可夫随机场的分割方

3、法是一种利用图像中像素间的空间相关性进行分割的方法，能够准确地描述每个像素所属类别与周围像素类别之间的重要依赖关系。国内外的研究者们已对MRF的图像分割方法进行了深入了研究，其最早是由S.German和D.German［3］提到的有关随机场概念的问题开始的，而后在1994年，Bouman等再次［4］提出了具有因果特性的马尔科夫先验模型，文献［5］利用马尔科夫随机场的MLL模型来描叙像素标号场的先验概率，文献［6］提出的HMRF-FCM图像分割算法实现了用均值场来逼近马尔科夫随机场，但是他的算法复杂度较高，计算效率较低。本文在综合分析已有文献算法的基础上，结合凝胶

4、图像的特点，提出了一种改进的MRF凝胶图像分割方法。改进算法首先利用非局部均值算法对图像进行预处理，根据像素点之间的相似度进行去噪，使其在结构和细节保护上具有一定优势，同时对各类灰度均值点密度聚类采用归一化加权方式引入权重3,使其后验概率的确定更加准确从而实现凝胶图像的分割。1算法介绍1.1马尔可夫随机场马尔可夫随机场方法［7］是建立在MRF模型与贝叶斯(Bayes)理论的基础上，提供了先验知识与不确定性描述联系的纽带，并根据观测图像的特点，利用估计理论和统计决策中某些最优准则确定分割问题的目标函数，以此来求解满足这些条件的最大可能分布，从而将分割问题转化为最优

5、化问题。基于MRF的分割方法是一种利用图像像素的空间相关性进行分割的方法，它先将图像的灰度信息进行分类，并且能够准确的描述每个像素所属分类与周围像素类别之间的相互依赖关系。1.2基于改进MRF的凝胶图像分割问题实现在马尔科夫随机场中，常用两个随机场来描述带分割的图像，假设X和Y是二维平面上的随机场，则令X是像素的标号场，一般用先验分布描述标号场的局部相关性；令Y为观测图像场，常以标号场为条件用分布函数来描述观测数据的分布。然后根据贝叶斯定理［8〜12］实现图像的分割，但是若依据传统公式进行计算，计算量很大，很难得到理想效果，所以我们假设观测场Y=y每一个成分是相

6、对独立于给定标记场X二x的，并且标记场中共有C个标记，这样的话，凝胶图像的分割问题可表示为：(1)其中P(X=x

7、Y=y)是标记场X=x在已知观测数据Y=y条件下的后验概率，P(yc

8、X=x)指的是在X=x的条件下观测场成分yc的概率分布，P(Y=y)则是已知观测数据Y=y的概率，在这里由于观测场Y二y指的是观测图像，是作为一个已知量存在的，因此在分析图像分割的最大后验概率过程可以忽略，P(Y=y)的影响。P(X=x)是X二x的先验概率，本文通过基于MRF分割模型MLL来获取标记场的先验概率，这里仅考虑二阶MLL模型，定义其能量函数如下：(2)其中V(xi,xj

9、)是势函数,当xi=xj,V(xi,xj)=-1；当xiHxj,V(xi,xj)=1；Ni是第i个位置的领域；？茁指的是空间位置函数[8],它的取值与隶属度函数有关，定义为？茁=1-0.5?滋i(xi),所以当隶属度？滋i(xi)值越小，说明其属于类c的程度越小，反之越大，也就是说在这里该点的分类很大程度上取决于它的领域状况。这样就得到了其先验概率：(3)通过前面的推导，仅有P(yc

10、X=x)作为未知量存在，由于一个类的特征数据通常被认为是正态分布的，在此我们将做进一步假设，令图像中属于同一类的像素服从高斯分布，那么在计算高斯分布的时候引入一个灰度密度的加权系数

11、来实现高斯分布灰度均值与