基于Contourlet域多状态HMT模型的含噪图像增强方法.pdf

(19)中华人民共和国国家知识产权局*CN102968771A*(12)发明专利申请(10)申请公布号CN102968771A(43)申请公布日2013.03.13(21)申请号201210508502.6(22)申请日2012.12.03(71)申请人北方民族大学地址750021宁夏回族自治区银川市西夏区文昌北路204号(72)发明人常霞高岳林黄永东马自萍(74)专利代理机构宁夏专利服务中心64100代理人赵明辉(51)Int.Cl.G06T5/00(2006.01)权利要求书权利要求书2页2页说明书说明书99页页附图附图66页(54)发明名称基于Contourlet域多状态HMT模型的含噪图像增强方法(57)摘要本发明涉及在Contourlet域上，利用多状态隐马尔科夫树HMT模型，对含噪图像进行增强的方法，尤其是基于Contourlet域多状态HMT模型的含噪图像增强方法。本发明方法具有以下优点：1.采用了非下采样Contourlet变换，可以有效捕捉到含噪图像中的线状和轮廓状几何信息，有利于图像的后续处理。2.引入了循环平移策略，可以有效避免因Contoutlet变换缺乏平移不变性，而易在增强后的图像中产生的人工失真。3.通过构建Contourlet域HMT模型的多种隐状态，进一步模拟出Contourlet域系数的近噪声和近强边缘成分，有利于获得细节特征突出，均匀区域平滑的增强结果。CN1029687ACN102968771A权利要求书1/2页1.一种基于Contourlet域多状态HMT模型的含噪图像增强方法，其特征在于，包括如下步骤：(1)输入大小为M×N的含噪图像im0(m，n)，其中1≤m≤M，1≤n≤N，M和N均为大于1的自然数，初始化计数量tem=0，行位移量u=1，列位移量v=1；(2)计数量递增，即tem=tem+1，对含噪图像im0(m，n)进行(u，v)步循环平移，得到循环平移后的含噪图像imtem(m，n)；(3)对循环平移后的含噪图像imtem(m，n)进行L层Contourlet变换，得到低频子带系-l数YA和高频方向子带系数其中0≤l≤L-1，1≤k≤kl，kl表示在尺度2上的高频方向子带的数目，i∈{1，2，3，...，Q}为Contourlet域系数四叉树的索引标记，Q为Contourlet域系数四叉树上待观察系数的数目，L为3～5，L为自然数；(4)采用期望最大化算法对步骤(3)得到的高频方向子带系数进行训练，得到确定的Contourlet域多状态HMT模型参数向量Θ，式中P0，k为Contourlet域系数四叉树根节点的隐状态概率密度向量，为高频方向子带系数的隐状态概率，S为隐状态变量，隐状态的类别数取为5，分别对应Contourlet域系数中的噪声、近噪声、弱边缘、近强边缘和强边缘成分，有{S|S=r，r∈{0，1，2，3，4}}，εl，k为尺度l上第k个高频方向子带系数的状态转移概率矩阵，σl，k和μl，k分别为尺度l上第k个高频方向子带系数的标准差和均值向量；(5)利用步骤(4)得到Contourlet域多状态HMT模型参数向量Θ，对高频方向子带系数进行增强；(6)对步骤(3)得到的低频子带系数YA和步骤(5)得到增强后的高频方向子带系数进行Contourlet逆变换，得到含噪图像imtem(m，n)的增强图像(7)对步骤(6)得到的含噪图像imtem(m，n)的增强图像进行(u，v)步逆向循环平移，得到图像并对其进行存储；(8)列位移量递增，即v=v+1，如果v<=8，转步骤(2)，否则转步骤(9)；(9)行位移量递增，即u=u+1，如果u<=8，转步骤(2)，否则转步骤(10)；(10)计算并得到输出含噪图像im0(m，n)的增强图像即2.如权利要求1所述的一种基于Contourlet域多状态HMT模型的含噪图像增强方法，其特征在于：其中步骤(5)所述的利用Contourlet域多状态HMT模型参数向量Θ，对高频方向子带系数进行增强，具体按如下过程进行：1)利用高频方向子带系数的隐状态概率将区分为五种类别，有2CN102968771A权利要求书2/2页2)利用增益规则公式，对获得了类别属性的高频方向子带系数进行初始修正，即：式中，Ai为所在高频方向子带系数的最大模值，设置a0=0，a1=1，a3=2，a2(·)为非线性增益函数，有：式中，参数t用于控制含噪图像增强的程度，设置参数t=6，e为自然常数，e=2.71828；3)利用高频方向子带系数的隐状态概率和步骤2)中得到的高频方向子带系数的初始修正结果计算增强后高频方向子带系数即式中，是尺度l上第k个高频方向子带的噪声系数的标准差，并且通过蒙特卡罗方法获得。3CN102968771A说明书1/9页基于Contourlet域