CN109358328B - 机动平台双基地前视sar的极坐标格式成像方法 - Google Patents
机动平台双基地前视sar的极坐标格式成像方法 Download PDFInfo
- Publication number
- CN109358328B CN109358328B CN201811310156.4A CN201811310156A CN109358328B CN 109358328 B CN109358328 B CN 109358328B CN 201811310156 A CN201811310156 A CN 201811310156A CN 109358328 B CN109358328 B CN 109358328B
- Authority
- CN
- China
- Prior art keywords
- wave number
- coordinate
- distance
- calculating
- polar
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 42
- 238000001228 spectrum Methods 0.000 claims abstract description 35
- 238000000034 method Methods 0.000 claims abstract description 25
- 238000012937 correction Methods 0.000 claims abstract description 22
- 238000005070 sampling Methods 0.000 claims abstract description 20
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 claims abstract description 19
- 238000012952 Resampling Methods 0.000 claims abstract description 12
- 238000013507 mapping Methods 0.000 claims abstract description 12
- 238000009825 accumulation Methods 0.000 claims abstract description 11
- 230000009466 transformation Effects 0.000 claims abstract description 8
- 230000006835 compression Effects 0.000 claims abstract description 4
- 238000007906 compression Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 17
- 238000005452 bending Methods 0.000 claims description 7
- 230000003595 spectral effect Effects 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims 1
- 238000011426 transformation method Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000001914 filtration Methods 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- 229910015234 MoCo Inorganic materials 0.000 description 2
- 239000000443 aerosol Substances 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 101710131169 Ribose-5-phosphate isomerase A 1 Proteins 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000926 separation 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
-
- 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/9058—Bistatic or multistatic SAR
-
- 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
- G01S13/9011—SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth
-
- 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/9043—Forward-looking SAR
-
- 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/9052—Spotlight mode
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开一种机动平台双基地前视聚束合成孔径雷达的极坐标格式成像方法,采用波数谱支撑域仿射变换的方法,将MP‑BFSAR的平行四边形波数谱支撑域变换为水平准矩形,然后通过距离向chirp‑z变换实现波数谱距离向极坐标到直角坐标的重采样,之后进行距离压缩和场景中心运动补偿,再在方位向上采用升采样映射积累的方法实现波数谱方位向极坐标到直角坐标的重采样;然后进行二维逆傅里叶变换,得到粗成像结果;最后通过散焦校正得到最终的成像结果;本发明所收集的波数谱的利用率相比于现有的方法明显提高。
Description
技术领域
本发明属于双基地前视合成孔径雷达(Bistatic Forward-looking SyntheticAperture Radar,BFSAR)成像领域,特别涉及一种机动平台双基地前视聚束合成孔径雷达极坐标格式算法。
背景技术
由于发射站和接收站的空间分离性,双基地SAR突破了SAR的前视约束,使接收平台能够在飞行方向上对前视地形进行成像,这种独特的特性使得BFSAR成为一种很有前景的传感技术。BFSAR能够提供全天时、全天候的前视高分辨率图像。
由于其双基地前视配置,MP-BFSAR的回波表现出严重的二维空间误差,其波数谱支撑区域通常是平行四边形而不是矩形,因此很难有效地使用所收集的波数谱;由于MP-BFSAR波束通常工作在模式,以期在许多应用中持续获得感兴趣区域的高分辨率图像,这导致回波的复杂时频特性,如多普勒混叠。
PFA算法(Parabola fitting algorithm:抛物线拟合算法)处理k空间(即波数域)中的回波。由于其独特的特性,如无模式引起的多普勒混叠、低平台轨迹限制和低计算复杂性等,所以PFA适用于处理SAR数据。目前PFA已广泛应用于单基地SAR和双基地SAR。
由于PFA采用平面波假设,但是实际波前存在弯曲,当成像场景较大时,需进行波前弯曲校正,对于大场景可能是无效的。
文献“Polar format algorithm for bistatic SAR,IEEETrans.Aerosp.Electron.Syst.,vol.40,no.4,pp.1147–1159,Oct.2004”和“Space-variant filtering for wavefront curvature correction in polar formattedbistatic SAR image,IEEE Trans.Aerosp.Electron.Syst.,vol.48,no.2,pp.940–950,Apr.2012”提出了k-空间坐标轴旋转的方法,将倾斜的平行四边形的波数谱支撑域转变成水平的平行四边形,然而,仅当支撑区域是矩形或准矩形时,k空间轴旋转才有效。因此,不适用于支撑区域是平行四边形的MP-BFSAR;文献“Polar format algorithm wavefrontcurvature compensation under arbitrary radar flight path,IEEE Geosci.RemoteSens.Lett.,vol.9,no.3,pp.526–530,May 2012”中,为了补偿由波前弯曲引起的相位误差,提出了一种空变波前弯曲补偿滤波器,它利用了极坐标到直角坐标重采样操作的内同性特性,然而,空变波前弯曲补偿滤波器是隐含的,其计算涉及二维插值,这意味着它在实践中不方便且效率低;文献“Space-variant filtering for wavefront curvaturecorrection in polar formatted bistatic SAR image,IEEETrans.Aerosp.Electron.Syst.,vol.48,no.2,pp.940–950,Apr.2012”中,提出了一种解析的空变波前弯曲补偿滤波器,然而,它是通过假设平台以匀速直线运动轨迹移动而得出的,因此它不适用于MP-BFSAR,因为MP-BFSAR的运动轨迹是高度非线性的。
发明内容
为解决上述技术问题,本发明提出一种机动平台双基地前视SAR的极坐标格式成像方法,采用KAM的方法,将平行四边形的波数谱支撑域变换到水平准矩形,使得收集的波数谱的利用率明显提高。
本发明采用的技术方案为:机动平台双基地前视SAR的极坐标格式成像方法,采用波数谱支撑域仿射变换(k-set affine mapping,KAM)的方法,将MP-BFSAR的平行四边形波数谱支撑域变换为水平准矩形,然后通过距离向chirp-z变换(Chirp-Z Transform,CZT)实现波数谱距离向极坐标到直角坐标的重采样,之后进行距离压缩和场景中心运动补偿,再在方位向上采用升采样映射积累(Upsampling mapping accumulation,UMA)的方法实现波数谱方位向极坐标到直角坐标的重采样;然后进行二维逆傅里叶变换(Inverse fastfourier transform,IFFT),得到粗成像结果;最后通过散焦校正得到最终的成像结果。
波数谱仿射变换(KAM),具体过程如下:
A1、计算回波时延梯度;
A2、计算回波多普勒梯度;
A3、根据回波时延梯度与回波多普勒梯度计算仿射矩阵及其逆矩阵;
A4、计算x坐标和y坐标方向的波数;
A5、根据x坐标和y坐标方向的波数以及仿射矩阵,计算KAM后波数分布范围;
A6、根据波数分布范围与仿射矩阵的逆矩阵,计算图像的坐标范围。
距离向CZT实现波数谱距离向极坐标到直角坐标的重采样;具体为:
B1、回波信号建模;
B2、对步骤B1建模的回波信号进行距离向CZT;
方位向升采样映射积累,具体包括以下过程:
C1、将步骤B2所得信号方位向升采样Nu倍;
C2、对经步骤C1之后的信号进行方位向信号积累。
方位向信号积累
C2、计算从η轴到kq轴的映射关系;
C3、将步骤C1得到的信号按照步骤C2的映射关系进行积累。
通过散焦校正得到最终的成像结果,具体包括以下过程:
D1、计算波前弯曲导致的相位误差;
D2、根据步骤D1的相位误差,计算计算波前弯曲导致的二次相位误差;
D3、根据二次相位误差构建散焦矫正滤波器;
D4、将粗成像结果划分为若干子图像;
D5、将个子图像乘上步骤D3构建的散焦矫正滤波器进行散焦矫正;
D6、将步骤D5得到的所有散焦校正后的子图像进行拼合,得到最终成像结果。
本发明的有益效果:本发明的方法,采用KAM的方法,将平行四边形的波数谱支撑域变换到水平准矩形,使得所收集的波数谱的利用率提高了约60%;由于CZT和UMA都能够通过傅里叶变换(Fast fourier transform,FFT)、IFFT和相位相乘高效实现,故本发明提出的PFA较为高效;本发明提出的解析的散焦校正滤波器相比现有的数值滤波器更加方便和高效,且本发明的散焦校正使得PFA算法能够适用于大的成像场景。
附图说明
图1为本发明方法的成像处理流程图;
图2为本发明所述机动平台双基地SAR回波录取几何模型;
图3为KAM前后的波数谱支撑域;
其中,图3(a)为KAM之前的波数谱支撑域;图3(b)为KAM之后的波数谱支撑域;
图4为点目标的成像结果;
图5为图4中点目标成像结果的放大等高图;
图6为面目标的成像结果。
具体实施方式
为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐释。
如图1所示为本发明的方案流程图,包括:采用波数谱支撑域仿射变换(k-setaffine mapping,KAM)的方法,将MP-BFSAR的平行四边形波数谱支撑域变换为水平准矩形,然后通过距离向chirp-z变换(Chirp-Z Transform,CZT)实现波数谱距离向极坐标到直角坐标的重采样,之后进行距离压缩和场景中心运动补偿,再在方位向上采用升采样映射积累(Upsampling mapping accumulation,UMA)的方法实现波数谱方位向极坐标到直角坐标的重采样;然后进行二维逆傅里叶变换(Inverse fast fourier transform,IFFT),得到粗成像结果;最后通过散焦校正得到最终的成像结果。
具体包括以下5个步骤:
1、波数谱仿射变换(KAM):
A1、回波时延梯度为
其中,c表示光速,分别表示方位零时刻发射站和接收站到场景中心视线方向单位矢量,rT0=(xT,yT,zT)、rR0=(xR,yR,zR)分别表示初始时刻发射站和接收站的位置矢量,RT0=|rT0|、RR0=|rR0|分别表示发射站和接收站在方位零时刻到场景中心的距离。
时延梯度在地面的投影为
Ggτ=PGτ
A2、计算多普勒梯度
回波多普勒梯度为
其中,λ表示载波波长,vT=(vTx,vTy,vTz)、vR=(vRx,vRy,vRz)分别表示发射站和接收站在方位零时刻的速度,I表示单位矩阵。
多普勒梯度在地面的投影为
Ggd=PGd
A3、计算仿射矩阵及其逆矩阵
对仿射矩阵求逆得
A4、计算x坐标和y坐标方向的波数
ΔR(η;x,y)=Rbi(η;x,y)-Rbi(0;x,y)表示差分距离,对差分距离作一阶泰勒近似:
ΔR(η;x,y)=RSAG(η;x,y)-RSAG(0;x,y)≈C10x+C01y,
其中
得到波数k在x轴和y轴的投影分别为kx=kC10,ky=kC01,其中
k=2π(fr+fc)/c
即,x坐标和y坐标方向的波数分别为
kx=kC10
ky=kC01
Rbi(η;x,y)=RT(η;x,y)+RR(η;x,y)表示双基地距离和,η表示慢时间,
RT(η;x,y)、RR(η;x,y)分别表示发射站和接收站距离历史,aT=(aTx,aTy,aTz)、aR=(aRx,aRy,aRz)分别表示发射站和接收站的加速度,rP=(x,y,0)表示场景内任意一个点目标的位置坐标。
A5、计算KAM后波数分布范围
其中,kxmin,kxmax分别表示kx的最小值和最大值,kymin,kymax分别表示ky的最小值和最大值;
KAM后距离向波数变量kp的分布范围为[kpmin,kpmax],其中
kpmin=min{kp1,kp2,kp3,kp4}
kpmax=max{kp1,kp2,kp3,kp4}
KAM后方位向波数变量kq的分布范围为[kqmin,kqmax],其中
kqmin=min{kq1,kq2,kq3,kq4}
kqmax=max{kq1,kq2,kq3,kq4}
A6、计算输出图像的坐标范围,具体如下:
其中,xmin,xmax分别表示原始场景的横坐标的最小值和最大值,ymin,ymax分别表示原始场景的纵坐标的最小值和最大值;
输出图像的距离向坐标范围为[pmin,pmax],其中
pmin=min{p1,p2,p3,p4}
pmax=max{p1,p2,p3,p4}
输出图像的方位向坐标范围为[qmin,qmax],其中
qmin=min{q1,q2,q3,q4}
qmax=max{q1,q2,q3,q4}
2、距离向CZT:
B1、回波信号建模
距离频域方位基带回波可建模为如下形式
其中,Kr表示发射脉冲调频率,τ表示快时间,c表示光速,λ表示载波波长。
B2、距离向CZT
对回波信号通过距离向CZT求回波在复平面单位圆如下位置的频谱
zi=exp{j(φ0+iΔφ)},i=0,1,2,...,Np
其中,
距离向CZT后信号可表示为
其中,kp表示KAM后的距离向波数变量。
3、匹配滤波(MF,Match filter)和场景中心运动补偿(MoCo to SC),即对距离向CZT后所得结果乘上如下相位函数:
MF和MoCo to SC之后的信号可近似表示为
其中,p,q分别表示输出图像的距离向和方位向坐标变量。
4、方位向升采样映射积累:
C1、将步骤3所得信号方位向升采样Nu倍,记所得信号为S3(kp,η(m);p,q),其中m表示升采样后η的采样序号
其中,Na表示升采样前η的采样点数,Nu一般取值为4或8;
C2、计算从η轴到kq轴的映射关系:
C3、方位向信号积累
将C1所得的信号在方位向按C2所得的映射关系进行积累:
其中,n、Nq分别表示信号在kq轴的采样点序号和采样点数。
至此,得到了在(kp,kq)域的信号,其可近似表示为
S4(kp,kq;p,q)≈exp{-j(kpp+kqq)}
C4、对C3所得信号进行二维逆傅里叶变换(IFFT),即可得到粗成像结果。
5、对于场景较大的情况,波前弯曲导致的边缘目标散焦效应不可忽略,需要对粗成像结果进行散焦校正:
D1、计算波前弯曲导致的相位误差,具体如下:
波前弯曲导致的相位误差为
ΦWC(k,η;x,y)=-kΔR(η;x,y)+kxx+kyy
D2、计算波前弯曲导致的二次相位误差
由隐函数定理和复合函数求导法则计算波前弯曲导致的二次相位误差为
D3、构建散焦校正滤波器
D4、子图像划分
将粗成像结果沿q轴划分为Nsub个子图像,再将每一个子图像通过关于q的傅里叶变换(Fast fourier transform,FFT)变换到(p,kq)域,其中Nsub为满足下式的一个正整数
其中qm由下面的方程解得
D5、子图像散焦校正
把每一个子图像乘上如下散焦校正滤波器
HRF(kq;p,qc)=exp{-jk2D(p,qc)}
其中,qc为子图像中心的q轴坐标,然后进行关于kq的IFFT,得到(p,q)域上的散焦校正后的子图像。
D6、子图像拼合
将D5所得的所有散焦校正后的子图像组合在一起,得到最终成像结果。
以下结合具体数据,对本发明的技术效果进行阐述。
如表1所示为本实施例采用的***参数,图2为本实施方式机动双基地SAR几何配置示意图,本实施方式中假定地面上有49个均匀分布的点目标,其x轴和y轴之上各个点目标之间的距离均为800米。
表1***参数
基于表1的***参数取值进行各项计算,具体为:
时延梯度在地面的投影为:Ggτ=PGτ=1.0×10-8[-0.2711,0.2634]T
回波多普勒梯度计算结果为:
多普勒梯度在地面的投影为:Ggd=PGd=[-0.0001,0.0458]T
耦合角计算结果为:
φc=|θgτ-θgd|=45.7°
仿射矩阵为:
对仿射矩阵求逆,得到:
对波数谱进行如下仿射变换
kp的最大值和最小值计算结果分别为:
kpmin=28.38rad/m、kpmax=30.99rad/m,
输出图像的距离向采样点数为:Np=3328。
Nu=8,升采样后η的采样序号为:
脉冲重复间隔为:PRI=1/PRF=4ms,
信号在kq轴的采样点数为:Nq=2224,
另外,由kqmax=-0.9044rad/m,kqmax=0.8413rad/m,得到kq轴的采样间隔为
信号在kq轴的采样序号为:
将粗聚焦图沿q轴划分为Nsub=8个子图像,并将每一个子图像通过关于q的FFT变换到(p,kq)域;
其中,D(p,qc)=7.31×10-6p2+1.02×10-5pqc+1.03×10-5q2,
p∈[-4.0020×103,3.9996×103]m;qc为各个子图像中心对应的q轴坐标,其取值集合为[-3500,-2500,-1500,-500,500,1500,2500,3500]m;
根据以上计算结果得到如图3、图4、图5、图6所示的仿真图像;KAM前后的波数谱支撑域仿真结果如图3所示,KAM之后,点目标的波数谱支撑域的利用效率由33%增加到了98%,由此可以看出,本发明所提出的方法可以有效提高波数谱的利用效率。本发明的点目标成像结果如图4所示;相应的点目标成像结果放大等高图如图5所示,可以看出点目标成像结果的等高图均匀规则分布,其峰值旁瓣比(PSLR)、积分旁瓣比(ISLR)与理论值-13.26dB、-9.96dB非常接近,说明本发明提出的算法对点目标回波的聚焦成像效果良好;面目标成像结果如图6所示,图6中可以看出面目标清晰可辨,说明本发明提出的算法对面目标回波的聚焦成像效果良好。综合图4、图5和图6可以看出,本发明所提出的方法可以有效地对机动平台双基地前视SAR目标回波进行聚焦成像。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (7)
1.机动平台双基地前视聚束合成孔径雷达的极坐标格式成像方法,其特征在于,
首先,通过波数谱仿射变换,将MP-BFSAR的平行四边形波数谱支撑域变换为水平准矩形;其次,通过距离向chirp-z变换实现波数谱距离向极坐标到直角坐标的重采样;再次,进行距离压缩和场景中心运动补偿;从次,波数谱方位向极坐标到直角坐标的重采样;然后,进行二维逆傅里叶变换,得到粗成像结果;最后,通过散焦校正得到最终的成像结果。
2.根据权利要求1所述的机动平台双基地前视聚束合成孔径雷达的极坐标格式成像方法,其特征在于,波数谱仿射变换,具体过程如下:
A1、计算回波时延梯度;
A2、计算回波多普勒梯度;
A3、根据回波时延梯度与回波多普勒梯度计算仿射矩阵及其逆矩阵;
A4、计算x坐标和y坐标方向的波数;
A5、根据x坐标和y坐标方向的波数以及仿射矩阵,计算KAM后波数分布范围;
A6、根据波数分布范围与仿射矩阵的逆矩阵,计算图像的坐标范围。
3.根据权利要求2所述的机动平台双基地前视聚束合成孔径雷达的极坐标格式成像方法,其特征在于,距离向CZT实现波数谱距离向极坐标到直角坐标的重采样;具体为:
B1、回波信号建模;
B2、对步骤B1建模的回波信号进行距离向CZT。
4.根据权利要求3所述的机动平台双基地前视聚束合成孔径雷达的极坐标格式成像方法,其特征在于,波数谱方位向极坐标到直角坐标的重采样,具体为:在方位向上采用升采样映射积累;包括以下步骤:
C1、将步骤B2所得信号方位向升采样Nu倍;
C2、对经步骤C1之后的信号进行方位向信号积累。
5.根据权利要求4所述的机动平台双基地前视聚束合成孔径雷达的极坐标格式成像方法,其特征在于,步骤C1采样后的信号为:S3(kp,η(m);p,q);
其中,kp表示波数谱支撑域仿射变换后距离向波数变量,η表示慢时间,p表示图像的距离向坐标,q表示图像的方位向坐标,m表示升采样后η的采样序号。
6.根据权利要求5所述的机动平台双基地前视聚束合成孔径雷达的极坐标格式成像方法,其特征在于,方位向信号积累具体为:
首先,计算从η轴到kq轴的映射关系;
然后,将步骤C1得到的信号按照映射关系进行积累。
7.根据权利要求6所述的机动平台双基地前视聚束合成孔径雷达的极坐标格式成像方法,其特征在于,通过散焦校正得到最终的成像结果,具体包括以下过程:
D1、计算波前弯曲导致的相位误差;
D2、根据步骤D1的相位误差,计算波前弯曲导致的二次相位误差;
D3、根据二次相位误差构建散焦矫正滤波器;
D4、将粗成像结果划分为若干子图像;
D5、将个子图像乘上步骤D3构建的散焦矫正滤波器进行散焦矫正;
D6、将步骤D5得到的所有散焦校正后的子图像进行拼合,得到最终成像结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811310156.4A CN109358328B (zh) | 2018-11-06 | 2018-11-06 | 机动平台双基地前视sar的极坐标格式成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811310156.4A CN109358328B (zh) | 2018-11-06 | 2018-11-06 | 机动平台双基地前视sar的极坐标格式成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109358328A CN109358328A (zh) | 2019-02-19 |
CN109358328B true CN109358328B (zh) | 2020-06-30 |
Family
ID=65344090
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811310156.4A Active CN109358328B (zh) | 2018-11-06 | 2018-11-06 | 机动平台双基地前视sar的极坐标格式成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109358328B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10495719B1 (en) * | 2013-09-24 | 2019-12-03 | Charles Alphonse Uzes | System for receiving communications |
CN109917361B (zh) * | 2019-04-02 | 2023-04-25 | 电子科技大学 | 一种基于双基雷达的三维未知场景成像方法 |
CN114167394A (zh) * | 2021-12-03 | 2022-03-11 | 东南大学 | 基于相位梯度误差估计技术的线谱特征增强方法及其*** |
CN114185047B (zh) * | 2021-12-09 | 2023-06-27 | 电子科技大学 | 一种基于最优极坐标变换的双基sar动目标重聚焦方法 |
CN114236543B (zh) * | 2021-12-16 | 2023-03-21 | 电子科技大学 | 一种机动平台双基前视sar轨迹设计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105223572A (zh) * | 2015-09-14 | 2016-01-06 | 北京航空航天大学 | 一种基于pfa算法的正前视双基sar成像处理方法 |
CN106990397A (zh) * | 2017-06-07 | 2017-07-28 | 电子科技大学 | 双基地前视sar非***距离徙动校正方法 |
CN107589421A (zh) * | 2017-10-31 | 2018-01-16 | 西安电子科技大学 | 一种阵列前视sar成像方法 |
CN108318880A (zh) * | 2018-01-26 | 2018-07-24 | 西安电子科技大学 | 一种运动误差参数化的极坐标sar成像方法 |
-
2018
- 2018-11-06 CN CN201811310156.4A patent/CN109358328B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105223572A (zh) * | 2015-09-14 | 2016-01-06 | 北京航空航天大学 | 一种基于pfa算法的正前视双基sar成像处理方法 |
CN106990397A (zh) * | 2017-06-07 | 2017-07-28 | 电子科技大学 | 双基地前视sar非***距离徙动校正方法 |
CN107589421A (zh) * | 2017-10-31 | 2018-01-16 | 西安电子科技大学 | 一种阵列前视sar成像方法 |
CN108318880A (zh) * | 2018-01-26 | 2018-07-24 | 西安电子科技大学 | 一种运动误差参数化的极坐标sar成像方法 |
Non-Patent Citations (2)
Title |
---|
"The Polar Imaging Algorithm for Forward-looking Bistatic SAR";Jin-ping SUN et al.;《7th European Conference on Synthetic Aperture Radar》;20080605;全文 * |
"多模式SAR成像及参数估计方法研究";吴玉峰;《中国博士学位论文全文数据库 信息科技辑》;20160315(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109358328A (zh) | 2019-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109358328B (zh) | 机动平台双基地前视sar的极坐标格式成像方法 | |
CN109143237B (zh) | 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法 | |
CN105842694B (zh) | 一种基于ffbp sar成像的自聚焦方法 | |
CN104391297A (zh) | 一种划分子孔径pfa雷达成像方法 | |
CN111999734B (zh) | 一种宽带斜视聚束sar两步成像方法 | |
CN111856461A (zh) | 基于改进pfa的聚束sar成像方法及其dsp实现 | |
CN114384520B (zh) | 一种机动平台对海面船只精细化雷达成像方法 | |
CN116299465B (zh) | 一种基于子空间时频映射的双基sar后向投影成像方法 | |
CN107607952A (zh) | 基于电磁涡旋波的三维合成孔径雷达成像方法 | |
CN105699970B (zh) | 一种用于合成孔径雷达大场景欺骗干扰的逆距离徙动方法 | |
CN115685200A (zh) | 一种高精度大前斜视sar成像运动补偿与几何校正方法 | |
CN111208514A (zh) | 一种曲线运动轨迹SAR的切比雪夫斜距模型和Chirp Scaling成像方法 | |
CN109143236B (zh) | 适用于复杂飞行轨迹的双基聚束sar大场景成像方法 | |
CN114325704A (zh) | 一种基于波数谱拼接的合成孔径雷达快速时域成像方法 | |
CN111722225A (zh) | 基于先验相位结构信息的双基sar两维自聚焦方法 | |
CN104931965B (zh) | 一种固定发射站双基正侧视sar成像方法 | |
CN111665501B (zh) | 基于改进cbp的mimo雷达二维成像方法 | |
Li et al. | Inverse-mapping filtering polar formation algorithm for high-maneuverability SAR with time-variant acceleration | |
CN111127334B (zh) | 基于rd平面像素映射的sar图像实时几何校正方法及*** | |
CN114185047B (zh) | 一种基于最优极坐标变换的双基sar动目标重聚焦方法 | |
CN115015920A (zh) | 一种基于距离空变频谱矫正的快速后向投影成像方法 | |
CN108469612B (zh) | 基于等效斜距的双基时变加速度前视sar成像方法 | |
Gui et al. | Dynamic ISAR imaging method for multiple moving vehicles based on OMP-CADMM | |
CN112859018B (zh) | 一种基于图像几何校正的视频sar成像方法 | |
Gong et al. | High resolution 3D InISAR imaging of space targets based on PFA algorithm with single baseline |
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 |