CN113779817B - 一种测量控制网基准稳定性分析方法 - Google Patents
一种测量控制网基准稳定性分析方法 Download PDFInfo
- Publication number
- CN113779817B CN113779817B CN202111333451.3A CN202111333451A CN113779817B CN 113779817 B CN113779817 B CN 113779817B CN 202111333451 A CN202111333451 A CN 202111333451A CN 113779817 B CN113779817 B CN 113779817B
- Authority
- CN
- China
- Prior art keywords
- poor
- station
- baseline
- stability
- measuring
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明公开了一种测量控制网基准稳定性分析方法。它包括如下步骤,步骤一:统计基线较差;步骤二:计算较差标称中误差;步骤三:计算基线标准化较差;步骤四:构造检验统计量;步骤五:假设检验,判断测站稳定性;步骤六:结合网型,确定平差起算基准;分析测量控制网网形,从上述步骤得到的稳定测站中选取2‑4个测站为起算基准,进行后续平差处理。本发明克服了基于相对测量的工程控制网建设方法中现有起算基准稳定性分析方法严重依赖经验和先验信息、偶然性强、理论不严密、成果不可靠、易导致分析结论错误等问题;具有理论严密、适应性广、操作性强、可靠性高的优点。
Description
技术领域
本发明涉及一种测量控制网基准稳定性分析方法,更具体地说它是一种测量控制网数据处理中起算基准稳定性分析的方法。
背景技术
测量控制网的建设是工程项目勘察、设计、施工、运维阶段的先导工程;按照相关规范要求,测量控制网需要定期复测,以保证控制成果的精度和现势性;根据测量原理不同,工程控制网测量可以分为基于绝对测量的方法和基于相对测量的方法;基于绝对测量的方法一般指GNSS精密单点定位法(Precise Point Positioning,PPP),包括后处理PPP和实时PPP等;PPP测量基于GNSS***精密星历、精密钟差等产品,实现控制点的实时/事后绝对定位,无需平差解算即得到控制成果,但存在数据处理模型复杂、收敛时间较长、可用数据处理平台不多等问题,目前一般用于科研领域;基于相对测量的方法一般指传统边角测量法、GNSS相对定位测量等,通过观测控制点间的基线、边长、角度等,经过平差处理,得到控制成果;基于相对测量的控制网建设方法具有精度高、数据处理模型简单、数据处理平台成熟等一系列优点,是目前控制网建设的主要方法,但在平差过程中对起算基准的要求比较严格,起算基准之间不兼容,或者起算基准存在变形,会导致基准误差传播至控制点成果,得到错误的变形结论;
传统的测量基准稳定性分析方法包括“加点法”、“减点法”等;以“加点法”为例,首先结合经验和先验信息选取2个起算基准点进行平差计算,根据前、后两期变形结论和点位中误差反复迭代剔除不稳定测站,得到稳定点组,再综合考虑网型确定起算基准,重新执行平差计算,得到最终成果;上述方法严重依赖首批基准点的稳定性,如果首批基准点选择不合适,或者网中变形控制点较多,极易导致错误的基准稳定性分析结论,存在偶然性强、理论不严密、成果不可靠等问题,已经成为测量控制网建设的制约因素,影响工程建设的安全性和可靠性,急需提出理论严密、可靠性高、操作性强的测量控制网基准稳定性分析方法,满足工程建设对于测量控制网精度和可靠性的迫切需求。
发明内容
本发明的目的是为了提供一种测量控制网基准稳定性分析方法,理论严密、适应性广、操作性强、可靠性高;克服了基于相对测量的工程控制网建设方法中现有起算基准稳定性分析方法严重依赖经验和先验信息、偶然性强、理论不严密、成果不可靠、易导致分析结论错误等问题;本发明针对的是起算点的稳定性分析方法,确定起算点稳定后才进行其他基准点的稳定性分析;克服了现有变形监测网直接认为起算点稳定,未考虑起算点是否兼容的问题。
为了实现上述目的,本发明的技术方案为:一种测量控制网基准稳定性分析方法,其特征在于:包括如下步骤,
步骤一:统计基线较差;
统计待检验测站本期与周边其他测站组成的所有基线长度与前期基线长度之差;
步骤二:计算较差标称中误差;
根据设备参数计算较差标称中误差,用于步骤三计算基线标准化较差;
步骤三:计算基线标准化较差;
采用基线较差归一化处理,获得基线的标准化较差;
步骤四:构造检验统计量;
构造检验统计量,以用于步骤五的假设检验;
步骤五:假设检验,判断测站稳定性;
若上述方法检验结果不符合实际情况,则采用下述方法,以进行两种方法之间的互相检核。首先计算测站稳定性指数,根据测站稳定性指数迭代判断测站稳定性,确定稳定测站,进入下一步;
步骤六:结合网形,确定平差起算基准;
分析测量控制网网形,从上述步骤得到的稳定测站中选取2-4个能控制整网的测站为起算基准,进行后续平差处理。
在上述技术方案中,在步骤二中,计算较差标称中误差,具体包括如下步骤,
根据所采用设备标称精度和基线长度,计算各条基线较差的标称中误差;基线较差的标称中误差根据误差传播定律计算:
在上述技术方案中,在步骤三中,各条基线前后两期较差除以对应标称中误差,称为基线较差归一化处理,其结果称为基线的标准化较差;
基线较差归一化处理的方法如下,
其中,U为检验统计量;为与某一测站相关的基线标准化较差的均值;为待检验测站相关基线条数;为基线标准化较差;为标准差,此处指标准正态分布的标准差,为0,无量纲;为分布的期望,此处指标准正态分布的期望,为1,无量纲。
上述技术方案中,在步骤五中,在实际情况中,若统计量不满足标准正态分布,则分析结果会出现不符合实际的情况,如结果显示全部点位都是稳定点或全部点位都是不稳定点;此时采用构造稳定性指数的方法进行检核:
2)计算每个测站的稳定性指数;
当所有测站稳定性指数均小于1且比较一致,则认为测站均比较稳定,未发生明显变形;
当部分测站稳定性指数大于1且数值出现分群现象,则剔除对应最大稳定性指数的测站及对应基线,针对剩余测站和基线重复迭代执行步骤1)、2),直到剩余测站稳定性指数均小于1且比较一致为止;认为剔除的测站发生形变,不是稳定测站,不可作为起算基准;剩余测站为稳定测站;
本发明针对基线进行统计分析,可以进行起算点的稳定性分析,也可用于后续其他基准点的稳定性分析;克服了现有监测基准点的稳定性分析方法针对坐标分析,仅能在引入起算基准之后用于其他基准点的稳定性分析的缺点。
本发明具有如下优点:
针对目前工程控制网建设中传统起算基准稳定性分析方法严重依赖经验和先验信息、偶然性强、理论不严密、成果不可靠、易导致分析结论错误等问题,本发明通过理论研究和大量实践总结,提出了基于数理统计原理的新的测量控制网数据处理中起算基准稳定性分析方法;本发明方法具有理论严密、适应性广、操作性强、可靠性高等优点,可以在工程控制网建设、地质灾害监测、地壳形变监测等领域推广应用;克服了传统算法的上述多项缺陷。
附图说明
图1为本发明算法流程图。
具体实施方式
下面结合附图详细说明本发明的实施情况,但它们并不构成对本发明的限定,仅作举例而已。同时通过说明使本发明的优点更加清楚和容易理解。
参阅附图1可知:本发明所述的一种新的测量控制网基准稳定性分析方法,具体内容如下:
步骤一:统计待检验测站本期与周边其他测站组成的所有基线长度与前期基线长度之差;
步骤二:根据所采用设备标称精度和基线长度,计算各条基线较差的标称中误差;基线较差的标称中误差根据误差传播定律计算:
步骤三:各条基线前后两期较差除以对应标称中误差,称为基线较差归一化处理,其结果称为基线的标准化较差;
步骤六:在实际数据处理过程中,由于样本较少,或者不稳定测站的干扰,上述基于假设检验的稳定性分析方法并不总是有效,可采用下述方法处理:
2)计算每个测站的稳定性指数;若所有测站稳定性指数均小于1且比较一致,则认为测站均比较稳定,未发生明显变形;若有部分测站稳定性指数大于1且数值出现分群现象,则剔除对应最大稳定性指数的测站及对应基线,针对剩余测站和基线重复迭代执行步骤1)、2),直到剩余测站稳定性指数均小于1且比较一致为止;认为剔除的测站发生形变,不可作为起算基准;剩余测站为稳定测站;
步骤七:分析测量控制网网形,从根据上述步骤得到的稳定测站中选取2-4个能控制整网的测站为起算基准,进行后续平差处理。
实例分析
采用本发明对某平面监测基准网进行起算基准稳定性分析,分析结果如下所示。
表1 某平面监测基准网起算基准稳定性分析结果
表1为某平面监测基准网起算基准稳定性分析情况;由表1可知,在3轮计算过程中,TN03、TN04、TN06三点依次对应每轮最大的值;剔除3点(即TN03、TN04、TN06)后,剩余点位对应的值较小且比较一致;根据上述分析结果,可将此平面监测基准网点分为两群,其中TN03、TN04、TN06为不稳定测站,其余测站为较稳定测站;以此分析结果,根据基准点几何分布,本期此网数据处理以TN01为固定点,以TN01—TN07为固定方向进行全网平差。
验证情况:上述实例本期复测成果与上年度复测成果比较见表2。
表2 某平面监测基准网复测成果较差(2019年-2018年)
其它未说明的部分均属于现有技术。
Claims (2)
1.一种测量控制网基准稳定性分析方法,其特征在于:包括如下步骤,
步骤一:统计基线较差;
统计待检验测站本期与周边其他测站组成的所有基线长度与前期基线长度之差;
步骤二:计算较差标称中误差;
在步骤二中,计算较差标称中误差,具体包括如下步骤,
根据所采用设备标称精度和基线长度,计算各条基线较差的标称中误差;基线较差的标称中误差根据误差传播定律计算:
步骤三:计算基线标准化较差;
采用基线较差归一化处理,获得基线的标准化较差;
在步骤三中,基线较差归一化处理的方法如下,
步骤四:构造检验统计量;
其中,U为检验统计量;为与某一测站相关的基线标准化较差的均值;为待检验测站相关基线条数;为基线标准化较差;为标准差,此处指标准正态分布的标准差,为0,无量纲;为分布的期望,此处指标准正态分布的期望,为1,无量纲;
步骤五:假设检验,判断测站稳定性;
若上述方法检验结果不符合实际情况,则采用下述方法,以进行两种方法之间的互相检核;首先计算测站稳定性指数,根据每条基线标准化较差,计算对应测站所有基线标准化较差的均方根值,称其为此测站的稳定性指数,根据测站稳定性指数迭代判断测站稳定性,确定稳定测站,进入下一步;
步骤六:结合网形,确定平差起算基准;
分析测量控制网网形,从上述步骤得到的稳定测站中选取2-4个测站为起算基准,进行后续平差处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111333451.3A CN113779817B (zh) | 2021-11-11 | 2021-11-11 | 一种测量控制网基准稳定性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111333451.3A CN113779817B (zh) | 2021-11-11 | 2021-11-11 | 一种测量控制网基准稳定性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113779817A CN113779817A (zh) | 2021-12-10 |
CN113779817B true CN113779817B (zh) | 2022-03-11 |
Family
ID=78956940
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111333451.3A Active CN113779817B (zh) | 2021-11-11 | 2021-11-11 | 一种测量控制网基准稳定性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113779817B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103134433A (zh) * | 2013-02-02 | 2013-06-05 | 青岛理工大学 | 一种利用位移监测鉴别边坡失稳致滑因子的方法 |
WO2015044930A1 (en) * | 2013-09-30 | 2015-04-02 | Koninklijke Philips N.V. | Device specific outlier rejection for stable optical shape sensing |
CN106289195A (zh) * | 2016-08-29 | 2017-01-04 | 长江空间信息技术工程有限公司(武汉) | 高原长距离大型工程测量控制***的建立方法 |
WO2017016472A1 (en) * | 2015-07-28 | 2017-02-02 | Huawei Technologies Co., Ltd. | Predicting network performance |
WO2018207524A1 (ja) * | 2017-05-07 | 2018-11-15 | 国立大学法人大阪大学 | 識別方法、分類分析方法、識別装置、分類分析装置および記憶媒体 |
CN109143266A (zh) * | 2018-07-16 | 2019-01-04 | 中铁大桥勘测设计院集团有限公司 | 一种gnss平面控制网已知点兼容性检验的方法 |
CN109520522A (zh) * | 2018-09-27 | 2019-03-26 | 湖北省水利水电规划勘测设计院 | 一种基于三维基线的控制点稳定性判定方法 |
CN109765592A (zh) * | 2019-02-27 | 2019-05-17 | 湖北省水利水电规划勘测设计院 | 一种基于方差协方差矩阵的变形监测网稳定性分析方法 |
CN112099069A (zh) * | 2020-08-31 | 2020-12-18 | 中国三峡建设管理有限公司 | 实测气象参数修正对流层经验模型的rtk算法及应用 |
CN113358017A (zh) * | 2021-06-02 | 2021-09-07 | 同济大学 | 一种多站协同处理的gnss高精度变形监测方法 |
CN113503865A (zh) * | 2021-07-07 | 2021-10-15 | 中国人民解放军战略支援部队信息工程大学 | 一种三维控制网建立方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11269303B2 (en) * | 2009-06-22 | 2022-03-08 | Johnson Controls Technology Company | Systems and methods for detecting changes in energy usage in a building |
US20120086598A1 (en) * | 2010-10-08 | 2012-04-12 | Canadian Space Agency | Apparatus and methods for driftless attitude determination and reliable localization of vehicles |
-
2021
- 2021-11-11 CN CN202111333451.3A patent/CN113779817B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103134433A (zh) * | 2013-02-02 | 2013-06-05 | 青岛理工大学 | 一种利用位移监测鉴别边坡失稳致滑因子的方法 |
WO2015044930A1 (en) * | 2013-09-30 | 2015-04-02 | Koninklijke Philips N.V. | Device specific outlier rejection for stable optical shape sensing |
WO2017016472A1 (en) * | 2015-07-28 | 2017-02-02 | Huawei Technologies Co., Ltd. | Predicting network performance |
CN106289195A (zh) * | 2016-08-29 | 2017-01-04 | 长江空间信息技术工程有限公司(武汉) | 高原长距离大型工程测量控制***的建立方法 |
WO2018207524A1 (ja) * | 2017-05-07 | 2018-11-15 | 国立大学法人大阪大学 | 識別方法、分類分析方法、識別装置、分類分析装置および記憶媒体 |
CN109143266A (zh) * | 2018-07-16 | 2019-01-04 | 中铁大桥勘测设计院集团有限公司 | 一种gnss平面控制网已知点兼容性检验的方法 |
CN109520522A (zh) * | 2018-09-27 | 2019-03-26 | 湖北省水利水电规划勘测设计院 | 一种基于三维基线的控制点稳定性判定方法 |
CN109765592A (zh) * | 2019-02-27 | 2019-05-17 | 湖北省水利水电规划勘测设计院 | 一种基于方差协方差矩阵的变形监测网稳定性分析方法 |
CN112099069A (zh) * | 2020-08-31 | 2020-12-18 | 中国三峡建设管理有限公司 | 实测气象参数修正对流层经验模型的rtk算法及应用 |
CN113358017A (zh) * | 2021-06-02 | 2021-09-07 | 同济大学 | 一种多站协同处理的gnss高精度变形监测方法 |
CN113503865A (zh) * | 2021-07-07 | 2021-10-15 | 中国人民解放军战略支援部队信息工程大学 | 一种三维控制网建立方法 |
Non-Patent Citations (4)
Title |
---|
A study on the model building for adjust baseline estimation using the thermal load calculation;Iwata T等;《支部学术研究发表会前刷集》;20081231;1-10 * |
GPS测量与高精度控制网成果的比较分析;***等;《测绘通报》;19980925(第09期);5-7、14 * |
南水北调中线工程施工测量控制网***研究;杨爱明等;《人民长江》;20101014(第19期);34-37 * |
平面控制网观测误差模拟方法研究;***等;《北京测绘》;20170908;75-78 * |
Also Published As
Publication number | Publication date |
---|---|
CN113779817A (zh) | 2021-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111460567B (zh) | 基于bim楼梯面层净高检查***、应用***的方法及工艺 | |
CN115077425B (zh) | 一种基于结构光三维视觉的产品检测设备及方法 | |
CN111291486A (zh) | 一种数控机床***组件可靠性评估方法 | |
CN106599367A (zh) | 一种航天器状态异常检测方法 | |
CN113779817B (zh) | 一种测量控制网基准稳定性分析方法 | |
CN110850358A (zh) | 基于逐步回归算法的电能表综合检定方法及*** | |
CN112131752B (zh) | 一种基于拟准检定的超强崩溃污染率抗差估计算法 | |
CN107247871B (zh) | 项目检测时间核查预警方法及装置 | |
CN111190200B (zh) | 一种动态环境下单频周跳检测与修复方法 | |
CN110765005B (zh) | 软件可靠性评估方法和装置 | |
CN111898924A (zh) | 一种装备测试性使用综合评价方法及*** | |
CN116339262A (zh) | 一种基于人工智能的数控加工生产质量监测*** | |
CN111060050B (zh) | 一种综合验证三坐标测量机检测软件的方法 | |
CN107798149B (zh) | 一种飞机维修性评估方法 | |
CN108073518A (zh) | 一种测试用例管理方法及装置 | |
JP6394787B2 (ja) | 製品検査装置、製品検査方法及びコンピュータプログラム | |
CN111157938B (zh) | 一种评价自动化检定***计量过程能力的方法和*** | |
CN113567033A (zh) | 一种静态扭矩检验方法 | |
CN106547690A (zh) | 基于新拟合判定标准下的软件可靠性建模方法及装置 | |
CN112329108A (zh) | 一种地铁车站优化抗浮验算方法及*** | |
CN106289492A (zh) | 一种砝码量值分量组合检定法 | |
CN107992287B (zh) | 一种***需求优先级排序结果的检查方法和装置 | |
CN111506967A (zh) | 一种强度试验约束点载荷异常的溯源方法 | |
CN106338770A (zh) | 一种炮检点数据互查方法和*** | |
CN109241114A (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 |