GDAL是应用广泛的空间数据处理库，可以处理几何、栅格数据，Python是一门简单易学的编程语言，常用来编写数据处理工具、脚本。本文讲解如何在Python环境中安装、配置、使用GDAL。本文示例中使用的GDAL版本为3.4.3

一、下载GDAL的whl包
可以通过链接：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”获取。

二、安装GDAL
在安装whl包前，需确保Python中已安装pip。这里将在cmd命令行中使用pip进行whl包的安装。
1.切换安装目录
使用“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.查看安装效果
通过PyCharm可以看到，GDAL已成功安装，如下图所示。

 

三、读取shapefile测试
使用GDAL读取shapefile，并打印其wkt坐标串、图形面积，代码及效果如下所示。

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))