CN110456397B - 一种多天线超短基线定位监测方法、装置及存储介质 - Google Patents
一种多天线超短基线定位监测方法、装置及存储介质 Download PDFInfo
- Publication number
- CN110456397B CN110456397B CN201910691238.6A CN201910691238A CN110456397B CN 110456397 B CN110456397 B CN 110456397B CN 201910691238 A CN201910691238 A CN 201910691238A CN 110456397 B CN110456397 B CN 110456397B
- Authority
- CN
- China
- Prior art keywords
- carrier phase
- pseudo
- correction
- range
- change rate
- 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.)
- Expired - Fee Related
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000012937 correction Methods 0.000 claims abstract description 188
- 238000012545 processing Methods 0.000 claims abstract description 14
- 230000008859 change Effects 0.000 claims description 72
- 230000008569 process Effects 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000012806 monitoring device Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 description 9
- 230000003068 static effect Effects 0.000 description 7
- 238000004891 communication Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000000819 phase cycle Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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
- 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/40—Correcting position, velocity or attitude
- G01S19/41—Differential correction, e.g. DGPS [differential GPS]
-
- 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
- G01S19/421—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system
-
- 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
- G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
- G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供一种多天线超短基线定位监测方法、装置及存储介质，方法为：通过预先部署的多个接收机接收卫星信号，将任意一个接收机作为基准站，得到基准站的原始观测数据，利用原始观测数据中的基准站精确位置计算得到载波相位差分改正数和伪距差分改正数；根据差分改正数和各个接收机的所在天线位置观测数据得到各个接收机的天线坐标；能够快速地切换基准站，得到接收机间相对位置，达到相互监测目的；根据各个接收机的天线所在位置对目标点进行定位；在不同方位部署多个接收机，通过基准站的定位信息得到载波相位差分改正数，确定多个接收机的天线坐标，从而对目标点进行定位监测，监测的范围较大，便于监测信息收集和处理。
Description
技术领域
本发明主要涉及卫星定位处理技术领域，具体涉及一种多天线超短基线定位监测方法、装置及存储介质。
背景技术
目前，形变监测装置多采用事后差分解算或者采用单天线RTK解算，采用事后差分解算可利用IGS发布的误差信息，虽然能提高一些误差修正，但通常需要较长的等待时间，事后差分解算或者采用单天线RTK解算不利于对形变的实时监测，对于滑坡等地质形变灾害来说实时性难以满足，且事后差分解算所能监测的区域较小，只能针对天线所在点或很小区域进行监测，不便于监测信息收集和处理。
发明内容
本发明所要解决的技术问题是针对现有技术的不足，提供一种多天线超短基线定位监测方法、装置及存储介质。
本发明解决上述技术问题的技术方案如下：一种多天线超短基线定位监测方法，包括如下步骤：
通过预先部署的多个接收机接收卫星信号，根据芯片输出协议将各个接收机接收的卫星信号进行解析，得到各个接收机对应的所在天线位置观测数据，将任意一个接收机作为基准站，所述基准站的所在天线位置观测数据作为原始观测数据；
在相邻两个观测时间点分别根据所述原始观测数据得到所述基准站的相邻两个载波相位差分改正数和相邻两个伪距差分改正数；
计算相邻两个观测时间点之间的间隔时间，根据所述相邻两个载波相位差分改正数以及所述间隔时间得到载波相位差分改正数变化率，并根据所述相邻两个伪距差分改正数以及所述间隔时间得到伪距差分改正数变化率；
根据所述载波相位差分改正数变化率和所述伪距差分改正数变化率得到共同观测卫星改正数；
根据所述共同观测卫星改正数和各接收机对应的所在天线位置观测数据分别得到各个接收机的天线坐标；
根据所述各个接收机的天线坐标对目标点进行定位监测。
本发明解决上述技术问题的另一技术方案如下：一种多天线超短基线定位监测装置，包括：
信号解析模块，用于通过预先部署的多个接收机接收卫星信号，根据芯片输出协议将各个接收机接收的卫星信号进行解析，得到各个接收机对应的所在天线位置观测数据，将任意一个接收机作为基准站，所述基准站的所在天线位置观测数据作为原始观测数据；
处理模块，用于在相邻两个观测时间点分别根据所述原始观测数据得到所述基准站的相邻两个载波相位差分改正数和相邻两个伪距差分改正数；
计算相邻两个观测时间点之间的间隔时间，根据所述相邻两个载波相位差分改正数以及所述间隔时间得到载波相位差分改正数变化率，并根据所述相邻两个伪距差分改正数以及所述间隔时间得到伪距差分改正数变化率；
根据所述载波相位差分改正数变化率和伪距差分改正数变化率得到共同观测卫星改正数；
根据所述共同观测卫星改正数和各接收机对应的所在天线位置观测数据分别得到各个接收机的天线坐标；
监测模块，用于根据所述各个接收机的天线坐标对目标点进行定位监测。
本发明解决上述技术问题的另一技术方案如下：一种多天线超短基线定位监测装置，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的嵌入式程序，当所述处理器执行所述嵌入式程序时，实现如上所述的多天线超短基线定位监测方法。
本发明解决上述技术问题的另一技术方案如下：一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，当所述计算机程序被处理器执行时，实现如上所述的多天线超短基线定位监测方法。
本发明的有益效果是：在不同方位部署多个接收机，将任意一个接收机作为基准站，能够快速地切换基准站，通过基准站的定位信息得到载波相位差分改正数，即得到接收机间相对位置，从而确定多个接收机的天线坐标，通过多个天线坐标对目标点进行定位监测，达到相互监测目的，监测的范围较大，降低成本，便于监测信息收集和处理。
附图说明
图1为本发明实施例提供的多天线超短基线定位监测方法的流程示意图；
图2为本发明实施例提供的多天线超短基线定位监测装置的模块框图；
图3为本发明实施例提供的多个接收机的布局图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
如图1所示，一种多天线超短基线形变监测方法，包括如下步骤：
通过预先部署的多个接收机接收卫星信号，根据芯片输出协议将各个接收机接收的卫星信号进行解析，得到各个接收机对应的所在天线位置观测数据，将任意一个接收机作为基准站，所述基准站的所在天线位置观测数据作为原始观测数据；
在相邻两个观测时间点分别根据所述原始观测数据得到所述基准站的相邻两个载波相位差分改正数和相邻两个伪距差分改正数；
计算相邻两个观测时间点之间的间隔时间，根据所述相邻两个载波相位差分改正数以及所述间隔时间得到载波相位差分改正数变化率，并根据所述相邻两个伪距差分改正数以及所述间隔时间得到伪距差分改正数变化率；
根据所述载波相位差分改正数变化率和所述伪距差分改正数变化率得到共同观测卫星改正数；
根据所述共同观测卫星改正数和各接收机对应的所在天线位置观测数据分别得到各个接收机的天线坐标；
根据所述各个接收机的天线坐标对目标点进行定位监测。
具体地，所述接收机上设有天线(GNSS天线)，天线和接收机是通过信号线连接。多个接收机同时挂载在CPU上，因此解算的位置是天线所在位置，形成的监测是对天线所在位置的。原始观测数据包括观测文件、星历文件、载波相位、伪距、多普勒等观测值。
应理解地，天线位置观测数据与接收机观测数据其实是同一种数据，但该数据解算出来的位置是天线所在位置即天线坐标。
上述实施例中，在不同方位部署多个接收机，将任意一个接收机作为基准站，能够快速地切换基准站，通过基准站的定位信息得到载波相位差分改正数，即得到接收机间相对位置，从而确定多个接收机的天线坐标，通过多个天线坐标对目标点进行定位监测，达到相互监测目的，监测的范围较大，降低成本，便于监测信息收集和处理。
可选地，作为本发明的一个实施例，所述得到所述基准站的相邻两个载波相位差分改正数和相邻两个伪距差分改正数的过程包括：
在相邻两个观测时间点分别得到所述基准站的第一原始观测数据和第二原始观测数据；
根据所述第一原始观测数据得到第一载波相位观测值和第一伪距观测值，根据所述第二原始观测数据得到第二载波相位观测值和第二伪距观测值；
从IGS情报采集系统中获取精密星历数据和误差改正数据，根据精密单点定位法、获取的所述精密星历数据和所述误差改正数据以及各个接收机的所在天线位置观测数据，对所述各个接收机进行位置标定，得到所述各个接收机的标定数据；
根据所述各个接收机的标定数据分别对所述第一载波相位观测值、所述第一伪距观测值、所述第二载波相位观测值和所述第二伪距观测值进行解算，分别得到第一载波相位差分改正数、第二伪距差分改正数、第二载波相位差分改正数和第二伪距差分改正数，其中，所述第一载波相位差分改正数和第二载波相位差分改正数为所述相邻两个载波相位差分改正数，所述第二载波相位差分改正数和第二载波相位差分改正数为所述相邻两个伪距差分改正数。
具体地，利用原始观测数据中的基准站精确位置来得到载波相位观测值和伪距观测值。
上述实施例中，由于互差分需要不同天线的精确位置，因此在布置天线后需要对每个天线位置进行精密单点定位(precise point positioning,PPP)，利用其中一台接收机(GNSS)的载波相位观测值和伪距观测值，采用高精度的卫星轨道和钟差产品，并通过模型改正或参数估计的方法精细考虑与卫星端、信号传播路径及接收机端有关误差对定位的影响，实现高精度定位。
可选地，作为本发明的一个实施例，所述根据所述相邻两个载波相位差分改正数以及所述间隔时间得到载波相位差分改正数变化率的过程包括：
将所述相邻两个载波相位差分改正数进行差值计算，得到载波相位数值差值；
将所述载波相位数值差值与所述间隔时间进行比值计算，得到载波相位差分改正数变化率；
所述根据所述相邻两个伪距差分改正数以及所述间隔时间得到伪距差分改正数变化率的过程包括：
将所述相邻两个伪距差分改正数进行差值计算，得到伪距相位数值差值；
将所述伪距相位数值差值与所述间隔时间进行比值计算，得到伪距相位差分改正数变化率。
上述实施例中，还包括如下步骤：
对卫星钟差进行修正。
在对载波相位观测值进行伪距平滑滤波处理前还包括载波相位观测值进行载波相位周跳探测及修复。
具体地，计算所述基准站和所述卫星的站星距离值利用精密坐标算法计算得到。
上述实施例中，伪距观测噪声约为30cm，而相位观测噪声为(1～2)mm，相差两个量级，利用接收机伪距信息，使得载波相位整周模糊度的固定提供一定的收索范围，减小在收索双差整周模糊度时所耗费的CPU资源，能够更有效的实现多接收机的定位解算，有效地提高各个接收机天线观测数据的解算。
上述实施例中，通过平滑滤波处理能够消除伪距噪声，得到相对准确的伪距值，将基准站的观测值与PPP标定位置进行解算，得到差分改正数，差分改正数变化率。通过选择不同接收机作为基准站，实现各监测接收机之间的相对高精度监测定位，达到整个面的形变监测。而且可将此部分接收机作为其他模组的基准站，实现监测站之间的相对监测，防止出现整个面的相变而不能有效监测的困境，实现设备的监测完整性。
可选地，作为本发明的一个实施例，根据所述载波相位差分改正数变化率和伪距差分改正数变化率得到共同观测卫星改正数的过程包括：
将当前观测时间点得到的载波相位差分改正数和伪距差分改正数以及所述载波相位差分改正数变化率和所述伪距差分改正数变化率下发至各个所述接收机中；
依次将各个所述接收机中的载波相位差分改正数变化率和所述伪距差分改正数变化率分别进行观测周及周内秒比对解算，得到解算结果，如果所述解算结果在设定范围内，则采用所述载波相位差分改正数变化率及所述伪距差分改正数变化率，否则舍弃该数据，进行下一轮差分解算；
将得到的伪距相位数值差值与所述载波相位差分改正数变化率及所述伪距差分改正数变化率相乘，将得到的乘积加上当前观测时间点得到的载波相位差分改正数和伪距差分改正数，得到所有接收机的共同观测卫星改正数。
具体地，在上述执行过程中，根据动静态载波相位差分算法来判断执行过程是否顺利，其过程为：
静态固定：载波相位双差静态算法固定载波整周模糊度。
动态判断：静态情况下根据设置的动态阈值R判断是否进入动态算法。
动态监测：载波相位双差动态算法固定载波整周模糊度。
形变监测：判断是否产生形变。
动态跟踪：利用动态算法跟踪形变。
静态判断：形变终止后根据设置的静态阈值R判断是否进入静态算法，其中，R为整周模糊度第一最佳解和第二最佳解公式比值，具体公式为
此算法用于实时得形变跟踪，具有跟踪性好，实时性强，精度高等优点。
上述实施例中，通过匹配算法对载波相位差分改正数的误差的公共量进行消除。即对接收的载波相位差分改正数的数据帧进行解析，匹配相同卫星的卫星号，匹配观测周数判断这一帧数据是否在同一个观测周内，匹配观测周内秒是否在同一个观测历元内，综合以上条件得到共同观测卫星改正数。
上述实施例中，对载波相位差分改正数的误差进行消除得到共同观测卫星改正数，能够使各天线实现相对高精度定位。
可选地，作为本发明的一个实施例，在下发所述载波相位差分改正数变化率和所述伪距差分改正数变化率前还包括步骤：
将所述基准站对应的卫星作为各个接收机共同的公共观测卫星；
将各个接收机与所述公共观测卫星的卫星号及卫星系统进行匹配，以确定不同接收机所需的载波相位差分改正数和伪距差分改正数。
可选地，作为本发明的一个实施例，所述根据所述共同观测卫星改正数得到各个接收机的天线坐标的过程包括：
根据第一式进行计算任一接收机的天线坐标，所述第一式为
A[x’k y’k z’k]-1＝[x’s y’s z’s]-1+[Δx’ Δy’ Δz’]-1，
其中，[x’k y’k z’k]-1为天线坐标，[x’s y’s z’s]-1为接收机天线所在位置解算坐标，所述接收机天线所在位置解算坐标从接收机对应的所在天线位置观测数据中获得，[Δx’ Δy’ Δz’]-1为共同观测卫星改正数，根据解算坐标消除误差。
应理解地，第一式仅表示解算示意，不代表实际解算程序。
上述实施例中是采用相对定位方式进行计算，上式中[Δx’ Δy’ Δz’]-1分量是准确的，即基准站与观测点站之间的坐标距离(除基准站外其他各接收站为观测点站)，因此，所得观测点站的监测点坐标相对于基准站坐标是严格准确的。接着，切换接收机作为基准站，实现对上一基准站的相对高精度定位。
上述实施例中，能够快速得到各个天线坐标，所得监测点坐标相对于基准站坐标是严格准确的，切换接收机作为基准站，实现对上一基准站的相对高精度定位。
由于在安装设备时需要获取天线所在的绝对、高精度的位置，需要对每个天线进行PPP定位。由于精密星历的延时发布，在后续计算基站站中使用精密单点定位(precisepoint positioning,PPP)技术，提高基准站坐标的精度，使得监测精度进一步提升，可以将观测数据传输到云端服务器，方便安装过后再对任意一接收机(GNSS)对应的天线进行定位解算。
防止某一区域整体的偏移形变，可以使用两组多天线短基线监测组件，进行相互监测，即实现组建整体模块化。在实际应用中可视具体情况选择数据链路，如将载波相位差分改正数上传至服务器端，各组监测点访问获取改正数；在短基线情况下使用有线链路，实时获取共同观测卫星改正数，提高监测精度；多检测点情况下采用无线广播方式等等。模块化组建方式使得监测点布置更加便捷，节约成本，有效提高监测精度，使得满足毫米级的微变监测功能。
图2为本发明实施例提供的监测装置的模块框图。
可选地，作为本发明的一个实施例，如图2所示，一种多天线超短基线形变监测装置，包括：
信号解析模块，用于通过预先部署的多个接收机接收卫星信号，根据芯片输出协议将各个接收机接收的卫星信号进行解析，得到各个接收机对应的所在天线位置观测数据，将任意一个接收机作为基准站，所述基准站的所在天线位置观测数据作为原始观测数据；
处理模块，用于在相邻两个观测时间点分别根据所述原始观测数据得到所述基准站的相邻两个载波相位差分改正数和相邻两个伪距差分改正数；
计算相邻两个观测时间点之间的间隔时间，根据所述相邻两个载波相位差分改正数以及所述间隔时间得到载波相位差分改正数变化率，并根据所述相邻两个伪距差分改正数以及所述间隔时间得到伪距差分改正数变化率；
根据所述载波相位差分改正数变化率和伪距差分改正数变化率得到共同观测卫星改正数；
根据所述共同观测卫星改正数和各接收机对应的所在天线位置观测数据分别得到各个接收机的天线坐标；
监测模块，用于根据所述各个接收机的天线坐标对目标点进行定位监测。
如图3所示，各个接收机通过处理模块与服务器端连接，将处理的定位监测数据上传至服务器或保存在存储介质中。所述信号解析模块根据芯片输出协议将所述基准站的卫星信号进行解析，处理模块与接收机的通讯方式可为串口通信、I2C通信、SPI通信，可减小对处理器硬件接口的要求。将得到所述基准站的原始观测数据传输至处理模块中，进行下一步解算处理，包括卫星匹配，差分解算定位，然后得到定位监测数据，该数据上传至服务器或保存在本地。由于差分改正数的大小会随着时间的推移而变化，而且监测站和基准站并不一定会在同一时刻对卫星信号进行观测，因此基准站需要将所有观测卫星数据计算上传服务器，将计算得到差分改正数及其变化率，进行相应的数据帧编码转化为二进制进行本地存储方便后续其他接收机数据差分解算。
根据所述各个接收机的天线坐标对目标点进行定位监测，具体为：确定各个接收机的天线坐标后，目标点向各个接收机发送目标信号，各个接收机通过天线坐标对该目标信号进行定位，从而得到目标点的定位信息。
可选地，作为本发明的另一个实施例，一种多天线超短基线形变监测装置，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的嵌入式程序，当所述处理器执行所述嵌入式程序时，实现如上所述的多天线超短基线形变监测方法。
可选地，作为本发明的另一个实施例，一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，当所述计算机程序被处理器执行时，实现如上所述的多天线超短基线定位监测方法。
将嵌入式程序解算结果储存，原始观测数据储存，用于其他定位分析及保存。
在本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。
作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本发明实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个嵌入式可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该嵌入式软件产品存储在一个存储介质中，包括若干指令用以使得一台嵌入式设备(可以是个人嵌入式，服务器，或者网络设备等)执行本发明各个实施例方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。
Claims (7)
1.一种多天线超短基线定位监测方法，其特征在于，包括如下步骤：
通过预先部署的多个接收机接收卫星信号，根据芯片输出协议将各个接收机接收的卫星信号进行解析，得到各个接收机对应的所在天线位置观测数据，将任意一个接收机作为基准站，所述基准站的所在天线位置观测数据作为原始观测数据；
在相邻两个观测时间点分别根据所述原始观测数据得到所述基准站的相邻两个载波相位差分改正数和相邻两个伪距差分改正数；
计算相邻两个观测时间点之间的间隔时间，根据所述相邻两个载波相位差分改正数以及所述间隔时间得到载波相位差分改正数变化率，并根据所述相邻两个伪距差分改正数以及所述间隔时间得到伪距差分改正数变化率；所述根据所述相邻两个载波相位差分改正数以及所述间隔时间得到载波相位差分改正数变化率的过程包括：
将所述相邻两个载波相位差分改正数进行差值计算，得到载波相位数值差值；
将所述载波相位数值差值与所述间隔时间进行比值计算，得到载波相位差分改正数变化率；
所述根据所述相邻两个伪距差分改正数以及所述间隔时间得到伪距差分改正数变化率的过程包括：
将所述相邻两个伪距差分改正数进行差值计算，得到伪距相位数值差值；
将所述伪距相位数值差值与所述间隔时间进行比值计算，得到伪距相位差分改正数变化率；
根据所述载波相位差分改正数变化率和所述伪距差分改正数变化率得到共同观测卫星改正数；根据所述载波相位差分改正数变化率和伪距差分改正数变化率得到共同观测卫星改正数的过程包括：
将当前观测时间点得到的载波相位差分改正数和伪距差分改正数以及所述载波相位差分改正数变化率和所述伪距差分改正数变化率下发至各个所述接收机中；
依次将各个所述接收机中的载波相位差分改正数变化率和所述伪距差分改正数变化率分别进行观测周及周内秒比对解算，得到解算结果，如果所述解算结果在设定范围内，则采用所述载波相位差分改正数变化率及所述伪距差分改正数变化率，否则舍弃该数据，进行下一轮差分解算；
将得到的载波相位数值差值与所述载波相位差分改正数变化率及所述伪距差分改正数变化率相乘，将得到的乘积加上当前观测时间点得到的载波相位差分改正数和伪距差分改正数，得到所有接收机的共同观测卫星改正数；
根据所述共同观测卫星改正数和各接收机对应的所在天线位置观测数据分别得到各个接收机的天线坐标；
根据所述各个接收机的天线坐标对目标点进行定位监测。
2.根据权利要求1所述的多天线超短基线定位监测方法，其特征在于，所述得到所述基准站的相邻两个载波相位差分改正数和相邻两个伪距差分改正数的过程包括：
在相邻两个观测时间点分别得到所述基准站的第一原始观测数据和第二原始观测数据；
根据所述第一原始观测数据得到第一载波相位观测值和第一伪距观测值，根据所述第二原始观测数据得到第二载波相位观测值和第二伪距观测值；
从IGS情报采集系统中获取精密星历数据和误差改正数据，根据精密单点定位法、获取的所述精密星历数据和所述误差改正数据以及各个接收机的所在天线位置观测数据，对所述各个接收机进行位置标定，得到所述各个接收机的标定数据；
根据所述各个接收机的标定数据分别对所述第一载波相位观测值、所述第一伪距观测值、所述第二载波相位观测值和所述第二伪距观测值进行解算，分别得到第一载波相位差分改正数、第二伪距差分改正数、第二载波相位差分改正数和第二伪距差分改正数，其中，所述第一载波相位差分改正数和第二载波相位差分改正数为所述相邻两个载波相位差分改正数，所述第二载波相位差分改正数和第二载波相位差分改正数为所述相邻两个伪距差分改正数。
3.根据权利要求1所述的多天线超短基线定位监测方法，其特征在于，在下发所述载波相位差分改正数变化率和所述伪距差分改正数变化率前还包括步骤：
将所述基准站对应的卫星作为各个接收机共同的公共观测卫星；
将各个接收机与所述公共观测卫星的卫星号及卫星系统进行匹配。
4.根据权利要求1至3任一项所述的多天线超短基线定位监测方法，其特征在于，所述根据所述共同观测卫星改正数得到各个接收机的天线坐标的过程包括：
根据第一式进行计算任一接收机的天线坐标，所述第一式为
[x’k y’k z’k]-1＝[x’s y’s z’s]-1+[Δx’ Δy’ Δz’]-1，
其中，[x’k y’k z’k]-1为天线坐标，[x’s y’s z’s]-1为接收机天线所在位置解算坐标，所述接收机天线所在位置解算坐标从接收机对应的所在天线位置观测数据中获得，[Δx’ Δy’ Δz’]-1为共同观测卫星改正数，根据解算坐标消除误差。
5.一种多天线超短基线定位监测装置，其特征在于，包括：
信号解析模块，用于通过预先部署的多个接收机接收卫星信号，根据芯片输出协议将各个接收机接收的卫星信号进行解析，得到各个接收机对应的所在天线位置观测数据，将任意一个接收机作为基准站，所述基准站的所在天线位置观测数据作为原始观测数据；
处理模块，用于在相邻两个观测时间点分别根据所述原始观测数据得到所述基准站的相邻两个载波相位差分改正数和相邻两个伪距差分改正数；
计算相邻两个观测时间点之间的间隔时间，根据所述相邻两个载波相位差分改正数以及所述间隔时间得到载波相位差分改正数变化率，并根据所述相邻两个伪距差分改正数以及所述间隔时间得到伪距差分改正数变化率；所述根据所述相邻两个载波相位差分改正数以及所述间隔时间得到载波相位差分改正数变化率的过程包括：
将所述相邻两个载波相位差分改正数进行差值计算，得到载波相位数值差值；
将所述载波相位数值差值与所述间隔时间进行比值计算，得到载波相位差分改正数变化率；
所述根据所述相邻两个伪距差分改正数以及所述间隔时间得到伪距差分改正数变化率的过程包括：
将所述相邻两个伪距差分改正数进行差值计算，得到伪距相位数值差值；
将所述伪距相位数值差值与所述间隔时间进行比值计算，得到伪距相位差分改正数变化率；
根据所述载波相位差分改正数变化率和所述伪距差分改正数变化率得到共同观测卫星改正数；根据所述载波相位差分改正数变化率和伪距差分改正数变化率得到共同观测卫星改正数的过程包括：
将当前观测时间点得到的载波相位差分改正数和伪距差分改正数以及所述载波相位差分改正数变化率和所述伪距差分改正数变化率下发至各个所述接收机中；
依次将各个所述接收机中的载波相位差分改正数变化率和所述伪距差分改正数变化率分别进行观测周及周内秒比对解算，得到解算结果，如果所述解算结果在设定范围内，则采用所述载波相位差分改正数变化率及所述伪距差分改正数变化率，否则舍弃该数据，进行下一轮差分解算；
将得到的载波相位数值差值与所述载波相位差分改正数变化率及所述伪距差分改正数变化率相乘，将得到的乘积加上当前观测时间点得到的载波相位差分改正数和伪距差分改正数，得到所有接收机的共同观测卫星改正数；
根据所述共同观测卫星改正数和各接收机对应的所在天线位置观测数据分别得到各个接收机的天线坐标；
监测模块，用于根据所述各个接收机的天线坐标对目标点进行定位监测。
6.一种多天线超短基线定位监测装置，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的嵌入式程序，其特征在于，当所述处理器执行所述嵌入式程序时，实现如权利要求1至4任一项所述的多天线超短基线定位监测方法。
7.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，当所述计算机程序被处理器执行时，实现如权利要求1至4任一项所述的多天线超短基线定位监测方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910691238.6A CN110456397B (zh) | 2019-07-29 | 2019-07-29 | 一种多天线超短基线定位监测方法、装置及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910691238.6A CN110456397B (zh) | 2019-07-29 | 2019-07-29 | 一种多天线超短基线定位监测方法、装置及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110456397A CN110456397A (zh) | 2019-11-15 |
| CN110456397B true CN110456397B (zh) | 2021-02-09 |
Family
ID=68483848
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910691238.6A Expired - Fee Related CN110456397B (zh) | 2019-07-29 | 2019-07-29 | 一种多天线超短基线定位监测方法、装置及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110456397B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110727008B (zh) * | 2019-11-19 | 2021-06-22 | 西南交通大学 | 一种gnss整周模糊度的快速解算方法 |
| CN112649818B (zh) * | 2020-11-20 | 2021-08-17 | 深圳市天弓导航科技有限公司 | 卫星导航接收机的检测方法、装置、终端设备及介质 |
| CN116482733B (zh) * | 2023-06-26 | 2023-09-08 | 银河航天(北京)通信技术有限公司 | 实现大范围应用载波相位技术的方法、装置及存储介质 |
| CN116577816B (zh) * | 2023-07-12 | 2023-09-19 | 银河航天(北京)网络技术有限公司 | 利用双差定位确定用户站位置的方法、装置及存储介质 |
Citations (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1995018978A1 (en) * | 1994-01-03 | 1995-07-13 | Trimble Navigation | A network for carrier phase differential gps corrections |
| WO2002061449A1 (en) * | 2000-12-15 | 2002-08-08 | Trimble Navigation Limited | Method and apparatus for correcting the delay error of gps signals |
| CN102414577A (zh) * | 2009-05-02 | 2012-04-11 | 天宝导航有限公司 | Gnss信号处理方法和设备 |
| JP2012112738A (ja) * | 2010-11-24 | 2012-06-14 | Nec Corp | 追尾装置及び追尾方法 |
| CN203011355U (zh) * | 2012-12-10 | 2013-06-19 | 陕西西北铁道电子有限公司 | 基于北斗兼容型的铁路桥梁状态监测系统 |
| CN104062895A (zh) * | 2014-06-26 | 2014-09-24 | 桂林电子科技大学 | 伪卫星时间同步及其定位方法 |
| CN104536027A (zh) * | 2015-01-28 | 2015-04-22 | 中国人民解放军国防科学技术大学 | 一种实时北斗精密相对定位的方法 |
| CN104536026A (zh) * | 2015-01-08 | 2015-04-22 | 中国航空无线电电子研究所 | 一种动态对动态实时测量系统 |
| CN104570024A (zh) * | 2014-12-22 | 2015-04-29 | 中国科学院国家授时中心 | 北斗天基高精度实时定位方法 |
| CN105319566A (zh) * | 2014-06-24 | 2016-02-10 | 北京大学 | 一种协同定位方法及装置 |
| CN105929424A (zh) * | 2015-11-16 | 2016-09-07 | 中国测绘科学研究院 | 一种bds/gps高精度定位方法 |
| CN107479066A (zh) * | 2017-09-14 | 2017-12-15 | 北方信息控制研究院集团有限公司 | 一种步进式移动地基增强方法 |
| CN107607032A (zh) * | 2017-10-23 | 2018-01-19 | 上海普适导航科技股份有限公司 | 一种gnss形变监测系统 |
| CN108317949A (zh) * | 2018-02-07 | 2018-07-24 | 桂林电子科技大学 | 一种rtk高精度差分定位形变监测系统及方法 |
| CN109085617A (zh) * | 2018-08-29 | 2018-12-25 | 桂林电子科技大学 | 一种gnss监测站的定位系统及定位方法 |
| CN109655846A (zh) * | 2019-01-30 | 2019-04-19 | 清华大学 | 一种多站差分后处理高精度时间同步方法及系统 |
| CN109782305A (zh) * | 2018-12-28 | 2019-05-21 | 中国电子科技集团公司第二十研究所 | 多基准条件下接收机优选策略生成rtk报文的方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001281317A (ja) * | 2000-03-30 | 2001-10-10 | Japan Radio Co Ltd | 移動体姿勢計測システムおよび複数アンテナ配置構造 |
| CN105044750A (zh) * | 2015-07-10 | 2015-11-11 | 上海北斗卫星导航平台有限公司 | 面向大众应用的高精度连续定位系统和方法 |
| CN105116427A (zh) * | 2015-08-31 | 2015-12-02 | 无锡伊佩克科技有限公司 | 一种基于预失真器的车载导航信道自适应系统 |
| CN105353392A (zh) * | 2015-10-30 | 2016-02-24 | 中国石油大学(华东) | 基于gnss多天线的动态载体精密定位方法 |
| CN106680854A (zh) * | 2017-01-17 | 2017-05-17 | 桂林电子科技大学 | 一种低成本高精度定位系统及方法 |
| CN109541663B (zh) * | 2018-11-12 | 2022-04-05 | 华东师范大学 | 一种gnss定位测姿多路径误差的纠正方法 |
-
2019
- 2019-07-29 CN CN201910691238.6A patent/CN110456397B/zh not_active Expired - Fee Related
Patent Citations (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1995018978A1 (en) * | 1994-01-03 | 1995-07-13 | Trimble Navigation | A network for carrier phase differential gps corrections |
| WO2002061449A1 (en) * | 2000-12-15 | 2002-08-08 | Trimble Navigation Limited | Method and apparatus for correcting the delay error of gps signals |
| CN102414577A (zh) * | 2009-05-02 | 2012-04-11 | 天宝导航有限公司 | Gnss信号处理方法和设备 |
| JP2012112738A (ja) * | 2010-11-24 | 2012-06-14 | Nec Corp | 追尾装置及び追尾方法 |
| CN203011355U (zh) * | 2012-12-10 | 2013-06-19 | 陕西西北铁道电子有限公司 | 基于北斗兼容型的铁路桥梁状态监测系统 |
| CN105319566A (zh) * | 2014-06-24 | 2016-02-10 | 北京大学 | 一种协同定位方法及装置 |
| CN104062895A (zh) * | 2014-06-26 | 2014-09-24 | 桂林电子科技大学 | 伪卫星时间同步及其定位方法 |
| CN104570024A (zh) * | 2014-12-22 | 2015-04-29 | 中国科学院国家授时中心 | 北斗天基高精度实时定位方法 |
| CN104536026A (zh) * | 2015-01-08 | 2015-04-22 | 中国航空无线电电子研究所 | 一种动态对动态实时测量系统 |
| CN104536027A (zh) * | 2015-01-28 | 2015-04-22 | 中国人民解放军国防科学技术大学 | 一种实时北斗精密相对定位的方法 |
| CN105929424A (zh) * | 2015-11-16 | 2016-09-07 | 中国测绘科学研究院 | 一种bds/gps高精度定位方法 |
| CN107479066A (zh) * | 2017-09-14 | 2017-12-15 | 北方信息控制研究院集团有限公司 | 一种步进式移动地基增强方法 |
| CN107607032A (zh) * | 2017-10-23 | 2018-01-19 | 上海普适导航科技股份有限公司 | 一种gnss形变监测系统 |
| CN108317949A (zh) * | 2018-02-07 | 2018-07-24 | 桂林电子科技大学 | 一种rtk高精度差分定位形变监测系统及方法 |
| CN109085617A (zh) * | 2018-08-29 | 2018-12-25 | 桂林电子科技大学 | 一种gnss监测站的定位系统及定位方法 |
| CN109782305A (zh) * | 2018-12-28 | 2019-05-21 | 中国电子科技集团公司第二十研究所 | 多基准条件下接收机优选策略生成rtk报文的方法 |
| CN109655846A (zh) * | 2019-01-30 | 2019-04-19 | 清华大学 | 一种多站差分后处理高精度时间同步方法及系统 |
Non-Patent Citations (5)
| Title |
|---|
| 《Design and Implementation of Dual-Mode Pseudorange Differential Positioning》;Chenglin Cai; Jia Hu;《2019 IEEE International Conference on Power, Intelligent Computing and Systems (ICPICS)》;20190714;283-286 * |
| 《Multi-monitor station "on the fly" ambiguity resolution: the impacts of satellite geometry and monitor station geometry》;H.Z. Abidin;《IEEE PLANS 92 Position Location and Navigation Symposium Record》;20020806;412-418 * |
| 《基于多参考站的分米级GPS伪距差分定位方法》;邓健;《东南大学学报(自然科学版)》;20100310;316-319 * |
| 《差分GPS基准站测量方法比较分析》;郁聪冲;《价值工程》;20160418;93-94 * |
| 《广域差分新方法的定位性能与差分网优化布局》;蔡成林; 李孝辉; 吴海涛;《宇航学报》;20090730;1404-1409 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110456397A (zh) | 2019-11-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110456397B (zh) | 一种多天线超短基线定位监测方法、装置及存储介质 | |
| CN107102346B (zh) | 一种基于北斗系统的多天线测姿方法 | |
| CN108931915B (zh) | 利用导航卫星的授时方法和装置、计算机可读存储介质 | |
| CN111045034B (zh) | 基于广播星历的gnss多系统实时精密时间传递方法及系统 | |
| JP2014510260A (ja) | 時計修正値を特定するための方法およびシステム | |
| JP5352422B2 (ja) | 測位装置及びプログラム | |
| CN108195283B (zh) | 基于无缝变换的gnss实时变形监测方法及系统 | |
| CN104714244A (zh) | 一种基于抗差自适应Kalman滤波的多系统动态PPP解算方法 | |
| CN112698563A (zh) | 卫星授时方法、装置、电子设备和存储介质 | |
| EP1844349A2 (en) | A system and method for monitoring and surveying movements of the terrain, large infrastructures and civil building works in general, based upon the signals transmitted by the gps navigation satellite system | |
| CN104459722B (zh) | 一种基于多余观测分量的整周模糊度可靠性检验方法 | |
| US20240255649A1 (en) | Network rtk anti-ionospheric disturbance positioning method, apparatus, system, device and storage medium | |
| CN113325446B (zh) | 一种多模共频gnss载波相位时间传递方法及系统 | |
| CN112146557A (zh) | 一种基于gnss的实时桥梁变形监测系统及方法 | |
| CN110824505B (zh) | Gnss卫星接收机的偏差估计方法及系统、定位方法及终端 | |
| CN109597099A (zh) | 判断基准站接收机是否移动的方法、oem板卡及接收机 | |
| EP4421534A1 (en) | Data preprocessing method, data preprocessing apparatus, and chip | |
| CN110736963A (zh) | 一种基于CSI的室内Wi-Fi定位方法、装置及存储介质 | |
| CN110068848B (zh) | 一种高性能rtk处理技术方法 | |
| CN116125371A (zh) | 一种卫星定向方法、装置、卫星导航芯片及存储介质 | |
| CN117853594A (zh) | 海底阵潜标位置标定方法、装置、存储介质及电子设备 | |
| CN112230254B (zh) | 一种gps载波相位多径误差的校正方法及装置 | |
| CN110058274B (zh) | 一种卫星导航系统间的时差监测方法及系统 | |
| CN114152961B (zh) | 一种导航系统的周跳处理方法和装置 | |
| CN115327595A (zh) | 基于gnss差分的航天器相对位置测量方法及其装置 |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210209 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |