gdal readraster和writeraster的用法

gdal是一个用于地理空间数据处理的开源库，其中包含了读写栅格数据的函数，其中readraster和writeraster是其中的两个函数，分别用于读取和写入栅格数据。

下面是它们的用法：

读取栅格数据：

from osgeo import gdal

# 打开栅格数据
dataset = gdal.Open('path/to/raster.tif', gdal.GA_ReadOnly)

# 获取栅格数据的基本信息
width = dataset.RasterXSize
height = dataset.RasterYSize
bands = dataset.RasterCount
projection = dataset.GetProjection()
geotransform = dataset.GetGeoTransform()

# 读取栅格数据
band = dataset.GetRasterBand(1)
data = band.ReadAsArray(0, 0, width, height)

# 关闭栅格数据
dataset = None

写入栅格数据：

from osgeo import gdal

# 创建栅格数据
driver = gdal.GetDriverByName('GTiff')
dataset = driver.Create('path/to/output.tif', width, height, bands, gdal.GDT_Float32)

# 设置栅格数据的基本信息
dataset.SetGeoTransform(geotransform)
dataset.SetProjection(projection)

# 写入栅格数据
band = dataset.GetRasterBand(1)
band.WriteArray(data)

# 关闭栅格数据
dataset = None

需要注意的是，读取和写入栅格数据时，需要根据实际情况进行修改。例如，读取时可以根据需要读取多个波段的数据，写入时可以根据需要设置不同的数据类型。