专利名称：：基于增强空域-变换域统计模型的纹理图像分割方法
技术领域：
：本发明涉及图像处理
技术领域：
，具体地说是一种纹理图像分割的方法，可用于对SAR图像的分割。技术背景一直以来，国内外学者提出了为数众多的纹理图像分割方法，主要有共生矩阵法、基于小波和Gabor滤波器的方法、马尔可夫随机场方法等。近年来，多尺度思想在图像分割领域得到了广泛应用，图像多尺度分解的优点在于，低分辨率的图像极大地降低了对图像处理的复杂度，同时为高分辨率图像的处理过程提供指导信息，从而大大减少图像处理对先验知识的依赖，从处理的图像特征角度看，多尺度分割方法可分为基于变换域和空间域两大类。基于变换域多尺度分割方法则先对图像做变换，如小波变换，然后分析变换域系数、建立合适的模型，最后运用建立的模型计算各个尺度上的初分割。基于空间域多尺度分割方法通常对图像下采样建立塔型结构，然后运用MRF(MarkovRandomField)方法对各个尺度上的特征建模，并得到多尺度的初分割结果。在基于变换域的多尺度统计图像建模方面，M.S.Crouse等人提出的小波域HMT模型能有效地描述小波系数在尺度间、尺度内的统计相关性，它是一种新的统计图像感知与识别方法。在图像分割领域，H.Choi等人提出了基于小波域HMT模型的多尺度图像分割方法HMTseg，实验表明HMTseg方法对均匀纹理分割性能良好，但对于非均匀纹理分割效果则有待提高。在基于空间域统计图像建模方面，Li和Gmy等人提出的2D-HMM(Two-DimensionalHiddenMarkovModel)对纹理有较强的描述能力，该模型先将纹理图像分割成大小相同的图像块，在考虑图像块内纹理特征的同时，更强调了对图像块间相关性的统计建模。2D-HMM捕捉了纹理图像宏观上的结构关系，应用到图像分割中效果良好。从捕捉纹理图像的宏观和微观特征出发，Lu和Carin等人将空间域的HMM模型和变换域的HMT模型相结合提出了HMM-HMT(HiddenMarkovModel-HiddenMarkovTree)模型。他们认为，HMTseg方法是基于图像块内的分割方法，该方法更侧重于考虑图像块内的统计特性，因此HMT模型不能准确地描述纹理图像的宏观统计规律性，尤其对结构型纹理图像HMT模型难以捕捉其纹理子结构间复杂的空间结构关系。为此，Lu和Carin等人对每一类纹理均联系一个HMM-HMT模型，然后通过HMM和HMT分别对图像块间的相关性和图像块内特征建模。HMM-HMT有效地描述了纹理图像块间的宏观和图像块内的微观统计特征，故该模型在一个尺度上得到了较可靠的初分割结果。然而，应用HMM-HMT模型计算各个尺度上的初分割时需由于要分别训练各个尺度上的HMM-HMT模型参数，因而极大地增加了计算复杂度；如果为了减小计算复杂度，只计算最粗尺度上的初分割而不计算其它各细尺度上的HMM-HMT模型参数，则会造成最终分割结果区域一致性差和边界保持不良的缺陷。
发明内容本发明的目的在于克服上述已有问题的缺陷，提出一种基于增强HMM-HMT和多状态加权HMT的纹理图像分割方法，以降低计算复杂度，提高图像分割的区域一致性和边界保持度。实现本发明目的的技术方案是在最粗尺度分割上，对现有HMM-HMT模型中对纹理宏观结构关系描述进行改进，提出增强HMM-HMT模型，以改善最粗尺度的初分割结果；在细尺度分割上，提出多状态加权隐马尔可夫树模型(Multi-StatesWeightedHiddenMarkovTree,MSWHMT)，以减少计算复杂度，提高各个细尺度上对区域间边界的感知能力，其具体实现步骤包括如下(1)输入纹理图像的类别数E，并将纹理图像初始分割成16X16大小互不重叠的粗图像块，对每一个图像块做二维离散小波变换；(2)对二维离散小波变换后的各个图像块内的系数，通过EM算法训练出各图像块的隐马尔可夫树模型皿丁参数//^7;1，/^，训代表平滑图像块，s代表奇异图像块；(3)对初始分割的各个图像块周围邻域的8个图像块进行标记，其中将中心块记为C,将与中心块直接相邻的4个块记为I类邻域块，将与中心块对角线方向的4个块记为I类邻域块；(4)采用类Viterbi算法训练各个初始分割图像块的块隐状态概率&=="、中心块c至l类邻域块的隐状态转移概率、，和中心块C至I类邻域块的隐状态转移概率A,"，构成增强HMM的参数，其中*={^力，代表平滑图像块，s代表奇异图像块；(5)输入HMT参数/fMr，,/fil化和增强HMM的参数;r,、《，、&，并利用这些参数计算各类纹理图像的似然"C|《)=wax{;rxa["x/〖^x(『》J"L(『。"rA，…,/t9*C^C^^I*C4^Ih^1、丑I"其中，^)和",(K皿)分别代表图像块《在I类、I类邻域块作用下它属于第f类纹理的似然，4(『3,)=^(『^1界"1,…，2，g表示第i个图像块的隐状态变量，K,为第i个图像块的小波系数，*={^^}，《={層丄，厨:，《，《，1}，,=1,2,...,￡;(6)从各类纹理的似然值中找出最大的一个似然值，用该最大似然所对应的纹理类别作为粗图像块的初分割结果；(7)依次将纹理图像分割成8X8、4X4和2X2大小互不重叠的细图像块，利用粗图像块上的得到的模型参数&、肌^1和//^1《，计算各细图像块上各节点/的似然值，并在各类纹理的似然值中找出最大的一个似然值，作为各细图像块的初分割结果，该似然值的计算公式如下、1=I"=《c)+《c其中，~1，2，...，￡，￡为纹理类别数，《为第f类纹理的增强HMM-HMT参数，似然和/;(旰)分别由参数i/MTl和i/MT;在四叉树上通过一次自下而上的扫描过程计算得到，其表达式为<formula>formulaseeoriginaldocumentpage8</formula>其中，T7表示由图像块内的小波系数构成的四叉树上的一个节点，~表示节点7/的块隐状态，^表示小波系数子树7/所包含的小波系数；(8)将粗尺度j上的初分割结果以1:4扩展，使之与细尺度j-l上的初分割的各个节点一一对应，然后判断粗尺度节点^是否为边界节点，确定出纹理图像的边界；(9)采用基于最大后验概率MAP的多尺度融合策略确定纹理图像边界节点的纹理类别，按照粗图像块的初分割结果确定纹理图像非边界节点的纹理类别，得到最终的纹理分割结果。本发明与现有的技术相比具有以下优点1、本发明提出的增强HMM-HMT模型对中心块与8邻域内的所有块建立直接联系，增强了模型对纹理空间结构的描述能力。实验表明，增强HMM-HMT模型在最粗尺度分割结果上得到了区域一致性好的初分割结果，相比HMM-HMT模型性能有所提高；2、本发明在计算各个细尺度上的初分割时，提出了多状态加权HMT模型。与增强HMM-HMT模型相比，多状态加权HMT模型只需要增强HMM-HMT模型在最粗尺度上的模型训练参数，因而多状态加权HMT模型极大地降低了时间复杂度；其次，多状态加权HMT模型舍弃了各个细尺度上块与邻域块间的联系，消除了块邻域成员对中心块的影响，不再通过描述块间的邻域关系来刻画该细尺度上的纹理空间结构信息，而采用基于块内的HMT模型来捕捉其细微特征，进而提高了该模型在各个细尺度上的边界保持度。图1是本发明的I类邻域块和I类邻域块结构图；图2是本发明的纹理分割方法流程图；图3是现有的由微纹理合成的2-6类测试纹理图；图4是用本发明与现有其它不同方法对图3微纹理的分割结果比较图；图5是现有由宏纹理合成的2-6类测试纹理图；图6是用本发明与现有其它不同方法对图5宏纹理的分割结果比较图；图7是现有由混合纹理合成的2-6类测试纹理图；图8是用本发明与现有其它不同方法对图7混合纹理的分割结果比较图。具体实施方式步骤l，输入纹理图像，以及纹理图像的类别数E。输入的纹理图像包括图3所示的由微纹理合成测试纹理图，其中图(a)是合成纹理图像mosaic7，图(b)是合成纹理图像mosaic8,图(c)是合成纹理图像mosaicl,图(d)是合成纹理图像mosaic4，图(e)是合成纹理图像mosaic9;图5所示的由宏纹理合成的测试纹理图，其中图(a)是合成纹理图像mosaic10，图(b)是合成纹理图像mosaicll,图(c)是合成纹理图像mosaic2，图(d)是合成纹理图像mosaic5，图(e)是合成纹理图像mosaic12;图7所示的由混合纹理合成的测试纹理图，其中图(a)是合成纹理图像mosaic13，图(b)是合成纹理图像mosaic14,图(c)是合成纹理图像mosaic3，图(d)是合成纹理图像mosaic6，图(e)是合成纹理图像mosaic15;这三种纹理图像的类别数均为2-6类；步骤2，对初始纹理图像分为16X16大小互不重叠的粗图像块，对每个图像块进行二维离散小波变换；步骤3,利用增强HMM-HMT模型计算粗尺度的初分割结果，具体步骤为(3.1)对二维离散小波变换后的各个图像块内的系数，通过EM算法训练出各图像块的隐马尔可夫树模型HMT参数ifMr，,/ZM7;，柳代表平滑图像块，^代表奇异图像块；(3.2)对初始分割的各个图像块周围邻域的8个图像块进行标记，其中将中心块记为C,将与中心块直接相邻的4个块记为I类邻域块，将与中心块对角线方向的4个块记为I类邻域块；其邻域结构图如图1所示；(3.3)采用类Viterbi算法训练各个初始分割图像块的块隐状态概率^=;7(5="、中心块C至I类邻域块的隐状态转移概率^和中心块C至I类邻域块的隐状态转移概率凡,,构成EHMM的参数，其中*={^^}，代表平滑图像块，^代表奇异图像块；(3.3.1)对于每一个图像块所联系的块隐状态变量，其初始值可任意赋值，并初始化增强HMM参数&，am,々m，*={w^}，附,"={儲力；(3.3.2)依据块隐状态变量的值对所有图像块分成i/iiT^和i/M7;2组，//M乙代表平滑图像块组HMT参数，ifM:代表奇异图像块组HMT参数，更新每一组的HMT参数，并计算每一个图像块^对块隐状态A的似然4(K,)，*,={^,4，/",2,…，2，2为图像块的总个数；(3.3.3)计算中心块C的后验概率^(yy和中心块C与其8邻域的联合状态概率^(^):<formula>formulaseeoriginaldocumentpage11</formula>(3.3.4)对(1)式求最大值，估计最可能的状态序列A，^，^,、,^,、，^，^，^L并对中心块C的隐状态变量赋值为化；(3.3.5)将每一个图像块均当作中心块C一次，并按步骤(3.3.3)和(3.3.4)估计出每个图像块的隐状态变量的值；(3.3.6)若参数收敛，则满足终止条件，训练停止；否则，转步骤(3.3.2);(3.4)输入1^丁参数/^7;1，/^^和增强1^1]\4的参数％、《，、&，并利用这些参数计算各类纹理图像的似然-<formula>formulaseeoriginaldocumentpage11</formula>其中,",)和《,(『Sl)分别代表图像块A在:类、n类邻域块作用下它属于第f类纹理的似然，"(f^^/K^IA"),、1,…，2，5,表示第i个图像块的隐状态变量，^为第i个图像块的小波系数，*={^^}，《={//^^丄,^^^7;，《,，~1，2,..，￡;(3.5)从各类纹理的似然值中找出最大的一个似然值，用该最大似然所对应的纹理类别作为粗图像块的初分割结果。步骤4，利用多状态加权HMT模型计算各细尺度的初分割结果，具体方法为依次将纹理图像分割成8X8、4X4和2X2大小互不重叠的细图像块，利用粗图像块上的得到的模型参数&、/fl《L和ffl《，计算各细图像块上各节点,'的似然值，并在各类纹理的似然值中找出最大的一个似然值，作为各细图像块的初分割结果，该似然值的计算公式如下其中，,=1，2,...，￡,￡为纹理类别数，《为第f类纹理的增强HMM-HMT参数,似然lw)和a^)分别由参数//M7:和//mt:在四叉树上通过一次自下而上的扫描过程计算得到,其表达式为、W)"(WI=力其中，T7表示由图像块内的小波系数构成的四叉树上的一个节点，~表示节点7y的块隐状态，(表示小波系数子树77所包含的小波系数。步骤5，将粗尺度j上的初分割结果以k4扩展，使之与细尺度j-l上的初分割的各个节点一一对应，然后判断粗尺度节点s是否为边界节点，确定出纹理图像的边界，具体步骤如下(5.1)若粗尺度y'上节点s的8邻域中存在与该节点s不同纹理类别的节点，则S为边界节点；(5.2)若粗尺度y'上节点s的纹理类别与其细尺度7—1上子节点f的纹理类别不一致，则还需要利用细尺度_/_2上的边界信息判断节点s是否为边界节点；(5.3)若同时满足如下两个子条件，则认为s是否为边界节点-第一个条件是对于粗尺度y上节点s的似然/,=(/:,/;，...,/，)，取/,的各分量中次小者为候选纹理，且候选纹理类别与细尺度上子节点^的纹理类别第二个条件是以细尺度>1上子节点f为中心划3x3大小的窗口，将该窗口联系一个窗口似然矢量/,>=(^''，...，/4^'',/5^'',f',…，C'),其中各分量，如表示节点f的8邻域窗口内序号为1的节点属于纹理类别/^^/,的似然，序号5对应窗口的中心子节点f，序号1-4和6-9对应子节点f的8邻域;同理,对粗尺度y做以上处理,并将以子节点s为中心的窗口联系一个窗口似然矢量^w-(r"',…,/rs/广'，/广s…，c')，则子节点纟的窗口似然矢量乙(。的方差比节点^的窗口似然矢量/,>的方差要大；即v"r"'，...w;，"，…，")>v"k/,v..,/:,/;,/6v..,")其中,v"K/;，…，/:，/;,/:,…,/;)定义为^(/;广.，/:,/;,/:,".，0=^1:(/;-")2j9'=i步骤6，采用基于最大后验概率MAP的多尺度融合策略确定纹理图像边界节点的纹理类别，按照粗图像块的初分割结果确定纹理图像非边界节点的纹理类别，得到最终的纹理分割结果；微纹理图像、宏纹理图像和混合纹理图像的分割结果分别如图4、图6、图8所示，其中(a)是分割的模板，(b)是HMT分割方法，(c)是结合边界信息的多尺度MAP融合HMT分割方法，(d)是多尺度MAP融合EHMM-HMT分割方法(e)是结合边界信息的多尺度MAP融合EHMM-HMT分割方法。本发明的效果可通过以下仿真实验进一步说明为了具体说明本发明的优势和特点，下面对该发明的实现过程进行详细介绍，并应用到纹理图像分割过程中去，得到的分割结果见附图。实验l:最粗尺度上的初分割结果实验1目的是测试并比较EHMM-HMT模型和HMM-HMT模型在最粗尺度上的初分割性能。我们将不同性质的纹理合成了3类测试纹理，分别进行了仿真实验。它们分别是均匀纹理、非均匀纹理及其混合纹理的和纹理。实验的各纹理均是Brodatz纹理图像库里的纹理样图，它包含112幅自然纹理图像，大小均为640x640，256级灰度。如图2所示，第一行由三幅4类测试纹理组成，其中mosaicl由均匀纹理D9、D75、D95、D68合成，mosaic2由非均匀纹理D42、D63、D91、D71合成，mosaic3由均匀非均匀混合纹理D95、D9、D68、D91合成;第二行由三幅5类测试纹理组成，其中mosaic4由均匀纹理D35、D21、D16、D18、D95合成，mosaic5由非均匀纹理D42、D63、D91、D71、D108合成，mosaic6由均匀非均匀混合纹理D75、D9、D68、D95、D42合成。表1给出了EHMM-HMT与HMM-HMT在最粗尺度上的初分割精度。表1EHMM-HMT与HMM-HMT在最粗尺度上初分割性能比较(正确率:％)模型均匀纹理非均匀纹理混合纹理mosaic1mosaic4mosaic2mosaic5mosaic3mosaic6画M-HMT95.0188.5787.9874.9093.1684.37EHMM-腿T95.7388.9689.1675.2092.2984.48由表1可见，对于大多数测试纹理，EHMM-HMT在最粗尺度的初分割性能比HMM-HMT要优。同时可见，由于EHMM-HMT对中心块与8邻域内的所有块建立直接联系，增强了模型对纹理空间结构的描述能力，进而提高了粗尺度的初分割性能。实验2:各细尺度的初分割结果本部分以实验1中的测试图像为例，分别运用EHMM-HMT、HMT、MSWHMT计算各细尺度上的初分割，然后给出了它们的边界检测精度和所耗时间。由于基于小波域的EHMM-HMT在最细尺度上剖分得到的图像块大小为2X2，此时基于块的HMT模型退化为IM模型，故在表4、表5、表6中未给出其结果，并将其边界检测精度和耗时分别标记为"X"和"/"。本实验的硬件测试平台是:Pentium4CPU，主频2.9GHz，内存l.OGB;软件平台为:WindowsXP操作系统和Matlab7.3.0。表2对均匀纹理在各尺度上的边界检测精度及耗时<table>tableseeoriginaldocumentpage14</column></row><table><table>tableseeoriginaldocumentpage15</column></row><table>表3对非均匀纹理在各尺度上的边界检测精度及耗时<table>tableseeoriginaldocumentpage15</column></row><table>表4对混合纹理在各尺度上的边界检测精度及耗时<table>tableseeoriginaldocumentpage15</column></row><table><table>tableseeoriginaldocumentpage16</column></row><table>表2、表3和表4分别给出了EHMM-HMT、HMT、MSWHMT在各细尺度上的初分割结果。结果表明，EHMM-HMT模型在更精细的尺度，比如8X8、4X4禾B2X2上重新估计EHMM-HMT参数不仅增加了算法的运行时间，而且得到的边界检测精度不尽如人意，在多数指标上逊于MSWHMT模型的结果。一方面，由于纹理特征通常以观察者所选取尺度的不同而变化，因此对图像分割而言，一幅图像针对某特定的应用总存在其最合理的尺度；另一方面，由于EHMM-HMT在细尺度上剖分的块随尺度变小而变小，而细尺度上较小的块却不能提供足够多的基本信息，于是基于较小的块计算出的似然可信度差，最终导致在此基础之上训练得到的EHMM参数并不能准确地描述纹理的宏观特征。分析MSWHMT模型我们发现，一方面，MSWHMT舍弃了图像块与邻域图像块间的联系，消除了块邻域成员对中心图像块的影响，不再通过描述图像块间的邻域关系来捕捉纹理空间结构信息，提高了MSWHMT对区域间边界的感知能力。另一方面，MSWHMT保留了描述不同纹理子结构的全局统计信息，实现了对各纹理子结构的宏观估计与局部描述的有机结合，在各细尺度上达到了减少区域边界误分的目的，因而在细尺度上对边界的检测比EHMM-HMT、HMT性能更优。权利要求1.一种基于增强空域-变换域统计模型的纹理图像分割方法，包括如下步骤(1)输入纹理图像的类别数E，并将纹理图像初始分割成16×16大小互不重叠的粗图像块，对每一个图像块做二维离散小波变换；(2)对二维离散小波变换后的各个图像块内的系数，通过EM算法训练出各图像块的隐马尔可夫树模型HMT参数HMTsm，HMTs，sm代表平滑图像块，s代表奇异图像块；(3)对初始分割的各个图像块周围邻域的8个图像块进行标记，其中将中心块记为C，将与中心块直接相邻的4个块记为I类邻域块，将与中心块对角线方向的4个块记为II类邻域块；(4)采用类Viterbi算法训练各个初始分割图像块的块隐状态概率πk＝p(B＝k)、中心块C至I类邻域块的隐状态转移概率αm，n和中心块C至II类邻域块的隐状态转移概率βm，n，构成增强HMM的参数，其中k＝{sm，s}，m，n＝{sm，s}，sm代表平滑图像块，s代表奇异图像块；(5)输入HMT参数HMTsm，HMTs和增强HMM的参数πk、αm，n、βm，n，并利用这些参数计算各类纹理图像的似然<mathsid="math0001"num="0001"><math><![CDATA[<mrow><mi>L</mi><mrow><mo>(</mo><mi>C</mi><mo>|</mo><msub><mi>θ</mi><mi>t</mi></msub><mo>)</mo></mrow><mo>=</mo><munder><mi>max</mi><mrow><msub><mi>k</mi><mn>1</mn></msub><mo>,</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><msub><mi>k</mi><mn>9</mn></msub></mrow></munder><mo>{</mo><msubsup><mi>π</mi><msub><mi>k</mi><mi>C</mi></msub><mi>t</mi></msubsup><mo>×</mo><msubsup><mi>α</mi><mrow><msub><mi>k</mi><mi>C</mi></msub><mo>→</mo><msub><mi>k</mi><mi>I</mi></msub></mrow><mi>t</mi></msubsup><mo>×</mo><msubsup><mi>β</mi><mrow><msub><mi>k</mi><mi>C</mi></msub><mo>→</mo><msub><mi>k</mi><mi>II</mi></msub></mrow><mi>t</mi></msubsup><mo>×</mo><msubsup><mi>l</mi><msub><mi>k</mi><mi>I</mi></msub><mi>t</mi></msubsup><mrow><mo>(</mo><msub><mi>W</mi><msub><mi>B</mi><mi>I</mi></msub></msub><mo>)</mo></mrow><mo>×</mo><msubsup><mi>l</mi><msub><mi>k</mi><mi>II</mi></msub><mi>t</mi></msubsup><mrow><mo>(</mo><msub><mi>W</mi><msub><mi>B</mi><mi>II</mi></msub></msub><mo>)</mo></mrow><mo>}</mo></mrow>]]></math></maths>其中，id="icf0002"file="A2009100233620002C2.tif"wi="12"he="4"top="178"left="45"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>和id="icf0003"file="A2009100233620002C3.tif"wi="12"he="4"top="179"left="65"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>分别代表图像块Bi在I类、II类邻域块作用下它属于第t类纹理的似然，<mathsid="math0002"num="0002"><math><![CDATA[<mrow><msub><mi>l</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>W</mi><msub><mi>B</mi><mi>i</mi></msub></msub><mo>)</mo></mrow><mo>=</mo><mi>p</mi><mrow><mo>(</mo><msub><mi>W</mi><msub><mi>B</mi><mi>i</mi></msub></msub><mo>|</mo><msub><mi>B</mi><mi>i</mi></msub><mo>=</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math>id="icf0004"file="A2009100233620002C4.tif"wi="49"he="5"top="187"left="82"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>i＝1，…，Q，Bi表示第i个图像块的隐状态变量，id="icf0005"file="A2009100233620002C5.tif"wi="4"he="4"top="197"left="99"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>为第i个图像块的小波系数，k＝{sm，s}，<mathsid="math0003"num="0003"><math><![CDATA[<mrow><msub><mi>θ</mi><mi>t</mi></msub><mo>=</mo><mo>{</mo><msubsup><mi>HMT</mi><mi>sm</mi><mi>t</mi></msubsup><mo>,</mo><msubsup><mi>HMT</mi><mi>s</mi><mi>t</mi></msubsup><mo>,</mo><msubsup><mi>π</mi><mi>k</mi><mi>t</mi></msubsup><mo>,</mo><msubsup><mi>α</mi><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow><mi>t</mi></msubsup><mo>,</mo><msubsup><mi>β</mi><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow><mi>t</mi></msubsup><mo>}</mo><mo>,</mo></mrow>]]></math>id="icf0006"file="A2009100233620002C6.tif"wi="58"he="4"top="205"left="48"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>t＝1，2，…，E；(6)从各类纹理的似然值中找出最大的一个似然值，用该最大似然所对应的纹理类别作为粗图像块的初分割结果；(7)依次将纹理图像分割成8×8、4×4和2×2大小互不重叠的细图像块，利用粗图像块上的得到的模型参数πk、HMTsmt和HMTst，计算各细图像块上各节点i的似然值，并在各类纹理的似然值中找出最大的一个似然值，作为各细图像块的初分割结果，该似然值的计算公式如下<mathsid="math0004"num="0004"><math><![CDATA[<mrow><msub><mi>l</mi><msub><mi>θ</mi><mi>t</mi></msub></msub><mrow><mo>(</mo><msubsup><mi>W</mi><mi>i</mi><mi>j</mi></msubsup><mo>|</mo><msub><mi>θ</mi><mi>t</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>W</mi><mi>i</mi><mi>j</mi></msubsup><mo>|</mo><msub><mi>θ</mi><mi>t</mi></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mi>π</mi><mi>sm</mi><mi>t</mi></msubsup><msubsup><mi>l</mi><mi>sm</mi><mi>t</mi></msubsup><mrow><mo>(</mo><msubsup><mi>W</mi><mi>i</mi><mi>j</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msubsup><mi>π</mi><mi>s</mi><mi>t</mi></msubsup><msubsup><mi>l</mi><mi>s</mi><mi>t</mi></msubsup><mrow><mo>(</mo><msubsup><mi>W</mi><mi>i</mi><mi>j</mi></msubsup><mo>)</mo></mrow></mrow>]]></math></maths>其中，t＝1，2，…，E，E为纹理类别数，θt为第t类纹理的增强HMM-HMT参数，似然lsmt(Wij)和lst(Wij)分别由参数HMTsmt和HMTst在四叉树上通过一次自下而上的扫描过程计算得到，其表达式为<mathsid="math0005"num="0005"><math><![CDATA[<mfencedopen='{'close=''><mtable><mtr><mtd><msubsup><mi>l</mi><mi>sm</mi><mi>t</mi></msubsup><mrow><mo>(</mo><msup><msub><mi>W</mi><mi>i</mi></msub><mi>j</mi></msup><mo>)</mo></mrow><mo>=</mo><mi>p</mi><mrow><mo>(</mo><msup><msub><mi>W</mi><mi>i</mi></msub><mi>j</mi></msup><mo>|</mo><msub><mi>k</mi><msubsup><mi>T</mi><mi>i</mi><mi>j</mi></msubsup></msub><mo>=</mo><mi>sm</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>l</mi><mi>s</mi><mi>t</mi></msubsup><mrow><mo>(</mo><msup><msub><mi>W</mi><mi>i</mi></msub><mi>j</mi></msup><mo>)</mo></mrow><mo>=</mo><mi>p</mi><mrow><mo>(</mo><msup><msub><mi>W</mi><mi>i</mi></msub><mi>j</mi></msup><mo>|</mo><msub><mi>k</mi><msubsup><mi>T</mi><mi>i</mi><mi>j</mi></msubsup></msub><mo>=</mo><mi>s</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math></maths>其中，Tij表示由图像块内的小波系数构成的四叉树上的一个节点，id="icf0009"file="A2009100233620003C3.tif"wi="4"he="5"top="124"left="33"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>表示节点Tij的块隐状态，Wij表示小波系数子树Tij所包含的小波系数；(8)将粗尺度j上的初分割结果以1∶4扩展，使之与细尺度j-1上的初分割的各个节点一一对应，然后判断粗尺度节点s是否为边界节点，确定出纹理图像的边界；(9)采用基于最大后验概率MAP的多尺度融合策略确定纹理图像边界节点的纹理类别，按照粗图像块的初分割结果确定纹理图像非边界节点的纹理类别，得到最终的纹理分割结果。2.根据权利要求1所述的纹理图像分割方法，其中步骤(4)的具体训练过程如下2a)对于每一个图像块所联系的块隐状态变量，其初始值可任意赋值，并初始化增强HMM参数&,am,y^，&=&附,4，》7，"={柳力；2b)依据块隐状态变量的值对所有图像块分成//m乙和//mt;2组，//m7^代表平滑图像块组HMT参数，/w《代表奇异图像块组HMT参数，更新每一组的HMT参数，并计算每一个图像块^,对块隐状态it的似然4(^,)，《="m,s}，"1，2，…，2，2为图像块的总个数；2c)计算中心块c的后验概率^(&)和中心块c与其8邻域的联合状态概率^(^):<formula>formulaseeoriginaldocumentpage4</formula>其中，<formula>formulaseeoriginaldocumentpage4</formula>2d)对(1)式求最大值，估计最可能的状态序列{、,*2,*3,*4,^,、,W~}，并对中心块c的隐状态变量赋值为k;2e)将每一个图像块均当作中心块C一次，并按步骤(2c)和(2d)估计出每个图像块的隐状态变量的值；2f)若参数收敛，则满足终止条件，训练停止；否则，转步骤(2b)。3.根据权利要求1所述的纹理图像分割方法，其中步骤(8)所述的判断粗尺度节点s是否为边界节点，确定出纹理图像的边界，具体实现过程如下3a)若粗尺度上节点s的8邻域中存在与该节点s不同纹理类别的节点，则s为边界节点；3b)若粗尺度7'上节点s的纹理类别与其细尺度户l上子节点f的纹理类别不一致，则还需利用细尺度户2上的边界信息判断s是否为边界节点；3c)若同时满足如下两个子条件，则认为s是否为边界节点第一个条件是对于粗尺度7上节点S的似然/,=(/乂,...,0,取/,的各分量中次小者为候选纹理，且候选纹理类别与细尺度上子节点f的纹理类别I^^一致；第二个条件是以细尺度>1上子节点^为中心划3x3大小的窗口，将该窗口联系一个窗口似然矢量c(c',…,cv广',c',…,c'),其中各分量，如r"'表示节点^的s邻域窗口内序号为i的节点属于纹理类别丄"^/,的似然，序号5对应窗口的中心子节点f,序号1-4和6-9对应子节点f的8邻域;同理，对粗尺度；做以上处理,并将以子节点s为中心的窗口联系一个窗口似然矢量Lw=(『'V..，d"''，fv..,C，则子节点f的窗口似然矢量/,)的方差比节点s的窗口似然矢量/^)的方差要大；即<formula>formulaseeoriginaldocumentpage5</formula>其中,v"k/;，…,/:,/;，/:,…,/;)定义为<formula>formulaseeoriginaldocumentpage5</formula>全文摘要本发明公开了一种基于增强空域-变换域统计模型的纹理图像分割方法，它属于图像处理
技术领域：
。主要解决现有纹理图像分割方法分割准确性差，计算复杂度高的问题。其分割步骤为(1)输入纹理图像和纹理类别数；(2)对纹理图像分成16×16的粗图像块，对每个图像块进行二维离散小波变换；(3)训练EHMM-HMT参数，计算粗图像块的似然值及粗尺度上的分割结果；(4)对纹理图像进行分成8×8、4×4和2×2的细图像块；(5)计算粗图像块的似然值及各细尺度上的分割结果；(6)结合边界信息的多尺度MAP融合。本发明具有纹理图像分割准确度高，计算复杂度低的优点，可用于微纹理和宏纹理图像的分割。文档编号G06T7/40GK101609557SQ20091002336公开日2009年12月23日申请日期2009年7月17日优先权日2009年7月17日发明者彪侯,芳刘,焦李成,爽王,缑水平,郝红侠,桦钟,陈蓉伟申请人:西安电子科技大学