python gdal写geotif

Python GDAL是一个用于处理地理空间数据的库，它提供了对数据集的读取、写入、转换等功能。gdal包括两个核心模块：gdal和osr。

使用Python GDAL库可以方便地写Geotiff文件。首先，需要导入gdal模块，并创建一个新的数据集：

from osgeo import gdal

# 创建新的Geotiff数据集
driver = gdal.GetDriverByName("GTiff")
dataset = driver.Create("new.tif", width, height, bands, datatype)

在上述代码中，"GTiff"表示创建的数据集为Geotiff格式。接下来，可以对数据集进行写入操作。例如，写入数据到某个波段：

# 获取指定波段
band = dataset.GetRasterBand(band_number)

# 将数据写入波段
band.WriteArray(data)

其中，"band_number"表示波段的索引号，从1开始。"data"为要写入的数据数组。

写入数据后，可以设置Geotiff文件的空间参考信息（空间坐标系统、投影、地理变换等）：

# 获取数据集的空间参考
srs = osr.SpatialReference()
srs.ImportFromEPSG(epsg_code)

# 设置地理变换信息
dataset.SetGeoTransform(geotransform)

# 设置投影信息
dataset.SetProjection(srs.ExportToWkt())

在上述代码中，"epsg_code"表示要设置的坐标系的EPSG码，"geotransform"是一个6个值的元组，包括地理变换信息。

最后，将数据集写入到Geotiff文件中：

# 将数据集写入文件
dataset.FlushCache()

使用Python GDAL库的gdal包，可以方便地创建和写入Geotiff文件，同时也可以设置空间参考信息。这样就可以使用GDAL库在Python中处理地理空间数据并输出为Geotiff格式。