GDAL工具应用–007–生成金字塔与生成COG

大家好，欢迎来到IT知识分享网。

参考：
https://gdal.org/download.html

1.生成金字塔

命令： gdaladdo.exe -r average D:\data\corp.tif -ro 2 4 8 16 32 64 128 256

其中的 -ro 参数表示按1/2, 1/4, 1/8, 1/16, 1/32, 1/64, 1/128, 1/256进行重采样

C:\Users\Administrator>gdaladdo.exe -r average D:\testData\corp\xiaomai.tif -ro 2 4 8 16 32 64 128 256 0...10...20...30...40...50...60...70...80...90...100 - done.

2.将栅格文件转换到COG

2.1投影变换法

第一种方法：投影变换到EPSG:4490坐标系，并在geotiff中构建tile

1.如果原始栅格本身是EPSG:4490坐标系，则只在geotiff中构建tile

命令：gdal_translate.exe D:/data/GF2DOM_4490.tif
D:/data/GF2DOM_4490_cog.tif -co TILED=YES -co COMPRESS=DEFLATE -co COPY_SRC_OVERVIEWS=YES

2.如果原始栅格本身不是EPSG:4490坐标系，则需要进行投影变换，并在结果geotiff中构建tile

命令：gdalwarp.exe -t_srs EPSG:4490 -r bilinear -wm 4096 -et 0 -wo USE_OPENCL=FALSE -co BIGTIFF=YES -co TILED=YES -co COMPRESS=DEFLATE -co COPY_SRC_OVERVIEWS=YES -multi D:/data/GF2DOM_4326.tif D:/data/GF2DOM_4490.tif

C:\Users\Administrator> gdalwarp.exe -t_srs EPSG:4490 -r bilinear -wm 4096 -et 0 -wo USE_OPENCL=FALSE -co BIGTIFF=YES  -co TILED=YES -co COMPRESS=DEFLATE -co COPY_SRC_OVERVIEWS=YES -multi D:/data/GF2DOM_4326.tif D:/data/GF2DOM_4490.tif 0...10...20...30...40...50...60...70...80...90...100 - done.

2.2Tile策略法

第二种方法：在结果geotiff中指定输出为COG，并且指定Tile策略为GoogleMapsCompatible参数

命令：gdalwarp.exe -t_srs EPSG:4326 -r bilinear -wm 4096 -et 0 -wo USE_OPENCL=FALSE -co BIGTIFF=YES -co TILING_SCHEME=GoogleMapsCompatible -co COMPRESS=DEFLATE -co OVERVIEWS=IGNORE_EXISTING -of COG -multi D:/data/test_5678.tif D:/data/test_4326.tif

其中：

-of COG 参数表示输出为COG

TILING_SCHEME=GoogleMapsCompatible参数表示Tile策略为GoogleMapsCompatible

C:\Users\Administrator>gdalwarp.exe -t_srs EPSG:4326 -r bilinear -wm 4096 -et 0 -wo USE_OPENCL=FALSE -co BIGTIFF=YES -co TILING_SCHEME=GoogleMapsCompatible -co COMPRESS=DEFLATE -co OVERVIEWS=IGNORE_EXISTING -of COG -multi D:/data/test_5678.tif D:/data/test_4326.tif 0...10...20...30...40...50...60...70...80...90...100 - done.