CN106802416A - 一种快速因式分解后向投影sar自聚焦方法 - Google Patents
一种快速因式分解后向投影sar自聚焦方法 Download PDFInfo
- Publication number
- CN106802416A CN106802416A CN201710093572.2A CN201710093572A CN106802416A CN 106802416 A CN106802416 A CN 106802416A CN 201710093572 A CN201710093572 A CN 201710093572A CN 106802416 A CN106802416 A CN 106802416A
- Authority
- CN
- China
- Prior art keywords
- formula
- phase error
- theta
- phi
- aperture
- 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
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/9004—SAR image acquisition techniques
- G01S13/9019—Auto-focussing of the SAR signals
-
- 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/9017—SAR image acquisition techniques with time 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
-
- 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/40—Means for monitoring or calibrating
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)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种快速因式分解后向投影SAR自聚焦方法,结合快速因式分解后向投影成像算法(FFBP)建立相位误差优化模型,并利用坐标下降方法和割线方法来求解该优化模型,有效的解决了相位误差补偿问题。与现有技术相比,本发明能够更加精确地估计相位误差以及得到良好聚焦的SAR图像,解决了现有FFBP自聚焦方法的运算量高、精度低、适应小的问题,从而实现结合FFBP成像的精确运动误差估计和良好聚焦。
Description
技术领域
本发明属于雷达技术领域,具体涉及一种基于最大图像锐度的快速因式分解后向投影SAR自聚焦方法的设计。
背景技术
合成孔径雷达(Synthetic Aperture Radar,SAR)是一种具有高分辨率的成像雷达,与光学传感器相比,SAR具有全天时全天候工作能力的独特优点。随着SAR技术的发展和提高,其分辨率越来越高,目前已接近或超光学成像的分辨率,因而被广泛的应用于地球遥感、海洋研究、资源勘探、灾情预报和军事侦察等领域。
基于时域的快速因式分解后向投影(FFBP)成像方法不存在近似处理,具有高精度成像能力和良好的相位保持特性,几乎适用于所有合成孔径雷达系统,包括双/多基地合成孔径雷达等具有复杂成像几何的合成孔径雷达系统;并且解决了传统后向投影(BP)成像方法的大的计算复杂度问题。
对于SAR的高分辨成像,自聚焦是关键步骤之一。FFBP成像算法的应用已经逐渐趋于成熟,但是对于FFBP的自聚焦还有待进一步的研究。传统的自聚焦方法,如相位梯度自聚焦、子视图相关等都是建立在图像域和距离压缩相位历史域的傅里叶变换关系上;然而这种变换关系对于FFBP却是不成立的。为了解决这个问题,文献Zhang Lei,Li Haolin,XingMengdao,Bao Zheng,“Integrating Autofocus Techniques With Fast FactorizedBack-projection For High-Resolution Spotlight SAR Imaging,”IEEE Geoscienceand remote sensingpp:1394-1398,2013引入虚拟极坐标系,使二者的傅里叶变换关系近似成立,同时对误差函数进行多项式拟合,提出子孔径图像偏移自聚焦,但是却受限于实际大场景且对高频误差进行估计时性能下降;文献Jan Torgrimsson,Patrik Dammert,HansHellsten,Lars M.H.Ulander,“Factorized Geometrical Autofocus for SyntheticApertuer Radar Processing,”IEEE Transactions on geoscience and remotesensing,2014和文献Jan Torgrimsson,Lars M.H Ulander,Patrik Dammert,HansHellsten,“Factorized geometrical autofocus:On the geometry search,”IEEETransactions on geoscience and remote sensing,2016从几何参数出发,建立评价图像聚焦性能的代价函数,利用优化求解的方法对几何参数进行估计,但是涉及多维参数优化求解,运算量大,速度慢,且当涉及复杂运动时估计精度降低,适用性小。
发明内容
本发明的目的是为了解决现有技术中针对快速因式分解后向投影成像方法中的运动误差补偿问题,提出了一种基于最大图像锐度的快速因式分解后向投影SAR自聚焦方法。
本发明的技术方案为:一种快速因式分解后向投影SAR自聚焦方法,包括以下步骤:
S0、系统参数初始化;
S1、快速因式分解后向投影成像;
S2、相位误差估计。
进一步地,步骤S0具体为:
以场景参考目标点O=[0 0 0]T为原点建立坐标系,场景中任意像素点位置记为P=[x y 0]T,雷达平台理想位置记为PA=[0 vt h]T;
其中,v为载机沿着y轴的理想飞行速度,h为载机沿着y轴的理想飞行高度,t为方位慢时间;
任意像素点距雷达平台的理想瞬时距离记为R(t)=||P-PA||;当飞行轨迹存在偏差时,雷达平台实际位置记为其中[Δx(t) Δy(t) Δz(t)]T表示载机位置偏差;
任意像素点距雷达平台的实际瞬时距离记为即:
将公式(1)展开可以得到:
其中ΔR表示距离误差;得到对应相位误差表达式:
式中λ表示雷达发射信号对应波长。
进一步地,步骤S1具体为:
根据因式分解的原理确定最佳的初始孔径长度l0以及每次合并的子图像个数I,即分解因子;
以各个初始子孔径的中心为原点建立极坐标系,划分出子图像的像素点坐标(r,θ)的取值范围,有其中,Rsrn为近端距离,Rsrf为远端距离,r表示像素点到子孔径坐标原点的距离,θ表示孔径矢量和距离矢量r之间的夹角,Θ表示积累角;极坐标系下的瞬时距离表示为:
其中,d表示雷达平台到对应坐标系原点的距离;
确定初始子图像的角度向分辨率
同时根据信号带宽确定距离向分辨率ρr;方
位向采样点个数为N,初始成像子孔径个数为K;通过公式(5)(6)对成像场景中的每一点(r,
θ)进行初始成像:
其中,w(θ)为天线方向函数,s(x,τ)为获取的回波,δ(·)为单位冲激函数,x为雷达在方位向的位置,τ为对应的快时间,下标q表示方位向采样点,k表示第k个子孔径,且有q∈[1,N],k∈[1,K],上标(1)表示第一级(初始)成像;
考虑确定雷达位置对应确定的相位误差,记飞行轨迹偏差导致的方位向相位误差为φ={φ1,φ2.....φN},代入公式(6),得到其中,j表示为虚数单位;则初始成像更新为:
对于i+1级某个成像网格中的某点(r,θ),计算其在第i级第k幅子图像中的位置(r′,θ′),将第k幅图像在该点的结果采用插值的方式累加至新的i+1级图像f(r,θ)(i+1)中,根据公式(8)-(10)对子孔径成像结果进行逐级合并:
rcos(θ)-r′cosθ′=dk (8)
此时,dk表示为当前子孔径对应的上一级合并时第k个子孔径距离该子孔径几何中心的距离;
每进行一级合并,将I幅子图像生成次一级子图像,第i+1幅子图像分辨率和第i级子图像分辨率存在以下的关系:
其中,ρθm为最终分辨率;得到最终图像:
其中,表示多次合并插值之后的值,表示第q个方位向采样点相位误差的估计值,用于补偿相位误差,则f(F)为最终补偿后的图像。
进一步地,步骤S2具体为:
记图像锐度为:
其中i表示第i个像素点,记最终得到的图像f(F)→f,对应的共轭表示为f*,则vi=fifi *;建立带有未知相位误差的最大图像锐度的优化模型来求解相位误差即是对φ的估计,最终使得图像良好聚焦,表示为:
使用坐标下降的迭代方法来求解最优解;在坐标下降的算法中,每一次迭代都使得优化变量按序得以更新,而保持其他变量固定不变;记为第i次迭代的第n个相位误差校正量,则第i+1次迭代表示为:
最终图像表示如下:
其中,x表示除了第n次采样点的其他所有采样点的后向投影之和,y表示对当前采样点校正后的后向投影;将公式(16)代入(13)得到:
其中,xi和yi为第i个像素点在公式(16)下的对应表示;xi *和yi *为xi和yi对应的共轭表示,
因此在坐标下降处理下将多维优化求解简化为单变量求极值问题由公式(13)得到关于φ的导数:
将公式(19)(20)代入公式(18)最终得到:
其中:
利用割线法对单变量优化问题进行求解,具体步骤如下:
T1、选择初始点(φ0)1,参数ε>0,γ>0,令(φ1)1=(φ0)1-γφ′((φ0)1),k:=0;
T2、若|φ′((φk)i)|>ε,则利用公式(23)和(24)更新αk,否则停止运算;其中,上标i表示第i个方位向采样点;
T3、通过公式(25)和(26)计算(φk+1)i,k=k+1,返回步骤T2;
公式(23)-(26)中的下标k表示对于第i个方位向采样点的第k次迭代求解;通过上述步骤,实现对相位误差进行精确估计并实现良好成像。
本发明的有益效果是:本发明结合快速因式分解后向投影成像算法(FFBP)建立相位误差优化模型,并利用坐标下降方法和割线方法来求解该优化模型,有效的解决了相位误差补偿问题。与现有技术相比,本发明能够更加精确地估计相位误差以及得到良好聚焦的SAR图像,解决了现有FFBP自聚焦方法的运算量高、精度低、适应小的问题,从而实现结合FFBP成像的精确运动误差估计和良好聚焦。
附图说明
图1为本发明提供的一种快速因式分解后向投影SAR自聚焦方法流程图。
图2为本发明实施例采用的快速因式分解后向投影成像具体实现结构示意图。
图3为本发明实施例采用的雷达系统结构示意图。
图4为本发明实施例的场景点目标分布图。
图5为本发明实施例的成像散焦结果图。
图6为本发明实施例的相位误差估计图。
图7为本发明实施例的成像聚焦结果图。
具体实施方式
下面结合附图对本发明的实施例作进一步的说明。
本发明提供了一种基于最大图像锐度的快速因式分解后向投影SAR自聚焦方法,如图1所示,包括以下步骤:
S0、系统参数初始化。
以场景参考目标点O=[0 0 0]T为原点建立坐标系,场景中任意像素点位置记为P=[x y 0]T,雷达平台理想位置记为PA=[0 vt h]T。
其中,v为载机沿着y轴的理想飞行速度,h为载机沿着y轴的理想飞行高度,t为方位慢时间。
任意像素点距雷达平台的理想瞬时距离记为R(t)=||P-PA||;当飞行轨迹存在偏差时,雷达平台实际位置记为其中[Δx(t) Δy(t) Δz(t)]T表示载机位置偏差。
任意像素点距雷达平台的实际瞬时距离记为即:
将公式(1)展开可以得到:
其中ΔR表示距离误差;得到对应相位误差表达式:
式中λ表示雷达发射信号对应波长。
S1、快速因式分解后向投影成像。
根据因式分解的原理确定最佳的初始孔径长度l0以及每次合并的子图像个数I,即分解因子。
以各个初始子孔径的中心为原点建立极坐标系,划分出子图像的像素点坐标(r,θ)的取值范围,有其中,Rsrn为近端距离,Rsrf为远端距离,r表示像素点到子孔径坐标原点的距离,θ表示孔径矢量和距离矢量r之间的夹角,Θ表示积累角;极坐标系下的瞬时距离表示为:
其中,d表示雷达平台到对应坐标系原点的距离。
确定初始子图像的角度向分辨率同时根据信号带宽确定距离向分辨率ρr;方位向采样点个数为N,初始成像子孔径个数为K;通过公式(5)(6)对成像场景中的每一点(r,θ)进行初始成像:
其中,w(θ)为天线方向函数,s(x,τ)为获取的回波,δ(·)为单位冲激函数,x为雷达在方位向的位置,τ为对应的快时间,下标q表示方位向采样点,k表示第k个子孔径,且有q∈[1,N],k∈[1,K],上标(1)表示第一级(初始)成像。
考虑确定雷达位置对应确定的相位误差,记飞行轨迹偏差导致的方位向相位误差为φ={φ1,φ2.....φN},代入公式(6),得到其中,j表示为虚数单位;则初始成像更新为:
对于i+1级某个成像网格中的某点(r,θ),计算其在第i级第k幅子图像中的位置(r′,θ′),将第k幅图像在该点的结果采用插值的方式累加至新的i+1级图像f(r,θ)(i+1)中,从图2所示的几何关系中可以得到公式(8)-(10):
rcos(θ)-r′cosθ′=dk (8)
此时,dk表示为当前子孔径对应的上一级合并时第k个子孔径距离该子孔径几何中心的距离。
根据公式(8)-(10)对子孔径成像结果进行逐级合并,每进行一级合并,将I幅子图像生成次一级子图像,第i+1幅子图像分辨率和第i级子图像分辨率存在以下的关系:
其中,ρθm为最终分辨率;得到最终图像:
其中,表示多次合并插值之后的值,表示第q个方位向采样点相位误差的估计值,用于补偿相位误差,则f(F)为最终补偿后的图像。
S2、相位误差估计。
记图像锐度为:
其中i表示第i个像素点,记最终得到的图像f(F)→f,对应的共轭表示为f*,则vi=fifi *;建立带有未知相位误差的最大图像锐度的优化模型来求解相位误差即是对φ的估计,最终使得图像良好聚焦,表示为:
使用坐标下降的迭代方法来求解最优解;在坐标下降的算法中,每一次迭代都使得优化变量按序得以更新,而保持其他变量固定不变;记为第i次迭代的第n个相位误差校正量,则第i+1次迭代表示为:
最终图像表示如下:
其中,x表示除了第n次采样点的其他所有采样点的后向投影之和,y表示对当前采样点校正后的后向投影;将公式(16)代入(13)得到:
其中,xi和yi为第i个像素点在公式(16)下的对应表示;xi *和yi *为xi和yi对应的共轭表示,
因此在坐标下降处理下将多维优化求解简化为单变量求极值问题由公式(13)得到关于φ的导数:
将公式(19)(20)代入公式(18)最终得到:
其中:
利用割线法对单变量优化问题进行求解,具体步骤如下:
T1、选择初始点(φ0)1,参数ε>0,γ>0,令(φ1)1=(φ0)1-γφ′((φ0)1),k:=0。
T2、若|φ′((φk)i)|>ε,则利用公式(23)和(24)更新αk,否则停止运算;其中,上标i表示第i个方位向采样点;
T3、通过公式(25)和(26)计算(φk+1)i,k=k+1,返回步骤T2。
公式(23)-(26)中的下标k表示对于第i个方位向采样点的第k次迭代求解;通过上述步骤,即可实现对相位误差进行精确估计并实现良好成像。
下面以一个具体实施例对本发明提供的一种基于最大图像锐度的快速因式分解后向投影SAR自聚焦方法作进一步描述:
S0、系统参数初始化。
本发明实施例采用的SAR几何结构如图3所示,系统仿真参数如下表所示:
仿真参数 | 取值 |
载频 | 9.6e9Hz |
带宽 | 200e6Hz |
脉冲重复频率 | 800Hz |
时间宽度 | 5e-6s |
相位误差 |
本发明实施例采用的目标场景如图4所示,图中的黑色圆点为布置于地面上的9个点目标,这9个点分布于x轴和y轴上,沿x方向(切航迹)间隔10m,沿y方向(沿航迹)间隔10m,平台沿y轴运动,且速度v=30m/s。
波束中心位于场景坐标原点处时记为零时刻,零时刻的雷达平台位置坐标为PA=(0 1000 1000)T,场景中心坐标为(0 0 0)T;场景中任一点目标的位置坐标为P[x y 0]T;在MATLAB平台上对成像区域内的任意目标计算其距离历程R,仿真出雷达回波数据,记为s0(t,τ)。
S1、对回波数据进行快速因式分解后向投影成像。
对回波数据s0(t,τ)按照下式进行距离压缩:
s′0(t,τ)=IFFT(FFT(s0(t,τ))×FFT(fτ(τ))) (27)
其中Kr为距离向调频斜率,取Kr=4.0e+13,fc为载频。
在距离压缩后的回波数据中加入相位误差Δφ,误差具体表达式如上表所示,记加入相位误差的回波数据为s(t,τ)。
取分解因子I为2,初始孔径长度为l0为0.0375m,确定初始子图像的角度向分辨率为0.0768rad,同时根据信号带宽确定距离向分辨率ρr为1.5m。在极坐标系下按照公式(4)计算每一个像素点距雷达平台的双程距离和,按照公式(5)(6)进行初始子图像成像。所有子孔径成像完成后按照公式(8)-(10)对子孔径成像结果进行逐级合并,下一级成像分辨率计算如公式(11)所示,直至得到最终成像结果,带有相位误差的成像结果如图5所示。
S2、相位误差估计。
根据公式(13)计算图像锐度,得到设置初始误差以及参数ε=0.001,γ=0.001。将ε,γ代入步骤T1中计算出将其代入公式(18)计算根据公式(23)和(24)更新参数(α0)1,代入(25)和(26)计算下一个重复以上步骤计算下一轮迭代的最终得到所有N个误差估计值相位误差估计结果如图6所示,误差补偿后成像结果如图7所示。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (4)
1.一种快速因式分解后向投影SAR自聚焦方法,其特征在于,包括以下步骤:
S0、系统参数初始化;
S1、快速因式分解后向投影成像;
S2、相位误差估计。
2.根据权利要求1所述的快速因式分解后向投影SAR自聚焦方法,其特征在于,所述步骤S0具体为:
以场景参考目标点O=[0 0 0]T为原点建立坐标系,场景中任意像素点位置记为P=[xy 0]T,雷达平台理想位置记为PA=[0 vt h]T;
其中,v为载机沿着y轴的理想飞行速度,h为载机沿着y轴的理想飞行高度,t为方位慢时间;
任意像素点距雷达平台的理想瞬时距离记为R(t)=||P-PA||;当飞行轨迹存在偏差时,雷达平台实际位置记为其中[Δx(t) Δy(t) Δz(t)]T表示载机位置偏差;
任意像素点距雷达平台的实际瞬时距离记为即:
将公式(1)展开可以得到:
其中ΔR表示距离误差;得到对应相位误差表达式:
式中λ表示雷达发射信号对应波长。
3.根据权利要求1所述的快速因式分解后向投影SAR自聚焦方法,其特征在于,所述步骤S1具体为:
根据因式分解的原理确定最佳的初始孔径长度l0以及每次合并的子图像个数I,即分解因子;
以各个初始子孔径的中心为原点建立极坐标系,划分出子图像的像素点坐标(r,θ)的取值范围,有其中,Rsrn为近端距离,Rsrf为远端距离,r表示像素点到子孔径坐标原点的距离,θ表示孔径矢量和距离矢量r之间的夹角,Θ表示积累角;极坐标系下的瞬时距离表示为:
其中,d表示雷达平台到对应坐标系原点的距离;
确定初始子图像的角度向分辨率同时根据信号带宽确定距离向分辨率ρr;方位向采样点个数为N,初始成像子孔径个数为K;通过公式(5)(6)对成像场景中的每一点(r,θ)进行初始成像:
其中,w(θ)为天线方向函数,s(x,τ)为获取的回波,δ(·)为单位冲激函数,x为雷达在方位向的位置,τ为对应的快时间,下标q表示方位向采样点,k表示第k个子孔径,且有q∈[1,N],k∈[1,K],上标(1)表示第一级成像;
考虑确定雷达位置对应确定的相位误差,记飞行轨迹偏差导致的方位向相位误差为φ={φ1,φ2.....φN},代入公式(6),得到其中,j表示为虚数单位;则初始成像更新为:
对于i+1级某个成像网格中的某点(r,θ),计算其在第i级第k幅子图像中的位置(r′,θ′),将第k幅图像在该点的结果采用插值的方式累加至新的i+1级图像f(r,θ)(i+1)中,根据公式(8)-(10)对子孔径成像结果进行逐级合并:
r cos(θ)-r′cosθ′=dk (8)
此时,dk表示为当前子孔径对应的上一级合并时第k个子孔径距离该子孔径几何中心的距离;
每进行一级合并,将I幅子图像生成次一级子图像,第i+1幅子图像分辨率和第i级子图像分辨率存在以下的关系:
其中,ρθm为最终分辨率;得到最终图像:
其中,表示多次合并插值之后的值;表示第q个方位向采样点相位误差的估计值,用于补偿相位误差,则f(F)为最终补偿后的图像。
4.根据权利要求3所述的快速因式分解后向投影SAR自聚焦方法,其特征在于,所述步骤S2具体为:
记图像锐度为:
其中i表示第i个像素点,记最终得到的图像f(F)→f,对应的共轭表示为f*,则建立带有未知相位误差的最大图像锐度的优化模型来求解相位误差即是对φ的估计,最终使得图像良好聚焦,表示为:
使用坐标下降的迭代方法来求解最优解;在坐标下降的算法中,每一次迭代都使得优化变量按序得以更新,而保持其他变量固定不变;记为第i次迭代的第n个相位误差校正量,则第i+1次迭代表示为:
最终图像表示如下:
其中,x表示除了第n次采样点的其他所有采样点的后向投影之和,y表示对当前采样点校正后的后向投影;将公式(16)代入(13)得到:
其中,xi和yi为第i个像素点在公式(16)下的对应表示;xi *和yi *为xi和yi对应的共轭表示,
因此在坐标下降处理下将多维优化求解简化为单变量求极值问题由公式(13)得到关于φ的导数:
将公式(19)(20)代入公式(18)最终得到:
其中:
利用割线法对单变量优化问题进行求解,具体步骤如下:
T1、选择初始点(φ0)1,参数ε>0,γ>0,令(φ1)1=(φ0)1-γφ′((φ0)1),k:=0;
T2、若|φ′((φk)i)|>ε,则利用公式(23)和(24)更新αk,否则停止运算;其中,上标i表示第i个方位向采样点;
T3、通过公式(25)和(26)计算(φk+1)i,k=k+1,返回步骤T2;
公式(23)-(26)中的下标k表示对于第i个方位向采样点的第k次迭代求解;通过上述步骤,实现对相位误差进行精确估计并实现良好成像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710093572.2A CN106802416B (zh) | 2017-02-21 | 2017-02-21 | 一种快速因式分解后向投影sar自聚焦方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710093572.2A CN106802416B (zh) | 2017-02-21 | 2017-02-21 | 一种快速因式分解后向投影sar自聚焦方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106802416A true CN106802416A (zh) | 2017-06-06 |
CN106802416B CN106802416B (zh) | 2020-04-07 |
Family
ID=58987517
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710093572.2A Active CN106802416B (zh) | 2017-02-21 | 2017-02-21 | 一种快速因式分解后向投影sar自聚焦方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106802416B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108205135A (zh) * | 2018-01-22 | 2018-06-26 | 西安电子科技大学 | 基于无插值融合快速后向投影的雷达视频成像方法 |
CN108614249A (zh) * | 2018-04-12 | 2018-10-02 | 北京航空航天大学 | 相位误差的估计方法、装置、补偿方法和系统 |
CN110095775A (zh) * | 2019-04-29 | 2019-08-06 | 西安电子科技大学 | 基于混合坐标系的颠簸平台sar快速时域成像方法 |
CN110554385A (zh) * | 2019-07-02 | 2019-12-10 | 中国航空工业集团公司雷华电子技术研究所 | 机动轨迹合成孔径雷达自聚焦成像方法、装置及雷达系统 |
CN111352108A (zh) * | 2020-02-28 | 2020-06-30 | 南昌大学 | 基于ffbp逆向处理的sar回波信号快速仿真方法 |
CN111537999A (zh) * | 2020-03-04 | 2020-08-14 | 云南电网有限责任公司电力科学研究院 | 一种稳健高效的分解投影自动聚焦方法 |
CN112946649A (zh) * | 2021-04-08 | 2021-06-11 | 电子科技大学 | 一种适用于任意子孔径长度的pfa成像方法 |
CN113640795A (zh) * | 2021-07-27 | 2021-11-12 | 北京理工大学 | 一种基于生成对抗网络的sar智能参数化自聚焦方法 |
CN113702974A (zh) * | 2021-09-01 | 2021-11-26 | 上海无线电设备研究所 | 一种机载/弹载合成孔径雷达图像快速优化方法 |
CN114578355A (zh) * | 2022-03-03 | 2022-06-03 | 西安电子科技大学 | 高超声速飞行器合成孔径雷达快速时域成像方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102854507A (zh) * | 2012-09-12 | 2013-01-02 | 电子科技大学 | 一种基于gpu后向投影双站合成孔径雷达成像方法 |
CN103728617A (zh) * | 2014-01-13 | 2014-04-16 | 电子科技大学 | 双基地合成孔径雷达时域快速成像方法 |
CN103869315A (zh) * | 2014-03-18 | 2014-06-18 | 电子科技大学 | 临近空间圆周合成孔径雷达快速后向投影成像方法 |
CN105093223A (zh) * | 2015-08-10 | 2015-11-25 | 中国人民解放军国防科学技术大学 | 双站前视sar的快速时域成像方法 |
-
2017
- 2017-02-21 CN CN201710093572.2A patent/CN106802416B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102854507A (zh) * | 2012-09-12 | 2013-01-02 | 电子科技大学 | 一种基于gpu后向投影双站合成孔径雷达成像方法 |
CN103728617A (zh) * | 2014-01-13 | 2014-04-16 | 电子科技大学 | 双基地合成孔径雷达时域快速成像方法 |
CN103869315A (zh) * | 2014-03-18 | 2014-06-18 | 电子科技大学 | 临近空间圆周合成孔径雷达快速后向投影成像方法 |
CN105093223A (zh) * | 2015-08-10 | 2015-11-25 | 中国人民解放军国防科学技术大学 | 双站前视sar的快速时域成像方法 |
Non-Patent Citations (1)
Title |
---|
李浩林 等: "一种适用于快速分解后向投影聚束SAR成像的自聚焦方法", 《航空学报》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108205135A (zh) * | 2018-01-22 | 2018-06-26 | 西安电子科技大学 | 基于无插值融合快速后向投影的雷达视频成像方法 |
CN108205135B (zh) * | 2018-01-22 | 2022-03-04 | 西安电子科技大学 | 基于无插值融合快速后向投影的雷达视频成像方法 |
CN108614249A (zh) * | 2018-04-12 | 2018-10-02 | 北京航空航天大学 | 相位误差的估计方法、装置、补偿方法和系统 |
CN108614249B (zh) * | 2018-04-12 | 2021-06-11 | 北京航空航天大学 | 相位误差的估计方法、装置、补偿方法和系统 |
CN110095775A (zh) * | 2019-04-29 | 2019-08-06 | 西安电子科技大学 | 基于混合坐标系的颠簸平台sar快速时域成像方法 |
CN110095775B (zh) * | 2019-04-29 | 2023-03-14 | 西安电子科技大学 | 基于混合坐标系的颠簸平台sar快速时域成像方法 |
CN110554385B (zh) * | 2019-07-02 | 2022-10-28 | 中国航空工业集团公司雷华电子技术研究所 | 机动轨迹合成孔径雷达自聚焦成像方法、装置及雷达系统 |
CN110554385A (zh) * | 2019-07-02 | 2019-12-10 | 中国航空工业集团公司雷华电子技术研究所 | 机动轨迹合成孔径雷达自聚焦成像方法、装置及雷达系统 |
CN111352108A (zh) * | 2020-02-28 | 2020-06-30 | 南昌大学 | 基于ffbp逆向处理的sar回波信号快速仿真方法 |
CN111537999A (zh) * | 2020-03-04 | 2020-08-14 | 云南电网有限责任公司电力科学研究院 | 一种稳健高效的分解投影自动聚焦方法 |
CN111537999B (zh) * | 2020-03-04 | 2023-06-30 | 云南电网有限责任公司电力科学研究院 | 一种稳健高效的分解投影自动聚焦方法 |
CN112946649A (zh) * | 2021-04-08 | 2021-06-11 | 电子科技大学 | 一种适用于任意子孔径长度的pfa成像方法 |
CN113640795A (zh) * | 2021-07-27 | 2021-11-12 | 北京理工大学 | 一种基于生成对抗网络的sar智能参数化自聚焦方法 |
CN113640795B (zh) * | 2021-07-27 | 2024-02-13 | 北京理工大学 | 一种基于生成对抗网络的sar智能参数化自聚焦方法 |
CN113702974A (zh) * | 2021-09-01 | 2021-11-26 | 上海无线电设备研究所 | 一种机载/弹载合成孔径雷达图像快速优化方法 |
CN113702974B (zh) * | 2021-09-01 | 2023-09-19 | 上海无线电设备研究所 | 一种机载/弹载合成孔径雷达图像快速优化方法 |
CN114578355B (zh) * | 2022-03-03 | 2022-10-21 | 西安电子科技大学 | 高超声速飞行器合成孔径雷达快速时域成像方法 |
CN114578355A (zh) * | 2022-03-03 | 2022-06-03 | 西安电子科技大学 | 高超声速飞行器合成孔径雷达快速时域成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106802416B (zh) | 2020-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106802416A (zh) | 一种快速因式分解后向投影sar自聚焦方法 | |
CN104931967B (zh) | 一种改进的高分辨率sar成像自聚焦方法 | |
CN104950305B (zh) | 一种基于稀疏约束的实波束扫描雷达角超分辨成像方法 | |
CN105549049B (zh) | 一种应用于gps导航的自适应卡尔曼滤波算法 | |
CN106249237B (zh) | 一种曲线轨迹下大斜视sar频域成像方法 | |
CN104007440B (zh) | 一种加速分解后向投影聚束合成孔径雷达成像方法 | |
CN105842694B (zh) | 一种基于ffbp sar成像的自聚焦方法 | |
CN103869311B (zh) | 实波束扫描雷达超分辨成像方法 | |
CN103364772B (zh) | 基于实数域广义多重信号分类算法的目标低仰角估计方法 | |
CN103777178B (zh) | 一种同步误差补偿方法、设备及系统 | |
CN107271993A (zh) | 一种基于最大后验的扫描雷达角超分辨成像方法 | |
Wang et al. | Extending Loffeld's bistatic formula for the general bistatic SAR configuration | |
CN109752710B (zh) | 一种基于稀疏贝叶斯学习的快速目标角度估计方法 | |
CN104316923A (zh) | 针对合成孔径雷达bp成像的自聚焦方法 | |
CN104713560A (zh) | 基于期望最大化的多源测距传感器空间配准方法 | |
CN105785358B (zh) | 一种方向余弦坐标系下带多普勒量测的雷达目标跟踪方法 | |
CN106291543A (zh) | 一种运动平台扫描雷达超分辨成像方法 | |
CN104251990A (zh) | 合成孔径雷达自聚焦方法 | |
CN106199600A (zh) | 基于多普勒估计的方位多通道合成孔径雷达成像方法 | |
CN104730500B (zh) | 一种合成孔径雷达残余距离徙动校正方法 | |
CN105699969A (zh) | 基于广义高斯约束的最大后验估计角超分辨成像方法 | |
CN107402380A (zh) | 一种实现多普勒波束锐化成像的快速自适应迭代方法 | |
CN105954742B (zh) | 一种球坐标系下带多普勒观测的雷达目标跟踪方法 | |
CN105022040A (zh) | 基于杂波数据联合拟合的阵元误差估计方法 | |
CN104793177B (zh) | 基于最小二乘法的麦克风阵列测向方法 |
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 |