CN112255599A - 实现卫星雷达波范围在地图上展示的方法、装置及设备 - Google Patents

实现卫星雷达波范围在地图上展示的方法、装置及设备 Download PDF

Info

Publication number
CN112255599A
CN112255599A CN202011109983.4A CN202011109983A CN112255599A CN 112255599 A CN112255599 A CN 112255599A CN 202011109983 A CN202011109983 A CN 202011109983A CN 112255599 A CN112255599 A CN 112255599A
Authority
CN
China
Prior art keywords
satellite
coordinate system
coordinate
vector
calculating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011109983.4A
Other languages
English (en)
Other versions
CN112255599B (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.)
Beijing Shi Guan Jin Yang Technology Development Co ltd
Original Assignee
Beijing Shi Guan Jin Yang Technology Development 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 Beijing Shi Guan Jin Yang Technology Development Co ltd filed Critical Beijing Shi Guan Jin Yang Technology Development Co ltd
Priority to CN202011109983.4A priority Critical patent/CN112255599B/zh
Priority claimed from CN202011109983.4A external-priority patent/CN112255599B/zh
Publication of CN112255599A publication Critical patent/CN112255599A/zh
Application granted granted Critical
Publication of CN112255599B publication Critical patent/CN112255599B/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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40—Means for monitoring or calibrating
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Remote Sensing (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本申请公开了一种实现卫星雷达波范围在地图上展示的方法、装置及设备，以实现对卫星仿真结果进行展示，该方法包括：获取卫星的轨道六根数参数，计算卫星在J2000坐标系下的实时位置及实时速度；根据实时位置、实时速度及朝向参数，计算在J2000坐标系下的欧拉角并转换为在WGS84坐标系下的第二向量坐标；根据实时位置计算卫星的经纬度；根据第二向量坐标以及经纬度，计算卫星在北东地坐标系下的第三向量坐标，从而计算卫星在北东地坐标系下的欧拉角；根据第三向量坐标、卫星在北东地坐标系下的欧拉角及卫星对应的雷达在北东地坐标系下的欧拉角，获得以雷达为起点，预设角度为圆周角度的射线；计算射线与地球的交点坐标；将交点坐标对应的区域在地图上展示。

Description

实现卫星雷达波范围在地图上展示的方法、装置及设备
技术领域
本申请涉及装备仿真技术领域，具体涉及一种实现卫星雷达波范围在地图上展示的方法、装置及设备。
背景技术
雷达卫星是载有合成孔径雷达的对地观测遥感卫星的统称。在现有技术中，可以通过对雷达卫星进行仿真，以模拟得到雷达卫星的观测范围。一般由于雷达卫星涉及的参数较多，在对雷达卫星进行仿真的过程中，缺乏一种可以直观对仿真得到的雷达卫星观测范围进行展示的方法。
发明内容
有鉴于此，本申请实施例提供一种实现卫星雷达波范围在地图上展示的方法、装置及设备，以实现对仿真得到的卫星雷达波范围在地图上进行展示。
为解决上述问题，本申请实施例提供的技术方案如下：
一种实现卫星雷达波范围在地图上展示的方法，所述方法包括：
获取卫星的轨道六根数参数，根据所述轨道六根数参数计算所述卫星在J2000坐标系下的实时位置以及实时速度；
根据所述实时位置、所述实时速度以及所述卫星的朝向参数，计算所述卫星在J2000坐标系下的欧拉角；
将所述卫星在J2000坐标系下的欧拉角转换为在J2000坐标系下的第一向量坐标；
将所述第一向量坐标转换为在WGS84坐标系下的第二向量坐标；
根据所述实时位置计算所述卫星对应的经纬度；
根据所述第二向量坐标以及所述经纬度，计算所述卫星在北东地坐标系下的第三向量坐标；
根据所述第三向量坐标，计算所述卫星在北东地坐标系下的欧拉角；
根据所述第三向量坐标、所述卫星在北东地坐标系下的欧拉角以及所述卫星对应的雷达在北东地坐标系下的欧拉角，获得以所述雷达为起点，预设角度为圆周角度的多条射线；
计算所述射线与地球的交点坐标；
将所述交点坐标对应的区域在地图上进行展示。
在一种可能的实现方式中，所述方法还包括：
如果所述交点坐标对应的区域经过经度180度，将所述交点坐标对应的区域沿所述经度180度进行切割，将切割后的区域在地图上进行展示。
在一种可能的实现方式中，所述将所述卫星在J2000坐标系下的欧拉角转换为在J2000坐标系下的第一向量坐标，包括：
利用以下公式，将所述卫星在J2000坐标系下的欧拉角转换为在J2000坐标系下的第一向量坐标：
X1＝cos(yaw1·π/180)*cos(pitch1·π/180)；
Y1＝sin(yaw1·π/180)*cos(pitch1·π/180)；
Z1＝sin(pitch1·π/180)*(-1)；
其中，yaw1为所述卫星在J2000坐标系下的欧拉角中的偏航角，pitch1为所述卫星在J2000坐标系下的欧拉角中的俯仰角；X1、Y1和Z1为J2000坐标系下的第一向量坐标。
在一种可能的实现方式中，所述根据所述实时位置计算所述卫星对应的经纬度，包括：
利用以下公式，根据所述实时位置计算所述卫星对应的经纬度：
L＝arctan(Y/X)；
Figure BDA0002728282510000031
HD＝Z/sinB-N(1-e2)；
其中，X、Y和Z为所述实时位置，N为法截弧的曲率半径，e为自然常数，L为所述卫星对应的经度、B为所述卫星对应的纬度、HD为所述卫星对应的大地高度。
在一种可能的实现方式中，所述根据所述第二向量坐标以及所述经纬度，计算所述卫星在北东地坐标系下的第三向量坐标，包括：
利用以下公式，根据所述第二向量坐标以及所述经纬度，计算所述卫星在北东地坐标系下的第三向量坐标：
X3＝-sin(B·π/180)*cos(L·π/180)-sin(B·π/180)*sin(L·π/180)*X2+cos(B·π/180)*Z2；
Y3＝-sin(L·π/180)*Y2+cos(L·π/180)*X2；
Z3＝-(cos(B·π/180)*cos(L·π/180)*Y2+cos(B·π/180)*sin(L·π/180)*X2+sin(B·π/180)*Z2)；
其中，X2、Y2和Z2为所述第二向量坐标；L为所述卫星对应的经度、B为所述卫星对应的纬度、HD为所述卫星对应的大地高度；X3、Y3和Z3为所述卫星在北东地坐标系下的第三向量坐标。
在一种可能的实现方式中，所述根据所述第三向量坐标，计算所述卫星在北东地坐标系下的欧拉角，包括：
利用以下公式，根据所述第三向量坐标，计算所述卫星在北东地坐标系下的欧拉角：
roll2＝0；
Figure BDA0002728282510000032
yaw2＝(atan2(Y3,X3))·180/π；
其中，X3、Y3和Z3为所述第三向量坐标，yaw2为所述卫星在北东地坐标系下的欧拉角中的偏航角，pitch2为所述卫星在北东地坐标系下的欧拉角中的俯仰角，roll2为所述卫星在北东地坐标系下的欧拉角中的翻滚角。
一种实现卫星雷达波范围在地图上展示的装置，所述装置包括：
获取单元，用于获取卫星的轨道六根数参数，根据所述轨道六根数参数计算所述卫星在J2000坐标系下的实时位置以及实时速度；
第一计算单元，用于根据所述实时位置、所述实时速度以及所述卫星的朝向参数，计算所述卫星在J2000坐标系下的欧拉角；
第一转换单元，用于将所述卫星在J2000坐标系下的欧拉角转换为在J2000坐标系下的第一向量坐标；
第二转换单元，用于将所述第一向量坐标转换为在WGS84坐标系下的第二向量坐标；
第二计算单元，用于根据所述实时位置计算所述卫星对应的经纬度；
第三计算单元，用于根据所述第二向量坐标以及所述经纬度，计算所述卫星在北东地坐标系下的第三向量坐标；
第四计算单元，用于根据所述第三向量坐标，计算所述卫星在北东地坐标系下的欧拉角；
获得单元，用于根据所述第三向量坐标、所述卫星在北东地坐标系下的欧拉角以及所述卫星对应的雷达在北东地坐标系下的欧拉角，获得以所述雷达为起点，预设角度为圆周角度的多条射线；
第五计算单元，用于计算所述射线与地球的交点坐标；
展示单元，用于将所述交点坐标对应的区域在地图上进行展示。
在一种可能的实现方式中，所述装置还包括：
切割单元，用于如果所述交点坐标对应的区域经过经度180度，将所述交点坐标对应的区域沿所述经度180度进行切割；
所述展示单元，具体用于将切割后的区域在地图上进行展示。
一种实现卫星雷达波范围在地图上展示的设备，包括：存储器，处理器，及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时，实现上述的实现卫星雷达波范围在地图上展示的方法。
一种计算机可读存储介质，所述计算机可读存储介质中存储有指令，当所述指令在终端设备上运行时，使得所述终端设备执行上述的实现卫星雷达波范围在地图上展示的方法。
由此可见，本申请实施例具有如下有益效果：
本申请实施例可以根据卫星的实时位置、实时速度以及朝向参数等，计算出卫星在J2000坐标系下的欧拉角，并通过在各个坐标系下的向量坐标最终转换为在北东地坐标系下的欧拉角。同时，结合卫星挂载的雷达在北东地坐标系下的欧拉角，可以计算得到雷达波的发射方向，从而计算出雷达波与地球的交点，将雷达波与地球的交点对应的区域在地图上进行展示，从而直观对雷达波的扫描范围在地图上进行展示。
附图说明
图1为本申请实施例提供的实现卫星雷达波范围在地图上展示的方法的流程图；
图2为本申请实施例提供的一种卫星雷达波范围在地图上展示的效果图；
图3为本申请实施例提供的另一种卫星雷达波范围在地图上展示的效果图；
图4为本申请实施例提供的实现卫星雷达波范围在地图上展示的装置的示意图。
具体实施方式
为使本申请的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本申请实施例作进一步详细的说明。
为了便于理解和解释本申请实施例提供的技术方案，下面将先对本申请的背景技术进行说明。
发明人在对传统的对雷达卫星进行仿真的方法进行研究后发现，可以通过对雷达卫星进行仿真，以模拟得到雷达卫星的观测范围。一般由于雷达卫星涉及的参数较多，在对雷达卫星进行仿真的过程中，缺乏一种可以直观对仿真得到的雷达卫星观测范围进行展示的方法。
基于此，本申请实施例提供了实现卫星雷达波范围在地图上展示的方法，可以根据卫星等装备的姿态参数、以及雷达相对挂载卫星等装备的姿态参数，使用向量法计算出雷达波发射波的向量方向。从而计算出静态的雷达投影效果，可以在地图上相对准确地画出雷达波效果，对雷达波扫描范围进行展示。进一步的，为了让工程师可以更加直观的可以在二维地图上查看雷达波探测范围，根据探测范围还可以调整装备挂载雷达的参数，来提高卫星实际使用的准确性和可靠性，对工程师能够完成仿真中装备的使用提供了更加直观和准确的依据。
为了便于理解本申请，下面结合附图对本申请实施例提供的一种实现卫星雷达波范围在地图上展示的方法进行说明。
参见图1，该图为本申请实施例提供的一种实现卫星雷达波范围在地图上展示的方法，如图2所示，该方法可以包括以下步骤：
步骤S101：获取卫星的轨道六根数参数，根据轨道六根数参数计算卫星在J2000坐标系下的实时位置以及实时速度。
卫星的轨道六根数参数，是描述在牛顿运动定律和牛顿万有引力定律的作用下的卫星在其开普勒轨道上运动时，确定其轨道所必要的六个参数。一般包括：半长轴、轨道离心率、轨道倾角、近地点辐角、升交点赤经、真近点角。
在实际应用中，用户可以根据仿真需求输入卫星的轨道六根数参数。同时还可以输入仿真开始时间以及仿真步长。则在获取卫星的轨道六根数参数、仿真开始时间以及仿真步长后，可以根据二体理论计算得到卫星在J2000坐标系下的实时位置以及实时速度。其中，J2000坐标系是根据世界协调时间UTC时间为2000年1月1日11:58：55.816时刻平赤道作为参考平面，平春分点作为主方向，组成的坐标系，即J2000平赤道平春分点坐标系，简称J2000坐标系。
步骤S102：根据实时位置、实时速度以及卫星的朝向参数，计算卫星在J2000坐标系下的欧拉角。
根据步骤S101计算得到的实时位置、实时速度，以及输入的卫星的朝向参数，例如对地，对日，对月，还是切向地等，可以求出卫星在J2000坐标系下的欧拉角。欧拉角是物体绕坐标系三个坐标轴的旋转角度，一般包括偏航角、俯仰角以及翻滚角。
步骤S103：将卫星在J2000坐标系下的欧拉角转换为在J2000坐标系下的第一向量坐标。
为了便于后续计算，可以将卫星在J2000坐标系下的欧拉角转换为在J2000坐标系下的第一向量坐标，第一向量坐标具有方向，利用向量坐标进行计算比直接使用欧拉角计算更为快速、准确。
在一种可能的实现方式中，可以利用以下公式，将卫星在J2000坐标系下的欧拉角转换为在J2000坐标系下的第一向量坐标：
X1＝cos(yaw1·π/180)*cos(pitch1·π/180)；
Y1＝sin(yaw1·π/180)*cos(pitch1·π/180)；
Z1＝sin(pitch1·π/180)*(-1)；
其中，yaw1为卫星在J2000坐标系下的欧拉角中的偏航角，pitch1为卫星在J2000坐标系下的欧拉角中的俯仰角。X1、Y1和Z1为J2000坐标系下的第一向量坐标，即X1、Y1和Z1分别为J2000坐标系下第一向量坐标对应于三个坐标轴的坐标。
步骤S104：将第一向量坐标转换为在WGS84坐标系下的第二向量坐标。
在步骤S103将卫星在J2000坐标系下的欧拉角转换为卫星在J2000坐标系下的第一向量坐标后，可以进一步将卫星在J2000坐标系下的第一向量坐标转换为WGS84坐标系下的第二向量坐标。其中，WGS84坐标系是一种国际上采用的地心坐标系。其坐标原点为地球质心，其地心空间直角坐标系的Z轴指向BIH(国际时间服务机构)1984.0定义的协议地球极(CTP)方向，X轴指向BIH 1984.0的零子午面和CTP赤道的交点，Y轴与Z轴、X轴垂直构成右手坐标系，称为1984年世界大地坐标系统。
步骤S105：根据实时位置计算卫星对应的经纬度。
根据步骤S101计算得到的实时位置还可以计算卫星对应的经纬度，卫星对应的经纬度包括卫星对应的经度、纬度以及大地高度。
在一种可能的实现方式中，可以利用以下公式，根据实时位置计算卫星对应的经纬度：
L＝arctan(Y/X)；
Figure BDA0002728282510000081
HD＝Z/sinB-N(1-e2)；
其中，X、Y和Z为实时位置，即X、Y和Z分别为卫星在J2000坐标系下实时位置对应于三个坐标轴的坐标。N为法截弧的曲率半径，e为自然常数，L为卫星对应的经度、B为卫星对应的纬度、HD为卫星对应的大地高度。
需要注意的是，本申请实施例不限定步骤S103-S104与步骤S105之间的顺序，可以先执行步骤S102-S104再执行步骤S105，也可以先执行步骤S105再执行步骤S102-S104，还可以步骤S102-S104与步骤S105并行执行。
步骤S106：根据第二向量坐标以及经纬度，计算卫星在北东地坐标系下的第三向量坐标。
在得到卫星在WGS84坐标系下的第二向量坐标以及卫星对应的经纬度之后，可以计算得到卫星在北东地坐标系下的第三向量坐标。其中，北东地坐标系简称为n坐标系，也叫做导航坐标系，是在导航时根据导航系统工作的需要而选取的用于导航解算的参考坐标系。北轴指向地球北；东轴指向地球东；地轴垂直于地球表面并指向下。
在一种可能的实现方式中，可以利用以下公式，根据第二向量坐标以及经纬度，计算卫星在北东地坐标系下的第三向量坐标：
X3＝-sin(B·π/180)*cos(L·π/180)-sin(B·π/180)*sin(L·π/180)*X2+cos(B·π/180)*Z2；
Y3＝-sin(L·π/180)*Y2+cos(L·π/180)*X2；
Z3＝-(cos(B·π/180)*cos(L·π/180)*Y2+cos(B·π/180)*sin(L·π/180)*X2+sin(B·π/180)*Z2)；
其中，X2、Y2和Z2为第二向量坐标，即卫星在WGS84坐标系下第二向量坐标对应于三个坐标轴的坐标。L为卫星对应的经度、B为卫星对应的纬度、HD为卫星对应的大地高度；X3、Y3和Z3为卫星在北东地坐标系下的第三向量坐标，即卫星在北东地坐标系下第三向量坐标对应于三个坐标轴的坐标。
步骤S107：根据第三向量坐标，计算卫星在北东地坐标系下的欧拉角。
由于卫星搭载的雷达一般使用北东地坐标系，则将卫星的相关参数最终均转换到北东地坐标系下。在步骤S106计算得到卫星在北东地坐标系下的第三向量坐标，可以进一步求出卫星在北东地坐标系下的欧拉角。
在一种可能的实现方式中，可以利用以下公式，根据第三向量坐标，计算卫星在北东地坐标系下的欧拉角，包括：
roll2＝0；
Figure BDA0002728282510000101
yaw2＝(atan2(Y3,X3))·180/π；
其中，X3、Y3和Z3为第三向量坐标，yaw2为卫星在北东地坐标系下的欧拉角中的偏航角，pitch2为卫星在北东地坐标系下的欧拉角中的俯仰角，roll2为卫星在北东地坐标系下的欧拉角中的翻滚角。
步骤S108：根据第三向量坐标、卫星在北东地坐标系下的欧拉角以及卫星对应的雷达在北东地坐标系下的欧拉角，获得以雷达为起点，预设角度为圆周角度的多条射线。
根据卫星在北东地坐标系下的位置(即第三向量坐标)，卫星在北东地坐标系下的欧拉角，以及输入的雷达在北东地坐标系下的欧拉角等参数，以卫星雷达位置为发射点，预设角度为圆周角度，发射多条射线。本申请实施例对预设角度不进行限定，可以根据实际情况设定，例如设定为0.2度。
步骤S109：计算射线与地球的交点坐标。
进一步可以根据线段与椭圆交点的求解方式，求得仿真发射的射线与地球的相交点，得到交点坐标。可以理解的是，多个交点坐标可以围成一个区域。
步骤S110：将交点坐标对应的区域在地图上进行展示。
本申请实施例则将交点坐标对应的区域在地图上进行展示。参见图2所示，示出了卫星雷达波范围在地图上展示的效果图。在实际应用中，可以根据交点坐标，以QGIS为工具向二维地图画矢量图层。从而实现雷达波对地面投影的可见性。对于装备仿真、工业仿真领域中，雷达波在二维地图的可见性展示效果有重要意义，本申请实施例提供了一种相对简单的二维地图雷达波投影成像方法，降低了计算雷达波的复杂度。
另外，在交点坐标对应的区域越过二维地图边界，即交点坐标对应的区域经过经度180度时，还可以优化雷达波范围的形状。
在一种可能的实现方式中，如果交点坐标对应的区域经过经度180度，将交点坐标对应的区域沿经度180度进行切割，将切割后的区域在地图上进行展示。
参见图3所示，示出了交点坐标对应的区域经过经度180度时，卫星雷达波范围在地图上展示的效果图。通过切割交点坐标对应的区域，将切割后的两个区域在二维地图的两侧分别显示。
这样，本申请实施例可以根据卫星的实时位置、实时速度以及朝向参数等，计算出卫星在J2000坐标系下的欧拉角，并通过在各个坐标系下的向量坐标最终转换为在北东地坐标系下的欧拉角。同时，结合卫星挂载的雷达在北东地坐标系下的欧拉角，可以计算得到雷达波的发射方向，从而计算出雷达波与地球的交点，将雷达波与地球的交点对应的区域在地图上进行展示，从而直观对雷达波的扫描范围在地图上进行展示。另外，本申请还可以在雷达波范围穿过经度180度时，优化雷达波形状，在地图上准确展示雷达波的扫描范围。
基于上述方法实施例提供的实现卫星雷达波范围在地图上展示的方法，本申请实施例还提供了一种实现卫星雷达波范围在地图上展示的装置，下面将结合附图对该实现卫星雷达波范围在地图上展示的装置进行说明。
参见图4，该图为本申请实施例提供的一种实现卫星雷达波范围在地图上展示的装置的结构示意图。如图4所示，该实现卫星雷达波范围在地图上展示的装置可以包括：
获取单元401，用于获取卫星的轨道六根数参数，根据所述轨道六根数参数计算所述卫星在J2000坐标系下的实时位置以及实时速度；
第一计算单元402，用于根据所述实时位置、所述实时速度以及所述卫星的朝向参数，计算所述卫星在J2000坐标系下的欧拉角；
第一转换单元403，用于将所述卫星在J2000坐标系下的欧拉角转换为在J2000坐标系下的第一向量坐标；
第二转换单元404，用于将所述第一向量坐标转换为在WGS84坐标系下的第二向量坐标；
第二计算单元405，用于根据所述实时位置计算所述卫星对应的经纬度；
第三计算单元406，用于根据所述第二向量坐标以及所述经纬度，计算所述卫星在北东地坐标系下的第三向量坐标；
第四计算单元407，用于根据所述第三向量坐标，计算所述卫星在北东地坐标系下的欧拉角；
获得单元408，用于根据所述第三向量坐标、所述卫星在北东地坐标系下的欧拉角以及所述卫星对应的雷达在北东地坐标系下的欧拉角，获得以所述雷达为起点，预设角度为圆周角度的多条射线；
第五计算单元409，用于计算所述射线与地球的交点坐标；
展示单元410，用于将所述交点坐标对应的区域在地图上进行展示。
在一种可能的实现方式中，所述装置还包括：
切割单元，用于如果所述交点坐标对应的区域经过经度180度，将所述交点坐标对应的区域沿所述经度180度进行切割；
所述展示单元，具体用于将切割后的区域在地图上进行展示。
在一种可能的实现方式中，所述第一转换单元具体用于：
利用以下公式，将所述卫星在J2000坐标系下的欧拉角转换为在J2000坐标系下的第一向量坐标：
X1＝cos(yaw1·π/180)*cos(pitch1·π/180)；
Y1＝sin(yaw1·π/180)*cos(pitch1·π/180)；
Z1＝sin(pitch1·π/180)*(-1)；
其中，yaw1为所述卫星在J2000坐标系下的欧拉角中的偏航角，pitch1为所述卫星在J2000坐标系下的欧拉角中的俯仰角；X1、Y1和Z1为J2000坐标系下的第一向量坐标。
在一种可能的实现方式中，所述第二计算单元具体用于：
利用以下公式，根据所述实时位置计算所述卫星对应的经纬度：
L＝arctan(Y/X)；
Figure BDA0002728282510000131
HD＝Z/sinB-N(1-e2)；
其中，X、Y和Z为所述实时位置，N为法截弧的曲率半径，e为自然常数，L为所述卫星对应的经度、B为所述卫星对应的纬度、HD为所述卫星对应的大地高度。
在一种可能的实现方式中，所述第三计算单元具体用于：
利用以下公式，根据所述第二向量坐标以及所述经纬度，计算所述卫星在北东地坐标系下的第三向量坐标：
X3＝-sin(B·π/180)*cos(L·π/180)-sin(B·π/180)*sin(L·π/180)*X2+cos(B·π/180)*Z2；
Y3＝-sin(L·π/180)*Y2+cos(L·π/180)*X2；
Z3＝-(cos(B·π/180)*cos(L·π/180)*Y2+cos(B·π/180)*sin(L·π/180)*X2+sin(B·π/180)*Z2)；
其中，X2、Y2和Z2为所述第二向量坐标；L为所述卫星对应的经度、B为所述卫星对应的纬度、HD为所述卫星对应的大地高度；X3、Y3和Z3为所述卫星在北东地坐标系下的第三向量坐标。
在一种可能的实现方式中，所述第四计算单元具体用于：
利用以下公式，根据所述第三向量坐标，计算所述卫星在北东地坐标系下的欧拉角：
roll2＝0；
Figure BDA0002728282510000141
yaw2＝(atan2(Y3,X3))·180/π；
其中，X3、Y3和Z3为所述第三向量坐标，yaw2为所述卫星在北东地坐标系下的欧拉角中的偏航角，pitch2为所述卫星在北东地坐标系下的欧拉角中的俯仰角，roll2为所述卫星在北东地坐标系下的欧拉角中的翻滚角。
另外，本申请实施例还提供了一种实现卫星雷达波范围在地图上展示的设备，包括：存储器，处理器，及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时，实现如上述实施例所述的实现卫星雷达波范围在地图上展示的方法。
本申请实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质中存储有指令，当所述指令在终端设备上运行时，使得所述终端设备执行如上述实施例所述的实现卫星雷达波范围在地图上展示的方法。
由此可见，本申请实施例可以根据卫星的实时位置、实时速度以及朝向参数等，计算出卫星在J2000坐标系下的欧拉角，并通过在各个坐标系下的向量坐标最终转换为在北东地坐标系下的欧拉角。同时，结合卫星挂载的雷达在北东地坐标系下的欧拉角，可以计算得到雷达波的发射方向，从而计算出雷达波与地球的交点，将雷达波与地球的交点对应的区域在地图上进行展示，从而直观对雷达波的扫描范围在地图上进行展示。
需要说明的是，本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统或装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
应当理解，在本申请中，“至少一个(项)”是指一个或者多个，“多个”是指两个或两个以上。“和/或”，用于描述关联对象的关联关系，表示可以存在三种关系，例如，“A和/或B”可以表示：只存在A，只存在B以及同时存在A和B三种情况，其中A，B可以是单数或者复数。字符“/”一般表示前后关联对象是一种“或”的关系。“以下至少一项(个)”或其类似表达，是指这些项中的任意组合，包括单项(个)或复数项(个)的任意组合。例如，a，b或c中的至少一项(个)，可以表示：a，b，c，“a和b”，“a和c”，“b和c”，或“a和b和c”，其中a，b，c可以是单个，也可以是多个。
还需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本申请。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下，在其它实施例中实现。因此，本申请将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (10)

1.一种实现卫星雷达波范围在地图上展示的方法，其特征在于，所述方法包括：
获取卫星的轨道六根数参数，根据所述轨道六根数参数计算所述卫星在J2000坐标系下的实时位置以及实时速度；
根据所述实时位置、所述实时速度以及所述卫星的朝向参数，计算所述卫星在J2000坐标系下的欧拉角；
将所述卫星在J2000坐标系下的欧拉角转换为在J2000坐标系下的第一向量坐标；
将所述第一向量坐标转换为在WGS84坐标系下的第二向量坐标；
根据所述实时位置计算所述卫星对应的经纬度；
根据所述第二向量坐标以及所述经纬度，计算所述卫星在北东地坐标系下的第三向量坐标；
根据所述第三向量坐标，计算所述卫星在北东地坐标系下的欧拉角；
根据所述第三向量坐标、所述卫星在北东地坐标系下的欧拉角以及所述卫星对应的雷达在北东地坐标系下的欧拉角，获得以所述雷达为起点，预设角度为圆周角度的多条射线；
计算所述射线与地球的交点坐标；
将所述交点坐标对应的区域在地图上进行展示。
2.根据权利要求1所述的方法，其特征在于，所述方法还包括：
如果所述交点坐标对应的区域经过经度180度，将所述交点坐标对应的区域沿所述经度180度进行切割，将切割后的区域在地图上进行展示。
3.根据权利要求1所述的方法，其特征在于，所述将所述卫星在J2000坐标系下的欧拉角转换为在J2000坐标系下的第一向量坐标，包括：
利用以下公式，将所述卫星在J2000坐标系下的欧拉角转换为在J2000坐标系下的第一向量坐标：
X1＝cos(yaw1·π/180)*cos(pitch1·π/180)；
Y1＝sin(yaw1·π/180)*cos(pitch1·π/180)；
Z1＝sin(pitch1·π/180)*(-1)；
其中，yaw1为所述卫星在J2000坐标系下的欧拉角中的偏航角，pitch1为所述卫星在J2000坐标系下的欧拉角中的俯仰角；X1、Y1和Z1为J2000坐标系下的第一向量坐标。
4.根据权利要求1所述的方法，其特征在于，所述根据所述实时位置计算所述卫星对应的经纬度，包括：
利用以下公式，根据所述实时位置计算所述卫星对应的经纬度：
L＝arctan(Y/X)；
Figure FDA0002728282500000021
HD＝Z/sinB-N(1-e2)；
其中，X、Y和Z为所述实时位置，N为法截弧的曲率半径，e为自然常数，L为所述卫星对应的经度、B为所述卫星对应的纬度、HD为所述卫星对应的大地高度。
5.根据权利要求1所述的方法，其特征在于，所述根据所述第二向量坐标以及所述经纬度，计算所述卫星在北东地坐标系下的第三向量坐标，包括：
利用以下公式，根据所述第二向量坐标以及所述经纬度，计算所述卫星在北东地坐标系下的第三向量坐标：
X3＝-sin(B·π/180)*cos(L·π/180)-sin(B·π/180)*sin(L·π/180)*X2+cos(B·π/180)*Z2；
Y3＝-sin(L·π/180)*Y2+cos(L·π/180)*X2；
Z3＝-(cos(B·π/180)*cos(L·π/180)*Y2+cos(B·π/180)*sin(L·π/180)*X2+sin(B·π/180)*Z2)；
其中，X2、Y2和Z2为所述第二向量坐标；L为所述卫星对应的经度、B为所述卫星对应的纬度、HD为所述卫星对应的大地高度；X3、Y3和Z3为所述卫星在北东地坐标系下的第三向量坐标。
6.根据权利要求1所述的方法，其特征在于，所述根据所述第三向量坐标，计算所述卫星在北东地坐标系下的欧拉角，包括：
利用以下公式，根据所述第三向量坐标，计算所述卫星在北东地坐标系下的欧拉角：
roll2＝0；
Figure FDA0002728282500000031
yaw2＝(atan2(Y3，X3))·180/π；
其中，X3、Y3和Z3为所述第三向量坐标，yaw2为所述卫星在北东地坐标系下的欧拉角中的偏航角，pitch2为所述卫星在北东地坐标系下的欧拉角中的俯仰角，roll2为所述卫星在北东地坐标系下的欧拉角中的翻滚角。
7.一种实现卫星雷达波范围在地图上展示的装置，其特征在于，所述装置包括：
获取单元，用于获取卫星的轨道六根数参数，根据所述轨道六根数参数计算所述卫星在J2000坐标系下的实时位置以及实时速度；
第一计算单元，用于根据所述实时位置、所述实时速度以及所述卫星的朝向参数，计算所述卫星在J2000坐标系下的欧拉角；
第一转换单元，用于将所述卫星在J2000坐标系下的欧拉角转换为在J2000坐标系下的第一向量坐标；
第二转换单元，用于将所述第一向量坐标转换为在WGS84坐标系下的第二向量坐标；
第二计算单元，用于根据所述实时位置计算所述卫星对应的经纬度；
第三计算单元，用于根据所述第二向量坐标以及所述经纬度，计算所述卫星在北东地坐标系下的第三向量坐标；
第四计算单元，用于根据所述第三向量坐标，计算所述卫星在北东地坐标系下的欧拉角；
获得单元，用于根据所述第三向量坐标、所述卫星在北东地坐标系下的欧拉角以及所述卫星对应的雷达在北东地坐标系下的欧拉角，获得以所述雷达为起点，预设角度为圆周角度的多条射线；
第五计算单元，用于计算所述射线与地球的交点坐标；
展示单元，用于将所述交点坐标对应的区域在地图上进行展示。
8.根据权利要求7所述的装置，其特征在于，所述装置还包括：
切割单元，用于如果所述交点坐标对应的区域经过经度180度，将所述交点坐标对应的区域沿所述经度180度进行切割；
所述展示单元，具体用于将切割后的区域在地图上进行展示。
9.一种实现卫星雷达波范围在地图上展示的设备，其特征在于，包括：存储器，处理器，及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时，实现如权利要求1-6任一项所述的实现卫星雷达波范围在地图上展示的方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有指令，当所述指令在终端设备上运行时，使得所述终端设备执行如权利要求1-6任一项所述的实现卫星雷达波范围在地图上展示的方法。
CN202011109983.4A 2020-10-16 实现卫星雷达波范围在地图上展示的方法、装置及设备 Active CN112255599B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011109983.4A CN112255599B (zh) 2020-10-16 实现卫星雷达波范围在地图上展示的方法、装置及设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011109983.4A CN112255599B (zh) 2020-10-16 实现卫星雷达波范围在地图上展示的方法、装置及设备

Publications (2)

Publication Number Publication Date
CN112255599A true CN112255599A (zh) 2021-01-22
CN112255599B CN112255599B (zh) 2024-08-30

Family

ID=

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115186044A (zh) * 2022-06-27 2022-10-14 北京连山科技股份有限公司 一种基于卫星临空的地图标注方法
CN115270522A (zh) * 2022-09-27 2022-11-01 北京未尔锐创科技有限公司 一种基于wgs84坐标仿真跟踪目标设备的方法及装置

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6611226B1 (en) * 2000-04-20 2003-08-26 Hughes Electronics Corp Satellite surveillance system and method
CN102967851A (zh) * 2012-11-14 2013-03-13 电子科技大学 一种双基地sar的空间同步方法
RU2483280C1 (ru) * 2012-01-10 2013-05-27 Антон Владимирович Чернявец Навигационный комплекс
CN103809178A (zh) * 2014-01-17 2014-05-21 西安空间无线电技术研究所 一种地球同步轨道合成孔径雷达实现覆盖区连续观测方法
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
CN106226768A (zh) * 2016-08-09 2016-12-14 北京空间飞行器总体设计部 超高分辨率敏捷sar卫星滑动聚束模式系统参数设计方法
CN107505615A (zh) * 2017-08-22 2017-12-22 西安电子科技大学 星载sar滑动聚束成像模式星体姿态设计方法
US20180299557A1 (en) * 2017-04-17 2018-10-18 Baidu Online Network Technology (Beijing) Co., Ltd Method and apparatus for updating maps
CN111060077A (zh) * 2019-12-20 2020-04-24 彭耿 基于稀疏控制点的遥感卫星图像定位方法
CN111580055A (zh) * 2020-06-03 2020-08-25 北京世冠金洋科技发展有限公司 一种绘制雷达辐射范围的方法及装置
CN111591472A (zh) * 2020-05-15 2020-08-28 北京世冠金洋科技发展有限公司 一种调整卫星姿态的方法和相关装置

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6611226B1 (en) * 2000-04-20 2003-08-26 Hughes Electronics Corp Satellite surveillance system and method
RU2483280C1 (ru) * 2012-01-10 2013-05-27 Антон Владимирович Чернявец Навигационный комплекс
CN102967851A (zh) * 2012-11-14 2013-03-13 电子科技大学 一种双基地sar的空间同步方法
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
CN103809178A (zh) * 2014-01-17 2014-05-21 西安空间无线电技术研究所 一种地球同步轨道合成孔径雷达实现覆盖区连续观测方法
CN106226768A (zh) * 2016-08-09 2016-12-14 北京空间飞行器总体设计部 超高分辨率敏捷sar卫星滑动聚束模式系统参数设计方法
US20180299557A1 (en) * 2017-04-17 2018-10-18 Baidu Online Network Technology (Beijing) Co., Ltd Method and apparatus for updating maps
CN107505615A (zh) * 2017-08-22 2017-12-22 西安电子科技大学 星载sar滑动聚束成像模式星体姿态设计方法
CN111060077A (zh) * 2019-12-20 2020-04-24 彭耿 基于稀疏控制点的遥感卫星图像定位方法
CN111591472A (zh) * 2020-05-15 2020-08-28 北京世冠金洋科技发展有限公司 一种调整卫星姿态的方法和相关装置
CN111580055A (zh) * 2020-06-03 2020-08-25 北京世冠金洋科技发展有限公司 一种绘制雷达辐射范围的方法及装置

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
夏伟;李杨寰;梁福来;周智敏;: "旋翼无人机载SAR杠杆臂误差补偿算法", 计算机仿真, no. 06 *
尚学刚: "基于STK的临近空间飞行器雷达探测仿真研究", 中国优秀硕士学位论文全文数据库 信息科技辑, no. 03 *
朱爱萍: "高分卫星轨道及其对地覆盖的可视化研究", 中国优秀硕士学位论文全文数据库工程科技Ⅱ辑, no. 02, 15 February 2020 (2020-02-15) *
杨震;高辰;牛文龙;孟新;李大林;: "基于坐标系转换与QTM编码的天球面探测覆盖栅格化分析方法研究", 空间科学学报, no. 01, 31 December 2014 (2014-12-31) *
林雪原;何友;: "三星定位性能分析及其与SINS的融合方法", 测绘学报, no. 02 *
江烨: "多模式星载合成孔径雷达成像算法研究", 中国优秀硕士学位论文全文数据库信息科技辑, no. 02 *
程春泉;张继贤;黄国满;马晶;: "地心直角坐标系下的航摄光学影像严密定位模型", 武汉大学学报(信息科学版), no. 05, 30 July 2015 (2015-07-30) *
连雯卓: "海洋遥感卫星组网观测仿真系统设计", 中国优秀硕士学位论文全文数据库基础科学辑, no. 09, pages 24 - 27 *
钟巍;高悠然;: "一种卫星对地覆盖区域计算方法", 现代雷达, no. 07, 25 July 2020 (2020-07-25) *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115186044A (zh) * 2022-06-27 2022-10-14 北京连山科技股份有限公司 一种基于卫星临空的地图标注方法
CN115270522A (zh) * 2022-09-27 2022-11-01 北京未尔锐创科技有限公司 一种基于wgs84坐标仿真跟踪目标设备的方法及装置

Similar Documents

Publication Publication Date Title
JP2516257B2 (ja) 航空機の複数のレ―ダ観測から航空機の位置を表示する方法
CN102565797A (zh) 一种针对聚束模式星载sar图像的几何校正方法
CN105067120A (zh) 星载微波辐射计观测亮温的动态滤波重采样方法及装置
JP2843666B2 (ja) 飛翔体の姿勢表示方法
CN103644918A (zh) 卫星对月探测数据定位处理方法
CN112991535B (zh) 高度信息增强墨卡托地图的三维空间态势表征方法和装置
CN103344958B (zh) 基于星历数据的星载sar高阶多普勒参数估算方法
CN103727937A (zh) 一种基于星敏感器的舰船姿态确定方法
Li et al. Autonomous orbit and attitude determination for Earth satellites using images of regular-shaped ground objects
CN106918317A (zh) 一种扫描式卫星对地侧摆观测覆盖带确定方法及装置
CN106897962A (zh) 大视场空间对地观测图像的圆盘投影和拼接方法
Zeng et al. Pan-location mapping and localization for the in-situ science exploration of Zhurong Mars rover
CN108225306B (zh) 基于遥感卫星凝视姿态的星敏感器安装布局方法
CN112255599A (zh) 实现卫星雷达波范围在地图上展示的方法、装置及设备
CN106250684B (zh) 基于地固系数据的卫星过境时间快速计算方法
CN112255599B (zh) 实现卫星雷达波范围在地图上展示的方法、装置及设备
CN104391311B (zh) 基于gps广播数据的星上无源定位方法
CN114021068B (zh) 太阳同步圆轨道卫星受晒因子计算方法、装置及电子设备
US20150325045A1 (en) Map drawing device
Golubaev Main parameters of meteoroid motion during the fall of the Chelyabinsk meteorite shower on February 15, 2013
Zhang et al. A new method of single celestial-body sun positioning based on theory of mechanisms
CN109581374B (zh) 单侧sar卫星成像形态仿真、动态仿真方法及系统
CN110096721B (zh) 一种航天器对地面复杂区域目标的观测可见性判断方法
CN102519471A (zh) 基于三角函数拟合的成像型地球敏感器地球扁率补偿方法
Tsou Using GIS to obtain Celestial fix under the framework of an ECDIS System

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