CN102136134A - 基于mrf先验的sar图像去斑方法 - Google Patents
基于mrf先验的sar图像去斑方法 Download PDFInfo
- Publication number
- CN102136134A CN102136134A CN 201110059882 CN201110059882A CN102136134A CN 102136134 A CN102136134 A CN 102136134A CN 201110059882 CN201110059882 CN 201110059882 CN 201110059882 A CN201110059882 A CN 201110059882A CN 102136134 A CN102136134 A CN 102136134A
- Authority
- CN
- China
- Prior art keywords
- sigma
- sar image
- pixel
- similarity
- priori
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Image Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种图像处理技术领域基于MRF先验的SAR图像去斑方法,包括如下步骤:1.输入一副待去斑的SAR图像;2.选取中心像素块;3.确定搜索窗;4.选取相似块;5.计算相似性权值;6.遍历搜索窗;7.计算待估计像素点的恢复值;8.遍历整幅SAR图像;9.恢复整幅SAR图像;10.迭代去斑;11.获得最终去斑结果。本发明充分利用了图像的先验信息,相对于目前空域滤波的去斑方法,能够更好的平滑斑点噪声,同时保持SAR图像的边缘和纹理细节,可用于对SAR图像的去斑处理。
Description
技术领域
本发明属于图像处理技术领域,更进一步涉及一种雷达图像处理技术的滤波方法,该方法使用马尔科夫(MRF)随机场提取先验信息,以贝叶斯非局部均值(BNL)滤波为去斑方法,可用于对合成孔径雷达(SAR)图像进行去斑处理。
背景技术
目前,基于空域滤波的去斑效果较好的算法有原始非局部均值方法,改进sigma滤波方法和PPB滤波方法。
Buades等人在“Buades,B.Coll,and J.M.Morel.A non-local algorithm for image denoising.Proc.IEEE CVPR,2005,2:60-65.”中提出了一种非局部均值去斑滤波方法。该方法不是用图像中单个像素的灰度值进行比较,而是对该像素周围的窗口中整个灰度的分布状况进行比较,根据灰度分布的相似性贡献权值。实践证明,使用该方法获得的灰度分布相似性权值不准确,导致图像结构信息如边缘、线性体、点等目标会在一定程度上被模糊或滤除。
J.S.Lee等人在“J.S.Lee,J.H.Wen,T.L.Ainsworth,K.S.Chen,and A.J.Chen,“Improved sigma filter for speckle filtering of SAR imagery,”IEEE Trans.Geosci.Remote Sens.,vol.47,no.1,Jan 2009,202-213.”中提出了改进sigma去斑滤波方法。该方法使用新的sigma范围选取窗口内的像素点,然后对选取的像素点采用最小均方误差(MMSE)估计中心点像素的值。实际使用中证明该方法在一定程度上压缩了强反射点目标的亮度,且不能很好平滑同质区域的噪声,导致去斑效果不理想。
C.A.Deledalle等人在“C.A.Deledalle,L.Denis,and F.Tupin.Iterative weighted maximum likelihood denoising with probabilistic patch-based weights.IEEE Trans.Image Processing,vol.18,no.12,Dec.2009,2661-2672.”中提出了一种在最大似然框架下迭代去斑的滤波方法。该方法在最大似然估计框架下得到加权平均公式,并定义块权值的计算公式,通过迭代方法逐步修改先验信息,最终收敛至最佳去斑结果。此方法在使用中需迭代25次,因此其复杂度很高,降低了该方法的实用性。
发明内容
本发明的目的在于克服现有技术的不足,先使用MRF模型提取图像的先验信息,再根据贝叶斯非局部均值滤波框架对图像去斑的方法,提出了一种基于MRF先验的BNL均值滤波器对SAR图像空域去斑方法,以实现对SAR图像降斑中边缘和平滑区域的兼顾,提高图像去斑效果。
本发明的具体步骤包括如下:
(1)输入一副待去斑的SAR图像;
(2)选取中心像素块:在待去斑的SAR图像中,逐行扫描选取一个像素点为待估计像素点,以待估计像素点为中心,以固定长度为块半径,选取一个正方形的中心像素块;
(3)确定搜索窗:以选取的待估计像素点为中心,以固定长度为搜索窗半径,选取一个正方形的搜索窗;
(4)选取相似块:在搜索窗中逐行扫描选取一个像素点作为相似点,以此相似点为中心,选取一个与中心像素块等大小的像素块作为相似块;
(5)计算相似性权值:
5a)计算相似性条件概率:依据条件概率密度函数,计算在步骤(4)所选取的相似块的条件下,中心像素块出现的条件概率;
5b)根据基于MRF先验的相似性度量公式,计算先验相似性:
其中,g(·)为先验相似性;β为平滑参数;ηs,k和ηt,k分别为us,k和ut,k的四邻域ui为ηs,k中的像素点;uj为ηt,k中的像素点
5c)将相似性条件概率与先验相似性相乘获得相似性权值;
(6)遍历搜索窗:判断搜索窗内的点是否搜索完,若搜索完则进行下一步骤,若没有搜索完则重复步骤(4)、(5),计算出整个搜索窗中所有的相似块与中心块的相似性权值;
(7)计算待估计像素点的恢复值:根据步骤(6)中求出的相似性权值,将获得的每个相似块的相似性权值和其对应的相似点的灰度值,加权求和得出待估计像素点的恢复值;
(8)遍历整幅SAR图像:判断整幅SAR图像中的点是否搜索完,若搜索完则进行下一步骤,若没有搜索完则重复步骤(2)、(3)、(4)、(5)、(6)、(7),计算出整幅SAR图像中所有点的恢复值;
(9)恢复整幅SAR图像:根据步骤(8)中获得的恢复值,恢复出整幅SAR图像;
(10)迭代去斑:使用最新恢复出的SAR图像计算相似性权值,将获得的相似性权值加权求和恢复出整幅图像;
(11)获得最终去斑结果:重复执行步骤(10),直到去斑图达到预想要求,中止迭代。
本发明与现有的技术相比具有以下优点:
1.本发明是在空域中直接对像素点进行处理的,实现过程简单,且可以并行实现。
2.本发明由于利用BNL均值滤波方法和SAR图像斑点模型以及MRF先验模型推导得到的一种新的相似性权值,能够准确地计算出SAR图像中以待估计像素点为中心的中心块和以搜寻窗内的相似点为中心的相似块之间的相似性,进而能准确的恢复出待估计像素点,较好的保持了图像中的线目标和点目标,并能很好的平滑同质区域的噪声,从而提高了图像的去斑效果。
3.本发明由于充分利用了图像中的先验信息,因此能够准确的计算相似性权值,从而能够准确计算待估计像素点的恢复值,因此,能够较好平滑斑点噪声,同时保持和恢复SAR图像的边缘和纹理细节,图像的去斑效果明显增强。
附图说明
图1为本发明的流程图;
图2为基于空域滤波的去斑方法的效果图;
图3为本发明的效果图。
具体实施方式
下面结合附图对本发明做进一步的描述。
参照图1,本发明的具体实施步骤如下:
步骤1,输入一副待去斑的SAR图像。
步骤2,选取中心像素块。
选取中心像素块主要是为了确定待估计像素点的领域信息。
1a)选取待估计像素点。在待去斑的SAR图像v中,逐行扫描选取一个像素点t作为待估计像素点;
1b)确定中心像素块。以待估计像素点为中心,以固定长度2~5像素点为块半径,选取一个大小为M×M正方形的中心像素块。本发明实施例中以像素点t为中心,选取一个大小为7×7的中心像素块vt。
步骤3,确定搜索窗。
由于搜索窗的大小直接决定了相似点的选取范围,它必须包含较多的相似点。确定搜索窗时,以选取的待估计像素点t为中心,以固定长度7~15像素点为搜索窗半径,选取大小为N×N的搜索窗。本发明实施例中选取大小为21×21的搜索窗。
步骤4,选取相似块。
在步骤2获得的搜索窗中逐行扫描选取一个像素点s作为相似点,以此相似点为中心,选取一个大小为M×M相似块us。相似块的大小一定要与中心相似块的大小相同,它描述了相似点s的领域信息。
步骤5,计算相似性权值。
5a)计算相似性条件概率。根据条件概率密度函数,计算在已知相似块us的条件下,中心像素vt块出现的条件概率。相似性条件概率部分决定了搜索窗中的相似点在加权平均中对恢复待估计相似点的贡献,根据强度和幅度的SAR图像,它有不同的计算公式。
适用于强度SAR图像的条件概率密度函数为:
适用于幅度SAR图像的条件概率密度函数为:
vt,k为vt中逐行扫描的第k个像素点
us,k为us中逐行扫描的第k个像素点
L为SAR图像视数
5b)计算先验相似性。先验相似性主要是充分利用了相似块与中心像素块之间的先验信息,根据它们之间的先验信息来计算它们之间的相似度。
计算先验相似性需使用先验信息,其相似性度量公式如下:
其中,g(·)为先验相似性;
β为平滑参数
ηs,k和ηt,k分别为us,k和ut,k的四邻域
ui为ηs,k中的像素点
uj为ηt,k中的像素点
5c)将相似性条件概率与先验相似性相乘获得相似性权值w(t,s)。相似性权值描述了相似点与待估计像素点之间的相似度。根据强度和幅度的SAR图像,它有不同的计算公式。
适用于强度SAR图像的相似性权值公式为:
其中,Z(t)为归一化因子:
适用于幅度SAR图像的相似性权值公式为:
其中,Z(t)为归一化因子:
步骤6,遍历搜索窗。判断搜索窗内的点是否搜索完,若搜索完则进行下一步骤,若没有搜索完则重复步骤4、5,计算出整个搜索窗中所有的相似块与中心块的相似性权值。
步骤7,计算待估计像素点的恢复值。根据步骤6中求出的相似性权值,将获得的每个相似块的相似性权值和其对应的相似点的灰度值,进行加权求和得出待估计像素点的恢复值。
其中,MBNLv(t)为加权平均后t点的恢复值
v(s)为原图中s处的灰度值
w(t,s)为像素点s相对于t的权值,且满足:0≤w(t,s)≤1,
步骤8,遍历整幅SAR图像。判断整幅SAR图像中的点是否搜索完,若搜索完则进行下一步骤,若没有搜索完则重复步骤2、3、4、5、6、7,计算出整幅SAR图像中所有点的恢复值。
步骤9,恢复整幅SAR图像:根据步骤8中获得的恢复值,获得整幅SAR图像的第一次恢复值u。
步骤10,迭代去斑。使用最新恢复出的SAR图像u,重复步骤2、3、4、5、6、7、8、9,恢复出整幅图像的更新值。
步骤11,获得最终去斑结果:重复执行步骤10,直到去斑图达到预想要求,中止迭代。
下面结合图2、3对本发明的效果做进一步的描述
图2(a)和(b)表示原始非局部均值滤波方法对两幅SAR图像的去斑效果图,图2(c)和(d)表示改进sigma滤波方法对两幅SAR图像的去斑效果图,图2(e)和(f)表示PPB滤波方法对两幅SAR图像的去斑效果图。
图3(a)和(b)表示输入待去斑的SAR图像,其中图3(a)表示待去斑的两视幅度SAR图像,图3(b)表示待去斑的四视强度SAR图像。图3(a)和(b)中,标号1、2、3、4和5标示出了待评价去斑效果的同质区域。图3(c)和(d)表示本发明两次迭代对两幅真实SAR图像去斑的效果图,为了证明我们滤波方法有较好的去斑效果,分别使用窗口大小为5×5的原始非局部均值滤波方法、改进sigma滤波方法和PPB滤波方法与本章滤波方法实验结果进行对比。将去斑后区域1、区域2、区域3、区域4和区域5对应的均值、标准差、等效视数ENL作为评价指标,比较各种空域滤波方法的优劣。其比较结果如表1所示:
表1各种去斑结果对比
从表1可以发现,本发明在均值、标准差和等效视数各方面都有较佳的表现,本方法的去斑效果显示出了同质区域斑点抑制能力和一定的保护点线等强反射目标的能力。同时,本发明与基于边缘、同质区域、点目标分类的空域滤波方法相比,在对细节信息,像线和纹理进行恢复时比其它方法表现出了更好的效果。本发明最明显的优点在于仅仅利用了图像中的先验信息,也能把点目标的亮度保持的相当好。从图3(c)和(d)中可以看出,本发明对同质区域的平滑性能很好,因此可能会因为过平滑而丢失细小纹理信息,但是还是能把大片区域的地物轮廓恢复出来。另外,如果对细小纹理感兴趣可以减少迭代次数,降低对同质区域的平滑力,就能保护好细小纹理结构信息。相对于其它的一些空域的去斑方法具有更好的性能,能够更好平滑斑点噪声同时保持SAR图像的边缘和纹理细节。
Claims (7)
1.一种基于MRF先验的SAR图像去斑方法,包括如下步骤:
(1)输入一副待去斑的SAR图像;
(2)选取中心像素块:在待去斑的SAR图像中,逐行扫描选取一个像素点为待估计像素点,以待估计像素点为中心,以固定长度为块半径,选取一个正方形的中心像素块;
(3)确定搜索窗:以选取的待估计像素点为中心,以固定长度为搜索窗半径,选取一个正方形的搜索窗;
(4)选取相似块:在搜索窗中逐行扫描选取一个像素点作为相似点,以此相似点为中心,选取一个与中心像素块等大小的像素块作为相似块;
(5)计算相似性权值:
5a)计算相似性条件概率:依据条件概率密度函数,计算在步骤(4)所选取的相似块的条件下,中心像素块出现的条件概率;
5b)根据基于MRF先验的相似性度量公式,计算先验相似性:
其中,g(·)为先验相似性;β为平滑参数;ηs,k和ηt,k分别为us,k和ut,k的四邻域ui为ηs,k中的像素点;uj为ηt,k中的像素点
5c)将相似性条件概率与先验相似性相乘获得相似性权值;
(6)遍历搜索窗:判断搜索窗内的点是否搜索完,若搜索完则进行下一步骤,若没有搜索完则重复步骤(4)、(5),计算出整个搜索窗中所有的相似块与中心块的相似性权值;
(7)计算待估计像素点的恢复值:根据步骤(6)中求出的相似性权值,将获得的每个相似块的相似性权值和其对应的相似点的灰度值,加权求和得出待估计像素点的恢复值;
(8)遍历整幅SAR图像:判断整幅SAR图像中的点是否搜索完,若搜索完则进行下一步骤,若没有搜索完则重复步骤(2)、(3)、(4)、(5)、(6)、(7),计算出整幅SAR图像中所有点的恢复值;
(9)恢复整幅SAR图像:根据步骤(8)中获得的恢复值,恢复出整幅SAR图像;
(10)迭代去斑:使用最新恢复出的SAR图像计算相似性权值,重复步骤2、3、4、5、6、7、8、9,恢复出整幅图像的更新值;
(11)获得最终去斑结果:重复执行步骤(10),直到去斑图达到预想要求,中止迭代。
2.根据权利要求1所述的一种基于MRF先验的SAR图像去斑方法,其特征在于:所述步骤(2)中的固定长度为2~5个像素点。
3.根据权利要求1所述的一种基于MRF先验的SAR图像去斑方法,其特征在于:所述步骤(3)中的固定长度为7~15个像素点。
6.根据权利要求1所述的一种基于MRF先验的SAR图像去斑方法,其特征在于:所述步骤5c)中适用于强度SAR图像的相似性权值公式为:
其中,Z(t)为归一化因子:
7.根据权利要求1所述的一种基于MRF先验的SAR图像去斑方法,其特征在于:所述步骤5c)中适用于幅度SAR图像的相似性权值公式为:
其中,Z(t)为归一化因子:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110059882 CN102136134B (zh) | 2011-03-11 | 2011-03-11 | 基于mrf先验的sar图像去斑方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110059882 CN102136134B (zh) | 2011-03-11 | 2011-03-11 | 基于mrf先验的sar图像去斑方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102136134A true CN102136134A (zh) | 2011-07-27 |
CN102136134B CN102136134B (zh) | 2012-12-26 |
Family
ID=44295909
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110059882 Expired - Fee Related CN102136134B (zh) | 2011-03-11 | 2011-03-11 | 基于mrf先验的sar图像去斑方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102136134B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102750675A (zh) * | 2012-05-31 | 2012-10-24 | 华中科技大学 | 一种斑点噪声污染图像的非局部均值滤波方法 |
CN103077499A (zh) * | 2013-01-09 | 2013-05-01 | 西安电子科技大学 | 基于相似块的sar图像预处理方法 |
CN104715458A (zh) * | 2015-03-23 | 2015-06-17 | 华中科技大学 | 一种双模非局部均值滤波方法 |
CN105761225A (zh) * | 2015-11-20 | 2016-07-13 | 广西慧云信息技术有限公司 | 一种优化图片显示效果的方法 |
CN116993626A (zh) * | 2023-09-26 | 2023-11-03 | 成都市晶林科技有限公司 | 一种基于时空域的红外图像降噪方法及*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101398487A (zh) * | 2007-09-28 | 2009-04-01 | 北京师范大学 | 一种降低合成孔径雷达影像斑点噪声的方法 |
EP2146315A1 (en) * | 2008-07-16 | 2010-01-20 | Galileian Plus s.r.l. | Method of filtering SAR images from speckle noise and related device. |
CN101639934A (zh) * | 2009-09-04 | 2010-02-03 | 西安电子科技大学 | 基于轮廓波域块隐马尔可夫模型sar图像去噪方法 |
CN101685158A (zh) * | 2009-07-02 | 2010-03-31 | 西安电子科技大学 | 基于隐马尔科夫树模型的sar图像去噪方法 |
-
2011
- 2011-03-11 CN CN 201110059882 patent/CN102136134B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101398487A (zh) * | 2007-09-28 | 2009-04-01 | 北京师范大学 | 一种降低合成孔径雷达影像斑点噪声的方法 |
EP2146315A1 (en) * | 2008-07-16 | 2010-01-20 | Galileian Plus s.r.l. | Method of filtering SAR images from speckle noise and related device. |
CN101685158A (zh) * | 2009-07-02 | 2010-03-31 | 西安电子科技大学 | 基于隐马尔科夫树模型的sar图像去噪方法 |
CN101639934A (zh) * | 2009-09-04 | 2010-02-03 | 西安电子科技大学 | 基于轮廓波域块隐马尔可夫模型sar图像去噪方法 |
Non-Patent Citations (2)
Title |
---|
《电子与信息学报》 20090331 宋珩,王世唏,计科峰,郁文贤 基于结构保持MRF模型的SAR图像去斑 第745-747页 1-7 第31卷, 第3期 * |
《电波科学学报》 20070430 吴艳,王霞,廖桂生 基于小波域隐马尔可夫混合模型的SAR图像降斑算法 第244-250页 1-7 第22卷, 第2期 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102750675A (zh) * | 2012-05-31 | 2012-10-24 | 华中科技大学 | 一种斑点噪声污染图像的非局部均值滤波方法 |
CN102750675B (zh) * | 2012-05-31 | 2014-08-27 | 华中科技大学 | 一种斑点噪声污染图像的非局部均值滤波方法 |
CN103077499A (zh) * | 2013-01-09 | 2013-05-01 | 西安电子科技大学 | 基于相似块的sar图像预处理方法 |
CN103077499B (zh) * | 2013-01-09 | 2015-04-08 | 西安电子科技大学 | 基于相似块的sar图像预处理方法 |
CN104715458A (zh) * | 2015-03-23 | 2015-06-17 | 华中科技大学 | 一种双模非局部均值滤波方法 |
CN104715458B (zh) * | 2015-03-23 | 2017-08-01 | 华中科技大学 | 一种双模非局部均值滤波方法 |
CN105761225A (zh) * | 2015-11-20 | 2016-07-13 | 广西慧云信息技术有限公司 | 一种优化图片显示效果的方法 |
CN116993626A (zh) * | 2023-09-26 | 2023-11-03 | 成都市晶林科技有限公司 | 一种基于时空域的红外图像降噪方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN102136134B (zh) | 2012-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103093441B (zh) | 基于变换域的非局部均值和双变量模型的图像去噪方法 | |
CN101661611B (zh) | 基于贝叶斯非局部均值滤波器的实现方法 | |
CN103077508B (zh) | 基于变换域非局部和最小均方误差的sar图像去噪方法 | |
CN101950414B (zh) | 自然图像非局部均值去噪方法 | |
CN104408700A (zh) | 基于形态学和pca的轮廓波红外与可见光图像融合方法 | |
CN102136134B (zh) | 基于mrf先验的sar图像去斑方法 | |
CN104103041B (zh) | 超声图像混合噪声自适应抑制方法 | |
CN109509163B (zh) | 一种基于fgf的多聚焦图像融合方法及*** | |
CN105894476A (zh) | 基于字典学习融合的sar图像降噪处理方法 | |
Teng et al. | Modified pyramid dual tree direction filter‐based image denoising via curvature scale and nonlocal mean multigrade remnant filter | |
CN103473755B (zh) | 基于变化检测的sar图像稀疏去噪方法 | |
CN103020918A (zh) | 基于形状自适应邻域均值的非局部均值去噪方法 | |
CN104657951A (zh) | 图像乘性噪声移除方法 | |
CN108460765B (zh) | 一种虹膜图像质量检测方法 | |
CN103077507A (zh) | 基于Beta算法的多尺度SAR图像降噪方法 | |
CN102222327A (zh) | 基于Treelet变换和最小均方误差估计的图像去噪方法 | |
CN103793889B (zh) | 基于字典学习和ppb算法的sar图像去斑方法 | |
CN104537624B (zh) | 基于ssim校正聚类稀疏表示的sar图像降斑方法 | |
CN103971345A (zh) | 一种基于改进双边滤波的图像去噪方法 | |
CN105469358A (zh) | 一种图像处理方法 | |
CN102314675A (zh) | 基于小波高频的贝叶斯去噪方法 | |
Wu et al. | Research on crack detection algorithm of asphalt pavement | |
Sun et al. | Color image denoising based on guided filter and adaptive wavelet threshold | |
CN107085839A (zh) | 基于纹理增强与稀疏编码的sar图像降斑方法 | |
CN103778615A (zh) | 基于区域相似性的多聚焦图像融合方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121226 Termination date: 20180311 |