GDAL python读取shp进行影像投影坐标和地理坐标的转换

最新推荐文章于 2025-07-21 09:12:43 发布

Fred240011

最新推荐文章于 2025-07-21 09:12:43 发布

阅读量1.1w

点赞数 4

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/Fred240011/article/details/81217472

本文介绍如何使用gdal将.shp文件中的投影坐标转换为地理坐标（经纬度），以便于在mongodb中建立空间索引。通过示例代码演示了整个转换过程。

之前只在网上找到了使用gdal进行tif格式影像的投影坐标与地理坐标转换的代码，但没有专门针对shp格式的数据的教程。正巧现在需要用mongodb做矢量入库，建立空间索引（2dshere）时必须用经纬度而不是投影坐标，就参考别人的代码做了一些修改，实现了shp的投影坐标转地理坐标。

参考博客地址：https://blog.youkuaiyun.com/theonegis/article/details/54427906

from osgeo import gdal
from osgeo import osr
from osgeo import ogr

#输入为待操作shp的路径，和一组待转换的投影坐标
def prj2geo(path,x,y):
    # 为了支持中文路径，请添加下面这句代码
    gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "NO")
    # 为了使属性表字段支持中文，请添加下面这句
    gdal.SetConfigOption("SHAPE_ENCODING", "")
    # 注册所有的驱动
    ogr.RegisterAll()
    # 数据格式的驱动
    driver = ogr.GetDriverByName('ESRI Shapefile')
    ds = driver.Open(path);
    layer0 = ds_simple.GetLayerByIndex(0);

    #或取到shp的投影坐标系信息
    prosrs = layer0.GetSpatialRef()
    geosrs = osr.SpatialReference()
    #设置输出坐标系为WGS84
    geosrs.SetWellKnownGeogCS("WGS84")

    ct = osr.CoordinateTransformation(prosrs, geosrs)
    coords = ct.TransformPoint(x, y)
    #输出为转换好的经纬度
    return coords[:2]