arcgis api for silverlight开发系列之三：定义坐标系

让我们用一个例子开始

<esri:Map x:Name="MyMap" Background="White">
<esri:Map.Extent>
<esri:Envelope XMin="661140" YMin="-1420246" XMax="3015668" YMax="1594451" >
<esri:Envelope.SpatialReference>
<esri:SpatialReference WKID="26777"/>
</esri:Envelope.SpatialReference>
</esri:Envelope>
</esri:Map.Extent>
<esri:Map.Layers>
<esri:ArcGISDynamicMapServiceLayer ID="DynamicLayer1"
Url="http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Demographics/ESRI_Census_USA/MapServer"/>
<esri:ArcGISDynamicMapServiceLayer ID="DynamicLayer2"
Url="http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Petroleum/KGS_OilGasFields_Kansas/MapServer"/>
</esri:Map.Layers>
</esri:Map>

在这个例子中，首先我们定义了地图展现的区域，接着定义了坐标系。26777是坐标系的ID，其对应的坐标系是：

26777

NAD_1927_StatePlane_Kansas_North_FIPS_1501

属于投影坐标系，更多的投影坐标系你可以在这里查看：http://resources.esri.com/help/9.3/arcgisserver/apis/javascript/arcgis/help/jshelp/pcs.htm

更多地地理坐标系你可以在这里查看：http://resources.esri.com/help/9.3/arcgisserver/apis/javascript/arcgis/help/jshelp/gcs.htm.

关于地理坐标系和大地坐标系的详细概念网上很多写的都很生涩，其实大家做开发不必专门去了解很多地图学的知识。我们只要理解，凡是以经纬度为单位的都是地理坐标系，因为它归根结底是一个椭球体，只不过各个国家为了反映该国家所在区域地球的真实形状，而采用不同的数学模型对本不是椭球体的地球进行椭球体化。而投影坐标系，是对地理坐标系按照某种方式投影到平面上的，所以可以认为它是一个平面坐标系，单位自然是米或千米。

地理坐标系

投影坐标系

那么我们在做开发的时候，尤其是web地图开发，两种坐标系至关重要

4326 GCS_WGS_1984

102100 WGS_1984_Web_Mercator_Auxiliary_Sphere

前者是wgs1984，属于地理坐标系，相信大家对它都有所耳闻，他就是大名鼎鼎的GPS采用的坐标系，也就是通过GPS拿到的坐标信息都是按这个坐标系给我们的经度和纬度。当然，如果你是做移动平台上的GPS，获得的经纬度也是按这个坐标系咯。
后则是目前在线地图采用的通用坐标系，属于投影坐标系，建议各位朋友看看我之前的一篇博客WGS 1984 Web-Mercator 对于在线地图服务的意义。相信你会理解的更多。

了解了这么多，大家应该明白一件事，如果你们采用GoogleMap做地图，然后想通过GPS将位置在地图上显示，不经过任何转换直接在GoogleMap上显示是不行的，因为他们的坐标系不统一。所以在显示之前就必须将GPS获取点进行坐标转换到WGS_1984_Web_Mercator，然后在GoogleMap上显示。arcgis api for silverlight中提供了两种坐标系的转换方法，首先要引入ESRI.ArcGIS.Client.Bing的dll，然后调用ESRI.ArcGIS.Client.Bing.Transform.GeographicToWebMercator方法就搞定。如果，你要自己写方法将两坐标系转可以借鉴下面的方法：

//经纬度转墨卡托
public MapPoint lonLat2Mercator(MapPoint lonLat)
{
MapPoint mercator = new MapPoint();
double x = lonLat.X * 20037508.34 / 180;
double y = Math.Log(Math.Tan((90 + lonLat.Y) * Math.PI / 360)) / (Math.PI / 180);
y = y * 20037508.34 / 180;
mercator.X = x;
mercator.Y = y;
return mercator;
}

//墨卡托转经纬度
public MapPoint Mercator2lonLat(MapPoint mercator)
{
MapPoint lonLat = new MapPoint();
double x = mercator.X / 20037508.34 * 180;
double y = mercator.Y / 20037508.34 * 180;
y = 180 / Math.PI * (2 * Math.Atan(Math.Exp(y * Math.PI / 180)) - Math.PI / 2);
lonLat.X = x;
lonLat.Y = y;
return lonLat;
}