Python库：gdal

---

GDAL（Geospatial Data Abstraction Library）是处理地理空间数据的开源库，广泛用于地理信息系统（GIS）和遥感领域。GDAL 提供了读取、写入、转换和处理栅格数据（如卫星图像、DEM等）和矢量数据（如Shapefile）的功能。

Python 通过 gdal 库提供了对 GDAL 的绑定，使得开发者可以在 Python 环境中使用 GDAL 的功能。

1. 安装 gdal 库

在安装 gdal 库之前，需要确保已经安装了 GDAL 库本身。然后可以通过 pip 安装 Python 的 gdal 绑定：

pip install gdal

2. 基本概念

在使用 gdal 库之前，了解一些基本概念是很有帮助的：

• 栅格数据：栅格数据是由像素组成的图像，每个像素都有一个值，通常表示某种测量值（如高度、温度、反射率等）。
• 矢量数据：矢量数据由点、线、面等几何对象组成，通常用于表示地理特征（如道路、建筑物、河流等）。
• 投影：地理数据通常需要投影到二维平面上，常见的投影方式包括 UTM、WGS84 等。
• 坐标系：地理数据通常使用特定的坐标系来表示位置，如 WGS84（经纬度）、UTM 等。

3. 读取栅格数据

gdal 库可以读取多种栅格数据格式，如 GeoTIFF、JPEG、PNG 等。以下是一个简单的示例，展示如何读取一个 GeoTIFF 文件并获取其基本信息。

from osgeo import gdal

# 打开栅格数据文件
dataset = gdal.Open('example.tif')

# 获取栅格数据的宽度和高度
width = dataset.RasterXSize
height = dataset.RasterYSize

# 获取栅格数据的波段数量
bands = dataset.RasterCount

# 获取栅格数据的投影信息
projection = dataset.GetProjection()

# 获取栅格数据的地理变换信息（仿射变换参数）
geotransform = dataset.GetGeoTransform()

# 获取第一个波段的数据
band = dataset.GetRasterBand(1)
data = band.ReadAsArray()

# 打印基本信息
print(f"Width: {
     width}, Height: {
     height}, Bands: {
     bands}")
print(f"Projection: {
     projection}")
print(f"GeoTransform: {
     geotransform}")
print(f"First band data shape: {
     data.shape}")

4. 读取矢量数据

gdal 库也可以读取矢量数据，如 Shapefile。以下是一个简单的示例，展示如何读取一个 Shapefile 并获取其几何对象和属性。

from osgeo import ogr

# 打开矢量数据文件
driver = ogr.GetDriverByName('ESRI Shapefile')
dataSour