GDAL自带的 rpc纠正和金字塔文件生成方法

GDAL纠正功能很强大，这里只做一些简单的纠正。
说明：首先切换到GDAL的bin目录下，然后运行下面命令

-multi 多线程

1、带DEM影像纠正的方法

gdalwarp -multi -rpc -to "RPC_DEM=D:\datas\dem\DEM-1000.tif" D:\datas\Rectify_TestData\datas\GF3_KAS_FSII_023906_W72.1_N6.3_20210223_L1A_HH_L10005498031intensity.tiff D:\datas\Rectify_TestData\datas_result\GF3exeWithDem_rectify.tiff

带入DEM影像，纠正结果可能不规则，如下图：

2、DEM为常数的纠正方法

gdalwarp -multi -rpc -to "RPC_HEIGHT=500" D:\datas\InDoc4Image\GF3_KAS_FSII_020490_W76.4_S8.9_20200701_L1A_VV_L10004898896intensity.tiff D:\datas\InDoc4Image_result\GF3exeWithDem500_rectify.tiff

输出结果是个完整的矩形，如下图：

还可以设置的参数：

-to "RPC_DEM_INTERPOLATION=CUBIC"
-to "RPC_DEM_MISSING_VALUE=0" 

-wo USE_OPENCL=FALSE

另外一个样例，可以设置输出影像的格式

gdalwarp -multi -tap -et 0.0 -r cubic -srcnodata 0 -rpc -overwrite -to "RPC_DEM_INTERPOLATION=CUBIC" -to "RPC_DEM_MISSING_VALUE=0" -to "RPC_DEM=D:\Temp\bug_gdalwarp\DEM.tif" -tr 30 30 -t_srs "+proj=aea +lat_1=44.75 +lat_2=55.75 +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs" D:\Temp\bug_gdalwarp\IMG.tif D:\Temp\bug_gdalwarp\IMG_ortho.tif

下面是直接纠正的命令，没用设置DEM

echo  %time% 
F:\gdal1112\gdalwarp.exe -rpc -t_srs WGS84 -r cubic -wm 1024 -et 0 -wo USE_OPENCL=FALSE F:\Data\711_214_26sep2006_p5\P5\banda.tif F:\gdal1112.tif

echo  %time% 
F:\gdal20\gdalwarp.exe -rpc -t_srs WGS84 -r cubic -wm 1024 -et 0 -wo USE_OPENCL=FALSE F:\Data\711_214_26sep2006_p5\P5\banda.tif F:\gdal21dev.tif

echo  %time% 

用RPC仿射六参数来纠正

gdalwarp.exe -rpc -to "RPC_AFFINE=-32.714672501057066 0.999199897235577 0.000158731686899 28.720843336473692 0.000589585516339 1.000068008511035" D:\rpctest\banda.tif D:\rpctest\banda_affine.tif

3、借助GDAL工具快速生成金字塔文件

只用将影像路径更换成待生成金字塔的影像，2 4 8 16 32 为金字塔级别列表（层数），后面可以继续乘以2进行累计

命令

gdaladdo.exe --config USE_RRD YES D:\te\Test16To8\GF6_PMS_E116.8_N40.2_20190415_L1A1119868086-PAN.tiff 2 4 8 16 32

4、gdalwarp命令解释

4.1、所有命令

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

4.2、命令含义

[-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]将目标裁减到切割线范围
<srcfile*>源文件路径，可以输入多个，每个文件使用空格隔开
生成文件路径，如果输出文件已经存在，支持拼接到现有输出文件中