重采样图像rasterio/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'！