convertCoordinates.py 文件源码

def batchConvertCoordinates(sourceCoordinates, sourceEPSG, targetEPSG):
    targetCoordinates = []
    sourceProjection = osr.SpatialReference()
    sourceProjection.ImportFromEPSG(sourceEPSG)
    destinationProjection = osr.SpatialReference()
    destinationProjection.ImportFromEPSG(targetEPSG)
    coordTrans = osr.CoordinateTransformation(sourceProjection, destinationProjection)
    for lon, lat in sourceCoordinates:
        point = ogr.Geometry(ogr.wkbPoint)
        point.AddPoint(lon,lat)
        point.Transform(coordTrans)
        targetCoordinates.append((point.GetX(), point.GetY()))
    return targetCoordinates