一、获取地理坐标或投影坐标
通过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)