小波变换在图像分割中的应用

存档日期：___＿_＿_______________存档编号:＿＿＿__________＿_____＿_＿徐州师范大学科文学院本科生毕业论文(设计)论文题目：小波变换在图像分割中的应用姓:王婷婷学号：078333106专业：自动化班级:0７自动化指导老师:李旭超科文学院教务处印刷摘要近年来,对图像分割的研究一直是图像技术研究的焦点。图像分割是一种很重要的图像分析技术,它的目的是把图像分为具有各种特性的区域并把感兴趣的部分提取出来。它融合了多个学科的成果,并且成功应用于工业、农业、医学、军事等领域，得到了广泛的应用。图像分割的方法有很多种,本文的题目是小波变换在图像分割中的应用。主要介绍了两种小波分割方法，即小波阈值分割方法和小波域马尔可夫随机场模型分割方法,其中通过了直方图、建立模型等手段对这两种方法做出具体的讨论。经过研究发现，区分真正的噪声和边缘是图像分割的难题之一，然而小波变换则可以解决这一问题，小波变换是一种时－频两域分析的工具。同时还利用Matlab分别对两种方法进行仿真,并得到了有效的结果。根据仿真结果我们可以看出不同分割方法的不同分割效果,从而更好地理解这些方法。图像分割是一个经典的问题，至今仍没有一种通用的解决方法。关键词:图像分割小波变换阈值马尔可夫随机场模型AbｓtractＲecentｌy，tｈeimagｅｓegmenｔａtiｏnissｔillａfocｕsｏfimaｇeprｏcessing.Itiｓasoiｍpｏrtanttechnolｏgyinthefieldofimageproceｓｓing，ｔhｅａiｍofwhicｈseparaｔesｔｈｅobjｅｃtintoｔheｒegｉonsｗithdiｆｆｅrentｃharacｔerｓ,ａnｄpickupｔheｉｎｔereｓtｉngregｉons.Iｔincｏrporａtｅsmultipｌeｄiｓciplｉｎｅｓ,itappliｅdintｏindustｒｙ,agriｃultｕrｅ，mｅdicｉｎe,miｌiｔarｙandotherfiｅｌdsｓuccessfully,ａｎdbeｅｎwｉｄｅlyused.Tｈｅrｅarｅaｌｏtofkiｎdsｏfiｍaｇesｅgmentatioｎ,ｔｈetopicoｆthｉspaｐerisＴheaｐpliｃatioｎofｗaveleｔtｒaｎsｆｏrmtotheimaｇeｓeｇmentａtｉoｎ.Itmaｉnｌｙiｎtroｄｕcｅstwokindsofwavｅｌｅtｄｅcｏｍpositiｏｎmethｏds,tｈeyaｒeｒespectivｅｌythewａveletthｒｅｓholdsｅｇｍｅntaｔionmetｈoｄａndthｅｗaveｌｅｔｄｏmaｉnｍaｒkovmodelsegmentａtｉｏnｍｅｔhｏｄwｉthｔheａirｐort，Ａmongthｅmｔｈrｏｕghthｅhｉstogｒam,mｏdｅlａｎddisｃussthemｅaｎsｏfthetｗoｍetｈoｄｓ.Ａftｅｒtｈｅstudyｆoｕnd,ｄisｔinguiｓhｔruenoiseanｄeｄgｅisｏneｏftheｄifficultproblemｓoｆtｈｅimａgeｓegmenｔatiｏn,howevｅrｗａveletｔｒａnsfoｒmcanresｏlvethisproblembｅtｔeｒ，becauｓewaｖeｌettransforｍisakｉｎｄoｆtiｍe－frｅqueｎｃｙaｎalyｓisｔoｏls．ＴｈｅntｈｅrearetwosummａtiｏnsfoｒｔwomｅthodswithMａtlab，tｈeresuｌtcａnshoｗtｈｅefｆeｃｔivｅphenomenon.Ａccordｉngtｏｔhesimuｌationｒesult，ｗecanseetｈｅdiffｅrｅntimagｅsｅgｍentatioｎeｆfectofｄiffｅｒｅntmethods,soweｃａnundｅrsｔandｔhesewａｙsofｉmagesegmentatioｎbｅｔｔｅｒ.Ｉｍagesegｍｅｎｔaｔioniｓaｃｌassｉcpｒoblem，ｂutｓtillcannoｔfinｄaｇｅnｅｒalsolutiｏｎtoit.Keｙwｏrd:imageｓegｍｅntationwaveｌｅttrａｎsfｏrｍtｈrｅｓhoｌdｗaｖeletｄomainmａrkovmｏdelwiｔhtheaｉrporｔ目录TOC\ｏ"１-3"＼h\z\uHYPＥＲLIＮK\ｌ"＿Toc29396１3０２"摘要ﻩPAGEREＦ_Ｔoc2９39613０２\ｈIHYPＥRLＩNK＼l"＿Toc2９396１30３"ＡbstｒａctﻩPAＧEREＦ_Toc2９396130３\hIＩＨYPERLINK\l＂＿Tｏｃ29396130４"1绪论ﻩＰＡGEREF_Toc2９3９6130４\h1HYPERＬＩNK1.1.1串行边界分割技术ﻩPAGEREF_Toｃ293９６1306＼h1HYPＥRLINK\l＂＿Ｔoｃ293９6130７"１.1.２串行区域分割技术 ＰAGEREＦ_Tｏc2９３961３0７＼h1HYPEＲLINK\l"_Ｔoc２93961308"１.１.3并行边界分割技术ﻩＰＡGＥRＥF_Toｃ２939６13０８\h1HＹPEＲLINK１．1．4并行区域分割技术ﻩPＡGEREF_Tｏｃ２939６１3０9＼h２HYPEＲLIＮK\ｌ"＿Ｔoc29３961310"１．1.5结合特定理论工具的分割技术 PAGEＲEF_Toc293961３10\h２HYＰＥRＬINK1.2频域图像分割 PＡGＥREＦ_Tｏc2９39６１3１１\h2HYPERＬINＫ１．3小波域图像分割 PAＧEREF_Ｔoｃ293961３12\ｈ3HYPERLINK＼l"＿Ｔoｃ2939６１３13"1.3.１图像分割的描述 PAGEREＦ_Ｔoc２9３9６1３1３\h3HYPEＲＬＩNK＼ｌ"_Ｔｏc29３9613１４"1．3.2图像分割的发展及现状ﻩPＡＧＥREF＿Toc２93961314\h5HYPEＲLINK＼l"_Toｃ２９３9６13１5"1.3.3基于小波变换的图像分割方法ﻩＰAGERＥＦ_Toc2９3９６1315＼h6ＨYPERLＩＮK\l"_Toc2939６13１６"1.４本文的主要工作ﻩＰＡGERＥF＿Toc293９6131６＼h7HＹＰERLINＫ２.1小波理论 PＡGＥREF_Toｃ293９６1318＼h8HYPERLＩNK2．2．1小波变换的概述ﻩPAGEREF_Tｏｃ29396１320\h8HＹPERLINＫ\l＂_Toc2９３9６1321＂2.2.2正交小波基的种类ﻩＰAGＥREＦ＿Ｔｏｃ293９61321\h9ＨYPERLＩNK＼l"_Tｏc２93９61322"2．2.3多分辨率分析ﻩＰAGEＲＥF_Ｔoc2９396132２＼h１0HＹPEＲLＩＮK２.２．4连续小波变换ﻩPAGEREF_Tｏｃ２９39６１323\ｈ11HYＰＥRＬINK2.2.5离散小波变换 PAGERＥF_Toｃ293961324\h1２ＨＹＰERLINＫ\ｌ"_Toc2９３9613２5"2.３马尔可夫随机场的基本理论 ＰAGEREF＿Toc293９61325＼ｈ１3HYPERLIＮK3图像分割中的小波阈值法ﻩPAGＥREF_Toc2９3９6１３26\h1５ＨYPERＬINＫ＼l＂＿Toc293961327"３.1小波阈值法的原理ﻩPAＧEＲEF＿Toc2939６1327＼h15HＹＰEＲLINK＼ｌ"_Ｔoc293961328"3.２图像直方图的多分辨率分析 PＡGEREF＿Ｔｏc293９6１32８＼h16HＹPERLINＫ３.4实验仿真ﻩPAGERＥF_Toｃ２9３961330\h1７ＨYＰERLINK\ｌ"_Tｏc2９396133１"４图像分割中的小波域马尔可夫随机场方法ﻩPAGEＲEF_Toｃ293９6１３31\h1９HYPERＬIＮK\l＂_Ｔoｃ293961332＂４.1图像分割中的马尔可夫随机场方法综述ﻩPAGEＲEＦ＿Ｔoc29３96１332\h19HYPＥRＬＩNＫ＼l＂_Ｔｏc２93９６１3３3＂４．1.1用马尔可夫随机场描述图像模型ﻩPAGEREF_Ｔｏc29396133３＼h１９HＹPＥRLINＫ\l"_Tｏｃ2９３9６1３34"４．1.２基于适当最优准则实现图像的分割 PＡGＥRＥＦ_Ｔoc2939６13３4＼h２1HＹPＥRLＩＮK4.２小波域马尔可夫随机场模型的分割算法 PＡGＥＲＥF＿Toc２９３961３3５\ｈ2２ＨYPＥRLIＮK\ｌ＂＿Ｔｏc29３９613３6"4．２.1小波域马尔可夫随机场模型的MAP准则的分割算法 PAGERＥＦ_Toｃ2９39６１3３6＼h22ＨYＰＥRLIＮK＼l"＿Ｔoc２９396133７"4.2．2小波域马尔可夫随机场模型多尺度概率值的计算ﻩＰＡＧＥRＥF_Toc２９3９613３7\h23HYＰERLIＮK４.2.3小波域马尔可夫随机场模型分割算法的具体实现ﻩPAGＥRＥＦ_Toc2939６133８＼h23HＹPERＬＩNK\ｌ"_Toc２9３96１339＂4.3小波域马尔可夫随机场模型的参数估计 PAＧEREF_Toc2939613３9\h24HYPＥＲLINK\l＂＿Toｃ293９613４0"HYＰＥRLＩNK\ｌ＂_Toｃ293９61３4１"４.4实验仿真ﻩＰAGＥREF_Ｔoc29396１３41\h2５ＨYPERLＩNＫ\l＂_Toｃ29396１3４2"5结论与展望ﻩPAGEREＦ_Tｏｃ29３961３42\ｈ2６HYPＥRLINK５.1论文的总结ﻩＰAGEREＦ_Tｏc29396１３43＼h26HYPＥRＬINK\ｌ"_Toc293９61344"５.2论文的展望ﻩPＡGERＥＦ_Ｔｏｃ２９3961344\ｈ2６HYPERＬINＫ致谢 PＡＧEREF＿Tｏc2939６1３４5\h28HYPＥRLINK参考文献ﻩPＡGEREF_Toc2939６１346\h29ＨYPＥＲＬINK\ｌ"_Ｔｏｃ2939６１3４6"附录ﻩ３1１绪论1.1空域图像分割空域是指图像平面本身，空域图像分割就是直接对图像的像素进行处理分割。研究者经过几十年的研究与努力,研究出了很多种空域图像分割方法。归纳起来大致包括:串行边界分割技术、串行区域分割技术、并行边界分割技术、并行区域分割技术、结合特定理论工具的分割技术等[1][1]LIYue-e,LIUQing-fang.TheApplicationofWaveletTransformtotheImageSegmentation[J].CollegeofPhysicsandElectronicsEngineering,ShanxiUniversity,Taiyuan030006,China,2009,32(4):566~571.1.1.1串行边界分割技术串行边界分割技术指通过顺序搜索边缘点，采用串行方式来对感兴趣目标的边界进行检测。主要有以下三个关键步骤：①首先确定一个顺序搜索的起始边界点；②然后在确定先前的搜索结果对下一边界点的检测和下一个结果的影响的前提之下，选择某种搜索策略,根据相应原则逐一检测新的边界点；③最后选定搜索终止的条件，结束整个搜索过程。串行边界分割技术所采取的策略主要有以下两种：①首先检测出边界点,然后再连接边界点；②以交叉结合的方式来进行边界点的检查和连接。１.1.２串行区域分割技术串行区域分割技术指通过对目标区域的直接检测，用串行方式来进行图像分割的技术。它的特点是将整个处理过程分解为的多个步骤依次进行,然后前续步骤的处理结果来决定对后续步骤的处理。结合了特定数学理论工具的一些图像分割方法也经常用串行区域分割的方式。串行区域分割技术有两种基本形式：①从单个像素出发,逐渐合并成所需的分割区域;②从全图出发，逐渐分裂成所要的分割区域。1.1.３并行边界分割技术并行边界分割技术指通过对感兴趣区域的边界进行检测，用并行方式来对图像进行分割技术。其过程主要有以下两个步骤:①检测感兴趣区域的边界点；②形成感兴趣区域的边界。对于步骤①可以采用各种微分算子来直接检测,也可以利用拟合方法与边缘模型进行间接检测。对于步骤②由于其过程较复杂,因此单纯利用微分算子不能形成闭合边界,需要通过一定的准则和数学工具将感兴趣的区域分离出来。1．1．4并行区域分割技术并行区域分割技术指通过对感兴趣区域的检测，用并行方式来进行图像分割的技术。在实际应用中,并行区域分割技术主要包括以下两大类：①特征空间聚类方法；②阈值化方法。1．１.5结合特定理论工具的分割技术①基于信息论的分割技术最近几年，出现了许多借助信息论中熵的概念的图像分割方法。这些方法利用信息论当中求熵的极值的方式来进行图像分割。例如:1D最大熵法、2Ｄ最大熵法、最大后验熵法、最小熵相关法、最大香农熵法、条件熵法等等。②基于小波分析和小波变换的分割技术小波变换是空间(时间）和频率的局域变换,通过平移、伸缩等运算对函数和信号进行多尺度的细化分析，有效的从信号中提取信息，从而解决了傅立叶变换不能解决的许多问题。近年来，在低频和高频分析时,有“变焦”特性的小波变换在图像分割中得到广泛应用。1.２频域图像分割频域图像分割法就是傅里叶变换的方法,也是一种最简单的图象分割的方法。傅里叶变换一直是信号处理领域中最完美、应用最广泛、效果最好的一种分析手段，只是傅里叶变换是一种单纯对频域的分析方法,它在频域的定位性上是最准确的（即频域分辨率是最高），而在时域没有任何分辨能力，也就是表示傅里叶变换所反映的是整个信号在所有时间下的整体频域的特征,无法提供任何一段局部时间段上的频率信息。在现实生活中,经常会出现一些非平稳信号，例如音乐和语音信号等，它们的频域特征都在随着时间的变化而变化，这些信号被称为时变信号。Ｇａbor为了研究信号在局部时间范围内的频域特征,而提出了非常著名的Gaｂor变换,后来随着不断研究发展为短时傅里叶变换（又称加窗傅里叶变换,简称SＴFT）。如今短时傅里叶变换在许多领域已经得到了广泛的应用，但是它的定义决定了其窗函数的形状和大小都与频率和时间没有关系而保持固定不变,这样在分析时变信号时是不方便的[2]章毓晋著.图像分割[M].北京:科学出版社,[2]章毓晋著.图像分割[M].北京:科学出版社,2001,78-90.[3]李世雄编著.小波理论及其应用[M].南开大学数学研究所计算教学年印,1992.小波变换继承并且发展了短时傅里叶变换的局部思想的优点，同时还克服了其窗口大小不能随着频率变化而变化,且缺少离散正交基的缺点，是一种比较理想的用来处理信号的数学工具。因此目前小波变换在图象分割中获得了广泛的应用。１.3小波域图像分割1．3.1图像分割的描述图像分割（iｍａｇｅsegmentaｔion)是一种非常重要的图像技术，它是图像分析与识别过程中的重要环节,其分割的结果直接关系到图像处理后期其他工作的质量，例如图像的测量、图像的识别等。一方面它是图像表达的基础，对于特征测量有着至关重要的作用;另一方面，由于图像分割及其基于分割的特征提取、目标表达和参数测量等将原始的图像形式转化为更加抽象紧凑的形式，从而使得更高层的图像分析和理解成为可能[4]彭玉华著.小波变换与工程应用[M].科学出版社,2000.。图像分割在不同的领域中都有重要的应用,所以有时会用不同的名称,比如目标识别(targetiｄentifｉcaｔioｎ）技术,阈值化(threshｏld)技术,目标检测(taｒｇetｄetection）技术，目标轮廓(objｅctdelineate）技术，图像区分或求差（ｉmagｅdiscrimｉnａtiｏn）技术等,虽然这些名称看起来是不同的技术,但其实它们本身或核心都是图像分割技术，由此可以看出图像分割技术的重要性[4]彭玉华著.小波变换与工程应用[M].科学出版社,2000.[5]崔屹.图象处理与分析-数学形态学方法及其应用[M].北京:科学出版社,2000.图1-1图像分割在图像工程中的位置其实图像分割简单地说,就是要把一幅数字图像分割成几个不同的区域，其中在同一区域内的具有在一定准则下的图像可以认为是有相同的颜色、纹理和灰度等性质，同时在任意相邻区域间的图像的性质是有显著的区别。在长时间的研究中，对于图像人们提出了很多不同的解释和表达,图像分割可以借助集合概念给出如下的定义:令整个图像的区域用集合表示,对于的分割可以看作是将分割成个满足下面五个条件的非空子集(子区域)：（1);(2)对所有的和,,都有;(3）对,有；(4)对,有;（5）对，是连通的区域。其中对于全部在集合中的元素,是一致性逻辑谓词,代表空集。在对图像进行研究和应用中,人们经常只是对图像中的某些部分或者是某些区域感兴趣,这些部分或者区域常被称为前景或目标（其他部分被称为背景),他们一般是对应图像中所特定的、具有独立性质的区域。为了辨识和分析目标，需要将他们从图像中分离并提取出来,只有这样才有可能对目标进一步的分析利用。概括起来说,图像分割就是指把图像分成各具特性的区域并提取感兴趣的目标，从而对其进行分析、应用的技术和过程。这里特性是像素的颜色、灰度、纹理等，预先定义的目标可以是对应的单个区域，同时也可以是对应的多个区域[6][6]陈武凡.小波分析及其在图像处理中的应用[M].北京:科学出版社,2002.1．3.2图像分割的发展及现状对于图像分割的研究已经有几十年的历史，人们对此一直都高度重视。早在50年代中期即计算机视觉理论体系形成以前,人们就已经开始了对图像分割的研究。能够找到一种通用的、普适的图像分割方法是人们几十年来不断追求的梦想，人们为此付出了许多艰辛的努力,但也取得了不少研究成果,提出了很多图像分割算法。然而随着研究的不断深入,人们逐渐意识到,在将三维景物投影为二维图像的过程中,丢失了深度以及不可见部分的信息，不同的视角下的同一物体的图像会有很大的不同，还有会因为前后物体的遮挡而丢失信息等;另外,在场景中的很多比如物体表面几何、光源、物理特性以及成像设备与物体之间的空间特性等不同的因素,都被综合成为单一的图像中像元的灰度值;还有在成像过程中会或多或少地引入一些噪声和畸变。这些问题都直接或间接地导致了图像分割问题是一种信息不足的不确定性问题,所以不可能存在一种通用的对任何图像都适用的、统一的图像分割方法[7][7]DonohoDL.WaveletShrinkageandW.V.D.ATen-minuteTour.ProgressinWavelet[J].AnalysisandApplication.1993:109-128.在过去的四十多年里,人们一直在高度重视图像分割的研究,至今为止，研究者提出了上千种不同类型的分割算法，如：小波分析法、水线法、匹配法、马尔可夫随机场模型法等，并且近年来每年都有上百篇相关研究成果的发表。但是,如今的方法大多是为了特定应用而设计的,具有很大的局限性和针对性，对图像分割的研究还是缺乏一个统一的理论体系。近年来,随着计算机技术的不断发展和成熟，现在的图像分割已经得到了广泛的应用,几乎出现在所有有关图像处理的领域,例如在文档处理，工业自动化,生产过程控制，在线产品检验,遥感和生物医学图像分析,保安监视，以及军事,农业工程,体育等诸多领域。然而在缺乏足够先验信息的前提下进行图像分割是一项很困难的任务，如医学上人脑组织的细微变化、监控系统的随机性等，人们为了克服这一困难，建立了大量的图像模型来完成分割的任务，如利用概率分布函数刻画图像的纹理特征而建立起来的统计图像模型等。相应地也研究出了各种分割算法，如基于区域的分割、合并分割算法等。基于图像模型的分割在易于实现的同时,准确的图像模型还可以提高分割的质量,因此对图像模型的研究，不仅是各种应用的需要，而且也是分割取得突破性进展的关键问题[8][8]陈武凡,小波分析及其在图像处理中的应用[M].北京:科学出版社,2002.1.3.３基于小波变换的图像分割方法小波变换是一种近年来得到广泛应用的数学工具，是空间（时间)和频率的局域变换，因而能有效地从信号中提取信息。小波变换从图像分割的角度看,有以下几个优点[9][9]LORENZThomas.Set_ValuedMapsforimageSegmentation[J].ComputingandVisualizationScience,2001(5):45-57.(1）通过选取合适的滤波器,小波变换可以最大程度地减小或消除所提取信息的不同特征之间的相关性;（２)小波变换实现上有快速算法(即Mallａｔ小波分解算法);(３)小波变换具有“变焦”特性，在高频段可以用低频率分辨率和高时间分辨率,在低频段可以用高频率分辨率和低时间分辨率;(４)小波分解可以覆盖整个频域，从而在数学上可以提供一个完备的描述。小波理论在图像处理方面的应用日益成熟，小波变换具有方向性、多分辨率特性、非冗余性、以及小波系数具有重拖尾和持续性的非高斯分布等特性，对于刻画图像的非平稳性,小波理论提供了有效的工具。通过利用小波变换的多分辨率特性,让特征场建立在一系列小波域上，使相应分辨率的标号场利用该尺度的特征场,这样有利于刻画图像的非平稳性,小波系数随着分辨率的不同而具有不同的特征向量，从而有利于反映特征场的本质特征。同时对每一小波系数引入有限的隐状态，来反映小波系数的非高斯分布的特性。描述小波系数的非高斯分布需要用两个状态的高斯混合分布,纹理的小波系数、边缘用较大状态的方差来表示，平滑区域的小波系数用较小状态的方差来表示。大多数真实图像都可以使用这种模型,但该模型不具有通用性,如背景为较平滑的区域，受光照影响的交通图像分割问题,物体有纹理则相应的小波系数较大，但阴影的小波系数介于背景与物体之间，所以很难用两个状态的小波系数来表示图像的特征[10][10]BesagJ.Spatialinteractionandstatisticalanalysisoflatticsystems[J].Journalofroyalstatisticalsociety,1974,serialB(6)：192-236.与小波系数服从非高斯分布与高斯分布相背离来比较,小波系数通过有限通用混合分布（FGM)来逼近小波系数的分布似乎更加合理。同时在标号场先验模型确定的基础上,HMT(即ＨiddeｎMarkovTree)模型没有具体表达式来给出标号场先验概率分布，认为父节点周围3６个小波系数的标号与标号有关，这与小波系数的紧支撑性不同并且与小波系数的局部相关性也相互矛盾。另外，认为小波系数的标号是相互独立的的说法与尺度间小波系数的持续性、紧支撑性相背离，由于父节点与子节点一般都具有相同的标号，而且父节点周围的三个节点小波系数的标号与子节点的标号有关联。标号场体现的不是独立的作用关系而是小波系数标号的相互关系。１.4本文的主要工作图像分割是数字图像处理中的重要技术之一，是图像处理的核心技术。但是图像模型的准确建立是影响图像分割质量好坏的重要因素。本文对小波变换在图像分割中的应用进行了研究,主要内容有：第一章首先介绍了空域图像分割、频域图像分割和小波域图像分割。详细地叙述了图像分割技术的研究背景，包含了图像分割在国内外的发展及现状和小波变换在图像分割中的方法,为全文研究的小波变换在图像分割中的应用提供理论基础。第二章首先详细阐述了小波理论、小波变换的理论和马尔可夫随机场理论。为第三章和第四章图像分割方法的提出奠定了理论基础。第三章首先介绍了小波阈值分割方法的原理,然后阐述了阈值分割的算法，最后分别对小波阈值分割方法和传统阈值分割法进行仿真实验，通过对比表现出小波阈值分割的优点,这也是本文的创新点。第四章首先介绍了小波域马尔可夫随机场模型的分割方法的综述、算法以及参数的估计,最后分别用kmeaｎs分割方法和本文的分割方法进行仿真实验，通过对比显现出本文分割方法的优势。第五章通过前四章对小波变换在图像分割中的应用的详细阐述，进一步进行总结,指出小波变换在图像分割中应用的广泛,但随着科技的日益发展,图像分割的技术与理论还需进一步的完善。2小波变换与马尔可夫随机场理论２．1小波理论小波(wａｖｅlet),是在有限的时间范围内变化,并且它的平均值为零的数学函数。小波函数的精确定义为：设为一平方可积函数,即,如果其傅里叶变换满足条件:（２-1）其中称为一个小波母函数或基本小波，称上面的公式为小波函数的可容许性条件。从小波的定义可以知道小波的特点是小,即它们在时域都具有紧支集或近似紧支集;还具有正负交替的波动性。2.2小波变换2．2.1小波变换的概述小波变换的基本思想是:先找到一个满足一定条件下的基本小波函数，通过对其进行伸缩和平移来构成小波函数族,其中每个在时域和频域上都具有比较好的局部性,再利用这个小波函数去表示和逼近所需研究的信号,这样有利于更好的分析和处理[11]王玲.多小波理论及其在图像处理中的应用研究[D][11]王玲.多小波理论及其在图像处理中的应用研究[D].西安电子科技大学博士学位论文,2000.7.一般我们把选取的基本小波函数用来表示,又称为母小波函数，而所选取的母小波函数又必须满足如下条件:(2－2)式(2-２）中的是的傅里叶变换。在满足允许条件的前提下,所选取的母小波函数都可看作是一个带通滤波器的脉冲响应．由上式可直接得出：（2-3）上式(2-2)和(2－３)是完全等价的，式(2-３)表明具有一定的震荡性。对母小波函数进行伸缩和平移以后便可得到如下小波函数族:（２-４)这里的称为小波函数(简称小波）。，,ａ、ｂ分别称为伸缩和平移因子，变量a反应函数的宽度或尺度，变量b指明函数沿着x轴的平移位置。一般情况下,母小波函数能量集中在原点,而小波函数能量集中在b点[8］。2．2.2正交小波基的种类（1)Haａr小波Ｈａar系的函数是不连续的,并且它在频域上随着的衰减速度只为，所以频域的局部性能不好,但是它的结构简单，经常运用到理论研究中。图2－1Haar小波函数图像(2)Ｍｏrlet小波Morｌet小波是高斯下的单频率复正弦函数:式中，表示虚数，常数。虽然Moｒｌet小波有解析表达式，但其不具有正交性的同时也不存在紧支集。Moｒleｔ小波的特点是能够提取信号中的幅值和相应的信息，并广泛应用于地球物理信号处理中。图2-３Morｌet小波函数图像(3)Ｍｅyer小波该小波在频域上具有紧支集和任意阶正则性，所以它在时域和频域都具有较好的局部性。图２-3Meyｅr小波２．2.3多分辨率分析Mａllａt在１9８7年将计算机视觉领域中的多尺度分析思想引入到小波分析中，并且提出了多分辨率分析概念。同时小波变换具有多分辨率性,多分辨率分析指满足下列性质的一系列闭子空间：（1）平移不变性:即,对所有的；(2）一致单调性:即；(３)伸缩规则性:即；（4)渐进完全性:即;（5）正交基存在性：即存在，使是的正交基，使,。由于离散化小波的信息量仍是冗余的,因此再次从数字计算机处理的角度考虑，人们仍然希望减小离散化小波的冗余量,直到得到一组正交基。这组正交基称为正交小波基。如何构成正交基,构造小波母函数，而解决这些问题的方法就是多分辨率分析理论。从多分辨率的定义可知,所有的闭子空间都由同一尺度函数伸缩后的平移系列展开成的尺度空间,被称为多分辨率分析中的尺度函数。多分辨率分析的一系列尺度空间是由同一尺度函数在不同尺度下展开的。所以多分辨率分析理论为小波变换提供了数学上的理论基础[12]付忠良.基于图像差距度量的阈值选取方法[J].北京:计算机研究与发展,2001,38(5):56-567.２[12]付忠良.基于图像差距度量的阈值选取方法[J].北京:计算机研究与发展,2001,38(5):56-567.２.2.4连续小波变换将任何空间里的函数在小波基下进行展开,我们称这种展开为函数的连续小波变换(简称为ＣWT),它的表达式为（2-5)从ＣＷＴ的定义可知，小波变换和傅里叶变换一样都是积分变换，被称为小波变换系数，因为小波基与傅里叶基不相同，所以小波变换和傅里叶变换的不同之处有很多。其中最重要的不同是，小波基有平移、尺度a两个参数,因此在小波基下把函数展开，就意味着将该时间函数投影到二维时间-尺度的相平面上。由于小波基本身的特点，把函数投影到小波变换域之后,可以更好的提取函数的一些本质特征。位移和尺度都连续变化的连续小波基函数与傅里叶基不同，形成了一组非正交的过度完全基。这表示其任意函数的小波展开系数间都有相关关系，如果用来描述两个基函数和的相关度大小，那么（２-6）其中表示时移半平面（由于a>0所以称半平面）、连续尺度的两个不同点间的CWT系数的相关关系,同时也称为再生核或者重建核,小波的选取决定了它的结构。CWT的缺点是系数具有较大的冗余量,但我们也可以利用它的冗余量来实现去噪和数据恢复[13]靳宏磊.多方法的图像分割和识别[D].[13]靳宏磊.多方法的图像分割和识别[D].东南大学博士学位论文,2000.2．2.5离散小波变换在实际应用中，尤其在计算机实现时,经常需要将连续小波离散化。对变换进行二进制离散是一种方便的形式,我们将离散化后的小波和所相应的小波变换分别称为二进小波和二进小波变换。与连续小波变换不同，因为离散小波变换在尺度-位移的相平面上，所对应是离散的点,所以被称为离散小波变换,把小波变换的连续相平面离散化。离散化小波变换表示函数的所有信息。任意函数都可表示为基本单元的加权和。在理想状态下,当离散后的小波基函数满足正交完备性条件时,所计算的小波变换系数没有任何冗余度，可以最大限度地压缩并减少计算量。可是在连续小波变换的概念出现以后，人们一直无法找到具有一定正则性的可作为空间的标准正交基的。所以人们就去研究具有一定冗余度、非正交的离散小波基离散栅格的取法和它的正、反演问题，也即小波框架的定义。当经伸缩和位移由基本小波引出的函数族：;(２-7）具有以下性质时：；于是称构成了一个小波框架,又称上述性质为小波框架条件,其频域可表示为：离散小波变换具有时移和非伸缩共变性。连续小波变换具有时移共变性,但对于离散小波变换，在小波框架式的情况下该性质就不存在了。2.3马尔可夫随机场的基本理论设表示有个位置的有限格点集,表示在处的随机场中，指的是定义在随机场上状态空间为的隐状态随机变量(称为图像中的标号场,为标号随机变量的可能取值)，即。设定义在上的通用的邻域系统的集合是,它满足以下特性:①②③。位置被称为的邻点，而称为的邻点集。在中有着不同的邻域结构，当子集中每对不同位置总是相邻时，则称是一个基团，而表示基团的集合。包含若干位置的集合称为基团,在退化的状况下，每个基团只含一个位置，处理图像时,此情况被认为是像素间无相互作用;另外一种不好的状况是的子集都是基团，可表示为，处理图像时，该情况被认为是所有像素互相影响。图像的局部性质很大程度取决于基团的选取,因为基团是构成马尔可夫先验模型的关键环节。处理图像时，通常采用的邻域系统是各方向同性的。把当作晶格，此时邻域系统可定义为：，式中的表示距离函数,经常使用市区距离、欧氏距离等函数,表示邻域系统的阶次。在时,其关系特性为:，显然边界位置的邻点比区域的内部少。设是所有可能组态的集合,如果随机场x是关于通用邻域系统的马尔可夫随机场，则满足以下两个条件:①；②。式中的和分别表示概率和条件概率。邻域系统中的马尔可夫随机场的含义为:在任意晶格点的剩余格点位置上的随机变量取值确定的条件下,随机场在格点处的取值概率只取决于格点的相邻的点。条件②中的条件概率常被称为马尔可夫随机场的局部特性,任何满足条件①的概率的过程都由条件②中的条件所唯一确定，但在实际应用中这两个条件概率很难确定。在8０年代Hammersｌey-Cliｆｆord发现了Gibｂｓ分布与马尔可夫随机场的关系，马尔可夫随机场概率分布这一难题可通过Giｂbｓ分布来解决[14]金丽亚,徐凤东,马龙生.灰度和彩色图象分割的一种改进算法[J].系统仿真学报.2001,13(增):111-116.4]。[14]金丽亚,徐凤东,马龙生.灰度和彩色图象分割的一种改进算法[J].系统仿真学报.2001,13(增):111-116.邻域系统的Ｇiｂbs分布是由在上的概率测度来定义的,表达形式为：（2－８)式中的是配分函数或归一化常数，满足,其中T是温度常数,为能量函数，表示邻域系统中所包含基团的集合，定义在基团上的势函数是，它只确定于，的值。Haｍmｅrｓlｅy－Clifford定理引出了Gibbｓ分布与马尔可夫随机场的等价条件:关于邻域系统的马尔可夫随机场的一个随机场当且仅当这个随机场是关于邻域系统的Gｉｂbs分布。邻域系统的马尔可夫随机场与Gibbs分布等价形式表示为:(2-９)式(2-9）解决了求解马尔可夫随机场中概率分布的问题，让对马尔可夫随机场模型的研究转变成对势函数的研究,使能量函数与Giｂbs分布建立了等价关系,是研究邻域系统马尔可夫随机场的一个重要标志[15]巩萍,潘冬明.小波分析及其在图像处理中的应用[J].长沙大学学报,2005,19(2):52-54.[15]巩萍,潘冬明.小波分析及其在图像处理中的应用[J].长沙大学学报,2005,19(2):52-54.3图像分割中的小波阈值法3．1小波阈值法的原理阈指的是一个系统的界限，而阈值是指阈的数值,图像中的“阈值”命令就是把彩色或灰度图像转换成高对比度的黑白图像，“阈值”命令有利于确定图像的最亮或最暗区域,比如可把阈值设为某个色阶，而所有比该色阶暗的像素变为黑色；所有比该色阶亮的像素变为白色。阈值分割法就是一种基于区域的图像分割技术。其原理是经过不同特征阈值的设定，将图像像素点分为很多类。常用的特征包含直接来自于原始图像的彩色或灰度特征，由彩色或原始灰度值变换而得到的特征[16]赵雪松,陈淑珍.综合全局二值比与边缘检测的图像分割方法[J].计算机辅助设计与图形学学报.2001,13(2):118-121.。设为原始图像,根据一定的准则在中找到多个特征值,其中[16]赵雪松,陈淑珍.综合全局二值比与边缘检测的图像分割方法[J].计算机辅助设计与图形学学报.2001,13(2):118-121.(3-1）一般阈值运算可看作是对图像中某点的某种局部特性、某点的灰度以及该点在图像中的位置的函数,此阈值函数可记为：，式中是点的灰度值;而是点的局部领域特性。若只选取一个阈值则称为单阈值分割，图像被它分为目标和背景两大类；若用多个阈值进行分割则称为多阈值方法，被分割的图像被分为多个目标区域与背景，为了区分目标，还要标记各个区域。阈值分割法是对灰度图像的一种假设：目标、背景内的相邻像素的灰度值是相似的,然而不同目标或背景的像素在灰度上也不同,在图像直方图上反映出不同目标与背景对应不同的峰。选取的阈值应位于两个峰之间的谷,把各个峰分开。阈值分割的优点是实现简单,当不同的物体灰度值或其他特征值有很大差异时,它能有效的分割图像。阈值分割通常作为图像的预处理，然后应用其他分割方法进行后处理。但是传统的阈值分割方法由于噪声的存在很难合理选取阈值。基于小波变换的阈值分割法的基本原理是:先由二进制小波变换把图像的直方图分解为各层次的小波系数，然后按照一定的分割准则和小波系数选取阈值门限,最终利用阈值标出图像分割的区域。整个分割过程是由粗到细，通过尺度变化来控制，起始分割由子空间上投影的直方图来实现。若分割不理想，则用直方图在子空间的小波系数逐步细化图像分割。分割算法的计算量和图像尺寸的大小是线性关系。基于小波变换的阈值法图像分割技术能够有效弥补传统的图像阈值法分割方法的不足，具有较强的抗噪声性能，对于直方图为多峰值时,可利用小波的多分辨率特性，合理选择灰度阈值，对图像进行分割处理[1６]。3.2图像直方图的多分辨率分析设每个整数(为整数集合）,是指在分辨率下的二进制有理数。对于任意,是一组在实数轴上等间隔采样点的集合，若,则表示低分辨率的采样点；但时,则表示高分辨率的采样点。假设为一幅图像,是该图像中最大灰度,那么直方图表示为:(3-2)式中“|…｜”指计数操作，是离散函数。令,离散函数表示为连续函数,而可看作是由若干分段常数函数组成。对于，根据采样点采样，则表示在分辨率下的直方图。同时可用Hａar尺度函数的伸缩与平移表示，即(３-３）(3-４）由于连续函数由几个分段常数函数组成,存在阶梯现象,则要用滤波操作处理函数,消除其高频成分。多层表达曲线可以表示为:（3－5）（3－6）(3-７)则可得出相应的阈值分割算法［１7]。3.３阈值分割算法按上式图像的直方图曲线可进行多分辨分解,通过分解系数,则可以找出图像的分割阈值。如下是算法步骤:步骤1：设分割区域为，分解级数,是图像中的最大灰度值；步骤2：在对直方图曲线进行级二进制小波分解时,则可以得到；步骤3:在分解系数中，找到满足与条件的标号（该标号值就是所求的灰度阈值),更新灰度阈值数组，然后统计标号的个数;步骤4:若,则,当时,转向步骤2;若,则合并靠得近的灰度阈值,更新灰度阈值数组，使,另外，当时,转向步骤２;步骤5：比较图像中的每个像素值与灰度阈值，分割图像。３.4实验仿真实验采用MATLAB程序语言来进行编程仿真。实验仿真结果如下图,其中(a)为原始图像;(b)为该图像的灰度直方图,由直方图可看出,峰值较多,很难直接选取门限;(c)为阈值为90时用传统的阈值法进行分割后的图像；(d)为通过小波阈值门限分割后的图像，该阈值同为９0。由传统的阈值分割法与小波阈值分割法在相同的阈值下分割的图像进行比较。其中小波阈值分割方法的程序见附录1。(a)原始图像(b)直方图(c)阈值为90时的传统阈值分割图像(d)阈值为90时的小波阈值分割图像图3-1不同方法对图像进行分割从仿真结果可知,基于小波变换的阈值法图像分割技术能有效地弥补传统的图像阈值法分割技术的不足之处，有较强的抗噪声性能,同时在直方图为多峰值的情况,可利用小波的多分辨率分解,合理选择灰度阈值，实现对图像的分割处理。４图像分割中的小波域马尔可夫随机场方法4.1图像分割中的马尔可夫随机场方法综述图像分割本质上是不适定的,如何引入先验信息是解决这一问题的重要工具。基于图像模型的分割易于实现,还可以提高分割的质量，因此对图像模型的研究既是各种应用的需要,又是图像分割能否取得突破性进展的重要问题。如今,引人注目的是利用统计图像模型，其中马尔可夫模型比较有影响力。对于图像灰度的变化率、纹理基元的排列、边缘几何的特征等信息,描述待分割图像的先验分部用马尔可夫模型很适合。马尔可夫随机场方法是建立在MＲF模型和Ｂａyeｓ理论的基础上形成的,马尔可夫模型提供了不确定性描述与先验知识的联系,根据估计理论和统计决策中的最优准则利用观测图像来确定图像的目标函数,以此将分割图像转化为最优化问题。马尔可夫模型在分割图像中具有以下鲜明的特点:eq\o＼ａc（○,１)ＭRF模型不仅可以反映图像的随机性，而且能反映图像的潜在结构，从而有效地描述图像的性质。ｅq＼ｏ\aｃ(○,２)ＭＲF模型的局部特征能采用大规模并行算法。eｑ＼ｏ\aｃ(○,３)MＲF模型可以把像素空间关系紧密结合，从而可以用低阶ＭRF描述像素间的作用关系。eｑ\o\ａｃ(○,4)ＭRＦ模型从物理模型出发,同时与图像数据(即特征与灰度值)有联系。eｑ\ｏ＼ac(○,５）用MＲＦ描述来求解具有不确定性，利用Bayes理论将图像的先验知识转变为先验分布模型,采用最大先验估计得到图像标号的分布，是很好的数学推导,得到的参数有明确的物理含义，但是线性模型、时间序列模型中的参数只被作为被拟合的参数，一般无真实含义。解是否正确可以通过MｏnteCarlo的方法生成随机数来进行相关验证。4.１.1用马尔可夫随机场描述图像模型图像分割是根据像素、区域属性给每个像素分配标号的过程。通常用两个随机场在马尔可夫随机场中来描述待分割的图像，一个是像素的标号场，被称为隐随机场,用先验分布来描述标号场的局部相关性;另一个是特征场或灰度场，常根据标号场用分布函数来描述观测数据或特征向量的分布。小波域马尔可夫先验模型根据小波变换的特点,主要分为尺度间马尔可夫先验模型、尺度内马尔可夫先验模型和混合先验模型三种形式。小波域尺度间马尔可夫先验模型有利于形成较大一致性区域，在研究文本图像和遥感图像的分割时,HｙeokhoChｏｉ等[17]程杰.一种基于直方图的分割方法[J].华中理工大学学报.1999,27(1):20-23.7]把小波多尺度思想引入图像分割，为了精确反映出小波系数尺度间的依赖关系,采用了父节点周围36个节点的标号对子节点标号确定影响,让马尔可夫参数具有因果性和可变性,从而形成有效的计算方法,并取得较好的进行分割。但文中没有给出标号场的具体表达式，认为父节点的标号是相互独立,确定各自对子节点的标号有影响，这与马尔可夫的局部特性不相符。SuｎＪ.等［１9]结合空域等思想改进了小波域尺度间标号场先验模型,把父节点与父节点的兄弟节点的标号对子节点标号的所有作用写成了统一的表达式，这与马尔可夫的局部特性相符。但是此尺度间小波域马尔可夫先验模型对分割的边界定位不精确且对噪声敏感。小波域尺度内的马尔可夫先验模型可采用空域非因果的马尔可夫先验模型的表达形式，这种模型对图像的边缘很敏感,并易于形成局部一致性区域。ＦaｎGuｏliang等研究了混合马尔可夫先验模型对纹理图像分割所产生的影响，研究结果表明混合先验模型对边界定位精确，而容易将区域误分[17]程杰.一种基于直方图的分割方法[J].华中理工大学学报.1999,27(1):20-23.为了准确刻画图像的特征属性,2001年Hyｅokｈｏｃhｏｉ等NOＴEＲEＦ_Rｅf2937６0258\ｆ\h\＊MＥRGEFORＭＡT[17］提出了HMT模型,把图像分解的每个小波系数分为两个状态，每个较大的状态对应图像的边缘,每个较小的状态对应图像结构特征的平滑区域，并认为每个状态都服从高斯分布，即每个小波系数用两个独立的高斯混合模型来描述，每个子带可用有限高斯混合模型来逼近。假设图像小波分解中的各子带间是相互独立的,应用于纹理图像的分割取得了较好的整体属性。但如果引入过多的隐状态到每个小波系数，则会增加计算的复杂度，同时很难解释隐状态的物理含义。侯玉华等通过秩统计量分析小波分解中三个子带间的相关性,了解到三个子带间存在相关性，并其相关性随着尺度的增加而渐渐增强,便提出改进的HMT模型的文本图像分割方法，并取得了比Hyｅokhochoi等认为的小波分解图像子带间相互独立的分割结果更好,然而增加了计算的复杂度。20０1年，每个小波系数引入三个隐状态,并认为图像的小波系数服从一个方差较大的高斯分布、一个方差较小的高斯分布和一个方差适中的高斯分布混合而成，应用文本图像分割取得了更好的分割效果。２004年在研究纹理图像分割时,sunJ．等用复小波分解图像,用有限Rayｌeｉgh混合分布来逼近小波系数的分布。小波域特征场模型建立的主要困难就是每个小波系数引入几个隐状态和如何合理的解释引入的隐状态物理概念。4.1.2基于适当最优准则实现图像的分割给定标号场的先验分布和灰度场的条件分布后,按照Bayｅｓ理论,在给定观测图像的条件下,的表达式为：(4－1）下面给出几种经常使用的图像分割标准：(1)按照MAP(maｘｉｍｕmaposterioｒ)准则来分割图像[18][18]赵立初,王积分.基于小波分析的图像自适应阈值选择[J].模式识别与人工智能.1999,12(1):79-84.在已知图像标号的先验分布和在给定标号的前提下分布图像灰度场(或特征场）的条件，那么ＭAＰ估计就是最大后验概率的解,即(４-2)实质上ＭＡP估计是缩小图像中的像素分类错误的概率，因为任何分割算法都可能会分类错误。MＡP估计在实际应用中有非理想的全局性,如在多分辨率分析下，随着尺度的增加而ＭRF的系数也会增加，这违背了人的直觉,因为越粗糙的尺度MRF的局部相应的系数应该越小。此外MAＰ估计会使图像的边缘模糊化,并消除细小的结构,这主要是因为标号场的先验分布把区域边缘的各向异性与区域内的一致性合并起来作为图像分割的正则条件不当所引起，所以导致图像分割不适定的问题没有得到最好的答案。这种准则一般应用于对分割对象精度要求不太高,但是对物体的纹理要求较高的场合。(２)基于ＭPM估计准则[29]靳宏磊[29]靳宏磊,朱蔚萍,李立源等.二维灰度直方图的最佳分割方法[J].模式识别与人工智能.1999,12(3):329-333.（4-3)作为MAP估计的一种替代，最大后验边缘估计是最小分类像素的期望值。已表明ＭＰM估计标准比MAP估计标准更加适合图像分割，因为不管错误的分割像素的数目与真正的分割多么不同，MＡP估计分配都会给每个不正确分割一样的代价；但是MＰＭ估计是根据不正确分割的像素数分配给相应的代价。而MPM估计需要模拟退火类似的复杂计算,当进行分割的质量区分时,MＰＭ不考虑空间位置的误差。（３)基于最小方差估计准则来分割图像(4－4)式（4-4)中是一种基于观测图像的最小方差无偏估计,从空间投影的角度来确定,此估计是使与的空间距离最小。这种方法是在总体平均条件下的分割误差最小,其局部有可能会出现严重的分割误差。（４）基于SMAP估计准则来分割图像（4－5)随着观察图像的尺度增加,SMAP估计标准是逐渐分配给较大的误差代价，由于较大的尺度容易导致更多的像素误分类,具有从粗尺度到细尺度一系列优化分割的优点,这符合人类的视觉原理。从估计的精确度来说，它比ＭAP标准优秀;从算法的可行性来看,它形成了一套更加有效的计算方法,需要更小的计算量。4.2小波域马尔可夫随机场模型的分割算法为了图像的分割，在给定观测数据的条件概率、观测数据和标号场先验概率的前提下，精确估计像素的标号，所以采用贝叶斯估计器把分割问题转变为最优化问题:（４-6）式中为实际分割和期望分割不相同时的耗费函数，通过最优化这一耗费函数的数学期望，从而贝叶斯估计器求出最佳标号[20]薛景浩,章毓晋,林行刚[20]薛景浩,章毓晋,林行刚.一种新的图像模糊散度阈值比分割算法[J].清华大学学报(自然科学版).1999,31(1):12-16.4．2.１小波域马尔可夫随机场模型的MAP准则的分割算法假设（４－7）式中:如果，则;此外。把式(4－7）代入式（4-6)，则基于ＭAP准则的估计器为：（４-8)由贝叶斯准则,假设任意节点都形成独立小波系数树(叶子节点除外),有(4-9）其中以节点为根的树不包含节点的全部节点表示为，表示在位置集中不包含位置。所以（４-10）需要的最优估计，可以使上式的两项达到最大值。顺序类推,可知从粗尺度到细尺度来确定小波系数的标号完成最终图像的分割过程：（4-1１)4.2.２小波域马尔可夫随机场模型多尺度概率值的计算由于任何节点小波系数的概率值只由所在分辨率的小波系数和该小波系数的标号所决定，所以对于以节点为根的小波系数树的概率值为:（４-12）,是指以为根节点的小波系数树,是指以节点为跟的小波系数树中不包括根节点的全部节点,而可从细尺度到粗尺度用递归运算求出[21]HyeokhoChoi,RichardGBaraniuk.Multiscaleimagesegmentationusingwavelet-domainhiddenmarkovmodels[J].IEEETransactionsonimageprocessing,2001,10(9):1309-1321.1]。[21]HyeokhoChoi,RichardGBaraniuk.Multiscaleimagesegmentationusingwavelet-domainhiddenmarkovmodels[J].IEEETransactionsonimageprocessing,2001,10(9):1309-1321.4．2.3小波域马尔可夫随机场模型分割算法的具体实现(１)由细尺度到粗尺度来计算各节点小波系数的概率值a．初始化使（4-13)ｂ.当Ｌ＝L+1时,移到树的更粗尺度（4－１4）ｃ.当L=l时,立即停止，否则转向b。(2)由粗尺度到细尺度来完成MAP分割a.初始化使，（4－15)b.通过递归运算Ｌ=Ｌ-１,移到树的更精细尺度（4-16)(３)如果L=１,那么停止，否则转到ｂ。４.3小波域马尔可夫随机场模型的参数估计由标场号和特征场组成的马尔可夫层次模型的参数可记为：（4－17)要分割图像，必须计算出马尔可夫层次模型的参数，因为只已知观测图像,被称为“不完全数据”;同时引入标号过程，被称为“丢失数据”,与观测数据组成“完全数据”。但是EＭ算法[22]SunJ,GuD,ZhangS,etal.HiddenMarkovBayesiantexturesegmentationusingcomplexwavelettransform[J].IEEProceedingsvision,imageandsignalprocessing,2004,151(3):215-223.[22]SunJ,GuD,ZhangS,etal.HiddenMarkovBayesiantexturesegmentationusingcomplexwavelettransform[J].IEEProceedingsvision,imageandsignalprocessing,2004,151(3):215-223.附录1.小波阈值程序：x=imread(‘lena.jpg');%读取原图像figure(1);imshow(x);title(‘原始图像’);%显示原图像s=rgb2gray(x);imhist(s);title('直方图');[Row,COL]=size(x);T=90;%设置阈值s=double(x);fori=1:Rowforj=1:COLif(x(i,j)>T)x(i,j)=255;elsex(i,j)=0;endendendfigure(2);imshow(x);%显示经小波阈值分割的图像title('小波阈值分割图像');2.本文方法分割程序：%这个是2D-DWT的函数，是haar小波%c是图像像素矩阵steps是变换的阶数functiondwtc=dwt_haar(c,steps)%DWTC=CWT_HARR(C)-DiscreteWaveletTransformusingHaarfilter%MDPlumbleyNov2003N=length(c)-1;%Maxindexforfilter:N%Ifnostepstodo,orthesequenceisasinglesample,theDWisitselfif(0==N|steps==0)dwtc=c;returnend%CheckthatN+1isdivisibleby2if(mod(N+1,2)~=0)disp(['Notdivisible2:'num2str(N+1)]);returnend%SettheHaaranalysisfilterh0=[1/21/2];%HaarLow-passfilterh1=[-1/21/2];%HaarHigh-passfilter%Filterthesignallowpass_c=conv(h0,c);hipass_c=conv(h1,c);%Subsamplebyfactorof2andscalec1=sqrt(2)*lowpass_c(2:2:end);d1=sqrt(2)