FDD-Fear Driven Development

祖国啊,我只是一个摆地摊的.

利用GDAL-python库为SHP格式添加Z值

Z值

z值通常表示一个点的高程值,高程值一般都包含在栅格数据中。本文偿试从栅格数据读出对应点的z值,再写入相应点中。

步骤

  1. 读取SHP格式,获得Geometry,再根据几何类型读出相应的点
  2. 读取栅格数据,跟据已获得相应的点,读取高程值
  3. 新建输出数据源、层。将原图层中的图层定义赋予新建的图层中
  4. 新建特定的Geometry对象,将有z值的点添加到Geometry对象中
  5. 新建Feature对象,并将原有图层的属性表、新建Geometry对象赋予新建Feature对象,

写得很啰嗦,要是有耐心看完wkb的格式,也可以直接操作wkb。

代码

* 在QGIS库中的QgsRasterLayer有个很简单的方法就可以提出z值—— rasterLayer.identify(point),也可以用GDAL-python的api读取z值

阅读更多
文章标签: float api
个人分类: Python
想对作者说点什么? 我来说一句

python 读取shp文件

2013年03月14日 2KB 下载

没有更多推荐了,返回首页

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭