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

代码关键点

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')
driver