python使用gdal实现影像重采样

gdal实现影像的重采样

对于tif格式的遥感影像而言，通常使用左上角那个像元真实的x坐标值，像素宽度，x像素旋转角度，原点y坐标，y像素旋转角度，像素高度（负值）共六个参数来表示整幅影像像素点的坐标。

当影像分辨率改变后，x坐标值和y坐标值度都不会改变，但像素高度和宽度会改变，而像素高度和宽度改变的话，影像的行列数也会改变。

1.gdal实现影像重采样至高分辨率

这里，将一幅影像的像元高度和宽度都变为原来的0.5倍，实现代码如下：

'''
Created on 2020年2月10日

@author: Sun Strong
'''
from osgeo import gdal
import os
path=r"C:\Users\Sun Strong\Desktop\gdal" 
os.chdir(path)
in_ds=gdal.Open("p047r027_7t20000730_z10_nn10.tif")
in_band=in_ds.GetRasterBand(1)
xsize=in_band.XSize
ysize=in_band.YSize
geotrans=list(in_ds.GetGeoTransform())
geotrans[1]/=2#像元宽度变为原来的0.5倍
geotrans[5]/=2#像元高度变为原来的0.5倍
#重采样后的影像
if os.path.exists('resampled.tif'):#如果已存在同名影像
    os