CN106556835B - 基于双基前视sar图像的目标定位方法 - Google Patents
基于双基前视sar图像的目标定位方法 Download PDFInfo
- Publication number
- CN106556835B CN106556835B CN201611092582.6A CN201611092582A CN106556835B CN 106556835 B CN106556835 B CN 106556835B CN 201611092582 A CN201611092582 A CN 201611092582A CN 106556835 B CN106556835 B CN 106556835B
- Authority
- CN
- China
- Prior art keywords
- receiver
- point
- scene
- center
- distance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
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/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
- 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)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于双基前视SAR图像的目标定位方法,以解决现有技术在广阔海域背景下,无法提取出具有标志性匹配点,最终导致匹配定位失效的问题。其实现方案是:利用相邻两合成孔径中心时刻接收机和发射机平台与目标的几何关系进行数学建模;根据模型寻找场景中心点的距离历程与接收机和发射机的几何关系,并结合校正到地平面的SAR图像得到目标与场景中心的精确位置关系,解算出目标相对接收机的位置信息,实现对目标的精确定位。本发明避免了从场景中提取参考点的操作,可实现在宽广海域背景下舰船目标的精确定位。
Description
技术领域
本发明属于雷达信号处理领域,特别涉及一种前视SAR图像的目标定位方法,可用于对海面上舰船目标SAR成像制导。
背景技术
合成孔径雷达SAR成像制导技术已经广泛应用于目标定位中,对于SAR目标定位技术的研究最早的研究工作可追溯到第二次世界大战以后,当时已有文献显示雷达图像作为定位的一种信息源。目前,实现SAR图像的目标定位的算法,主要有以下几种:
1.多项式对应关系的方法
该方法的原理是在SAR的测绘区域中找出一些位置确知的参考点,通过建立多项式图像相应点之间的坐标关系来预测图像上其他任意像素的地理位置。由于SAR图像的集合变形由多种因素引起,其变化规律十分复杂。为此,把SAR图像的总体变形看作是平移、缩放、旋转、仿射、偏扭、弯曲以及更高次的基本变形的综合作用结果,难以用严格的数字表达式来描述,而是用一个适当的多项式来描述校正前后。该方法对各种类型传感器图像的校正都是适用的,对参加计算的同名点的要求为在影像上为明显的地物点,易于判读,均匀分布。
这种方法的优点是操作简单且原理清晰,缺点是此方法回避了SAR成像的空间几何过程,直接对图像变形的本身进行数字模拟,在实际应用上具有内在的明显局限性,特别是在如海面、沙漠这些宽阔的场景背景下,很难找到用于图像定位的可靠参考点,无法取得匹配点。
2.空空定位理论
SAR的目标定位技术发展的第二阶段是在八十年代,在这一阶段SAR目标图像的集合构想数学模型,即共线方程模型,以该模型进行像素定位。他是建立在图像坐标与地面坐标的严格数学变换关系的基础上的,是对成像空间几何形态的直接描述。该方法的定位过程需要有地面高程信息,可以改正因地形起伏而引起的投影差。
该方法的优点是精度比多项式对应关系方法高,但是缺点也十分明显,共线方程需要数字高程信息,计算量比较大。利用共线方程还必须知道SAR平台的姿态参数等外方位元素,而这些参数是不稳定的,他们是时间的若干阶函数,这会导致该方法对目标的定位精度下降。可以看出,共线方程需要的数据量和数据类别很复杂,而且对于SAR成像而言,共线方程也只是一个虚构的中心投影构想模型,并没有利用SAR成像图像,并不能真实的反映出SAR的距离投影的集合形态,若再增加附加的约束条件,将使方程的可计算性变得很差。
上述提及的两种定位算法都通过对目标区域的图像进行参考点提取,与先前得到的基准图进行匹配得到匹配点,进而通过数学模型进行数字模拟反解出雷达平台与目标的相对位置,为平台轨迹修正提供参量。但当定位目标为宽广海面上的舰船时,图像背景被广阔海域占据,很难提取目标区域的参考点,进而无法提取匹配点,最终导致根据匹配点的定位失效。
发明内容
本发明的目的在于提出一种基于双基前视SAR图像几何定位方法,以解决上述现有技术在广阔海域背景下,无法提取出具有标志性匹配点,导致匹配定位失效的问题。
本发明的技术方案是:利用双基前视SAR在相邻两个合成孔径中心时刻的平台坐标与场景中心点的相对位置关系,进行几何建模,寻找场景中心点的距离历程与接收机和发射机平台的几何关系。并利用双基各平台的波束指向角,解算出场景中心距接收机的距离。根据得到的地距SAR图像,求解出目标与场景中心点的相对位置关系,建立起目标与接收机的几何关系,实现接收机对目标的定位。其实现步骤如下:
(1)根据双基前视SAR在相邻两个合成孔径中心时刻的平台坐标与场景中心点S的相对位置关系,建立双基平台与场景中心的几何模型;
(2)根据步骤(1)建立的几何模型,并利用双基收发平台的波束指向角和双基收发平台与场景中心点S的三角关系,解算得到Br时刻接收机平台与场景中心点S的距离Rrb;
(3)根据步骤(1)建立的几何模型,求解目标点P与场景中心点S的距离向距离Δx和方位向距离Δy;
(4)在步骤(1)建立的几何模型中设置辅助点G并结合步骤(2)的结果,计算辅助点G与接收机的距离Rrg;
(5)根据步骤(3)和步骤(4)的结果,求解目标点P相对于接收机的位置信息,即目标点P距接收机的距离Rrp,目标点距接收机的方位角θ,目标点距接收机的俯仰角φ:
本发明具有如下优点:
本发明通过建立前后相邻两个合成孔径中心时刻接收机和发射机的实际运行轨迹与目标点P、场景中心点S、辅助点G之间的几何模型,可得到目标点P相对接收机的距离、方位角和俯仰角,实现目标定位。
本发明的定位方法由于有SAR平台运行轨迹的介入,具有较好的实际意义,并且避免了从场景中提取参考点的操作,可在成像背景为宽广海面等难以提取参考点的场景中,对目标进行全程定位,弥补了利用场景参考点进行目标定位方法的不足。
附图说明
图1是本发明的实现流程图;
图2是本发明中建立的双基雷达定位模型图;
图3是本发明中的接收机相对目标点的几何关系示意图;
图4是本发明仿真使用的场景布置图;
图5是用本发明在接收机距目标的不同处仿真定位误差结果图。
具体实施方式
参照图1,本发明的实现步骤如下:
步骤1,建立接收机、发射机与场景中心点S的几何模型。
以目标场景中心点S为原心,建立北天东坐标系,如图2所示,其中x轴正方向指向正北方向,y轴正方向指向正东方向,z轴正方向指向垂直于地面方向;设接收机沿y轴正向从第一合成孔径中心Ar运动到第二合成孔径中心Br,与场景中心S构成了三角形ΔSArBr,为接收机在Ar时相对场景中心点S的擦地角,表示接收机在Br时相对场景中心点S的擦地角;同时发射机沿着与y轴不同指向的y'轴正向从第一合成孔径中心At运动到第二合成孔径中心Bt,与场景中心点S构成了三角形ΔSAtBt,设α'为角∠BtAtS,β'为角∠AtBtS。接收机在Ar时接收机与场景中心点S的距离为Rra,接收机在Br时接收机与场景中心点S的距离Rrb,发射机在At时发射机与场景中心点S的距离Rta,发射机在Bt时发射机与场景中心点S的距离为Rtb。
步骤2,求解当接收机位于Br时,接收机与场景中心点S的距离Rrb。
(2a)设Rra,Rta之和为Ra;Rrb,Rtb之和为Rb,即:
Ra=Rra+Rta <1>
Rb=Rrb+Rtb <2>
(2b)利用接收机在Ar、Br位置与场景中心点S构成的三角形ΔSArBr的正弦性质和余弦性质,得到关系式:
在<1>、<2>、<3>、<4>四个方程式中,Ra、Rb、α'、β'、为已知量,可通过几何校正后的SAR图像得到,v为接收机的运行速度,Δtm为相邻合成孔径中心的时间差,由雷达平台技术参数确定为已知量;Rra、Rta、Rrb、Rtb为未知量。
(2c)根据多项式理论,通过联立式<1>、式<2>、式<3>和式<4>这四个方程,求解出所述的四个未知量,得到接收机与场景中心点S的距离Rrb。
步骤3,确定目标点P与场景中心点S的位置关系。
(3a)从几何校正后的SAR图像中分别得到方位向分辨率ρa和目标点P相对场景中心点S的方位向间隔xn,计算目标点P相对于场景中心点S的方位向距离:Δx=ρa×xn;
(3b)从几何校正后的SAR图像中分别得到距离向分辨率ρr和目标点P相对场景中心点S的距离向间隔yn,计算目标点P相对于场景中心点S的距离向距离:Δy=ρr×yn。
步骤4,计算目标点P相对接收机的位置信息Rrp、θ和φ。
(4a)在几何校正后的SAR图像中按距离向和方位向按一定得间隔划出距离线和方位线,设置目标点P所在距离线与场景中心点S所在方位线的交点G作为辅助点,则G具有目标点P的距离向信息和场景中心点S的方位向信息;
(4b)利用辅助点G、目标点P和接收机第二合成孔径中心Br,构造三角形ΔPGBr利用三角形ΔPGBr的余弦性质得到辅助点G相对接收机的距离:
(4c)利用上述计算得到的Rrb、Δx、Δy和Rrg结果,结合图3所示的三角形ΔPGBr和三角形ΔSGBr,根据正弦性质和余弦性质,计算目标点P相对接收机的位置信息:
其中,Rrp是目标点P距接收机的距离,θ是目标点距接收机的方位角,φ是目标点距接收机的俯仰角。
本发明的效果通过以下仿真进一步说明:
1.仿真条件
1.1)设置仿真实验参数,如表1所示
表1接收机与发射机的雷达信息和运动参数
1.2)设置目标点P距场景中心点S的方位向距离Δx和距离向距离Δy均为100m如图4所示。
2.仿真实验:
仿真一:在不考虑双基“空间误差”等外部误差的情况下,验证本发明理论模型的正确性:
选取目标点P与接收机的几个典型距离进行模型验证,将利用模型计算的目标点P与接收机的距离结果与实际距离对比,距离误差为实际距离的10-9数量级,即成功验证本发明理论模型的正确性,结果如表2所示。
表2典型距离接收机对目标测距结果
实际接收机目标距离 | 测得的接收机目标距离 | 距离误差 |
69.062km | 69.062km | -0.0724m |
49.772km | 49.772km | -0.1004m |
35.035km | 35.035km | -0.1426m |
19.625km | 19.625km | -0.2542m |
9.9171km | 9.9165km | -0.5017m |
4.3782km | 4.377km | -1.1308m |
从表2可见,在不考虑外部误差的影响下,本方法能够对目标进行全程定位,辅助惯性导航***INS进行制导。
仿真二:考虑空间同步问题,在表1的数据中加入(±50m,±50m)的地平面
SAR图像空间误差,进一步验证本发明模型的有效性:
选取目标点P与接收机实际距离分别为69.06km、35.02km、9.90km、4.30km处进行仿真,结果如图5,其中:
图(5a)是接收机距目标69.02km处定位误差仿真结果图;
图(5b)是接收机距目标35.02km处定位误差仿真结果图;
图(5c)是接收机距目标9.90km处定位误差仿真结果图;
图(5d)是接收机距目标4.30km处定位误差仿真结果图。
从图5可见,各个仿真距离处的最大距离定位误差均不超过±50m,且随着目标点P与接收机距离越来越小,测量距离误差成收敛趋势;俯仰角和方位角误差均在10-3到10-2数量级,且随着目标点P与接收机距离越来越小,俯仰角和方位角的误差最大值仍不超过0.02度,满足精确定位精度要求,进一步说明了本方法的有效性。
Claims (1)
1.一种基于双基前视SAR图像的目标定位方法,包括如下步骤:
(1)根据双基前视SAR在相邻两个合成孔径中心时刻的平台坐标与场景中心点S的相对位置关系,建立双基平台与场景中心的几何模型:
(1a)是以目标场景中心点S为原心,建立北天东坐标系,即x轴正方向指向正北方向,y轴正方向指向正东方向,z轴正方向指向垂直于地面方向;
(1b)设接收机平台沿y轴正向运动,接收机相邻两个合成孔径中心为Ar和Br,发射机平台沿着与y轴不同指向的y'轴正向运动,发射机相邻两个合成孔径中心为At和Bt;接收机在Ar且发射机在At时,SAR平台与场景中心点S的距离和为Ra,接收机在Br且发射机在Bt时,SAR平台与场景中心点S的距离和为Rb;
(2)根据步骤(1)建立的几何模型,并利用双基收发平台的波束指向角和双基收发平台与场景中心点S的三角关系,解算接收机第二合成孔径中心Br与场景中心点S的距离Rrb:
其中,Rra表示接收机在Ar时接收机与场景中心点S的距离,Rta表示发射机在At时发射机与场景中心点S的距离,Ar表示接收机第一合成孔径中心,At表示发射机第一合成孔径中心,Ra为Rra与Rta之和;Rrb表示接收机第二合成孔径中心Br与场景中心点S的距离,Rtb表示发射机第二合成孔径中心Bt距场景中心点S的距离,Br表示接收机第二合成孔径中心,Bt表示发射机第二合成孔径中心,Rb为Rrb与Rtb之和;为发射机第一合成孔径中心Ar相对场景中心点S的擦地角,表示接收机第二合成孔径中心Br相对场景中心点S的擦地角,v为接收机的运行速度,Δtm为相邻合成孔径中心的时间差;α为∠BtAtS,β为∠AtBtS;
(3)根据步骤(1)建立的几何模型,求解目标点P与场景中心点S的距离向距离Δx和方位向距离Δy:
Δx=ρa×xn,Δy=ρr×yn
其中,ρa表示方位向分辨率,ρr表示距离向分辨率,xn表示目标点P相对场景中心点S的方位向间隔,yn表示目标点P相对场景中心点S的距离向间隔;
(4)在步骤(1)建立的几何模型中设置辅助点G并结合步骤(2)的结果,计算辅助点G与接收机的距离Rrg:
(4a)设置目标点P所在距离线与场景中心点S所在的方位线的交点G作为辅助点,构造三角形ΔPGBr,其中Br表示第二合成孔径中心时刻接收机的位置;
(4b)利用三角形ΔPGBr的余弦性质得到:
其中,Δy表示目标点P与场景中心点S的距离向距离,Rrb表示接收机第二合成孔径中心Br与场景中心点S的距离,表示接收机第二合成孔径中心Br相对场景中心点S的擦地角;
(5)根据步骤(3)和步骤(4)的结果,求解目标点P相对于接收机的位置信息,即目标点P距接收机的距离Rrp,目标点距接收机的方位角θ,目标点距接收机的俯仰角φ:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611092582.6A CN106556835B (zh) | 2016-12-01 | 2016-12-01 | 基于双基前视sar图像的目标定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611092582.6A CN106556835B (zh) | 2016-12-01 | 2016-12-01 | 基于双基前视sar图像的目标定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106556835A CN106556835A (zh) | 2017-04-05 |
CN106556835B true CN106556835B (zh) | 2019-03-26 |
Family
ID=58445667
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611092582.6A Active CN106556835B (zh) | 2016-12-01 | 2016-12-01 | 基于双基前视sar图像的目标定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106556835B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107092014B (zh) * | 2017-06-27 | 2020-05-01 | 西安电子科技大学 | 一种弹载双基前视sar舰船目标定位的优化方法 |
CN116564159A (zh) * | 2023-07-07 | 2023-08-08 | 中国科学院长春光学精密机械与物理研究所 | 一种光电测控设备跟踪操作仿真训练***及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412310A (zh) * | 2013-08-26 | 2013-11-27 | 电子科技大学 | 双基地前视合成孔径雷达地面动目标检测方法与成像方法 |
CN104678417A (zh) * | 2015-02-11 | 2015-06-03 | 中国地质大学(北京) | 一种基于导航卫星反射信号的目标定位方法 |
CN104865573A (zh) * | 2015-06-12 | 2015-08-26 | 西安电子科技大学 | 一种双基前视高机动平台sar扩展场景成像优化方法 |
CN104898120A (zh) * | 2015-06-12 | 2015-09-09 | 西安电子科技大学 | 一种基于回波模拟的双基前视高机动平台sar成像方法 |
-
2016
- 2016-12-01 CN CN201611092582.6A patent/CN106556835B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412310A (zh) * | 2013-08-26 | 2013-11-27 | 电子科技大学 | 双基地前视合成孔径雷达地面动目标检测方法与成像方法 |
CN104678417A (zh) * | 2015-02-11 | 2015-06-03 | 中国地质大学(北京) | 一种基于导航卫星反射信号的目标定位方法 |
CN104865573A (zh) * | 2015-06-12 | 2015-08-26 | 西安电子科技大学 | 一种双基前视高机动平台sar扩展场景成像优化方法 |
CN104898120A (zh) * | 2015-06-12 | 2015-09-09 | 西安电子科技大学 | 一种基于回波模拟的双基前视高机动平台sar成像方法 |
Non-Patent Citations (2)
Title |
---|
Fast raw data simulator of extended scenes for bistatic forward-looking synthetic aperture radar with constant acceleration;Ziqiang Meng 等;《2015 IEEE Signal Processing and Signal Processing Education Workshop》;20151231;237-242 |
基于斜距等效的弹载双基前视SAR相位空变校正方法;孟自强 等;《电子与信息学报》;20160331;第38卷(第3期);613-621 |
Also Published As
Publication number | Publication date |
---|---|
CN106556835A (zh) | 2017-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108614258B (zh) | 一种基于单水声信标距离量测的水下定位方法 | |
WO2022127532A1 (zh) | 一种激光雷达与imu的外参标定方法、装置及设备 | |
CN106990424B (zh) | 一种双天线gps测姿方法 | |
ES2530264T3 (es) | Procedimiento para la determinación de las coordenadas geográficas de píxeles en imágenes SAR | |
CN107044852B (zh) | 不平整状态下全站仪测量方法 | |
CN108061889A (zh) | Ais与雷达角度***偏差的关联方法 | |
CN108225185A (zh) | 一种车载扫描***检校方法 | |
CN104076348A (zh) | 一种雷达超视距基线无源协同定位方法 | |
CN111811538B (zh) | 一种基于常见星体的光电经纬仪定向标校方法 | |
CN109738902B (zh) | 一种基于同步信标模式的水下高速目标高精度自主声学导航方法 | |
CN103954953A (zh) | 一种基于数据驱动的机载激光雷达盲源误差补偿方法 | |
CN105242252B (zh) | 基于图像匹配的下降轨聚束sar雷达定位方法 | |
CN110441760B (zh) | 一种基于先验地形图的大范围海底地形图拓展构图方法 | |
CN106988312B (zh) | 基于北斗定向定位技术的机械设备归心改正方法及*** | |
CN108458698A (zh) | 全站仪不整平测量的未知点三维坐标计算方法 | |
CN106556835B (zh) | 基于双基前视sar图像的目标定位方法 | |
CN103744081A (zh) | 一种机载圆迹合成孔径雷达高精度三维成像补偿方法 | |
CN108957433A (zh) | 一种双基地雷达高精度测距装置及方法 | |
CN105067276B (zh) | 一种发动机推力线测量方法 | |
CN105929400A (zh) | 基于精确测量的恒定prf下变航速机载sar数据补偿方法 | |
CN110275139B (zh) | 一种基于旋转式基元复用的超短基线定位***及方法 | |
CN102207380A (zh) | 一种高精度横轴倾斜误差补偿方法 | |
Yu et al. | Automatic extrinsic self-calibration of mobile LiDAR systems based on planar and spherical features | |
CN103925904A (zh) | 一种基于对称测线的超短基线安装角度偏差无偏估计方法 | |
CN106353756A (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 |