GDAL空间参考坐标——Python地理数据处理学习分享

一、获取地理坐标或投影坐标

通过GetSpatialRef()函数来获取图层的坐标系统

# 调用模板
from osgeo import ogr
from osgeo import osr
# 打印SRS(未投影)
ds=ogr.Open('D:\\Python_study\\osgeopy-data\\US\\states_48.shp')
print(ds.GetLayer().GetSpatialRef())

利用osr模块，创建一个空的SpatialReference对象，利用ImportFromEPSG（）函数导入投影。

# 使用标准的EPSG代码
sr=osr.SpatialReference()
sr.ImportFromEPSG(26912)
print(sr)

利用GetAttrValue()函数获取各种值。

其中PROJCS为投影和基准信息

# 获取SRS的基准面和投影
jz=sr.GetAttrValue('PROJCS')
print(jz)

二、自定义投影

创建一个空的SpatialReference对象，设置一个投影（采取SetProjCS()函数），设置一个基准（采用SetWellKnownGeogCS()函数），利用SetACEA()函数创建空间参考。

SetACEA（a,b,c,d,e,f)函数中a\b\c\d\e\f六个参数分别为标准平行线1、标准平行线2、中央纬线、中央经线、东移假定值、北移假值。

# 自定义投影
sr=osr.SpatialReference()
sr.SetProjCS('USGS Albers')
sr.SetWellKnownGeogCS('NAD83')
sr.SetACEA(29.5,45.5,23,-96,0,0)
sr.Validate()
print(sr)

三、创建一个带坐标系统的图层

# 创建带坐标的图层
import os 
from osgeo import ogr
from osgeo import osr
sr=osr.SpatialReference()
sr.ImportFromEPSG(26912)
driver=ogr.GetDriverByName('ESRI Shapefile')
ds=driver.CreateDataSource('D:\\Python_study\\python\\5c\\test.shp')
lyr=ds.CreateLayer('counties',sr,ogr.wkbPolygon)
print(ds.GetLayer().GetSpatialRef())

四、重投影一个数据

from osgeo import osr
from osgeo import ogr
sr=osr.SpatialReference()
sr.ImportFromEPSG(3857)
ds=ogr.Open('D:\\Python_study\\osgeopy-data\\US',1)
in_lyr=ds.GetLayer('us_volcanos')
out_lyr=ds.CreateLayer('us_volcanos_aea',sr,ogr.wkbPoint)
out_lyr.CreateFields(in_lyr.schema)
out_feat=ogr.Feature(out_lyr.GetLayerDefn())
for in_feat in in_lyr:
	geom=in_feat.geometry().Clone()
	geom=TransformTo(sr)
	out_feat.SetGeometry(geom)
	for i in range(in_feat.GetFieldCount()):
		out_feat.SetField(i,in_feat.GetField(i))
	out_lyr.CreateFeature(out_feat)