coordTrans = osr.CoordinateTransformation(sourceSR, targetSR)
geom.Transform(coordTrans)
但是这段代码很让人蛋疼!在windows里面跑不通。老外的论坛里面有讨论,说在linux里面没问题,windows死活不行,哎。。。
另外还有几个要注意的地方:
要在适当的时候编辑Geometry,投影变换之后最好就不要再动了吧。
对一个数据源DataSource里面的所有Geometry做投影变换,你得一个一个来。用个循环吧。
将你的投影写入.prj文件,其实很简单。首先MorphToESRI(),转成字符串,然后开个文本文件往里面写就行了。例如:
targetSR.MorphToESRI()
file = open('test.prj', 'w')
file.write(targetSR.ExportToWkt())
ffile.close()
以上就是python gdal教程之:几何形状geometry与投影projection的内容,更多相关内容请关注PHP中文网(HdhCmsTestgxlcms测试数据)!
查看更多关于pythongdal教程之:几何形状geometry与投影projection的详细内容...
声明:本文来自网络,不代表【好得很程序员自学网】立场,转载请注明出处:http://www.haodehen.cn/did86315