CN107229048A - 一种高分宽幅sar动目标速度估计与成像方法 - Google Patents
一种高分宽幅sar动目标速度估计与成像方法 Download PDFInfo
- Publication number
- CN107229048A CN107229048A CN201710415901.0A CN201710415901A CN107229048A CN 107229048 A CN107229048 A CN 107229048A CN 201710415901 A CN201710415901 A CN 201710415901A CN 107229048 A CN107229048 A CN 107229048A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- moving
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 41
- 239000011159 matrix material Substances 0.000 claims description 26
- 238000005070 sampling Methods 0.000 claims description 13
- 238000012546 transfer Methods 0.000 claims description 10
- 230000005540 biological transmission Effects 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 5
- 230000001131 transforming effect Effects 0.000 claims description 5
- 230000009467 reduction Effects 0.000 claims description 4
- 241001269238 Data Species 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 2
- 230000003044 adaptive effect Effects 0.000 abstract description 2
- 238000003672 processing method Methods 0.000 abstract description 2
- 238000000034 method Methods 0.000 description 49
- 238000011045 prefiltration Methods 0.000 description 12
- 238000005516 engineering process Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S13/589—Velocity or trajectory determination systems; Sense-of-movement determination systems measuring the velocity vector
-
- 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
- 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/9004—SAR image acquisition 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/904—SAR modes
-
- 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/904—SAR modes
- G01S13/9041—Squint mode
-
- 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/904—SAR modes
- G01S13/9047—Doppler beam sharpening mode
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动目标速度估计与成像方法,在高分宽幅方位多通道模型下,首先利用空时自适应处理方法对多通道距离压缩后的数据进行杂波抑制;然后采用吕氏分布方法,估计动目标距离向和方位向速度;而后利用动目标距离向速度抑制动目标方位模糊;最后,利用后向投影算法(BP)及动目标距离向和方位向速度实现高分辨成像。
Description
技术领域
本发明属于雷达技术领域,特别涉及一种高分宽幅SAR体制下的动目标检测与成像技术,即GMTI技术。
背景技术
高分宽幅SAR(HRWS,High Resolution Wide Swath,高分宽幅;SAR,SyntheticAperture Radar,合成孔径雷达)是一种可以同时获取高分辨率和宽测绘带的SAR成像技术,可以应用到土地规划、海洋监测、地图绘制等要求大面积高精度的成像中。高分宽幅下的动目标成像可以应用到交通运输中检测车辆及战场侦察信息等,具有重要作用。其中,动目标速度的准确估计是其高质量成像的关键前提,因此构建高性能的动目标速度估计方法及在此基础上的动目标成像是HRWS SAR成像应用的重要基础。
目前已公开发表的HRWS SAR动目标成像相关论文中,最具有代表性的包括Baumgartner总结HRWS条件下利用穷举搜索法来估计动目标速度(见文献1:Baumgartner SV,Krieger G.Simultaneous high-resolution wide-swath SAR imaging and groundmoving target indication:Processing approaches and system concepts[J].IEEEJournal of Selected Topics in Applied Earth Observations and Remote Sensing,2015,8(11):5015-5029.);ShuangXi-Zhang等人提出在方位多通道条件下,基于最小熵的方法来估计动目标速度(见文献2:Zhang S X,Xing M D,Xia X G,et al.A novel movingtarget imaging algorithm for HRWS SAR based on local maximum-likelihoodminimum entropy[J].IEEE Transactions on Geoscience and Remote Sensing,2014,52(9):5333-5348.);XiangyuWang进一步总结了穷举法估计动目标速度的应用(见文献3:Wang X,Wang R,Li N,et al.A velocity estimation method of moving target forSAR high-resoluton wide-swath mode[C]//Geoscience and Remote SensingSymposium(IGARSS),2016IEEE International.IEEE,2016:6819-6822.)。然而,上述动目标速度估计方法中,穷举法需要对每一个假设的动目标速度,完成后续处理看最优检测结果,存在很多冗余的处理步骤,导致运算量大,而且因为假设的速度有限,可能导致估计精度低;而另一类基于最小熵的方法,由于需要用最优化的方法搜索满足最小熵的最优解,算法的运算量也很大,为而且得到的最优解也和算法收敛标准的设置有关,差的收敛标准也会直接导致速度估计精度降低。因此这两类动目标速度估计的方法都存在运算量大和潜在的精度低的问题,会给后续动目标成像处理带来影响,导致成像质量变差。
发明内容
本发明为克服现有HRWS SAR动目标成像中动目标速度估计方法运算量大的不足,而且弥补现有方法一般只能估计动目标距离向速度的缺陷,抑制了动目标方位向速度对成像结果的影响;通过将吕氏分布(LVD)方法应用到HRWS SAR动目标成像中,可以高效估计动目标距离向和方位向速度,利用估计得到的动目标速度信息结合多通道处理方法得到高分辨率动目标成像。
本发明采用的技术方案是:一种高分宽幅SAR动目标速度估计与成像方法,包括:
S1、初始化方位多通道***参数,通过仿真得到多通道回波数据;
所述初始化方位多通道***参数具体为:发射信号脉冲宽度,记为Tr;***脉冲重复频率,记为PRF;信号采样频率,记为Fs;通道数为N,N为偶数;发射信号带宽,记为Br;信号多普勒带宽为Bdop;平台速度为Vs;发射信号波长,记为λ;通道间隔为d;雷达中心频率,记为f0;
根据初始化的方位多通道***参数,通过仿真得到N个通道回波数据;所述每个通道回波数据均为Na×Nr维的矩阵;Na,Nr分别为方位向和距离向采样点数;
S2、将步骤S1得到的各通道的回波数据,分别进行距离压缩,得到各通道压缩后的回波数据;
S3、将步骤S2得到各通道回波数据进行两两校准,得到各个通道校准后的回波数据;
S4、对步骤S3得到的各通道校准后的回波数据抑制杂波,得到N/2组回波数据;
S5、选取一组步骤S4得到的回波数据,通过吕氏分布方法,估计动目标速度;
S6、根据步骤S5得到动目标速度,对S4杂波抑制后的数据,进行多通道重建、抑制动目标方位模糊,得到方位模糊抑制后的数据;
S7、根据步骤S6得到的方位模糊抑制后的数据,以及步骤S5估计的动目标速度,实现动目标聚焦成像。
进一步地,步骤S1所述单个通道的回波数据表达式为:
se,i(τ,t)=Aexp{-j4πf0Ri(t)/c}exp{jπKr(τ-2Ri(t)/c)2},i=1,…,N;
其中,A是复系数;i表示通道序号;τ,t分别为距离快时间和方位慢时间;c为电磁波速度;Kr=Br/Tr为距离调频斜率;j为虚数单位。
进一步地,所述步骤S2具体为:首先利用傅里叶变换将各通道回波数据逐行变换到距离频域;然后构建参考函数,并将该参考函数变换到频域;将各通道变换到距离频域的回波数据与变换到频域的参考函数相乘,最后经反傅里叶变换得到N个通道各自距离压缩后的回波数据。
进一步地,所述步骤S3具体为:将步骤S2得到的各通道回波数据进行方位傅里叶变换,根据傅里叶变换的时移性质,构建时延式,将变换到方位频域的各通道回波数据乘以时延式,得到各通道两两校准后的回波数据。
更进一步地,所述时延式为:
exp{-j2πf d/2/Vs};
其中,d/2是通道间隔,f表示方位向频率,Vs是平台速度。
进一步地,所述步骤S4具体为:对步骤S3得到的各通道校准后的回波数据,进行方位向和距离向离散采样,得到每个通道采样信号对通道校准后的两个通道都在(a,r)位置采样,组成二维的向量
其中,a,r分别表示方位向和距离向采样单元序号,r=1,...,Nr,a=1,...,Na;
分别选取待估计距离门的相邻L个距离门的数据作为样本,得到像素点(a,r)位置的杂波协方差矩阵R(a,r);
其中,上标H表示共轭转置;
根据计算得到的杂波协方差矩阵,得到像素点(a,r)位置的权矢量w(a,r);
其中,为动目标空域导引矢量,为空间频率,θ,分别为俯仰角和方位角;
根据计算出的权矢量,对像素点(a,r)位置的采样信号进行线性滤波,得到杂波抑制后的结果sstap(a,r);
进一步地,所述步骤S5具体为:
选取一组步骤S4得到的杂波抑制后的回波数据,根据吕氏分布方法,得到位于域的信号
其中,A表示信号幅值,q是一个恒定的时延参数,h为一个标量因子, 域表示以为横坐标,以为纵坐标的频域坐标系,表示多普勒中心频率变量,表示线性调频斜率变量,fd、γd分别表示实际估计出来的目标回波信号的多普勒中心频率值和线性调频斜率值;
根据fd和γd值,求解出动目标速度信息v=(vr,va)
其中,vr表示动目标距离向速度,va表示动目标方位向速度,(x,y)为动目标初始位置,X0为场景中心到平台轨迹中心的水平距离。
进一步地,所述步骤S6具体为:根据步骤S5得到的动目标距离向速度,计算动目标的多普勒频率;
其中,fa(-PRF/2<fa<PRF/2)是单通道模糊的方位向频率,b(=0,1,2)是划分的多普勒频带标号,PRF为脉冲重复频率,λ为波长,N是通道数,round表示取整;
然后构建多通道传递函数矩阵H(f),
其中,N表示通道数,f为方位频率;
对多通道传递函数矩阵H(f)求逆得到重建矩阵P(f),
P(f)=H-1(f);
将杂波抑制后的数据sstap沿方位向进行隔位补零,并变换到方位频域得到而后和重建矩阵P(f)相乘,经过方位向逆傅里叶变换可以得到重建结果为srec:
其中,为逆傅里叶变换算子。
更进一步地,所述传递函数矩阵H(f)中的第(x,y)个元素Hx,y通过下式计算得到:
其中,dy表示第y个通道和参考通道之间的间隔,y=1,..,N,fx=f+(x-1)PRF,x=1,..,N。
进一步地,所述步骤S7具体包括:根据步骤S5得到的动目标距离向和方位向速度;计算得到动目标在整个合成孔径时间内的运动轨迹;
根据计算得到的运动轨迹结合平台的运动轨迹得到动目标的斜距史;
根据动目标的斜距史以及步骤S6得到的数据,采用BP成像方法对动目标进行聚焦成像。
本发明的有益效果:本发明的一种高分宽幅SAR动目标速度估计与成像方法,将LVD方法应用在高分宽幅SAR条件下,进行动目标速度估计,利用估计得到的动目标距离向速度进行重建,抑制动目标方位模糊;利用估计得到的动目标距离向和方位向速度改进成像方法,抑制动目标成像方位散焦,最终得到聚焦效果更好的动目标成像;并且本申请只需要杂波抑制后的单个通道数据就可以估计动目标速度,有效提高了处理效率,大大减少了运算量。
附图说明
图1为本发明的处理流程框图。
图2为本发明实施例提供的经本申请方法处理得到的杂波抑制结果示意图;
其中,图2(a)为杂波抑制前,图2(b)为杂波抑制后。
图3为本发明实施例提供的经本申请方法处理获取的动目标二维成像结果示意图;
其中,图3(a)为成像结果散焦;图3(b)为重聚焦成像结果。
图4为本发明实施例提供的经本申请方法处理得到的动目标成像结果剖面图;
其中,图4(a)为成像结果散焦;图4(b)为重聚焦成像结果。
具体实施方式
为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐释。
为了方便描述本发明的内容,首先作以下术语定义:
定义1、HRWS多通道模型动目标斜距
第i个通道在方位慢时间tm时刻斜距为
其中,R0是最短斜距,vr是动目标距离向速度,Vs是平台速度,di是第i个通道与参考通道的间隔。
定义2、空时自适应处理(STAP)方法
STAP处理过程可以表示为:
sstap(a,r)=wH(a,r)sr(a,r) (2)
其中,w(a,r)为权矢量,sr为距离压缩后的数据,a,r分别表示距离向和方位向采样单元序号,sr(a,r)表示多个通道在像素点(a,r)的采样信号组成的向量最优权矢量w(a,r)由下式计算:
其中,为动目标空域导引矢量;为空间频率;θ,分别为俯仰角和方位角;
杂波协方差矩阵需要从多个距离门样本中估计得到,像素点(a,r)位置的杂波协方差矩阵R(a,r)估计式如下
其中,L为距离门样本数。
定义3、吕氏分布LVD方法
LVD估计的表达式可以表示为
其中,Ak表示幅值,q是一个恒定的时延参数,h为一个标量因子,为了使得估计的值令h=1,依据以下吕氏分布参考文献,qh=1,则设置q=1,分别表示第k个目标回波信号的多普勒中心频率和线性调频斜率,可以通过下式计算
其中,(x,y)为动目标初始位置,X0为场景中心到平台轨迹中心的水平距离。
具体的吕氏分布可参考文献“Lv X,Bi G,Wan C,et al.Lv's distribution:principle,implementation,properties,and performance[J].IEEE Transactions onSignal Processing,2011,59(8):3576-3591.”
定义4、动目标多普勒频率偏移
动目标重建的方位频率为
其中,fa(-PRF/2<fa<PRF/2)是单通道模糊的方位向频率,b(=0,1,2)是划分的多普勒频带,PRF为脉冲重复频率,λ为波长,N是通道数,round表示取整。
定义5、多通道重建
多通道重建矩阵为P(f)=H-1(f),其中H(f)为传递函数矩阵,形式如下
其中,N表示通道数,f为方位频率,传递函数矩阵H(f)中的第(x,y)个元素Hx,y通过下式计算:
其中,fx=f+(x-1)PRF,x=1,..,N,y=1,..,N,dy表示第y个通道和参考通道之间的间隔,各个通道间隔同为d,则三通道的间隔关系为d=(0,d,2d)。
如图1所示为本发明方法的流程图,本发明的技术方案为:一种高分宽幅SAR动目标速度估计与成像方法,包括:
S1、初始化方位多通道***参数,通过仿真得到多通道回波数据;
初始化方位多通道***参数包括:发射信号脉冲宽度,记为Tr;***脉冲重复频率,记为PRF;信号采样频率,记为Fs;通道数为N;发射信号带宽,记为Br;信号多普勒带宽为Bdop;平台速度为Vs;发射信号波长,记为λ;通道间隔为d;雷达中心频率,记为f0。
通过仿真得到N个通道各自的回波数据为:
se,i(τ,t)=Aexp{-j4πf0Ri(t)/c}exp{jπKr(τ-2Ri(t)/c)2},
其中,A是复系数;i表示通道序号,i=1,…,N;τ,t分别为距离快时间和方位慢时间;c为电磁波速度;Kr=Br/Tr为距离调频斜率;j为虚数单位,每个通道的回波数据维度均为:Na×Nr;
其中,Na,Nr分别为方位向和距离向采样点数。
S2、距离压缩
将步骤S1中得到的各通道回波数据se,i(τ,t),首先逐行变换到距离频域,然后构建参考函数并将该参考函数变换到频域后与变换到距离频域的回波信号相乘,最后经逆傅里叶变换(IFFT)得到各通道距离压缩后的回波数据
sr,i(τ,t)=Aexp{-j4πf0R(t)/c}。
S3、通道配准
对步骤S2得到的距离压缩结果sr(τ,t)进行傅里叶变换(FFT)得到方位频域数据,然后利用傅里叶变换时移性质,构建时延式exp{-j2πf d/2/Vs},
其中,d/2是通道间隔,f表示方位向频率,Vs是平台速度。
将变换到方位频域的N/2个通道的回波数据和时延式相乘得到两两校准后的N个通道的距离压缩数据
S4、STAP技术抑制杂波
对步骤S3得到的各通道校准后的回波数据,进行方位向和距离向离散采样,得到每个通道采样信号其中a,r分别表示方位向和距离向采样单元序号,r=1,...,Nr,a=1,...,Na,对通道校准后的两个通道都在(a,r)位置采样,组成二维的向量分别选取待估计距离门的相邻L个距离门的数据作为样本,参考式(4)估计得到像素点(a,r)位置的杂波协方差矩阵R(a,r),
其中,上标H表示共轭转置;表示像素点(a,l)处的采样信号;
然后参考式(3)相应求出像素点(a,r)位置的权矢量w(a,r),
参考式(2)对配准后的两通道数据进行线性滤波得到杂波抑制后的结果sstap(a,r),
因为这里是分别对N/2组配准的数据进行杂波抑制,所以最终会得到杂波抑制后的数据为sstap,k,k=1,…,N/2。
S5、LVD方法估计动目标速度
选取一组步骤S4得到的杂波抑制后的数据,参考式(5)进行LVD变换处理,得到位于域的信号
其中,A表示信号幅值,q是一个恒定的时延参数,h为一个标量因子, 域表示以为横坐标,以为纵坐标的频域坐标系,表示多普勒中心频率变量,表示线性调频斜率变量,fd、γd分别表示实际估计出来的目标回波信号的多普勒中心频率值和线性调频斜率值;
可以通过图像中的位置确定信号的fd和γd值,结合式(7)即可求解出动目标速度信息v=(vr,va)。
S6、多通道重建抑制动目标方位模糊
利用步骤S5估计的动目标速度vr,参考式(8)计算出定义4中动目标的多普勒频率fa,b,
然后参考式(9)构建多通道传递函数矩阵H(f),
对多通道传递函数矩阵H(f)求逆得到重建矩阵P(f),P(f)=H-1(f)
将步骤S4得到的杂波抑制后的数据sstap沿方位向进行隔位补零并变换到方位频域得到而后和重建矩阵P(f)相乘,经过方位向逆傅里叶变换可以得到重建结果为:
其中,为逆傅里叶变换算子,srec维度为(N/2*Na)×Nr。
S7、抑制动目标方位散焦得到高分辨成像
对步骤S6得到的重建数据做成像处理,首先根据步骤S5利用LVD方法估计得到的动目标速度v=(vr,va),计算出动目标在整个合成孔径时间内的运动轨迹Ptar(t)=vt,结合平台的运动轨迹PS(t)即可准确得到动目标的斜距史R(t)=|PS(t)-Ptar(t)|,最后利用BP成像方法对动目标进行聚焦成像。
如图2(a)所示,杂波抑制前动目标信息淹没在杂波中无法识别,如图2(b)所示,对配准后的数据进行STAP杂波抑制,杂波抑制后得到的动目标距离压缩结果,说明杂波已被有效抑制。图3(a)反映了动目标方位向速度对成像结果方位散焦的影响,如图4(a)中实线所示,传统动目标估计方法忽略了动目标方位向速度,直接成像导致最终的重建结果模糊抑制效果变差,成像方位向散焦。本发明的方法在利用LVD方法估计出动目标的两个方向的速度后,首先是利用距离向速度来重建抑制方位向模糊,然后在聚焦成像时结合动目标速度v=(vr,va)通过步骤S7得到完全聚焦的成像,对比结果如图3(b)以及图4(b)所示剖面结果,可以看到聚焦后动目标方位散焦被抑制。
现有HRWS SAR动目标成像中,动目标速度估计方法主要是穷举法和最小熵估计的方法。穷举法需要假设动目标速度,完成后续处理后看检测效果,再假设动目标速度如此循环取最终检测效果最好的情况,运算量大且估计精度较低;而最小熵方法运算量为也存在运算量大的问题。本发明将LVD方法应用到HRWS SAR中,只需要杂波抑制后的单个通道数据就可以估计动目标速度,有效提高了处理效率,运算量为并且可以同时估计得到动目标的距离向和方位向速度,在有效重建抑制动目标方位模糊的前提下,还可以解决动目标成像的方位散焦问题。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (10)
1.一种高分宽幅SAR动目标速度估计与成像方法,其特征在于,包括:
S1、初始化方位多通道***参数,通过仿真得到多通道回波数据;
所述初始化方位多通道***参数具体为:发射信号脉冲宽度,记为Tr;***脉冲重复频率,记为PRF;信号采样频率,记为Fs;通道数为N,N为偶数;发射信号带宽,记为Br;信号多普勒带宽为Bdop;平台速度为Vs;发射信号波长,记为λ;通道间隔为d;雷达中心频率,记为f0;
根据初始化的方位多通道***参数,通过仿真得到N个通道回波数据;所述每个通道回波数据均为Na×Nr维的矩阵;Na,Nr分别为方位向和距离向采样点数;
S2、将步骤S1得到的各通道的回波数据,分别进行距离压缩,得到各通道压缩后的回波数据;
S3、将步骤S2得到各通道回波数据进行两两校准,得到各个通道校准后的回波数据;
S4、对步骤S3得到的各通道校准后的回波数据抑制杂波,得到N/2组回波数据;
S5、选取一组步骤S4得到的回波数据,通过吕氏分布方法,估计动目标速度;
S6、根据步骤S5得到动目标速度,对S4杂波抑制后的数据,进行多通道重建、抑制动目标方位模糊,得到方位模糊抑制后的数据;
S7、根据步骤S6得到的方位模糊抑制后的数据,以及步骤S5估计的动目标速度,实现动目标聚焦成像。
2.根据权利要求1所述的一种高分宽幅SAR动目标速度估计与成像方法,其特征在于,步骤S1所述单个通道的回波数据表达式为:
se,i(τ,t)=Aexp{-j4πf0Ri(t)/c}exp{jπKr(τ-2Ri(t)/c)2},i=1,…,N;
其中,A是复系数;i表示通道序号;τ,t分别为距离快时间和方位慢时间;c为电磁波速度;Kr=Br/Tr为距离调频斜率;j为虚数单位。
3.根据权利要求1所述的一种高分宽幅SAR动目标速度估计与成像方法,其特征在于,所述步骤S2具体为:首先利用傅里叶变换将各通道回波数据逐行变换到距离频域;然后构建参考函数,并将该参考函数变换到频域;将各通道变换到距离频域的回波数据与变换到频域的参考函数相乘,最后经逆傅里叶变换得到N个通道各自距离压缩后的回波数据。
4.根据权利要求1所述的一种高分宽幅SAR动目标速度估计与成像方法,其特征在于,所述步骤S3具体为:将步骤S2得到的各通道回波数据进行方位傅里叶变换,根据傅里叶变换的时移性质,构建时延式,将变换到方位频域的各通道回波数据乘以时延式,得到各通道两两校准后的回波数据。
5.根据权利要求4所述的一种高分宽幅SAR动目标速度估计与成像方法,其特征在于,所述时延式为:
exp{-j2πf d/2/Vs};
其中,d/2是通道间隔,f表示方位向频率,Vs是平台速度。
6.根据权利要求1所述的一种高分宽幅SAR动目标速度估计与成像方法,其特征在于,所述步骤S4具体为:对步骤S3得到的各通道校准后的回波数据,进行方位向和距离向离散采样,得到每个通道采样信号对通道校准后的两个通道都在(a,r)位置采样,组成二维的向量
其中,a,r分别表示方位向和距离向采样单元序号,r=1,...,Nr,a=1,...,Na;
分别选取待估计距离门的相邻L个距离门的数据作为样本,得到像素点(a,r)位置的杂波协方差矩阵R(a,r);
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>L</mi>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mi>r</mi>
<mo>-</mo>
<mi>L</mi>
<mo>/</mo>
<mn>2</mn>
<mo>,</mo>
<mi>l</mi>
<mo>&NotEqual;</mo>
<mi>r</mi>
</mrow>
<mrow>
<mi>r</mi>
<mo>+</mo>
<mi>L</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<msub>
<mover>
<mi>s</mi>
<mo>~</mo>
</mover>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mover>
<mi>s</mi>
<mo>~</mo>
</mover>
<mi>r</mi>
<mi>H</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中,上标H表示共轭转置;表示像素点(a,l)处的采样信号;
根据计算得到的杂波协方差矩阵R(a,r),得到像素点(a,r)位置的权矢量w(a,r);
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>R</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>p</mi>
</mrow>
<mrow>
<msup>
<mi>p</mi>
<mi>H</mi>
</msup>
<msup>
<mi>R</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>p</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
其中,为动目标空域导引矢量,为空间频率,θ,分别为俯仰角和方位角;
根据计算出的权矢量w(a,r),对像素点(a,r)位置的采样信号进行线性滤波,得到杂波抑制后的结果sstap(a,r);
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>a</mi>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>w</mi>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>s</mi>
<mo>~</mo>
</mover>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
7.根据权利要求1所述的一种高分宽幅SAR动目标速度估计与成像方法,其特征在于,所述步骤S5具体为:
选取一组步骤S4得到的杂波抑制后的回波数据,根据吕氏分布方法,得到位于域的信号
其中,A表示信号幅值,q是一个恒定的时延参数,h为一个标量因子, 域表示以为横坐标,以为纵坐标的频域坐标系,表示多普勒中心频率变量,表示线性调频斜率变量,fd、γd分别表示实际估计出来的目标回波信号的多普勒中心频率值和线性调频斜率值;
根据fd和γd值,由下式求解出动目标速度信息v=(vr,va)
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>f</mi>
<mi>d</mi>
</msup>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mrow>
<msub>
<mi>&lambda;R</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>V</mi>
<mi>s</mi>
</msub>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>v</mi>
<mi>a</mi>
</msub>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>v</mi>
<mi>r</mi>
</msub>
<msub>
<mi>X</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>xX</mi>
<mn>0</mn>
</msub>
</mrow>
<msubsup>
<mi>R</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&gamma;</mi>
<mi>d</mi>
</msup>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mrow>
<msub>
<mi>&lambda;R</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<mn>2</mn>
<msub>
<mi>v</mi>
<mi>y</mi>
</msub>
<msub>
<mi>V</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>v</mi>
<mi>a</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>v</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<mfrac>
<msubsup>
<mi>X</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>R</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
其中,vr表示动目标距离向速度,va表示动目标方位向速度,(x,y)为动目标初始位置,X0为场景中心到平台轨迹中点的水平距离。
8.根据权利要求1所述的一种高分宽幅SAR动目标速度估计与成像方法,其特征在于,所述步骤S6具体为:根据步骤S5得到的动目标距离向速度,计算动目标的多普勒频率fa,b;
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<mo>&lsqb;</mo>
<mi>r</mi>
<mi>o</mi>
<mi>u</mi>
<mi>n</mi>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>v</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mo>&CenterDot;</mo>
<mi>P</mi>
<mi>R</mi>
<mi>F</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>b</mi>
<mo>-</mo>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</mfrac>
<mo>&rsqb;</mo>
<mi>P</mi>
<mi>R</mi>
<mi>F</mi>
<mo>;</mo>
</mrow>
其中,fa(-PRF/2<fa<PRF/2)是单通道的方位向频率,b(=0,1,2)是多普勒频带标号,PRF为脉冲重复频率,λ为波长,N是通道数,round表示取整;
然后根据下式构建多通道传递函数矩阵H(f):
<mrow>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>d</mi>
<mi>N</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<msub>
<mi>d</mi>
<mi>N</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>f</mi>
<mi>N</mi>
</msub>
<mo>,</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>f</mi>
<mi>N</mi>
</msub>
<mo>,</mo>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>f</mi>
<mi>N</mi>
</msub>
<mo>,</mo>
<msub>
<mi>d</mi>
<mi>N</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
其中,N表示通道数,f为方位频率;
对多通道传递函数矩阵H(f)求逆得到重建矩阵P(f)
P(f)=H-1(f);
将步骤S4得到的杂波抑制后的数据sstap沿方位向进行隔位补零,并变换到方位频域得到而后将和重建矩阵P(f)相乘,并经过方位向逆傅里叶变换得到多通道重建后的数据srec:
其中,为逆傅里叶变换算子。
9.根据权利要求8所述的一种高分宽幅SAR动目标速度估计与成像方法,其特征在于,所述传递函数矩阵H(f)中的第(x,y)个元素Hx,y通过下式计算得到:
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<msubsup>
<mi>&pi;d</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&lambda;R</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mo>&CenterDot;</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mfrac>
<msub>
<mi>d</mi>
<mi>y</mi>
</msub>
<mrow>
<mn>2</mn>
<msub>
<mi>V</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>f</mi>
<mi>x</mi>
</msub>
<mo>&rsqb;</mo>
<mo>;</mo>
</mrow>
其中,dy表示第y个通道和参考通道之间的间隔,y=1,..,N,fx=f+(x-1)PRF,x=1,..,N。
10.根据权利要求1所述的一种高分宽幅SAR动目标速度估计与成像方法,其特征在于,所述步骤S7具体包括:根据步骤S5得到的动目标距离向和方位向速度;计算得到动目标在整个合成孔径时间内的运动轨迹;
根据计算得到的运动轨迹结合平台的运动轨迹得到动目标的斜距史;
根据动目标的斜距史以及步骤S6得到的数据,采用BP成像方法对动目标进行聚焦成像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710415901.0A CN107229048B (zh) | 2017-06-06 | 2017-06-06 | 一种高分宽幅sar动目标速度估计与成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710415901.0A CN107229048B (zh) | 2017-06-06 | 2017-06-06 | 一种高分宽幅sar动目标速度估计与成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107229048A true CN107229048A (zh) | 2017-10-03 |
CN107229048B CN107229048B (zh) | 2020-07-31 |
Family
ID=59934851
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710415901.0A Expired - Fee Related CN107229048B (zh) | 2017-06-06 | 2017-06-06 | 一种高分宽幅sar动目标速度估计与成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107229048B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107831480A (zh) * | 2017-10-13 | 2018-03-23 | 西安电子科技大学 | 弹载雷达和差通道稳健自适应杂波抑制方法 |
CN108490417A (zh) * | 2018-05-22 | 2018-09-04 | 电子科技大学 | 一种精确的sar动目标参数估计方法 |
CN108549081A (zh) * | 2018-05-02 | 2018-09-18 | 北京空间飞行器总体设计部 | 一种高轨合成孔径雷达动目标速度检测方法 |
CN109725313A (zh) * | 2019-03-01 | 2019-05-07 | 中国科学院电子学研究所 | 一种sar海浪成像方法、***、电子设备和介质 |
CN109975804A (zh) * | 2019-03-04 | 2019-07-05 | 广东工业大学 | 一种多平台星座sar融合相参成像方法 |
CN110109102A (zh) * | 2019-04-04 | 2019-08-09 | 电子科技大学 | 一种sar运动目标检测与速度估计的方法 |
CN110261855A (zh) * | 2019-07-29 | 2019-09-20 | 上海无线电设备研究所 | 一种sar图像的近海岸地物杂波及其方位模糊抑制方法 |
CN110568434A (zh) * | 2019-07-15 | 2019-12-13 | 电子科技大学 | 一种多通道匀加速轨迹sar动目标二维速度估计方法 |
CN110646774A (zh) * | 2019-09-30 | 2020-01-03 | 中国人民解放军战略支援部队信息工程大学 | 基于乘积变尺度周期吕分布的机动目标相参检测方法及装置 |
CN111965643A (zh) * | 2020-07-10 | 2020-11-20 | 西安电子科技大学 | 一种斜视sar bp图像中运动舰船目标重聚焦方法 |
CN112098988A (zh) * | 2020-08-10 | 2020-12-18 | 西安电子科技大学 | 一种基于子空间投影的动目标径向速度估计方法 |
CN113484862A (zh) * | 2021-08-04 | 2021-10-08 | 电子科技大学 | 一种自适应的高分宽幅sar清晰重构成像方法 |
CN114488025A (zh) * | 2022-01-18 | 2022-05-13 | 西安电子科技大学 | 基于平滑吕变换的空间目标三维成像方法 |
CN117169887A (zh) * | 2023-11-03 | 2023-12-05 | 武汉能钠智能装备技术股份有限公司 | 基于方向确定的sar地面移动目标定位方法 |
-
2017
- 2017-06-06 CN CN201710415901.0A patent/CN107229048B/zh not_active Expired - Fee Related
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107831480B (zh) * | 2017-10-13 | 2021-06-01 | 西安电子科技大学 | 弹载雷达和差通道稳健自适应杂波抑制方法 |
CN107831480A (zh) * | 2017-10-13 | 2018-03-23 | 西安电子科技大学 | 弹载雷达和差通道稳健自适应杂波抑制方法 |
CN108549081B (zh) * | 2018-05-02 | 2020-05-08 | 北京空间飞行器总体设计部 | 一种高轨合成孔径雷达动目标速度检测方法 |
CN108549081A (zh) * | 2018-05-02 | 2018-09-18 | 北京空间飞行器总体设计部 | 一种高轨合成孔径雷达动目标速度检测方法 |
CN108490417A (zh) * | 2018-05-22 | 2018-09-04 | 电子科技大学 | 一种精确的sar动目标参数估计方法 |
CN109725313A (zh) * | 2019-03-01 | 2019-05-07 | 中国科学院电子学研究所 | 一种sar海浪成像方法、***、电子设备和介质 |
CN109975804A (zh) * | 2019-03-04 | 2019-07-05 | 广东工业大学 | 一种多平台星座sar融合相参成像方法 |
CN110109102A (zh) * | 2019-04-04 | 2019-08-09 | 电子科技大学 | 一种sar运动目标检测与速度估计的方法 |
CN110568434A (zh) * | 2019-07-15 | 2019-12-13 | 电子科技大学 | 一种多通道匀加速轨迹sar动目标二维速度估计方法 |
CN110568434B (zh) * | 2019-07-15 | 2022-06-03 | 电子科技大学 | 一种多通道匀加速sar动目标二维速度估计方法 |
CN110261855A (zh) * | 2019-07-29 | 2019-09-20 | 上海无线电设备研究所 | 一种sar图像的近海岸地物杂波及其方位模糊抑制方法 |
CN110646774A (zh) * | 2019-09-30 | 2020-01-03 | 中国人民解放军战略支援部队信息工程大学 | 基于乘积变尺度周期吕分布的机动目标相参检测方法及装置 |
CN110646774B (zh) * | 2019-09-30 | 2021-10-22 | 中国人民解放军战略支援部队信息工程大学 | 基于乘积变尺度周期吕分布的机动目标相参检测方法及装置 |
CN111965643A (zh) * | 2020-07-10 | 2020-11-20 | 西安电子科技大学 | 一种斜视sar bp图像中运动舰船目标重聚焦方法 |
CN111965643B (zh) * | 2020-07-10 | 2022-11-04 | 西安电子科技大学 | 一种斜视sar bp图像中运动舰船目标重聚焦方法 |
CN112098988A (zh) * | 2020-08-10 | 2020-12-18 | 西安电子科技大学 | 一种基于子空间投影的动目标径向速度估计方法 |
CN112098988B (zh) * | 2020-08-10 | 2023-02-10 | 西安电子科技大学 | 一种基于子空间投影的动目标径向速度估计方法 |
CN113484862A (zh) * | 2021-08-04 | 2021-10-08 | 电子科技大学 | 一种自适应的高分宽幅sar清晰重构成像方法 |
CN113484862B (zh) * | 2021-08-04 | 2023-10-17 | 电子科技大学 | 一种自适应的高分宽幅sar清晰重构成像方法 |
CN114488025A (zh) * | 2022-01-18 | 2022-05-13 | 西安电子科技大学 | 基于平滑吕变换的空间目标三维成像方法 |
CN114488025B (zh) * | 2022-01-18 | 2024-06-07 | 西安电子科技大学 | 基于平滑吕变换的空间目标三维成像方法 |
CN117169887A (zh) * | 2023-11-03 | 2023-12-05 | 武汉能钠智能装备技术股份有限公司 | 基于方向确定的sar地面移动目标定位方法 |
CN117169887B (zh) * | 2023-11-03 | 2024-04-19 | 武汉能钠智能装备技术股份有限公司 | 基于方向确定的sar地面移动目标定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107229048B (zh) | 2020-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107229048A (zh) | 一种高分宽幅sar动目标速度估计与成像方法 | |
CN111142105B (zh) | 复杂运动目标isar成像方法 | |
CN109143237B (zh) | 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法 | |
CN104597447B (zh) | 一种子孔径SAR大斜视改进Omega‑K成像方法 | |
Saeedi et al. | Synthetic aperture radar imaging using nonlinear frequency modulation signal | |
CN110632594B (zh) | 一种长波长星载sar成像方法 | |
CN104020471B (zh) | 一种基于分块处理的sar实时成像方法及*** | |
CN109270528B (zh) | 基于全解析距离模型的一站固定式双站sar成像方法 | |
CN104698431B (zh) | 基于模糊分量doa估计的多通道sar方位解模糊方法 | |
CN106842210A (zh) | 一种新的多子阵合成孔径声纳快速成像算法 | |
CN111505639A (zh) | 基于变重频采样模式的合成孔径雷达宽幅稀疏成像方法 | |
CN114545411B (zh) | 一种基于工程实现的极坐标格式多模高分辨sar成像方法 | |
CN104777479A (zh) | 基于多核dsp的前侧视sar实时成像方法 | |
CN105301589A (zh) | 高分辨宽测绘带sar地面运动目标成像方法 | |
CN103969644A (zh) | 一种多通道调频连续波sar成像方法 | |
Yu et al. | Ground moving target motion parameter estimation using Radon modified Lv's distribution | |
CN113759372A (zh) | 弹载大斜视小孔径多通道sar的成像方法 | |
CN111487614A (zh) | 基于子孔径的曲线航迹弹载sar波前重建成像方法及*** | |
CN104155653B (zh) | 一种基于特征距离子空间的sar后向投影成像方法 | |
CN109143235B (zh) | 一种双基前视合成孔径雷达地面运动目标检测方法 | |
CN111722225A (zh) | 基于先验相位结构信息的双基sar两维自聚焦方法 | |
CN104698448B (zh) | 运动平台下基于流形分离的共形阵列稳健估角方法 | |
CN103217677A (zh) | 一种基于联合检测量的单通道sar动目标检测方法 | |
CN110441771B (zh) | 基于方位时间重采样的高速机动飞行sar成像方法 | |
CN115453530B (zh) | 一种基于参数化模型的双基sar滤波反投影两维自聚焦方法 |
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 | ||
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: 20200731 |