? 在生產(chǎn)中，會(huì)遇到將我們既有的坐標(biāo)系轉(zhuǎn)換到高德或者百度坐標(biāo)系下，從而可以使用百度高德的基礎(chǔ)地圖。亦或者將獲得的百度高德數(shù)據(jù)轉(zhuǎn)換到普通的84坐標(biāo)系下。本文基于Arcpy腳本實(shí)現(xiàn)了對(duì)shp或者gdb數(shù)據(jù)的坐標(biāo)系之間的相互轉(zhuǎn)換。

1、從常用坐標(biāo)系到百度高德坐標(biāo)系的轉(zhuǎn)換

思路如下：
1）讀取原始數(shù)據(jù)的坐標(biāo)系；
2）利用Arcpy將原始坐標(biāo)轉(zhuǎn)換成WGS84經(jīng)緯度；
3）對(duì)于線和面，利用Arcpy對(duì)數(shù)據(jù)進(jìn)行修復(fù)；
4）逐個(gè)取點(diǎn)，將WGS84經(jīng)緯度轉(zhuǎn)換成高德或者百度坐標(biāo)系。

# -*- coding:UTF-8 -*-
import arcpy, os
import ConvertLag_Lat as CVLL
import sys
reload(sys)
sys.setdefaultencoding('utf8')

#輸入工作空間
in_features = arcpy.GetParameterAsText(0)
#輸出的空間坐標(biāo)系
outCSStr=arcpy.GetParameterAsText(1)
#輸出工作空間
output_folder = arcpy.GetParameterAsText(2)

outCS=arcpy.SpatialReference(4326)
arcpy.BatchProject_management(in_features, output_folder, outCS)
arcpy.env.workspace = output_folder
transType=0
if(outCSStr=='wgs84'):
  arcpy.AddMessage("transform succeed")
elif(outCSStr=='gaode'):
  arcpy.AddMessage("target coords gaode")
  transType='3'
elif(outCSStr=='baidu'):
  transType='6'
  arcpy.AddMessage("target coords baidu")
arcpy.AddMessage("start transforming......")
featureclasses = arcpy.ListFeatureClasses()
for fc in featureclasses:
  cur= arcpy.UpdateCursor(fc)
  SR= arcpy.Describe(fc).spatialReference
  des = arcpy.Describe(fc)
  arcpy.RepairGeometry_management(fc)
  fc_count=int(arcpy.GetCount_management(fc).getOutput(0))
  arcpy.AddMessage(fc_count)
  arcpy.SetProgressor("step", "transform " fc,0, fc_count, 1)
  if des.shapeType.upper()=='POINT':
    for r in cur:
      geom = r.getValue("SHAPE")
      r.setValue("SHAPE",CVLL.offsetPoint(geom,transType))
      cur.updateRow(r)
      arcpy.SetProgressorPosition()
    del r,cur 
  elif des.shapeType.upper()=='POLYLINE':   
    for r in cur:
      geom = r.getValue("SHAPE")
      r.setValue("SHAPE",CVLL.offsetLine(geom,SR,transType))  
      cur.updateRow(r)
      arcpy.SetProgressorPosition()
    del r,cur 
  elif des.shapeType.upper()=='POLYGON':   
    for r in cur:
      geom = r.getValue("SHAPE")
      r.setValue("SHAPE",CVLL.offsetPolygon(geom,SR,transType))
      cur.updateRow(r)
      arcpy.SetProgressorPosition()
    del r,cur
  arcpy.AddMessage(fc " transform finished")

2、百度高德與WGS84互轉(zhuǎn)

import arcpy, os
import ConvertLag_Lat as CVLL
import sys
reload(sys)
sys.setdefaultencoding('utf8')

#input path of folder where all shapefiles reside in
in_fold = arcpy.GetParameterAsText(0)
output_folder = arcpy.GetParameterAsText(1)
outCSStr= arcpy.GetParameterAsText(2)
outCS=arcpy.SpatialReference(4326)

arcpy.env.workspace = in_fold
#copy data
featureclasses = arcpy.ListFeatureClasses()
for fc in featureclasses:
  arcpy.Copy_management(fc,output_folder "\\" fc)
arcpy.env.workspace=output_folder

transType='0'
if(outCSStr=='mars2baidu'):
  transType='1'

elif(outCSStr=='baidu2mars'):
  transType='2'

elif(outCSStr=='mars2baidu'):
  transType='3'

elif(outCSStr=='mars2wgs84'):
  transType='4'

elif(outCSStr=='baidu2wgs'):
  transType='5'

elif(outCSStr=='wgs2baidu'):
  transType='6'


arcpy.AddMessage("transfor type is " str(transType))
featureclasses = arcpy.ListFeatureClasses()
for fc in featureclasses:
  cur= arcpy.UpdateCursor(fc)
  SR= arcpy.SpatialReference(4326)
  des = arcpy.Describe(fc)
  arcpy.RepairGeometry_management(fc)
  fc_count=int(arcpy.GetCount_management(fc).getOutput(0))
  arcpy.SetProgressor("step", "transform " fc,0, fc_count, 1)
  if des.shapeType.upper()=='POINT':
    for r in cur:
      geom = r.getValue("SHAPE")
      r.setValue("SHAPE",CVLL.offsetPoint(geom,transType))
      cur.updateRow(r)
      arcpy.SetProgressorPosition()
    del r,cur 
  elif des.shapeType.upper()=='POLYLINE':   
    for r in cur:
      geom = r.getValue("SHAPE")
      r.setValue("SHAPE",CVLL.offsetLine(geom,SR,transType))  
      cur.updateRow(r)
      arcpy.SetProgressorPosition()
    del r,cur 
  elif des.shapeType.upper()=='POLYGON':   
    for r in cur:
      geom = r.getValue("SHAPE")
      r.setValue("SHAPE",CVLL.offsetPolygon(geom,SR,transType))
      cur.updateRow(r)
      arcpy.SetProgressorPosition()
    del r,cur
  arcpy.AddMessage(fc " transform finished")

不足

這里有幾個(gè)方面的不足
1）我使用的是git上給出的坐標(biāo)轉(zhuǎn)換方式，經(jīng)測(cè)試轉(zhuǎn)換后有米級(jí)的誤差；
2）這里用的是取點(diǎn)在重組的方式，因此對(duì)于弧，多部件要素的支持不行，只能支持對(duì)簡(jiǎn)單的點(diǎn)線面進(jìn)行轉(zhuǎn)換。
這里給出源碼，希望有需要的人可以繼續(xù)優(yōu)化完善。