重新采样图像栅格/ gdal，Python

Question

如何使用双线性插值对单个波段的GeoTIFF进行重新采样？

import os

import rasterio
from rasterio.enums import Resampling
from rasterio.plot import show,show_hist
import numpy as np

if __name__ == "__main__":
    

    input_Dir = 'sample.tif'
    #src = rasterio.open(input_Dir) 
    #show(src,cmap="magma")

    upscale_factor = 2

    with rasterio.open(input_Dir) as dataset:
        
        # resample data to target shape
        data = dataset.read(
            out_shape=(
                dataset.count,
                int(dataset.height * upscale_factor),
                int(dataset.width * upscale_factor)
            ),
            resampling=Resampling.bilinear
        )

        # scale image transform
        transform = dataset.transform * dataset.transform.scale(
            (dataset.width / data.shape[-1]),
            (dataset.height / data.shape[-2])
        )
        show(dataset,cmap="magma",transform=transform)

我尝试了以下代码，我的输出如下：

我正在尝试实现以下输出：

Answer 1

一种选择是使用GDAL python绑定。然后，您可以在内存中执行重新采样（或者如果需要，可以保存图像）。假设旧的栅格分辨率为0.25x0.25，而您将重新采样为0.10x0.10：

from osgeo import gdal

input_Dir = 'sample.tif'

ds = gdal.Translate('', input_Dir, xres=0.1, yres=0.1, resampleAlg="bilinear", format='vrt')

如果要保存图像，则将输出文件路径而不是第一个参数的空字符串放入，并将格式更改为“ tif”！