主权项 |
一种基于方向波变换和模糊理论的图像增强方法,包括如下步骤:(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>μ</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>μ</mi><mrow><mi>m</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>μ</mi><mi>c</mi></msub><mo>-</mo><msup><mrow><mo>(</mo><msup><msub><mi>μ</mi><mi>c</mi></msub><mn>2</mn></msup><mo>-</mo><msubsup><mi>μ</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>≤</mo><msub><mi>μ</mi><mrow><mi>m</mi><mi>n</mi></mrow></msub><mo>≤</mo><msub><mi>μ</mi><mi>c</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>μ</mi><mi>c</mi></msub><mo>+</mo><msup><mrow><mo>(</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>μ</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>μ</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>μ</mi><mi>c</mi></msub><mo>≤</mo><msub><mi>μ</mi><mrow><mi>m</mi><mi>n</mi></mrow></msub><mo>≤</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<b<1,c为控制增强强度的参数,取值为0<c<1,sigm函数定义为:<img file="FDA0000795794210000023.GIF" wi="439" he="143" />y为函数sigm的输入值;(4)将增强后的低频信息和高频信息进行方向波逆变换,分别得到135度、90度及45度增强的图像;(5)将三个增强图像进行加权平均,获得最终的增强图像。 |