大同推广型网站开发,黄骅市网站建设价格,代写文章价格表,沈阳图书设计公司本文转载于CSDN作者rsyaoxin这是一篇关于坐标和投影的「神文」截止目前浏览量已达2698239是相关文章中不可打破的神话...文末有本文作者推荐的两款坐标转换的小工具下载链接回想一下#xff0c;接触遥感专业也有几个年头了#xff0c;而现在越来越偏离遥感了#xff0c;突然… 本文转载于CSDN作者rsyaoxin这是一篇关于坐标和投影的「神文」截止目前浏览量已达2698239是相关文章中不可打破的神话...文末有本文作者推荐的两款坐标转换的小工具下载链接回想一下接触遥感专业也有几个年头了而现在越来越偏离遥感了突然想着把自己脑中的遥感知识整理出来。首先想到的便是坐标系和投影我想这个东西困扰着80%以上的测绘、遥感和GIS领域的从业人员吧经常有人问我自己曾经也很迷糊什么大地坐标系啊、地心坐标系啊、高斯投影啊、UTM投影啊搞得头都大了可是这玩意又不能不玩毕竟空间信息是遥感数据的灵魂如果定位不准影像就只能当风景画看了。好吧那就啃吧基础的东西还是要的经过一番探索再加上去年专门做过坐标转换一些工作算是清晰明了些了在此就把我的一些见解分享给大家欢迎各位菜鸟和专家前来吐槽。好了咱们言归正传开始坐标系探险记。1、什么是坐标系 或许很多人心中明白这个概念却又不甚明白。什么是坐标系人们描述空间中的一个点或者一个位置通常会采用坐标这个概念。可是这个坐标该怎么计算呢它的参考是谁呢如果中国定义中南海坐标是(1921,1949,2012)美国佬也定义他们白宫的坐标是(1921,1949,2012)哪天黑马说咱们打它中南海一炮吧那就可能把他自己烤熟。所以必须要建立一个合理的坐标系(⊙﹏⊙b汗~~这个理由好牵强啊)。 坐标系由原点和坐标轴组成。坐标系种类很多我们大家在数学中想必都学过笛卡尔坐标系、极坐标系、球面坐标系和柱面坐标系吧在地学领域用到最多的是平面坐标系、空间直角坐标系(前两者属于笛卡尔坐标系)和球面坐标系。比如说一个点坐标是(-2850017.4724690744.5233237959.973)就是指空间直角坐标而我们经常看到的Google Earth上的点的坐标(37°20′17″N112°33′20″E)就是指球面坐标。 测量学上坐标系怎样定义的呢地球是一个不规则的类椭球怎样用严格的数学方式表示它应该是测绘学家们所追求的高峰。为了表示地球上每一个点的位置是不是要建立一个统一的世界坐标系呢建立坐标系是不是要确定坐标原点和坐标轴呢那怎样建立呢测量学家们把地球当作一个规则的椭球来处理这下就好办多了椭球中心就原点呗长轴短轴就作为坐标轴这样一个坐标系就出来了啊。可是有人会问那这个椭球怎么表示呢嘛嘛的我也想问原来这是科学家们利用天文观测得到的而且不同的科学家得到的椭球还不一样比如一个叫Krasovsky的人搞了个克拉索夫斯基椭球还有人搞了什么IUGG-1975、WGS-84、GRS80椭球他们主要特点是长半轴和扁率不同。可是又有人会问(Y的问题好多啊)搞这么多椭球干嘛有一个不就行啦。呵呵我们知道地球坑坑洼洼的用严密的椭球来表示肯定有误差有的国家为了使自己的国家与椭球面吻合(最好大家都站在椭球面上)这样根据各自的情况就定义了不同的参考椭球比如北京54坐标系就采用了苏联老大哥的克拉索夫斯基椭球。可是问题又来了怎样才算吻合得好呢肯定会有人站在椭球面上有人站在椭球面下真头疼。此时测量学家们引入了大地基准面来衡量椭球与大地的吻合度。大地基准面是由大地水准面而来是指平均海平面延伸到大陆得到的一个封闭曲面。比如在建立北京54坐标系时专家们肯定会选择与中国的大地水准面吻合比较好的椭球。此时的椭球称为参考椭球建立的坐标系称为参心坐标系我国的北京54和西安80坐标系都是参心坐标系是一种局部范围的坐标系。然而这种坐标系对于全球定位来说极其不便误差很大所以山姆大叔率先针对GPS系统设计了全球大地坐标系WGS-84坐标系统这时的大地原点不再是参考椭球的中心而是地球的质心。WGS84椭球体的相关参数和WGS84坐标系的坐标轴指向请参考相关专业书籍。我国现有的国家2000坐标系也是一种全球大地坐标系其与WGS-84坐标系稍微有点差异。下面是几种常见坐标系的椭球参数其中北京54坐标系和西安80坐标系是参心坐标系而WGS-84坐标系与国家2000坐标系是地心坐标系坐标原点是地球质心。 好啦这样大家明白了坐标系的定义了吧首先需要定义参考椭球体有了参考椭球还需要大地基准面(全球大地坐标系就不要了)然后需要定义坐标系原点和坐标轴的指向。这样一个坐标系就建立了以后找妹子就方便多了全球定位吧哈哈2、为什么要投影 大家会想有了坐标就行了为什么还要搞个让人迷糊的投影呵呵前面我们讲到的是以椭球体为参考来进行空间定位一点都不直观如果哪天你和妹纸约会妹纸说她在(-2850017.4724690744.5233237959.973)或者(112°E38°N)尼玛坑爹这到底在哪个国家离哥哥我有多远啊不知道啊不至于拿个尺子去测吧。这时候泡妞高手们想出了一个办法把球面投影到一个平面用一个平面坐标(x,y)来表示地面点的位置两点之间求距离是不是很容易啊这时候你会发现那个妹子不就是隔壁那妞嘛200米不到(囧哈哈)。当然投影最大的目的不是方便把妹纸而是地图。所以投影就是把球面坐标转化为平面坐标也就是3D到2D的转换。 投影有很多种按性质分比如等角投影等积投影等距投影任意投影等。大家都知道球面展开成平面肯定是一个不严密(也可说不完美)的过程会有不同程度的变形。如何选择呢比如在航海上就需要等角投影如果方向错了就会差很多我猜如果哥伦布那时知道这些就不会跑到美洲还以为到了印度吧。如果需要丈量面积那就要选择等积投影了。3、测量坐标有哪几种表示方式 常见的测量坐标包括大地坐标(B、L、H)、空间直角坐标(X、Y、Z)、平面坐标(x、y、H)。具体参考相关教程。4、坐标转换 终于讲到重点了各种坐标怎样转换是大家最关心的。首先有一点要牢记 同一参考椭球下大地坐标与空间直角坐标之间的转换是严密的(数学关系对应)它们与平面坐标的转换是不严密的需要做投影转换(想想也明白把球面展成平面那可是难住了好多科学家呀)。而不同参考椭球之间的坐标转换永远都是非严密的。 坐标转换原理同一椭球下的转换 同一椭球下大地坐标(B、L、H)与空间直角坐标(X、Y、Z)之间的转换是严密的其公式为 而大地坐标(B、L、H)与空间直角坐标(X、Y、Z)向平面直角坐标的转换属于非严密的需要进行球面到平面的投影选择通常将空间直角坐标转换为大地坐标然后在大地坐标和平面直角坐标之间采用高斯正算和反算公式进行计算。不同椭球下的转换 不同参考椭球下的坐标转换实质是基准的转换。如空间定位技术所采用的全球基准与地面网所采用的局部基准间的转换。通常的转换模型有布尔莎-沃尔夫模型和莫洛金斯基模型。这两种模型都常用且非常相似布尔莎模型在进行全球或者较大范围内较为常用但是莫洛金斯基模型可以克服布尔莎模型中旋转参数与平移参数相关性高的问题。 两个坐标系的转换通常有三维七参数模型和二维四参数模型。 布尔莎模型又称为七参数转换或者七参数赫尔默特变换。该模型共采用7个参数分别为三个平移参数(ΔX、ΔY、ΔZ)和三个旋转参数(ωx、 ωy、ωz)和一个尺度参数k。 上式是一个WGS84下的空间直角坐标转换到CGCS2000下的空间直角坐标的布尔莎模型有七个未知参数简单的求解只需要3个公共点就可以了如果要得到严密解就需要更多的公共点进行最小二乘平差解算。而对于大地坐标可以转成空间直角坐标再解算也可以直接利用布尔莎模型。5、遥感制图的投影和坐标系选择 关于遥感影像的坐标系和投影大家经常会听到所谓的地理坐标系(Geographic coordinate system)和投影坐标系(Projection coordinate system)的概念。打开Arcgis的坐标系文件夹也会发现这两个子文件夹。地理坐标系是以经纬度表示的是一种球面坐标系而投影坐标系是米或者千米等为单位的是一种以xy表示的平面坐标系。下面是Arcgis中两种坐标系的定义(a) 地理坐标系 (b)投影坐标系 上图中的左图a是WGS-84下的地理坐标系我们可以看出它定义了一个WGS-84椭球体然后是WGS84下的基准面。有了Spheroid和Datum就可以使用地理坐标系了。右图b是西安80坐标系下的投影坐标系投影方式是高斯-克吕格投影。可以看出它除了Gauss-Kruger投影参数外还定义了Xian-1980地理坐标系。 有人会问为什么投影坐标系要带一个地理坐标系呢我们知道投影本质就是把球面投影到平面那么投影就必须定义是用哪个球面坐标系来投影。比如UTM投影坐标系我们熟知的是WGS84坐标系来投影其实我们也可以用国家2000坐标系来投影就我所知在月球探索过程中人们也用到了UTM投影而这时用来投影的坐标系就是月球坐标系了。我国常用地图投影的选择 在我国大中比例尺(≥1:50万)的地图一般选择高斯-克吕格投影高斯-克吕格投影有3°和6°分带两种1:2.5万1:50万比例尺地形图采用经差6度分带1:1万比例尺的地形图采用经差3度分带。小比例尺的地图(比如全中国地图)应采用等角正轴割圆锥投影也就是兰勃特投影(Lambert Conformal Conic)但是小比例尺海图多选择等角正轴圆柱投影也就是墨卡托投影。另外我们平时下载的遥感影像大都是WGS-84下的地理坐标系或者UTM投影坐标系。 下面就一个例子具体说下地理坐标系到投影坐标系的转换(投影的选择) 比如北京的中心经纬度大致在(117°E40°N)那么其在北京54投影坐标系下6°分带就应该选择Beijing 1954 GK Zone 20.prj或者Beijing 1954 GK Zone 20N.prj(后者会在横坐标前加上带号)中央经线带号(比如20)*6-3。而在西安80投影坐标系下6°分带应该选择Xian 1980 GK Zone 20.prj。北京地区的遥感影像的UTM投影一般会选择WGS 1984 UTM Zone 50N因为Zone为20*6-3的中央经线为117左114右120正好覆盖北京而东半球要加30(全球360/6共60个Zone)所以是50NN表示北半球。此外需要说明的是跨带投影会造成较大的变形。很多人下载全球的植被图、夜间灯光图等想要裁剪出自己需要的区域原始数据往往是WGS84下的地理坐标系而矢量数据多是UTM或者TM等投影坐标系这时需要将矢量数据转换到地理坐标系再去裁剪如果反过来将全球图去投影会发现变形得很难看。6、推荐一些好用的坐标转换工具(1) 坐标转换工具用过很多坐标转换工具也自己写过现在给大家推荐一款GPS工具箱发现它很容易用精度也不错这是下载地址在文末。可能大家平时会用的比较多的是Coord这个软件但是我发现它结果老是不对下面是几个软件的比较第一个是GPS工具箱第二个是我自己写的一个小软件第三个是Coord结果明显不对不知道是不是我下载的Coord有问题。(2)开源坐标和投影转换库Proj4。它的功能主要有经纬度坐标与地理坐标的转换坐标系的转换包括基准变换等。这是官方网站http://trac.osgeo.org/proj/可以下载到源码和文档源码是C所写需要编译才能集成到你的程序使用。公众号对话框发送613获取两款坐标转换工具的下载链接- END -现有测绘成果转换到2000国家大地坐标系统指南2000国家大地坐标系转换技术培训(含DLG,DOM,DEM,DRG转换方法及后处理)一篇非常详细的坐标系统学习PPT值得学习坐标问题不再头疼墨卡托、高斯-克吕格、UTM...这些投影的小知识你真的知道吗
