发明名称 基于空间引导滤波的图像细节增强方法
摘要 本发明公开了一种基于空间引导滤波的图像细节增强方法,其特征是按如下步骤进行:首先,利用图像边缘检测算子对源图像提取边缘特征响应图,并进行归一化;然后,对于不同的灰度区间分别建立二值化的空间指示图,并对每个空间指示图进行高斯卷积,得到空间滤波图,并计算每个空间滤波图的权值;其次,计算累加图,并对累加图进行引导图像滤波得到空间引导图;最后,求解基础图像和残差图像,建立基于空间引导滤波的图像细节增强模型对源图像进行图像细节增强。本发明能有效提升对图像细节的增强效果。
申请公布号 CN104318525A 申请公布日期 2015.01.28
申请号 CN201410552788.7 申请日期 2014.10.17
申请人 合肥工业大学 发明人 汪萌;郝世杰;杨勋;韦希明
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 安徽省合肥新安专利代理有限责任公司 34101 代理人 何梅生
主权项 一种基于空间引导滤波的图像细节增强方法,其特征是按如下步骤进行:步骤1,对于分辨率为m×n的源图像I,I∈R<sup>m×n</sup>,利用图像边缘检测算子提取所述源图像I的边缘响应图,并将所述边缘响应图除以255进行归一化,获得归一化边缘响应矩阵|▽I|,|▽I|∈R<sup>m×n</sup>;m为源图像I的长度;n为源图像I的宽度;步骤2,将所述归一化边缘响应矩阵|▽I|的灰度范围区间[0,1]均匀划分成k个子区间Ω<sub>i</sub>,利用式(1)分别对所述k个子区间Ω<sub>i</sub>中的每个子区间建立相对应的空间指示图M(i),M(i)∈R<sup>m×n</sup>,1≤i≤k;<img file="FDA0000588493580000011.GIF" wi="1043" he="172" />式(1)中:(x,y)表示所述归一化边缘响应矩阵|▽I|中元素的坐标,并对应于所述源图像I中像素的位置;1≤x≤m,1≤y≤n;|▽I(x,y)|∈Ω<sub>i</sub>表示所述归一化边缘响应矩阵|▽I|中第x行第y列元素的取值位于第i个子区间Ω<sub>i</sub>中;M<sub>(x,y)</sub>(i)表示第i个空间指示图M(i)中第x行第y列元素的取值;步骤3、利用高斯卷积依次对每一个空间指示图进行滤波处理,得到k个空间滤波图I<sub>Gauss</sub>(i),I<sub>Gauss</sub>(i)∈R<sup>m×n</sup>;步骤4、依次统计所述归一化边缘响应矩阵|▽I|中所有元素落在每个子区间中元素的个数,并利用式(2)计算每个空间滤波图的权重;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><msub><mi>h</mi><mi>i</mi></msub><mi>mn</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000588493580000012.GIF" wi="666" he="142" /></maths>式(2)中,h<sub>i</sub>表示所述归一化边缘响应矩阵|▽I|中的元素落在第i个子区间Ω<sub>i</sub>中元素的个数;w(i)表示落在第i个空间滤波图的权重;步骤5、利用式(3)计算累加图S<sub>a</sub>,S<sub>a</sub>∈R<sup>m×n</sup>;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>S</mi><mi>a</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>I</mi><mi>Gauss</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000588493580000013.GIF" wi="676" he="155" /></maths>步骤6、将所述源图像I作为引导图,采用引导图像滤波方法对所述累加图S<sub>a</sub>进行引导图像滤波处理,得到空间引导图S,S∈R<sup>m×n</sup>;步骤7、将所述源图像I作为引导图,采用引导图像滤波方法对所述源图像I自身进行引导图像滤波处理,得到基础图像I<sub>b</sub>,I<sub>b</sub>∈R<sup>m×n</sup>;步骤8、利用式(4)所示的图像细节增强模型,对所述源图像I进行图像细节增强处理,获得细节增强图像I<sub>o</sub>,I<sub>o</sub>∈R<sup>m×n</sup>;I<sub>o</sub>=I<sub>b</sub>+S<sub>0</sub>·S⊙I<sub>r</sub>   (4)式(4)中:S<sub>0</sub>为滤波强度;S<sub>0</sub>为标量;⊙为哈达玛乘积符号,表示两个矩阵对应相乘;I<sub>r</sub>表示残差图像,并有I<sub>r</sub>=I‑I<sub>b</sub>。
地址 230009 安徽省合肥市包河区屯溪路193号