python将tif遥感影像行列数变为1024*1024而像素大小不变

如果你需要将tif遥感影像的行列数变为1024x1024而不改变像素大小，可以使用Python的GDAL库来实现。以下是一个示例代码：

from osgeo import gdal

# 打开遥感影像
ds = gdal.Open('your_image.tif')

# 获取遥感影像的行列数和波段数
rows = ds.RasterYSize
cols = ds.RasterXSize
bands = ds.RasterCount

# 获取遥感影像的投影和地理变换信息
proj = ds.GetProjection()
trans = ds.GetGeoTransform()

# 计算新的列数和行数
new_cols = 1024
new_rows = 1024

# 计算每个像素在新的影像中的大小
pixel_size_x = trans[1] * cols / new_cols
pixel_size_y = trans[5] * rows / new_rows

# 计算新的地理变换信息
new_trans = (trans[0], pixel_size_x, trans[2], trans[3], trans[4], -pixel_size_y)

# 创建输出遥感影像
driver = gdal.GetDriverByName('GTiff')
out_ds = driver.Create('your_output_image.tif', new_cols, new_rows, bands, gdal.GDT_Float32)

# 设置输出遥感影像的投影和地理变换信息
out_ds.SetProjection(proj)
out_ds.SetGeoTransform(new_trans)

# 读取每个波段的数据，并将影像的行列数变为1024x1024
for b in range(bands):
    band = ds.GetRasterBand(b+1)
    data = band.ReadAsArray(0, 0, cols, rows)
    out_band = out_ds.GetRasterBand(b+1)
    out_band.WriteArray(data)

# 关闭数据集
ds = None
out_ds = None

这个示例代码将遥感影像的行列数都变为1024x1024，但是像素大小保持不变。注意，新的列数和行数需要与原始影像的宽高比例相同，否则会导致影像形状发生变化。新的地理变换信息的计算方法是根据新的列数和行数计算每个像素在新的影像中的大小，并且保持原始影像的左上角坐标不变。