【实用的GDAL】使用 gdalwarp 影像裁剪和重采样

主要功能

Image reprojection and warping utility

影像重采样和裁剪

全部参数

gdalwarp [--help-general] [--formats]
    [-s_srs srs_def] [-t_srs srs_def] [-ct string] [-to "NAME=VALUE"]* [-vshift | -novshift]
    [[-s_coord_epoch epoch] | [-t_coord_epoch epoch]]
    [-order n | -tps | -rpc | -geoloc] [-et err_threshold]
    [-refine_gcps tolerance [minimum_gcps]]
    [-te xmin ymin xmax ymax] [-te_srs srs_def]
    [-tr xres yres] [-tap] [-ts width height]
    [-ovr level|AUTO|AUTO-n|NONE] [-wo "NAME=VALUE"] [-ot Byte/Int16/...] [-wt Byte/Int16]
    [-srcnodata "value [value...]"] [-dstnodata "value [value...]"]
    [-srcalpha|-nosrcalpha] [-dstalpha]
    [-r resampling_method] [-wm memory_in_mb] [-multi] [-q]
    [-cutline datasource] [-cl layer] [-cwhere expression]
    [-csql statement] [-cblend dist_in_pixels] [-crop_to_cutline]
    [-if format]* [-of format] [-co "NAME=VALUE"]* [-overwrite]
    [-nomd] [-cvmd meta_conflict_value] [-setci] [-oo NAME=VALUE]*
    [-doo NAME=VALUE]*
    srcfile* dstfile

常用参数介绍

[-s_srs]原始影像的坐标系，如未指定，默认从影像中查找坐标系
[-t_srs]设置目标坐标系
[-tr]定义生成结果的分辨率，命令行中直接写 x、y 方向的分辨率，中间使用空格隔开
[-r]重采样方法，可用的方法有：near， bilinear，cubic，cubicspline，lanczos，average，rms，mode，max，min，med，q1，q3，sum
[-of]选择输出格式，从 GDAL2.3 开始，如未指定，则从扩展名中猜测格式（以前是 GTiff ）。
[-co]许多命令有一个或多个选项，此参数控制有关命令的详细选项。例：生成影像压缩方式的选择：
-co COMPRESS=[JPEG/ LZW /PACKBITS/DEFLATE/CCITTRLE/CCITTFAX3/CCITTFAX4/LZMA/ZSTD/LERC/LERC_DEFLATE/LERC_ZSTD/WEBP/JXL/NONE]（通常选择 LZW 方式）
[-cutline]使用 OGR 支持数据源的混合切割线
[-crop_to_cutline]将目标裁减到切割线范围
源文件路径，可以输入多个，每个文件使用空格隔开
生成文件路径，如果输出文件已经存在，支持拼接到现有输出文件中

更多参数使用方法：GDAL文档

使用例子

影像裁剪

在命令行中输入：

gdalwarp -tr 10 10 -co COMPRESS=LZW -crop_to_cutline -cutline /your_file.geojson /image_need_to_be_cropped.tif /save_image.tif

影像重采样

gdalwarp -tr x y -r near -co COMPRESS=LZW /your_images.tif* /save_image.tif

2022.12.22补充

此命令类似于可执行文件，之前写的是在 Linux 系统中运行，如果需要在 Win 下使用，可以在 gdal 的 pip 安装文件中找到相关程序，我的可执行文件路径：
D:\anaconda3\envs\pytorch\Lib\site-packages\osgeo\gdalwarp.exe
Linux 中如果不能直接调用 gdalwarp 命令，可以使用
/usr/local/lib/python3.6/dist-packages/osgeo/gdalwarp ，若之后要经常使用，则可以在 .bashrc 文件中建立快捷命令

echo "alias gdalwarp='/usr/local/lib/python3.6/dist-packages/osgeo/gdalwarp'" >> ~/.bashrc
source ~/.bashrc