ArcGIS 墨卡托投影出现偏移

最近项目中接入了其他图商的底图,发现他们的地图是gcj02mc的,坐标中的点位都是米制单位,但是我在地图选点时在百度反向地理编码出现了位置偏移问题,反向地理编码出的位置和原位置隔了大概20km左右。

仔细分析下为什么会出现偏移：

1.情况介绍:

先介绍下我们公司的经纬度情况 (原来的)
1.底图 : 百度墨卡托投影(bd09mc)
2.定位: 取的是系统自带的经纬度(wgs84) IOS
3.百度反向地理编码,百度POI检索: bd09

修改后:
1.底图 : gcj02 Mercator
2.定位: 取的是系统自带的经纬度(wgs84) IOS
3.百度反向地理编码,百度POI检索: bd09

2.问题分析:

使用场景: 地图取点-> 反向地理编码
原来的场景下: bd09mc->bd09
现有场景: gcj02 Mercator->gcj02ll->bd09

完整完成这个流程后发现点位情况发生了偏移

场景	取点	转换
原场景	13276655.077991,3997641.462104	119.264924, 33.947265
现场景	13276655.077991,3997641.462104	119.266222,33.767619

在这能很明显的展示出来左边情况,大概偏移了有20km。
我这边使用的是百度投影转换算法,这边有一个 (在线转换工具)[http://www.site-digger.com/tools/mct2latlng.html]
仔细观察后发现图商使用的是Arcgis中自带的web墨卡托投影算法,我们之前项目中采用的是百度米制算法,仔细了解得知墨卡托投影算法还有很多种。Arcgis中具体的解释:https://desktop.arcgis.com/zh-cn/arcmap/latest/map/projections/mercator.htm ,Google Web 墨卡托,在此之前我已经详细介绍过了我项目里的百度米制算法,这里讲下如何直接使用Arcgis自带方法将web墨卡托坐标转换成地理坐标系。

小tip:

BD09 MC是在标准Web墨卡托的基础上进行GCJ-02加偏之后，再加上百度自身的加偏算法，也就是在Web墨卡托的基础之上进行了两次加偏，该坐标系的坐标值为Web墨卡托格式，单位为米。

3.解决方案:

使用arcgis中的api直接转换坐标

//经纬度转web墨卡托  
+(CLLocationCoordinate2D) LonLatToWebMercator:(CLLocationCoordinate2D) lonLat
{
    // 使用arcgis中的方法直接转换
    AGSPoint* point_2  = [AGSPoint pointWithX:lonLat.longitude y:lonLat.latitude spatialReference:[AGSSpatialReference spatialReferenceWithWKID:4326 WKT:nil]];
    
    AGSSpatialReference *srs = [AGSSpatialReference spatialReferenceWithWKID:3857];
    AGSGeometryEngine *geoEng=[AGSGeometryEngine defaultGeometryEngine];
    
    AGSPoint *mappoint2 = (AGSPoint *)[geoEng projectGeometry:point_2 toSpatialReference:srs];
    
    CLLocationCoordinate2D point;
    point.longitude = mappoint2.x;
    point.latitude = mappoint2.y;
    return point;
}

//web墨卡托转经纬度 
+(CLLocationCoordinate2D) webMercatorTolonLat:(CLLocationCoordinate2D) mercator;
{
     // 使用arcgis中的方法直接转换
    AGSPoint* point_2  = [AGSPoint pointWithX:mercator.longitude y:mercator.latitude spatialReference:[AGSSpatialReference spatialReferenceWithWKID:3857 WKT:nil]];
    
    AGSSpatialReference *srs = [AGSSpatialReference spatialReferenceWithWKID:4326];
    AGSGeometryEngine *geoEng=[AGSGeometryEngine defaultGeometryEngine];
    
    AGSPoint *mappoint2 = (AGSPoint *)[geoEng projectGeometry:point_2 toSpatialReference:srs];
    
    CLLocationCoordinate2D point;
    point.longitude = mappoint2.x;
    point.latitude = mappoint2.y;
    return point;
}

其中 WKID

EPSG:4326
大地坐标系，WGS84
EPSG:3857
投影坐标系，墨卡托投影

也可以设置成102100,从ArcGIS官网得知:102100实际上也是 Web Mercator 的WKID，是Esri早年规定的，现在已经被 3857 取代。

之前也没有系统的学习ArcGIS内容,大部分做的是电商有关的项目,也是到了这家公司之后才开始全面接触ArcGIS,因为之前和后台都是约定的百度米制坐标,也没有关注相关的问题,后面还是要加强学习这方面的书籍,在这里记录下遇到的问题希望和大家一起学习。