用python读取删格图像中一部分，并进行重采样

Question

用python，使用 GDAL 读取 2012_test.tif 的第 1000 行至11000 行和第 1000 列至11000 列，压缩成 2500*2500 尺寸，显示出来，并存储到新数据集。要求具体代码。

Answer 1

from osgeo import gdal, gdalconst

# 输入和输出影像文件路径
input_raster = "2012_test.tif"
output_raster = "resampled_image.tif"

# 打开影像数据集
dataset = gdal.Open(input_raster, gdalconst.GA_ReadOnly)

# 定义裁剪区域
start_col, start_row, num_cols, num_rows = 1000, 1000, 10000, 10000

# 读取部分影像数据
band = dataset.GetRasterBand(1)  # 获取第一个波段
data = band.ReadAsArray(start_col, start_row, num_cols, num_rows)

# 计算重采样后的大小
resampled_cols, resampled_rows = 2500, 2500

# 创建输出数据集
driver = gdal.GetDriverByName('GTiff')
out_data = driver.Create(output_raster, resampled_cols, resampled_rows, 1, gdal.GDT_Float32)

# 计算重采样的缩放比例
resample_x = num_cols / resampled_cols
resample_y = num_rows / resampled_rows

# 重采样数据
out_data.GetRasterBand(1).WriteArray(data[::resample_y, ::resample_x])

# 释放资源
out_data.FlushCache()
out_data = None
dataset = None