CN102779231A - 基于邻近的坐标转换参数计算方法 - Google Patents
基于邻近的坐标转换参数计算方法 Download PDFInfo
- Publication number
- CN102779231A CN102779231A CN201210199071XA CN201210199071A CN102779231A CN 102779231 A CN102779231 A CN 102779231A CN 201210199071X A CN201210199071X A CN 201210199071XA CN 201210199071 A CN201210199071 A CN 201210199071A CN 102779231 A CN102779231 A CN 102779231A
- Authority
- CN
- China
- Prior art keywords
- reference mark
- control point
- coordinate conversion
- adjacent control
- coordinate
- 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
Images
Landscapes
- Processing Or Creating Images (AREA)
Abstract
本发明公开了一种基于邻近的无缝地图坐标转换参数计算方法,属于大地测量、数字地图和地理信息技术领域,该发明主要针对大范围无缝坐标转换参数计算中,待求点周围的控制点如何自动、合理选择,最终实现无缝坐标转换参数合理、高效自动计算。本发明首先基于已知控制点构建Delaunay三角网,然后基于三角网自动提取每个控制点的邻近控制点集,计算参数回代后计算各点的残差并统计中误差,从而对粗差及可疑控制点进行过滤。然后加入图廓点或格网等求参数点,基于控制点构建的Delaunay三角网自动计算各点所对应的邻近控制点集,基于对应的邻近控制点集解决出各点所需的坐标转换参数。访方法可以批量、高效、无缝地计算广大范围内的坐标转换参数,在大地坐标转换、地图数据坐标转换以及地理信息***领域能发挥很大的作用,并获得很好的经济效益。
Description
技术领域
本发明涉及大地测量、数字地图和地理信息技术领域,尤其涉及不同椭球之间坐标转换参数的计算方法。
背景技术
不同国家和地区使用不同的椭球及投影***,由此派生出不同的空间坐标系,进尔导致空间点位的坐标及地图数据坐标系的不同。如我国先后建立了北京54、西安80、2000国家坐标系。由于坐标系的不同,导致测绘成果存在多种坐标形式,如果不通过坐标转换,不同坐标系的测绘成果之间无法直接利用。随着应用的深入,需要高效自动地实现地图数据的无缝坐标转换,地图数据无缝坐标转换的关键是图廓点或格网点坐标转换参数的求解,由于坐标转换参数属于保密数据,有时需要根据控制点数据解算转换参数。对于大范围坐标转换,需要逐图廓点或格网点(以下简称待求点)解算出对应的坐标转换参数(7参数、4参数等),对于大范围或大批量的转换参数计算,如何基于待求点的位置自动合理地选择周围相关的参数计算控制点,是实现大范围无缝转换参数自动计算的关键。
解决问题
该发明主要针对大范围无缝坐标转换参数计算中,待求点周围的控制点如何自动、合理选择,最终实现无缝坐标转换参数合理、高效自动计算。
发明内容
本发明的目的是为了解决大地测量、数字地图、地理信息技术领域中不同测绘成果的自动无缝坐标转换参数计算问题。
本发明所采用的技术方案与步骤如下:
(1)Delaunay三角网生成。基于平面坐标计算所有控制点集构成的Delaunay三角网,剔除最***小于阈值的三角形(如最小角度小于5度),如附图1所示;
(2)邻近控制点集提取。将存在三角形边直接相连的控制点称为邻近控制点,通过控制点集的Delaunay三角网,逐点查找每个控制点对应的邻近控制点集;
(3)控制点检核。参与坐标转换参数计算的控制点可能达不到精度或存在粗差,需要进行过滤与筛选。先逐点查找控制点的邻近点集,然后基于控制点的邻近点集求解出坐标转换参数,回代计算出各个控制点的目标坐标系坐标,求出各点残差,并统计中误差,按照一定阈值筛选可疑控制点。
(4)自动搜索待求点的邻近控制点集。逐三角形判断外接圆是否包含待求点,成立则说明三顶点对应的控制点与待求点存在邻近关系,建立所有待求点所对应的邻近控制点集,将其作为解算待求点坐标转换参数的基础。对于远离控制点的待求点,找出距离该点最近的三个控制点,作为邻近控制点集。
(5)计算待求点的坐标转换参数。基于待求点的邻近控制点集采用4参数或七参数计算模型,解算出待求点处对应的坐标转换参数。
(6)
附图说明
图1是控制点构成Delaunay三角网示意图;
图2是图廓点或格网点邻近控制点集计算查找示意图;
图3是完整计算流程图。
有益效果
通过建立控制点之间的Delaunay三角网,动态***待求点,自动合理地基于拓扑邻近关系查找该待求点周围的相关控制点,计算得到的参数不仅考虑参数计算必要的控制点数,同时邻近点之间自动保持重叠,也将待求点与坐标转换参数建立了一一对应关系。
具体实施方式
下面结合附图和实例,对本发明实施步骤详细说明如下:
如附图1,图中已知控制点共有11个,待求点有140个,计算每个待求点处的坐标转换七参数。
步骤一:基于平面坐标计算所有控制点组成的Delaunay三角网,如图1所示。对于凸壳边上的小角度三角形可以在算法中设定最小角度阈值自动将其过滤掉。图中2-8-10三角形就被过滤掉,特殊的需要可以手动去除或添加。
步骤二:控制点检核。参与参数计算的控制点可能达不到精度或存在粗差,所以需要进行过滤与筛选。首先提取每个控制点的相邻控制点(相邻近控制点之间存在直接Delaunay三角形边相连),然后基于该控制点及邻近控制点列立误差方程,基于最小二乘法解算出转换参数。再基于该参数,通过源控制点坐标计算目标坐标系坐标,求出该控制点的符合差。采用此方法,计算出所有控制点的残差,并统计其中误差,将中误差3倍作为筛选可疑控制点的依据,剔除可疑控制点后再重复步骤二,达到要求为至。如图1中,1号点自动与邻近点3、4、8构成邻近点集,7号点与邻近点2、3、4、6、11、10、8组成邻近点集,并通过同名点的已知源坐标与目标坐标,计算转换参数,回代求出残差,并统计邻近点集的中误差。
步骤三:计算每个图廓点或格网点所对应的邻近控制点集。邻近控制点是指至少存在一个以该控制点为顶点的Delaunay三角形的外接圆包围该图廓点。图2中,A分别位于三角形6-7-11以及7-11-10对应的外接圆内,所以A的邻近控制点集为{6,7,10,11}。对于***没有邻近控制点的格网点,选择距离该点最近的三个控制点,作为邻近控制点,图2中,B点邻近控制点集为{8,7,10}。
步骤四:基于邻近控制点计算所有格网点的坐标转换参数。找出邻近控制点的源坐标与目标坐标,采用四参数或七参数计算模型,列立误差方程,并解算出该格网点的从源坐标到目标坐标的转换参数。
对应的流程图如图3所示。
Claims (3)
1.基于邻近的坐标转换参数计算方法,包括以下几个步骤,其特征在于:
步骤1:基于平面坐标计算所有控制点组成的Delaunay三角网,如图1所示。对于凸壳边上的小角度三角形可以在算法中设定最小角度阈值自动将其过滤掉。
步骤2:控制点检核。参与参数计算的控制点可能达不到精度或存在粗差,所以需要进行过滤与筛选。首先提取每个控制点的相邻控制点(相邻近控制点之间存在直接Delaunay三角形边相连),然后基于该控制点及邻近控制点列立误差方程,基于最小二乘法解算出转换参数。再基于该参数,通过源控制点坐标计算目标坐标系坐标,求出该控制点的符合差。采用此方法,计算出所有控制点的残差,并统计其中误差,将中误差3倍作为筛选可疑控制点的依据,剔除可疑控制点后再重复步骤2,达到要求为止。
步骤3:计算每个图廓点或格网点所对应的邻近控制点集。邻近控制点是指至少存在一个以该控制点为顶点的Delaunay三角形的外接圆包围该图廓点。
步骤4:基于邻近控制点计算所有格网点的坐标转换参数。找出邻近控制点的源坐标与目标坐标,采用四参数或七参数计算模型,列立误差方程,并解算出该格网点的从源坐标到目标坐标的转换参数。
2.如权利要求1所述的基于邻近的坐标转换参数计算方法,其特征在于:基于控制点和拓扑邻近控制点列立误差方程,基于最小二乘法解算出转换参数。再基于该参数,通过源控制点坐标计算目标坐标系坐标,求出该控制点的符合差。计算出所有控制点的残差,并统计其中误差,将中误差3倍作为筛选可疑控制点的依据,剔除可疑控制点后再重复步骤2,达到要求为止。
3.如权利要求1所述的基于邻近的坐标转换参数计算方法,其特征在于:计算每个图廓点或格网点所对应的邻近控制点集。邻近控制点是指至少存在一个以该控制点为顶点的Delaunay三角形的外接圆包围该图廓点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210199071.XA CN102779231B (zh) | 2012-06-18 | 2012-06-18 | 基于邻近的坐标转换参数计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210199071.XA CN102779231B (zh) | 2012-06-18 | 2012-06-18 | 基于邻近的坐标转换参数计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102779231A true CN102779231A (zh) | 2012-11-14 |
CN102779231B CN102779231B (zh) | 2015-11-25 |
Family
ID=47124141
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210199071.XA Expired - Fee Related CN102779231B (zh) | 2012-06-18 | 2012-06-18 | 基于邻近的坐标转换参数计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102779231B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103500530A (zh) * | 2013-04-07 | 2014-01-08 | 郑州大学 | 基于道路网空间分布的自动初始比例尺地图制图方法 |
CN104111061A (zh) * | 2014-06-30 | 2014-10-22 | 中国电力工程顾问集团中南电力设计院 | 一种基础资料缺乏地区测量控制点成果获取方法 |
CN105160197A (zh) * | 2015-09-23 | 2015-12-16 | 湖北省基础地理信息中心 | 一种综合性地理空间数据坐标转换方法及*** |
CN106847069A (zh) * | 2017-03-01 | 2017-06-13 | 王俊锋 | 一种鉴定地图平面数学基础的方法 |
CN108681595A (zh) * | 2018-02-08 | 2018-10-19 | 中国人民解放军61540部队 | 地理pdf地图混合平面坐标***模型及其构建方法 |
CN108871291A (zh) * | 2017-05-12 | 2018-11-23 | 北京中天易观信息技术有限公司 | 一种空中三角测量中大数据匹配点智能压缩方法 |
CN109191579A (zh) * | 2018-07-02 | 2019-01-11 | 广东省地质测绘院 | 坐标转换的方法、装置、计算机设备和存储介质 |
CN109947877A (zh) * | 2018-07-25 | 2019-06-28 | 武汉光庭信息技术股份有限公司 | 一种提高gis移动终端地图定位精度的方法及*** |
CN110542410A (zh) * | 2019-09-03 | 2019-12-06 | 东南大学 | 一种基于格网模型的大地测量成果转换方法 |
CN110941587A (zh) * | 2019-11-18 | 2020-03-31 | 内蒙古师范大学 | 基于附加文件的空间矢量数据存储方法及坐标系转换*** |
CN111723167A (zh) * | 2019-05-10 | 2020-09-29 | 腾讯科技(深圳)有限公司 | 地图坐标转换方法、装置、可读存储介质和计算机设备 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1595067A (zh) * | 2004-06-23 | 2005-03-16 | 武汉虹信通信技术有限责任公司 | 导航***中的一种快速地理坐标投影转换方法 |
-
2012
- 2012-06-18 CN CN201210199071.XA patent/CN102779231B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1595067A (zh) * | 2004-06-23 | 2005-03-16 | 武汉虹信通信技术有限责任公司 | 导航***中的一种快速地理坐标投影转换方法 |
Non-Patent Citations (4)
Title |
---|
WEI ZENG 等: "《Voronoi Diagrams in Science and Engineering (ISVD), 2011 Eighth International Symposium on》", 30 June 2011, article "《Global Surface Remeshing Using Symmetric Delaunay Triangulation in Uniformization Spaces》", pages: 160-169 * |
崔涛 等: "《游动九参数法用于大区域GPS坐标转换的研究与实践》", 《测绘与空间地理信息》, vol. 32, no. 1, 28 February 2009 (2009-02-28) * |
梁洪宝 等: "《基于Delaunay三角网的大范围GPS坐标转换研究》", 《测绘工程》, vol. 20, no. 6, 31 December 2011 (2011-12-31) * |
毛泽卫: "《土地利用信息中空间数据的质量控制》", 《万方学位论文全文数据库》, 30 November 2011 (2011-11-30) * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103500530B (zh) * | 2013-04-07 | 2015-11-25 | 郑州大学 | 基于道路网空间分布的自动初始比例尺地图制图方法 |
CN103500530A (zh) * | 2013-04-07 | 2014-01-08 | 郑州大学 | 基于道路网空间分布的自动初始比例尺地图制图方法 |
CN104111061A (zh) * | 2014-06-30 | 2014-10-22 | 中国电力工程顾问集团中南电力设计院 | 一种基础资料缺乏地区测量控制点成果获取方法 |
CN104111061B (zh) * | 2014-06-30 | 2016-09-21 | 中国电力工程顾问集团中南电力设计院有限公司 | 一种基础资料缺乏地区测量控制点成果获取方法 |
CN105160197A (zh) * | 2015-09-23 | 2015-12-16 | 湖北省基础地理信息中心 | 一种综合性地理空间数据坐标转换方法及*** |
CN106847069A (zh) * | 2017-03-01 | 2017-06-13 | 王俊锋 | 一种鉴定地图平面数学基础的方法 |
CN106847069B (zh) * | 2017-03-01 | 2019-02-26 | 王俊锋 | 一种鉴定地图平面数学基础的方法 |
CN108871291A (zh) * | 2017-05-12 | 2018-11-23 | 北京中天易观信息技术有限公司 | 一种空中三角测量中大数据匹配点智能压缩方法 |
CN108681595B (zh) * | 2018-02-08 | 2021-10-29 | 中国人民解放军61540部队 | 地理pdf地图混合平面坐标***模型及其构建方法 |
CN108681595A (zh) * | 2018-02-08 | 2018-10-19 | 中国人民解放军61540部队 | 地理pdf地图混合平面坐标***模型及其构建方法 |
CN109191579A (zh) * | 2018-07-02 | 2019-01-11 | 广东省地质测绘院 | 坐标转换的方法、装置、计算机设备和存储介质 |
CN109947877A (zh) * | 2018-07-25 | 2019-06-28 | 武汉光庭信息技术股份有限公司 | 一种提高gis移动终端地图定位精度的方法及*** |
CN109947877B (zh) * | 2018-07-25 | 2022-04-29 | 武汉光庭信息技术股份有限公司 | 一种提高gis移动终端地图定位精度的方法及*** |
CN111723167A (zh) * | 2019-05-10 | 2020-09-29 | 腾讯科技(深圳)有限公司 | 地图坐标转换方法、装置、可读存储介质和计算机设备 |
CN110542410A (zh) * | 2019-09-03 | 2019-12-06 | 东南大学 | 一种基于格网模型的大地测量成果转换方法 |
CN110941587A (zh) * | 2019-11-18 | 2020-03-31 | 内蒙古师范大学 | 基于附加文件的空间矢量数据存储方法及坐标系转换*** |
CN110941587B (zh) * | 2019-11-18 | 2023-04-18 | 内蒙古师范大学 | 基于附加文件的空间矢量数据存储方法及坐标系转换*** |
Also Published As
Publication number | Publication date |
---|---|
CN102779231B (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102779231A (zh) | 基于邻近的坐标转换参数计算方法 | |
CN105933294B (zh) | 网络用户定位方法、装置及终端 | |
CN104821013A (zh) | 基于大地坐标系数字高程模型的地表面积提取方法及*** | |
CN104469676A (zh) | 一种移动终端的定位方法及*** | |
CN105427317A (zh) | 一种适用于多视角自动化配准多站地面激光点云数据的方法 | |
CN108241712A (zh) | 一种地图数据处理方法和装置 | |
CN102930509A (zh) | 一种机载激光点云数据的智能化滤波方法 | |
CN109410330A (zh) | 一种基于bim技术无人机航拍建模方法 | |
CN103927788A (zh) | 基于城市竖向规划的建筑地物分层dem模型制作方法 | |
CN103425801A (zh) | 一种面向数字地球的数据融合方法和装置 | |
CN106846251B (zh) | 基于建筑物房顶矢量的正射影像镶嵌线网络自动选择方法 | |
CN105355042A (zh) | 一种基于出租车gps的道路网络提取方法 | |
CN105427177A (zh) | 一种基于gis的农地四至自动计算方法 | |
CN104270785A (zh) | 一种基于地理栅格聚合的无线网络区域问题定位方法 | |
CN102013114B (zh) | 一种基于Microstation v8i的城市快速建模方法 | |
CN113537141B (zh) | 一种堤坝管涌和滑坡病害快速检测方法及*** | |
CN107423373A (zh) | 一种城市级三维建筑模型索引方法 | |
CN106951453A (zh) | 一种快速更新和数据共享的地理实体编码方法 | |
CN106844757A (zh) | 一种天线方位角异常数据挖掘方法 | |
CN105335478B (zh) | 构建城市土地空间立体调查数据语义关联的方法和装置 | |
CN104915539A (zh) | 一种基于遥感分类获取流域下垫面信息的方法 | |
CN105184854B (zh) | 针对地下空间扫描点云成果数据的快速建模方法 | |
CN106875330B (zh) | 一种平面模型转球面模型的方法 | |
CN114881466A (zh) | 一种基于多源数据的人口空间分区拟合方法 | |
CN113032404A (zh) | 一种面向地表覆盖数据的图斑变化轨迹提取方法与*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151125 Termination date: 20180618 |
|
CF01 | Termination of patent right due to non-payment of annual fee |