CN107918126B - 基于多特征自动分割的多通道近海岸模糊杂波抑制方法 - Google Patents
基于多特征自动分割的多通道近海岸模糊杂波抑制方法 Download PDFInfo
- Publication number
- CN107918126B CN107918126B CN201711144523.3A CN201711144523A CN107918126B CN 107918126 B CN107918126 B CN 107918126B CN 201711144523 A CN201711144523 A CN 201711144523A CN 107918126 B CN107918126 B CN 107918126B
- Authority
- CN
- China
- Prior art keywords
- clutter
- interference
- channel
- fuzzy
- view
- 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
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000001629 suppression Effects 0.000 title claims abstract description 24
- 230000011218 segmentation Effects 0.000 title claims abstract description 9
- 238000001514 detection method Methods 0.000 claims abstract description 19
- 239000011159 matrix material Substances 0.000 claims abstract description 19
- 238000012549 training Methods 0.000 claims abstract description 10
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 7
- 238000009499 grossing Methods 0.000 claims abstract description 5
- 230000008569 process Effects 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 6
- 230000004927 fusion Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 238000012937 correction Methods 0.000 claims description 2
- 238000012804 iterative process Methods 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 2
- 230000002401 inhibitory effect Effects 0.000 abstract 2
- 238000004088 simulation Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010587 phase diagram Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9029—SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/354—Extracting wanted echo-signals
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提出一种基于多特征自动分割的近海岸模糊杂波位置确定与抑制方法,主要解决现有技术在方位多通道SAR体制下对近海岸模糊杂波位置确定与抑制准确性差的问题。其方案是:首先建立多视干涉有限混合杂波模型;然后结合杂波模型利用期望最大化算法实现SAR图像自动分割,并利用马尔科夫随机场平滑技术剔除孤立错分点;最后根据获得的模糊杂波位置和训练样本构造杂波协方差矩阵实现方位模糊杂波抑制。本发明在不依赖于精确的系统参数和方位模糊源位置信息的前提下,能够获得精确的方位模糊位置信息和良好的杂波抑制性能,可应用于近海岸区域慢速动目标的检测。
Description
技术领域
本发明属于雷达技术领域,涉及近海岸模糊杂波的抑制,可应用于近海岸区域慢速动目标的检测。
背景技术
对于近海岸区域而言,临近强陆地杂波的方位模糊分量或者“鬼影”是明显可见的,并且成为运动目标或舰船检测的主要干扰源。因此,抑制方位模糊成为了近海岸区域运动目标或舰船检测的首要任务。根据通道自由度划分,方位模糊抑制方法可以分为单通道和多通道方法。无论是单通道方法还是多通道方法,在后处理技术当中大部分的方位模糊抑制方法的前提在于方位模糊位置的确定。其中:
单通道方法方位模糊位置的确定严重依赖于参数的精确程度,而且在某些特殊情况下还可能会失效。多通道图像域模糊杂波抑制方法要求精确获得模糊杂波训练样本,否则杂波抑制性能受限。束宇翔,何嘉懿,廖桂生等人在多通道SAR近岸水面区域模糊杂波抑制方法(电子与信息学报,2014,36(5):1030-1035.)一文中提出利用SAR图像的干涉相位特征和幅度特征联合进行模糊杂波位置确定,根据幅度和相位直方图选取对应的分割门限进行处理,该方法由于要求干涉相位和幅度直方图具有明显可区分的多峰特性,且对应的门限参数需要人工选择,因而适用性受到了制约。
发明内容
本发明的目的在于,针对上述现有技术的不足,在方位多通道SAR体制下,提出一种基于多特征自动分割的近海岸模糊杂波抑制方法,以避免人工选择过程,利用杂波的多特征对杂波进行自动分类,扩大其适用性。
本发明的技术思路是:利用双通道SAR图像的多视干涉幅度和相位特征,采用期望最大化EM算法实现SAR图像自动分割;根据马尔科夫随机场MRF平滑技术剔除孤立错分点,并从分割结果中获得模糊杂波位置;利用获得的模糊杂波样本构造模糊杂波协方差矩阵进行模糊杂波抑制。其实现步骤包括如下:
本发明的技术方案包括如下步骤:
(1)获取沿航迹多视SAR图像的干涉幅度η、干涉相位ψ和干涉复图像;
(2)获得沿航迹多视SAR图像的归一化干涉幅度和干涉相位的联合概率密度函数:
其中,P表示杂波类数,αp表示第p类杂波所占的比例,表示参数为Θp={np,ρp,θp,σp,1,σp,2}的第p类杂波的归一化干涉幅度和干涉相位联合概率密度函数,n表示多视数,ρ表示两个通道输出间的实相干系数,θ表示杂波多视干涉相位均值,和分别表示第1通道和第2通道的功率,Ψ为干涉相位随机变量,ε为归一化幅度随机变量。
(3)对SAR干涉复图像进行降采样处理,得到降采样后的干涉复图像;
(4)利用最大期望算法EM对降采样后的干涉复图像进行杂波分类,得到不同参数特征的杂波集合;
(5)利用杂波有限混合模型和步骤(5)获得的不同参数特征的杂波集合构造皮尔逊卡方拟合检测统计量χ2;
(6)对统计量χ2进行检测:若检测统计量χ2满足则结束分类过程,执行步骤(7),否则,继续对分类结果进行分类,其中表示卡方分布χ2(Ac-1)的上α分位点,1-α表示置性水平,Ac表示样本总体所有可能取值组成的全体所划分成的互不相交子集的个数;
(7)用MRF方法对分类结果进行平滑,得到平滑后的杂波分类结果;
(8)根据SAR图像分类结果可以进一步获得模糊杂波的位置,记模糊杂波的位置集合为Ωamb。由此可以获取模糊杂波多通道的训练样本数据集合 其中x(i,j)=[x1(i,j),...,xM(i,j)]T表示SAR图像像素点(i,j)处的通道数为M的方位多通道数据,xm(i,j)表示第m个方位通道SAR图像像素点(i,j)处的复数值,其中m=1,...,M,上标T表示矩阵转置操作;
(9)利用平滑之后的杂波分类结果构造杂波抑制矢量:
本发明与现有技术相比具有以下优点:
第一,本发明利用EM等算法自动对杂波进行分类,可以将近海岸处各类杂波分类完好,相比于现有技术,不会将陆地与海岸杂波错分为模糊杂波。
第二,本发明在不依赖于精确的雷达系统参数和方位模糊源位置信息的前提下,相比于现有技术,可以获得精确的方位模糊位置和良好的杂波抑制性能。
附图说明
图1为本发明使用的场景图;
图2为本发明的实现流程图;
图3为用本发明仿真过程中的星载SAR数据幅度图;
图4为用本发明仿真过程中的星载SAR数据干涉相位图;
图5为用本发明对SAR图像杂波自动分类结果图;
图6为用本发明对模糊杂波位置的检测结果图;
图7为用本发明对方位模糊杂波的抑制结果图。
具体实施方式
下面结合附图对本发明做更清楚、更完整的描述。
参照图1,本发明使用的场景有M个通道,第m个通道的位置坐标为(-dm,0),其中dm为第m个通道与第一个通道之间的距离,m=2,...,M,杂波点的位置为(X0,R0)。
参照图2,本发明的实施步骤如下:
步骤1,获取沿航迹多视SAR图像的干涉幅度η、干涉相位ψ和干涉复图像。
(1a)获取干涉复图像:
先获得多视SAR干涉复图像中像素点(i,j)处的复数值In(i,j)为:
其中,z1(k)和z2(k)分别表示在以像素点(i,j)为中心一个空域多视窗内,单视SAR复图像通道1和通道2的第k个独立像素点复数值,n表示多视数,即独立样本个数,上标“*”表示复共轭操作;
再利用各像素点的复数值In(i,j)求得干涉复图像;
(1b)获取干涉幅度η和干涉相位ψ:
(1b1)获得n视样本协方差矩阵:
其中,上标H表示矩阵复共轭转置操作,T表示矩阵转置操作,*表示共轭操作;(1b2)将(1b1)中的次对角元素作为多视干涉复数值:令In(i,j)=ηexp(jψ),则得到多视干涉幅度η和干涉相位ψ如下:
ψ=angle{In(i,j)},
其中angle{·}表示取辐角操作。
步骤2,获得沿航迹多视SAR图像的归一化干涉幅度和干涉相位的联合概率密度函数gε,Ψ(η,ψ)。
(2a)获取样本协方差矩阵期望R为:
(2b)根据(2a)所得结果,获得随机变量干涉相位Ψ=ψ和归一化幅度ε=η/(σ1σ2)的联合概率密度,为:
其中,Γ(·)表示伽马函数,Kn-1(·)表示n-1阶修正的贝塞尔函数,参数集合Θ={n,ρ,θ,σ1,σ2},该集合中的参数通过文献C.H.Gierull.Statistical analysis ofmultilook SAR interferograms for CFAR detection of ground moving targets[J]的方法估计得出;
(2c)利用(2b)所得结果,获得沿航迹多视SAR图像的,T类杂波的归一化干涉幅度和归一化干涉相位的联合概率密度函数:
步骤3,对SAR干涉复图像进行降采样处理,得到降采样后的干涉复图像。
步骤4,利用最大期望算法EM对降采样后的干涉复图像进行杂波分类,得到不同参数特征的杂波集合。
其中P为分类总类数,上标m为第m次迭代;
(4c)对第m次迭代分类结果分别计算各类的参数:
其中,表示第p类杂波第m+1次迭代时所占的比例,为第m+1次迭代时的参数集合,为m+1次迭代时的第p类杂波的多视数,为m+1次迭代时通道1和通道2之间的实相干系数,为m+1次迭代时的干涉相位均值,和为m+1次迭代时通道1和通道2的功率,N表示图像总的像素点数,#{·}表示取元素个数操作,当迭代次数大于10次或者前后两次迭代的杂波所占的比例相对变化小于10-3时,迭代过程终止。
步骤5,根据文献E.L.Lehmann and J.P.Romano.Testing statisticalhypotheses 3rd edition(Springer texts in statistics)[M]中的方法利用有限杂波混合模型和步骤4获得的不同参数特征的杂波集合,构造皮尔逊卡方拟合检测统计量χ2:
其中,为干涉相位取值范围所划分为不相交子集的个数,为归一化干涉幅度取值范围所划分为不相交子集的个数,fPHk,Mgk为所有检验杂波样本中第PHk,Mgk类杂波样本的个数,ePHk,Mgk为第PHk,Mgk类杂波样本发生的概率。
步骤6,对统计量χ2进行检测。
(6c)由卡方分布得出样本个数门限κ,将最坏类样本个数κp与门限κ进行比较:
若最坏类pworse的样本个数κp小于等于门限κ,则将该最坏类与剩余类中参数差异最小的类融合,并根据融合结果更新参数集合Θp,返回到步骤4;
若最坏类pworse的样本个数κp大于门限κ,则进行随机对半分裂,并根据分裂或者融合结果更新参数集合Θp,返回到步骤4。
步骤7,用马尔科夫随机场MRF方法对分类结果进行平滑,得到平滑后的杂波分类结果。
步骤9,利用平滑之后的杂波分类结果构造杂波抑制矢量:
表示利用模糊杂波的训练数据集合估计的模糊杂波样本协方差矩阵,其中,L表示模糊杂波的训练数据集合的元素个数,上标H表示矩阵复共轭转置操作,(·)-1为矩阵求逆操作,表示径向速度为vr的目标方位多通道导向矢量,dm是第m个方位通道与参考通道间的等效相位中心间距,m=2,...,M,va为雷达平台的恒定速度,λ为波长;
步骤10,将杂波抑制矢量wamb加权于SAR图像上像素点(i,j)处的方位多通道数据x(i,j),即wamb乘以x(i,j)后可以获得模糊杂波抑制后的SAR图像像素点(i,j)处的复数值y(i,j):
本发明的效果可以通过以下仿真进一步说明:
一.仿真条件:
仿真场景的散射系数根据文献J.Kim,M.Younis,P.Prats-Iraola,et al.,"Firstspaceborne demonstration of digital beamforming for azimuth ambiguitysuppression,"中的SAR图像得到,中间区域为陆地,陆地两边为平静水面。
仿真参数:卫星等效雷达速度7500m/s,波长0.0313m,天线方位孔径10m,均匀划分为2个通道,全孔径发射子孔径同时接收,相邻等效自发自收相位中心间距d2=2.5m,PRF=1200Hz,主杂波的杂噪比CNR=25dB。
二.仿真内容
仿真1,用本发明对星载SAR场景接收的数据成像,结果如图3和4所示,其中:
图3为星载SAR图像数据的单视幅度图。
图4为星载SAR图像数据的多视干涉相位图。
从图3和图4可以看出,仿真星载SAR图像中包含了陆地主杂波、左右一次陆地模糊杂波以及平静的海面。同时根据图4所示的3×3的空域多视窗多视干涉相位图可以发现,左右一次模糊杂波两通道间的干涉相位均值的理论计算值θamb=±2π·PRF·d2/va≈±2.5,其与实际值相符。
仿真2,对仿真1所得的结果进行杂波分类,结果如图5.
由图5可以看出,本发明的方法可以将海杂波、陆地主杂波、左右一次陆地模糊杂波三类杂波被很好地分割开,说明了本发明方法可以有效进行SAR图像杂波分类。
仿真3,对仿真1所得的结果方位模糊杂波位置检测,结果如图6,其中:
图6(a)为左一次模糊检测结果,图6(b)为右一次模糊检测结果;
由图6可以进一步看出,本发明的方法可以将左右一次陆地模糊杂波同其它完好的分类开。
仿真4,对图5中分类后的杂波进行杂波抑制,结果如图7;
图7(a)为左一次模糊抑制结果,图7(b)为右一次模糊抑制结果;
从图7可以看出,本发明方法对模糊杂波的抑制的性能良好。
上述仿真实验结果表明:本发明的方法可以有效地将杂波分类开,获得模糊杂波位置,并对杂波可以进行有效的抑制,该实验结果验证了本发明的有效性和正确性。
以上描述仅是本发明的一个具体实例,并未构成对本发明的任何限制,显然对于本领域的专业人员来说,在了解了本发明内容和原理后,都可能在不背离本原理、结构的情况下,进行形式和细节上的各种修改和改变,但是这些基于本发明的思想的修正和改变仍在本发明的权利要求保护范围之内。
Claims (5)
1.基于多特征自动分割的多通道近海岸模糊杂波抑制方法,其特征在于,包括:
(1)获取沿航迹多视SAR图像的干涉幅度η、干涉相位ψ和干涉复图像;
(2)获得沿航迹多视SAR图像的归一化干涉幅度和干涉相位的联合概率密度函数:
其中,P表示杂波类数,αp表示第p类杂波所占的比例,表示参数为Θp={np,ρp,θp,σp,1,σp,2}的第p类杂波的归一化干涉幅度和干涉相位联合概率密度函数,n表示多视数,ρ表示两个通道输出间的实相干系数,θ表示杂波多视干涉相位均值,σp,1和σp,2分别表示第p类杂波在第1通道和第2通道的幅度,Ψ为干涉相位随机变量,ε为归一化幅度随机变量;
(3)对SAR干涉复图像进行降采样处理,得到降采样后的干涉复图像;
(4)利用最大期望算法EM对降采样后的干涉复图像进行杂波分类,得到不同参数特征的杂波集合;
(5)利用杂波有限混合模型和步骤(4)获得的不同参数特征的杂波集合构造皮尔逊卡方拟合检测统计量χ2;
(6)对统计量χ2进行检测:若检测统计量χ2满足则结束分类过程,执行步骤(7),否则,继续对分类结果进行分类,其中表示卡方分布χ2(Ac-1)的上α分位点,1-α表示置性水平,Ac表示样本总体所有可能取值组成的全体所划分成的互不相交子集的个数;
(7)用马尔科夫随机场MRF方法对分类结果进行平滑,得到平滑后的杂波分类结果;
(8)根据杂波分类结果获得模糊杂波的位置,记模糊杂波的位置集合为Ωamb。由此可以获取模糊杂波多通道的训练样本数据集合其中x(i,j)=[x1(i,j),…,xM(i,j)]T表示SAR图像像素点(i,j)处的通道数为M的方位多通道数据,xm(i,j)表示第m个方位通道SAR图像像素点(i,j)处的复数值,其中m=1,...,M,上标T表示矩阵转置操作;
(9)利用平滑之后的杂波分类结果构造杂波抑制矢量:
表示利用模糊杂波的训练数据集合估计的模糊杂波样本协方差矩阵,其中,L表示模糊杂波的训练数据集合的元素个数,上标H表示矩阵复共轭转置操作,(·)-1为矩阵求逆操作,表示径向速度为vr的目标方位多通道导向矢量,dm是第m个方位通道与参考通道间的等效相位中心间距,m=2,...,M,va为雷达平台的恒定速度,λ为波长;
(10)将杂波抑制矢量wamb加权于SAR图像上像素点(i,j)处的方位多通道数据x(i,j),即wamb乘以x(i,j)后可以获得模糊杂波抑制后的SAR图像像素点(i,j)处的复数值y(i,j):
2.根据权利要求1中所述的方法,其特征在于步骤(1)中获取沿航迹多视SAR图像的干涉幅度η、干涉相位ψ和干涉复图像,按如下步骤进行:
(1a)获取干涉复图像:
先获得多视SAR干涉复图像中像素点(i,j)处的复数值In(i,j)为:
其中,z1(k)和z2(k)分别表示在以像素点(i,j)为中心一个空域多视窗内,单视SAR复图像通道1和通道2的第k个独立像素点复数值,n表示多视数,即独立样本个数,上标“*”表示复共轭操作;
再利用各像素点的复数值In(i,j)求得干涉复图像;
(1b)获取干涉幅度η和干涉相位ψ:
(1b1)获得n视样本协方差矩阵:
其中,上标H表示矩阵复共轭转置操作,T表示矩阵转置操作,*表示共轭操作;
令In(i,j)=ηexp(jψ),则得到多视干涉幅度η和干涉相位ψ如下:
ψ=angle{In(i,j)},
其中angle{·}表示取辐角操作。
3.根据权利要求1中所述的方法,其特征在于步骤(2)中获得沿航迹多视SAR图像的归一化干涉幅度和干涉相位的联合概率密度函数,按如下步骤进行:
(2a)获取样本协方差矩阵期望R为:
(2b)根据(2a)所得结果,获得随机变量干涉相位Ψ=ψ和归一化幅度ε=η/(σ1σ2)的联合概率密度,为:
其中,Γ(·)表示伽马函数,Kn-1(·)表示n-1阶修正的贝塞尔函数,参数集合Θ={n,ρ,θ,σ1,σ2};
(2c)根据(2b)所得结果,获得沿航迹多视SAR图像的P类杂波的归一化干涉幅度和归一化干涉相位的联合概率密度函数:
4.根据权利要求1中所述的方法,其特征在于步骤(4)中利用最大期望算法EM对降采样后的干涉复图像进行杂波分类,按如下步骤进行:
其中P为分类总类数,上标m为第m次迭代;
(4c)对第m次迭代分类结果分别计算各类的参数:
(6c)由卡方分布可得出样本个数门限κ,将最坏类样本个数κp与门限κ进行比较:
若最坏类pworse的样本个数κp小于等于门限κ,则将该最坏类与剩余类中参数差异最小的类r融合,并根据融合结果更新参数集合Θr,此时杂波类数减1,返回到步骤4;
若最坏类pworse的样本个数κp大于门限κ则最坏类pworse进行随机对半分裂,获得两个新的杂波类t和h,并更新参数集合Θt和Θh,此时杂波类数加1,返回到步骤4。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711144523.3A CN107918126B (zh) | 2017-11-17 | 2017-11-17 | 基于多特征自动分割的多通道近海岸模糊杂波抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711144523.3A CN107918126B (zh) | 2017-11-17 | 2017-11-17 | 基于多特征自动分割的多通道近海岸模糊杂波抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107918126A CN107918126A (zh) | 2018-04-17 |
CN107918126B true CN107918126B (zh) | 2021-05-18 |
Family
ID=61896624
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711144523.3A Active CN107918126B (zh) | 2017-11-17 | 2017-11-17 | 基于多特征自动分割的多通道近海岸模糊杂波抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107918126B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108445466B (zh) * | 2018-05-04 | 2021-09-28 | 西安电子科技大学 | 一种基于下雨腐蚀算法的雷达相邻目标分割方法 |
CN109975804B (zh) * | 2019-03-04 | 2023-04-07 | 广东工业大学 | 一种多平台星座sar融合相参成像方法 |
CN111812598B (zh) * | 2020-07-30 | 2022-03-22 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种基于时域频域多特征的地海杂波分类方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102928840A (zh) * | 2012-10-29 | 2013-02-13 | 中国人民解放军空军装备研究院侦察情报装备研究所 | 多通道sar地面慢动目标检测方法及装置 |
CN104851097A (zh) * | 2015-05-19 | 2015-08-19 | 西安电子科技大学 | 基于目标形状与阴影辅助的多通道sar-gmti方法 |
CN106569212A (zh) * | 2016-11-09 | 2017-04-19 | 西安空间无线电技术研究所 | 一种多通道sar‑gmti距离模糊杂波抑制方法 |
US9791563B1 (en) * | 2014-01-08 | 2017-10-17 | National Technology & Engineering Solutions Of Sandia, Llc | Joint synthetic aperture radar plus ground moving target indicator from single-channel radar using compressive sensing |
-
2017
- 2017-11-17 CN CN201711144523.3A patent/CN107918126B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102928840A (zh) * | 2012-10-29 | 2013-02-13 | 中国人民解放军空军装备研究院侦察情报装备研究所 | 多通道sar地面慢动目标检测方法及装置 |
US9791563B1 (en) * | 2014-01-08 | 2017-10-17 | National Technology & Engineering Solutions Of Sandia, Llc | Joint synthetic aperture radar plus ground moving target indicator from single-channel radar using compressive sensing |
CN104851097A (zh) * | 2015-05-19 | 2015-08-19 | 西安电子科技大学 | 基于目标形状与阴影辅助的多通道sar-gmti方法 |
CN106569212A (zh) * | 2016-11-09 | 2017-04-19 | 西安空间无线电技术研究所 | 一种多通道sar‑gmti距离模糊杂波抑制方法 |
Non-Patent Citations (2)
Title |
---|
SAR Clutter Suppression Using Probability Density Skewness;R. Fogler et al.;《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》;19940430;第3卷(第2期);第622-626页 * |
多通道SAR 近岸水面区域模糊杂波抑制方法;束宇翔 等;《电子与信息学报》;20140531;第36卷(第5期);第1030-1035页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107918126A (zh) | 2018-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103744068B (zh) | 双通道调频连续波sar系统的动目标检测成像方法 | |
Santi et al. | Spatial resolution improvement in GNSS-based SAR using multistatic acquisitions and feature extraction | |
CN110146857B (zh) | 一种颠簸平台sar三维运动误差估计方法 | |
CN108469608B (zh) | 一种运动平台雷达多普勒质心精确估计方法 | |
CN107918126B (zh) | 基于多特征自动分割的多通道近海岸模糊杂波抑制方法 | |
US8816896B2 (en) | On-board INS quadratic correction method using maximum likelihood motion estimation of ground scatterers from radar data | |
CN104730505A (zh) | 一种多通道sar地面目标检测定位方法和装置 | |
US8937849B1 (en) | Auto-focus for circular synthetic aperture sonar | |
CN109557533B (zh) | 一种基于模型的联合跟踪与识别方法 | |
CN114545411B (zh) | 一种基于工程实现的极坐标格式多模高分辨sar成像方法 | |
Ma et al. | Minimum spanning tree co-registration approach for time-series Sentinel-1 TOPS data | |
CN113589287B (zh) | 合成孔径雷达稀疏成像方法、装置、电子设备及存储介质 | |
CN104502895A (zh) | 一种稳健的基于多通道sar特征投影的杂波抑制方法及装置 | |
CN104977585A (zh) | 一种鲁棒的运动声纳目标检测方法 | |
CN109061640B (zh) | 一种用于顺轨干涉sar海流反演的方位模糊抑制方法 | |
WO2018037533A1 (ja) | レーダ装置 | |
Newey et al. | Detection and imaging of moving targets with limit SAR data | |
JP2009236720A (ja) | 移動目標検出装置 | |
CN114089333B (zh) | 一种基于直升机平台的sar振动误差估计及补偿方法 | |
CN108020835B (zh) | 一种星载合成孔径雷达sar地面运动目标指示gmti的强杂波抑制方法及装置 | |
CN105572648A (zh) | 一种合成孔径雷达回波数据距离徙动校正方法和装置 | |
WO2023186386A1 (en) | Detecting and suppressing ambiguities in synthetic aperture radar data and images | |
Wang et al. | High-coherence oriented image formation algorithm based on adaptive elevation ramp fitting for GNSS-based InBSAR systems | |
Wang et al. | Seafloor terrain detection from acoustic images utilizing the fast two-dimensional CMLD-CFAR | |
CN108983192B (zh) | 基于gps辐射源的雷达运动目标参数估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |