CN105223572B - 一种基于pfa算法的正前视双基sar成像处理方法 - Google Patents
一种基于pfa算法的正前视双基sar成像处理方法 Download PDFInfo
- Publication number
- CN105223572B CN105223572B CN201510582465.7A CN201510582465A CN105223572B CN 105223572 B CN105223572 B CN 105223572B CN 201510582465 A CN201510582465 A CN 201510582465A CN 105223572 B CN105223572 B CN 105223572B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msubsup
- mfrac
- signal
- 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.)
- Expired - Fee Related
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/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
-
- 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
- G01S13/9043—Forward-looking SAR
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)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种基于PFA算法的正前视双基SAR成像处理方法,包括读入正前视双基SAR模式回波仿真数据及相应的成像参数、对回波仿真数据在场景中心点进行时域的完全距离徙动补偿、将解线频调后的回波仿真数据在方位向进行去斜处理、将去斜处理后的回波信号进行坐标转换并离散化、将离散化后的信号利用sinc插值依次进行距离向插值和方位向插值、将经过二维插值解耦后的信号进行二维Fourier变换,得到聚焦后的信号。本发明通过对场景中心点进行二维去斜处理,并在处理过程中近似认为微波传输满足远场条件,即波前近似为平面,可以对系统的运动误差具有较好的抵制能力。
Description
技术领域
本发明属于信号处理领域,特别涉及一种基于PFA(Polar Format Algorithm)算法的正前视双基SAR(Synthetic Aperture Radar)成像处理方法。
背景技术
合成孔径雷达(SAR)是一种先进的天基信息获取的有效载荷,是获取信息的一种重要手段。相比光学传感器,合成孔径雷达(SAR)能克服云雾雨雪和夜暗条件的限制对地面目标成像,实现全天时、全天候、高分辨率、大幅面对地观测,在海洋观测、精细农业、环境保护、灾害监测、资源勘查、地质测绘、政府公共决策等方面有着广泛的应用。许多情况下,能够起到其它空间遥感手段无法替代的作用。
然而,传统的单基地机载SAR用于军事侦察时,几乎没有隐蔽性可言,同时,单基地SAR在正前视的情况下,由于多普勒调频率为零,导致目标在整个合成孔径时间内的带宽为零,从而使得单基地SAR无法对正前方的目标进行成像。双基SAR作为多平台SAR中的一种,几何构型非常灵活,由于收发平台的分置,它的距离分辨率和多普勒分辨率的大小和方向由收发双方共同决定,因此可以通过合理设置双基构型,使得某一平台前视时,距离分辨率方向和多普勒分辨率方向仍存在夹角(即仍能构成一个二维分辨单元),从而可实现前视二维成像。
但是,正前视双基SAR的成像处理与传统的单基SAR模式有着很多的不同,由于双基SAR回波距离历程由收发天线与目标的相对运动共同决定,其距离历程较单站SAR复杂,单基地SAR的经典频域算法或混合域算法均无法适用于双基SAR的成像,同时,正前视双基SAR回波还具有二维空变性,给成像处理带来了更大的难度,因此,如何合理、高效、精确的对正前视双基SAR的回波数据进行成像处理,显得尤为重要。
发明内容
本发明的目的是为了解决上述问题,针对正前视双基SAR回波数据量具有较强的二维空变性等特点,基于传统的单基地SAR的PFA成像算法,提出了一种正前视双基SAR成像处理方法。
一种基于PFA算法的正前视双基SAR成像处理方法,包括以下步骤:
步骤一:读入正前视双基SAR模式回波仿真数据S(τ,t;RS)以及相应的成像参数,其中,τ为快时间、t为慢时间、RS为发射平台和接收平台到目标的距离和,具体包括:方位向采样点数Na,距离向采样点数Nr,信号采样率fc,信号带宽Bw,脉冲宽度Tp,脉冲重复频率PRF,零时刻发射平台的位置T(xt,yt,zt)和接收平台的位置R(xr,yr,zr),参考点P的位置(x0,y0,0),发射平台速度接收平台速度信号波长λ,光速c,合成孔径时间Tsyn;
步骤二:对回波仿真数据S(τ,t;r)在场景中心点进行时域的完全距离徙动补偿,即解线频调(Dechirp)处理,得到解线频调后的信号S'(τ,t;r);
步骤三:将解线频调后的回波仿真数据S'(τ,t;r)在方位向进行去斜处理(RVP),得到去斜处理后的信号Sr1(τ,t;r);
步骤四:设定正前视双基SAR的成像坐标系,将去斜处理后的回波信号Sr1(τ,t;r)进行坐标转换并离散化,得到离散后的信号Sr2(i,m),其中i表示距离向采样点序列,m表示方位向采样点序列;
步骤五:将离散化后的信号Sr2(i,m)利用sinc插值依次进行距离向插值和方位向插值,得到二维插值后的信号Sr3(i,m);
步骤六:将经过二维插值解耦后的信号Sr3(i,m)进行二维Fourier变换,得到聚焦后的信号Sr4(i',m')。
本发明的优点在于:
(1)本发明通过对场景中心点进行二维去斜处理,并在处理过程中近似认为微波传输满足远场条件,即波前近似为平面,可以对系统的运动误差具有较好的抵制能力。
附图说明
图1是本发明的方法流程图。
图2是本发明所设定的成像坐标系。
图3是本发明步骤五中距离向插值效果示意图。
图4是本发明步骤五中方位向插值效果示意图。
图5是仿真场景设计示意图。
图6是正前视双基SAR点目标成像处理结果。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明是一种基于PFA算法的正前视双基SAR成像处理方法,具体流程如图1所示,包括以下步骤:
步骤一:读入正前视双基SAR模式回波仿真数据S(τ,t;RS)以及相应的成像参数,其中,τ为快时间、t为慢时间、RS为发射平台和接收平台到目标的距离和,具体包括:方位向采样点数Na,距离向采样点数Nr,信号采样率fc,信号带宽Bw,脉冲宽度Tp,脉冲重复频率PRF,零时刻发射平台的位置T(xt,yt,zt)和接收平台的位置R(xr,yr,zr),其中,xt、yt和zt为发射机在直角坐标系中的位置,xr、yr和zr为接收机在直角坐标系中的位置。参考点P的位置(x0,y0,0),其中x0,y0为参考点在直角坐标系下的位置。发射平台速度接收平台速度信号波长λ,光速c;
步骤二:对回波仿真数据S(τ,t;r)在场景中心点进行时域的完全距离徙动补偿,即解线频调(Dechirp)处理,得到解线频调后的信号S'(τ,t;r);
具体为:将回波仿真数据S(τ,t;r)乘以补偿滤波函数H1(τ,t),得到解线频调后的信号S'(τ,t;r),补偿滤波函数H1(τ,t)为:
其中:Tsyn是合成孔径时间,tc为合成孔径中心时刻,γ为距离向调频率。
步骤三:将解线频调后的回波仿真数据S'(τ,t;r)在方位向进行去斜处理,得到去斜处理后的信号Sr1(τ,t;r);
具体为:
(a)对回波仿真数据S'(τ,t;r)进行去斜处理,去斜处理后的信号在方位向进行傅里叶变换,得到距离多普勒域信号;
(b)距离多普勒域信号再乘以Dechirp参考函数H2(fτ,t),得到方位向去斜处理后的信号S′i(τ,t;r),Dechirp参考函数H2(fτ,t)为:
其中:fτ为快时间对应的频率,γ为距离向调频率;
(c)方位向去斜处理后的信号S′i(τ,t;r)在方位向进行傅里叶变换,得到去斜处理后的时域信号Sr1(τ,t;r)。
步骤四:设定正前视双基SAR的成像坐标系,将去斜处理后的回波信号Sr1(τ,t;r)进行坐标转换并离散化,得到离散后的信号Sr2(i,m),其中i表示距离向采样点序列,m表示方位向采样点序列;
具体为:
(a)设定正前视双基SAR的成像坐标系,场景中心为原点,YS轴指向距离向方向,ZS轴指向上方,XS轴为原点与YS轴、ZS轴以右手定则确定。图2所示为直角坐标系和所设定的成像坐标系之间的转换示意图,其中,X,Y,Z为直角坐标系下的三轴,Ω为YS轴与X轴夹角,为Ys轴与中心时刻发射机与场景中心斜距在地面的投影方向夹角,为Ys轴与中心时刻接收机与场景中心斜距在地面的投影方向夹角。
(b)将去斜处理后的回波信号Sr1(τ,t;r)进行坐标转换并离散化,得到离散后的信号Sr2(i,m)。
步骤五:将离散化后的信号Sr2(i,m)利用sinc插值依次进行距离向插值和方位向插值,得到二维插值后的信号Sr3(i,m);
具体为:
(a)对离散化后的信号Sr2(i,m)进行距离向插值,得到经过距离向插值后的信号,图3表示距离向插值后的效果。其中,目标插值序列的选取如下:
KYsref(i)=max(KYs(1,m))+ΔKYsrefi,(i=1,...,Nr,,m=1,...,Na) (3)
式中:
式中:Ω为YS轴与X轴夹角,i为离散化后距离向采样点序列i,m为离散化后方位向采样点序列用。
(b)对经过距离向插值后的信号进行方位向插值,得到经过方位向插值后的信号,即二维插值后的信号Sr3(i,m),图4表示方位向插值后的效果。其中,目标插值序列的选取如下:
KXsref(m)=KXsref(1)+ΔKXsrefm,(m=1,...,Na) (4)
式中:
步骤六:将经过二维插值解耦后的信号Sr3(i,m)进行二维Fourier变换,得到聚焦后的信号Sr4(i',m')。
实施例:
本发明的一种基于PFA算法的正前视双基SAR成像处理方法,具体为:
步骤一:读入正前视双基SAR模式回波仿真数据S(τ,t;RS)以及相应的成像参数,包括:方位向采样点数Na=4096,距离向采样点数Nr=4096,信号采样率fc=15GHz,信号带宽Bw=180MHz,脉冲宽度τ=3.5us,脉冲重复频率PRF=500Hz,零时刻发射平台的位置T=(-40000,0,6000)和接收平台的位置R=(0,-50000,6000),参考点P的位置(0,0,0),发射平台速度接收平台速度信号波长λ=0.02m,光速c=3×108m/s,合成孔径时间Tsyn=6s;
步骤二:对回波仿真数据S(τ,t;r)在场景中心点进行时域的完全距离徙动补偿,即解线频调(Dechirp)处理,得到解线频调后的信号S'(τ,t;r);方法如公式(1)所示;
步骤三:将解线频调后的回波仿真数据S'(τ,t;r)在方位向进行去斜处理,得到去斜处理后的信号Sr1(τ,t;r);
(a)对解线频调后的信号在方位向进行傅里叶变换;
(b)再乘以Dechirp参考函数;
方法如公式(4)所示;
(c)对去斜处理hoist的信号在方位向进行傅里叶逆变换;
步骤四:设定正前视双基SAR的成像坐标系,将去斜处理后的回波信号Sr1(τ,t;r)进行坐标转换并离散化,得到离散后的信号Sr2(i,m),其中i表示距离向采样点序列,m表示方位向采样点序列;:
(a)依据图2所设的成像坐标系统对去斜后的信号进行坐标转换;
(b)将转换后的信号进行离散化处理,得到离散后的信号
步骤五:将离散化后的信号Sr2(i,m)利用sinc插值依次进行距离向插值和方位向插值,得到二维插值后的信号Sr3(i,m);
(a)对离散化后的信号进行距离向插值;
目标插值序列的选取如公式(3)所示;
(b)对距离向插值后的信号进行方位向插值;
目标插值序列的选取如公式(4)所示;
步骤六:将经过二维插值解耦后的信号Sr3(i,m)进行二维Fourier变换,得到聚焦后的信号Sr4(i',m')。
为了说明该方法的有效性,进行如下仿真实验,场景设置为3×3点阵如图5所示,设置点目标O为场景中心点。成像参数如表1所示。
表1成像参数
利用本发明的方法进行成像处理并对成像质量进行评估。成像结果如图6所示,其中垂直方向为方位向,水平方向为距离向。表2详细的给出了各点目标在距离向的峰值旁瓣比和积分旁瓣比以及距离向和方位向的分辨率信息。评估结果表明各点目标的峰值旁瓣比和积分旁瓣比都和理论值较为接近,从而验证可本发明的有效性和准确性。
表2点目标距离向成像质量评估结果
本发明提供了一种基于PFA算法的正前视双基SAR成像处理方法。该方法主要用于接收机正前视的双基SAR模式,针对其SAR回波数据具有二维空变性的特点,提出了一种正前视双基SAR成像处理方法。
Claims (3)
1.一种基于PFA算法的正前视双基SAR成像处理方法,包括以下步骤:
步骤一:读入正前视双基SAR模式回波仿真数据S(τ,t;RS)以及相应的成像参数,其中,τ为快时间、t为慢时间、RS为发射平台和接收平台到目标的距离和,具体包括:方位向采样点数Na,距离向采样点数Nr,信号采样率fc,信号带宽Bw,脉冲宽度Tp,脉冲重复频率PRF,零时刻发射平台的位置T(xt,yt,zt)和接收平台的位置R(xr,yr,zr),其中,xt、yt和zt为发射机在直角坐标系中的位置,xr、yr和zr为接收机在直角坐标系中的位置;参考点P的位置(x0,y0,0),其中x0,y0为参考点在直角坐标系下的位置;发射平台速度接收平台速度信号波长λ,光速c;
步骤二:对回波仿真数据S(τ,t;r)在场景中心点进行时域的完全距离徙动补偿,即解线频调Dechirp处理,得到解线频调后的信号S'(τ,t;r);
将回波仿真数据S(τ,t;r)乘以补偿滤波函数H1(τ,t),得到解线频调后的信号S'(τ,t;r),补偿滤波函数H1(τ,t)为:
<mrow>
<msub>
<mi>H</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>r</mi>
<mi>e</mi>
<mi>c</mi>
<mi>t</mi>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
</mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>s</mi>
<mi>y</mi>
<mi>n</mi>
</mrow>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<mi>exp</mi>
<mo>{</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mi>c</mi>
</msub>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
</mrow>
<mi>c</mi>
</mfrac>
<mo>}</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mi>&pi;</mi>
<mi>&gamma;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
<mi>c</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:Tsyn是合成孔径时间,tc为合成孔径中心时刻,γ为距离向调频率;
步骤三:将解线频调后的回波仿真数据S'(τ,t;r)在方位向进行去斜处理,得到去斜处理后的信号Sr1(τ,t;r);
(a)对回波仿真数据S'(τ,t;r)进行去斜处理,去斜处理后的信号在方位向进行傅里叶变换,得到距离多普勒域信号;
(b)距离多普勒域信号再乘以Dechirp参考函数H2(fτ,t),得到方位向去斜处理后的信号Si'(τ,t;r),Dechirp参考函数H2(fτ,t)为:
<mrow>
<msub>
<mi>H</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&tau;</mi>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mi>&pi;</mi>
<mfrac>
<msubsup>
<mi>f</mi>
<mi>&tau;</mi>
<mn>2</mn>
</msubsup>
<mi>&gamma;</mi>
</mfrac>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:fτ为快时间对应的频率,γ为距离向调频率;
(c)方位向去斜处理后的信号Si'(τ,t;r)在方位向进行傅里叶变换,得到去斜处理后的时域信号Sr1(τ,t;r);
步骤四:设定正前视双基SAR的成像坐标系,将去斜处理后的回波信号Sr1(τ,t;r)进行坐标转换并离散化,得到离散后的信号Sr2(i,m),其中i表示距离向采样点序列,m表示方位向采样点序列;
步骤五:将离散化后的信号Sr2(i,m)利用sinc插值依次进行距离向插值和方位向插值,得到二维插值后的信号Sr3(i,m);
步骤六:将经过二维插值解耦后的信号Sr3(i,m)进行二维Fourier变换,得到聚焦后的信号Sr4(i',m')。
2.根据权利要求1所述的一种基于PFA算法的正前视双基SAR成像处理方法,所述的步骤四具体包括:
(a)设定正前视双基SAR的成像坐标系,场景中心为原点,YS轴指向距离向方向,ZS轴指向上方,XS轴为原点与YS轴、ZS轴以右手定则确定;
(b)将去斜处理后的回波信号Sr1(τ,t;r)进行坐标转换并离散化,得到离散后的信号Sr2(i,m)。
3.根据权利要求1所述的一种基于PFA算法的正前视双基SAR成像处理方法,所述的步骤五具体包括:
具体为:
(a)对离散化后的信号Sr2(i,m)进行距离向插值,得到经过距离向插值后的信号,目标插值序列的选取如下:
KYsref(i)=max(KYs(1,m))+ΔKYsrefi (3)
式中:i=1,...,Nr;m=1,...,Na;
<mrow>
<msub>
<mi>&Delta;K</mi>
<mrow>
<mi>Y</mi>
<mi>s</mi>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>Y</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>N</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>Y</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>r</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>Y</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mo>-</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&Omega;</mi>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>T</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<msubsup>
<mi>x</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>y</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>z</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<msubsup>
<mi>x</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>y</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>z</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&Omega;</mi>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mi>T</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<msubsup>
<mi>x</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>y</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>z</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<msubsup>
<mi>x</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>y</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>z</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中:Ω为YS轴与X轴夹角,i为离散化后距离向采样点序列i,m为离散化后方位向采样点序列用;
(b)对经过距离向插值后的信号进行方位向插值,得到经过方位向插值后的信号,即二维插值后的信号Sr3(i,m),目标插值序列的选取如下:
KXsref(m)=KXsref(1)+ΔKXsrefm (4)
式中:m=1,...,Na;
<mrow>
<msub>
<mi>&Delta;K</mi>
<mrow>
<mi>X</mi>
<mi>s</mi>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>Y</mi>
<mi>s</mi>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>X</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<msub>
<mi>N</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>Y</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<msub>
<mi>N</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>Y</mi>
<mi>s</mi>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>X</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>Y</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>X</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&Omega;</mi>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>T</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<msubsup>
<mi>x</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>y</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>z</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<msubsup>
<mi>x</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>y</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>z</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>+</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&Omega;</mi>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mi>T</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<msubsup>
<mi>x</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>y</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>z</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<msubsup>
<mi>x</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>y</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>z</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
2
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510582465.7A CN105223572B (zh) | 2015-09-14 | 2015-09-14 | 一种基于pfa算法的正前视双基sar成像处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510582465.7A CN105223572B (zh) | 2015-09-14 | 2015-09-14 | 一种基于pfa算法的正前视双基sar成像处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105223572A CN105223572A (zh) | 2016-01-06 |
CN105223572B true CN105223572B (zh) | 2017-08-29 |
Family
ID=54992629
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510582465.7A Expired - Fee Related CN105223572B (zh) | 2015-09-14 | 2015-09-14 | 一种基于pfa算法的正前视双基sar成像处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105223572B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107229050B (zh) * | 2017-05-11 | 2020-04-21 | 西北工业大学 | 一种基于极坐标格式的雷达成像优化方法 |
CN108318880B (zh) * | 2018-01-26 | 2021-10-08 | 西安电子科技大学 | 一种运动误差参数化的极坐标sar成像方法 |
CN109143237B (zh) * | 2018-08-31 | 2020-10-09 | 电子科技大学 | 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法 |
CN109358328B (zh) * | 2018-11-06 | 2020-06-30 | 电子科技大学 | 机动平台双基地前视sar的极坐标格式成像方法 |
CN112346057B (zh) * | 2019-08-06 | 2023-07-14 | 北京华航无线电测量研究所 | 一种图像处理过程中的二维插值模块 |
CN110515080B (zh) * | 2019-08-30 | 2022-12-02 | 西安电子科技大学 | 基于半径角度插值的高超声速机动平台sar成像方法 |
CN111158056B (zh) * | 2019-12-26 | 2022-05-24 | 北京遥测技术研究所 | 一种基于稀疏阵列的安检装置及方法 |
CN113341412B (zh) * | 2021-05-26 | 2024-06-11 | 南京航空航天大学 | 一种基于改进Sinc插值技术的PFA成像方法 |
CN116148856B (zh) * | 2023-04-17 | 2023-07-04 | 南京邮电大学 | 一种sar动目标二维自聚焦成像处理方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7511656B2 (en) * | 2006-02-10 | 2009-03-31 | Raytheon Company | Spotlight synthetic aperture radar (SAR) system and method for generating a SAR map in real-time using a modified polar format algorithm (PFA) approach |
CN103969628B (zh) * | 2014-04-08 | 2017-01-04 | 南京航空航天大学 | 一种基于压缩感知的合成孔径雷达pfa成像算法 |
CN104391297B (zh) * | 2014-11-17 | 2017-09-15 | 南京航空航天大学 | 一种划分子孔径pfa雷达成像方法 |
-
2015
- 2015-09-14 CN CN201510582465.7A patent/CN105223572B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105223572A (zh) | 2016-01-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105223572B (zh) | 一种基于pfa算法的正前视双基sar成像处理方法 | |
CN103869311B (zh) | 实波束扫描雷达超分辨成像方法 | |
CN109143237B (zh) | 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法 | |
CN107238824B (zh) | 基于先验dem数据的星载sar图像几何精校正方法 | |
CN105259552A (zh) | 一种基于非线性调频信号的合成孔径雷达成像方法和装置 | |
CN103487803B (zh) | 迭代压缩模式下机载扫描雷达成像方法 | |
CN104833972B (zh) | 一种双基地调频连续波合成孔径雷达频率变标成像方法 | |
CN102749621B (zh) | 一种双基地合成孔径雷达频域成像方法 | |
CN103983974B (zh) | 双站调频连续波合成孔径雷达成像方法 | |
CN105487074B (zh) | 一种双基地合成孔径雷达数值距离多普勒成像方法 | |
CN103630902B (zh) | 一种用于分时接收阵列sar的相位校正成像方法 | |
CN102004250B (zh) | 基于频域展开的星机联合双基地合成孔径雷达成像方法 | |
CN104122549B (zh) | 基于反卷积的雷达角超分辨成像方法 | |
CN103235309B (zh) | 临近空间慢速平台sar成像方法 | |
CN104330779B (zh) | 机载合成孔径雷达运动误差补偿方法 | |
CN105093224A (zh) | 一种大斜视合成孔径雷达成像处理方法 | |
CN103969628B (zh) | 一种基于压缩感知的合成孔径雷达pfa成像算法 | |
CN107544068A (zh) | 一种基于频域bp的图像域宽带合成方法 | |
CN103308913A (zh) | 一种高速飞行器载双天线前视sar解模糊算法 | |
CN106093932A (zh) | 一种扫描波束高分辨率的雷达散射计 | |
CN108710111A (zh) | 一种机载双基前视sar方位相位的二维空变校正方法 | |
CN103543451A (zh) | 一种基于压缩感知的多径虚像抑制sar后处理方法 | |
CN107153191B (zh) | 一种针对隐形飞机的双基isar成像检测方法 | |
CN104820221B (zh) | 多基合成孔径雷达的目标三维定位方法 | |
CN109143236B (zh) | 适用于复杂飞行轨迹的双基聚束sar大场景成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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: 20170829 Termination date: 20200914 |