前言

前段时间数据工程测绘部领导找到我,需要一个坐标转换工具。自2018年7月1号起,国家全面启动CGCS2000坐标系,以前的54坐标系,80坐标系都要转成2000坐标。测绘部曾花5万块钱临时购买了一个转换工具,但后期更高精度坐标转换还需要费用。网上也有一些坐标转换的小工具,但是,只能支持单点或简单txt格式的转换。

GIS专业的我,对于坐标转换原理、思路还是以前在学校的时候学的,于是,我开始网上搜索,布尔莎模型、二维7参数模型、二维四参数模型,高斯正算、高斯反算、最小二乘法。。。结合国家下发的《80T2000坐标转换》方案,以及网上能找到的资料,历时2个月,坐标转换工具终于完成,自己也小有成就感。

工具包括二维四参数转换,和布尔莎七参数转换等。转换结果跟网上多个转换工具比较,结果高度一致,精度远高出国家要求转换精度。由于,要经常验证各种数据精度,基于让我自己使用便捷的原则,哪里不爽改哪里,经过不断完善,参数计算转换使用便捷,验证方便。现将四参数转换过程、实现思路、结果验证,以及遇到的问题分享出来,供大家学习!布尔莎七参数转换,将在后期整理后发出。

界面

先上个二维四参数计算的界面,示例中数据为54转80,平均残差0.04,如果是80转2000,转换精度更高。

该参数计算还支持导入参数验证残差,导入参数0,即可验证不使用转换参数的情况,残差在50m左右。

遇到的问题

1.x y 的困惑

当我第一次把工具交到数据工程部领导时,还比较满意,但是也提了两点,其中一点,他说到“我干了测试这么多年,始终分不清XY,主要是很多工具对XY的定义不一样!”,在他的建议下,我把XY坐标换成了,北坐标和东坐标,这样一来,测绘人员就可以很清楚的做区分了。

2.都是带号的错

刚开始计算的转换参数跟南方CASS工具有一定的差别,纵偏移-500m左右,但是我带入四参数模型发现残差很小,满足转换精度。我把南方CASS计算的参数代入,残差却很大。我用Coord工具验证,计算结果跟我计算的几乎完全一致。

百思不得其解,后来经过反复实验,才发现是带号的问题。南方CASS工具计算的时候去掉了带号,我把带号去掉后,转换参数就一致了。在后面的计算及转换中,我都将带号做了处理。

参数计算

四参数转换模型比较简单,如下:

但是要根据模型,求得四参数的值,就需要将模型变形成下面这种形式:

将该模型用计算机表示,代码如下:

然后利用要用最小二乘法求解,即可。最小二乘法即几次矩阵变换,转置,求逆等,形式比较固定,大家可以网上去查。

坐标转换

计算出四参数后,带入模型就可以进行坐标转换了,该模型适合小区域、大比例尺下坐标转换。ArcGIS并不支持四参数转换,需要将Geometry中坐标带入模型逐点转换。

小结

二维四参数坐标转换的核心是构造出求解矩阵,然后最小二乘法求解。先写这么多了,后期敬请关注!

基于二维四参数模型的坐标转换相关推荐

  1. HALCON示例程序class_2dim_unsup.hdev基于二维直方图的聚类分类

    HALCON示例程序class_2dim_unsup.hdev基于二维直方图聚类分类的颜色分类 示例程序源码(加注释) 读入图片 read_image (Image, 'patras') 关闭窗口 d ...

  2. 基于二维码的设备巡检系统相关技术

    基于二维码的设备巡检的设计与实现 作者:不染心 时间:2022/5/28 项目地址: https://mbd.pub/o/author-aWaVlmpkYw==/work 文章目录 基于二维码的设备巡 ...

  3. 坐标系转换--二维四参数大地坐标系转换模型变换关系理解

    二维四参数空间直角坐标系转换模型变换关系理解 [ X 2 Y 2 ] = [ T x T y ] + ( 1 + m ) [ cos ⁡ α − sin ⁡ α sin ⁡ α cos ⁡ α ] [ ...

  4. AAAI 2020 | 时间可以是二维的吗?基于二维时间图的视频内容片段检测

    作者 | 彭厚文.傅建龙 来源 | 微软研究院AI头条 编者按:当时间从一维走向二维,时序信息处理问题中一种全新的建模思路由此产生.根据这种新思路及其产生的二维时间图概念,微软亚洲研究院提出一种新的解 ...

  5. 时间可以是二维的?基于二维时间图的视频内容片段检测 | AAAI 2020

    作者 | 彭厚文.傅建龙 来源 | 微软研究院AI头条(ID: MSRAsia) 编者按:当时间从一维走向二维,时序信息处理问题中一种全新的建模思路由此产生.根据这种新思路及其产生的二维时间图概念,微 ...

  6. 图像分割——基于二维灰度直方图的阈值处理

    前言 像素灰度值仅仅反映了像素灰度级的幅值大小,并没有反映出像素与邻域的空间相关信息. 二维灰度直方图的概念 二维灰度直方图:像素的灰度值分布和邻域的平均灰度值分布构成的二维直方图 二维直方图的值N( ...

  7. Matlab实现基于二维伽马函数的光照不均匀图像自适应校正算法

    Matlab程序: %% % 基于二维伽马函数的光照不均匀图像自适应校正算法 % %% clc,close all; tic; im=imread('你的图片.jpg'); figure; imsho ...

  8. 二维四象限,就是从非此即彼的二分法里解放出来

    二维四象限,就是从非此即彼的二分法里解放出来,用两个对立统一的重要属性作为依据,画出四象限图,分别讨论情况,逐个解决问题.用二维四象限来分析问题,会让思维更完整.更辩证.风险管理模型.乔哈里视窗理论. ...

  9. 基于二维码的室内定位技术(一)——原理

    作者介绍:周语馨 from 南京大学 to 英特尔亚太研发有限公司 504849766@qq.com 原文地址: https://zhou-yuxin.github.io 哎,不知道怎么说呢.自从朱富 ...

最新文章

  1. QCon讲师对对碰——洪小军采访梁宇鹏:就是爱Golang
  2. 山东大学和哈工大的教师招聘条件对比,心里要有点数
  3. svmtrain和svmpredict简介(转)
  4. 纯虚函数能为private吗?
  5. 可访问范围 与 visual 或 abstract “修饰符”
  6. ajax异步加载和cmd,异步传输Ajax(JQ)
  7. mbot编程机器人怎么连接蓝牙_台式机蓝牙怎么连接
  8. Atitit.100% 多个子元素自适应布局属性
  9. (综述,讲得很好)基于3DMM的三维人脸重建技术总结
  10. Linux一键安装Aria2+Yaaw+FileManager实现BT磁力下载,并在线查看/观看
  11. 常识之外:全表扫描为何产生大量 db file sequential read 单块读?
  12. 微信小程序如何更换头像
  13. xxl-job任务调度平台
  14. fedora安装视频播放器
  15. 我的2011--虚荣、挣扎、总结和转变的一年
  16. 【欧拉计划第 5 题】最小公倍数 Smallest multiple
  17. 视觉平衡与物理平衡_设计中的构图平衡,对称与反对称
  18. Ultra Librarian下载与安装使用教程
  19. 品牌热搜拆解 | 奈雪的茶上市、索尼道歉……
  20. MyBatis动态拼接SQL

热门文章

  1. 电脑丢失dll文件一键修复需要什么软件?快速修复dll文件的方法
  2. Python获取当前时间(time模块)
  3. 网络编程_UDP socket
  4. windows下go安装及govendor设置
  5. 9.19 万兴科技笔试
  6. 可搜索加密机制研究与进展
  7. java多线程下载文件(断点下载、进度展示、网速展示)
  8. Team Leader 你不再只是编码, 来炖一锅石头汤吧
  9. 记录一次nginx配置php-fpm 404
  10. 浏览器F12查看项目源码