python调用gdal的代码

在Python中，GDAL（Geospatial Data Abstraction Library）是一个用于处理各种地理空间数据的库，包括图像和矢量文件。以下是一个简单的例子，展示如何使用GDAL读取和显示GeoTIFF图像：

from osgeo import gdal

# 打开GeoTIFF文件
ds = gdal.Open('path_to_your_file.tif')

# 获取数据集的一些基本信息
band_count = ds.RasterCount   # 图像有多少条光谱带
width = ds.RasterXSize       # 图像宽度
height = ds.RasterYSize      # 图像高度
projection = ds.GetProjection()  # 图像投影信息
geotransform = ds.GetGeoTransform()  # GPS坐标转换矩阵

# 访问并显示第一幅图层的数据
band1 = ds.GetRasterBand(1)
image_data = band1.ReadAsArray(0, 0, width, height)  # 读取像素值
print(image_data)

# 关闭数据集
ds = None