python gdal使用教程

GDAL (Geospatial Data Abstraction Library) 是一个用于处理地理空间数据的开源库，它支持多种栅格和矢量地理空间数据格式。在 Python 中，我们可以使用 GDAL 库来读取、写入和处理地理空间数据。

以下是一个简单的 Python GDAL 使用教程，帮助你入门：

1. 安装 GDAL

在开始之前，需要先安装 GDAL。在 Linux 上，可以使用包管理器安装 GDAL，例如：

sudo apt-get install gdal-bin

在 Windows 上，可以下载预编译的 GDAL 安装程序并安装。

2. 安装 Python GDAL 模块

在 Python 中使用 GDAL，需要安装 GDAL 模块。你可以使用 pip 命令来安装 GDAL 模块：

pip install GDAL

3. 读取地理空间数据

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

4. 写入地理空间数据

使用 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 使用教程，希望对你有所帮助。