Python空间数据处理----矢量数据在不同空间参考下的转换

最新推荐文章于 2024-06-11 16:17:59 发布

x**3

最新推荐文章于 2024-06-11 16:17:59 发布

阅读量2.9k

点赞数 6

分类专栏： Python与GIS数据处理

本文链接：https://blog.csdn.net/qq_36181130/article/details/83543544

版权

本文介绍了Python中进行空间数据处理时，如何利用GDAL/OGR库进行矢量数据在不同空间参考之间的转换。通过使用SpatialReference、CoordinateTransformation等函数，详细阐述了从EPSG导入空间参考、创建数据源、添加图层、转换几何坐标以及保存转换后的数据的步骤。

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

流程图如下：
需要用到的函数：
（1）spatialref=osr.SpatialReference()
该函数是SpatialReference类的一个构造函数，它构造一个 SpatialReference对象，另外，它接受wkt格式的字符串来构建SpatialReference对象，如果没有接受，可以通过SpatialReference对象调用成员函数ImportFromWkt()、ImportFromEPSG()来构建。
（2）spatialref.ImportFromEPSG(4326)
该函数从EPSG中导入一个空间参考。
（3）trans=osr.CoordinateTransformation(spatialref_source,spatialref_target)
该函数是CoordinateTransformation（坐标系统转换）类的构造函数，用来创建坐标系统转换对象，后续可以通过该函数创建的坐标系统转换对象来对数据的空间参考进行转换。
（4）driver=ogr.GetDriverByName(“ESRI Shapefile”)
注册驱动，可以理解为：这一个函数的作用是告诉程序需要在磁盘中创建多少个格式的文件，如：对于 .shp文件，这一函数则告诉程序需要在文件中创建 .shp、.shx、.prj、.dbf等文件，创建完之后，可以将数据写入相应的文件中。该函数返回一个Driver类的对象。