CN112346028B - 分布式InSAR卫星方位向空间同步在轨测试方法及*** - Google Patents
分布式InSAR卫星方位向空间同步在轨测试方法及*** Download PDFInfo
- Publication number
- CN112346028B CN112346028B CN202011197192.1A CN202011197192A CN112346028B CN 112346028 B CN112346028 B CN 112346028B CN 202011197192 A CN202011197192 A CN 202011197192A CN 112346028 B CN112346028 B CN 112346028B
- Authority
- CN
- China
- Prior art keywords
- satellite
- azimuth
- orbit
- calculating
- threshold
- 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
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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供了一种分布式InSAR卫星方位向空间同步在轨测试方法及***,方法包括:步骤1:通过主星轨道及姿态信息计算方位向波束指向AP01;步骤2:计算方位向波束指向AP02;步骤3:判断|AP01‑AP02|是否超出门限;步骤4:通过辅星轨道及姿态信息计算方位向波束指向AP11;步骤5:计算多普勒中心AP12;步骤6:判断|AP11‑AP12|是否超出门限;步骤7:选取|AP01‑AP02|、|AP11‑AP12|同时处于门限范围内的测试样本对,则方位向空间同步测试结果为|AP12‑AP02|/AW×100%,AW为主星方位向波束宽度;步骤8:计算各组测试样本结果,对各组计算值取平均,进一步降低随机误差,并获取最终测试结果。本发明可有效获取分布式InSAR卫星在轨方位向空间同步测试结果。
Description
技术领域
本发明涉及航天***技术领域,具体地,涉及一种分布式InSAR卫星方位向空间同步在轨测试方法及***。
背景技术
干涉合成孔径雷达(InSAR)是获取高精度地面高程模型(DSM)的重要遥感手段。它利用沿垂直航向分布的两部SAR天线以不同的视角观测同一地区,并将获取的两幅复SAR图像进行干涉处理,求取主副雷达天线相位中心与目标之间的斜距差,进而获得观测区域的DSM。分布式卫星InSAR***将两部SAR安装在编队飞行的两颗卫星上同时对地观测,它可以克服重复航过InSAR面临的时间去相干及基线精度低等问题,可获取高精度的DSM。由于SAR发射天线和接收天线分置在不同的卫星平台上,要实现单星SAR成像以及InSAR干涉测量,收发天线波束必须能同时覆盖同一地面照射区域,实现双星SAR波束的空间同步,以保证辅星能够接收到足够的回波能量并与主星保持较好的相干性。
本发明针对最大相干法,提出了一种方位向空间同步指标的在轨测试方法。采用最大相干法的分布式InSAR卫星,主辅星SAR均以正侧视模块工作,SAR方位向波束宽度相同、波束指向均处于各自的全零多普勒面内,如图2所示。最大相干法下的方位向空间同步要求主辅星分别按照各自的姿态导引规律将波束照射在同一个垂直于卫星飞行方向的测绘带内,该模块下主辅星回波多普勒中心频率残余误差较小,且基本相同,保证了较高的多普勒去相干指标,对InSAR***获取高精度DSM 产品提供有力支持。
公开号为CN106054185的专利文献“一种基于分布式POS的机载双天线InSAR 基线计算方法”解决的是机载双天线InSAR基线的计算方法,而本发明中主要解决的是分布式InSAR卫星空间同步在轨运行的测试方法,两者的应用方向、适用范围、技术途径上均有明显差异。
基于卫星编队InSAR空间同步对***性能影响的分析,武汉大学学报(信息科学版),200710。主要差异:本专利主要解决的是分布式InSAR卫星空间同步在轨运行的测试方法,而该论文主要论证了编队条件下InSAR高程测量对空间同步的要求并分析了空间同步对InSAR***的高程测量精度和分辨率的影响,应用方向、适用范围、技术途径上均有明显差异。
干涉SAR卫星编队波束同步方法,中国空间科学技术,201005。主要差异:本专利主要解决的是分布式InSAR卫星空间同步在轨运行的测试方法,而该论文主要提出了通过卫星编队姿态导引完成波束同步的工程实现途径,应用方向、适用范围、技术途径上均有明显差异。
一种实用化的星机双基地SAR空间同步方法,电子与信息学报,200806。主要差异:本专利主要解决的是分布式InSAR卫星空间同步在轨运行的测试方法,而该论文主要提出卫星和飞机收发平台波束指向的空间同步实现途径,应用方向、适用范围、技术途径上均有明显差异。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种分布式InSAR卫星方位向空间同步在轨测试方法及***。
根据本发明的一个方面,提供了一种分布式InSAR卫星方位向空间同步在轨测试方法,包括如下步骤:
步骤1:通过主星轨道及姿态信息计算方位向波束指向AP01;
步骤2:利用观测热带雨林地区的回波数据计算多普勒中心,并计算方位向波束指向AP02;
步骤3:判断|AP01-AP02|是否超出门限,若处于门限值内可进行下步工作,若超门限值则该段数据不满足测试要求,更换样本;
步骤4:通过辅星轨道及姿态信息计算方位向波束指向AP11;
步骤5:利用观测热带雨林地区的回波数据计算多普勒中心AP12;
步骤6:判断|AP11-AP12|是否超出门限,若处于门限值内可进行下步工作,若超门限值则该段数据不满足测试要求,更换样本;
步骤7:选取|AP01-AP02|、|AP11-AP12|同时处于门限范围内的测试样本对,则方位向空间同步测试结果为|AP12-AP02|/AW×100%,AW为主星方位向波束宽度;
步骤8:计算各组测试样本结果,对各组计算值取平均,进一步降低随机误差,并获取最终测试结果。
AP01=θcosα-ψsinα。
优选地,所述步骤2中,AP02的计算方法是,设多普勒中心频率为fd,卫星速度为v,SAR载波频率为fc,光速为c,则
优选地,所述步骤2中,分布式卫星测试时过热带雨林地区需尽量一致,要求多星测试区域相差不超过10%,否则测试精度下降严重;有效测试样本对越多测试精度越高,本专利方法要求远视角下有效样本Sfar、中视角下有效样本Smiddle、近视角下有效样本Snear至少各2组才能保证本测试方法精度。
优选地,所述步骤3中,判断|AP1-AP0|是否超出门限,主辅星对同一雨林地区的回波数据同时处于门限内则为有效样本,若主辅星任一数据超门限则为无效样本;
门限值的设置与SAR工作频段、卫星姿态导引规律和产品生产水平因素相关,可由使用者根据实际情况自行选取。
根据本发明的另一个方面,提供了一种分布式InSAR卫星方位向空间同步在轨测试***,包括:
模块M1:通过主星轨道及姿态信息计算方位向波束指向AP01;
模块M2:利用观测热带雨林地区的回波数据计算多普勒中心,并计算方位向波束指向AP02;
模块M3:判断|AP01-AP02|是否超出门限,若处于门限值内可进行下步工作,若超门限值则该段数据不满足测试要求,更换样本;
模块M4:通过辅星轨道及姿态信息计算方位向波束指向AP11;
模块M5:利用观测热带雨林地区的回波数据计算多普勒中心AP12;
模块M6:判断|AP11-AP12|是否超出门限,若处于门限值内可进行下步工作,若超门限值则该段数据不满足测试要求,更换样本;
模块M7:选取|AP01-AP02|、|AP11-AP12|同时处于门限范围内的测试样本对,则方位向空间同步测试结果为|AP12-AP02|/AW×100%,AW为主星方位向波束宽度;
模块M8:计算各组测试样本结果,对各组计算值取平均,进一步降低随机误差,并获取最终测试结果。
AP01=θcosα-ψsinα。
优选地,所述模块M2中,AP02的计算方法是,设多普勒中心频率为fd,卫星速度为v,SAR载波频率为fc,光速为c,则
优选地,所述模块M2中,分布式卫星测试时过热带雨林地区需尽量一致,要求多星测试区域相差不超过10%,否则测试精度下降严重;有效测试样本对越多测试精度越高,本专利方法要求远视角下有效样本Sfar、中视角下有效样本Smiddle、近视角下有效样本Snear至少各2组才能保证本测试方法精度。
优选地,所述模块M3中,判断|AP1-AP0|是否超出门限,主辅星对同一雨林地区的回波数据同时处于门限内则为有效样本,若主辅星任一数据超门限则为无效样本;
门限值的设置与SAR工作频段、卫星姿态导引规律和产品生产水平因素相关,可由使用者根据实际情况自行选取。
与现有技术相比,本发明具有如下的有益效果:
1、本发明可用于分布式InSAR卫星的在轨测试工作,提出一种最大相干法下方位向空间同步指标的在轨测试手段;
2、本发明可有效获取最大相干法分布式InSAR卫星在轨方位向空间同步测试结果;
3、本发明具备在轨测试数据的有效性检查功能,可以通过设置门限筛选出由于地形地貌变化造成的无效数据,提升了测试算法的有效性及可靠性。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为分布式InSAR卫星方位向空间同步在轨测试方法的处理步骤流程图;
图2为方位向空间同步示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
本发明提供了一种分布式InSAR卫星方位向空间同步在轨测试方法,如图1所示,在本实施例中,具体实现步骤如下:
(1)通过主星轨道及姿态信息计算方位向波束指向AP01;
(2)利用观测热带雨林地区的回波数据计算多普勒中心,并计算方位向波束指向AP02;
(3)判断|AP01-AP02|是否超出门限,若处于门限值内可进行下步工作,若超门限值则该段数据不满足测试要求,更换样本;
(4)通过辅星轨道及姿态信息计算方位向波束指向AP11;
(5)利用观测热带雨林地区的回波数据计算多普勒中心AP12;
(6)判断|AP11-AP12|是否超出门限,若处于门限值内可进行下步工作,若超门限值则该段数据不满足测试要求,更换样本;
(7)选取|AP01-AP02|、|AP11-AP12|同时处于门限范围内的测试样本对,则方位向空间同步测试结果为|AP12-AP02|/AW×100%,AW为主星方位向波束宽度。
(8)计算各组测试样本结果,对各组计算值取平均,进一步降低随机误差,并获取最终测试结果。
需要特殊说明的是,门限值的设置与SAR工作频段、卫星姿态导引规律、产品生产水平等诸多因素相关,可由使用者根据实际情况自行选取。
以德国TanDEM-X***为例,主要信息如表1所示:
表1 TanDEM-X卫星***信息
以巴西境内赤道附近的亚马逊热带雨林地区为例,指向测量精度(即AP01、AP11精度)为0.01°;多普勒中心估计精度为30Hz,折算至AP02、AP12的计算精度约0.0025°;因此门限值可设置为0.015°,本发明应用在德国TanDEM-X***上,单组测试样本的测试精度可达(0.0025°+0.0025°)/0.33°×100%=1.5%,通过对多组有效样本的测试结果取平均,可进一步降低随机误差,测试精度有望进一步提高。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (10)
1.一种分布式InSAR卫星方位向空间同步在轨测试方法,其特征在于,包括:
步骤1:通过主星轨道及姿态信息计算方位向波束指向AP01;
步骤2:利用观测热带雨林地区的回波数据计算多普勒中心,并利用主星观测的数据计算方位向波束指向AP02;
步骤3:判断|AP01-AP02|是否超出门限,若处于门限值内可进行下步工作,若超门限值则该回波数据不满足测试要求,更换样本;
步骤4:通过辅星轨道及姿态信息计算方位向波束指向AP11;
步骤5:利用观测热带雨林地区的回波数据计算多普勒中心,并利用辅助观测的数据计算方位向波束指向AP12;
步骤6:判断|AP11-AP12|是否超出门限,若处于门限值内可进行下步工作,若超门限值则该回波数据不满足测试要求,更换样本;
步骤7:选取|AP01-AP02|、|AP11-AP12|同时处于门限范围内的测试样本对,则方位向空间同步测试结果为|AP12+AP02|/AW×100%,AW为主星方位向波束宽度;
步骤8:计算各组测试样本结果,对各组计算值取平均,进一步降低随机误差,并获取最终测试结果;
所述步骤2中,分布式卫星测试时过热带雨林地区一致,要求多星测试区域相差不超过10%。
4.根据权利要求3所述的分布式InSAR卫星方位向空间同步在轨测试方法,其特征在于,有效测试样本对越多测试精度越高,要求远视角下有效样本Sfar、中视角下有效样本Smiddle、近视角下有效样本Snear至少各2组才能保证测试精度。
5.根据权利要求4所述的分布式InSAR卫星方位向空间同步在轨测试方法,其特征在于,所述步骤3中,判断|AP1-AP0|是否超出门限,主星对同一雨林地区的回波数据同时处于门限内则为有效样本,若主星任一数据超门限则为无效样本;门限值的设置与SAR工作频段、卫星姿态导引规律和产品生产水平因素相关,可由使用者根据实际情况自行选取。
6.一种分布式InSAR卫星方位向空间同步在轨测试***,其特征在于,包括:
模块M1:通过主星轨道及姿态信息计算方位向波束指向AP01;
模块M2:利用观测热带雨林地区的回波数据计算多普勒中心,并利用主星观测的数据计算方位向波束指向AP02;
模块M3:判断|AP01-AP02|是否超出门限,若处于门限值内可进行下步工作,若超门限值则该回波数据不满足测试要求,更换样本;
模块M4:通过辅星轨道及姿态信息计算方位向波束指向AP11;
模块M5:利用观测热带雨林地区的回波数据计算多普勒中心,并利用辅助观测的数据计算方位向波束指向AP12;
模块M6:判断|AP11-AP12|是否超出门限,若处于门限值内可进行下步工作,若超门限值则该回波数据不满足测试要求,更换样本;
模块M7:选取|AP01-AP02|、|AP11-AP12|同时处于门限范围内的测试样本对,则方位向空间同步测试结果为|AP12+AP02|/AW×100%,AW为主星方位向波束宽度;
模块M8:计算各组测试样本结果,对各组计算值取平均,进一步降低随机误差,并获取最终测试结果;
所述模块M2中,分布式卫星测试时过热带雨林地区一致,要求多星测试区域相差不超过10%。
9.根据权利要求8所述的分布式InSAR卫星方位向空间同步在轨测试***,其特征在于,有效测试样本对越多测试精度越高,要求远视角下有效样本Sfar、中视角下有效样本Smiddle、近视角下有效样本Snear至少各2组才能保证测试精度。
10.根据权利要求9所述的分布式InSAR卫星方位向空间同步在轨测试***,其特征在于,所述模块M3中,判断|AP1-AP0|是否超出门限,主辅星对同一雨林地区的回波数据同时处于门限内则为有效样本,若主辅星任一数据超门限则为无效样本;门限值的设置与SAR工作频段、卫星姿态导引规律和产品生产水平因素相关,可由使用者根据实际情况自行选取。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011197192.1A CN112346028B (zh) | 2020-10-30 | 2020-10-30 | 分布式InSAR卫星方位向空间同步在轨测试方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011197192.1A CN112346028B (zh) | 2020-10-30 | 2020-10-30 | 分布式InSAR卫星方位向空间同步在轨测试方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112346028A CN112346028A (zh) | 2021-02-09 |
CN112346028B true CN112346028B (zh) | 2022-10-25 |
Family
ID=74355251
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011197192.1A Active CN112346028B (zh) | 2020-10-30 | 2020-10-30 | 分布式InSAR卫星方位向空间同步在轨测试方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112346028B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116859344B (zh) * | 2023-08-28 | 2023-11-03 | 中国电子科技集团公司第十四研究所 | 一种面向相干最优的能谱自适应分布式InSAR空间同步方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101464512A (zh) * | 2009-01-21 | 2009-06-24 | 电子科技大学 | 一种星机双基地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 |
CN110531328A (zh) * | 2019-06-25 | 2019-12-03 | 西安空间无线电技术研究所 | 一种星载环扫雷达波束中心指向误差估计方法 |
CN110907932A (zh) * | 2019-11-26 | 2020-03-24 | 上海卫星工程研究所 | 分布式InSAR卫星测高精度影响因素分析方法及*** |
CN111381256A (zh) * | 2020-03-10 | 2020-07-07 | 上海卫星工程研究所 | 主动遥感卫星天线相位中心偏移误差计算的方法和*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6825806B2 (en) * | 2002-06-03 | 2004-11-30 | The Boeing Company | Satellite methods and structures for improved antenna pointing and wide field-of-view attitude acquisition |
-
2020
- 2020-10-30 CN CN202011197192.1A patent/CN112346028B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101464512A (zh) * | 2009-01-21 | 2009-06-24 | 电子科技大学 | 一种星机双基地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 |
CN110531328A (zh) * | 2019-06-25 | 2019-12-03 | 西安空间无线电技术研究所 | 一种星载环扫雷达波束中心指向误差估计方法 |
CN110907932A (zh) * | 2019-11-26 | 2020-03-24 | 上海卫星工程研究所 | 分布式InSAR卫星测高精度影响因素分析方法及*** |
CN111381256A (zh) * | 2020-03-10 | 2020-07-07 | 上海卫星工程研究所 | 主动遥感卫星天线相位中心偏移误差计算的方法和*** |
Non-Patent Citations (5)
Title |
---|
InSAR编队卫星波束同步姿态策略分析;槐超等;《上海航天》;20150425(第02期);全文 * |
分布式SAR卫星偏航控制及精度要求;陈筠力等;《航天器工程》;20080115(第01期);全文 * |
分布式星载SAR***时间同步和波束同步误差分析;李航舰等;《雷达学报》;20170607(第02期);全文 * |
星载InSAR编队波束同步策略;贺东雷等;《宇航学报》;20090930(第05期);全文 * |
最大相干波束同步方法及其性能验证研究;张永俊等;《电子学报》;20120815(第08期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112346028A (zh) | 2021-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rabus et al. | The shuttle radar topography mission—a new class of digital elevation models acquired by spaceborne radar | |
EP2535735B1 (en) | Acquisition of sar images for computing a height or a digital elevation model by interferometric processing | |
Tian et al. | Review on interferometric ISAR 3D imaging: concept, technology and experiment | |
CN102169173B (zh) | 一种倾斜地球同步轨道合成孔径雷达模糊度分析方法 | |
US8791859B2 (en) | High precision radio frequency direction finding system | |
CN105068065A (zh) | 星载激光测高仪在轨检校方法及*** | |
CN110823191B (zh) | 混合基线双天线斜视干涉sar洋流测量性能确定方法及*** | |
Alvarez-Salazar et al. | Mission design for NISAR repeat-pass Interferometric SAR | |
CN113281747A (zh) | 一种基于sar成像的卫星 | |
KR102151362B1 (ko) | 극좌표변환을 이용한 항공기기반 영상복원장치 및 이를 이용한 영상복원방법 | |
CN112346028B (zh) | 分布式InSAR卫星方位向空间同步在轨测试方法及*** | |
Poe et al. | Geolocation error analysis of the special sensor microwave imager/sounder | |
Berardino et al. | On the Time-Domain Airborne SAR Focusing in the Presence of Strong Azimuth Variations of the Squint Angle | |
US11018705B1 (en) | Interference mitigation, target detection, location and measurement using separable waveforms transmitted from spatially separated antennas | |
Fasano et al. | Formation geometries for multistatic SAR tomography | |
CN112327262B (zh) | 分布式InSAR卫星SAR波束指向一致性在轨校准方法及*** | |
Liu et al. | A necessary model to quantify the scanning loss effect in spaceborne iGNSS-R ocean altimetry | |
RU2499279C1 (ru) | Способ оценки высоты полета летательного аппарата по радиолокационным изображениям земной поверхности | |
Kohlhase et al. | Interferometric baseline performance estimations for multistatic synthetic aperture radar configurations derived from GRACE GPS observations | |
Zhao et al. | Error analysis for the baseline estimation and calibration of distributed InSAR satellites | |
Qiu et al. | MEO-SAR in-orbit Elevation Antenna Pattern Determination Using Nano Calibration Satellite | |
Lim et al. | Phase statistical model and correction in imagery of ground based synthetic aperture radar (GBSAR) for land deformation monitoring | |
Hensley et al. | Lunar topographic mapping using a new high resolution mode for the GSSR radar | |
Wei et al. | Geolocation analysis for near-space synthetic aperture radar based on Rang-Doppler model | |
Buck et al. | A status update of investigations into the Wavemill concept |
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 |