对GEE下载时分块影像进行镶嵌(GDAL)
前言
GDAL自带的镶嵌功能比较单一,只能将图像简单拼接到一起,不能实现直方图匀色以及羽化等功能,但是对GEE导出的分块影像进行镶嵌较为适合。
优点:
- 使用代码,镶嵌多个影像较为简便
- GDAL较为稳定
- 镶嵌时NAN的区域不会改为0或者其他自定义的值,仍为NAN
缺点:
- 速度较ENVI慢
- 镶嵌进度不明显
- 镶嵌过程中占用大量硬盘空间,超过100G数据不建议使用该方法
具体代码
一、GDAL镶嵌命令
gdal_merge [--help] [--help-general]
[-o <out_filename>] [-of <out_format>] [-co <NAME>=<VALUE>]...
[-ps <pixelsize_x> <pixelsize_y>] [-tap] [-separate] [-q] [-v] [-pct]
[-ul_lr <ulx> <uly> <lrx> <lry>] [-init "<value>[ <value>]..."]
[-n <nodata_value>] [-a_nodata <output_nodata_value>]
[-ot <datatype>] [-createonly] <input_file> [<input_file>]...
仅列举常用参数,具体参数参阅https://gdal.org/programs/gdal_merge.html
-o <out_filename>
输出文件的名称,如果尚不存在则将创建该文件(默认为“out.tif”)。
-of <format>
选择输出格式。从 GDAL 2.3 开始,如果未指定,则从扩展名猜测格式(之前为 GTiff)。使用短格式名称。
-ot <type>
强制输出图像带具有驱动程序支持的特定数据类型,可能是以下之一: Byte 、 Int8 、 UInt16 、 Int16 、 UInt32 、 Int32 、 UInt64 、 Int64 、 Float32 、 Float64 、 CInt32 、 CFloat32 或 CFloat64 。
-ps <pixelsize_x> <pixelsize_y>
用于输出文件的像素大小。如果未指定,将使用第一个输入文件的分辨率。
-ul_lr <ulx> <uly> <lrx> <lry>
输出文件的范围。如果未指定,将使用所有输入文件的聚合范围。
-v
完成镶嵌操作后生成详细输出。
-n <nodata_value>
忽略与此像素值合并的文件中的像素。
-a_nodata <output_nodata_value>
Assign a specified nodata value to output bands.
将指定的无数据值分配给输出波段。
-init <"value(s)">
使用这些值预初始化输出图像带。但是,它未在输出文件中标记为 nodata 值。如果仅给出一个值,则在所有频段中使用相同的值。(常用于生成单值图像)
<input_file> [<input_file>]
输入图像路径,可以为列表格式
二、代码示例
# 引用函数库
from osgeo import gdal
import os
import glob
inputdir = r'输入文件的文件夹地址'
outputdir = r'输出文件的文件地址,需要加上文件名.tif'
files = glob.glob(os.path.join(inputdir,'*.tif'))
paths = []
for file in files:
filepath = os.path.join(inputdir,file)
paths.append(filepath)
paths = ' '.join(paths)
print(paths)
command = f"python gdal_merge.py -o {outputdir} {paths}"
os.system(command)
三、注意事项
- 使用本方法主要因为可以自动化执行,如果遇到单幅图像ENVI打开速度很慢的情况,请先到GEE中验证时间筛选范围以及其他限制条件,同时由于均值镶嵌易产生位数较多的数据,所以中值镶嵌较均值镶嵌占用小,且不易出现ENVI无法打开的情况。
- 该代码仅为gdal中的一种镶嵌方法,且由于使用命令行命令进行镶嵌,所以无法进行多线程运算,速度较慢,占用空间较慢,如果条件不足请勿选择本方法。
- 可通过command字符串自行修改命令中的相关参数。