CN111766581B - 高精度聚束式双基sar空间同步角度计算方法及装置 - Google Patents

高精度聚束式双基sar空间同步角度计算方法及装置 Download PDF

Info

Publication number
CN111766581B
CN111766581B CN202010667096.2A CN202010667096A CN111766581B CN 111766581 B CN111766581 B CN 111766581B CN 202010667096 A CN202010667096 A CN 202010667096A CN 111766581 B CN111766581 B CN 111766581B
Authority
CN
China
Prior art keywords
zynq
coordinates
unit
coordinate system
target area
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
CN202010667096.2A
Other languages
English (en)
Other versions
CN111766581A (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.)
Chengdu Huirong Guoke Microsystem Technology Co ltd
Original Assignee
Chengdu Huirong Guoke Microsystem Technology Co ltd
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 Chengdu Huirong Guoke Microsystem Technology Co ltd filed Critical Chengdu Huirong Guoke Microsystem Technology Co ltd
Priority to CN202010667096.2A priority Critical patent/CN111766581B/zh
Publication of CN111766581A publication Critical patent/CN111766581A/zh
Application granted granted Critical
Publication of CN111766581B publication Critical patent/CN111766581B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/9094—Theoretical aspects
    • 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/418—Theoretical aspects

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)
  • Radar Systems Or Details Thereof (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明提供一种基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法及装置，所述方法包括：记录高精度聚束式双基SAR成像的目标区域中心GPS坐标，实时获取运动平台当前的GPS坐标；结合地球模型参数，Zynq PS单元将GPS坐标系下目标区域中心GPS坐标和运动平台当前的GPS坐标转换成地心地固坐标系下的目标区域中心坐标和运动平台坐标；Zynq PL单元将地心地固坐标系下的目标区域中心坐标和运动平台坐标结合空间同步算法，计算出相控阵天线的离轴角度因子和旋转角度因子；由Zynq PS单元调用C语言数学函数库计算结果。根据本发明的方案，确保了计算结果的精度足够高以及数据计算足够快。

Description

高精度聚束式双基SAR空间同步角度计算方法及装置
技术领域
本发明涉及雷达成像处理领域，尤其涉及一种基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法及装置。
背景技术
双基合成孔径雷达系统通过分置的收发运动平台，可以实现接收平台前视区域的高分辨率成像，并且由于其收发分置的架构，使得该类型的雷达系统具有抗干扰性能好、通过发远收近的系统配置能够获得隐蔽性强、抗截获性好等优点，是当前雷达领域的研究热点。在早匹配制导、自主着陆等领域都有很好的应用前景。
双基SAR的两个运动平台位置分离，通常是发射平台斜视目标区域，接收平台前视目标区域，接收场景的回波信号。聚束式双基SAR工作模式和普通推扫式的双基SAR工作模式相似。聚束式双基SAR的特点在于发射平台和接收平台在运动的过程中，会实时控制它们的相控阵天线波束指向角始终指向同一片目标区域，持续地获取该目标区域的图像信息。在聚束式双基SAR的这种工作模式下，由于能够提前知道目标区域的坐标，因此，其波束方法与现有技术中普遍采用的方法不同。
在现有技术中，双基SAR系统的姿态信息和位置信息一般是由组合惯导系统以一定的帧率，例如100次/秒转发至双基SAR系统控制板上的FPGA芯片中，继而由FPGA根据空间同步算法计算出相应的相控阵天线波束指向角度，再发送至相控阵天线，完成一次指向角控制。然而，由于FPGA本身的定点计算特性，其计算结果的精度和动态范围不能同时兼顾，特别是在空间同步计算时，会涉及到对数值较大的参数进行计算，例如计算地球半径，在这种情况下，计算精度不足，会使得最后的天线波束指向角误差难以接受。
发明内容
为解决上述技术问题，本发明提出了一种基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法及装置，所述方法及装置，用以解决现有技术中对于聚束式双基SAR空间同步角度计算时，计算精度低、准确性差、计算速度慢的问题。
根据本发明的第一方面，提供一种基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法，所述方法包括以下步骤：
步骤S101：选定高精度聚束式双基SAR成像的目标区域，并记录目标区域中心GPS坐标(lonp,latp,hgtp)；
步骤S102：Zynq PS单元通过与组合惯导系统互连的UART接口实时获取运动平台当前的GPS坐标(lons,lats,hgts)；
步骤S103：结合地球模型参数，Zynq PS单元将GPS坐标系下所述目标区域中心GPS坐标和所述运动平台当前的GPS坐标转换成地心地固坐标系下的目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)；
步骤S104：Zynq PS单元将目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)通过与Zynq PL单元互连的AXI接口发送到Zynq PL单元中；
步骤S105：Zynq PL单元将目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)结合空间同步算法，计算出相控阵天线的离轴角度因子α和旋转角度因子φ；
步骤S106：Zynq PL单元将相控阵天线的离轴角度因子α和旋转角度因子φ通过AXI接口发送到Zynq PS单元，由Zynq PS单元调用C语言数学函数库计算结果。
进一步地，所述步骤S103：结合地球模型参数，Zynq PS单元将GPS坐标系下所述目标区域中心GPS坐标和所述运动平台当前的GPS坐标转换成地心地固坐标系下的目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)，包括：
Figure BDA0002580736260000031
Figure BDA0002580736260000032
Figure BDA0002580736260000033
xs1＝(Rgs+hgts)×cos(lats)×cos(lons) (4)
ys1＝(Rgs+hgts)×cos(lats)×sin(lons) (5)
zs1＝((1-ε2)×Rgs+hgts)sin(lats) (6)
xp1＝(Rgp+hgtp)×cos(latp)×cos(lonp) (7)
yp1＝(Rgp+hgtp)×cos(latp)×sin(lonp) (8)
zp1＝((1-ε2)×Rgp+hgtp)sin(latp) (9)
其中，Ra＝6378137米，为地球长轴；Rb＝6356752米，为地球短轴；(xs1,ys1,zs1)为运动平台在地心地固坐标系下的坐标，(xp1,yp1,zp1)为目标区域中心在地心地固坐标系下的坐标。
进一步地，所述步骤S105：Zynq PL单元将目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)结合空间同步算法，计算出相控阵天线的离轴角度因子α和旋转角度因子φ，包括：
步骤S1051：Zynq PL单元将地心地固坐标系下的坐标转换至运动平台坐标系，
xps2＝xp1-xs1
yps2＝yp1-ys1 (10)
zps2＝zp1-zs1
(xps2,yps2,zps2)为以运动平台为坐标中心的运动平台坐标系下目标区域中心的坐标；
步骤S1052：Zynq PL单元将运动平台坐标系下的目标区域中心的坐标旋转到北东天坐标系下，获得目标区域中心在北东天坐标系下的坐标(xps5,yps5,zps5)，计算方式如下：
Figure BDA0002580736260000041
Figure BDA0002580736260000042
xps5＝zps4
yps5＝xps4 (13)
zps5＝yps4
步骤S1053：Zynq PL单元将北东天坐标系下的目标区域中心的坐标转换成以运动平台为中心的惯性坐标系下的坐标(xps8,yps8,zps8)：
Figure BDA0002580736260000043
Figure BDA0002580736260000044
Figure BDA0002580736260000045
其中，at，βt，γt分别表示北东天坐标系下的运动平台的航向角、俯仰角和滚转角；
步骤S1054：在Zynq PL单元中计算出相控阵天线的离轴角度因子α和旋转角度因子φ：
Figure BDA0002580736260000046
φ＝zps8/yps8。
进一步地，所述步骤S106：Zynq PL单元将相控阵天线的离轴角度因子α和旋转角度因子φ通过AXI接口发送到Zynq PS单元，由Zynq PS单元调用C语言数学函数库计算结果，包括：
步骤S1061：Zynq PL单元将相控阵天线的离轴角度因子α和旋转角度因子φ通过AXI接口发送到Zynq PS单元；
步骤S1062：由Zynq PS单元调用C语言数学函数库中的反余弦和反正切函数，完成离轴角和旋转角的计算：
θL＝acos(a)
θX＝atan(φ)。
根据本发明第二方面，提供一种基于Zynq的高精度聚束式双基SAR空间同步角度计算的装置，所述装置包括：
第一坐标获取模块：选定高精度聚束式双基SAR成像的目标区域，并记录目标区域中心GPS坐标(lonp,latp,hgtp)；
第二坐标获取模块：Zynq PS单元通过与组合惯导系统互连的UART接口实时获取运动平台当前的GPS坐标(lons,lats,hgts)；
第一坐标转换模块：结合地球模型参数，Zynq PS单元将GPS坐标系下所述目标区域中心GPS坐标和所述运动平台当前的GPS坐标转换成地心地固坐标系下的目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)；
传输模块：Zynq PS单元将目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)通过与Zynq PL单元互连的AXI接口发送到Zynq PL单元中；
第一计算模块：Zynq PL单元将目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)结合空间同步算法，计算出相控阵天线的离轴角度因子α和旋转角度因子φ；
第二计算模块：Zynq PL单元将相控阵天线的离轴角度因子α和旋转角度因子φ通过AXI接口发送到Zynq PS单元，由Zynq PS单元调用C语言数学函数库计算结果。
进一步地，所述第一坐标转换，包括：
转换子模块：用于计算
Figure BDA0002580736260000061
Figure BDA0002580736260000062
Figure BDA0002580736260000063
xs1＝(Rgs+hgts)×cos(lats)×cos(lons) (4)
ys1＝(Rgs+hgts)×cos(lats)×sin(lons) (5)
zs1＝((1-ε2)×Rgs+hgts)sin(lats) (6)
xp1＝(Rgp+hgtp)×cos(latp)×cos(lonp) (7)
yp1＝(Rgp+hgtp)×cos(latp)×sin(lonp) (8)
zp1＝((1-ε2)×Rgp+hgtp)sin(latp) (9)
其中，Ra＝6378137米，为地球长轴；Rb＝6356752米，为地球短轴；(xs1,ys1,zs1)为运动平台在地心地固坐标系下的坐标，(xp1,yp1,zp1)为目标区域中心在地心地固坐标系下的坐标。
进一步地，所述第一计算模块，包括：
第一坐标转换子模块：Zynq PL单元将地心地固坐标系下的坐标转换至运动平台坐标系，
xps2＝xp1-xs1
yps2＝yp1-ys1 (10)
zps2＝zp1-zs1
(xps2,yps2,zps2)为以运动平台为坐标中心的运动平台坐标系下目标区域中心的坐标；
第二坐标转换子模块：Zynq PL单元将运动平台坐标系下的目标区域中心的坐标旋转到北东天坐标系下，获得目标区域中心在北东天坐标系下的坐标(xps5,yps5,zps5)，计算方式如下：
Figure BDA0002580736260000064
Figure BDA0002580736260000065
xps5＝zps4
yps5＝xps4 (13)
zps5＝yps4
第三坐标转换子模块：Zynq PL单元将北东天坐标系下的目标区域中心的坐标转换成以运动平台为中心的惯性坐标系下的坐标(xps8,yps8,zps8)：
Figure BDA0002580736260000071
Figure BDA0002580736260000072
Figure BDA0002580736260000073
其中，at，βt，γt分别表示北东天坐标系下的运动平台的航向角、俯仰角和滚转角；
第一计算子模块：在Zynq PL单元中计算出相控阵天线的离轴角度因子α和旋转角度因子φ：
Figure BDA0002580736260000074
φ＝zps8/yps8。
进一步地，所述第二计算模块，包括：
第一数据发送子模块：Zynq PL单元将相控阵天线的离轴角度因子α和旋转角度因子φ通过AXI接口发送到Zynq PS单元；
第二数据计算子模块：由Zynq PS单元调用C语言数学函数库中的反余弦和反正切函数，完成离轴角和旋转角的计算：
θL＝acos(α)
θX＝atan(φ)。
根据本发明第三方面，提供一种基于Zynq的高精度聚束式双基SAR空间同步角度计算的系统，包括：
处理器，用于执行多条指令；
存储器，用于存储多条指令；
其中，所述多条指令，用于由所述存储器存储，并由所述处理器加载并执行如前所述的一种基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法。
根据本发明第四方面，提供一种计算机可读存储介质，所述存储介质中存储有多条指令；所述多条指令，用于由处理器加载并执行如前所述的一种基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法。
根据本发明的上述方案，相对于现有技术，本发明充分利用软硬件协同设计的优势，将聚束式双基SAR空间同步角度计算算法中适合由软件执行的步骤放在Zynq ps中运行，适合以硬件逻辑执行步骤放在Zynq pl中运行，既确保了计算结果的精度足够高，也确保了数据计算足够快，可以充分发挥软硬件协同设计的优势。
上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，并可依照说明书的内容予以实施，以下以本发明的较佳实施例并配合附图详细说明如后。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解，本发明提供如下附图进行说明。在附图中：
图1为本发明一个实施方式的基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法流程图；
图2为本发明一个实施方式的基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法的交互图；
图3为本发明一个实施方式的Zynq PL单元将地心地固坐标系下的坐标进行转换的实现方式框图；
图4为本发明一个实施方式的Zynq PL单元将运动平台坐标系下的目标区域中心的坐标进行转换的实现方式框图；
图5为本发明一个实施方式的Zynq PL单元将北东天坐标系下的目标区域中心的坐标进行转换的实现方式框图；
图6为本发明一个实施方式的Zynq PL单元计算相控阵天线的离轴角度因子α和旋转角度因子φ的实现方式框图；
图7为本发明一个实施方式的基于Zynq的高精度聚束式双基SAR空间同步角度计算的装置结构框图。
具体实施方式
Zynq是Xilinx公司推出的一款集成了ARM-Cortex-A9双核(Processing System，PS)以及最多可达500多万个逻辑门的可编程逻辑单元(Programmable Logic，PL)的高端芯片。在芯片内部，Zynq PS单元与Zynq PL单元通过基于axi协议的高速接口互联。Zynq PS单元，其本质是一个arm处理单元，包含丰富的外设接口，例如spi、uart、Internet等。Zynq PS的编程实现方式和普通的arm处理器相似，例如用c语言编程，由指令代码串执行。Zynq PS单元可以执行双精度浮点型运算，有足够高的数值计算精度。Zynq PL单元本质上是一个FPGA逻辑处理单元，具有常见的FPGA内部资源。Zynq PL单元在执行算法时具有足够高的速度。
首先结合图1-2说明本发明的基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法。图1为基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法流程图，图2为基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法的交互图。如图1-2所示，所述方法包括以下步骤：
步骤S101：选定高精度聚束式双基SAR成像的目标区域，并记录目标区域中心GPS坐标(lonp,latp,hgtp)；
本实施例中采用WGS-84坐标系，其中，lon代表经度，lat代表纬度，hgt代表高度；lonp，latp，hgtp分别为目标区域中心在坐标系下的各个轴上的坐标；
步骤S102：Zynq PS单元通过与组合惯导系统互连的UART接口实时获取运动平台当前的GPS坐标(lons,lats,hgts)；
本实施例中采用WGS-84坐标系，其中，lon代表经度，lat代表纬度，hgt代表高度；lons，lats，hgts分别为运动平台在坐标系下的各个轴上的坐标；
步骤S103：结合地球模型参数，Zynq PS单元将GPS坐标系下所述目标区域中心GPS坐标和所述运动平台当前的GPS坐标转换成地心地固坐标系下的目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)；
所述地心地固坐标系，其原点O与地球质心重合，Z轴指向地球北极，X轴指向格林尼治平均子午面与赤道的交点，Y轴垂直于XOZ平面构成右手坐标系。
步骤S104：Zynq PS单元将目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)通过与Zynq PL单元互连的AXI接口发送到Zynq PL单元中；
步骤S105：Zynq PL单元将目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)结合空间同步算法，计算出相控阵天线的离轴角度因子α和旋转角度因子φ；
步骤S106：Zynq PL单元将相控阵天线的离轴角度因子α和旋转角度因子φ通过AXI接口发送到Zynq PS单元，由Zynq PS单元调用C语言数学函数库计算结果。
所述步骤S103：结合地球模型参数，Zynq PS单元将GPS坐标系下所述目标区域中心GPS坐标和所述运动平台当前的GPS坐标转换成地心地固坐标系下的目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)，包括：
Figure BDA0002580736260000111
Figure BDA0002580736260000112
Figure BDA0002580736260000113
xs1＝(Rgs+hgts)×cos(lats)×cos(lons) (4)
ys1＝(Rgs+hgts)×cos(lats)×sin(lons) (5)
zs1＝((1-ε2)×Rgs+hgts)sin(lats) (6)
xp1＝(Rgp+hgtp)×cos(latp)×cos(lonp) (7)
yp1＝(Rgp+hgtp)×cos(latp)×sin(lonp) (8)
zp1＝((1-ε2)×Rgp+hgtp)sin(latp) (9)
其中，Ra＝6378137米，为地球长轴；Rb＝6356752米，为地球短轴；(xs1,ys1,zs1)为运动平台在地心地固坐标系下的坐标，(xp1,yp1,zp1)为目标区域中心在地心地固坐标系下的坐标。
所述步骤S105：Zynq PL单元将目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)结合空间同步算法，计算出相控阵天线的离轴角度因子α和旋转角度因子φ，包括：
步骤S1051：Zynq PL单元将地心地固坐标系下的坐标转换至运动平台坐标系，
xps2＝xp1-xs1
yps2＝yp1-ys1 (10)
zps2＝zp1-zs1
运动平台坐标系仍处在地心地固坐标系下，所述运动平台坐标系的原点为运动平台在地心地固坐标系下运动平台的位置；
(xps2,yps2,zps2)为以运动平台为坐标中心的运动平台坐标系下目标区域中心的坐标；
本实施例中，Zynq PL单元将地心地固坐标系下的坐标进行转换的实现框图如图3所示。
步骤S1052：Zynq PL单元将运动平台坐标系下的目标区域中心的坐标旋转到北东天坐标系下，获得目标区域中心在北东天坐标系下的坐标(xps5,yps5,zps5)，计算方式如下：
Figure BDA0002580736260000121
Figure BDA0002580736260000122
xps5＝zps4
yps5＝xps4 (13)
zps5＝yps4
本实施例中，Zynq PL单元将运动平台坐标系下的目标区域中心的坐标进行转换的实现框图如图4所示。
本实施例中，使用Xilinx公司IP库中的cordic ip核完成正余弦运算。
步骤S1053：Zynq PL单元将北东天坐标系下的目标区域中心的坐标转换成以运动平台为中心的惯性坐标系下的坐标(xps8,yps8,zps8)：
Figure BDA0002580736260000123
Figure BDA0002580736260000124
Figure BDA0002580736260000125
其中，at，βt，γt分别表示北东天坐标系下的运动平台的航向角、俯仰角和滚转角；
本实施例中，Zynq PL单元将北东天坐标系下的目标区域中心的坐标进行转换的实现框图如图5所示。
步骤S1054：在Zynq PL单元中计算出相控阵天线的离轴角度因子α和旋转角度因子φ：
Figure BDA0002580736260000131
φ＝zps8/yps8
本实施例中，Zynq PL单元计算相控阵天线的离轴角度因子和旋转角度因子φ的实现框图如图6所示。
所述步骤S106：Zynq PL单元将相控阵天线的离轴角度因子α和旋转角度因子φ通过AXI接口发送到Zynq PS单元，由Zynq PS单元调用C语言数学函数库计算结果，包括：
步骤S1061：Zynq PL单元将相控阵天线的离轴角度因子α和旋转角度因子φ通过AXI接口发送到Zynq PS单元；
步骤S1062：由Zynq PS单元调用C语言数学函数库中的反余弦和反正切函数，完成离轴角和旋转角的计算：
θL＝acos(α)
θX＝atan(φ)。
本实施例中，部分数值较大的运算在Zynq PS中执行，部分非线性的三角函数运算由Xilinx公司公开的IP库中的cordic ip核完成。
本发明实施例进一步给出一种基于Zynq的高精度聚束式双基SAR空间同步角度计算的装置，如图7所示，所述装置包括：
第一坐标获取模块：选定高精度聚束式双基SAR成像的目标区域，并记录目标区域中心GPS坐标(lonp,latp,hgtp)；
第二坐标获取模块：Zynq PS单元通过与组合惯导系统互连的UART接口实时获取运动平台当前的GPS坐标(lons,lats,hgts)；
第一坐标转换模块：结合地球模型参数，Zynq PS单元将GPS坐标系下所述目标区域中心GPS坐标和所述运动平台当前的GPS坐标转换成地心地固坐标系下的目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)；
传输模块：Zynq PS单元将目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)通过与Zynq PL单元互连的AXI接口发送到Zynq PL单元中；
第一计算模块：Zynq PL单元将目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)结合空间同步算法，计算出相控阵天线的离轴角度因子α和旋转角度因子φ；
第二计算模块：Zynq PL单元将相控阵天线的离轴角度因子α和旋转角度因子φ通过AXI接口发送到Zynq PS单元，由Zynq PS单元调用C语言数学函数库计算结果。
本发明实施例进一步给出一种基于Zynq的高精度聚束式双基SAR空间同步角度计算的系统，包括：
处理器，用于执行多条指令；
存储器，用于存储多条指令；
其中，所述多条指令，用于由所述存储器存储，并由所述处理器加载并执行如前所述的一种基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法。
本发明实施例进一步给出一种计算机可读存储介质，所述存储介质中存储有多条指令；所述多条指令，用于由处理器加载并执行如前所述的基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法。
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
在本发明所提供的几个实施例中，应该理解到，所揭露的系统，装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。
上述以软件功能单元的形式实现的集成的单元，可以存储在一个计算机可读取存储介质中。上述软件功能单元存储在一个存储介质中，包括若干指令用以使得一台计算机装置(可以是个人计算机，实体机服务器，或者网络云服务器等，需安装Windows或者Windows Server操作系统)执行本发明各个实施例所述方法的部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(Read-Only Memory，ROM)、随机存取存储器(RandomAccess Memory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。

Claims (8)

1.一种基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法，其特征在于，包括以下步骤：
步骤S101：选定高精度聚束式双基SAR成像的目标区域，并记录目标区域中心GPS坐标(lonp,latp,hgtp)；
步骤S102：Zynq PS单元通过与组合惯导系统互连的UART接口实时获取运动平台当前的GPS坐标(lons,lats,hgts)；
步骤S103：结合地球模型参数，Zynq PS单元将GPS坐标系下所述目标区域中心GPS坐标和所述运动平台当前的GPS坐标转换成地心地固坐标系下的目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)；
步骤S104：Zynq PS单元将目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)通过与Zynq PL单元互连的AXI接口发送到Zynq PL单元中；
步骤S105：Zynq PL单元将目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)结合空间同步算法，计算出相控阵天线的离轴角度因子α和旋转角度因子φ；
步骤S106：Zynq PL单元将相控阵天线的离轴角度因子α和旋转角度因子φ通过AXI接口发送到Zynq PS单元，由Zynq PS单元调用C语言数学函数库计算结果；
所述步骤S105，包括：
步骤S1051：Zynq PL单元将地心地固坐标系下的坐标转换至运动平台坐标系，
Figure FDA0003411557570000011
(xps2,yps2,zps2)为以运动平台为坐标中心的运动平台坐标系下目标区域中心的坐标；
步骤S1052：Zynq PL单元将运动平台坐标系下的目标区域中心的坐标旋转到北东天坐标系下，获得目标区域中心在北东天坐标系下的坐标(xps5,yps5,zps5)，计算方式如下：
Figure FDA0003411557570000021
Figure FDA0003411557570000022
Figure FDA0003411557570000023
步骤S1053：Zynq PL单元将北东天坐标系下的目标区域中心的坐标转换成以运动平台为中心的惯性坐标系下的坐标(xps8,yps8,zps8)：
Figure FDA0003411557570000024
Figure FDA0003411557570000025
Figure FDA0003411557570000026
其中，at，βt，γt分别表示北东天坐标系下的运动平台的航向角、俯仰角和滚转角；
步骤S1054：在Zynq PL单元中计算出相控阵天线的离轴角度因子α和旋转角度因子φ：
Figure FDA0003411557570000027
φ＝zps8/yps8。
2.如权利要求1所述的基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法，其特征在于，所述步骤S103：结合地球模型参数，Zynq PS单元将GPS坐标系下所述目标区域中心GPS坐标和所述运动平台当前的GPS坐标转换成地心地固坐标系下的目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)，包括：
Figure FDA0003411557570000031
Figure FDA0003411557570000032
Figure FDA0003411557570000033
xs1＝(Rgs+hgts)×cos(lats)×cos(lons) (4)
ys1＝(Rgs+hgts)×cos(lats)×sin(lons) (5)
zs1＝((1-ε2)×Rgs+hgts)sin(lats) (6)
xp1＝(Rgp+hgtp)×cos(latp)×cos(lonp) (7)
yp1＝(Rgp+hgtp)×cos(latp)×sin(lonp) (8)
zp1＝((1-ε2)×Rgp+hgtp)sin(latp) (9)
其中，Ra＝6378137米，为地球长轴；Rb＝6356752米，为地球短轴；(xs1,ys1,zs1)为运动平台在地心地固坐标系下的坐标，(xp1,yp1,zp1)为目标区域中心在地心地固坐标系下的坐标。
3.如权利要求1所述的基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法，其特征在于，所述步骤S106：Zynq PL单元将相控阵天线的离轴角度因子α和旋转角度因子φ通过AXI接口发送到Zynq PS单元，由Zynq PS单元调用C语言数学函数库计算结果，包括：
步骤S1061：Zynq PL单元将相控阵天线的离轴角度因子α和旋转角度因子φ通过AXI接口发送到Zynq PS单元；
步骤S1062：由Zynq PS单元调用C语言数学函数库中的反余弦和反正切函数，完成离轴角和旋转角的计算：
θL＝acos(α)
θX＝atan(φ)。
4.一种基于Zynq的高精度聚束式双基SAR空间同步角度计算的装置，其特征在于，所述装置包括：
第一坐标获取模块：选定高精度聚束式双基SAR成像的目标区域，并记录目标区域中心GPS坐标(lonp,latp,hgtp)；
第二坐标获取模块：Zynq PS单元通过与组合惯导系统互连的UART接口实时获取运动平台当前的GPS坐标(lons,lats,hgts)；
第一坐标转换模块：结合地球模型参数，Zynq PS单元将GPS坐标系下所述目标区域中心GPS坐标和所述运动平台当前的GPS坐标转换成地心地固坐标系下的目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)；
传输模块：Zynq PS单元将目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)通过与Zynq PL单元互连的AXI接口发送到Zynq PL单元中；
第一计算模块：Zynq PL单元将目标区域中心坐标(xp1,yp1,zp1)和运动平台坐标(xs1,ys1,zs1)结合空间同步算法，计算出相控阵天线的离轴角度因子α和旋转角度因子φ；
第二计算模块：Zynq PL单元将相控阵天线的离轴角度因子α和旋转角度因子φ通过AXI接口发送到Zynq PS单元，由Zynq PS单元调用C语言数学函数库计算结果；
所述第一计算模块，包括：
第一坐标转换子模块：Zynq PL单元将地心地固坐标系下的坐标转换至运动平台坐标系，
Figure FDA0003411557570000041
(xps2,yps2,zps2)为以运动平台为坐标中心的运动平台坐标系下目标区域中心的坐标；
第二坐标转换子模块：Zynq PL单元将运动平台坐标系下的目标区域中心的坐标旋转到北东天坐标系下，获得目标区域中心在北东天坐标系下的坐标(xps5,yps5,zps5)，计算方式如下：
Figure FDA0003411557570000051
Figure FDA0003411557570000052
Figure FDA0003411557570000053
第三坐标转换子模块：Zynq PL单元将北东天坐标系下的目标区域中心的坐标转换成以运动平台为中心的惯性坐标系下的坐标(xps8,yps8,zps8)：
Figure FDA0003411557570000054
Figure FDA0003411557570000055
Figure FDA0003411557570000056
其中，at，βt，γt分别表示北东天坐标系下的运动平台的航向角、俯仰角和滚转角；
第一计算子模块：在Zynq PL单元中计算出相控阵天线的离轴角度因子α和旋转角度因子φ：
Figure FDA0003411557570000057
φ＝zps8/yps8。
5.如权利要求4所述的基于Zynq的高精度聚束式双基SAR空间同步角度计算的装置，其特征在于，所述第一坐标转换，包括：
转换子模块：用于计算：
Figure FDA0003411557570000061
Figure FDA0003411557570000062
Figure FDA0003411557570000063
xs1＝(Rgs+hgts)×cos(lats)×cos(lons) (4)
ys1＝(Rgs+hgts)×cos(lats)×sin(lons) (5)
zs1＝((1-ε2)×Rgs+hgts)sin(lats) (6)
xp1＝(Rgp+hgtp)×cos(latp)×cos(lonp) (7)
yp1＝(Rgp+hgtp)×cos(latp)×sin(lonp) (8)
zp1＝((1-ε2)×Rgp+hgtp)sin(latp) (9)
其中，Ra＝6378137米，为地球长轴；Rb＝6356752米，为地球短轴；(xs1,ys1,zs1)为运动平台在地心地固坐标系下的坐标，(xp1,yp1,zp1)为目标区域中心在地心地固坐标系下的坐标。
6.如权利要求5所述的基于Zynq的高精度聚束式双基SAR空间同步角度计算的装置，其特征在于，所述第二计算模块，包括：
第一数据发送子模块：Zynq PL单元将相控阵天线的离轴角度因子α和旋转角度因子φ通过AXI接口发送到Zynq PS单元；
第二数据计算子模块：由Zynq PS单元调用C语言数学函数库中的反余弦和反正切函数，完成离轴角和旋转角的计算：
θL＝acos(α)
θX＝atan(φ)。
7.一种基于Zynq的高精度聚束式双基SAR空间同步角度计算的系统，其特征在于，包括：
处理器，用于执行多条指令；
存储器，用于存储多条指令；
其中，所述多条指令，用于由所述存储器存储，并由所述处理器加载并执行如权利要求1-3之任一项所述基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法。
8.一种计算机可读存储介质，其特征在于，所述存储介质中存储有多条指令；所述多条指令，用于由处理器加载并执行如权利要求1-3之任一项所述基于Zynq的高精度聚束式双基SAR空间同步角度计算的方法。
CN202010667096.2A 2020-07-13 2020-07-13 高精度聚束式双基sar空间同步角度计算方法及装置 Active CN111766581B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010667096.2A CN111766581B (zh) 2020-07-13 2020-07-13 高精度聚束式双基sar空间同步角度计算方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010667096.2A CN111766581B (zh) 2020-07-13 2020-07-13 高精度聚束式双基sar空间同步角度计算方法及装置

Publications (2)

Publication Number Publication Date
CN111766581A CN111766581A (zh) 2020-10-13
CN111766581B true CN111766581B (zh) 2022-02-08

Family

ID=72725015

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010667096.2A Active CN111766581B (zh) 2020-07-13 2020-07-13 高精度聚束式双基sar空间同步角度计算方法及装置

Country Status (1)

Country Link
CN (1) CN111766581B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2816371A1 (en) * 2013-12-13 2014-12-24 Institute of Electronics, Chinese Academy of Sciences Method and device for steering attitude of satellite carrying synthetic aperture radar

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4172632A (en) * 1976-01-21 1979-10-30 Holmes Lawrence Jr Method and apparatus producing three-dimensional shadow images
CN102967851B (zh) * 2012-11-14 2014-03-26 电子科技大学 一种双基地sar的空间同步方法
CN104049241B (zh) * 2014-05-29 2017-01-04 电子科技大学 目标位置坐标未知的双基地合成孔径雷达的空间同步方法
CN111125873B (zh) * 2019-11-13 2023-09-29 北京航天长征飞行器研究所 一种飞行器天线波束选择控制器

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2816371A1 (en) * 2013-12-13 2014-12-24 Institute of Electronics, Chinese Academy of Sciences Method and device for steering attitude of satellite carrying synthetic aperture radar

Also Published As

Publication number Publication date
CN111766581A (zh) 2020-10-13

Similar Documents

Publication Publication Date Title
AU2018426350B2 (en) Systems and methods for identifying grids of geographical region in map
CN109000654B (zh) 定位方法、装置、设备及存储介质
CN107004961A (zh) 二维天线系统、用于定位目标的方法和设备
Groves The complexity problem in future multisensor navigation and positioning systems: A modular solution
US11069080B1 (en) Collaborative airborne object tracking systems and methods
CN111680596B (zh) 基于深度学习的定位真值校验方法、装置、设备及介质
US20200349727A1 (en) Obstacle detection method, mobile platform, and computer readable storage medium
CN113834482A (zh) 定位方法、装置、电子设备及存储介质
CN115079168B (zh) 基于激光雷达与毫米波雷达融合的建图方法、装置及设备
CN114637329A (zh) 一种飞机多机密集编队队形重构方法及系统
CN115164936A (zh) 高精地图制作中用于点云拼接的全局位姿修正方法及设备
RU2513900C1 (ru) Способ и устройство определения координат объектов
CN111766581B (zh) 高精度聚束式双基sar空间同步角度计算方法及装置
Soloviev et al. Reconfigurable Integration Filter Engine (RIFE) for Plug-and-Play Navigation
CN117214932A (zh) 卫星编队的波束位置的确定方法、装置、设备及存储介质
WO2021056503A1 (zh) 可移动平台的定位方法、装置、可移动平台及存储介质
Koks Using rotations to build aerospace coordinate systems
CN114518585B (zh) 基于单颗卫星的目标定位方法、装置、终端及存储介质
CN113297780B (zh) 载体遮挡与多径信号的高速仿真方法、装置、介质及设备
CN112230194B (zh) 一种基于平移阵列的解模糊方法、设备及存储介质
CN115147482A (zh) 一种位姿初始化方法
CN115079102A (zh) 可灵活配置的雷达监测态势显示方法及装置
JP2008241321A (ja) Ｒｃｓ制御装置
Sulaj et al. Examples of real-time UAV data processing with cloud computing
CN113960614A (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