CN106643670B - 一种无人机航摄站点坐标求解装置及方法 - Google Patents

一种无人机航摄站点坐标求解装置及方法 Download PDF

Info

Publication number
CN106643670B
CN106643670B CN201710105351.2A CN201710105351A CN106643670B CN 106643670 B CN106643670 B CN 106643670B CN 201710105351 A CN201710105351 A CN 201710105351A CN 106643670 B CN106643670 B CN 106643670B
Authority
CN
China
Prior art keywords
gps
ground
track
representing
gps data
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
CN201710105351.2A
Other languages
English (en)
Other versions
CN106643670A (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.)
Xi'an Dadi Surveying And Mapping Co ltd
Original Assignee
Xi'an Dadi Surveying And Mapping 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 Xi'an Dadi Surveying And Mapping Co ltd filed Critical Xi'an Dadi Surveying And Mapping Co ltd
Priority to CN201710105351.2A priority Critical patent/CN106643670B/zh
Publication of CN106643670A publication Critical patent/CN106643670A/zh
Application granted granted Critical
Publication of CN106643670B publication Critical patent/CN106643670B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
    • G01C11/04—Interpretation of pictures
    • G01C11/30—Interpretation of pictures by triangulation
    • G01C11/34—Aerial triangulation
    • 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
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Multimedia (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种无人机航摄站点坐标求解装置，包括无人机搭载机构和地面基准站，无人机搭载机构上设置有机载GPS接收机和航摄数据采集单元，地面基准站上设置有地面GPS接收机和计算机，所述机载GPS接收机包括第一微控制器、飞机GPS天线、飞机GPS芯片和飞机GPS数据存储器；所述航摄数据采集单元包括航摄相机和存储卡；所述地面GPS接收机包括第二微控制器、地面GPS天线、地面GPS芯片和地面GPS数据存储器。本发明还公开了一种无人机航摄站点坐标求解方法。本发明通过GPS飞行轨迹约束自由网下飞行轨迹，采用事后计算，不需要使用信号线进行连接，能够避免航摄相机和GPS设备之间因为信号影响产生的同步误差。

Description

一种无人机航摄站点坐标求解装置及方法
技术领域
本发明属于航空摄影测量技术领域，具体涉及一种无人机航摄站点坐标求解装置及方法。
背景技术
近年来，轻小型测绘无人机以其成本低、灵活机动强、获取的影像分辨率高等特点，已广泛应用在应急测绘数据获取、海岛礁测绘、困难地区测绘等领域，是我国测绘航摄的重要手段和装备。GPS辅助空中三角形测量技术是提高无人机的作业效率的常用手段，为了实现GPS辅助空中三角形测量，就必须能够实现GPS设备与航测相机之间的同步，在现有的无人机航空摄影技术中，为了实现航摄相机与GPS设备的同步所采用的主要解决方案有以下三种：一、利用Y型线，将航摄相机和GPS设备并联到一起，但是由于GPS设备得到的信号时刻并不是航摄相机曝光的时刻，由于快门打开的时间并不是确定的，周围环境的改变相应的快门曝光时间也会有变化，造成GPS设备获取的时间精度误差大，从而影响到曝光点的坐标精度；二、通过航摄相机闪光灯的信号实现航摄相机和GPS设备的同步，这种方式在一定程度上提高了同步的精度，使得时间误差仅仅存在于航摄相机快门曝光的误差，有较高的精度；三、从航摄相机CCD处直接获得航摄相机曝光瞬间的信号，将此信号供给GPS设备，这种方式是精度最好的方式，能够直接获取航摄相机感光瞬间的信号，从而精确求得摄站点的坐标，但是这种方式需要专业的航摄相机才能够实现对于一般测绘公司来说应用成本过高。
为了解决上述问题，提出了一种航摄相机和GPS设备独立工作的无人机航摄站点坐标求解装置，通过GPS设备测量的飞行轨迹约束航摄相机测量得到的飞行相似轨迹，可算出航摄相机曝光点的坐标，不需要使用信号线进行连接，能够避免航摄相机和GPS设备之间因为信号影响产生的同步误差，同时也能降低硬件使用时的难度，使用效果好。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种无人机航摄站点坐标求解装置及方法，其结构简单、设计合理，通过GPS设备测量的GPS飞行轨迹约束航摄相机测量得到的自由网下飞行轨迹，采用事后计算的方法来获取摄站点坐标，不需要使用信号线进行连接，能够避免航摄相机和GPS设备之间因为信号影响产生的同步误差，实用性强，使用效果好，便于推广使用。
为解决上述技术问题，本发明采用的技术方案是：一种无人机航摄站点坐标求解装置，其特征在于：包括无人机搭载机构和地面基准站，所述无人机搭载机构上设置有机载GPS接收机和航摄数据采集单元，所述地面基准站上设置有地面GPS接收机和用于进行数据处理的计算机，所述机载GPS接收机包括第一微控制器和飞机GPS天线，所述飞机GPS天线与第一微控制器之间接有飞机GPS芯片，所述第一微控制器接有飞机GPS数据存储器，所述飞机GPS数据存储器与计算机通过通信接口连接；所述航摄数据采集单元包括用于采集测区影像信息的航摄相机和用于存储航摄相机信息的存储卡，所述存储卡与计算机通过通信接口连接；所述地面GPS接收机包括第二微控制器和地面GPS天线，所述地面GPS天线和第二微控制器之间接有地面GPS芯片，所述第二微控制器接有地面GPS数据存储器，所述地面GPS数据存储器与计算机通过通信接口连接。
上述的一种无人机航摄站点坐标求解装置，其特征在于：所述飞机GPS天线和飞机GPS芯片采用屏蔽线连接，所述地面GPS天线和地面GPS芯片采用屏蔽线连接。
上述的一种无人机航摄站点坐标求解装置，其特征在于：所述通信接口为USB接口。
利用上述的装置进行航摄站点坐标求解的方法，其特征在于，该方法包括以下步骤：
步骤一、GPS测量的飞行轨迹获取：
步骤101、机载GPS数据获取：
步骤1011、获取无人机与卫星之间的伪距：利用飞机GPS天线和飞机GPS芯片获取卫星发出的与无人机位置相关的GPS信号：
Figure BDA0001232960340000031
其中i表示卫星的个数且i≥3，
Figure BDA0001232960340000032
表示飞机GPS天线与第i颗卫星之间的伪距，C表示电磁波的传播速度，dδu表示飞机GPS芯片的时钟偏差；
Figure BDA0001232960340000033
表示第i颗卫星的时钟偏差；
Figure BDA0001232960340000034
表示第i颗卫星的星历引起的机载GPS数据误差；dρuion表示高空大气层中电离层效应引起的机载GPS数据偏差；dρutrop表示高空大气层中对流层时延引起的机载GPS数据偏差；dMu表示多路径效应引起的机载GPS数据偏差；vu表示机载GPS接收机的噪声值，
Figure BDA0001232960340000035
表示计算得到的无人机与第i颗卫星之间的真实距离，其中，
Figure BDA0001232960340000036
(xu，yu，zu)表示飞机GPS天线的坐标位置，
Figure BDA0001232960340000037
表示第i颗卫星的坐标位置；
步骤1012、将步骤1011中获得的GPS信号存储在飞机GPS数据存储器中；
步骤102、地面GPS数据获取：
步骤1021、获取地面基准站与卫星之间的伪距：利用地面GPS天线和地面GPS芯片获取卫星发出的与地面基准站位置相关的距离信号：
Figure BDA0001232960340000041
其中，
Figure BDA0001232960340000042
表示测量到的地面基准站与第i颗卫星之间的伪距，
Figure BDA0001232960340000043
表示计算得到的地面基准站到第i颗卫星的真实距离；dδb表示地面GPS芯片的时钟偏差；
Figure BDA0001232960340000044
表示第i颗卫星的时钟偏差；
Figure BDA0001232960340000045
表示第i颗卫星的星历引起的地面GPS数据误差；dρbion表示电离层效应引起的地面GPS数据偏差；dρbtrop表示对流层时延引起的地面GPS数据偏差；dMb表示多路径效应引起的地面GPS数据偏差；vb表示地面GPS接收机噪声值；
步骤1022、将步骤2011中获得的数据存储在地面GPS数据存储器中；
步骤103、伪距修正：
步骤1031、获取伪距修正值：计算机调用地面GPS数据存储器中的数据，进行计算：
Figure BDA0001232960340000046
其中
Figure BDA0001232960340000047
表示地面基准站测得的地面基准站与第i颗卫星之间的伪距修正值；
步骤1032、改正无人机与卫星之间的伪距：
当无人机距离地面基准站的距离小于1000km时，
Figure BDA0001232960340000048
dρuion≈dρbion，dρutrop≈dρbtrop，设Δdρ＝C(dδu-dδb)+(dMu-dMb)+(vu-vb)，因此，
Figure BDA0001232960340000051
步骤104、获得GPS飞行轨迹：求解得到：GPS飞行轨迹U(t)＝(xu,yu,zu)，
步骤二、自由网情况下飞行轨迹的获取：
步骤201、影像信息获取：利用航摄相机对相应的像点摄像得到像点照片，并存储在存储卡中，计算机利用立体设备求得像点照片的像点坐标(xj,yj)，通过工作人员在测区地面采集测区地面点影像信息(Xj,Yj,Zj)，其中，j＝1,2，3；
步骤202、计算机对影像信息进行空三处理：
Figure BDA0001232960340000052
其中，
Figure BDA0001232960340000061
其中，
Figure BDA0001232960340000062
表示拍摄像点照片时无人机机头所指方向和地面的夹角，ω表示拍摄像点照片时无人机机头所指方向在地面上的投影和正北方向的夹角、κ表示拍摄像点照片时无人机机翼的倾斜角，f表示拍摄像点照片时航摄相机的焦距；
步骤203、获得自由网下飞行轨迹：求解得到自由网下飞行轨迹：
Figure BDA0001232960340000063
设S(t)＝(XS,YS,ZS)，S(t)表示自由网情况下摄站点的坐标。
步骤三、空间直线参数提取：
步骤301、提取GPS飞行轨迹空间直线参数：在GPS飞行轨迹U(t)上任取两个轨迹点
Figure BDA0001232960340000064
和
Figure BDA0001232960340000065
Figure BDA0001232960340000066
表示th时刻的GPS飞行轨迹点，
Figure BDA0001232960340000067
表示th+1时刻的GPS飞行轨迹点，利用向量
Figure BDA0001232960340000068
描述GPS飞行轨迹：
Figure BDA0001232960340000069
Figure BDA00012329603400000610
表示单位向量，
Figure BDA00012329603400000611
表示垂直于单位向量
Figure BDA00012329603400000612
的法向量；
步骤302、提取自由网下飞行轨迹空间直线参数：在自由网下飞行轨迹S(t)上任取两个轨迹点
Figure BDA00012329603400000613
和
Figure BDA00012329603400000614
Figure BDA00012329603400000615
表示tj时刻的自由网下飞行轨迹点，
Figure BDA00012329603400000616
表示tj+1时刻的自由网下飞行轨迹点，利用向量
Figure BDA00012329603400000617
描述自由网下飞行轨迹：
Figure BDA00012329603400000618
Figure BDA00012329603400000619
表示轨迹点
Figure BDA00012329603400000620
和轨迹点
Figure BDA00012329603400000621
在空间内的单位向量，
Figure BDA00012329603400000622
表示垂直于单位向量
Figure BDA00012329603400000623
的法向量；
步骤303、求解相似变换参数：相似变换参数
Figure BDA0001232960340000071
Figure BDA0001232960340000072
其中μ表示自由网下的轨迹转换到GPS轨迹时的尺度因子，θ表示自由网下的轨迹转换到GPS轨迹时的旋转矩，q表示自由网下的轨迹转换到GPS轨迹时的平移参数，利用最小二乘法求解，得到T(μ,θ,q)；
步骤四、曲线拟合得到摄站点坐标：
步骤401、自由网下飞行轨迹空间坐标转换：对自由网情况下摄站点的坐标(XS,YS,ZS)进行坐标转换，
Figure BDA0001232960340000073
其中S(t)'＝(X′S,Y′S,Z′S)表示相似转换后的自由网下飞行轨迹空间坐标；
步骤402、建立曲线拟合模型：建立曲线拟合模型
Figure BDA0001232960340000074
其中U(t)'表示曲线拟合后得到的GPS飞行轨迹点坐标，Nh,k(t)为K阶B样条基函数，
Figure BDA0001232960340000075
th为GPS飞行轨迹点的对应时间，Nh,1(t)表示1阶B样条基函数在第h段区间上的值；
步骤403：求摄站点坐标的准确值：根据最小二乘优化公式求摄站点坐标P的准确值，得到：
Figure BDA0001232960340000076
本发明与现有技术相比具有以下优点：
1、本发明的结构简单、设计合理，实现及使用操作方便。
2、本发明通过设置飞机GPS数据存储器将测量过程中得到的机载GPS数据进行保存，通过设置地面GPS数据存储器将测量过程中得到的地面GPS数据进行保存，测量结束后，对测量得到的机载GPS数据和地面GPS数据进行事后计算，不需要使用信号线进行连接，避免了实时差分中存在的由于数据传输以及整周跳变对测量精度的影响。
3、本发明对航摄相机和GPS设备之间的曲线拟合采用事后计算，避免因为信号影响产生的同步误差，提高了摄站点坐标精度，同时也能降低硬件使用时的难度，使用效果好。
4、本发明以地面GPS数据作为差分基准站数据，通过对机载GPS数据和地面GPS数据进行差分计算，得到GPS飞行轨迹，对航摄相机采集的影像数据进行空三处理，得到自由网下飞行轨迹，通过对两条飞机飞行轨迹进行曲线拟合，通过GPS飞行轨迹约束自由网下飞行轨迹，得到摄站点坐标的准确值，提高了摄站点坐标精度。
综上所述，本发明结构简单、设计合理，通过GPS设备测量的GPS飞行轨迹约束航摄相机测量得到的自由网下飞行轨迹，采用事后计算，不需要使用信号线进行连接，并且能够避免航摄相机和GPS设备之间因为信号影响产生的同步误差，实用性强，使用效果好，便于推广使用。
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明无人机航摄站点坐标求解装置的电路框图。
图2为本发明无人机航摄站点坐标求解方法的流程框图。
附图标记说明:
1—飞机GPS数据存储器； 2—飞机GPS天线；
3—飞机GPS芯片； 4—第一微控制器； 5—存储卡；
6—航摄相机； 7—地面GPS天线； 8—地面GPS芯片；
9—第二微控制器； 10—地面GPS数据存储器；
11—计算机。
具体实施方式
如图1所示的一种无人机航摄站点坐标求解装置，包括无人机搭载机构和地面基准站，所述无人机搭载机构上设置有机载GPS接收机和航摄数据采集单元，所述地面基准站上设置有地面GPS接收机和用于进行数据处理的计算机11，所述机载GPS接收机包括第一微控制器4和飞机GPS天线2，所述飞机GPS天线2与第一微控制器4之间接有飞机GPS芯片3，所述第一微控制器4接有飞机GPS数据存储器1，所述飞机GPS数据存储器1与计算机11通过通信接口连接；所述航摄数据采集单元包括用于采集测区影像信息的航摄相机6和用于存储航摄相机6信息的存储卡5，所述存储卡5与计算机11通过通信接口连接；所述地面GPS接收机包括第二微控制器9和地面GPS天线7，所述地面GPS天线7和第二微控制器9之间接有地面GPS芯片8，所述第二微控制器9接有地面GPS数据存储器10，所述地面GPS数据存储器10与计算机11通过通信接口连接。
本实施例中，所述飞机GPS天线2和飞机GPS芯片3采用屏蔽线连接，所述地面GPS天线7和地面GPS芯片8采用屏蔽线连接。
采用屏蔽线连接可以有效屏蔽外界对导线中传输信号的干扰。
本实施例中，所述通信接口为USB接口。
如图1和如图2所示，本发明一种无人机航摄站点坐标求解方法，包括以下步骤：
步骤一、GPS测量的飞行轨迹获取：
步骤101、机载GPS数据获取：
步骤1011、获取无人机与卫星之间的伪距：利用飞机GPS天线2和飞机GPS芯片3获取卫星发出的与无人机位置相关的GPS信号：
Figure BDA0001232960340000091
其中i表示卫星的个数且i≥3，
Figure BDA0001232960340000092
表示飞机GPS天线2与第i颗卫星之间的伪距，C表示电磁波的传播速度，dδu表示飞机GPS芯片3的时钟偏差；
Figure BDA0001232960340000093
表示第i颗卫星的时钟偏差；
Figure BDA0001232960340000094
表示第i颗卫星的星历引起的机载GPS数据误差；dρuion表示高空大气层中电离层效应引起的机载GPS数据偏差；dρutrop表示高空大气层中对流层时延引起的机载GPS数据偏差；dMu表示多路径效应引起的机载GPS数据偏差；vu表示机载GPS接收机的噪声值，
Figure BDA0001232960340000101
表示计算得到的无人机与第i颗卫星之间的真实距离，其中，
Figure BDA0001232960340000102
(xu，yu，zu)表示飞机GPS天线2的坐标位置，
Figure BDA0001232960340000103
表示第i颗卫星的坐标位置；
需要说明的是，步骤1011中，
Figure BDA0001232960340000104
dρuion、dρutrop、dMu和
Figure BDA0001232960340000105
的单位均为米，C的单位为米/秒，dδu和
Figure BDA0001232960340000106
单位均为ppm；vu的单位为分贝。
步骤1012、将步骤1011中获得的GPS信号存储在飞机GPS数据存储器1中；
步骤102、地面GPS数据获取：
步骤1021、获取地面基准站与卫星之间的伪距：利用地面GPS天线7和地面GPS芯片8获取卫星发出的与地面基准站位置相关的距离信号：
Figure BDA0001232960340000107
其中，
Figure BDA0001232960340000108
表示测量到的地面基准站与第i颗卫星之间的伪距，
Figure BDA0001232960340000109
表示计算得到的地面基准站到第i颗卫星的真实距离；dδb表示地面GPS芯片8的时钟偏差；
Figure BDA00012329603400001010
表示第i颗卫星的时钟偏差；
Figure BDA00012329603400001011
表示第i颗卫星的星历引起的地面GPS数据误差；dρbion表示电离层效应引起的地面GPS数据偏差；dρbtrop表示对流层时延引起的地面GPS数据偏差；dMb表示多路径效应引起的地面GPS数据偏差；vb表示地面GPS接收机噪声值；
需要说明的是，步骤1021中，
Figure BDA0001232960340000111
dρbion、dρbtrop和dMb的单位均为米，dδb的单位为ppm；vb的单位为分贝。
步骤1022、将步骤2011中获得的数据存储在地面GPS数据存储器10中；
步骤103、伪距修正：
步骤1031、获取伪距修正值：计算机11调用地面GPS数据存储器10中的数据，进行计算：
Figure BDA0001232960340000112
其中
Figure BDA0001232960340000113
表示地面基准站测得的地面基准站与第i颗卫星之间的伪距修正值；
步骤1032、改正无人机与卫星之间的伪距：
当无人机距离地面基准站的距离小于1000km时，
Figure BDA0001232960340000114
dρuion≈dρbion，dρutrop≈dρbtrop，设Δdρ＝C(dδu-dδb)+(dMu-dMb)+(vu-vb)，因此，
Figure BDA0001232960340000115
步骤104、获得GPS飞行轨迹：求解得到GPS飞行轨迹U(t)＝(xu,yu,zu)，
步骤二、自由网情况下飞行轨迹的获取：
步骤201、影像信息获取：利用航摄相机6对相应的像点摄像得到像点照片，并存储在存储卡5中，计算机11利用立体设备求得像点照片的像点坐标(xj,yj)，通过工作人员在测区地面采集测区地面点影像信息(Xj,Yj,Zj)，其中，j＝1,2，3；
步骤202、计算机11对影像信息进行空三处理：
Figure BDA0001232960340000121
其中，
Figure BDA0001232960340000122
Figure BDA0001232960340000123
表示拍摄像点照片时无人机机头所指方向和地面的夹角，ω表示拍摄像点照片时无人机机头所指方向在地面上的投影和正北方向的夹角、κ表示拍摄像点照片时无人机机翼的倾斜角，
Figure BDA0001232960340000124
ω和κ三个角度表示像点照片的空间姿态，f表示拍摄像点照片时航摄相机6的焦距。
步骤203、获得自由网下飞行轨迹：求解得到自由网下飞行轨迹
Figure BDA0001232960340000125
设S(t)＝(XS,YS,ZS)，S(t)表示自由网情况下摄站点的坐标。
步骤三、空间直线参数提取：
步骤301、提取GPS飞行轨迹空间直线参数：在GPS飞行轨迹U(t)上任取两个轨迹点
Figure BDA0001232960340000131
和
Figure BDA0001232960340000132
Figure BDA0001232960340000133
表示th时刻的GPS飞行轨迹点，
Figure BDA0001232960340000134
表示th+1时刻的GPS飞行轨迹点，利用向量
Figure BDA0001232960340000135
描述GPS飞行轨迹：
Figure BDA0001232960340000136
Figure BDA0001232960340000137
表示单位向量，
Figure BDA0001232960340000138
表示垂直于单位向量
Figure BDA0001232960340000139
的法向量；
步骤302、提取自由网下飞行轨迹空间直线参数：在自由网下飞行轨迹S(t)上任取两个轨迹点
Figure BDA00012329603400001310
和
Figure BDA00012329603400001311
Figure BDA00012329603400001312
表示tj时刻的自由网下飞行轨迹点，
Figure BDA00012329603400001313
表示tj+1时刻的自由网下飞行轨迹点，利用向量
Figure BDA00012329603400001314
描述自由网下飞行轨迹：
Figure BDA00012329603400001315
Figure BDA00012329603400001316
表示轨迹点
Figure BDA00012329603400001317
和轨迹点
Figure BDA00012329603400001318
在空间内的单位向量，
Figure BDA00012329603400001319
表示垂直于单位向量
Figure BDA00012329603400001320
的法向量；
步骤303、求解相似变换参数：相似变换参数
Figure BDA00012329603400001321
Figure BDA00012329603400001322
其中μ表示自由网下的轨迹转换到GPS轨迹时的尺度因子，θ表示自由网下的轨迹转换到GPS轨迹时的旋转矩，q表示自由网下的轨迹转换到GPS轨迹时的平移参数，利用最小二乘法求解，得到T(μ,θ,q)；
步骤四、曲线拟合得到摄站点坐标：
步骤401、自由网下飞行轨迹空间坐标转换：对自由网情况下摄站点的坐标(XS,YS,ZS)进行坐标转换，
Figure BDA00012329603400001323
其中S(t)'＝(X′S,Y′S,Z′S)表示相似转换后的自由网下飞行轨迹空间坐标；
步骤402、建立曲线拟合模型：建立曲线拟合模型
Figure BDA00012329603400001324
其中U(t)'表示曲线拟合后得到的GPS飞行轨迹点坐标，Nh,k(t)为K阶B样条基函数，
Figure BDA0001232960340000141
th为GPS飞行轨迹点的对应时间，Nh,1(t)表示1阶B样条基函数在第h段区间上的值；
步骤403：求摄站点坐标的准确值：根据最小二乘优化公式求摄站点坐标P的准确值，得到：
Figure BDA0001232960340000142
具体实施时，将飞机GPS数据存储器1、飞机GPS天线2、飞机GPS芯片3、第一微控制器4、航摄相机6和存储卡5都搭设在无人机上。通过航摄相机6依据相关的航测规范采集测区影像信息，存储卡5用于存储测区影像信息。利用飞机GPS天线2和飞机GPS芯片3接收卫星信号，并通过第一微控制器4将卫星信号存储在飞机GPS数据存储器1上，利用地面GPS天线7和地面GPS芯片8在地面上接收卫星信号，并通过第二微控制器9将卫星信号存储在地面GPS数据存储器10上，飞机GPS芯片3和地面GPS芯片8采集GPS数据的数据更新率均为10Hz，使得机载GPS数据和地面GPS数据频率对应。GPS数据采集和测区影像信息采集之间不需要进行数据通讯，完全独立工作进行数据采集，能够避免测区影像信息采集和GPS数据采集之间因为信号影响产生的同步误差。
在GPS数据采集和测区影像信息采集结束以后，计算机11调用地面GPS数据存储器10、飞机GPS数据存储器1和存储卡5中的数据，对测量得到的机载GPS数据和地面GPS数据进行事后计算，不需要使用信号线进行连接，避免了实时差分中存在的由于数据传输以及整周跳变对测量精度的影响。将地面站作为差分基准站，将测得的机载GPS数据与已知的地面GPS数据进行比较，确定误差，得出准确的改正值，得到GPS飞行轨迹U(t)＝(xu,yu,zu)，即为事后差分GPS的工作过程。
计算机11对影像信息进行空三处理，得到自由网情况下摄站点的坐标S(t)＝(XS,YS,ZS)，S(t)即为自由网情况下飞行轨迹。由于GPS飞行轨迹U(t)和自由网情况下飞行轨迹S(t)在空间内是两条相似的三维曲线，计算得到两条三维曲线的相似变换参数T(μ,θ,q)，对GPS飞行轨迹U(t)和自由网情况下飞行轨迹S(t)进行空间直线参数提取，对所得的直线进行相似变换，求取尺度因子μ、旋转矩阵初值θ和平移参数q，其中μ为尺度因子，θ为3阶正交矩阵旋转矩阵，q为3×1矩阵。根据相似变换参数T(μ,θ,q)，求出自由网情况下飞行轨迹S(t)转换到GPS飞行轨迹U(t)上的转换坐标S(t)'＝(X′S,Y′S,Z′S)。采用B样条曲线Nh,k(t)进行曲线拟合，时间th在B样条曲线Nh,k(t)上进行滑动，当曲线拟合后得到的GPS飞行轨迹点坐标U(t)'与相似转换后的自由网下飞行轨迹空间坐标S(t)'距离最近的时候，得到摄站点坐标P，将摄站点坐标P用非线性最小二乘法优化，最终得到相应的摄站点坐标P的准确值。
以上所述，仅是本发明的实施例，并非对本发明作任何限制，凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化，均仍属于本发明技术方案的保护范围内。

Claims (4)

1.一种无人机航摄站点坐标求解方法，其特征在于：采用的无人机航摄站点坐标求解装置包括无人机搭载机构和地面基准站，所述无人机搭载机构上设置有机载GPS接收机和航摄数据采集单元，所述地面基准站上设置有地面GPS接收机和用于进行数据处理的计算机(11)，所述机载GPS接收机包括第一微控制器(4)和飞机GPS天线(2)，所述飞机GPS天线(2)与第一微控制器(4)之间接有飞机GPS芯片(3)，所述第一微控制器(4)接有飞机GPS数据存储器(1)，所述飞机GPS数据存储器(1)与计算机(11)通过通信接口连接；所述航摄数据采集单元包括用于采集测区影像信息的航摄相机(6)和用于存储航摄相机(6)信息的存储卡(5)，所述存储卡(5)与计算机(11)通过通信接口连接；所述地面GPS接收机包括第二微控制器(9)和地面GPS天线(7)，所述地面GPS天线(7)和第二微控制器(9)之间接有地面GPS芯片(8)，所述第二微控制器(9)接有地面GPS数据存储器(10)，所述地面GPS数据存储器(10)与计算机(11)通过通信接口连接；
该方法包括以下步骤：
将飞机GPS数据存储器(1)、飞机GPS天线(2)、飞机GPS芯片(3)、第一微控制器(4)、航摄相机(6)和存储卡(5)都搭设在无人机上；通过航摄相机(6)依据相关的航测规范采集测区影像信息，存储卡(5)用于存储测区影像信息；利用飞机GPS天线(2)和飞机GPS芯片(3)接收卫星信号，并通过第一微控制器(4)将卫星信号存储在飞机GPS数据存储器(1)上，利用地面GPS天线(7)和地面GPS芯片(8)在地面上接收卫星信号，并通过第二微控制器(9)将卫星信号存储在地面GPS数据存储器(10)上，飞机GPS芯片(3)和地面GPS芯片(8)采集GPS数据的数据更新率均为10Hz，使得机载GPS数据和地面GPS数据频率对应；GPS数据采集和测区影像信息采集之间不需要进行数据通讯，完全独立工作进行数据采集，能够避免测区影像信息采集和GPS数据采集之间因为信号影响产生的同步误差；
在GPS数据采集和测区影像信息采集结束以后，计算机(11)调用地面GPS数据存储器(10)、飞机GPS数据存储器(1)和存储卡(5)中的数据，对测量得到的机载GPS数据和地面GPS数据进行事后计算，不需要使用信号线进行连接，避免了实时差分中存在的由于数据传输以及整周跳变对测量精度的影响；将地面站作为差分基准站，将测得的机载GPS数据与已知的地面GPS数据进行比较，确定误差，得出准确的改正值，得到GPS飞行轨迹U(t)＝(xu,yu,zu)，即为事后差分GPS的工作过程；
计算机(11)对影像信息进行空三处理，得到自由网情况下摄站点的坐标S(t)＝(XS,YS,ZS)，S(t)即为自由网情况下飞行轨迹；由于GPS飞行轨迹U(t)和自由网情况下飞行轨迹S(t)在空间内是两条相似的三维曲线，计算得到两条三维曲线的相似变换参数T(μ,θ,q)，对GPS飞行轨迹U(t)和自由网情况下飞行轨迹S(t)进行空间直线参数提取，对所得的直线进行相似变换，求取尺度因子μ、旋转矩阵初值θ和平移参数q，其中μ为尺度因子，θ为3阶正交矩阵旋转矩阵，q为3×1矩阵；根据相似变换参数T(μ,θ,q)，求出自由网情况下飞行轨迹S(t)转换到GPS飞行轨迹U(t)上的转换坐标S(t)'＝(X′S,Y′S,Z′S)；采用B样条曲线Nh,k(t)进行曲线拟合，时间th在B样条曲线Nh,k(t)上进行滑动，当曲线拟合后得到的GPS飞行轨迹点坐标U(t)'与相似转换后的自由网下飞行轨迹空间坐标S(t)'距离最近的时候，得到摄站点坐标P，将摄站点坐标P用非线性最小二乘法优化，最终得到相应的摄站点坐标P的准确值。
2.按照权利要求1所述的一种无人机航摄站点坐标求解方法，其特征在于：所述飞机GPS天线(2)和飞机GPS芯片(3)采用屏蔽线连接，所述地面GPS天线(7)和地面GPS芯片(8)采用屏蔽线连接。
3.按照权利要求1所述的一种无人机航摄站点坐标求解方法，其特征在于：所述通信接口为USB接口。
4.一种利用如权利要求1所述的装置进行航摄站点坐标求解的方法，其特征在于，该方法包括以下步骤：
步骤一、GPS测量的飞行轨迹获取：
步骤101、机载GPS数据获取：
步骤1011、获取无人机与卫星之间的伪距：利用飞机GPS天线(2)和飞机GPS芯片(3)获取卫星发出的与无人机位置相关的GPS信号：
Figure FDA0003882790350000031
其中i表示卫星的个数且i≥3，
Figure FDA0003882790350000032
表示飞机GPS天线(2)与第i颗卫星之间的伪距，C表示电磁波的传播速度，dδu表示飞机GPS芯片(3)的时钟偏差；
Figure FDA0003882790350000033
表示第i颗卫星的时钟偏差；
Figure FDA0003882790350000034
表示第i颗卫星的星历引起的机载GPS数据误差；dρuion表示高空大气层中电离层效应引起的机载GPS数据偏差；dρutrop表示高空大气层中对流层时延引起的机载GPS数据偏差；dMu表示多路径效应引起的机载GPS数据偏差；vu表示机载GPS接收机的噪声值，
Figure FDA0003882790350000035
表示计算得到的无人机与第i颗卫星之间的真实距离，其中，
Figure FDA0003882790350000036
(xu，yu，zu)表示飞机GPS天线(2)的坐标位置，
Figure FDA0003882790350000037
表示第i颗卫星的坐标位置；
步骤1012、将步骤1011中获得的GPS信号存储在飞机GPS数据存储器(1)中；
步骤102、地面GPS数据获取：
步骤1021、获取地面基准站与卫星之间的伪距：利用地面GPS天线(7)和地面GPS芯片(8)获取卫星发出的与地面基准站位置相关的距离信号：
Figure FDA0003882790350000041
其中，
Figure FDA0003882790350000042
表示测量到的地面基准站与第i颗卫星之间的伪距，
Figure FDA0003882790350000043
表示计算得到的地面基准站到第i颗卫星的真实距离；dδb表示地面GPS芯片(8)的时钟偏差；
Figure FDA0003882790350000044
表示第i颗卫星的时钟偏差；
Figure FDA0003882790350000045
表示第i颗卫星的星历引起的地面GPS数据误差；dρbion表示电离层效应引起的地面GPS数据偏差；dρbtrop表示对流层时延引起的地面GPS数据偏差；dMb表示多路径效应引起的地面GPS数据偏差；vb表示地面GPS接收机噪声值；
步骤1022、将步骤2011中获得的数据存储在地面GPS数据存储器(10)中；
步骤103、伪距修正：
步骤1031、获取伪距修正值：计算机(11)调用地面GPS数据存储器(10)中的数据，进行计算：
Figure FDA0003882790350000046
其中
Figure FDA0003882790350000047
表示地面基准站测得的地面基准站与第i颗卫星之间的伪距修正值；
步骤1032、改正无人机与卫星之间的伪距：
当无人机距离地面基准站的距离小于1000km时，
Figure FDA0003882790350000048
dρuion≈dρbion，dρutrop≈dρbtrop，设Δdρ＝C(dδu-dδb)+(dMu-dMb)+(vu-vb)，因此，
Figure FDA0003882790350000051
步骤104、获得GPS飞行轨迹：求解得到GPS飞行轨迹U(t)＝(xu,yu,zu)，
步骤二、自由网情况下飞行轨迹的获取：
步骤201、影像信息获取：利用航摄相机(6)对相应的像点摄像得到像点照片，并存储在存储卡(5)中，计算机(11)利用立体设备求得像点照片的像点坐标(xj,yj)，通过工作人员在测区地面采集测区地面点影像信息(Xj,Yj,Zj)，其中，j＝1，2，3；
步骤202、计算机(11)对影像信息进行空三处理：
Figure FDA0003882790350000052
其中，
Figure FDA0003882790350000061
Figure FDA0003882790350000062
表示拍摄像点照片时无人机机头所指方向和地面的夹角，ω表示拍摄像点照片时无人机机头所指方向在地面上的投影和正北方向的夹角、κ表示拍摄像点照片时无人机机翼的倾斜角，f表示拍摄像点照片时航摄相机(6)的焦距；
步骤203、获得自由网下飞行轨迹：求解得到自由网下飞行轨迹
Figure FDA0003882790350000063
设S(t)＝(XS,YS,ZS)，S(t)表示自由网情况下摄站点的坐标；
步骤三、空间直线参数提取：
步骤301、提取GPS飞行轨迹空间直线参数：在GPS飞行轨迹U(t)上任取两个轨迹点
Figure FDA0003882790350000064
和
Figure FDA0003882790350000065
表示th时刻的GPS飞行轨迹点，
Figure FDA0003882790350000066
表示th+1时刻的GPS飞行轨迹点，利用向量
Figure FDA0003882790350000067
描述GPS飞行轨迹：
Figure FDA0003882790350000068
Figure FDA0003882790350000069
表示单位向量，
Figure FDA00038827903500000610
表示垂直于单位向量
Figure FDA00038827903500000611
的法向量；
步骤302、提取自由网下飞行轨迹空间直线参数：在自由网下飞行轨迹S(t)上任取两个轨迹点
Figure FDA00038827903500000612
和
Figure FDA00038827903500000613
表示tj时刻的自由网下飞行轨迹点，
Figure FDA00038827903500000614
表示tj+1时刻的自由网下飞行轨迹点，利用向量
Figure FDA00038827903500000615
描述自由网下飞行轨迹：
Figure FDA00038827903500000616
Figure FDA00038827903500000617
表示轨迹点
Figure FDA00038827903500000618
和轨迹点
Figure FDA00038827903500000619
在空间内的单位向量，
Figure FDA00038827903500000620
表示垂直于单位向量
Figure FDA00038827903500000621
的法向量；
步骤303、求解相似变换参数：相似变换参数
Figure FDA0003882790350000071
Figure FDA0003882790350000072
其中μ表示自由网下的轨迹转换到GPS轨迹时的尺度因子，θ表示自由网下的轨迹转换到GPS轨迹时的旋转矩，q表示自由网下的轨迹转换到GPS轨迹时的平移参数，利用最小二乘法求解，得到T(μ,θ,q)；
步骤四、曲线拟合得到摄站点坐标：
步骤401、自由网下飞行轨迹空间坐标转换：对自由网情况下摄站点的坐标(XS,YS,ZS)进行坐标转换，
Figure FDA0003882790350000073
其中S(t)'＝(X′S,Y′S,Z′S)表示相似转换后的自由网下飞行轨迹空间坐标；
步骤402、建立曲线拟合模型：建立曲线拟合模型
Figure FDA0003882790350000074
其中U(t)'表示曲线拟合后得到的GPS飞行轨迹点坐标，Nh,k(t)为K阶B样条基函数，
Figure FDA0003882790350000075
th为GPS飞行轨迹点的对应时间，Nh,1(t)表示1阶B样条基函数在第h段区间上的值；
步骤403：求摄站点坐标的准确值：根据最小二乘优化公式求摄站点坐标P的准确值，得到：
Figure FDA0003882790350000076
CN201710105351.2A 2017-02-26 2017-02-26 一种无人机航摄站点坐标求解装置及方法 Active CN106643670B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710105351.2A CN106643670B (zh) 2017-02-26 2017-02-26 一种无人机航摄站点坐标求解装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710105351.2A CN106643670B (zh) 2017-02-26 2017-02-26 一种无人机航摄站点坐标求解装置及方法

Publications (2)

Publication Number Publication Date
CN106643670A CN106643670A (zh) 2017-05-10
CN106643670B true CN106643670B (zh) 2022-12-09

Family

ID=58847892

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710105351.2A Active CN106643670B (zh) 2017-02-26 2017-02-26 一种无人机航摄站点坐标求解装置及方法

Country Status (1)

Country Link
CN (1) CN106643670B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107818697B (zh) * 2017-10-30 2020-11-10 中国科学院遥感与数字地球研究所 基于地形高程的非水平航线设计方法、终端及存储介质
CN107894620A (zh) * 2017-11-02 2018-04-10 中国科学院电子学研究所 一种航空编码混场源电磁探测系统
CN109612447B (zh) * 2018-12-29 2020-12-15 湖南璇玑信息科技有限公司 遥感影像地图数据的增强定位变换模型的构建方法、增强定位方法及增强定位服务器
WO2021007855A1 (zh) * 2019-07-18 2021-01-21 深圳市大疆创新科技有限公司 基站、像控点定位方法、电子设备和计算机可读介质
CN111536871B (zh) * 2020-05-07 2022-05-31 武汉大势智慧科技有限公司 一种多时相摄影测量数据体积变化量的精确计算方法
CN115442743B (zh) * 2022-11-07 2023-03-24 上海特金信息科技有限公司 基于rid信号数据融合的定位方法、装置、设备

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101201248A (zh) * 2007-12-07 2008-06-18 中国科学院武汉岩土力学研究所 基于无人机的航空近景摄影位移测量系统及其方法
CN102937443A (zh) * 2012-01-13 2013-02-20 唐粮 一种基于无人机的目标快速定位系统及方法
CN105182384A (zh) * 2015-08-24 2015-12-23 桂林电子科技大学 一种双模实时伪距差分定位系统和伪距改正数据生成方法
CN106093987A (zh) * 2016-08-10 2016-11-09 深圳供电局有限公司 一种应用于无人机的低成本差分gps系统及其实现方法
CN206479174U (zh) * 2017-02-26 2017-09-08 西安大地测绘股份有限公司 一种无人机航摄站点坐标求解装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6564146B1 (en) * 2000-01-24 2003-05-13 The United States Of America As Represented By The Secretary Of The Navy Tracking system for providing position information

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101201248A (zh) * 2007-12-07 2008-06-18 中国科学院武汉岩土力学研究所 基于无人机的航空近景摄影位移测量系统及其方法
CN102937443A (zh) * 2012-01-13 2013-02-20 唐粮 一种基于无人机的目标快速定位系统及方法
CN105182384A (zh) * 2015-08-24 2015-12-23 桂林电子科技大学 一种双模实时伪距差分定位系统和伪距改正数据生成方法
CN106093987A (zh) * 2016-08-10 2016-11-09 深圳供电局有限公司 一种应用于无人机的低成本差分gps系统及其实现方法
CN206479174U (zh) * 2017-02-26 2017-09-08 西安大地测绘股份有限公司 一种无人机航摄站点坐标求解装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
GPS PPK技术在像控点测量中的应用及精度分析;郭保;《测绘与空间地理信息》;20141130;第37卷(第11期);全文 *
GPS数据事后处理方法研究;胥辉旗 等;《海军航空工程学院学报》;20060331;第21卷(第2期);第245页右栏第1段-第246页右栏倒数第1段，图1 *
利用事后广域差分GPS方法处理北极GPS数据;鄂栋臣 等;《武汉测绘科技大学学报》;20001031;第25卷(第5期);全文 *

Also Published As

Publication number Publication date
CN106643670A (zh) 2017-05-10

Similar Documents

Publication Publication Date Title
CN106643670B (zh) 一种无人机航摄站点坐标求解装置及方法
US11221216B2 (en) Placement table for unmanned aerial vehicle, surveying method, surveying device, surveying system and program
Hol Sensor fusion and calibration of inertial sensors, vision, ultra-wideband and GPS
CN105974940B (zh) 适用于飞行器的目标跟踪方法
EP3454008A1 (en) Survey data processing device, survey data processing method, and survey data processing program
CN101241011B (zh) 激光雷达平台上高精度定位、定姿的装置和方法
CN106454879A (zh) 基于无人机航拍确定无线基站工参的方法及系统
CN104268935A (zh) 一种基于特征的机载激光点云与影像数据融合系统及方法
CN107449444A (zh) 一种多星图姿态关联的星敏感器内参数标定方法及其装置
CN111649737B (zh) 一种面向飞机精密进近着陆的视觉-惯性组合导航方法
CN104237922A (zh) 一种gnss/imu一体化无人机测绘的方法和系统
CN108562279A (zh) 一种无人机测绘方法
KR20190049086A (ko) 무인기 고속비행 중 정밀위치 영상 획득 장치 및 그를 이용한 정밀위치 획득 방법
WO2021007855A1 (zh) 基站、像控点定位方法、电子设备和计算机可读介质
CN110986888A (zh) 一种航空摄影一体化方法
CN109341685B (zh) 一种基于单应变换的固定翼飞机视觉辅助着陆导航方法
CN116123998A (zh) 多站点基于视频采集对空中炸点进行实时测量方法
CN109975836B (zh) 解算ccd图像地面位置方法、装置、电子设备及介质
CN105424034A (zh) 一种船载全天时星光惯导组合导航系统
CN106353756A (zh) 基于图像匹配的下降轨聚束合成孔径雷达定位方法
KR100469801B1 (ko) 실시간 항측 기준점 측량시스템 및 방법
CN107607933B (zh) 一种基于俯仰角迭代的远程无人机高精度定位方法
CN110927765B (zh) 激光雷达与卫星导航融合的目标在线定位方法
Mostafa Boresight calibration of integrated inertial/camera systems
CN109946682A (zh) 基于ICESat/GLAS的GF3数据基线估计方法

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Chen Gang

Inventor after: Zhang Fujun

Inventor after: Zhang Lili

Inventor after: Mao Zhihui

Inventor after: Yue Bo

Inventor after: Qi Jianzhu

Inventor after: Zhao Ying

Inventor after: Luo Ping

Inventor after: Zhang Jiao

Inventor before: Zhang Lili

Inventor before: Mao Zhihui

Inventor before: Chen Gang

Inventor before: Yue Bo

Inventor before: Qi Jianzhu

Inventor before: Zhao Ying

Inventor before: Zhang Fujun

Inventor before: Luo Ping

Inventor before: Zhang Jiao