Python+OGR库学习（四）：重投影shp文件并另存，属性表保持不变

最新推荐文章于 2023-10-10 09:17:54 发布

Learner_P

最新推荐文章于 2023-10-10 09:17:54 发布

阅读量3.1k

点赞数 3

本文链接：https://blog.csdn.net/weixin_42924891/article/details/85469881

版权

本文介绍了如何使用Python结合OGR库进行shp文件的重投影操作，详细阐述了代码关键点，包括转换参数设定、几何对象坐标转换、属性表保持不变的方法。通过定义转换关系，打开输入文件，读取图层并创建输出文件，遍历并转换每个要素，最终保存结果。此外，还补充了在处理过程中应注意的细节，如字段属性的读取和空间参考文件的创建。

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

代码关键点

1、首先要定义好转换参数
2、主要操作对象是要素，需要提前创建好输出文件，然后遍历所有要素，对每一个几何对象进行坐标转换
3、输出文件的字段属性定义需要从输入文件读取

代码思路

1、导入相关包，切换路径，注册驱动
2、定义转换关系
3、打开输入文件，读取到图层，创建输出文件，定义字段属性
4、遍历要素，循环读取要素、转换投影、创建对应新要素写入几何，并添加要素到图层
5、定义空间参考文件.prj

代码

#! usr/bin/env python3
# -*- coding:utf-8 -*-
#导入相关包
import ogr,os,osr
#定义投影转换关系
inosr = osr.SpatialReference()
inosr.ImportFromEPSG(4269)

outosr = osr.SpatialReference()
outosr.ImportFromEPSG(26912)

trans = osr.CoordinateTransformation(inosr,outosr)
#读取矢量文件，获取图层
os.chdir(r'F:\Python+gdal\7weeks数据\7weeks数据\ospy_data2\ospy_data2'