CN115685200A - 一种高精度大前斜视sar成像运动补偿与几何校正方法 - Google Patents
一种高精度大前斜视sar成像运动补偿与几何校正方法 Download PDFInfo
- Publication number
- CN115685200A CN115685200A CN202211283419.3A CN202211283419A CN115685200A CN 115685200 A CN115685200 A CN 115685200A CN 202211283419 A CN202211283419 A CN 202211283419A CN 115685200 A CN115685200 A CN 115685200A
- Authority
- CN
- China
- Prior art keywords
- wave number
- azimuth
- motion compensation
- center
- domain
- 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.)
- Pending
Links
- 230000033001 locomotion Effects 0.000 title claims abstract description 100
- 238000003384 imaging method Methods 0.000 title claims abstract description 50
- 238000012937 correction Methods 0.000 title claims abstract description 39
- 208000004350 Strabismus Diseases 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000001228 spectrum Methods 0.000 claims abstract description 54
- 238000013507 mapping Methods 0.000 claims description 12
- 230000006835 compression Effects 0.000 claims description 11
- 238000007906 compression Methods 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 10
- 230000001131 transforming effect Effects 0.000 claims description 7
- 230000014509 gene expression Effects 0.000 claims description 5
- 230000021615 conjugation Effects 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 28
- 230000007547 defect Effects 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 39
- 238000004364 calculation method Methods 0.000 description 8
- 238000012545 processing Methods 0.000 description 8
- 230000003595 spectral effect Effects 0.000 description 7
- 230000005012 migration Effects 0.000 description 6
- 238000013508 migration Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 238000001914 filtration Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000005526 G1 to G0 transition Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000006243 chemical reaction Methods 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
- 238000001514 detection method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000007480 spreading Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- -1 OXY forms Chemical group 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000012292 cell migration Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation 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
- 230000003287 optical effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及机载合成孔径雷达成像领域,公开了一种高精度大前斜视SAR成像运动补偿与几何校正方法,通过对原始回波进行多普勒中心估计,并将方位波数谱搬移至基带,将数据进行距离向傅里叶变换,并进行一阶运动补偿,然后根据二维波数域参考函数,完成场景中心的完全聚焦和相同LOS斜距目标的一致压缩。通过沿视线方向的二维波数谱扩展Stolt插值,实现距离‑方位解耦,提高了波数谱利用率,避免传统波数域成像算法中需要方位补零去混叠、图像有黑边、旁瓣倾斜等缺陷;也能实现一阶一致补偿与二阶空变运动补偿,通过方位去斜将图像聚焦至波数域,避免空间域大量补零,将带有几何畸变的斜平面图像通过几何校正转换至地平面图像,实现高精度聚焦。
Description
技术领域
本发明涉及机载合成孔径雷达成像领域,具体涉及一种高精度大前斜视SAR成像运动补偿与几何校正方法。
背景技术
合成孔径雷达(SAR)是一种主动微波遥感设备,具有全天时、全天候对感兴趣区域进行高分辨率成像和运动目标检测定位能力。SAR图像反映的是被辐射场景的后向散射特性,可作为光学图像的有效补充。高分辨率机载SAR成像的意义是显而易见的,在军事方面可以对战场进行实时侦察,对典型目标进行精细化描述,为指挥系统及时地提供情报。SAR在地图测绘、地质勘探、灾害评估、农业普查和植被分类等民用领域也具有独特的优势,具有广泛的应用前景和发展潜力。
机载大前斜视SAR是指雷达波束的视线方向远远偏离飞行航迹的法向,前斜角(波束视线与航迹法线的夹角)可达60°甚至80°,如图2所示。传统的正侧视SAR只能对载机侧方区域进行成像,无法获取载机前斜视区域的信息。大前斜最主要的优势是能够提前探测目标,可以有较长的机动时间,而且前斜角越大,成像范围越大,可获取目标多角度散射特性,大幅提升雷达的侦察识别能力,因此大前斜SAR在战场侦察、目标识别、对地精确打击、提高战机生存能力等军事领域有着重要应用。上世纪80年代,美国就开始了大前斜SAR的研究,Raython公司研制的AN/APG-77有源相控阵雷达在F/A-22战斗机上验证了前斜SAR的可行性。
大前斜SAR的成像算法研究是大前斜SAR的核心问题,主要难点在于目标的距离徙动很大,导致目标的二维耦合非常严重,空变性较强。同时由于SAR平台运动的非理想特性以及运动参数测量误差,SAR运动补偿也是大前斜高分辨SAR成像必须解决的问题。另外,由于大前斜SAR特殊的成像几何,斜距平面和地距平面除了存在俯仰向的转动还存在方位向的旋转,地距图像的几何畸变较为严重,影响目标高精度定位。因此大前斜SAR成像技术的研究,可进一步提升雷达的战场侦察、提前感知和检测识别能力,为军事发展提供支撑。
经过国内外学者多年的努力,目前常用的斜视成像算法有后向投影(BackProjection,BP)算法、距离-多普勒(Range-Doppler,RD)算法、ω-k(Omega-K)类算法、Chirp Scaling(CS)算法、非线性CS(Non-linear CS,NLCS)算法、极坐标格式算法(PolarFormat Algorithm,PFA)等。时域类的BP以及改进的快速FFBP等算法在大幅宽、大场景成像中的计算效率依然偏低,且较难与运动补偿相结合,对惯导设备精度要求极高。添加二次距离压缩(Second Range Compression,SRC)处理步骤的RD成像算法也无法适用大斜视角SAR成像,CS算法利用调频变标实现SAR距离徙动校正,在大斜视角下也无法适用,相关的扩展算法如NLCS考虑了调频率的方位空变问题,能够适应60°左右斜视角,但都是对距离徙动进行了高阶近似,计算较为复杂且仍存有残余徙动,也没有考虑运动误差的影响。PFA算法可支持小场景的大前斜视成像,但是大场景高分辨条件下平面波前假设无法满足,空变滤波补偿又需要分块计算目标真实位置,计算量较大,且块内统一补偿仍存在近似误差。Omega-K算法在速度不空变的条件下是一种高精确性成像算法,扩展的Omega-K算法还支持二阶运动补偿,但是其成像平面均为斜平面,成像系方位与航迹一致,大斜视条件下有效场景为倾斜状,黑边较多,利用效率不高,且旁瓣为倾斜状,难以结合自聚焦算法;另外Omega-K算法将图像聚焦到二维空间域,当波束照射区域大于孔径长度支撑域时(例如条带模式),需要大量补零避免图像方位混叠,极大增加了处理量。
综上所述,现有大斜视成像技术主要存在以下缺点:
1)难以兼顾运算效率与成像精度,除Omega-K算法外均存在一定程度的近似,即使对距离历程进行高阶展开,在大场景高分辨条件下仍存在空变性,难以适用75°以上斜视角成像。
2)多数算法难以与运动补偿相结合,并考虑运动补偿的空变性。
3)Omega-K,CS类算法成像平面为斜距平面,图像方位向与航向一致,将目标校正到垂直航迹的最短斜距处,大斜视时方位补零点数较多,增大计算量,且图像黑边较大,需要补零避免方位混叠。
4)以视线(Line Of Sight,LOS)方向建立坐标系的图像多为斜平面图像,不利于目标高精度定位。
发明内容
有鉴于此,本发明提供一种高精度大前斜视SAR成像运动补偿与几何校正方法,解决由大前斜成像几何引起的大距离徙动、二维耦合严重、空变性强、图像难以聚焦等难点,并与测量惯导数据相结合,实现二阶距离空变运动补偿,在此基础上将存在俯仰和方位二维旋转的斜平面图像校正到地平面,实现图像的高精度聚焦与无畸变成像。
一种高精度大前斜视SAR成像运动补偿与几何校正方法,包括如下步骤:
步骤一:对原始回波进行多普勒中心估计得到斜平面斜视角,并构造波数谱中心校正函数;基于波数谱中心校正函数与原始回波,将方位波数谱搬移至基带,得到方位基带回波;
步骤二:对方位基带回波进行距离向傅里叶变换,并根据雷达参考斜距处的一阶运动补偿函数,完成一阶运动补偿;
步骤三:将一阶运动补偿后的数据变换到二维波数域,并根据二维波数域的景中心参考函数,完成场景中心聚焦与一致压缩;
步骤四:对一致压缩后的二维波数谱开展沿视线方向的扩展Stolt插值,得到沿视线和垂直视线方向的二维波数谱,实现距离-方位解耦;
步骤五:将扩展Stolt插值后的回波变换到二维空间域,并将原始回波的运动误差映射到视线方向,完成视线方向空变的二阶运动补偿;
步骤六:根据二阶运动补偿后的数据以及方位向的二次参考函数,进行空间域去斜,变换到方位波数域,完成目标聚焦,得到斜平面聚焦SAR图像;
步骤七:根据目标在沿视线方向的斜平面与地平面的几何映射关系,通过二维插值几何校正,得到无畸变的地平面大前斜视SAR图像。
进一步地,步骤一中通过将构造的波数谱中心校正函数与基带原始回波的二维空间域表达式相乘,将方位波数谱变换至基带,得到方位基带回波;其中波数谱中心校正函数Hazi_b(X,fdc)=exp{-jKxcX},Kxc为方位中心波数,Kxc=2πfdc/v,fdc为原始回波的方位多普勒中心频率,v为载机飞行速度。
进一步地,步骤二中对方位基带回波进行距离向傅里叶变换,并与参考成像中心处的一阶运动补偿函数相乘,完成一阶运动补偿;一阶运动补偿函数为Hmoco_1st(X,Kr)=exp{j2Kr·Rerr_1st(X)},其中Kr=2π(fc+fr)/c为距离向波数,c为光速,fc为中心频率,fr为发射LFM信号的距离向频率,-fs/2≤fr≤fs/2,fs为距离向采样率。
进一步地,步骤三中在一阶运动补偿后,将数据沿方位向做傅里叶变换,得到场景中心的二维波数谱,将二维波数谱与二维波数域参考函数相乘,完成场景中心聚焦与一致压缩;其中二维波数域参考函数为:
进一步地,步骤五中通过将沿LOS方向进行扩展的Stolt插值后的回波数据与方位空间域扩展函数相乘,然后做二维逆傅里叶变换到二维空间域,并将原始回波域的运动误差映射到LOS方向,完成LOS方向空变的二阶运动补偿;其中方位空间域扩展函数Hazi_ext(Ka;Bc)为,
进一步地,步骤六中通过在空间域将二阶运动补偿后的数据与二次相位参考函数相乘,在LOS⊥方向完成去斜;其中二次去斜参考函数为
其中,a为LOS⊥方向空变的目标位置,b为LOS方向空变的目标位置。
与现有技术相比,本发明的有益效果为:
1、本发明通过沿视线方向的二维波数谱扩展Stolt插值,实现距离-方位解耦,易与运动补偿结合,提高了波数谱利用率,避免传统波数域成像算法中需要方位补零去混叠、图像有黑边、旁瓣倾斜等缺陷;详细推导了原始运动误差在LOS坐标系下的变换形式,实现一阶一致补偿与二阶空变运动补偿;通过方位去斜将图像聚焦至波数域,避免空间域大量补零;将带有几何畸变的斜平面图像通过几何校正转换至地平面图像,实现高精度聚焦。
2、本发明可以解决由波束控制误差、天线安装误差、载机非理想运动误差引起的成像中心偏差问题,并与当前基于PFA算法的实时处理架构相结合,避免实时处理中的时间滞后性问题,实现当帧估计当帧补偿的实时处理能力,,可用于机载聚束SAR、条带SAR实时成像领域。
附图说明
为了更清楚地说明本申请实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1是实施例2中高精度大前斜视SAR成像运动补偿与几何校正方法流程图;
图2是实施例2中大前斜视SAR斜平面成像几何示意图;
图3是实施例2中大前斜视SAR地平面成像几何示意图;
图4为实施例2中两个点目标仿真结果图;
图5为实施例2中斜视角72°实测数据处理结果图;
图6为图5中局部放大示意图。
具体实施方式
下面结合附图对本申请实施例进行详细描述。
以下通过特定的具体实例说明本申请的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本申请的其他优点与功效。显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。本申请还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本申请的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
实施例1
一种高精度大前斜视SAR成像运动补偿与几何校正方法,包括如下步骤:
步骤一:对原始回波进行多普勒中心估计得到斜平面斜视角θs,并构造波数谱中心校正函数;基于波数谱中心校正函数与原始回波,将方位波数谱搬移至基带,得到方位基带回波;
主要是针对大斜视带来的方位波数谱混叠现象,对原始回波进行多普勒中心估计,构造波数谱中心校正函数,将二维波数谱中心搬移到零点,解决方位波数谱混叠问题。
步骤二:对方位基带回波进行距离向傅里叶变换,并根据雷达参考斜距处的一阶运动补偿函数,完成一阶运动补偿;
步骤三:将一阶运动补偿后的数据变换到二维波数域,并根据二维波数域的景中心参考函数,完成场景中心的完全聚焦和相同LOS斜距目标的一致压缩;
步骤四:对一致压缩后的二维波数谱开展沿视线方向的扩展Stolt插值,得到沿视线和垂直视线方向的二维波数谱,实现距离-方位解耦;完成距离徙动校正,同时保留了方位相位历史域数据,便于开展距离空变的二阶运动补偿;
步骤五:将扩展Stolt插值后的回波变换到二维空间域,并将原始回波的运动误差映射到视线方向,完成视线方向空变的二阶运动补偿;
步骤六:根据二阶运动补偿后的数据以及方位向的二次参考函数,进行空间域去斜,变换到方位波数域,将图像聚焦至波数域,得到聚焦良好的LOS斜平面图像;
步骤七:根据目标在沿视线方向的斜平面与地平面的几何映射关系,通过二维插值几何校正,得到无畸变的地平面大前斜视SAR图像。
在本实施例中,针对大斜视带来的方位波数谱混叠现象,对原始回波进行多普勒中心估计,构造波数谱中心校正函数,将二维波数谱中心搬移到零点,解决方位波数谱混叠问题;将回波数据变换到距离频域,完成参考斜距处的一阶运动补偿,然后根据场景中心位置建立二维波数域参考函数,完成场景中心的完全聚焦和相同LOS斜距目标的一致压缩。通过沿视线方向的二维波数谱扩展Stolt插值,实现距离-方位解耦,易与运动补偿结合,提高了波数谱利用率,避免传统波数域成像算法中需要方位补零去混叠、图像有黑边、旁瓣倾斜等缺陷;也能实现一阶一致补偿与二阶空变运动补偿,通过方位去斜将图像聚焦至波数域,避免空间域大量补零,将带有几何畸变的斜平面图像通过几何校正转换至地平面图像,实现高精度聚焦。
实施例2
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并结合图1-3,对本发明方法及效果作进一步的详细说明。
本实施例中,大前斜视SAR成像几何模型如图2、图3所示,载机飞行高度为H,雷达零时刻位置为O,载机以水平速度v沿X方向飞行。雷达工作在条带模式,雷达零时刻波束中心与地面交于点C,也是参考的场景中心。Y轴与X轴垂直,OXY构成成像斜平面,作用距离为R0,θs为斜平面斜视角,C点到载机航迹的最近距离为Yc,方位位置为Xc,波束与Z轴的夹角为P点是地平面任意一点目标,坐标为(Xp,Yp),OgXgYg构成地平面,θg为地平面斜视角。b轴为视线LOS方向,a轴为垂直视线LOS⊥方向,Oab位于成像斜平面。某一时刻载机偏离理想航迹的误差为(Δx,Δy,Δz)。
步骤一:对原始回波进行多普勒中心估计,得到斜平面斜视角θs,构造波数谱中心校正函数Hazi_b(X,fdc),与原始回波相乘,将方位波数谱变换至基带,避免频谱混叠。
雷达发射线性调频(LFM)信号,解调后的基带原始回波的二维空间域表达式为:
其中,K为线性调频信号调频斜率,λ为波长,r为距离向采样位置,j是复数虚部,c为光速,Rp(X)为X处雷达到目标的瞬时斜距。
对S(X,r)进行多普勒中心估计,得到斜平面斜视角θs。多普勒中心估计是估计回波的方位多普勒中心频率fdc,由于前斜角较大,fdc通常会跨越多个PRF,需要给定斜视角指令值计算PRF模糊数。本实施例采用测量方位样本间的相位变化来估计信号的多普勒频率。
根据fdc可得到斜平面斜视角由于原始回波方位波数谱不在基带,直接做方位FFT会使频谱折叠,方位谱坐标不连续,不利于后续插值,因此需要将S(X,r)与波数谱中心校正函数Hazi_b(X,fdc)相乘,将方位波数谱搬移到基带,结果为S1(X,r):
S1(X,r)=S(X,r)·Hazi_b(X,fdc)
Hazi_b(X,fdc)=exp{-jKxcX}
方位中心波数为Kxc,Kxc=2πfdc/v。
步骤二:对方位基带回波S1(X,r)进行距离向傅里叶变换,并与参考成像中心处的一阶运动补偿函数相乘,完成一阶运动补偿。参考成像中心为成像时设定的图像中心,如图2中的C点。通过该中心的坐标与载机位置计算一阶运动补偿函数,成像后,该点位于图像的中心。
因为载机的非理想直线运动引入了运动误差,使得理想的二维匹配滤波失配,影响目标的聚焦性能。此处只考虑运动误差的距离空变性,即认为沿视线方向同一距离的目标具有相同的运动误差,只需要沿方位向计算一阶运动误差并补偿。
场景中心处目标(Xc,Yc)的真实距离单元徙动(RCM)为
其中,X=m·v/PRF,m∈[-NX/2,NX/2-1]是载机理想的方位位置,m是方位离散采样序号,NX是方位采样点数,Δx(m),Δy(m),Δz(m)是m时刻的载机运动误差。
一阶运动补偿函数为Hmoco_1st(X,Kr)=exp{j2Kr·Rerr_1st(X)},其中,Kr=2π(fc+fr)/c为距离向波数,c为光速,fc为中心频率,fr为发射LFM信号的距离向频率,-fs/2≤fr≤fs/2,fs为距离向采样率。一阶运动补偿在距离向波数域对参考斜距处的目标进行了包络补偿和相位补偿,使得距离空变的残余运动误差不超过一个距离单元,后续只需要进行二阶相位补偿。
步骤三:将数据沿方位向做傅里叶变换,变换到二维波数域,并与二维波数域参考函数相乘,完成场景中心聚焦与一致压缩。该步骤在二维波数域完成了场景中心目标的匹配滤波,使场景中心完全聚焦。具体流程为:
经过一阶运动补偿后,利用POSP(驻定相位原理)可得场景中心的二维波数谱为
因此二维波数域参考函数为
P*(Kr)是P(Kr)的复数共轭。
步骤四:对景中心聚焦的二维波数谱开展沿视线方向的扩展Stolt插值,得到沿视线和垂直视线方向的二维波数谱,实现距离-方位解耦。
对于场景中任一点目标P,在斜平面OXY坐标系的坐标为(Xp,Yp),其回波经过上述步骤后,表达式为
随着斜视角的增大,波数谱支撑区的斜拉程度变大。传统的Stolt变换是将均匀的Kr-KX谱插值成均匀的KX-KY谱,这种方式下矩形支撑区内的有效频谱将会减少,降低频谱利用率,或者存在大量空白区域,降低运算效率数据。本发明采用沿视线的插值方式,将均匀的Kr-KX谱插值成均匀的Kb-Ka谱,其中Kb为视线(LOS)方向,即斜视角方向,Ka为垂直视线(LOS⊥)方向,与Kb形成右手法则,且Ka-Kb均位于斜平面。这种插值方式等效于沿着波数谱支撑区的斜拉方向进行,可以最大程度保留波数谱支撑区,提高频谱利用率。另外,对于大前斜视SAR,由于天线方向图效应,其波束照射区域在OXY平面是倾斜的条状,采用传统的插值方式,空间域图像只在视线方向有目标,其余均为波束无法覆盖的黑色区域,不利于判读解译,而且X方向的场景支撑域远大于合成孔径长度,需要额外补零避免混叠,降低了运算效率。采用沿LOS方向插值的方式图像黑边较少,波数谱支撑区近似矩形,二维旁瓣倾斜方向互相垂直,有利于后续的自聚焦处理。
沿视线方向的二维Stolt插值的波数映射关系为:
同样地,由Ka,Kb转换到Kr,KX的映射关系为
上述插值方式同时实现了RCMC和方位匹配滤波两个作用,但是RCMC后的相位历史域数据无法与常规运动补偿算法相结合,进行二阶距离空变相位补偿,因此需要采用扩展的沿LOS方向的Stolt插值方法,将均匀的Kr-KX谱插值成均匀的Kb_ext-Ka谱,映射关系如下:
由于直接进行二维插值计算较为复杂,本发明中分解为两个一维插值,即先将Kr-KX谱插值为Kr-Ka谱,再将Kr-Ka谱插值为Ka-Kb_ext谱,可以避免复杂的二次根式计算。
将Kr-Kx谱插值为Kr-Ka谱的方法为:每次选择一固定的Kr(n),根据待插值的Ka坐标计算待插值的KX_interp坐标:
其中n=0,1,L,Nr-1为距离向序号,Nr为距离向采样点数。
将向量S3(KX,Kr(n))插值成S3(KX_interp,Kr(n)),插值核可选为sinc核函数。
遍历所有n,得到均匀采样的S3(Ka,Kr)。
将Kr-Ka谱插值为Ka-Kb_ext谱的方法为:每次选择一固定的Ka(ma),根据待插值的Kb_ext坐标计算待插值的Kr_interp坐标:
其中ma=0,1,L,Na-1为LOS⊥方向序号,Na为LOS⊥向点数;
将向量S3(Kr,Ka(ma))插值成S3(Kr_interp,Ka(ma)),插值核可选为sinc核函数。遍历所有ma,得到均匀采样的S4(Ka,Kb_ext),即
步骤五:将沿LOS方向进行扩展的Stolt插值后的回波S4(Ka,Kb_ext)与方位空间域扩展函数Hazi_ext(Ka;Bc)相乘,做二维逆傅里叶变换到二维空间域,并将原始回波域的运动误差映射到LOS方向,完成LOS方向空变的二阶运动补偿。
点目标P的回波经过一阶运动补偿以及距离向匹配滤波后的表达式为
沿方位向做傅里叶变换后为
SMF(Kr,KX)
=∫SMF(Kr,X)exp(-jKXX)dX
=∫exp{-j[2Kr(Rp+ΔRerr_res)]}exp(-jKXX)dX
根据POSP原理,将上式中的积分相位对X求偏导并取零点
代入相位积分式,得到二维波数谱为
经过参考函数相乘后,回波变为
沿LOS方向的扩展Stolt插值后,回波变为
化简后为
其中
Yn=Yp-Yc,Xn=Xp-Xc
Bn=Yncosθs+Xnsinθs,An=-Ynsinθs+Xncosθs
(An,Bn)是目标P在以C为原点的Cab坐标系下的坐标。残余运动误差相位为
数据在LOS⊥方向的空间域表达式为
上式中,第一个相位项为扩展的Stolt插值引入的二次相位,与目标的LOS坐标Bn有关,在后续方位匹配滤波中能够补偿。第二个相位项为关于目标目标LOS向坐标Bn的线性相位项,在逆傅里叶变换后将目标压缩至正确的位置,第三项为残余运动误差,是二阶运动补偿重点关注的内容。
第一项
第二项中
在二阶运动补偿中,通常认为运动误差沿LOS⊥是空不变的,即可用Ap=0代入上式,有
分母
分子
因此
最终
通常残余运动误差ΔRerr_res远小于1个分辨单元,上式中关于Kb_ext的线性相位引起的距离徙动可忽略,因此在a位置处的空间域运动误差与X方向位于的运动误差相同。由于方位是离散采样,可以通过插值方式获取非整数网格上的运动误差。一种ΔRerr_res(a)的近似表达为
b为LOS方向空变的目标位置,Δbg(a)为轨迹偏离误差在地面LOS方向的投影,Δz(a)为轨迹偏离误差在高程方向的投影。
经过二阶运动补偿后的数据为
第六步:将二阶运动补偿后的数据在LOS⊥空间域与去斜函数相乘,变换到LOS⊥波数域,得到斜平面聚焦SAR图像。
由于Omega-K类算法均为将图像聚焦到二维空间域,而在实际雷达数据录取时,雷达波束照射范围通常大于孔径长度(空间域支撑区),使得图像在支撑区中产生混叠。常用避免混叠的方法是在空间域补零,扩展空间域支撑区,但是对于一些远距离成像场景,补零点数将远大于方位采样点数,造成运算量急剧增大。
针对这种情况,本文通过在空间域将二阶运动补偿后的数据与二次相位参考函数相乘,在LOS⊥方向完成去斜,接着沿LOS⊥方向做傅里叶变换,完成图像聚焦。
二阶运动补偿后的数据为
回波沿LOS⊥方向为近似线性调频信号,斜率为-Krc/(Bc+b),去斜后只剩关于目标位置An的线性相位,进行a方向的傅里叶变换后,结果为
其中Akn为目标在波数域的聚焦位置,有
因此,目标的斜平面波数域坐标Akn不仅是空间域坐标An的函数,还与LOS向位置Bn有关,存在一定的缩放,需要在后续的几何校正中补偿。
第七步:根据目标在沿LOS方向的斜平面与地平面的几何映射关系,通过二维插值完成几何校正,得到无畸变的地平面大前斜视SAR图像。
如图2和图3所示,目标在斜平面Cab的坐标为(An,Bn),斜平面OXY的坐标为(Xp,Yp),地平面Cagbg的坐标为(Agn,Bgn),地平面OgXgYg的坐标为(Xgp,Ygp),
有
Xp=Agn·cosθg+Bgn·sinθg+Xc
其中
因此
Bn=(Yp-Yc)·cosθs+(Xp-Xc)·sinθs
An=-(Yp-Yc)·sinθs+(Xp-Xc)·cosθs
根据坐标转换关系,通过二维插值将斜平面图像(Akn,Bn)变换到地平面图像(Agn,Bgn),完成几何校正,最终结果为
S7(ag,bg)=interp{S6(Ka,b)}。
表1示出了本发明仿真试验中涉及的参数设置,场景中按田字型布放9个点目标,相邻目标距离向间距为1931.9m,方位向间距为1488.9m。沿视线方向添加最大幅度为5m的正弦运动误差,该运动误差已超过多个距离分辨单元。场景方位幅宽大于孔径方位投影长度,本发明将图像聚焦到方位波数域,无需补零即可避免图像混叠。
表1仿真参数设置表
中心频率 | 10GHz | 带宽 | 150MHz |
载机速度 | 250m/s | 采样率 | 187.5MHz |
飞行高度 | 10km | PRF | 750Hz |
作用距离 | 80km | 方位分辨率 | 1m |
合成孔径时间 | 23s | 地平面斜视角 | 75° |
垂直航向幅宽 | 1km | 方位波束宽度 | 2.2° |
表2点目标仿真结果
图4示出了点目标1和点目标5的二维仿真结果图,仿真结果相关数据如表2,其中点目标1位于场景边缘,点目标5位于场景中心。各点目标经过二阶运动补偿后聚焦良好,位置准确,图像无几何畸变。由于设置为条带模式,场景边缘目标合成孔径角小,分辨率降低。
图5、图6示出了实测数据处理结果,前斜角约为72°,分辨率为3m。其中图5为大幅宽全场景图,图6为图5中建筑区白色框内局部放大子图,可见图像聚焦良好,全图无几何畸变。
以上仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以权利要求的保护范围为准。
Claims (6)
1.一种高精度大前斜视SAR成像运动补偿与几何校正方法,其特征在于,包括如下步骤:
步骤一:对原始回波进行多普勒中心估计得到斜平面斜视角,并构造波数谱中心校正函数;基于波数谱中心校正函数与原始回波,将方位波数谱搬移至基带,得到方位基带回波;
步骤二:对方位基带回波进行距离向傅里叶变换,并根据雷达参考斜距处的一阶运动补偿函数,完成一阶运动补偿;
步骤三:将一阶运动补偿后的数据变换到二维波数域,并根据二维波数域的景中心参考函数,完成场景中心聚焦与一致压缩;
步骤四:对一致压缩后的二维波数谱开展沿视线方向的扩展Stolt插值,得到沿视线和垂直视线方向的二维波数谱,实现距离-方位解耦;
步骤五:将扩展Stolt插值后的回波变换到二维空间域,并将原始回波的运动误差映射到视线方向,完成视线方向空变的二阶运动补偿;
步骤六:根据二阶运动补偿后的数据以及方位向的二次参考函数,进行空间域去斜,变换到方位波数域,完成目标聚焦,得到斜平面聚焦SAR图像;
步骤七:根据目标在沿视线方向的斜平面与地平面的几何映射关系,通过二维插值几何校正,得到无畸变的地平面大前斜视SAR图像。
2.根据权利要求1所述的高精度大前斜视SAR成像运动补偿与几何校正方法,其特征在于,步骤一中通过将构造的波数谱中心校正函数与基带原始回波的二维空间域表达式相乘,将方位波数谱变换至基带,得到方位基带回波;其中波数谱中心校正函数Hazi_b(X,fdc)=exp{-jKxcX},Kxc为方位中心波数,Kxc=2πfdc/v,fdc为原始回波的方位多普勒中心频率,v为载机飞行速度。
3.根据权利要求1所述的高精度大前斜视SAR成像运动补偿与几何校正方法,其特征在于,步骤二中对方位基带回波进行距离向傅里叶变换,并与参考成像中心处的一阶运动补偿函数相乘,完成一阶运动补偿;一阶运动补偿函数为Hmoco_1st(X,Kr)=exp{j2Kr·Rerr_1st(X)},其中Kr=2π(fc+fr)/c为距离向波数,c为光速,fc为中心频率,fr为发射LFM信号的距离向频率,-fs/2≤fr≤fs/2,fs为距离向采样率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211283419.3A CN115685200A (zh) | 2022-10-20 | 2022-10-20 | 一种高精度大前斜视sar成像运动补偿与几何校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211283419.3A CN115685200A (zh) | 2022-10-20 | 2022-10-20 | 一种高精度大前斜视sar成像运动补偿与几何校正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115685200A true CN115685200A (zh) | 2023-02-03 |
Family
ID=85065964
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211283419.3A Pending CN115685200A (zh) | 2022-10-20 | 2022-10-20 | 一种高精度大前斜视sar成像运动补偿与几何校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115685200A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116148856A (zh) * | 2023-04-17 | 2023-05-23 | 南京邮电大学 | 一种sar动目标二维自聚焦成像处理方法 |
CN117784077A (zh) * | 2024-02-27 | 2024-03-29 | 南京奥联智驾科技有限公司 | 一种基于频率积累的弱小目标检测方法、终端及介质 |
-
2022
- 2022-10-20 CN CN202211283419.3A patent/CN115685200A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116148856A (zh) * | 2023-04-17 | 2023-05-23 | 南京邮电大学 | 一种sar动目标二维自聚焦成像处理方法 |
CN117784077A (zh) * | 2024-02-27 | 2024-03-29 | 南京奥联智驾科技有限公司 | 一种基于频率积累的弱小目标检测方法、终端及介质 |
CN117784077B (zh) * | 2024-02-27 | 2024-05-24 | 南京奥联智驾科技有限公司 | 一种基于频率积累的弱小目标检测方法、终端及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20240319364A1 (en) | Segmented aperture imaging and positioning method of multi-rotor unmanned aerial vehicle-borne synthetic aperture radar | |
CN106249237B (zh) | 一种曲线轨迹下大斜视sar频域成像方法 | |
Frey et al. | Focusing of airborne synthetic aperture radar data from highly nonlinear flight tracks | |
Li et al. | Image coregistration in SAR interferometry | |
CN104007440B (zh) | 一种加速分解后向投影聚束合成孔径雷达成像方法 | |
CN104597447B (zh) | 一种子孔径SAR大斜视改进Omega‑K成像方法 | |
CN115685200A (zh) | 一种高精度大前斜视sar成像运动补偿与几何校正方法 | |
CN105759263B (zh) | 一种高分辨率大场景下的星载斜视sar雷达成像方法 | |
US5923278A (en) | Global phase unwrapping of interferograms | |
CN110488288B (zh) | 一种机载sar高分辨层析成像方法 | |
CN109738894B (zh) | 一种用于大视场合成孔径雷达的高斜视多角度成像方法 | |
CN109116321B (zh) | 一种星载干涉成像高度计的相位滤波方法及高度测量方法 | |
CN103454638B (zh) | 一种圆迹合成孔径雷达三维层析成像方法 | |
CN102866393B (zh) | 一种基于pos与dem数据的sar多普勒参数估计方法 | |
CN103576130B (zh) | 一种进动锥体的三维成像方法 | |
CN102608576B (zh) | 一种用于大前斜合成孔径雷达回波图像的几何校正方法 | |
CN114545411B (zh) | 一种基于工程实现的极坐标格式多模高分辨sar成像方法 | |
Zuo et al. | Unified coordinate system algorithm for terahertz video-SAR image formation | |
CN109358328A (zh) | 机动平台双基地前视sar的极坐标格式成像方法 | |
CN111208514A (zh) | 一种曲线运动轨迹SAR的切比雪夫斜距模型和Chirp Scaling成像方法 | |
CN114325704A (zh) | 一种基于波数谱拼接的合成孔径雷达快速时域成像方法 | |
CN112882029A (zh) | 基于快速后向投影的机动轨迹大斜视波束指向sar成像方法 | |
CN111880179A (zh) | 一种弹载弧线俯冲大斜视tops sar的成像方法 | |
Gong et al. | High resolution 3D InISAR imaging of space targets based on PFA algorithm with single baseline | |
CN113671497B (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 |