赞
踩
目录
GDAL是应用广泛的空间数据处理库,可以处理几何、栅格数据,Python是一门简单易学的编程语言,常用来编写数据处理工具、脚本。本文讲解如何在Python环境中安装、配置、使用GDAL。本文示例中使用的GDAL版本为3.4.3
下载链接:https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal,下载所需要版本GDAL的whl包,如下图所示,这里需要注意,下载的whl包版本需要与自己所使用的Python版本、位数一致。我这里使用的是64位版的Python3.8,因此我这里下载的是:GDAL-3.4.3-cp38-cp38-win_amd64.whl。
需要本教程所用whl包及示例代码、示例数据,请在微信公众号“GIS工具乐园”,回复“GDAL343”获取。
图1 下载GDAL
在安装whl包前,需确保Python中已安装pip。这里将在cmd命令行中使用pip进行whl包的安装。
使用“cd /d”命令将cmd的工作目录切换至Python安装路径下的“Scripts”目录,例如:cd /d D:\Program Files\Python38\Scripts,如下图所示。
图2 切换安装目录
使用“pip install”命令进行whl包安装,例如:pip install D:\GDAL-3.4.3-cp38-cp38-win_amd64.whl,如下图所示。
图3 安装GDAL
通过PyCharm可以看到,GDAL已成功安装,如下图所示。
图4 安装效果
使用GDAL读取shapefile,并打印其wkt坐标串、图形面积,代码及效果如下所示。关注微行公众号:“GIS工具乐园”,持续更新开发技巧分享哦。
- strVectorFile =r"C:Users\HX\Desktop\test\示例数据\示例数据.shp"
-
-
- dataSource = ogr.Open(strVectorFile, 0)
- layer = dataSource.GetLayerByIndex(0)
- layer.ResetReading()
- featureNum = layer.GetFeatureCount()
- i = 0
- while i < featureNum:
- i += 1
- feature = layer.GetNextFeature()
- geom = feature.GetGeometryRef()
- wkt = geom.ExportToWkt()
- area = geom.Area()
- print("图形" + str(i) + "的wkt:" + wkt)
- print("图形" + str(i) + "的面积:" + str(area))
图5 使用GDAL读取shp
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。