CN107765244A - 基于机载双天线InSAR基线测量方法和装置 - Google Patents
基于机载双天线InSAR基线测量方法和装置 Download PDFInfo
- Publication number
- CN107765244A CN107765244A CN201710837745.7A CN201710837745A CN107765244A CN 107765244 A CN107765244 A CN 107765244A CN 201710837745 A CN201710837745 A CN 201710837745A CN 107765244 A CN107765244 A CN 107765244A
- Authority
- CN
- China
- Prior art keywords
- phase center
- auxiliary antenna
- calculating
- deformation
- antenna
- 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
Links
- 238000000691 measurement method Methods 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 37
- 238000005259 measurement Methods 0.000 claims abstract description 27
- 230000033001 locomotion Effects 0.000 claims abstract description 25
- 238000004804 winding Methods 0.000 claims description 45
- 238000004364 calculation method Methods 0.000 claims description 28
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000009466 transformation Effects 0.000 claims description 12
- 238000005452 bending Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 7
- 239000002243 precursor Substances 0.000 claims description 7
- 239000000463 material Substances 0.000 claims description 3
- 230000009977 dual effect Effects 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 10
- 238000003384 imaging method Methods 0.000 abstract description 8
- 238000012546 transfer Methods 0.000 abstract description 4
- 230000007812 deficiency Effects 0.000 abstract 1
- 230000008859 change Effects 0.000 description 9
- 230000007547 defect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
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
-
- 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/9023—SAR image post-processing techniques combined with interferometric techniques
-
- 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)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种基于机载双天线InSAR基线测量方法，对初始信息进行标校；基于标校的初始信息计算分布式POS相对角运动，获取辅天线相位中心处的绕曲变形角；通过绕曲变形角计算辅天线相位中心处的绕曲变形位置矢量；通过计算得到的辅天线相位中心处绕曲变形位置矢量，计算形变后的柔性基线。该方法克服了由于柔性基线动态变化过程难以精确建模而导致直接利用主POS至子IMU传递对准得到的位置计算基线产生较大误差的不足，且该方法具有计算精度高的特点，提高了机载双天线InSAR基线测量精度，提高了InSAR成像的精度。本发明还公开了一种基于机载双天线InSAR基线测量装置。
Description
技术领域
本发明涉及航空航天技术领域，特别是涉及一种基于机载双天线InSAR基线测量方法和装置。
背景技术
干涉合成孔径雷达(Interferometric Synthetic Aperture Radar，InSAR)技术通过将合成孔径雷达与微波干涉相结合，以不同视角下两幅或者多幅雷达复图像所形成的干涉相位为信息源进行反演而得到地表三维地形及其变化信息。双天线InSAR由主天线和辅天线组成，主辅天线相位中心之间的空间位移矢量(基线)的测量精度是影响InSAR干涉成像精度的关键因素之一，必须对其进行精确测量。
位置和姿态测量系统(Position and Orientation System，POS)是目前获取SAR天线运动参数的主要手段，主要由惯性测量单元(Inertial Measurement Unit，IMU)、全球定位系统(Global Positioning System，GPS)、POS导航计算机和后处理软件组成，IMU通常与SAR天线固定连接，用于测量三维线加速度和三维角速度，通过捷联解算和滤波算法将惯性量测信息同GPS量测信息进行融合，从而连续、实时地获取天线相位中心的位置、速度和姿态等全面的运动参数。对于机载双天线InSAR系统而言，由于工作过程中载机机体(包括机翼)的挠曲变形(尤其辅天线相位中心处存在较大变形)，天线相位中心之间的基线存在柔性变化，使得单个POS系统无法实现基线的准确测量。分布式POS由一个主POS和多个IMU(子IMU)组成，主POS由主IMU和GPS组成且与InSAR主天线固定连接，安装于机腹位置。主POS的主IMU和GPS通过惯性/卫星组合导航实现主天线相位中心运动信息的测量；子IMU与InSAR辅天线固定连接，并通过主POS至子IMU传递对准实现辅天线相位中心运动信息的测量。但是，由于主辅天线之间基线的柔性变化过程难以精确描述，使得主POS的位置和速度信息无法准确传递到子IMU，导致基线的测量误差，进而影响InSAR干涉成像精度。
发明内容
基于此，有必要针对传统技术存在的问题，提供一种基于机载双天线InSAR基线测量方法和装置，能够解决双天线InSAR的主辅天线之间基线随时间变化，难以对其精确建模而导致基线测量误差的问题，。本公开克服了由于柔性基线动态变化过程难以精确建模而导致直接利用主POS至子IMU传递对准得到的位置计算基线产生较大误差的不足。本公开具有计算精度高的特点，提高了机载双天线InSAR基线测量精度，提高了InSAR成像的精度。
第一方面，本发明实施例提供了一种基于机载双天线InSAR基线测量方法，所述方法包括：对初始信息进行标校；基于标校的初始信息计算分布式POS相对角运动，获取辅天线相位中心处的绕曲变形角；通过绕曲变形角计算辅天线相位中心处的绕曲变形位置矢量；通过计算得到的辅天线相位中心处绕曲变形位置矢量，计算形变后的柔性基线。
在其中一个实施例中，所述对初始信息进行标校包括：利用激光全站仪测量点坐标的方式获得测量机载双天线InSAR基线所需要的初始标校信息；其中，包括主POS坐标系与悬臂梁坐标系oBxByBzB之间坐标转换矩阵主POS坐标系与柔性基线形变前子IMU坐标系之间坐标转换矩阵悬臂梁坐标系oBxByBzB的oBxB轴方向的单位矢量eBx在主POS坐标系下的投影eBx＝[ex,ey,ez]T，ex、ey和ez分别表示eBx在轴、轴和轴的投影，形变前柔性基线长度r，形变前柔性基线r在主POS坐标系下的投影rx、ry和rz分别表示r在轴、轴和轴的投影。
在其中一个实施例中，所述基于标校的初始信息计算分布式POS相对角运动，获取辅天线相位中心处的绕曲变形角包括：利用标校的所述初始信息以及所述分布式POS的主POS和子IMU输出的高精度姿态矩阵，计算主POS和子IMU之间的相对角运动，并计算辅天线相位中心处的挠曲变形角σ。
在其中一个实施例中，所述利用标校的所述初始信息以及所述分布式POS的主POS和子IMU输出的高精度姿态矩阵，计算主POS和子IMU之间的相对角运动，并计算辅天线相位中心处的挠曲变形角包括：利用t时刻主POS和子IMU输出的姿态矩阵以及标校获得的主POS坐标系与柔性基线形变前子IMU坐标系之间坐标转换矩阵计算形变前子IMU坐标系与形变后子IMU坐标系之间的坐标转换矩阵
其中，公式为：
其中，为主POS在t时刻输出的姿态矩阵，为子IMU在t时刻输出的姿态矩阵，Tij为矩阵的第i行、第j列元素(i,j＝1,2,3)；
利用计算t时刻辅天线相位中心处的挠曲变形角σ；
其中，公式为：
其中，σ的取值范围满足且当ex(T32-T23)<0时，σ<0，当ex(T32-T23)>0时，σ>0。
在其中一个实施例中，所述通过绕曲变形角计算辅天线相位中心处的绕曲变形位置矢量包括：利用标校的所述初始信息以及获得的所述辅天线相位中心处的挠曲变形角σ，通过将柔性基线形变问题等效为悬臂梁自由端大挠度弯曲问题，根据工程力学理论建立辅天线相位中心处挠曲变形角σ与挠曲变形位置矢量△r之间的解析关系，计算得到辅天线相位中心处的挠曲变形位置矢量△r。
在其中一个实施例中，所述通过将柔性基线形变问题等效为悬臂梁自由端大挠度弯曲问题，根据工程力学理论建立辅天线相位中心处挠曲变形角σ与挠曲变形位置矢量△r之间的解析关系，计算得到辅天线相位中心处的挠曲变形位置矢量包括：建立t时刻辅天线相位中心处的挠曲变形角σ与挠曲变形位置矢量△r之间解析关系式；
其中，公式为：
其中，△rB为挠曲变形位置矢量△r在悬臂梁坐标系oBxByBzB下的矢量表示， 和分别为挠曲形变△r在悬臂梁坐标系oBxB轴、oByB轴和oBzB轴的投影，EI为柔性基线的弯曲刚度，根据柔性基线结构和材料确定，为常数，P为辅天线相位中心处的集中载荷，子IMU和辅天线的重量和，E(μ)为第二类完全椭圆积分，E(φ,μ)为第二类椭圆积分，F(μ)为第一类完全椭圆积分，F(φ,μ)为第一类完全椭圆积分；再将所述辅天线相位中心处的挠曲变形角σ与所述挠曲变形位置矢量△r之间解析关系式进行近似简化操作；
其中，公式为：
其中，α为指数系数，其取值为
在其中一个实施例中，所述通过计算得到的辅天线相位中心处绕曲变形位置矢量，计算形变后的柔性基线包括：通过计算得到的所述辅天线相位中心处挠曲变形位置矢量△rB，计算导航坐标系下形变后的柔性基线r′n(t)；
其中，公式为：
第二方面，本发明实施例提供了一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现上述第一方面的基于机载双天线InSAR基线测量方法。
第三方面，本发明实施例提供了一种包含指令的计算机程序产品，当该计算机程序产品在计算机上运行时，使得计算机执行上述第一方面所述的方法。
第四方面，本发明实施例还提供了一种基于机载双天线InSAR基线测量装置，所述装置包括：标校模块，对初始信息进行标校；获取模块，用于基于标校的初始信息计算分布式POS相对角运动，获取辅天线相位中心处的绕曲变形角；第一计算模块，用于通过绕曲变形角计算辅天线相位中心处的绕曲变形位置矢量；第二计算模块，用于通过计算得到的辅天线相位中心处绕曲变形位置矢量，计算形变后的柔性基线。
本发明提供的一种基于机载双天线InSAR基线测量方法和装置，对初始信息进行标校；基于标校的初始信息计算分布式POS相对角运动，获取辅天线相位中心处的绕曲变形角；通过绕曲变形角计算辅天线相位中心处的绕曲变形位置矢量；通过计算得到的辅天线相位中心处绕曲变形位置矢量，计算形变后的柔性基线。该方法克服了由于柔性基线动态变化过程难以精确建模而导致直接利用主POS至子IMU传递对准得到的位置计算基线产生较大误差的不足。本公开具有计算精度高的特点，提高了机载双天线InSAR基线测量精度，提高了InSAR成像的精度。
附图说明
图1为本发明一个实施例中的一种基于机载双天线InSAR基线测量方法的流程示意图；
图2为本发明一个实施例中的一种基于机载双天线InSAR基线测量方法中的柔性基线形变及坐标系示意图；以及
图3为本发明一个实施例中的一种基于机载双天线InSAR基线测量装置的结构示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例对本发明基于机载双天线InSAR基线测量方法和装置进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
如图1所示，为一个实施例中的一种基于机载双天线InSAR基线测量方法的流程示意图。具体包括以下步骤：
步骤102，对初始信息进行标校。本实施例中，对初始信息进行标校包括：利用激光全站仪测量点坐标的方式获得测量机载双天线InSAR基线所需要的初始标校信息。
其中，包括主POS坐标系与悬臂梁坐标系oBxByBzB之间坐标转换矩阵主POS坐标系与柔性基线形变前子IMU坐标系之间坐标转换矩阵悬臂梁坐标系oBxByBzB的oBxB轴方向的单位矢量eBx在主POS坐标系下的投影eBx＝[ex,ey,ez]T，ex、ey和ez分别表示eBx在轴、轴和轴的投影，形变前柔性基线长度r，形变前柔性基线r在主POS坐标系下的投影rx、ry和rz分别表示r在轴、轴和轴的投影。
步骤104，基于标校的初始信息计算分布式POS相对角运动，获取辅天线相位中心处的绕曲变形角。
本实施例中，基于标校的初始信息计算分布式POS相对角运动，获取辅天线相位中心处的绕曲变形角包括：利用标校的所述初始信息以及分布式POS的主POS和子IMU输出的高精度姿态矩阵，计算主POS和子IMU之间的相对角运动，并计算辅天线相位中心处的挠曲变形角σ。
进一步地，利用标校的初始信息以及分布式POS的主POS和子IMU输出的高精度姿态矩阵，计算主POS和子IMU之间的相对角运动，并计算辅天线相位中心处的挠曲变形角包括：
利用t时刻主POS和子IMU输出的姿态矩阵以及标校获得的主POS坐标系与柔性基线形变前子IMU坐标系之间坐标转换矩阵计算形变前子IMU坐标系与形变后子IMU坐标系之间的坐标转换矩阵
其中，公式为：
其中，为主POS在t时刻输出的姿态矩阵，为子IMU在t时刻输出的姿态矩阵，Tij为矩阵的第i行、第j列元素(i,j＝1,2,3)；
利用计算t时刻辅天线相位中心处的挠曲变形角σ；
其中，公式为：
其中，σ的取值范围满足且当ex(T32-T23)<0时，σ<0，当ex(T32-T23)>0时，σ>0。
步骤106，通过绕曲变形角计算辅天线相位中心处的绕曲变形位置矢量。
本实施例中，通过绕曲变形角计算辅天线相位中心处的绕曲变形位置矢量包括：利用标校的初始信息以及获得的辅天线相位中心处的挠曲变形角σ，通过将柔性基线形变问题等效为悬臂梁自由端大挠度弯曲问题，根据工程力学理论建立辅天线相位中心处挠曲变形角σ与挠曲变形位置矢量△r之间的解析关系，计算得到辅天线相位中心处的挠曲变形位置矢量△r。
进一步地，通过将柔性基线形变问题等效为悬臂梁自由端大挠度弯曲问题，根据工程力学理论建立辅天线相位中心处挠曲变形角σ与挠曲变形位置矢量△r之间的解析关系，计算得到辅天线相位中心处的挠曲变形位置矢量包括：
建立t时刻辅天线相位中心处的挠曲变形角σ与挠曲变形位置矢量△r之间解析关系式；
其中，公式为：
其中，△rB为挠曲变形位置矢量△r在悬臂梁坐标系oBxByBzB下的矢量表示， 和分别为挠曲形变△r在悬臂梁坐标系oBxB轴、oByB轴和oBzB轴的投影，EI为柔性基线的弯曲刚度，根据柔性基线结构和材料确定，为常数，P为辅天线相位中心处的集中载荷，子IMU和辅天线的重量和，E(μ)为第二类完全椭圆积分，E(φ,μ)为第二类椭圆积分，F(μ)为第一类完全椭圆积分，F(φ,μ)为第一类完全椭圆积分；
再将辅天线相位中心处的挠曲变形角σ与所述挠曲变形位置矢量△r之间解析关系式进行近似简化操作；
其中，公式为：
其中，α为指数系数，其取值为
步骤108，通过计算得到的辅天线相位中心处绕曲变形位置矢量，计算形变后的柔性基线。
本实施例中，通过计算得到的辅天线相位中心处绕曲变形位置矢量，计算形变后的柔性基线包括：通过计算得到的辅天线相位中心处挠曲变形位置矢量△rB，计算导航坐标系下形变后的柔性基线r′n(t)；
其中，公式为：
本发明提供的一种基于机载双天线InSAR基线测量方法，对初始信息进行标校；基于标校的初始信息计算分布式POS相对角运动，获取辅天线相位中心处的绕曲变形角；通过绕曲变形角计算辅天线相位中心处的绕曲变形位置矢量；通过计算得到的辅天线相位中心处绕曲变形位置矢量，计算形变后的柔性基线。该方法克服了由于柔性基线动态变化过程难以精确建模而导致直接利用主POS至子IMU传递对准得到的位置计算基线产生较大误差的不足，且该方法具有计算精度高的特点，提高了机载双天线InSAR基线测量精度，提高了InSAR成像的精度。
为了更清楚的理解与应用本公开提出的基于机载双天线InSAR基线测量方法，进行以下示例。需要说明的是，本公开的保护范围不局限以下示例。
如图2所示，为一个实施例中的一种基于机载双天线InSAR基线测量方法中的柔性基线形变及坐标系示意图。其中，主系统包括主POS和InSAR主天线，子系统包括子IMU和InSAR辅天线。oBxByBzB为悬臂梁坐标系，oB点位于梁的约束端中心位置，oBxB轴沿悬臂梁延伸方向，oBzB轴垂直于初始梁平面，oByB轴与oBxB轴、oBzB轴构成右手坐标系。为主POS坐标系，为柔性基线形变前子IMU坐标系，为柔性基线形变后子IMU坐标系。柔性基线形变前表示为r，柔性基线形变后表示为r′，辅天线相位中心处柔性基线的挠曲形变位置矢量为△r，σ为辅天线相位中心处的挠曲变形角。
综上所述，结合图1和图2，本公开的一种基于机载双天线InSAR基线测量方法的原理：针对机载双天线InSAR干涉天线相位中心之间的柔性基线难以精确建模，使得直接利用分布式POS的主POS和子IMU位置信息计算基线存在较大误差的问题，本发明首先利用激光全站仪测量点坐标的方式获得测量机载双天线InSAR基线所需要的初始标校信息；其次，利用主POS和子IMU的高精度姿态信息计算主POS和子IMU的相对角运动，进而得到辅天线位置处的挠曲变形角；再次，将柔性基线形变问题等效为悬臂梁自由端大挠度弯曲问题，得到辅天线相位中心处挠曲变形角与挠曲变形位置矢量之间解析关系，进而计算辅天线相位中心处的挠曲变形位置矢量；最后，利用得到的辅天线相位中心处的柔性基线挠曲变形位置矢量计算形变后的柔性基线。
进一步地，本公开提出的一种基于机载双天线InSAR基线测量方法与现有技术相比的优点在于：本公开克服了由于柔性基线动态变化过程难以精确建模而导致直接利用主POS至子IMU传递对准得到的位置计算基线产生较大误差的难题，不需要建立柔性基线随时间变化的动态模型，且不直接利用主POS和子IMU的位置信息计算基线，而是利用分布式POS的主POS和子IMU高精度姿态信息，通过计算主POS和子IMU之间相对角运动获得辅天线成像中心处的挠曲变形位置矢量，进而实现柔性基线的精确测量，提高了机载双天线InSAR基线测量精度，提高了InSAR干涉成像精度。
基于同一发明构思，还提供了一种基于机载双天线InSAR基线测量装置。由于此装置解决问题的原理与前述一种基于机载双天线InSAR基线测量方法相似，因此，该装置的实施可以按照前述方法的具体步骤时限，重复之处不再赘述。
如图7所示，为一个实施例中的一种基于机载双天线InSAR基线测量装置的结构示意图。该基于机载双天线InSAR基线测量装置10包括：标校模块200、获取模块400、第一计算模块600和第二计算模块800。
其中，标校模块200对初始信息进行标校；获取模块400用于基于标校的初始信息计算分布式POS相对角运动，获取辅天线相位中心处的绕曲变形角；第一计算模块600用于通过绕曲变形角计算辅天线相位中心处的绕曲变形位置矢量；第二计算模块800用于通过计算得到的辅天线相位中心处绕曲变形位置矢量，计算形变后的柔性基线。
本发明提供的一种基于机载双天线InSAR基线测量装置，通过标校模块200对初始信息进行标校；再通过获取模块400基于标校的初始信息计算分布式POS相对角运动，获取辅天线相位中心处的绕曲变形角；继而通过第一计算模块600通过绕曲变形角计算辅天线相位中心处的绕曲变形位置矢量；最终通过第二计算模块800通过计算得到的辅天线相位中心处绕曲变形位置矢量，计算形变后的柔性基线。该装置克服了由于柔性基线动态变化过程难以精确建模而导致直接利用主POS至子IMU传递对准得到的位置计算基线产生较大误差的不足。本公开具有计算精度高的特点，提高了机载双天线InSAR基线测量精度，提高了InSAR成像的精度。
本发明实施例还提供了一种计算机可读存储介质。该计算机可读存储介质上存储有计算机程序，该程序被图1或图2中处理器执行。
本发明实施例还提供了一种包含指令的计算机程序产品。当该计算机程序产品在计算机上运行时，使得计算机执行上述图1或2的方法。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)或随机存储记忆体(Random AccessMemory，RAM)等。
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。
Claims (9)
1.一种基于机载双天线InSAR基线测量方法，所述方法包括：
对初始信息进行标校；
基于标校的初始信息计算分布式POS相对角运动，获取辅天线相位中心处的绕曲变形角；
通过绕曲变形角计算辅天线相位中心处的绕曲变形位置矢量；
通过计算得到的辅天线相位中心处绕曲变形位置矢量，计算形变后的柔性基线。
2.根据权利要求1所述的方法，其中，所述对初始信息进行标校包括：利用激光全站仪测量点坐标的方式获得测量机载双天线InSAR基线所需要的初始标校信息；
其中，包括主POS坐标系与悬臂梁坐标系oBxByBzB之间坐标转换矩阵主POS坐标系与柔性基线形变前子IMU坐标系之间坐标转换矩阵悬臂梁坐标系oBxByBzB的oBxB轴方向的单位矢量eBx在主POS坐标系下的投影eBx＝[ex,ey,ez]T，ex、ey和ez分别表示eBx在轴、轴和轴的投影，形变前柔性基线长度r，形变前柔性基线r在主POS坐标系下的投影rx、ry和rz分别表示r在轴、轴和轴的投影。
3.根据权利要求1所述的方法，其中，所述基于标校的初始信息计算分布式POS相对角运动，获取辅天线相位中心处的绕曲变形角包括：利用标校的所述初始信息以及所述分布式POS的主POS和子IMU输出的高精度姿态矩阵，计算主POS和子IMU之间的相对角运动，并计算辅天线相位中心处的挠曲变形角σ。
4.根据权利要求3所述的方法，其中，所述利用标校的所述初始信息以及所述分布式POS的主POS和子IMU输出的高精度姿态矩阵，计算主POS和子IMU之间的相对角运动，并计算辅天线相位中心处的挠曲变形角包括：
利用t时刻主POS和子IMU输出的姿态矩阵以及标校获得的主POS坐标系与柔性基线形变前子IMU坐标系之间坐标转换矩阵计算形变前子IMU坐标系与形变后子IMU坐标系之间的坐标转换矩阵
其中，公式为：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>C</mi>
<msubsup>
<mi>b</mi>
<mi>s</mi>
<mo>&prime;</mo>
</msubsup>
<msub>
<mi>b</mi>
<mi>s</mi>
</msub>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>C</mi>
<msub>
<mi>b</mi>
<mi>m</mi>
</msub>
<msub>
<mi>b</mi>
<mi>s</mi>
</msub>
</msubsup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msubsup>
<mi>C</mi>
<msubsup>
<mi>b</mi>
<mi>s</mi>
<mo>&prime;</mo>
</msubsup>
<mi>n</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>C</mi>
<msub>
<mi>b</mi>
<mi>m</mi>
</msub>
<mi>n</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>13</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>23</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mn>31</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>32</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>33</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中，为主POS在t时刻输出的姿态矩阵，为子IMU在t时刻输出的姿态矩阵，Tij为矩阵的第i行、第j列元素(i,j＝1,2,3)；
利用计算t时刻辅天线相位中心处的挠曲变形角σ；
其中，公式为：
<mrow>
<mi>&sigma;</mi>
<mo>=</mo>
<mn>2</mn>
<mi>arccos</mi>
<mrow>
<mo>(</mo>
<msqrt>
<mfrac>
<mrow>
<msubsup>
<mi>e</mi>
<mi>x</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>T</mi>
<mn>11</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>e</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>T</mi>
<mn>22</mn>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>e</mi>
<mi>x</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>e</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>T</mi>
<mn>33</mn>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msubsup>
<mi>e</mi>
<mi>x</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>e</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</msqrt>
<mo>)</mo>
</mrow>
</mrow>
其中，σ的取值范围满足且当ex(T32-T23)<0时，σ<0，当ex(T32-T23)>0时，σ>0。
5.根据权利要求1所述的方法，其中，所述通过绕曲变形角计算辅天线相位中心处的绕曲变形位置矢量包括：利用标校的所述初始信息以及获得的所述辅天线相位中心处的挠曲变形角σ，通过将柔性基线形变问题等效为悬臂梁自由端大挠度弯曲问题，根据工程力学理论建立辅天线相位中心处挠曲变形角σ与挠曲变形位置矢量Δr之间的解析关系，计算得到辅天线相位中心处的挠曲变形位置矢量Δr。
6.根据权利要求5所述的方法，其中，所述通过将柔性基线形变问题等效为悬臂梁自由端大挠度弯曲问题，根据工程力学理论建立辅天线相位中心处挠曲变形角σ与挠曲变形位置矢量Δr之间的解析关系，计算得到辅天线相位中心处的挠曲变形位置矢量包括：
建立t时刻辅天线相位中心处的挠曲变形角σ与挠曲变形位置矢量Δr之间解析关系式；
其中，公式为：
<mrow>
<msub>
<mi>&Delta;r</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<msub>
<mi>r</mi>
<msub>
<mi>x</mi>
<mi>B</mi>
</msub>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;r</mi>
<msub>
<mi>y</mi>
<mi>B</mi>
</msub>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;r</mi>
<msub>
<mi>z</mi>
<mi>B</mi>
</msub>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mi>r</mi>
<mo>&CenterDot;</mo>
<msqrt>
<mfrac>
<mrow>
<mn>2</mn>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&sigma;</mi>
</mrow>
<mrow>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>&mu;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>&phi;</mi>
<mo>,</mo>
<mi>&mu;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mn>2</mn>
<mo>&lsqb;</mo>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>&mu;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>&phi;</mi>
<mo>,</mo>
<mi>&mu;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>&CenterDot;</mo>
<msqrt>
<mfrac>
<mrow>
<mi>E</mi>
<mi>I</mi>
</mrow>
<mi>P</mi>
</mfrac>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中，ΔrB为挠曲变形位置矢量Δr在悬臂梁坐标系oBxByBzB下的矢量表示， 和分别为挠曲形变Δr在悬臂梁坐标系oBxB轴、oByB轴和oBzB轴的投影，EI为柔性基线的弯曲刚度，根据柔性基线结构和材料确定，为常数，P为辅天线相位中心处的集中载荷，子IMU和辅天线的重量和，E(μ)为第二类完全椭圆积分，E(φ,μ)为第二类椭圆积分，F(μ)为第一类完全椭圆积分，F(φ,μ)为第一类完全椭圆积分；
再将所述辅天线相位中心处的挠曲变形角σ与所述挠曲变形位置矢量Δr之间解析关系式进行近似简化操作；
其中，公式为：
<mrow>
<msub>
<mi>&Delta;r</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<msub>
<mi>r</mi>
<msub>
<mi>x</mi>
<mi>B</mi>
</msub>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;r</mi>
<msub>
<mi>y</mi>
<mi>B</mi>
</msub>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;r</mi>
<msub>
<mi>z</mi>
<mi>B</mi>
</msub>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mi>r</mi>
<mo>&CenterDot;</mo>
<msqrt>
<mrow>
<mn>2</mn>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&sigma;</mi>
<mo>/</mo>
<mrow>
<mo>(</mo>
<mfrac>
<msup>
<mi>Pr</mi>
<mn>2</mn>
</msup>
<mrow>
<mi>E</mi>
<mi>I</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>r</mi>
<mi> </mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&sigma;</mi>
<mo>-</mo>
<mn>2</mn>
<msup>
<mi>r&sigma;</mi>
<mn>2</mn>
</msup>
<mo>/</mo>
<mn>3</mn>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msup>
<mi>Pr</mi>
<mn>2</mn>
</msup>
<mrow>
<mi>E</mi>
<mi>I</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>&alpha;</mi>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中，α为指数系数，其取值为
7.根据权利要求1所述的方法，其中，所述通过计算得到的辅天线相位中心处绕曲变形位置矢量，计算形变后的柔性基线包括：通过计算得到的所述辅天线相位中心处挠曲变形位置矢量ΔrB，计算导航坐标系下形变后的柔性基线r′n(t)；
其中，公式为：
<mrow>
<msup>
<mi>r</mi>
<mrow>
<mo>&prime;</mo>
<mi>n</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>C</mi>
<msub>
<mi>b</mi>
<mi>m</mi>
</msub>
<mi>n</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mo>&lsqb;</mo>
<msub>
<mi>r</mi>
<msub>
<mi>b</mi>
<mi>m</mi>
</msub>
</msub>
<mo>+</mo>
<msubsup>
<mi>C</mi>
<mi>B</mi>
<msub>
<mi>b</mi>
<mi>m</mi>
</msub>
</msubsup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>&Delta;r</mi>
<mi>B</mi>
</msub>
<mo>&rsqb;</mo>
<mo>.</mo>
</mrow>
8.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1-7任一项所述的方法。
9.一种基于机载双天线InSAR基线测量装置，所述装置包括：
标校模块，对初始信息进行标校；
获取模块，用于基于标校的初始信息计算分布式POS相对角运动，获取辅天线相位中心处的绕曲变形角；
第一计算模块，用于通过绕曲变形角计算辅天线相位中心处的绕曲变形位置矢量；
第二计算模块，用于通过计算得到的辅天线相位中心处绕曲变形位置矢量，计算形变后的柔性基线。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710837745.7A CN107765244B (zh) | 2017-09-18 | 2017-09-18 | 基于机载双天线InSAR基线测量方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710837745.7A CN107765244B (zh) | 2017-09-18 | 2017-09-18 | 基于机载双天线InSAR基线测量方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107765244A true CN107765244A (zh) | 2018-03-06 |
| CN107765244B CN107765244B (zh) | 2020-10-27 |
Family
ID=61265848
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710837745.7A Active CN107765244B (zh) | 2017-09-18 | 2017-09-18 | 基于机载双天线InSAR基线测量方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107765244B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108801166A (zh) * | 2018-05-29 | 2018-11-13 | 北京航空航天大学 | 基于悬臂梁理论的光纤光栅机翼形变测量建模及标定方法 |
| CN110068817A (zh) * | 2019-05-07 | 2019-07-30 | 中国科学院电子学研究所 | 一种基于激光测距和InSAR的地形测图方法、仪器和系统 |
| CN110297220A (zh) * | 2019-07-19 | 2019-10-01 | 西安科技大学 | 一种双天线InSAR系统基线矢量的测量方法 |
| CN111006578A (zh) * | 2019-12-26 | 2020-04-14 | 东南大学 | 一种基于gnss双天线的高铁桥墩变形监测的方法与装置 |
| CN111289013A (zh) * | 2020-04-03 | 2020-06-16 | 北京航空航天大学 | 一种分离复合多源扰动影响的机翼挠曲角速度估计方法 |
| CN113970305A (zh) * | 2021-06-16 | 2022-01-25 | 广西大学 | 一种通过挠度测量压杆轴向位移的方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102621565A (zh) * | 2012-04-17 | 2012-08-01 | 北京航空航天大学 | 一种机载分布式pos的传递对准方法 |
| CN106054185A (zh) * | 2016-05-23 | 2016-10-26 | 北京航空航天大学 | 一种基于分布式POS的机载双天线InSAR基线计算方法 |
| CN106289246A (zh) * | 2016-07-25 | 2017-01-04 | 北京航空航天大学 | 一种基于位置和姿态测量系统的柔性杆臂测量方法 |
| CN106679612A (zh) * | 2017-02-10 | 2017-05-17 | 北京航空航天大学 | 一种基于惯性测量匹配的非线性挠曲变形估计方法 |
-
2017
- 2017-09-18 CN CN201710837745.7A patent/CN107765244B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102621565A (zh) * | 2012-04-17 | 2012-08-01 | 北京航空航天大学 | 一种机载分布式pos的传递对准方法 |
| CN106054185A (zh) * | 2016-05-23 | 2016-10-26 | 北京航空航天大学 | 一种基于分布式POS的机载双天线InSAR基线计算方法 |
| CN106289246A (zh) * | 2016-07-25 | 2017-01-04 | 北京航空航天大学 | 一种基于位置和姿态测量系统的柔性杆臂测量方法 |
| CN106679612A (zh) * | 2017-02-10 | 2017-05-17 | 北京航空航天大学 | 一种基于惯性测量匹配的非线性挠曲变形估计方法 |
Non-Patent Citations (3)
| Title |
|---|
| 唐晓青 等: "考虑基线抖动的双天线干涉SAR原始回波仿真", 《电子与信息学报》 * |
| 房建成 等: "机载分布式POS传递对准建模与仿真", 《中国惯性技术学报》 * |
| 朱庄生 等: "分布式POS传递对准对InSAR干涉测量影响的分析", 《中国惯性技术学报》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108801166A (zh) * | 2018-05-29 | 2018-11-13 | 北京航空航天大学 | 基于悬臂梁理论的光纤光栅机翼形变测量建模及标定方法 |
| CN108801166B (zh) * | 2018-05-29 | 2020-05-26 | 北京航空航天大学 | 基于悬臂梁理论的光纤光栅机翼形变测量建模及标定方法 |
| CN110068817A (zh) * | 2019-05-07 | 2019-07-30 | 中国科学院电子学研究所 | 一种基于激光测距和InSAR的地形测图方法、仪器和系统 |
| CN110297220A (zh) * | 2019-07-19 | 2019-10-01 | 西安科技大学 | 一种双天线InSAR系统基线矢量的测量方法 |
| CN110297220B (zh) * | 2019-07-19 | 2021-04-20 | 西安科技大学 | 一种双天线InSAR系统基线矢量的测量方法 |
| CN111006578A (zh) * | 2019-12-26 | 2020-04-14 | 东南大学 | 一种基于gnss双天线的高铁桥墩变形监测的方法与装置 |
| CN111006578B (zh) * | 2019-12-26 | 2021-03-23 | 东南大学 | 一种基于gnss双天线的高铁桥墩变形监测的方法与装置 |
| CN111289013A (zh) * | 2020-04-03 | 2020-06-16 | 北京航空航天大学 | 一种分离复合多源扰动影响的机翼挠曲角速度估计方法 |
| CN111289013B (zh) * | 2020-04-03 | 2021-09-03 | 北京航空航天大学 | 一种分离复合多源扰动影响的机翼挠曲角速度估计方法 |
| CN113970305A (zh) * | 2021-06-16 | 2022-01-25 | 广西大学 | 一种通过挠度测量压杆轴向位移的方法 |
| CN113970305B (zh) * | 2021-06-16 | 2023-07-28 | 广西大学 | 一种通过挠度测量压杆轴向位移的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107765244B (zh) | 2020-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107765244B (zh) | 基于机载双天线InSAR基线测量方法和装置 | |
| CN110834333B (zh) | 一种机器人手眼标定方法及存储介质 | |
| CN108801166B (zh) | 基于悬臂梁理论的光纤光栅机翼形变测量建模及标定方法 | |
| CN107728182B (zh) | 基于相机辅助的柔性多基线测量方法和装置 | |
| CN108845335A (zh) | 一种基于图像和导航信息的无人机地面目标定位方法 | |
| CN114018274B (zh) | 车辆定位方法、装置及电子设备 | |
| CN108398130B (zh) | 挠曲形变测量网络的分布式pos传递对准建模方法和装置 | |
| CN108375383B (zh) | 多相机辅助的机载分布式pos柔性基线测量方法和装置 | |
| CN111043979A (zh) | 基于光纤光栅传感器阵列的板状结构形变检测装置和方法 | |
| CN108387227A (zh) | 机载分布式pos的多节点信息融合方法及系统 | |
| CN110297220B (zh) | 一种双天线InSAR系统基线矢量的测量方法 | |
| CN116540285B (zh) | 惯性辅助的gnss双天线定向方法、装置与电子设备 | |
| CN106767673A (zh) | 一种卫星高精度光学敏感载荷的指向测量方法 | |
| JP2017009555A (ja) | カメラのキャリブレーション装置、カメラのキャリブレーション方法およびカメラのキャリブレーションを行うためのプログラム | |
| CN108458709A (zh) | 基于视觉辅助测量的机载分布式pos数据融合方法和装置 | |
| CN104807477A (zh) | 一种基于靶标控制点的卫星ccd阵列影像几何检校方法 | |
| CN109141411A (zh) | 定位方法、定位装置、移动机器人及存储介质 | |
| CN114754798A (zh) | 一种陀螺误差特性参数在轨辨识与标定方法 | |
| Gu et al. | Airborne distributed POS flexible baseline measurement method based on MCLS | |
| JP5772265B2 (ja) | 地磁気測定装置、オフセット決定方法、及びオフセット決定プログラム | |
| CN107764268B (zh) | 一种机载分布式pos传递对准的方法和装置 | |
| CN116026293B (zh) | 激光gnss-rtk全站仪坐标转换方法及装置 | |
| Lu et al. | Analysis and application of geometric dilution of precision based on altitude-assisted INS/SAR integrated navigation | |
| CN115950450A (zh) | 二维多普勒测速仪的标定方法、装置、计算机设备和介质 | |
| CN112230194B (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 |