convertCoordinates.py 文件源码

def convertCoordinates(coordinate, sourceEPSG, targetEPSG):
    sourceProjection = osr.SpatialReference()
    sourceProjection.ImportFromEPSG(sourceEPSG)
    destinationProjection = osr.SpatialReference()
    destinationProjection.ImportFromEPSG(targetEPSG)
    coordTrans = osr.CoordinateTransformation(sourceProjection, destinationProjection)
    point = ogr.Geometry(ogr.wkbPoint)
    lon, lat = coordinate
    point.AddPoint(lon,lat)
    point.Transform(coordTrans)
    return point.GetX(), point.GetY()