发明名称 结合可见光图像信息的SAR图像斑点噪声抑制方法
摘要 一种结合可见光图像信息的SAR图像斑点噪声抑制方法,其步骤为:(1)进行像素点区域属性分类,利用光学图像信息将像素点划分为均匀区域、缓变区域和边缘细节区域三类;(2)根据像素点所属类别,分别对SAR图像分类的斑点噪声进行抑制处理。本发明利用可见光图像将SAR图像分为均匀区域、缓变区域和边缘细节区域三类图像区域,然后在不同图像区域采用不同的噪声抑制运算;其优点在于:(1)利用不存在斑点噪声的可见光图像信息能够准确地对区域进行分类;(2)均匀区域和缓变区域采用了最适合的降斑运算,能够获取最好的降斑效果;(3)边缘细节区域采用了保持边缘处理,因此能够获取最好的边缘保持性能。
申请公布号 CN100544400C 申请公布日期 2009.09.23
申请号 CN200710120813.4 申请日期 2007.08.27
申请人 北京航空航天大学 发明人 徐华平;周芳
分类号 H04N5/217(2006.01)I;H04N5/14(2006.01)I 主分类号 H04N5/217(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 代理人 贾玉忠;卢 纪
主权项 1、一种结合可见光图像信息的SAR图像斑点噪声抑制方法,其特征在于步骤如下:(1)进行像素点区域属性分类,利用光学图像信息将像素点划分为均匀区域、缓变区域和边缘细节区域三类;(2)根据像素点所属类别,分别对SAR图像分类的斑点噪声进行抑制处理;所述步骤(1)中利用光学图像信息进行区域属性分类的步骤如下:a.对光学图像f(i,j)进行梯度运算得到梯度图f<sub>i</sub>(i,j)<maths num="0001"><![CDATA[<math><mrow><msub><mi>f</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msqrt><msup><mrow><mo>[</mo><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>[</mo><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup></msqrt></mrow></math>]]></maths>其中,各像素点的关系如表1所示; <tables num="0001"><table><tgroup cols="2"><colspec colname="c001" colwidth="50%" /><colspec colname="c002" colwidth="50%" /><thead><row><entry morerows="1">f(i,j)</entry><entry morerows="1">f(i,j+1)</entry></row></thead><tbody><row><entry morerows="1">f(i+1,j)</entry><entry morerows="1"></entry></row></tbody></tgroup></table></tables>b.根据两个梯度阈值—均匀阈值T<sub>u</sub>和边缘阈值T<sub>e</sub>,判定当前像素f(i,j)所属的区域类别,即<img file="C200710120813C00022.GIF" wi="901" he="224" />均匀阈值T<sub>u</sub>和边缘阈值T<sub>e</sub>根据先验知识获得;所述步骤(2)中对SAR图像分类进行噪声抑制处理的步骤为:以当前像素点g(i,j)为中心,取出大小为W×W的窗口R内所有像素点,根据光学图像分类结果进行判断:a.如果窗口内所有像素点属于均匀区域,则以窗口内所有像素点的局域均值作为当前像素点降噪后的取值,进行均值滤波,即<maths num="0002"><![CDATA[<math><mrow><mover><mi>g</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msup><mi>W</mi><mn>2</mn></msup></mfrac><munder><mi>&Sigma;</mi><mrow><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>&Element;</mo><mi>R</mi></mrow></munder><mi>g</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow></mrow></math>]]></maths>b.如果窗口内的像素点仅属于缓变区域或属于缓变和均匀区域,则以当前像素点散射特性的最小均方误差估计值作为当前像素点降噪后的取值,进行Lee滤波,即<maths num="0003"><![CDATA[<math><mrow><mover><mi>g</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>&mu;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><msup><mi>&alpha;</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>&beta;</mi><mn>2</mn></msup></mrow><mrow><msup><mi>&alpha;</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><msup><mi>&beta;</mi><mn>4</mn></msup></mrow></mfrac><mo>&CenterDot;</mo><mrow><mo></mo><mo>[</mo><mi>g</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mi>&mu;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo><mo></mo></mrow></mrow></math>]]></maths>其中,<maths num="0004"><![CDATA[<math><mrow><mi>&beta;</mi><mo>=</mo><mfrac><mn>1</mn><msqrt><mi>L</mi></msqrt></mfrac><mo>,</mo></mrow></math>]]></maths>L为图像视数,由所需的噪声抑制效果确定,<maths num="0005"><![CDATA[<math><mrow><mi>&alpha;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>&sigma;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mrow><mi>&mu;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mi>&mu;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msup><mi>W</mi><mn>2</mn></msup></mfrac><munder><mi>&Sigma;</mi><mrow><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>&Element;</mo><mi>R</mi></mrow></munder><mi>g</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msup><mi>&sigma;</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msup><mi>W</mi><mn>2</mn></msup><mo>-</mo><mn>1</mn></mrow></mfrac><munder><mi>&Sigma;</mi><mrow><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>&Element;</mo><mi>R</mi></mrow></munder><msup><mrow><mrow><mo>[</mo><mi>g</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>-</mo><mi>&mu;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo></mrow></mrow><mn>2</mn></msup><mo>;</mo></mrow></math>]]></maths>c.如果窗口内的像素点有属于边缘细节区域的,则保留当前像素点的原值。
地址 100083北京市海淀区学院路37号
您可能感兴趣的专利