CN104749570B - 一种移不变机载双基合成孔径雷达目标定位方法 - Google Patents

一种移不变机载双基合成孔径雷达目标定位方法 Download PDF

Info

Publication number
CN104749570B
CN104749570B CN201510151647.9A CN201510151647A CN104749570B CN 104749570 B CN104749570 B CN 104749570B CN 201510151647 A CN201510151647 A CN 201510151647A CN 104749570 B CN104749570 B CN 104749570B
Authority
CN
China
Prior art keywords
scene
distance
target
point
echo data
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
Application number
CN201510151647.9A
Other languages
English (en)
Other versions
CN104749570A (zh
Inventor
杨建宇
钟徐琦
武俊杰
黄钰林
杨海光
李中余
杨晓波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201510151647.9A priority Critical patent/CN104749570B/zh
Publication of CN104749570A publication Critical patent/CN104749570A/zh
Application granted granted Critical
Publication of CN104749570B publication Critical patent/CN104749570B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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/904—SAR modes

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)
  • Radar Systems Or Details Thereof (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明公开了一种移不变机载双基合成孔径雷达目标定位方法，包括以下步骤：S1、计算双基SAR的回波数据；S2、对回波数据进行距离向脉冲压缩；S3、对场景中心点进行多普勒质心估计；S4、场景中心点定位；S5、对回波数据进行聚焦；S6、场景中所有目标点定位；S7、目标场景地距投影。本发明实现了在双基SAR下对成像目标点的精确定位，同时利用邻近像素点的关系式，快速导出整个图像的地理位置，并将其投影到地理坐标系中，极大地减少了计算量；可以广泛应用于双基SAR图像几何校正、目标定位、运动目标探测等领域。

Description

一种移不变机载双基合成孔径雷达目标定位方法
技术领域
本发明属于雷达信号处理技术领域，特别涉及一种移不变机载双基合成孔径雷达目标定位方法。
背景技术
SAR(合成孔径雷达)是一种全天时、全天候的现代高分辨率微波遥感成像雷达，在军事侦察、地形测绘、植被分析、海洋及水文观测、环境及灾害监视、资源勘探以及地壳微变检测等领域，SAR发挥了越来越重要的作用。
双基SAR由于收发分置而有着很多突出的优点，它能获取目标的非后向散射信息，具有作用距离远、隐蔽性和抗干扰性强等特点。另外，由于双基SAR接收机不含大功率器件，其功耗低、体积小、重量轻，便于多种类型的飞机携带，造价较低。总之，双基SAR作为一种空间对地观测的新手段，在民用和军用领域都有着广阔的发展空间。
在双基SAR中，不同距离单元的目标与收发站的相对位置关系不同，因而导致不同距离单元的目标有不同的多普勒频移；由于两个根号的影响，相比单基SAR目标定位，双基SAR目标定位所需要求解的多普勒质心公式与双基距离和公式都要更加复杂，且无法求解出精确解。
目前国际上针对双基SAR目标定位的研究，公开文献有：An Improved TargetLocation Algorithm for General Airborne Bistatic SAR System,Xin Zhou,Yong Liand Ya Li,Publication Year:2013,Page(s):1-4.中，给出了一种双基SAR目标定位方法，但它只考虑对特定目标进行定位，而忽略了SAR图像中不同目标间的位置关系。因此，对场景中所有目标点的定位都要进行同样复杂的计算过程，效率极低。
相关文献：S.Wenfeng,C.An,and Z.Changyao,“Range-doppler approach forcalibration and location of air-borne SAR image,”in Radar,2006.CIE’06.International Conference on,Oct 2006,pp.1-4.中，针对单基SAR的配置和特点，设计了一种利用回波信号的距离和多普勒质心信息的目标定位方法。而双基SAR与单基SAR在几何配置上有很大不同，因此该方法无法应用于双基SAR目标定位中。
发明内容
本发明的目的在于克服现有技术的不足，提供一种实现了在双基SAR下对成像目标点的精确定位，可以应用于双基SAR图像几何校正、目标定位、运动目标探测等领域的移不变机载双基合成孔径雷达目标定位方法。
本发明的目的是通过以下技术方案来实现的：一种移不变机载双基合成孔径雷达目标定位方法，包括以下步骤：
S1、计算双基SAR的回波数据；
S2、对回波数据进行距离向脉冲压缩；
S3、对场景中心点进行多普勒质心估计：选择回波数据条带的距离向中心距离门为场景中心所在距离门，取出该距离门所在列并对其进行多普勒质心估计；
S4、场景中心点定位：根据场景中心点所在距离单元建立双基距离和方程，再根据S3中估计出的多普勒质心建立多普勒频率方程，然后求解场景中心点位置；
S5、对回波数据进行聚焦：利用所求得到场景中心点在直角坐标系中的位置与接收站、发射站的相对位置对回波数据进行校正并聚焦；
S6、场景中所有目标点定位：利用场景中心点和它在图像中的像素位置推导邻近像素点在所建立的三维地理坐标系下的坐标；
S7、目标场景地距投影：利用S6中求得的像素与坐标位置对应关系，将图像投影到地理坐标系中，得到地理坐标系下的真实图像。
进一步地，所述的步骤S1具体实现方法为：令场景中心点被波束中心照射时刻，发射平台固定，位置记为(xT,yT,zT)，其中，xT、yT和hT分别为发射站的x轴、y轴和z轴坐标；接收站位置记为(0,0,hR)，其中，0、0和hR分别为接收站的x轴、y轴和z轴坐标；接收站与发射站速度记为V，并沿y轴运动；建立了以接收站正下方为原点、高度为Z轴、速度方向为Y轴的三维坐标系：
将方位时间向量记为：
其中，PRI为脉冲重复间隔，Na为目标回波方位点数；
双基距离历史和为：Rb(t；x,y)＝RT(t；x,y)+RR(t；x,y)
其中，t为方位时间，RT(t；x,y)和RR(t；x,y)分别为发射站和接收站的距离历史：
从而得到回波数据的表达式为：
其中，A0是散射系数的幅度，ωr(·)为距离向包络，ωa(·)方位向包络，τ是快时间变量，fc是载波频率，c是光速，Kr是距离向调频率，Ta是合成孔径时间，t0是目标点(x,y)的波束中心穿越时刻。
进一步地，所述的步骤S2具体实现方法为：利用发射的Chirp信号作为参考函数对回波数据进行距离向脉压，Chirp信号的表达式为：
S(τ)＝A0wr(τ)exp(jπKrτ2) (4)
取其反向共轭，建立距离向脉压信号的表达式：
S*(-τ)＝A0wr(-τ)exp(-jπKr(-τ)2) (5)
将回波数据的距离向数据与式(5)分别进行FFT后，在频域上相乘，然后进行IFFT得到距离向脉压后的回波数据。
进一步地，所述的步骤S3具体实现方法为：从回波数据表达式能够得出，多普勒质心fdc就是方位向频谱的中心，且方位向频谱相对于多普勒质心fdc为中心左右对称；
以估计的接收天线斜视角所计算出来的场景中心点位置为初始位置，计算该初始位置的多普勒质心并利用为迭代初值，采用以下步骤对场景中心点的多普勒质心fdc进行精确估计：
S3.1：分别计算小于的宽度频谱能量大于的宽度频谱能量
S3.2：比较与的大小，如果小于则转到步骤S3.3，否则转到步骤S3.6；
S3.3：按多普勒频谱分辨率增大即
S3.4：分别计算小于的宽度频谱能量大于的宽度频谱能量
S3.5：比较与的大小，如果小于则i＝i+1并转到步骤S3.3，否则输出结束；
S3.6：按多普勒频谱分辨率减小即
S3.7：分别计算小于的宽度频谱能量大于的宽度频谱能量
S3.8：比较与的大小，如果大于则i＝i+1并转到步骤S3.6，否则输出结束；
其中，i≥1，Δfdc为多普勒频谱分辨率，fa为方位向采样频率，Na为方位向采样点数。
进一步地，所述的步骤S4具体实现方法为：
已知当t＝0时，双基距离和方程为
其中，(xdc,ydc)为场景中心点的坐标；
多普勒质心方程为：
双基距离和已知为Rb，根据式(6)推导出：
其中，a＝R2-(yT-yR)2，b＝2yT(yT-yR)+2(xT-xR)xdc+k，
由于两个根号的影响，式(7)与式(8)无法解出精确解，因此采用搜索的方式来估计场景中心点的位置，具体步骤如下：
S4.1：由接收机与发射机的位置和目标点的双基距离和(式6)，得到目标点所在椭圆曲线方程为y＝fellipse(x)；
S4.2：由测量的天线斜视角计算出目标点初始X轴坐标计算出初始Y轴坐标
S4.3：根据计算出多普勒质心
S4.4：令椭圆曲线上左侧点的X轴坐标为Y轴坐标为ydc_left＝fellipse(xdc_left)，代入式(7)中，计算其多普勒质心为fdc_left；
S4.5：令椭圆曲线上右侧点的X轴坐标为Y轴坐标为ydc_right＝fellipse(xdc_right)，代入式(7)中，计算出多普勒质心为fdc_right；
S4.6：比较|fdc-fdc_left|与|fdc-fdc_right|的大小，如果|fdc-fdc_left|<|fdc-fdc_right|，则转到步骤S4.7；否则转到S4.10；
S4.7：令xdc＝xdc_left，fD＝fdc_left；
S4.8：令xleft＝xdc-Δx，则左侧点的位置坐标为(xleft,fellipse(xleft))，计算其多普勒质心，记为f'dc_left；
S4.9：比较|fdc-f'dc_left|与|fdc-fD|的大小，如果|fdc-f'dc_left|<＝|fdc-fD|，则转到S4.7；否则转到S4.13；
S4.10：令xdc＝xdc_right，fD＝fdc_right；
S4.11：令xright＝xdc+Δx，则右侧点的位置坐标为(xright,fellipse(xright))，计算其多普勒质心，记为f'dc_right；
S4.12：比较|fdc-f'dc_right|与|fdc-fD|的大小，如果|fdc-f'dc_right|<＝|fdc-fD|，则转到S4.10；否则转到S4.13；
S4.13：得到场景中心点坐标为(xdc,ydc)，结束；
其中，Δx为根据精度要求设定出来的搜索步长。
具体地，所述的步骤S5具体实现方法为：
S5.1、把步骤S4中得到场景中心点位置代入距离徒动公式中，得到回波数据中距离徒动校正量与多普勒频率的关系函数为：
式(9)中的未知系数由以下方程组求得：
其中，
S5.2、将回波数据进行方位向FFT变换，转换到距离多普勒域；
S5.3、在距离多普勒域进行距离徒动校正。
进一步地，所述的步骤S6具体实现方法为：通过场景中心点的地理坐标和它图像中的像素位置，推导出邻近像素点对应的坐标位置，具体推导方法包括以下步骤：
S6.1、已知某个像素点(i,j)对应的坐标为(xi,j,yi,j)，其中(i,j)分别为图像中距离向与方位向位置，它距离向的相邻点分别为(xi-1,j,yi-1,j)和(xi+1,j,yi+1,j)，首先，由波束模型可知：
对于近点，可得
其中，θc为接收天线斜视角；
最终得出：
其中，
S6.2、采用与步骤S61相同方法计算出(xi+1,j,yi+1,j)；
S6.3、它方位向上相邻点分别为(xi,j-1,yi,j-1)、(xi,j+1,yi,j+1)；
对于点(xi,j+1,yi,j+1)来说，发射机与接收机的新位置分别为(0,v0*PRI,HR)、(xT,yT+vy*PRI,HT)，已知关系式为：
其中θc为接收天线斜视角；
最终得出：
其中，
S6.4、采用与步骤S63相同方法计算出(xi,j-1,yi,j-1)；
S6.5、根据步骤S61～S64的计算方法得到整个图像中的每个像素所对应的地理坐标位置。
进一步地，所述的步骤S7中图像投影到地理坐标系中的具体投影关系为：
本发明的有益效果是：在单基SAR定位的基础上，利用估计场景中心点的双基距离和与多普勒质心信息，并通过所设计的算法来计算出场景中心的精确位置坐标；同时，又导出场景中心点与邻近像素点的位置关系，利用位置关系快速计算出图像上所有像素点所对应的地理坐标位置；实现了在双基SAR下对成像目标点的精确定位，同时利用邻近像素点的关系式，快速导出整个图像的地理位置，并将其投影到地理坐标系中，极大地减少了计算量；可以广泛应用于双基SAR图像几何校正、目标定位、运动目标探测等领域。
附图说明
图1为本发明的目标定位方法流程图；
图2是本发明具体实施例中采用的移不变模式双基SAR系统结构图；
图3为本发明具体实施例中采用的目标场景布置图；
图4为本发明具体实施例中回波距离向脉压后的二维时域图；
图5为本发明的多普勒质心估计方法流程图；
图6是本发明的场景中心点定位方法流程图；
图7是本发明具体实施例中对回波数据进行聚焦后得到的图像；
图8是本发明具体实施例中对图3中9个点目标进行定位投影后的结果图。
具体实施方式
下面结合如图和具体实施例对本发明的技术方案进行进一步描述。
如图1所示，一种移不变机载双基合成孔径雷达目标定位方法，包括以下步骤：
S1、计算双基SAR的回波数据：令场景中心点被波束中心照射时刻，发射平台固定，位置记为(xT,yT,zT)，其中，xT、yT和hT分别为发射站的x轴、y轴和z轴坐标；接收站位置记为(0,0,hR)，其中，0、0和hR分别为接收站的x轴、y轴和z轴坐标；接收站与发射站速度记为V，并沿y轴运动；建立了以接收站正下方为原点、高度为Z轴、速度方向为Y轴的三维坐标系：
将方位时间向量记为：
其中，PRI为脉冲重复间隔，Na为目标回波方位点数；
双基距离历史和为：Rb(t；x,y)＝RT(t；x,y)+RR(t；x,y)
其中，t为方位时间，RT(t；x,y)和RR(t；x,y)分别为发射站和接收站的距离历史：
从而得到回波数据的表达式为：
其中，A0是散射系数的幅度，ωr(·)为距离向包络，ωa(·)方位向包络，τ是快时间变量，fc是载波频率，c是光速，Kr是距离向调频率，Ta是合成孔径时间，t0是目标点(x,y)的波束中心穿越时刻。
本实施例的移不变模式双基SAR系统结构图如图2所示，本实施例的接收机坐标为(0,0,4)km、速度为(0,200,0)m/s，发射机坐标为(-12,2,6)km、速度为(0,200,0)m/s，目标场景中心坐标为(6,2,0)km，目标场景图如图3所示，其具体参数如表一所示。
表一
参数 符号 数值
载频 fc 9.65GHz
发射站零时刻位置 (xT,yT,hT) (-12km,2km,6km)
接收站零时刻位置 (xR,yR,hR) (0,0,4km)
平台运动速度 V0 200m/s
发射信号带宽 Br 200MHz
发射信号时宽 Tr 10us
脉冲采样频率 PRF 1500Hz
S2、对回波数据进行距离向脉冲压缩：利用发射的Chirp信号作为参考函数对回波数据进行距离向脉压，Chirp信号的表达式为：
S(τ)＝A0wr(τ)exp(jπKrτ2) (4)
取其反向共轭，建立距离向脉压信号的表达式：
S*(-τ)＝A0wr(-τ)exp(-jπKr(-τ)2) (5)
将回波数据的距离向数据与式(5)分别进行FFT后，在频域上相乘，然后进行IFFT得到距离向脉压后的回波数据，本实施例进行距离向脉冲压缩后的回波数据如图4所示。
S3、对场景中心点进行多普勒质心估计：选择回波数据条带的距离向中心距离门为场景中心所在距离门，本实施例的条带的中心所在距离门为4097，取出该距离门所在列并利用所提供的质心对其进行多普勒质心估计；具体实现方法为：
从回波数据表达式能够得出，多普勒质心fdc就是方位向频谱的中心，且方位向频谱相对于多普勒质心fdc为中心左右对称；
以估计的接收天线斜视角所计算出来的场景中心点位置为初始位置，计算该初始位置的多普勒质心并利用为迭代初值，采用以下步骤对场景中心点的多普勒质心fdc进行精确估计，如图5所示：
S3.1：分别计算小于的宽度频谱能量大于的宽度频谱能量
S3.2：比较与的大小，如果小于则转到步骤S3.3，否则转到步骤S3.6；
S3.3：按多普勒频谱分辨率增大即
S3.4：分别计算小于的宽度频谱能量大于的宽度频谱能量
S3.5：比较与的大小，如果小于则i＝i+1并转到步骤S3.3，否则输出结束；
S3.6：按多普勒频谱分辨率减小即
S3.7：分别计算小于的宽度频谱能量大于的宽度频谱能量
S3.8：比较与的大小，如果大于则i＝i+1并转到步骤S3.6，否则输出结束；
其中，i≥1，Δfdc为多普勒频谱分辨率，fa为方位向采样频率，Na为方位向采样点数，本实施例的计算方法得到的多普勒质心fdc值为1.7203kHz。
S4、场景中心点定位：根据场景中心点所在距离单元建立双基距离和方程，再根据S3中估计出的多普勒质心建立多普勒频率方程，然后求解场景中心点位置；具体实现方法为：
已知当t＝0时，双基距离和方程为
其中，(xdc,ydc)为场景中心点的坐标；
利用S3中得到的多普勒质心fdc，得到多普勒质心方程为：
双基距离和已知为Rb(Rb的计算方法为：本实施例计算得到的双机距离和为26.457km)。根据式(6)推导出：
其中，a＝R2-(yT-yR)2，b＝2yT(yT-yR)+2(xT-xR)xdc+k，
由于两个根号的影响，式(7)与式(8)无法解出精确解，因此采用搜索的方式来估计场景中心点的位置，如图6所示，其具体步骤如下：
S4.1：由接收机与发射机的位置和目标点的双基距离和(式6)，得到目标点所在椭圆曲线方程为y＝fellipse(x)；
S4.2：由测量的天线斜视角计算出目标点初始X轴坐标计算出初始Y轴坐标
S4.3：根据计算出多普勒质心
S4.4：令椭圆曲线上左侧点的X轴坐标为Y轴坐标为ydc_left＝fellipse(xdc_left)，代入式(7)中，计算其多普勒质心为fdc_left；
S4.5：令椭圆曲线上右侧点的X轴坐标为Y轴坐标为ydc_right＝fellipse(xdc_right)，代入式(7)中，计算出多普勒质心为fdc_right；
S4.6：比较|fdc-fdc_left|与|fdc-fdc_right的大小，如果|fdc-fdc_left|<|fdc-fdc_right|，则转到步骤S4.7；否则转到S4.10；
S4.7：令xdc＝xdc_left，fD＝fdc_left；
S4.8：令xleft＝xdc-Δx，则左侧点的位置坐标为(xleft,fellipse(xleft))，计算其多普勒质心，记为f'dc_left；
S4.9：比较|fdc-f'dc_left|与|fdc-fD|的大小，如果|fdc-f'dc_left|<＝|fdc-fD|，则转到S4.7；否则转到S4.13；
S4.10：令xdc＝xdc_right，fD＝fdc_right；
S4.11：令xright＝xdc+Δx，则右侧点的位置坐标为(xright,fellipse(xright))，计算其多普勒质心，记为f'dc_right；
S4.12：比较|fdc-f'dc_right|与|fdc-fD|的大小，如果|fdc-f'dc_right|<＝|fdc-fD|，则转到S4.10；否则转到S4.13；
S4.13：得到场景中心点坐标为(xdc,ydc)，结束；
其中，Δx为根据精度要求设定出来的搜索步长。利用本发明所提供的上述搜索算法，求解双基距离和方程与多普勒频率方程，从而解得场景中心点坐标为(xdc,ydc)，值为(6000.0 2000.2 0.0)m，与真实位置(6000 2000 0)m差0.2m。
S5、对回波数据进行聚焦：利用所求得到场景中心点在直角坐标系中的位置与接收站、发射站的相对位置对回波数据进行校正并聚焦，包括以下步骤：
S5.1、把步骤S4中得到场景中心点位置代入距离徒动公式中，得到回波数据中距离徒动校正量与多普勒频率的关系函数为：
式(9)中的未知系数由以下方程组求得：
其中，
S5.2、将回波数据进行方位向FFT变换，转换到距离多普勒域；
S5.3、在距离多普勒域进行距离徒动校正，本实施例得到聚焦后的图像如图7所示。
S6、场景中所有目标点定位：利用场景中心点和它在图像中的像素位置推导邻近像素点在所建立的三维地理坐标系下的坐标；其具体实现方法为：
通过场景中心点的地理坐标和它图像中的像素位置，推导出邻近像素点对应的坐标位置，具体推导方法包括以下步骤：
S6.1、已知某个像素点(i,j)对应的坐标为(xi,j,yi,j)，其中(i,j)分别为图像中距离向与方位向位置，它距离向的相邻点分别为(xi-1,j,yi-1,j)和(xi+1,j,yi+1,j)，首先，由波束模型可知：
对于近点，可得
其中，θc为接收天线斜视角；
最终得出：
其中，
S6.2、采用与步骤S61相同方法计算出(xi+1,j,yi+1,j)；
S6.3、它方位向上相邻点分别为(xi,j-1,yi,j-1)、(xi,j+1,yi,j+1)；
对于点(xi,j+1,yi,j+1)来说，发射机与接收机的新位置分别为(0,v0*PRI,HR)、(xT,yT+vy*PRI,HT)，已知关系式为：
其中θc为接收天线斜视角；
最终得出：
其中，
S6.4、采用与步骤S63相同方法计算出(xi,j-1,yi,j-1)；
S6.5、根据步骤S61～S64的计算方法得到整个图像中的每个像素所对应的地理坐标位置。
S7、目标场景地距投影：利用S6中求得的像素与坐标位置对应关系，将图像投影到地理坐标系中，得到地理坐标系下的真实图像；其具体投影关系为：
本实施例中对图3中9个点的定位投影后的图像如图8所示，其定位结果如表二所示。
表二
目标 真实坐标 像素位置 定位坐标 定位误差
1 (5900,1900)m (3590,3916) (5902.0,1900.1)m 2.00m
2 (6000,1900)m (3344,4096) (6000.0,1899.9)m 0.07m
3 (6100,1900)m (3097,4277) (6098.2,1899.8)m 1.81m
4 (5900,2000)m (4339,3916) (5902.0,1999.9)m 2.00m
5 (6000,2000)m (4094,4096) (6000.0,1999.9)m 0.07m
6 (6100,2000)m (3848,4277) (6098.2,1999.9)m 1.8m
7 (5900,2100)m (5090,3916) (5902.0,2100.1)m 2.00m
8 (6000,2100)m (4845,4096) (6000.0,2100.1)m 0.07m
9 (6100,2100)m (4599,4277) (6098.2,2100.0)m 1.8m
从图8和表二的结果来看，本发明提供的方法可以很好的实现移不变模式下双基地SAR的精确定位。
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。

Claims (8)

1.一种移不变机载双基合成孔径雷达目标定位方法，其特征在于，包括以下步骤：
S1、计算双基SAR的回波数据；
S2、对回波数据进行距离向脉冲压缩；
S3、对场景中心点进行多普勒质心估计：选择回波数据条带的距离向中心距离门为场景中心所在距离门，取出该距离门所在列并对其进行多普勒质心估计；
S4、场景中心点定位：根据场景中心点所在距离单元建立双基距离和方程，再根据S3中估计出的多普勒质心建立多普勒频率方程，然后求解场景中心点位置；
S5、对回波数据进行聚焦：利用所求得到场景中心点在直角坐标系中的位置与接收站、发射站的相对位置对回波数据进行校正并聚焦；
S6、场景中所有目标点定位：利用场景中心点和它在图像中的像素位置推导邻近像素点在所建立的三维地理坐标系下的坐标；
S7、目标场景地距投影：利用S6中求得的像素与坐标位置对应关系，将图像投影到地理坐标系中，得到地理坐标系下的真实图像。
2.根据权利要求1所述的移不变机载双基合成孔径雷达目标定位方法，其特征在于，所述的步骤S1具体实现方法为：令场景中心点被波束中心照射时刻，发射站固定，位置记为(xT,yT,hT)，其中，xT、yT和hT分别为发射站的x轴、y轴和z轴坐标；接收站位置记为(0,0,hR)，其中，0、0和hR分别为接收站的x轴、y轴和z轴坐标；接收站与发射站速度记为V，并沿y轴运动；建立了以接收站正下方为原点、高度为Z轴、速度方向为Y轴的三维地理坐标系：
将方位时间向量记为：
其中，PRI为脉冲重复间隔，Na为目标回波方位点数；
双基距离历史和为：Rb(t；x,y)＝RT(t；x,y)+RR(t；x,y)
其中，t为方位时间，RT(t；x,y)和RR(t；x,y)分别为发射站和接收站的距离历史：
<mrow> <msub> <mi>R</mi> <mi>T</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>;</mo> <mo>,</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <msub> <mi>x</mi> <mi>T</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>y</mi> <mo>+</mo> <mi>v</mi> <mi>t</mi> <mo>-</mo> <msub> <mi>y</mi> <mi>T</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msubsup> <mi>h</mi> <mi>T</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>R</mi> <mi>R</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>;</mo> <mo>,</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <msub> <mi>x</mi> <mi>R</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>y</mi> <mo>+</mo> <mi>v</mi> <mi>t</mi> <mo>-</mo> <msub> <mi>y</mi> <mi>R</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msubsup> <mi>h</mi> <mi>R</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
从而得到回波数据的表达式为：
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>s</mi> <mi>r</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>,</mo> <mi>&amp;tau;</mi> <mo>;</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>A</mi> <mn>0</mn> </msub> <msub> <mi>&amp;omega;</mi> <mi>r</mi> </msub> <mrow> <mo>(</mo> <mi>&amp;tau;</mi> <mo>-</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>;</mo> <mo>,</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> <msub> <mi>&amp;omega;</mi> <mi>a</mi> </msub> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>t</mi> <mo>-</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> </mrow> <msub> <mi>T</mi> <mi>a</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mi>exp</mi> <mo>&amp;lsqb;</mo> <mo>-</mo> <mi>j</mi> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mi>c</mi> </msub> <mfrac> <mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>;</mo> <mo>,</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>&amp;times;</mo> <mi>exp</mi> <mo>&amp;lsqb;</mo> <msub> <mi>j&amp;pi;K</mi> <mi>r</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>&amp;tau;</mi> <mo>-</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>;</mo> <mo>,</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
其中，A0是散射系数的幅度，ωr(·)为距离向包络，ωa(·)为方位向包络，τ是快时间变量，fc是载波频率，c是光速，Kr是距离向调频率，t0是目标点(x,y)的波束中心穿越时刻。
3.根据权利要求2所述的移不变机载双基合成孔径雷达目标定位方法，其特征在于，所述的步骤S2具体实现方法为：利用发射的Chirp信号作为参考函数对回波数据进行距离向脉压，Chirp信号的表达式为：
S(τ)＝A0wr(τ)exp(jπKrτ2) (4)
取其反向共轭，建立距离向脉压信号的表达式：
S*(-τ)＝A0wr(-τ)exp(-jπKr(-τ)2) (5)
将回波数据的距离向数据与式(5)分别进行FFT后，在频域上相乘，然后进行IFFT得到距离向脉压后的回波数据。
4.根据权利要求3所述的移不变机载双基合成孔径雷达目标定位方法，其特征在于，所述的步骤S3具体实现方法为：从回波数据表达式能够得出，多普勒质心fdc就是方位向频谱的中心，且方位向频谱相对于多普勒质心fdc为中心左右对称；
以估计的接收天线斜视角所计算出来的场景中心点位置为初始位置，计算该初始位置的多普勒质心并利用为迭代初值，采用以下步骤对场景中心点的多普勒质心fdc进行精确估计：
S3.1：分别计算小于的宽度频谱能量大于的宽度频谱能量
S3.2：比较与的大小，如果小于则转到步骤S3.3，否则转到步骤S3.6；
S3.3：按多普勒频谱分辨率增大即
S3.4：分别计算小于的宽度频谱能量大于的宽度频谱能量
S3.5：比较与的大小，如果小于则i＝i+1并转到步骤S3.3，否则输出结束；
S3.6：按多普勒频谱分辨率减小即
S3.7：分别计算小于的宽度频谱能量大于的宽度频谱能量
S3.8：比较与的大小，如果大于则i＝i+1并转到步骤S3.6，否则输出结束；
其中，i≥1，Δfdc为多普勒频谱分辨率，fa为方位向采样频率，Na为方位向采样点数。
5.根据权利要求4所述的移不变机载双基合成孔径雷达目标定位方法，其特征在于，所述的步骤S4具体实现方法为：
已知当t＝0时，双基距离和方程为
<mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>y</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>T</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>T</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msubsup> <mi>h</mi> <mi>T</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> <mo>+</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>R</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>R</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msubsup> <mi>h</mi> <mi>R</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
其中，(xdc,ydc)为场景中心点的坐标；
多普勒质心方程为：
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>f</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mover> <mi>v</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>T</mi> </msub> <mo>*</mo> <msub> <mover> <mi>R</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>T</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>T</mi> </msub> <mo>*</mo> <mi>&amp;lambda;</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <msub> <mover> <mi>v</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>R</mi> </msub> <mo>*</mo> <msub> <mover> <mi>R</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>R</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>R</mi> </msub> <mo>*</mo> <mi>&amp;lambda;</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>V</mi> <mrow> <mi>T</mi> <mi>x</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>T</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>V</mi> <mrow> <mi>T</mi> <mi>y</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>T</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&amp;lambda;</mi> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>T</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>T</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msubsup> <mi>h</mi> <mi>T</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <msub> <mi>V</mi> <mrow> <mi>R</mi> <mi>x</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>R</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>V</mi> <mrow> <mi>R</mi> <mi>y</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>R</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&amp;lambda;</mi> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>R</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>R</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msubsup> <mi>h</mi> <mi>R</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
双基距离和已知为Rb，根据式(6)推导出：
<mrow> <msub> <mi>y</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>T</mi> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>R</mi> </msub> <mo>)</mo> </mrow> <mi>b</mi> <mo>-</mo> <mn>2</mn> <msub> <mi>y</mi> <mi>T</mi> </msub> <msubsup> <mi>R</mi> <mi>b</mi> <mn>2</mn> </msubsup> <mo>&amp;PlusMinus;</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> <msqrt> <mrow> <msup> <mi>b</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>4</mn> <mi>a</mi> <mi>c</mi> </mrow> </msqrt> </mrow> <mrow> <mn>2</mn> <mi>a</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
其中，b＝2yT(yT-yR)+2(xT-xR)xdc+k，
由于两个根号的影响，式(7)与式(8)无法解出精确解，采用搜索的方式来估计场景中心点的位置，具体步骤如下：
S4.1：由接收站与发射站的位置和目标点的双基距离和式(6)，得到目标点所在椭圆曲线方程为y＝fellipse(x)；
S4.2：由测量的天线斜视角计算出目标点初始X轴坐标计算出初始Y轴坐标
S4.3：根据计算出多普勒质心
S4.4：令椭圆曲线上左侧点的X轴坐标为Y轴坐标为ydc_left＝fellipse(xdc_left)，代入式(7)中，计算其多普勒质心为fdc_left；
S4.5：令椭圆曲线上右侧点的X轴坐标为Y轴坐标为ydc_right＝fellipse(xdc_right)，代入式(7)中，计算出多普勒质心为fdc_right；
S4.6：比较|fdc-fdc_left|与|fdc-fdc_right|的大小，如果|fdc-fdc_left|<|fdc-fdc_right|，则转到步骤S4.7；否则转到S4.10；
S4.7：令xdc＝xdc_left，fD＝fdc_left；
S4.8：令xleft＝xdc-Δx，则左侧点的位置坐标为(xleft,fellipse(xleft))，计算其多普勒质心，记为f'dc_left；
S4.9：比较|fdc-f'dc_left|与|fdc-fD|的大小，如果|fdc-f'dc_left|<＝|fdc-fD|，则转到S4.7；否则转到S4.13；
S4.10：令xdc＝xdc_right，fD＝fdc_right；
S4.11：令xright＝xdc+Δx，则右侧点的位置坐标为(xright,fellipse(xright))，计算其多普勒质心，记为f'dc_right；
S4.12：比较|fdc-f'dc_right|与|fdc-fD|的大小，如果|fdc-f'dc_right|<＝|fdc-fD|，则转到S4.10；否则转到S4.13；
S4.13：得到场景中心点坐标为(xdc,ydc)，结束；
其中，Δx为根据精度要求设定出来的搜索步长。
6.根据权利要求5所述的移不变机载双基合成孔径雷达目标定位方法，其特征在于，所述的步骤S5具体实现方法为：
S5.1、把步骤S4中得到场景中心点位置代入距离徒动公式中，得到回波数据中距离徒动校正量与多普勒频率的关系函数：
<mrow> <msub> <mi>&amp;phi;</mi> <mrow> <mi>r</mi> <mi>c</mi> <mi>m</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>t</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <mn>4</mn> <mi>&amp;pi;</mi> <mo>&amp;lsqb;</mo> <msub> <mi>R</mi> <mi>x</mi> </msub> <mfrac> <mrow> <mfrac> <msub> <mi>f</mi> <mi>c</mi> </msub> <mi>c</mi> </mfrac> <mo>-</mo> <mfrac> <msub> <mi>a</mi> <mn>1</mn> </msub> <msub> <mi>V</mi> <mi>x</mi> </msub> </mfrac> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>cf</mi> <mi>t</mi> </msub> <mo>+</mo> <mn>2</mn> <msub> <mi>a</mi> <mn>1</mn> </msub> <msub> <mi>f</mi> <mi>c</mi> </msub> </mrow> <mrow> <mn>2</mn> <msub> <mi>V</mi> <mi>x</mi> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> <msqrt> <mrow> <msubsup> <mi>f</mi> <mi>c</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>cf</mi> <mi>t</mi> </msub> <mo>+</mo> <mn>2</mn> <msub> <mi>a</mi> <mn>1</mn> </msub> <msub> <mi>f</mi> <mi>c</mi> </msub> </mrow> <mrow> <mn>2</mn> <msub> <mi>V</mi> <mi>x</mi> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> <msub> <mi>cos&amp;theta;</mi> <mrow> <mi>s</mi> <mi>q</mi> <mi>x</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>x</mi> </msub> <mfrac> <mrow> <msub> <mi>a</mi> <mn>1</mn> </msub> <msub> <mi>sin&amp;theta;</mi> <mrow> <mi>s</mi> <mi>q</mi> <mi>x</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>V</mi> <mi>x</mi> </msub> <mi>c</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <msub> <mi>a</mi> <mn>0</mn> </msub> <mi>c</mi> </mfrac> <mo>&amp;rsqb;</mo> <msub> <mi>f</mi> <mi>r</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
式(9)中的未知系数由以下方程组求得：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>s</mi> <mi>q</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mi>a</mi> <mi> </mi> <mi>sin</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>k</mi> <mn>3</mn> </msub> <msqrt> <mrow> <mn>5</mn> <msubsup> <mi>k</mi> <mn>3</mn> <mn>2</mn> </msubsup> <mo>-</mo> <mn>3</mn> <msub> <mi>k</mi> <mn>2</mn> </msub> <msub> <mi>k</mi> <mn>4</mn> </msub> </mrow> </msqrt> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mi>x</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>k</mi> <mn>2</mn> </msub> <msub> <mi>k</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mn>5</mn> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&amp;theta;</mi> <mrow> <mi>s</mi> <mi>q</mi> <mi>x</mi> </mrow> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mn>2</mn> <msub> <mi>k</mi> <mn>4</mn> </msub> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&amp;theta;</mi> <mrow> <mi>s</mi> <mi>q</mi> <mi>x</mi> </mrow> </msub> <msub> <mi>sin&amp;theta;</mi> <mrow> <mi>s</mi> <mi>q</mi> <mi>x</mi> </mrow> </msub> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mi>x</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mn>3</mn> <msub> <mi>k</mi> <mn>2</mn> </msub> <msub> <mi>V</mi> <mi>x</mi> </msub> <msub> <mi>sin&amp;theta;</mi> <mrow> <mi>s</mi> <mi>q</mi> <mi>x</mi> </mrow> </msub> </mrow> <msub> <mi>k</mi> <mn>3</mn> </msub> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>a</mi> <mn>0</mn> </msub> <mo>=</mo> <mfrac> <msub> <mi>k</mi> <mn>0</mn> </msub> <mn>2</mn> </mfrac> <mo>-</mo> <msub> <mi>R</mi> <mi>x</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>a</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <msub> <mi>k</mi> <mn>1</mn> </msub> <mn>2</mn> </mfrac> <mo>-</mo> <msub> <mi>V</mi> <mi>x</mi> </msub> <msub> <mi>sin&amp;theta;</mi> <mrow> <mi>s</mi> <mi>q</mi> <mi>x</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
其中，
S5.2、将回波数据进行方位向FFT变换，转换到距离多普勒域；
S5.3、在距离多普勒域进行距离徒动校正。
7.根据权利要求6所述的移不变机载双基合成孔径雷达目标定位方法，其特征在于，所述的步骤S6具体实现方法为：通过场景中心点的地理坐标和它图像中的像素位置，推导出邻近像素点对应的坐标位置，具体推导方法包括以下步骤：
S6.1、已知某个像素点(i,j)对应的坐标为(xi,j,yi,j)，其中(i,j)分别为图像中距离向与方位向位置，它距离向的相邻点分别为(xi-1,j,yi-1,j)和(xi+1,j,yi+1,j)，首先，由波束模型可知：
<mrow> <mfrac> <mrow> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> </mrow> </mfrac> <mo>=</mo> <msub> <mi>tan&amp;theta;</mi> <mi>c</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
对于近点，可得
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>R</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>R</mi> <mo>=</mo> <msqrt> <mrow> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> <mn>2</mn> </msubsup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>R</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msubsup> <mi>h</mi> <mi>R</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> <mo>+</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>T</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>T</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msubsup> <mi>h</mi> <mi>T</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <msub> <mi>tan&amp;theta;</mi> <mi>c</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
其中，θc为接收天线斜视角；
最终得出：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <msub> <mi>b</mi> <mn>0</mn> </msub> <mo>&amp;PlusMinus;</mo> <msqrt> <mrow> <msup> <msub> <mi>b</mi> <mn>0</mn> </msub> <mn>2</mn> </msup> <mo>-</mo> <mn>4</mn> <msub> <mi>a</mi> <mn>2</mn> </msub> <msub> <mi>c</mi> <mn>0</mn> </msub> </mrow> </msqrt> </mrow> <mrow> <mn>2</mn> <msub> <mi>a</mi> <mn>2</mn> </msub> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <msub> <mi>tan&amp;theta;</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
其中，
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>a</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>4</mn> <msup> <mrow> <mo>(</mo> <mi>R</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>R</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msup> <mi>tan</mi> <mn>2</mn> </msup> <msub> <mi>&amp;theta;</mi> <mi>c</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mn>4</mn> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>T</mi> </msub> <mo>+</mo> <mo>(</mo> <mrow> <msub> <mi>y</mi> <mi>T</mi> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>R</mi> </msub> </mrow> <mo>)</mo> <msub> <mi>tan&amp;theta;</mi> <mi>c</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>b</mi> <mn>0</mn> </msub> <mo>=</mo> <mn>8</mn> <msub> <mi>tan&amp;theta;</mi> <mi>c</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>R</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>R</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <msub> <mi>tan&amp;theta;</mi> <mi>c</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>y</mi> <mi>R</mi> </msub> <mo>&amp;rsqb;</mo> <mo>+</mo> <mn>4</mn> <mo>&amp;lsqb;</mo> <msub> <mi>C</mi> <mn>0</mn> </msub> <mo>-</mo> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>T</mi> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>R</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <msub> <mi>tan&amp;theta;</mi> <mi>c</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mi>T</mi> </msub> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>T</mi> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>R</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>tan&amp;theta;</mi> <mi>c</mi> </msub> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>c</mi> <mn>0</mn> </msub> <mo>=</mo> <mn>4</mn> <msup> <mrow> <mo>(</mo> <mi>R</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>R</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>{</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <msub> <mi>tan&amp;theta;</mi> <mi>c</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>y</mi> <mi>R</mi> </msub> <mo>&amp;rsqb;</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msubsup> <mi>h</mi> <mi>R</mi> <mn>2</mn> </msubsup> <mo>}</mo> <mo>-</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>T</mi> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>R</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>tan&amp;theta;</mi> <mi>c</mi> </msub> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>C</mi> <mn>0</mn> </msub> <mo>&amp;rsqb;</mo> </mrow> <mn>2</mn> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>C</mi> <mn>0</mn> </msub> <mo>=</mo> <msubsup> <mi>x</mi> <mi>T</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>y</mi> <mi>T</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>h</mi> <mi>T</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>y</mi> <mi>R</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>h</mi> <mi>R</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mi>R</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>R</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
S6.2、采用与步骤S6.1相同方法计算出(xi+1,j,yi+1,j)；
S6.3、它方位向上相邻点分别为(xi,j-1,yi,j-1)、(xi,j+1,yi,j+1)；
对于点(xi,j+1,yi,j+1)来说，接收站与发射站的新位置分别为(0,v0*PRI,hR)、(xT,yT+vy*PRI,hT)，已知关系式为：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>R</mi> <mo>=</mo> <msqrt> <mrow> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>2</mn> </msubsup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>v</mi> <mn>0</mn> </msub> <mo>*</mo> <mi>P</mi> <mi>R</mi> <mi>I</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msubsup> <mi>h</mi> <mi>R</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> <mo>+</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>T</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>T</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>y</mi> </msub> <mo>*</mo> <mi>P</mi> <mi>R</mi> <mi>I</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msubsup> <mi>h</mi> <mi>T</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mrow> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>v</mi> <mn>0</mn> </msub> <mo>*</mo> <mi>P</mi> <mi>R</mi> <mi>I</mi> </mrow> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mfrac> <mo>=</mo> <msub> <mi>tan&amp;theta;</mi> <mi>c</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
最终得出：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <msub> <mi>b</mi> <mn>1</mn> </msub> <mo>&amp;PlusMinus;</mo> <msqrt> <mrow> <msubsup> <mi>b</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <mn>4</mn> <msub> <mi>a</mi> <mn>3</mn> </msub> <msub> <mi>c</mi> <mn>1</mn> </msub> </mrow> </msqrt> </mrow> <mrow> <mn>2</mn> <msub> <mi>a</mi> <mn>3</mn> </msub> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <msub> <mi>tan&amp;theta;</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>v</mi> <mn>0</mn> </msub> <mi>P</mi> <mi>R</mi> <mi>I</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> </mrow>
其中，
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>a</mi> <mn>3</mn> </msub> <mo>=</mo> <mn>4</mn> <msup> <mi>R</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msup> <mi>tan</mi> <mn>2</mn> </msup> <msub> <mi>&amp;theta;</mi> <mi>c</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mn>4</mn> <msup> <mrow> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mi>T</mi> </msub> <mo>+</mo> <msub> <mi>y</mi> <mi>T</mi> </msub> <msub> <mi>tan&amp;theta;</mi> <mi>c</mi> </msub> <mo>&amp;rsqb;</mo> </mrow> <mn>2</mn> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>b</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>4</mn> <mo>&amp;lsqb;</mo> <msub> <mi>C</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>2</mn> <msub> <mi>v</mi> <mn>0</mn> </msub> <msub> <mi>y</mi> <mi>T</mi> </msub> <mi>P</mi> <mi>R</mi> <mi>I</mi> <mo>&amp;rsqb;</mo> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mi>T</mi> </msub> <mo>+</mo> <msub> <mi>y</mi> <mi>T</mi> </msub> <msub> <mi>tan&amp;theta;</mi> <mi>c</mi> </msub> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>c</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>4</mn> <msup> <mi>R</mi> <mn>2</mn> </msup> <msubsup> <mi>h</mi> <mi>R</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <msub> <mi>C</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>2</mn> <msub> <mi>v</mi> <mn>0</mn> </msub> <msub> <mi>y</mi> <mi>T</mi> </msub> <mi>P</mi> <mi>R</mi> <mi>I</mi> <mo>&amp;rsqb;</mo> </mrow> <mn>2</mn> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>C</mi> <mn>1</mn> </msub> <mo>=</mo> <msubsup> <mi>x</mi> <mi>T</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>T</mi> </msub> <mo>+</mo> <msub> <mi>v</mi> <mn>0</mn> </msub> <mi>P</mi> <mi>R</mi> <mi>I</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msubsup> <mi>h</mi> <mi>T</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>v</mi> <mn>0</mn> </msub> <mi>P</mi> <mi>R</mi> <mi>I</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <msubsup> <mi>h</mi> <mi>R</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>R</mi> <mn>2</mn> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
S6.4、采用与步骤S6.3相同方法计算出(xi,j-1,yi,j-1)；
S6.5、根据步骤S6.1～S6.4的计算方法得到整个图像中的每个像素所对应的地理坐标位置。
8.根据权利要求7所述的移不变机载双基合成孔径雷达目标定位方法，其特征在于，所述的步骤S7中图像投影到地理坐标系中的具体投影关系为：
<mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> <mo>&amp;DoubleLeftRightArrow;</mo> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>)</mo> <mo>.</mo> </mrow>
CN201510151647.9A 2015-04-01 2015-04-01 一种移不变机载双基合成孔径雷达目标定位方法 Active CN104749570B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510151647.9A CN104749570B (zh) 2015-04-01 2015-04-01 一种移不变机载双基合成孔径雷达目标定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510151647.9A CN104749570B (zh) 2015-04-01 2015-04-01 一种移不变机载双基合成孔径雷达目标定位方法

Publications (2)

Publication Number Publication Date
CN104749570A CN104749570A (zh) 2015-07-01
CN104749570B true CN104749570B (zh) 2018-06-01

Family

ID=53589546

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510151647.9A Active CN104749570B (zh) 2015-04-01 2015-04-01 一种移不变机载双基合成孔径雷达目标定位方法

Country Status (1)

Country Link
CN (1) CN104749570B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105204020B (zh) * 2015-10-09 2017-07-14 电子科技大学 基于粒子群优化的移不变双基前视sar飞行模式设计方法
CN108267734B (zh) * 2016-12-30 2024-08-09 北京行易道科技有限公司 雷达
CN109283503B (zh) * 2017-07-23 2022-08-23 北京遥感设备研究所 一种均匀场景下的雷达回波信号多普勒质心确定方法
CN108132467B (zh) * 2017-12-23 2018-11-30 成都汇蓉国科微系统技术有限公司 基于增强型adc的dsp+fpga双基前视sar成像方法及成像装置
CN109765556A (zh) * 2018-12-29 2019-05-17 成都航天科工微电子系统研究院有限公司 一种基于级数反演的双基sar快速几何校正方法及装置
CN110361733B (zh) * 2019-07-01 2021-07-16 西安电子科技大学 一种基于时频联合重采样的中轨sar大斜视成像方法
CN112163454A (zh) * 2020-08-27 2021-01-01 中国海洋大学 基于rd谱增强的高频地波雷达杂波智能分类与定位方法
CN112560295B (zh) * 2021-02-20 2021-05-18 北京理工大学 一种用于被动合成孔径定位的卫星等效速度计算方法
CN116827458B (zh) * 2023-08-24 2023-12-05 北京理工大学 一种基于多角度的频偏估计及干扰源定位方法
CN117934617B (zh) * 2024-03-22 2024-05-28 中国科学院空天信息创新研究院 一种星载sar实时处理快速几何定位方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6747593B1 (en) * 2003-06-26 2004-06-08 Raytheon Company Generalized clutter tuning for bistatic radar systems
CN101000376A (zh) * 2007-01-08 2007-07-18 清华大学 双基地合成孔径雷达的双门限恒虚警运动目标检测方法
CN102565797A (zh) * 2011-12-21 2012-07-11 北京航空航天大学 一种针对聚束模式星载sar图像的几何校正方法
CN102890277A (zh) * 2012-09-19 2013-01-23 电子科技大学 移不变双基地合成孔径雷达距离徙动成像方法
CN103376447A (zh) * 2012-04-24 2013-10-30 中国人民解放军海军航空工程学院 一种实现非合作双多基地雷达目标三维定位的方法
CN103412300A (zh) * 2013-06-28 2013-11-27 南京航空航天大学 双基雷达目标三维定位方法
CN103412310A (zh) * 2013-08-26 2013-11-27 电子科技大学 双基地前视合成孔径雷达地面动目标检测方法与成像方法
CN103869314A (zh) * 2014-03-18 2014-06-18 电子科技大学 双飞移变双基地前视合成孔径雷达成像方法
CN103885062A (zh) * 2014-04-08 2014-06-25 电子科技大学 双基前视sar动目标成像方法与动目标速度估计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9335409B2 (en) * 2013-03-20 2016-05-10 Raytheon Company Bistatic inverse synthetic aperture radar imaging

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6747593B1 (en) * 2003-06-26 2004-06-08 Raytheon Company Generalized clutter tuning for bistatic radar systems
CN101000376A (zh) * 2007-01-08 2007-07-18 清华大学 双基地合成孔径雷达的双门限恒虚警运动目标检测方法
CN102565797A (zh) * 2011-12-21 2012-07-11 北京航空航天大学 一种针对聚束模式星载sar图像的几何校正方法
CN103376447A (zh) * 2012-04-24 2013-10-30 中国人民解放军海军航空工程学院 一种实现非合作双多基地雷达目标三维定位的方法
CN102890277A (zh) * 2012-09-19 2013-01-23 电子科技大学 移不变双基地合成孔径雷达距离徙动成像方法
CN103412300A (zh) * 2013-06-28 2013-11-27 南京航空航天大学 双基雷达目标三维定位方法
CN103412310A (zh) * 2013-08-26 2013-11-27 电子科技大学 双基地前视合成孔径雷达地面动目标检测方法与成像方法
CN103869314A (zh) * 2014-03-18 2014-06-18 电子科技大学 双飞移变双基地前视合成孔径雷达成像方法
CN103885062A (zh) * 2014-04-08 2014-06-25 电子科技大学 双基前视sar动目标成像方法与动目标速度估计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
临近空间慢速平台SAR地面动目标检测与成像;杨海光等;《电子科技大学学报》;20141130;第43卷(第6期);838-844 *
非移变模式机载双基地SAR 成像算法及实验结果;龚镇强等;《电子与信息学报》;20090131;第31卷(第1期);173-176 *

Also Published As

Publication number Publication date
CN104749570A (zh) 2015-07-01

Similar Documents

Publication Publication Date Title
CN104749570B (zh) 一种移不变机载双基合成孔径雷达目标定位方法
CN102393518B (zh) 一种适用于大斜视角的机载sar成像方法
CN103869311B (zh) 实波束扫描雷达超分辨成像方法
CN104833972B (zh) 一种双基地调频连续波合成孔径雷达频率变标成像方法
CN108427115B (zh) 合成孔径雷达对运动目标参数的快速估计方法
CN102749621B (zh) 一种双基地合成孔径雷达频域成像方法
CN109143237B (zh) 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法
CN105487074B (zh) 一种双基地合成孔径雷达数值距离多普勒成像方法
CN104898119B (zh) 一种基于相关函数的动目标参数估计方法
CN105182340B (zh) 一种双基地前视sar运动补偿方法
CN108872971B (zh) 一种基于运动单阵列的目标定位方法与装置
CN103336278A (zh) 多视角观测下前视三维sar成像方法
CN105137424A (zh) 一种杂波背景下实波束扫描雷达角超分辨方法
CN111273290B (zh) 基于预成像曲线轨迹的三维sar成像方法
CN114545411B (zh) 一种基于工程实现的极坐标格式多模高分辨sar成像方法
CN108710111A (zh) 一种机载双基前视sar方位相位的二维空变校正方法
CN105044667A (zh) 一种运动目标的双星跟踪方法、装置和系统
CN104820221B (zh) 多基合成孔径雷达的目标三维定位方法
CN109298417B (zh) 一种基于雷达信号处理的建筑物内部结构探测方法及装置
CN105158745A (zh) 移变双基前视合成孔径雷达距离徙动校正方法
Congfeng et al. Direct solution for fixed source location using well-posed TDOA and FDOA measurements
CN109143236B (zh) 适用于复杂飞行轨迹的双基聚束sar大场景成像方法
CN112285707B (zh) 基于gps导航信号的无源多基地高分辨成像方法
CN104931965A (zh) 一种固定发射站双基正侧视sar成像方法
CN116359921A (zh) 基于加速轨迹双基前视合成孔径雷达的快速时域成像方法

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