GDAL (Geospatial Data Abstraction Library) 是一个用于处理地理空间数据的开源库，它支持多种栅格和矢量地理空间数据格式。在 Python 中，我们可以使用 GDAL 库来读取、写入和处理地理空间数据。
以下是一个简单的 Python GDAL 使用教程，帮助你入门：
- 安装 GDAL
在开始之前，需要先安装 GDAL。在 Linux 上，可以使用包管理器安装 GDAL，例如：
sudo apt-get install gdal-bin
在 Windows 上，可以下载预编译的 GDAL 安装程序并安装。
- 安装 Python GDAL 模块
在 Python 中使用 GDAL，需要安装 GDAL 模块。你可以使用 pip 命令来安装 GDAL 模块：
pip install GDAL
- 读取地理空间数据
使用 GDAL 可以读取多种地理空间数据格式，例如 GeoTIFF、ESRI Shapefile 等。以下是一个简单的示例代码，演示如何使用 GDAL 读取 GeoTIFF 文件：
from osgeo import gdal
# 打开 GeoTIFF 文件
dataset = gdal.Open('path/to/file.tif')
# 获取像素宽度和高度
cols = dataset.RasterXSize
rows = dataset.RasterYSize
# 读取指定位置的像素值
pixel_value = dataset.ReadAsArray(0, 0, cols, rows)
# 输出像素值
print(pixel_value)
- 写入地理空间数据
使用 GDAL 可以将数据写入多种地理空间数据格式。以下是一个简单的示例代码，演示如何使用 GDAL 将数据写入 GeoTIFF 文件：
from osgeo import gdal, gdal_array
# 设置输出文件名和驱动程序
output_filename = 'path/to/output.tif'
driver = gdal.GetDriverByName('GTiff')
# 设置输出文件的大小和投影信息
cols = 128
rows = 128
projection = gdal.osr.SpatialReference()
projection.ImportFromEPSG(4326)
# 创建输出文件
output_dataset = driver.Create(output_filename, cols, rows, 1, gdal_array.NumericTypeCode.Uint8)
output_dataset.SetProjection(projection.ExportToWkt())
# 写入数据
data = gdal_array.numpy.zeros((rows, cols), gdal_array.NumericTypeCode.Uint8)
output_dataset.GetRasterBand(1).WriteArray(data)
# 关闭输出文件
output_dataset = None
以上是一个简单的 Python GDAL 使用教程，希望对你有所帮助。
