中添加2000坐标系_坐标系那些事儿(下)

3296124878666cb1409c39fc44873340.png

3、带号与中央子午线的对应关系

分了那么多投影带,为了方便区分我们得给它们编个号,也就是带号。怎么编带号呢?下面我就开始本文最难的计算部分了。 从0度经线(通过英国伦敦格林尼治天文台原址的那条经线称为0度经线,也叫本初子午线)开始算,每隔6度为一带,依次编号为1、2、3……。在实际工作中,我们需要非常明确每一带的带号和它的中央子午线的换算方法。其实根本不用算,我们就从第1个投影带开始数,第1带的中央子午线是3度经线,第2带的中央子午线是9度经线,第3带的中央子午线是15度经线,是不是已经能列出公式了,中央子午线=带号✖️6-3。把这个公式换一下形式就能得到:带号=(中央子午线+3)÷ 6。是不是很简单? 我们再看3度带。3度带是在6度带的基础上形成的。我们还从第1带开始算,3度带第1带的中央子午线和6度带第1带的是重合的,也就是说3度带第1带并不是从0度经线开始算起的,而是从1.5度经线开始的,那么3度带第1带的中央子午线就是3度经线,第2带的中央子午线就是6度经线,第3带的中央子午线就是9度经线,是不是已经看出规律了?3度带的中央子午线和带号的换算方法就是:中央子午线=带号✖️3, 带号=中央子午线÷ 3。是不是6度带还简单?

d1af12597bf495dfcc213f3c980de098.png

我国领土大概覆盖到73度到135度经线,按换算方法计算一下,就是6度带的第13带到第23带,3度带的第25带到第45带。所以只看带号就能判断投影是6度带还是3度带。 4、高斯平面直角坐标系 和数学坐标系一样,投影完的高斯平面直角坐标系也需要有横轴、纵轴和原点,我们把每个带的中央子午线的投影作为纵轴,把地球赤道的投影作为横轴,两轴交点作为原点。需要注意的是,测绘中的高斯平面直角坐标系和数学直角坐标系的默认的x轴和y轴是反着的,我们把横轴叫做y轴,把纵轴叫做x轴,上北下南左西右东,所以我们把x坐标叫做北坐标,把y坐标叫做东坐标。

48a946c11341cae544a01904faf5d861.png

其次,有一个向东平移500公里的事儿需要解释一下。我国位于北半球,因此x坐标都是正数。但在中央子午线左边(西边)的y坐标就是负数了,负数在平时的工作中很不方便。粗略的算一下,在赤道上1秒角度的距离大概是30米,那么3度(6度带中央子午线左边的一半)的距离就是30✖️3600✖️3 = 324000,也就是最小的y坐标是-324000。为了避免出现负数,同时让y坐标都是6位数,我们就统一给所有位置的y坐标都加上500000,也就是把坐标向东平移500公里。 最后,为了从坐标上看出点位位于哪个投影带,需要在东坐标的前面加上2位带号。比方说第38度带的y坐标是512345.123,加上带号就是38512345.123,相当于把坐标向东平移了38000000。 有了上面对高斯平面直角坐标系的了解,看到一个点的坐标,就能对它的位置有大概的了解。举个例子,一个点坐标:x=4045678.456,y=38543210.123,x坐标是7位数的正数,是在北半球,粗略的计算,纬度应该在37.5度左右,y坐标是8位数,说明前两位是带号,第38度带,那肯定是3度带投影,那中央子午线就是38✖️3 = 114,剩下的543210.123,说明这个点的位置是在中央子午线114度的东边,因为它大于500公里。

四、坐标系转换

坐标系转换这个问题现在比较火,我们单位从事规划的同事也经常会找我们帮着把已有的坐标或地形图转换到2000国家大地坐标系下。为什么现在比以前多了呢?因为从2018年7月1日起全面使用2000国家大地坐标系,这也意味着1954北京坐标系和1980西安坐标系正式退出历史舞台。那手头儿现有的那些北京54和西安80的成果也不能全扔掉从头再来啊,那成本就太高了。我们还有方法把它们转换到2000国家大地坐标系下。最稳妥的方法就是计算转换参数,找至少3个控制点,这些控制点能大致将你要转换坐标系的区域包含住,它们同时具有两个坐标系下精度可靠的坐标,然后通过专业软件拟合计算出转换参数,然后不是马上用参数进行坐标系转换,而是检查转换参数是否精确。怎么检查呢?这里就涉及到误差和限差了。 比方说,你要把一些测绘成果由北京54坐标系转到2000坐标系,有10个控制点同时有这两个坐标系下的坐标,计算转换参数时用了1~5号点5个控制点,6~10号控制点当做检查点,不参与转换参数的计算。参数计算出来了,下面要检查参数是否精确了。先用这套参数把1~5号控制点的北京54坐标转换成2000坐标,然后和他们原本就有的2000坐标相减统计误差,这叫内符合精度;然后用这套参数把6~10号控制点的北京54坐标转换成2000坐标,再和他们原本就有的2000坐标相减统计误差,这叫外符合精度。这两种精度数值就是这套参数转换坐标系的误差的体现,如果这个误差满足规范中的限差要求,就可以继续下一步坐标转换工作了,否则就得查找原因,重新计算,直到满足限差要求。 话又说回来,这些控制点的坐标成果可不是那么容易就能得到的,而且转换参数也是属于保密数据。如果需要准确可靠的2000国家大地坐标系的测绘资料,那还是交给我们测绘专业的人来做吧,毕竟这确实是个瓷器活儿。当然,如果精度要求不高的话,在有必要的基准数据的情况下,也可以在GIS软件中做一下配准,具体的配准方法可以点开下面的链接,有具体的操作步骤。 ArcGIS 地理配准步骤(底图校正)详细讲解 五、注意事项 了解了以上这些知识点,应付跟坐标系相关的日常工作应该是够用了。除此之外,还有几点需要特别注意的。 1、高斯投影与坐标系转换的区别 高斯投影与坐标系转换最大的区别就在于转换后有没有误差。前者只是根据严密的高斯投影正反算公式进行计算,把经纬度坐标转成高斯平面直角坐标叫高斯正算,反之,把高斯平面直角坐标转成经纬度坐标叫高斯反算,只是同一坐标下的坐标的两种形式的互换而已,就像一个人的大名和小名的关系一样,不管称呼范闲叫范闲也好、叫他小闲闲也罢,范闲还是那个范闲,互换不需要转换参数,不存在误差。每一种坐标系都可以有经纬度和高斯平面直角坐标,就像很多人都有大名和小名一样。而坐标系转换就不同,它是指两种坐标系之间的转换,比方说要从北京54坐标系转到2000国家大地坐标系,转换过程中必然会引入误差,只要满足限差要求,就可以使用。我举个例子,也许不是很恰当,这就好比是把低版本word的doc文件转换成高本版word的docx文件,虽然看起来没什么大的变化,但其实里面多少会有些小变化,比方说之前页眉页脚是居中的,转换后就不居中了,公式没法再编辑了等等,但那都不是大毛病,可以忍的话就可以继续用,不然就再简单编辑一下就可以了。 2、GIS软件需要prj文件 现在好多同事都在用使用GIS软件进行国土空间规划工作,有时候会遇到两套数据没法叠加到一起,或者软件提示缺失投影信息,那是因为prj文件出了问题。GIS软件跟CAD不一样,CAD只是固定的显示平面直角坐标,而且在CAD里不需要体现坐标系信息,而GIS软件一般是可以动态的显示坐标信息的,只要不涉及坐标系转换,如果想在经纬度坐标和高斯平面直角坐标系之间进行切换,在GIS软件中设置一下即可,要实现这个功能就需要每个加载到GIS软件中的数据具有准确的投影信息,一般都是存储在文件名是和数据文件同名、扩展名是prj的文件里,就是投影的英文projection的前三个字母,其中包括椭球参数、投影的中央子午线、北平移量、东平移量等等信息,有个这些信息GIS软件就可以动态的显示数据,什么叫动态呢?比方说,你有两个shp格式的数据,都是2000国家大地坐标系的,一个是某市的行政区划的范围线,中央子午线是114,坐标有带号,另一个是这个市内的交通主干道,中央子午线是115,坐标没带号。这两个数据完全不是一个投影带的,一个是标准高斯投影,一个是独立坐标系,但是只要prj信息完整,就能在GIS软件中套合到一起,准确的显示两者的相对位置关系。所以,prj文件很重要,拷数据的时候一定要一并拷走,文件虽小,重要性巨大。 3、应用测绘成果前务必确定坐标系 同一个点的不同坐标系下的坐标值会差很多,少的差几米,多的甚至差上百米,所以在拿到测绘成果时首先要做的是确定它是什么坐标系的,是不是规划成果最终要的坐标系,如果不是可根据具体工作需要和已有数据情况选择严密的坐标系转换或者粗略配准。否则,后果不堪设想! 4、尽量不要手动编辑测绘成果 有时候会把从多个途径获取的数据放在一起查看,但有的数据又没有准确的坐标信息,这种情况下尽量不要手动平移、旋转和缩放测绘成果,因为它的坐标信息是准确的,是我们费了九牛二虎之力才绘制成的,你这么简单的一挪基本上就浪费了我们一半的劳动成果,再做个缩放旋转,另一半也没了。实在逼不得已要改,那就先做个备份,之后发现错误了还有回头路。 好了,关于坐标系的一些事情就跟大家分享完了,希望能对大家的日常工作有一点点帮助。如果你觉得这篇文章确实能帮到你和你的同事,请分享、收藏、点“在看”。这篇文章末尾应该是能留言了,有问题可以给我留言,我一定抽时间答复。 多谢关注,下次再见!
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值