CN106897985B - 一种基于可见度分类的多角度sar图像融合方法 - Google Patents
一种基于可见度分类的多角度sar图像融合方法 Download PDFInfo
- Publication number
- CN106897985B CN106897985B CN201710042610.1A CN201710042610A CN106897985B CN 106897985 B CN106897985 B CN 106897985B CN 201710042610 A CN201710042610 A CN 201710042610A CN 106897985 B CN106897985 B CN 106897985B
- Authority
- CN
- China
- Prior art keywords
- pixel
- visibility
- image
- angle
- window
- 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.)
- Active
Links
- 238000007500 overflow downdraw method Methods 0.000 title claims abstract description 35
- 238000003384 imaging method Methods 0.000 claims abstract description 105
- 230000004927 fusion Effects 0.000 claims abstract description 59
- 230000011218 segmentation Effects 0.000 claims abstract description 21
- 238000000034 method Methods 0.000 claims description 50
- 238000012545 processing Methods 0.000 claims description 34
- 230000001427 coherent effect Effects 0.000 claims description 29
- 238000001914 filtration Methods 0.000 claims description 21
- 238000010586 diagram Methods 0.000 claims description 17
- 238000006731 degradation reaction Methods 0.000 claims description 12
- 230000015556 catabolic process Effects 0.000 claims description 11
- 230000000873 masking effect Effects 0.000 claims description 10
- 238000009825 accumulation Methods 0.000 claims description 7
- 230000000007 visual effect Effects 0.000 claims description 6
- 238000002156 mixing Methods 0.000 abstract 1
- 238000001514 detection method Methods 0.000 description 37
- 230000000694 effects Effects 0.000 description 9
- 238000011160 research Methods 0.000 description 9
- 238000007499 fusion processing Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 5
- 238000003706 image smoothing Methods 0.000 description 4
- 230000002411 adverse Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 230000001094 effect on targets Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明提供一种全新的多角度SAR图像融合成像算法——一种基于可见度分类的多角度SAR图像融合方法。本发明根据图像像素的叠掩区与非叠掩区二值分割结果,计算各方位视角图像的非叠掩区像素Inon_overlap的多视角可见度LI,通过可见度指标自动实现图像像素的精细分类;此外本发明区分了不同视角观测结果对图像像素的贡献,且仅利用有贡献的角度进行融合,避免了简单融合方法不能获得最优融合结果的缺点;同时本发明原理简单,融合质量佳,且可见度指标计算简单,分类精度高,融合算法工程实现简单。
Description
技术领域
本发明属于成像雷达探测领域,同时还属于图像融合处理领域,尤其涉及一种基于可见度分类的多角度SAR图像融合方法。
背景技术
合成孔径雷达(Synthetic Aperture Radar,以下简称SAR)作为主动式微波传感器,充分利用了脉冲压缩技术、合成孔径原理以及信号相干处理方法,能够以真实小孔径天线获取距离向和方位向二维高分辨的雷达图像,可实现全天时、全天候的目标观测,可透过地表和植被获取目标信息,可跟踪测量目标的散射特性,还可实现对机动目标的成像识别。这些特点使其一跃成为现代军事侦察的主要工具,其获取的高分辨雷达图像极大地推进了军事活动的信息化进程并深刻改变了现代战争的形态。
但是,经典SAR的直线观测几何仅能获取目标在较小方位角度内的散射信息,其“距离-多普勒”成像原理以及特殊的侧视成像探测模式,导致输出图像不可避免地存在目标遮挡、探测阴影(盲区)和目标叠掩等固有问题,给系统对目标检测与识别带来了不利影响。例如,SAR图像中由于遮挡造成的阴影区域成为探测的盲区,无法提供该部分的目标信息,而叠掩现象使得不同空间位置的目标散射信息相互叠加在同一个分辨单元内,严重时甚至会对感兴趣的目标散射信息产生遮盖,导致系统对目标检测和识别的正确率大大降低。
在军事侦察中,准确而全面地了解战场环境和获取重点目标的信息非常重要,上述问题让成像效果大打折扣,大大降低了SAR在城市、山地等复杂场景下成像结果的实用性和可读性,严重制约了SAR的探测性能。为进一步提升SAR在城市、山地等复杂场景下的探测性能及应用效能,提高SAR全面、准确获取目标信息的能力,多角度SAR成像技术应运而生。它通过融合多个角度下的探测信息,可以实现对目标电磁散射特征的完整描述,有效拓展被探测目标的空间谱支撑区,一定程度上克服传统SAR直线观测中存在的目标遮挡、探测阴影(盲区)和目标叠掩等问题,对于战场军事侦察和城市测绘具有重要的现实意义和极高的实际应用价值。
作为当前SAR成像领域的热点研究方向之一,国内外众多研究机构和学者已成功开展了一系列基于星载、机载和地基平台的多角度SAR成像理论和试验研究,并获得了一批良好的成像探测结果,充分展示了多角度SAR技术在对地观测方面的巨大应用潜力。但是,目前的多角度SAR成像技术研究所公开发表的成像探测结果,尚未解决山地、城市等地形起伏剧烈,雷达图像叠掩概率高、探测盲区大的成像场景的有效处理。
应用于未来战争环境下的合成孔径雷达不仅要求成像性能好,而且要求叠掩概率低、探测盲区小;现有多角度SAR成像处理技术在减小成像叠掩和探测盲区方面能力不足,如何提高多角度SAR成像处理技术消除成像叠掩和探测盲区的能力已成为重大现实问题。
本质上讲,多角度SAR成像需对多个角度的目标信息进行融合处理。根据融合处理方式的不同,多角度SAR成像理论主要可分为图像层融合的成像理论和信号层融合的成像理论。其中,图像层融合的多角度SAR成像理论主要在图像域实现目标信息融合,以图像配准、图像特征提取和图像融合方法等关键问题为研究重点。研究结果表明,利用不同角度的SAR图像可以提高目标参数的估计精度,并在一定程度上克服遮挡和多次散射对目标反演的影响,从原理上证明了多角度SAR成像的正确性。目前基于图像层的多角度SAR成像理论研究更为普遍,发展也相对成熟。
信号层融合的多角度SAR成像理论的关键问题是如何实现不同观测角度回波数据的融合,以及如何对融合后的数据进行统一的成像处理。由于多角度SAR观测角度差异较大,各角度SAR数据的空间不连续性和空间采样的非均匀性对基于信号层融合的成像处理带来了很大的挑战,传统基于傅里叶变换的方法会导致主瓣展宽、旁瓣升高等散焦问题,因此必须开展符合上述数据特点的特殊成像处理算法。同时,受目标散射特性非各向同性的影响,基于经典理论的信号层融合多角度SAR成像理论适用范围有限,亟待开发和完善基于新的信号处理理论的信号层融合多角度SAR成像理论及算法。
图像层融合的多角度SAR成像算法利用现有SAR成像算法获取目标在多个角度上的成像结果,通过图像融合技术,实现对目标信息的融合。该类算法可结合现有SAR成像理论和图像融合领域的研究成果,对系统成像条件要求较低,具有较强的兼容性和鲁棒性。但是,如何实现对目标特征信息的有效提取和不同角度下目标图像信息的准确融合仍然是该算法需要重点研究的问题。目前,在城市高精度测绘问题中,为克服单一视角SAR成像中建筑物遮挡和叠掩对高价值目标(如道路、桥梁、车辆等)探测的影响,多使用图像层融合算法。
信号层融合的多角度SAR成像算法首先在信号层对回波信号进行融合,然后再对融合后的数据进行成像处理,最终获得目标的成像结果。由于此类算法对接收数据的相干性要求较高,期望通过全相干累积的方法获得目标的最优分辨率和高信噪比图像,但是由于实际目标并不满足各向同性散射条件,这一目标实现难度很大。相比图像层融合的多角度SAR成像,其对系统的观测几何、平台稳定性和信号的采样要求更高。目前,信号层融合处理的成像算法仍是多角度SAR成像算法研究的一个难点。
可见,针对多角度SAR融合成像算法的研究,更适于综合利用信号层融合的相干处理(各视角或各子孔径单独相干成像)以及图像层融合的非相干处理(不同视角或不同子孔径非相干融合)来获得更优、更准确的成像探测结果。
发明内容
为解决上述多角度SAR融合成像算法在减小成像叠掩和探测盲区方面能力不足的问题,本发明提供一种全新的多角度SAR图像融合成像算法——一种基于可见度分类的多角度SAR图像融合方法。本发明适用于对山地、城市等高地形起伏的目标区域实施合成孔径成像探测,可有效消除单一视角SAR成像探测结果中叠掩、阴影等现象,克服其对目标检测与识别的不利影响,确保复杂场景下对敌军事目标的成像侦察能力不受影响。经该算法融合处理的图像,大大提升了图像的可读性和可懂度,非常有利于图像应用时的目标判读和理解。
一种基于可见度分类的多角度SAR图像融合方法,其特征在于,包括以下步骤:
步骤1:确定成像平面,并采用时域快速后向投影算法在与成像平面相对应的成像网格上完成N个单一视角SAR数据的相干累积处理,得到N幅单一视角SAR图像,其中N至少为3;
步骤2:对各幅单一视角SAR图像实施去相干斑滤波处理;
步骤3:对各幅相干斑滤波处理后的单一视角SAR图像分别进行大尺度降质处理,得到N幅粗分辨率的降质图像Ic;
步骤4:采用滑动窗口法,以降质图像Ic各像素为中心,计算滑动窗口内的图像能量,得到能量图Ec;
步骤5:根据能量图Ec和单一视角SAR图像的统计分布特性,确定各幅单一视角SAR图像的分割阈值Lc并完成二值分割;其中,如果任一像素点的像素值大于Lc,则该像素属于叠掩区,否则属于非叠掩区;
步骤6:计算各单一视角SAR图像的非叠掩区像素Inon_overlap的可见度LI,并依据可见度LI进行非叠掩区像素的精细分类;具体过程是:
假设可见度LI初始值为0,且处于非叠掩区的像素为可见,则如果成像平面任一点处于任一方位的单一视角SAR图像的非叠掩区,则该像素点的可见度LI增加1;其中,可见度取值为0到N;统计成像平面上各像素的可见度LI,按可见度LI的大小将具有N个不同的方位视角的成像平面上各像素划分为N+1个可见度等级;
步骤7:对可见度等级相同的各个像素进行融合,得到最终输出的融合图像Iopt。
一种基于可见度分类的多角度SAR图像融合方法,步骤7所述的对可见度等级相同的各个像素进行融合的方法为线性均值融合,具体为:分别计算各个可见度等级中包含的所有像素点像素值的平均值,取各等级中的平均值作为融合图像Iopt相同可见度等级包含像素点的像素值。
一种基于可见度分类的多角度SAR图像融合方法,步骤7所述的平均值采用滑动窗口法计算:
其中,为图像的平均值,I(i,j)为像素点的位置,nw为矩形滑动窗口内总的像素数,i=1,2,3,......nw;j=1,2,3,......nw。
一种基于可见度分类的多角度SAR图像融合方法,步骤7所述的对可见度等级相同的各个像素进行融合的方法为最大对比度融合,具体为:分别计算各个可见度等级中包含的所有像素点的对比度,取各等级中对比度最大的像素点对应的像素值作为融合图像Iopt相同可见度等级包含像素点的像素值。
一种基于可见度分类的多角度SAR图像融合方法,步骤7所述的对比度采用滑动窗口法计算:
Imax=max[I(i,j)],i=1,2,3,......nw;j=1,2,3,......nw
Imin=min[I(i,j)],i=1,2,3,......nw;j=1,2,3,......nw
为滑动窗口内的图像I的局部平均值,Imax和Imin分别为各个可见度等级中包含的所有像素值的最大值和最小值,nw为矩形滑动窗口内总的像素数;I(i,j)为像素点的位置。
一种基于可见度分类的多角度SAR图像融合方法,步骤7所述的对可见度等级相同的各个像素进行融合的方法为最小图像熵融合,具体为:分别计算各个可见度等级中包含的所有像素点的图像熵,取各等级中图像熵最小的像素点对应的像素值作为融合图像Iopt相同可见度等级包含像素点的像素值。
一种基于可见度分类的多角度SAR图像融合方法,步骤7所述的图像熵采用滑动窗口法计算:
P(n)为各图像像素取值出现的概率,nw为矩形滑动窗口内总的像素数,L=nw 2为矩形滑动窗口内的总像素数,I(i,j)为像素点的位置,j=1,2,3,......nw。
一种基于可见度分类的多角度SAR图像融合方法,步骤4所述的采用滑动窗口法,以各像素为中心计算窗口内图像能量,得到能量图Ec具体方法为:
假设滑动窗口W的尺寸为nw×nw,其表达式为:
n=[0 1 2......nw]T
式中α为0~1之间的窗口参数;
滑动窗口W在降质图像Ic上逐像素移动,以每个像素为中心计算窗口内的能量值,并将该能量值作为能量图对应的窗口中心像素点的像素值,直至滑动窗口W的中心遍历所有像素,得到能量图Ec:
其中,nw为奇数。
一种基于非可见度分类的多角度SAR图像融合方法,所述窗口为矩形窗、汉宁窗、海明窗或布拉克曼窗。
一种基于可见度分类的多角度SAR图像融合方法,步骤5所述的阈值Lc的确定方法为:
其中,μ和σ分别为能量图Ec的均值和方差;Pfa为恒虚警概率,取值范围是10%~20%。
有益效果:
1、本发明根据图像像素的叠掩区与非叠掩区二值分割结果,计算各方位视角图像的非叠掩区像素Inon_overlap的多视角可见度LI,通过可见度指标自动实现图像像素的精细分类;此外本发明区分了不同视角观测结果对图像像素的贡献,且仅利用有贡献的角度进行融合,避免了简单融合方法不能获得最优融合结果的缺点;同时本发明原理简单,融合质量佳,且可见度指标计算简单,分类精度高,融合算法工程实现简单;
2、本发明适用于对山地、城市等高地形起伏的目标区域实施合成孔径成像探测,可有效消除单一视角SAR成像探测结果中叠掩、阴影等现象,克服其对目标检测与识别的不利影响,确保复杂场景下对敌军事目标的成像侦察能力不受影响。经本发明方法融合处理的图像,大大提升了图像的可读性和可懂度,非常有利于图像应用时的目标判读和理解。
3、本发明可批量完成各像素的融合处理,运算效率较高,算法易于并行实现。
附图说明
图1为本发明基于可见度分类的多角度SAR图像融合方法流程图;
图2为本发明单一视角SAR图像;
图3为本发明大尺度降质处理后的粗分辨率的降质图像Ic;
图4为本发明单一视角SAR图像的能量图;
图5为本发明能量图经恒虚警检测的叠掩、非叠掩二值分割结果;
图6为本发明非叠掩区目标的可见度分类图;
图7为本发明可见度指标统计分布图;
图8为本发明基于可见度分类的多角度SAR图像融合成像结果。
具体实施方式
下面结合附图和实施例对本发明的基于非叠掩区目标可见度分类的多角度SAR非相干融合成像做进一步详细描述。
实施例一:单一平台SAR多视角融合成像探测。携带SAR载荷的单一平台沿设计的曲线飞行路径,环绕感兴趣的目标区域实施成像探测。由于目标散射特性以及成像空间几何随观测视角的剧烈变化,导致全孔径相干累积处理不可能实现最优的成像探测效果,必须考虑多视角SAR融合成像处理。本发明的基于可见度分类的多角度SAR图像融合方法可用于解决此时的融合成像处理问题。如图1所示,为本发明基于可见度分类的多角度SAR图像融合方法流程图,具体步骤为:
步骤1:确定成像平面,并采用时域快速后向投影算法在与成像平面相对应的成像网格上完成N个子孔径SAR数据的相干累积处理,得到N幅子孔径相干合成图像,这些图像分别对应不同的方位中心视角,其中N至少为3;
步骤2:采用经典的Lee滤波算法对获得的各子孔径相干合成图像实施去相干斑滤波处理;具体实施过程是:
相干斑是一种不相关的乘性噪声,可以表示为:
v=un
其中,v表示SAR图像观测到的灰度值,即步骤1得到一系列子孔径相干合成图像,u为真实目标的雷达反射系数,n为相干斑噪声,通常假设相干斑在统计上独立于信号。假定先验均值和方差可以通过计算匀质区域内的均值和方差得到,将上式用一阶泰勒级数展开,可以得到线性模型:
其中为SAR图像观测到的灰度值的均值,为真实目标的雷达反射系数的均值,为相干斑噪声的均值;
以最小均方差为基础,在固定窗口中计算均值和方差,可以得到滤波公式:
其中为经相干斑滤波处理后的子孔径相干合成图像,即步骤3的输入数据;
k为加权系数。是固定窗口内的像素平均值,为所取窗口内像素方差,为所取窗口内噪声方差。
此外,Lee滤波时,把图像分为三类区域:第一类是均匀区域,其中的相干斑噪声可以简单地用均值滤波平滑掉;第二类是不均匀区域,在去除噪声的同时应该尽量保持其中的结构和纹理信息;第三类是包含分离点目标的区域,滤波器应尽可能地保留原始的观察值。因此亦可采用另一种方法,即如下增强的Lee滤波算子对获得的各子孔径相干合成图像实施去相干斑滤波处理:
其中,NL是图像视数,对于单视复图像,NL=1。统计滤波方法利用局部统计信息对去除相干斑噪声后的数据进行估计,能够自适应的平滑掉均匀区域内的相干斑噪声,同时一定程度上保留图像的纹理信息和边界信息。
步骤3:对各幅相干斑滤波处理后的子孔径相干合成图像的结果实施大尺度降质处理,得到粗分辨率的降质图像Ic。具体实施过程是:
对于二维图像,其尺度降质的二维高斯核函数G(x,y,σ2)定义为:
图像的尺度空间表征可由原始图像与二维高斯核函数的卷积得到,即
式中(x,y)代表图像上像素的位置;为经相干斑滤波处理后的子孔径相干合成图像,即步骤3的输入数据;Ic为降质图像;σ2是尺度空间的空间尺度因子。σ2越小,图像平滑的范围越小,降质的效果越差;σ2越大,图像平滑的范围越大,降质的效果越好;通过选择空间尺度因子σ2,可以得到不同尺度下降质的序列SAR图像,空间尺度因子σ2越大,降质的图像越粗糙。
步骤4:基于降质图像Ic,采用滑动窗口法,以各像素为中心计算窗口内图像能量,得到能量图Ec。其实现过程是:
假设滑动窗口W的尺寸为nw×nw,其取值为经典的窗函数,可以是矩形窗、汉宁窗、海明窗、布拉克曼窗等升余弦窗口,其表达式为:
n=[0 1 2......nw]T
式中α为0~1之间的窗口参数。
滑动窗口W在降质图像Ic上逐像素移动,以每个像素为中心计算窗口内的能量值,并将该能量值作为能量图对应的窗口中心像素点的像素值,直至滑动窗口W的中心遍历所有像素,于是可以得到能量图Ec:
上式中,为了方便计算,取nw为奇数。
步骤5:根据能量图Ec和SAR图像的统计分布特性,并假设图像各像素实施恒虚警概率为Pfa,确定各幅子孔径相干合成图像的分割阈值Lc并完成二值分割;其中,如果任一像素点的像素值大于Lc,则该像素属于叠掩区,否则属于非叠掩区;如图5所示,为本发明能量图经恒虚警检测的叠掩、非叠掩二值分割结果;其实现过程是:
能量图Ec所得实际上是SAR降质图像的幅度,根据SAR图像统计特性理论,幅度图像往往满足瑞利分布;
由于恒虚警概率检测(CFAR)技术本质上是一种根据杂波统计特性在给定的虚警概率条件下自适应获取检测阈值的技术。关键是确定杂波分布的模型和制定的虚警概率,然后计算出分割阈值。假设p(x)为雷达杂波分布模型的概率密度函数,其概率分布函数为F(x)。显然,F(x)在[0,+∞)上为递增函数,假设指定的虚警概率为Pfa,则通过计算下述方程可以得到分割阈值Lc:
其中,μ和σ为能量图Ec的均值和方差;Pfa为恒虚警概率,取值范围是10%~20%;可见,通过计算能量图像Ec的均值和标准差,根据上式即可确定恒虚警概率Pfa条件下的分割阈值,对图像像素实施二值化分割:若I>Lc,则像素属于叠掩区,否则,则属于非叠掩区。
步骤6:根据叠掩区与非叠掩区二值分割的结果,计算各方位子孔径图像的非叠掩区像素Inon_overlap的多视角可见度LI,并依据可见度指标LI实施非叠掩区像素的精细分类。如图6所示,为本发明非叠掩区目标的可见度分类图;图7为本发明可见度指标统计分布图;具体过程是:
假设可见度LI初始值为0,总共有N个不同的方位子孔径图像,定义处于非叠掩区的像素为可见,若成像平面某一点处于某方位子孔径图像的非叠掩区,则其可见度LI增加1。其中,可见度取值为区间[0,N]之间的整数。统计成像平面上各像素的可见度LI,按其指标的大小进行可见度精细分类,对于具有N个不同的方位子孔径图像的情况,成像平面上各像素总共可以划分为N+1个可见度等级;
步骤7:根据非叠掩区像素可见度分类结果,对可见度等级LI相同的各类像素实施线性均值融合或最大对比度融合或最小图像熵融合,从而实现多角度SAR图像的融合处理,得到最终输出的融合图像Iopt,如图8所示,为本发明基于可见度分类的多角度SAR图像融合成像结果,其中:
线性均值融合具体为:分别计算各个可见度等级中包含的所有像素点像素值的平均值,取各等级中的平均值作为自身可见度等级包含像素点的像素值。
最大对比度融合具体为:分别计算各个可见度等级中包含的所有像素点的对比度,取各等级中对比度最大的像素点对应的像素值作为自身可见度等级包含像素点的像素值。
最小图像熵融合具体为:分别计算各个可见度等级中包含的所有像素点的图像熵,取各等级中图像熵最小的像素点对应的像素值作为自身可见度等级包含像素点的像素值。
具体实现如下:
其中,∪I表示各个可见度等级LI的所有像素;mean(∪I)表示取各个可见度等级LI中包含的所有像素点像素值的平均值;max(Icontrast),∪I表示取各等级中对比度最大的像素点对应的像素值作为自身可见度等级包含像素点的像素值;min(Ientropy),∪I表示取各等级中图像熵最小的像素点对应的像素值作为自身可见度等级包含像素点的像素值。
其中,图像均值图像对比度Icontrast以及图像熵Ientropy的具体计算方法如下:
I∈∪I
Imax=max[I(i,j)],i=1,2,3,......nw;j=1,2,3,......nw
Imin=min[I(i,j)],i=1,2,3,......nw;j=1,2,3,......nw
为图像的平均值,Imax和Imin分别为各个可见度等级中包含的所有像素值的最大值和最小值,nw为矩形滑动窗口内总的像素数;(i,j)为像素点的位置;P(n)为各图像像素取值出现的概率,L=nw 2为矩形滑动窗口内的总像素数。
实施例二:多平台SAR同时多视角融合成像探测。携带SAR载荷的多个探测平台沿设计的飞行路径,同时对感兴趣的目标区域实施协同式成像探测。由于各成像平台与目标之间空间几何关系完全不同,导致各平台成像探测时对应的目标散射特性各不相同,而且图像的叠掩、阴影等图像特征和现象也不尽相同,难以实现不同平台不同视角成像结果的相干累积,难以得到最优的成像探测效果。此时,必须考虑多视角SAR融合成像处理,实现不同成像平台、不同视角的成像探测结果的融合,提升对感兴趣的目标区域的高质量成像侦察。本发明的基于可见度分类的多角度SAR图像融合方法可用于解决此时的多平台、多角度融合成像处理问题。如图1所示,为本发明基于可见度分类的多角度SAR图像融合方法流程图,具体步骤为:
步骤1:确定成像平面,并采用时域快速后向投影算法在与成像平面相对应的成像网格上完成N个单一视角SAR数据的相干累积处理,得到N幅单一视角SAR图像,这些图像分别对应不同的方位中心视角,其中N至少为3;如图2所示,为本发明单一视角SAR图像;
步骤2:采用经典的Lee滤波算法对获得的各单一视角SAR图像实施去相干斑滤波处理;具体实施过程是:
相干斑是一种不相关的乘性噪声,可以表示为:
v=un
其中,v表示SAR图像观测到的灰度值,即步骤1得到一系列单一视角SAR图像,u为真实目标的雷达反射系数,n为相干斑噪声,通常假设相干斑在统计上独立于信号。假定先验均值和方差可以通过计算匀质区域内的均值和方差得到,将上式用一阶泰勒级数展开,可以得到线性模型:
其中为SAR图像观测到的灰度值的均值,为真实目标的雷达反射系数的均值,为相干斑噪声的均值;
以最小均方差为基础,在固定窗口中计算均值和方差,可以得到滤波公式:
其中为经相干斑滤波处理后的单一视角SAR图像,即步骤3的输入数据;
k为加权系数。是固定窗口内的像素平均值,为所取窗口内像素方差,为所取窗口内噪声方差。
此外,Lee滤波时,把图像分为三类区域:第一类是均匀区域,其中的相干斑噪声可以简单地用均值滤波平滑掉;第二类是不均匀区域,在去除噪声的同时应该尽量保持其中的结构和纹理信息;第三类是包含分离点目标的区域,滤波器应尽可能地保留原始的观察值。因此亦可采用另一种方法,即如下增强的Lee滤波算子对获得的各单一视角SAR图像实施去相干斑滤波处理:
其中,NL是图像视数,对于单视复图像,NL=1。统计滤波方法利用局部统计信息对去除相干斑噪声后的数据进行估计,能够自适应的平滑掉均匀区域内的相干斑噪声,同时一定程度上保留图像的纹理信息和边界信息。
步骤3:对各幅相干斑滤波处理后的单一视角SAR图像实施大尺度降质处理,得到粗分辨率的降质图像Ic。如图3所示,为本发明大尺度降质处理后的粗分辨率的降质图像Ic;具体实施过程是:
对于二维图像,其尺度降质的二维高斯核函数G(x,y,σ2)定义为:
图像的尺度空间表征可由原始图像与二维高斯核函数的卷积得到,即
式中(x,y)代表图像上像素的位置;为经相干斑滤波处理后的单一视角SAR图像,即步骤3的输入数据;Ic为降质图像;σ2是尺度空间的空间尺度因子。σ2越小,图像平滑的范围越小,降质的效果越差;σ2越大,图像平滑的范围越大,降质的效果越好;通过选择空间尺度因子σ2,可以得到不同尺度下降质的序列SAR图像,空间尺度因子σ2越大,降质的图像越粗糙。
步骤4:基于降质图像Ic,采用滑动窗口法,以各像素为中心计算窗口内图像能量,得到能量图Ec。如图4所示,为本发明单一视角SAR图像的能量图Ec;其实现过程是:
假设滑动窗口W的尺寸为nw×nw,其取值为经典的窗函数,可以是矩形窗、汉宁窗、海明窗、布拉克曼窗等升余弦窗口,其表达式为:
n=[0 1 2......nw]T
式中α为0~1之间的窗口参数。
滑动窗口W在降质图像Ic上逐像素移动,以每个像素为中心计算窗口内的能量值,并将该能量值作为能量图对应的窗口中心像素点的像素值,直至滑动窗口W的中心遍历所有像素,于是可以得到能量图Ec:
上式中,为了方便计算,取nw为奇数。
步骤5:根据能量图Ec和SAR图像的统计分布特性,并假设图像各像素实施恒虚警概率为Pfa,确定各幅单一视角SAR图像的分割阈值Lc并完成二值分割;其中,如果任一像素点的像素值大于Lc,则该像素属于叠掩区,否则属于非叠掩区;如图5所示,为本发明能量图经恒虚警检测的叠掩、非叠掩二值分割结果;其实现过程是:
能量图Ec所得实际上是SAR降质图像的幅度,根据SAR图像统计特性理论,幅度图像往往满足瑞利分布;
由于恒虚警概率检测(CFAR)技术本质上是一种根据杂波统计特性在给定的虚警概率条件下自适应获取检测阈值的技术。关键是确定杂波分布的模型和制定的虚警概率,然后计算出分割阈值。假设p(x)为雷达杂波分布模型的概率密度函数,其概率分布函数为F(x)。显然,F(x)在[0,+∞)上为递增函数,假设指定的虚警概率为Pfa,则通过计算下述方程可以得到分割阈值Lc:
其中,μ和σ为能量图Ec的均值和方差;Pfa为恒虚警概率,取值范围是10%~20%;可见,通过计算能量图像Ec的均值和标准差,根据上式即可确定恒虚警概率Pfa条件下的分割阈值,对图像像素实施二值化分割:若I>Lc,则像素属于叠掩区,否则,则属于非叠掩区。
步骤6:根据叠掩区与非叠掩区二值分割的结果,计算各单一视角SAR图像的非叠掩区像素Inon_overlap的多视角可见度LI,并依据可见度指标LI实施非叠掩区像素的精细分类。如图6所示,为本发明非叠掩区目标的可见度分类图;图7为本发明可见度指标统计分布图;具体过程是:
假设可见度LI初始值为0,总共有N个不同的单一视角SAR图像,定义处于非叠掩区的像素为可见,若成像平面某一点处于某单一视角SAR图像的非叠掩区,则其可见度LI增加1。其中,可见度取值为区间[0,N]之间的整数。统计成像平面上各像素的可见度LI,按其指标的大小进行可见度精细分类,对于具有N个不同的单一视角SAR图像的情况,成像平面上各像素总共可以划分为N+1个可见度等级;
步骤7:根据非叠掩区像素可见度分类结果,对可见度等级LI相同的各类像素实施线性均值融合或最大对比度融合或最小图像熵融合,从而实现多角度SAR图像的融合处理,得到最终输出的融合图像Iopt,如图8所示,为本发明基于可见度分类的多角度SAR图像融合成像结果,其中:
线性均值融合具体为:分别计算各个可见度等级中包含的所有像素点像素值的平均值,取各等级中的平均值作为自身可见度等级包含像素点的像素值。
最大对比度融合具体为:分别计算各个可见度等级中包含的所有像素点的对比度,取各等级中对比度最大的像素点对应的像素值作为自身可见度等级包含像素点的像素值。
最小图像熵融合具体为:分别计算各个可见度等级中包含的所有像素点的图像熵,取各等级中图像熵最小的像素点对应的像素值作为自身可见度等级包含像素点的像素值。
具体实现如下:
其中,∪I表示各个可见度等级LI的所有像素;mean(∪I)表示取各个可见度等级LI中包含的所有像素点像素值的平均值;max(Icontrast),∪I表示取各等级中对比度最大的像素点对应的像素值作为自身可见度等级包含像素点的像素值;min(Ientropy),∪I表示取各等级中图像熵最小的像素点对应的像素值作为自身可见度等级包含像素点的像素值。
其中,图像均值I、图像对比度Icontrast以及图像熵Ientropy的具体计算方法如下:
I∈∪I
Imax=max[I(i,j)],i=1,2,3,......nw;j=1,2,3,......nw
Imin=min[I(i,j)],i=1,2,3,......nw;j=1,2,3,......nw
为图像的平均值,Imax和Imin分别为各个可见度等级中包含的所有像素值的最大值和最小值,nw为矩形滑动窗口内总的像素数;(i,j)为像素点的位置;P(n)为各图像像素取值出现的概率,L=nw 2为矩形滑动窗口内的总像素数。
当然,本发明还可有其他多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (10)
1.一种基于可见度分类的多角度SAR图像融合方法,其特征在于,包括以下步骤:
步骤1:确定成像平面,并采用时域快速后向投影算法在与成像平面相对应的成像网格上完成N个单一视角SAR数据的相干累积处理,得到N幅单一视角SAR图像,其中N至少为3;
步骤2:对各幅单一视角SAR图像实施去相干斑滤波处理;
步骤3:对各幅相干斑滤波处理后的单一视角SAR图像分别进行大尺度降质处理,得到N幅粗分辨率的降质图像Ic;
步骤4:采用滑动窗口法,以降质图像Ic各像素为中心,计算滑动窗口内的图像能量,得到能量图Ec;
步骤5:根据能量图Ec和单一视角SAR图像的统计分布特性,确定各幅单一视角SAR图像的分割阈值Lc并完成二值分割;其中,如果任一像素点的像素值大于Lc,则该像素属于叠掩区,否则属于非叠掩区;
步骤6:计算各单一视角SAR图像的非叠掩区像素Inon_overlap的可见度LI,并依据可见度LI进行非叠掩区像素的精细分类;具体过程是:
假设可见度LI初始值为0,且处于非叠掩区的像素为可见,则如果成像平面任一点处于任一方位的单一视角SAR图像的非叠掩区,则该像素点的可见度LI增加1;其中,可见度取值为0到N;统计成像平面上各像素的可见度LI,按可见度LI的大小将具有N个不同的方位视角的成像平面上各像素划分为N+1个可见度等级;
步骤7:对可见度等级相同的各个像素进行融合,得到最终输出的融合图像Iopt。
2.如权利要求1所述的一种基于可见度分类的多角度SAR图像融合方法,其特征在于,步骤7所述的对可见度等级相同的各个像素进行融合的方法为线性均值融合,具体为:分别计算各个可见度等级中包含的所有像素点像素值的平均值,取各等级中的平均值作为融合图像Iopt相同可见度等级包含像素点的像素值。
3.如权利要求2所述的一种基于可见度分类的多角度SAR图像融合方法,其特征在于,步骤7所述的平均值采用滑动窗口法计算:
其中,为图像的平均值,I(i,j)为像素点的位置,nw为矩形滑动窗口边长包含的像素数,i=1,2,3,......nw;j=1,2,3,......nw。
4.如权利要求1所述的一种基于可见度分类的多角度SAR图像融合方法,其特征在于,步骤7所述的对可见度等级相同的各个像素进行融合的方法为最大对比度融合,具体为:分别计算各个可见度等级中包含的所有像素点的对比度,取各等级中对比度最大的像素点对应的像素值作为融合图像Iopt相同可见度等级包含像素点的像素值。
5.如权利要求4所述的一种基于可见度分类的多角度SAR图像融合方法,其特征在于,步骤7所述的对比度采用滑动窗口法计算:
Imax=max[I(i,j)],i=1,2,3,......nw;j=1,2,3,......nw
Imin=min[I(i,j)],i=1,2,3,......nw;j=1,2,3,......nw
为滑动窗口内的图像I的局部平均值,Imax和Imin分别为各个可见度等级中包含的所有像素值的最大值和最小值,nw为矩形滑动窗口边长包含的像素数;I(i,j)为像素点的位置。
6.如权利要求1所述的一种基于可见度分类的多角度SAR图像融合方法,其特征在于,步骤7所述的对可见度等级相同的各个像素进行融合的方法为最小图像熵融合,具体为:分别计算各个可见度等级中包含的所有像素点的图像熵,取各等级中图像熵最小的像素点对应的像素值作为融合图像Iopt相同可见度等级包含像素点的像素值。
7.如权利要求6所述的一种基于可见度分类的多角度SAR图像融合方法,其特征在于,步骤7所述的图像熵采用滑动窗口法计算:
P(n)为各图像像素取值出现的概率,nw为矩形滑动窗口边长包含的像素数,L=nw 2为矩形滑动窗口内的总像素数,I(i,j)为像素点的位置,j=1,2,3,......nw。
8.如权利要求1所述的一种基于可见度分类的多角度SAR图像融合方法,其特征在于,步骤4所述的采用滑动窗口法,以各像素为中心计算窗口内图像能量,得到能量图Ec具体方法为:
假设滑动窗口W的尺寸为nw×nw,其表达式为:
n=[0 1 2 ...... nw]T
式中α为0~1之间的窗口参数;
滑动窗口W在降质图像Ic上逐像素移动,以每个像素为中心计算窗口内的能量值,并将该能量值作为能量图对应的窗口中心像素点的像素值,直至滑动窗口W的中心遍历所有像素,得到能量图Ec:
其中,nw为奇数。
9.如权利要求8所述的一种基于可见度分类的多角度SAR图像融合方法,其特征在于,所述窗口为矩形窗、汉宁窗、海明窗或布拉克曼窗。
10.如权利要求1所述的一种基于可见度分类的多角度SAR图像融合方法,其特征在于,步骤5所述的阈值Lc的确定方法为:
其中,μ和σ分别为能量图Ec的均值和方差;Pfa为恒虚警概率,取值范围是10%~20%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710042610.1A CN106897985B (zh) | 2017-01-20 | 2017-01-20 | 一种基于可见度分类的多角度sar图像融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710042610.1A CN106897985B (zh) | 2017-01-20 | 2017-01-20 | 一种基于可见度分类的多角度sar图像融合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106897985A CN106897985A (zh) | 2017-06-27 |
CN106897985B true CN106897985B (zh) | 2019-10-29 |
Family
ID=59198162
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710042610.1A Active CN106897985B (zh) | 2017-01-20 | 2017-01-20 | 一种基于可见度分类的多角度sar图像融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106897985B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107590790B (zh) * | 2017-09-21 | 2021-04-13 | 长沙全度影像科技有限公司 | 一种基于对称边缘填充的简单透镜边缘区域去模糊方法 |
CN107861121A (zh) * | 2017-09-28 | 2018-03-30 | 北京华航无线电测量研究所 | 一种基于高斯模糊和cfar检测的危险物品检测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102044072A (zh) * | 2010-11-29 | 2011-05-04 | 北京航空航天大学 | 基于统计模型的sar图像融合处理方法 |
CN103871040A (zh) * | 2014-03-12 | 2014-06-18 | 北京理工大学 | 基于多角度导航卫星双基地合成孔径雷达图像融合方法 |
CN105809649A (zh) * | 2016-03-03 | 2016-07-27 | 西安电子科技大学 | 基于变分多尺度分解的sar图像与可见光图像融合方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9275445B2 (en) * | 2013-08-26 | 2016-03-01 | Disney Enterprises, Inc. | High dynamic range and tone mapping imaging techniques |
-
2017
- 2017-01-20 CN CN201710042610.1A patent/CN106897985B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102044072A (zh) * | 2010-11-29 | 2011-05-04 | 北京航空航天大学 | 基于统计模型的sar图像融合处理方法 |
CN103871040A (zh) * | 2014-03-12 | 2014-06-18 | 北京理工大学 | 基于多角度导航卫星双基地合成孔径雷达图像融合方法 |
CN105809649A (zh) * | 2016-03-03 | 2016-07-27 | 西安电子科技大学 | 基于变分多尺度分解的sar图像与可见光图像融合方法 |
Non-Patent Citations (2)
Title |
---|
Multi-angle SAR non-coherent image fusion algorithm based on HIS statistic characteristics;Da Ran et al.;《Eighth International Conference on Digtial Image Processing》;20160831;第1-6页 * |
多角度合成孔径雷达成像技术研究进展;冉达 等;《装备学院学报》;20160831;第27卷(第4期);第86-92页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106897985A (zh) | 2017-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106910177B (zh) | 一种局域图像指标最优化的多角度sar图像融合方法 | |
CN108387896B (zh) | 一种基于探地雷达回波数据的自动汇聚成像方法 | |
US8125370B1 (en) | Polarimetric synthetic aperture radar signature detector | |
US8193967B2 (en) | Method and system for forming very low noise imagery using pixel classification | |
Zeng et al. | Multiangle BSAR imaging based on BeiDou-2 navigation satellite system: Experiments and preliminary results | |
JP5932643B2 (ja) | 同一領域上で取得されたsar画像から得られるインターフェログラムのフィルタリング処理方法 | |
Chen et al. | A 3D reconstruction strategy of vehicle outline based on single-pass single-polarization CSAR data | |
CN103605119B (zh) | 一种星载合成孔径雷达条带模式方位模糊抑制方法 | |
CN113960595A (zh) | 一种地表形变监测方法及系统 | |
CN113570632B (zh) | 基于高时相星载sar序贯图像的运动小目标检测方法 | |
Rossi et al. | High-resolution InSAR building layovers detection and exploitation | |
CN116047519A (zh) | 一种基于合成孔径雷达干涉测量技术的选点方法 | |
Jiang et al. | Delineation of built-up land change from SAR stack by analysing the coefficient of variation | |
Dungan et al. | Wide-area wide-angle SAR focusing | |
CN112415515A (zh) | 一种机载圆迹sar对不同高度目标分离的方法 | |
CN113064162A (zh) | 应用于机场跑道异物检测雷达系统的检测方法、装置 | |
Newey et al. | Detection and imaging of moving targets with limit SAR data | |
Girón et al. | Nonparametric edge detection in speckled imagery | |
CN114998365A (zh) | 一种基于极化干涉sar的地物分类方法 | |
CN106897985B (zh) | 一种基于可见度分类的多角度sar图像融合方法 | |
CN106910178B (zh) | 一种基于色调统计特性分类的多角度sar图像融合方法 | |
Li et al. | Building detection in sar images based on bi-dimensional empirical mode decomposition algorithm | |
Zhao et al. | SAR deception jamming identification via differential feature enhancement | |
Cetin et al. | Evaluation of a regularized SAR imaging technique based on recognition-oriented features | |
Garren | SAR focus theory of complicated range migration signatures due to moving targets |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: 101416 Beijing City Bayi Yanqi Town Huairou District Road No. Patentee after: China People's Liberation Army Strategic Support Unit Space Engineering University Address before: 3380 box 101416, Huairou District, Beijing, China Patentee before: PLA 'S Equipment College |