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


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

参考博客地址:https://blog.csdn.net/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]

转载自:https://blog.csdn.net/Fred240011/article/details/81217472

You may also like...