2000坐标系转换为wgs84坐标系

记录一下如何进行2000坐标系转化为wgs84坐标系的过程。

在处理矢量数据转换的过程中，有几个关键步骤确保了数据的有效性和可用性。

首先，在转换之前，通过检查并创建目标保存目录，保证了正确的目录结构，这对于避免转换失败或产生错误输出至关重要。

接着，通过osr.SpatialReference对象定义目标坐标系，并使用ImportFromEPSG()方法设置EPSG代码，确保了转换后的数据能够在地理信息系统中正确显示。

此外，通过selectFields参数可以选择性地保留某些字段，而geometryType参数则定义了输出文件中的几何类型，这些选项使用户可以根据具体需求定制输出文件。

重投影选项reproject=True指示GDAL/OGR进行坐标系转换，这对于整合不同坐标系下的数据至关重要。

最后，通过dim参数可以控制输出数据的维度，例如仅保留X和Y坐标，这对于某些特定应用场景是必要的。这些步骤共同构成了一个完整的矢量数据转换流程，确保了数据的准确性和实用性。

具体代码如下：

from osgeo import ogr, gdal
from osgeo import osr
import os

def VectorTranslate(
        shapeFilePath,
        saveFolderPath,
        format="GeoJSON",
        accessMode=None,
        dstSrsESPG=4326,
        selectFields=None,
        geometryType="POLYGON",
        dim="XY",
):
    """
    转换矢量文件，包括坐标系，名称，格式，字段，类型，纬度等。
    :param shapeFilePath: 要转换的矢量文件
    :param saveFolderPath: 生成矢量文件保存目录
    :param format: 矢量文件格式
    :param accessMode:None代表creation,'update','append','overwrite'
    :param dstSrsESPG: 目标坐标系EPSG代码，4326是wgs84地理坐标系
    :param selectFields: 需要保留的字段列表如果都保留，则为None
    :param geometryType: 几何类型,"POLYGON","POINT"。。。
    :param dim: 新矢量文件坐标纬度,建议查阅官方API。
    :return:
    """
    if os.path.exists(saveFolderPath) ==0:
        os.makedirs(saveFolderPath)
    ogr.RegisterAll()
    gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "YES")
    data = ogr.Open(shapeFilePath)
    layer = data.GetLayer()
    spatial = layer.GetSpatialRef()
    layerName = layer.GetName()
    data.Destroy()
    dstSRS = osr.SpatialReference()
    dstSRS.ImportFromEPSG(int(dstSrsESPG))
    if format == "GeoJSON":
        destDataName = layerName + ".geojson"
        destDataPath = os.path.join(saveFolderPath, destDataName)
    elif format == "ESRI Shapefile":
        destDataName = os.path.join(saveFolderPath, layerName)
        flag = os.path.exists(destDataName)
        os.makedirs(destDataName) if not flag else None
        destDataPath = os.path.join(destDataName, layerName + ".shp")
    else:
        print("不支持该格式！")
        return
    options = gdal.VectorTranslateOptions(
        format=format,
        accessMode=accessMode,
        srcSRS=spatial,
        dstSRS=dstSRS,
        reproject=True,
        selectFields=selectFields,
        layerName=layerName,
        geometryType=geometryType,
        dim=dim
    )
    gdal.VectorTranslate(
        destDataPath,
        srcDS=shapeFilePath,
        options=options
    )
    return destDataPath

if __name__ == '__main__':
    shapeFilePath = r'G:\\vector.shp'
    saveFolderPath = r'G:\\84'
    VectorTranslate(
        shapeFilePath,
        saveFolderPath,
        format="ESRI Shapefile",
        accessMode=None,
        dstSrsESPG=4326,
        selectFields=None,
        geometryType="POLYGON",
        dim="XY",
    )

在上述代码的主程序部分，我们指定了一个Shapefile文件路径和保存目录，并调用了 VectorTranslate 函数将其转换为ESRI Shapefile格式，并重新投影到WGS84坐标系。

在实际应用中，这段代码可以被用来批量处理大量的矢量数据文件，从而简化数据准备过程。例如，当需要将来自不同来源的Shapefile文件统一到一个标准坐标系中时，这样的脚本可以极大地提高工作效率。