发明名称 一种基于格林函数的图像局部重建方法
摘要 本发明涉及一种基于格林函数的图像局部重建方法。本发明首先采用稳定场来描述图像局部纹理,建立并求解图像局部纹理的稳定场方程,提出基于格林函数法的图像局部重建模型;利用格林函数的定义及物理意义,实现上述重建模型;对上述计算中涉及的梯度计算,提出“单向计算梯度分量”模式,并利用“就近原则”对参与梯度分量单向计算的相邻点设置优先级。本发明对每一个缺损点都有一个计算与赋值的过程,且梯度计算模式的设置,使得缺损边缘特征能够及时全面地被反映在梯度及格林函数计算中,这些使得本发明在图像准确方面取得了良好结果。且由于本发明的方法没有使用迭代计算,重建效率也比较高。
申请公布号 CN104732502A 申请公布日期 2015.06.24
申请号 CN201510075202.7 申请日期 2015.02.12
申请人 杭州电子科技大学 发明人 叶学义;齐珍珍;惠舒芸;陈华华;赵治栋;何志伟
分类号 G06T5/20(2006.01)I 主分类号 G06T5/20(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 杜军
主权项 一种基于格林函数的图像局部重建方法,其特征在于:首先,对图像局部纹理的缺损区域建立稳定场模型及方程:‑ΔI(r)=f,(r∈U);其中,I(r)为描述图像局部纹理的函数;若为灰度图像,I(r)的值表示对应空间点的灰度值;若为彩色图像,I(r)表示其中任一颜色通道的值;r表示图像局部纹理区域的空间坐标,U表示整个图像局部纹理区域;Δ为拉普拉斯算子;对于稳定场,f表示场源,在图像重建中表示缺损区域周围的已知区域像素信息;林函数法求解上述方程得到图像局部重建模型:<img file="FDA0000672103190000011.GIF" wi="472" he="126" />其中,为格林函数,<img file="FDA0000672103190000012.GIF" wi="49" he="65" />为已知区域的像素值,i表示已知区域图像点序号;然后,采用如下计算方法进行格林函数的计算:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>G</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><msub><mi>r</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>g</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><msub><mi>r</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>k</mi><mrow><mo>(</mo><msub><mi>r</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>20</mn></munderover><mi>g</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><msub><mi>r</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>k</mi><mrow><mo>(</mo><msub><mi>r</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000672103190000013.GIF" wi="669" he="269" /></maths>其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>g</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><msub><mi>r</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msqrt><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><mo>|</mo><mo>|</mo><msub><mrow><mo>&dtri;</mo><mi>I</mi></mrow><mn>0</mn></msub><mrow><mo>(</mo><msub><mi>r</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>|</mo><mo>&times;</mo><msub><mi>d</mi><msub><mi>rr</mi><mi>i</mi></msub></msub><mo>&times;</mo><mi>cos</mi><mi>&theta;</mi><mo>|</mo></mrow><mrow><msub><mi>I</mi><mn>0</mn></msub><mrow><mo>(</mo><msub><mi>r</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>,</mo></mrow>]]></math><img file="FDA0000672103190000014.GIF" wi="962" he="253" /></maths><img file="FDA0000672103190000015.GIF" wi="704" he="195" />g为已知区域像素点r<sub>i</sub>对缺损点r的影响;<img file="FDA0000672103190000016.GIF" wi="169" he="92" />为<img file="FDA0000672103190000017.GIF" wi="150" he="73" />的模,<img file="FDA0000672103190000018.GIF" wi="160" he="75" />为I<sub>0</sub>(r<sub>i</sub>)的梯度,I<sub>0</sub>(r<sub>i</sub>)为已知区域像素点r<sub>i</sub>的值;<img file="FDA00006721031900000110.GIF" wi="77" he="75" />为缺损点和已知点之间的距离,θ为缺损点和已知点之间连线与I<sub>0</sub>(r<sub>i</sub>)的梯度的夹角,k(r<sub>i</sub>)为权值;计算<img file="FDA0000672103190000019.GIF" wi="154" he="76" />采用单向计算梯度分量的梯度计算模式,即对缺损点周围的已知像素点梯度分量的计算,选择前向或者后向差分法来计算;对参与梯度分量单向计算的相邻点设置优先级,距离越近的相邻点的优先级越高。
地址 310018 浙江省杭州市江干区下沙高教园区2号大街