发明名称 基于方向波变换和模糊理论的图像增强方法
摘要 本发明公开了一种基于方向波变换和模糊理论的图像增强算法,主要解决传统图像增强算法只对高频信息进行增强处理,而忽略了对包含大部分图像能量和信息的低频信息处理导致图像增强效果差的问题。其实现步骤为:(1)输入原始图像,对其分别进行45度、90度及135度三个方向的方向波变换,得到对应的低频信息和高频信息;(2)对低频信息进行模糊增强处理;(3)对高频信息进行非线性抑噪增强处理;(4)将增强后的低频信息和高频信息进行方向波逆变换,分别得到45度、90度及135度增强的图像;(5)将三个增强图像进行加权平均,获得最终的增强图像。本发明具有方向选择性好,细节增强效果明显,对比度高的优点,可用于低对比度图像的增强处理。
申请公布号 CN103310414B 申请公布日期 2016.02.10
申请号 CN201310071636.0 申请日期 2013.03.06
申请人 西安电子科技大学 发明人 白静;焦李成;王爽;周华吉;高艺菡;李阳阳;马文萍;马晶晶
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种基于方向波变换和模糊理论的图像增强方法,包括如下步骤:(1)输入原始图像,对它分别进行135度、90度及45度三个方向的方向波变换,得到对应的低频信息和高频信息:(L<sub>1</sub>,H<sub>11</sub>,H<sub>21</sub>,H<sub>31</sub>,H<sub>41</sub>,H<sub>51</sub>,H<sub>61</sub>,H<sub>71</sub>);(L<sub>2</sub>,H<sub>12</sub>,H<sub>22</sub>,H<sub>32</sub>,H<sub>42</sub>,H<sub>52</sub>,H<sub>62</sub>,H<sub>72</sub>);(L<sub>3</sub>,H<sub>13</sub>,H<sub>23</sub>,H<sub>33</sub>,H<sub>43</sub>,H<sub>53</sub>,H<sub>63</sub>,H<sub>73</sub>);其中,L<sub>1</sub>、L<sub>2</sub>、L<sub>3</sub>分别表示135度、90度及45度方向波变换对应的低频信息;H<sub>11</sub>,H<sub>21</sub>,H<sub>31</sub>,H<sub>41</sub>,H<sub>51</sub>,H<sub>61</sub>,H<sub>71</sub>表示135度方向波变换对应的高频信息;H<sub>12</sub>,H<sub>22</sub>,H<sub>32</sub>,H<sub>42</sub>,H<sub>52</sub>,H<sub>62</sub>,H<sub>72</sub>表示90度方向波变换对应的高频信息;H<sub>13</sub>,H<sub>23</sub>,H<sub>33</sub>,H<sub>43</sub>,H<sub>53</sub>,H<sub>63</sub>,H<sub>73</sub>表示45度方向波变换对应的高频信息;(2)对低频信息进行模糊增强处理:(2a)通过隶属度函数将低频信息转换至模糊域:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&mu;</mi><mrow><mi>m</mi><mi>n</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>L</mi><mi>k</mi></msub><mo>-</mo><mi>n</mi></mrow><mrow><mi>m</mi><mo>-</mo><mi>n</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000795794210000011.GIF" wi="308" he="142" /></maths>其中,μ<sub>mn</sub>为求出的模糊域,L<sub>k</sub>为方向波低频系数,k=1,2,3,m=max(L<sub>k</sub>),n=min(L<sub>k</sub>);(2b)对步骤(2a)中转换所得的低频系数,通过如下函数进行增强处理:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>T</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>&mu;</mi><mrow><mi>m</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&mu;</mi><mi>c</mi></msub><mo>-</mo><msup><mrow><mo>(</mo><msup><msub><mi>&mu;</mi><mi>c</mi></msub><mn>2</mn></msup><mo>-</mo><msubsup><mi>&mu;</mi><mrow><mi>m</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>)</mo></mrow><mn>0.5</mn></msup></mrow></mtd><mtd><mrow><mn>0</mn><mo>&le;</mo><msub><mi>&mu;</mi><mrow><mi>m</mi><mi>n</mi></mrow></msub><mo>&le;</mo><msub><mi>&mu;</mi><mi>c</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&mu;</mi><mi>c</mi></msub><mo>+</mo><msup><mrow><mo>(</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&mu;</mi><mi>c</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&mu;</mi><mrow><mi>m</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mn>0.5</mn></msup></mrow></mtd><mtd><mrow><msub><mi>&mu;</mi><mi>c</mi></msub><mo>&le;</mo><msub><mi>&mu;</mi><mrow><mi>m</mi><mi>n</mi></mrow></msub><mo>&le;</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000795794210000012.GIF" wi="1260" he="183" /></maths>其中,T<sub>1</sub>(μ<sub>mn</sub>)为低频信息的增强图像,μ<sub>c</sub>为方向波低频系数的隶属度;(3)对高频信息进行非线性抑噪增强处理:(3a)设置阈值T对高频信息进行噪声的抑制,其中,<img file="FDA0000795794210000013.GIF" wi="365" he="173" />σ为噪声方差,由中值估计器得到:<img file="FDA0000795794210000021.GIF" wi="543" he="143" />函数abs(·)表示绝对值,函数median(·)表示中值,H<sub>i</sub>为高频信息,i=1,2,3,4,5,6,7,N为高频系数个数,j为分解尺度;噪声抑制过程是当高频系数值大于或等于T时,该高频系数值不变;当高频系数值小于T时,该高频系数值取零;(3b)通过如下增益函数对抑制噪声后的高频信息进行增强处理:f(H′<sub>ij</sub>)=a[sigm(c(H′<sub>ij</sub>‑b))‑sigm(‑c(H′<sub>ij</sub>+b))]其中,f(H′<sub>ij</sub>)为高频信息的增强图像,H′<sub>ij</sub>为噪声抑制后的高频信息,i=1,2,3,j=1,2,3,4,5,6,7,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>a</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mi>s</mi><mi>i</mi><mi>g</mi><mi>m</mi><mrow><mo>(</mo><mi>c</mi><mo>(</mo><mn>1</mn><mo>-</mo><mi>b</mi><mo>)</mo></mrow><mo>)</mo><mo>-</mo><mi>s</mi><mi>i</mi><mi>g</mi><mi>m</mi><mrow><mo>(</mo><mo>-</mo><mi>c</mi><mo>(</mo><mn>1</mn><mo>+</mo><mi>b</mi><mo>)</mo></mrow><mo>)</mo></mrow></mfrac></mrow>]]></math><img file="FDA0000795794210000022.GIF" wi="843" he="158" /></maths>为中间值,b为控制整个高频信息增强范围的参数,取值为0&lt;b&lt;1,c为控制增强强度的参数,取值为0&lt;c&lt;1,sigm函数定义为:<img file="FDA0000795794210000023.GIF" wi="439" he="143" />y为函数sigm的输入值;(4)将增强后的低频信息和高频信息进行方向波逆变换,分别得到135度、90度及45度增强的图像;(5)将三个增强图像进行加权平均,获得最终的增强图像。
地址 710071 陕西省西安市太白南路2号