基于osgeo获取矢量文件坐标系和矢量文件坐标系变换

John H.

已于 2022-09-29 17:16:56 修改

阅读量528

点赞数

分类专栏： OSGEO 文章标签： python

于 2022-03-24 09:14:08 首次发布

本文链接：https://blog.csdn.net/weixin_42778485/article/details/123702266

版权

本文介绍了如何利用Python的osgeo模块获取矢量文件的坐标系信息，并详细阐述了进行坐标系转换的方法，包括坐标系的理解、坐标转换的原理以及实际操作步骤。

摘要由CSDN通过智能技术生成

from osgeo import ogr, gdal, gdal_array, ogr, osr



def obtain_corrdinate_inform(_shapefile):
    '''获取矢量文件的坐标系，返回的是EPSG编码'''
    _ds = ogr.Open(_shapefile)
    _lyr = _ds.GetLayer()
    _spatial = _lyr.GetSpatialRef()
    return int(_spatial.GetAttrValue('AUTHORITY',1))


def spatial_reference_change(_need_to_change_shapefile, _changed_corrdinate_shapefile, _target_coordinate_shapefile):
    '''坐标系变换
    _target_coordinate_shapefile目标坐标系矢量文件
    _need_to_change_shapefile需要转换坐标系的矢量文件
    _chagned_corrdinate_shapefile坐标系转换后的矢量文件'''
    # _name = os.path.basename(_need_to_change_shapefile).split('.')[0]
    # _changed_corrdinate_shapefile = os.path.join(os.path.dirname(_need_to_change_shapefile), (_name + '_changed_coordinate.shp'))
    shpdriver = ogr.GetDriverByName('ESRI Shapefile')
    if os.path.exists(_changed_corrdinate_shapefile):