发明名称 贝叶斯非局部均值图像修复方法
摘要 本发明公开了一种贝叶斯非局部图像修复方法,主要解决现有基于样本的非局部均值修复方法中相似块搜索不准确及参数值由经验确定的问题。其实现过程是:(1)对待修复图像I,确定待修复区域Ω及其边界δ;(2)找出边界上优先级最高的修复块<img file="dda00001915900500011.GIF" wi="71" he="43" />并对<img file="dda00001915900500012.GIF" wi="53" he="43" />使用贝叶斯框架建模;(3)通过自适应阈值的对搜寻区域进行预选取。(4)搜索与修复块<img file="dda00001915900500013.GIF" wi="53" he="44" />最相似的m个样本块,将其加权均值作为修复块的填充块<img file="dda00001915900500014.GIF" wi="83" he="43" />(5)更新置信度项及待修复区域,重复步骤(1)~(5)直至待修复区域内的所有点被修复。本发明可用于修复图像破损区域,修复图像划痕和图像中文本的去除。
申请公布号 CN102800077A 申请公布日期 2012.11.28
申请号 CN201210253269.1 申请日期 2012.07.20
申请人 西安电子科技大学 发明人 钟桦;焦李成;朱波;王桂婷;侯彪;王爽;张小华;田小林
分类号 G06T5/50(2006.01)I 主分类号 G06T5/50(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种贝叶斯非局部均值图像修复方法,包括如下步骤:(1)对于输入的待修复图像I,确定待修复区域Ω和待修复区域的边界δ;(2)利用如下公式,计算出中心点在待修复区域的边界δ上的所有修复块的优先级P(p):P(p)=C(p)·D(p),其中,D(p)为数据项,C(p)为置信度项,表示图像像素点的可信程度,C(p)初始化为C(p)=0,p∈Ω,C(p)=1,p∈I-Ω;(3)以优先级最高的修复块<img file="FDA00001915900200011.GIF" wi="53" he="44" />的中心点<img file="FDA00001915900200012.GIF" wi="30" he="52" />为中心,选取大小为M×M的邻域作为该修复块的搜索区域,定义该区域内以点<img file="FDA00001915900200013.GIF" wi="24" he="52" />为中心的块<img file="FDA00001915900200014.GIF" wi="49" he="44" />为样本块;(3.1)对修复块<img file="FDA00001915900200015.GIF" wi="53" he="44" />和样本块<img file="FDA00001915900200016.GIF" wi="49" he="44" />利用贝叶斯框架建模,分别计算修复块的均值<img file="FDA00001915900200017.GIF" wi="56" he="46" />和样本块的均值<img file="FDA00001915900200018.GIF" wi="63" he="45" />并计算出它们的均值差<img file="FDA00001915900200019.GIF" wi="93" he="45" />即<img file="FDA000019159002000110.GIF" wi="284" he="45" />(3.2)根据均值差<img file="FDA000019159002000111.GIF" wi="68" he="45" />服从<img file="FDA000019159002000112.GIF" wi="191" he="56" />的高斯分布的特性,定义一个自适应阈值:t=λσ<sub>0</sub>,其中,σ<sub>0</sub>为<img file="FDA000019159002000113.GIF" wi="68" he="45" />的标准差,λ=1.65,u<sub>0</sub>为<img file="FDA000019159002000114.GIF" wi="68" he="45" />的均值,<img file="FDA000019159002000115.GIF" wi="49" he="56" />为<img file="FDA000019159002000116.GIF" wi="68" he="45" />的方差;(3.3)选择搜索区域内所有<img file="FDA000019159002000117.GIF" wi="249" he="57" />的点,作为预选取后新的搜索区域;(4)计算新的搜索区域内的修复块<img file="FDA000019159002000118.GIF" wi="53" he="44" />与样本块<img file="FDA000019159002000119.GIF" wi="49" he="44" />的相似性距离:<img file="FDA000019159002000120.GIF" wi="484" he="63" />其中,<img file="FDA000019159002000121.GIF" wi="89" he="56" />为2范数;(5)根据相似性距离<img file="FDA000019159002000122.GIF" wi="195" he="57" />服从自由度为n的卡方分布X<sup>2</sup>(n)的特性,当n≥25时,分位点<img file="FDA000019159002000123.GIF" wi="441" he="106" />选取相似性距离位于分位点β左侧的最相似的m个样本块,作为修复块<img file="FDA000019159002000124.GIF" wi="53" he="44" />的集合;(6)根据如下公式,计算集合内的样本块<img file="FDA000019159002000125.GIF" wi="49" he="44" />和修复块<img file="FDA000019159002000126.GIF" wi="53" he="44" />的相似性权值:<img file="FDA000019159002000127.GIF" wi="499" he="117" />其中,Z为归一化参数,σ<sup>2</sup>为修复块<img file="FDA000019159002000128.GIF" wi="53" he="43" />的方差,N为修复块<img file="FDA000019159002000129.GIF" wi="53" he="43" />中像素值已知的 点的个数;(7)根据相似性权值,将集合内的全部样本块的加权均值,作为填充块Ψ<sub>0</sub>,并用该填充块对修复块<img file="FDA00001915900200021.GIF" wi="53" he="44" />进行填充修复;(8)当修复块<img file="FDA00001915900200022.GIF" wi="53" he="44" />完成修复后,更新待修复区域,并用修复块<img file="FDA00001915900200023.GIF" wi="53" he="44" />中心点<img file="FDA00001915900200024.GIF" wi="30" he="51" />的置信度更新已完成修复的点的置信度C(p):<img file="FDA00001915900200025.GIF" wi="277" he="62" /><img file="FDA00001915900200026.GIF" wi="249" he="56" />其中,<img file="FDA00001915900200027.GIF" wi="111" he="62" />为修复块<img file="FDA00001915900200028.GIF" wi="53" he="44" />中心点<img file="FDA00001915900200029.GIF" wi="30" he="52" />的置信度,∩表示‘与’关系;(9)重复步骤(1)~(8),直至待修复区域内的所有点被修复。
地址 710071 陕西省西安市太白南路2号