arcgis for jsapi开发：坐标系、经纬度与平面坐标的互换

arcgis for jsapi开发：坐标系、经纬度与平面坐标的互换

练武不练功，到老一场空。学习GIS开发，坐标系是绕不过去的，此乃基本概念，一定要搞清楚。某从事WEB GIS开发差不多2年，一直对这个坐标系感到模糊。这两天因为求解一个问题，方始有点认知。

GIS世界里，坐标系有很多，而且描述得十分抽象，让我们这些没有天赋，缺乏空间想象力的傻叉晕头转向。不过，对于arcgis for jsapi开发人员来说，经常接触的大概只有2种坐标系：

一、4326地理坐标系 通常，我们说地图上的某个点，有经纬度，往往指的就是这种坐标系。

这种坐标系，属于地理坐标系。因为地球是个球体，所以，这是一种球面坐标系。

全称：GCS_WGS_1984，ID为4326，一般GPS、遥感影像、外业测绘数据等就采用它。当然，还有其他地理坐标系，如beijing54，xian80等等，4326只是其中之一。

var sr4326 = new SpatialReference({wkid: 4326 });

二、102113投影坐标系

102113，全称：GCS_WGS_1984_web_mercator_auxiliary_sphere，又称为平面直角坐标系，笛卡尔平面坐标系，Google在线地图所采用。

如上所述，4326地理坐标系是球面坐标系，而我们用的2维地图，是平面的，所以一般用的是平面坐标系。球面坐标系如何换算到平面坐标系，采用投影法。投影法有很多，如米勒投影、墨卡托投影、高斯投影，等等。arcgis中，似乎一般就是墨卡托投影。墨卡托这个名字挺好听，有逼格，我喜欢。102113就是墨卡托投影法搞出来的，故而也可以将它称为墨卡托投影坐标系？

所谓的墨卡托投影是这样的：假设有一个圆锥体或圆柱体，刚好将地球包起来（就是切于地球的赤道），然后地球球心会发光，将地球的经纬线投影到圆锥体或圆柱体上面，这样就得到了地球的平面坐标。圆柱可以看作是特殊的圆锥，即顶点在无限远处。

使用这种方式得到的投影地图在两极会拉长，但经纬度是平行的直线。如图所示：

var sr102113 = new SpatialReference({wkid: 102113 });

因为我们加载的底图一般都用它，也可以直接引用地图对象的：

var sr102113 = map.spatialReference;

如果不确定地图对象的坐标系，可以直接看看地图的地图服务：

百度地图用的是102100（GCS_WGS_1984_web_mercator），应该也是一种墨卡托投影坐标系。

三、坐标系转换 为啥要转换？

当然要转换了。

var pt = new Point(longitude, latitude, new SpatialReference({wkid: 4326 });//最后一个参数，注明要使用何种坐标系来创建本对象

那么，经纬度如何转化为平面坐标系中的点坐标呢？转化出来的点，其单位是什么呢？墨卡托投影，出来的点的坐标，单位是米。所以经纬度换算成平面坐标后，数值都很大：

XMin: 7750103.171888566YMin: 1252344.2714235506XMax: 1.693722247553803E7YMax: 7973910.790707026Spatial Reference: 102113 (3785)

上代码：

var _f = 6378137, _p = 0.017453292519943;function lngLatToXY(pt){//经纬度转化为平面坐标 var lng = pt[0];//经度 var lat = pt[1];//纬度 if (lat > 89.999999){ lat = 89.999999; } else if (lat < -89.999999){ lat = -89.999999; } var c = lat * _p; x = lng * _p * _f; y = _f / 2 * Math.log((1 + Math.sin(c)) / (1 - Math.sin(c))); return [x,y];}

参考资料中有经纬度和平面坐标的互换代码，好像与我给的代码有点不同。我这个代码是同事给的，他是我司GIS的权威和专家。我没仔细比较，先记录下来。

四、实际应用 有辆车，已知当前的坐标、车速和方向，求一段时间后，这辆车预计会到达哪里？

1、这里就要用到坐标转换了。因为经纬度的单位是度分秒，可不是米，假设15分钟后，车子走了10公里，你说它会到哪里？在赤道附近，经纬度一度都约为111公里；但在南北极，经度之间的距离是0！所以要将经纬度先转成平面坐标，再参加运算。

2、本问题是已知圆心（车辆当前位置）、半径（车辆所走的距离）、角度（方向），求圆弧上的点坐标。arcgis地图中，以正北方向为0度，顺时针起算。也就是说，y轴的正方向为0度。有：

圆心坐标(a,b)，半径r的圆的参数方程（根据三角比的定义可以得出）为：(以y轴正方向为0度，设 t 为圆上某一点的度数) x = a + r＊sin t y = b + r＊cos t 所以圆上任一一点的坐标为 P(a + r＊sin t, b + r＊cos t)，（其中 t 为圆上某一点的度数）

function drawSpeed(car){ var t = car; var sp = 转成平面坐标([t.longitude,t.latitude]);//起点 var fp = 获得终点的平面坐标(sp,t.speed,t.direction);//终点 var sr102113 = map.spatialReference; var line ={geometry:{ "paths":[[sp,fp]], "spatialReference":sr102113}, "symbol":{"color":_cbd, "width":1, "type":"esriSLS", "style":"esriSLSSolid"} }; return new Graphic(line);}

幸甚至哉。

参考资料：​​​利用米勒投影法实现经纬度和平面坐标的相互转换​​​​Web墨卡托坐标与WGS84坐标互转​​

版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。