GDAL分幅数据切瓦片。不拼接直接切

gdal提供很多现成的工具比如gdal2tiles是可以生成瓦片。但是在处理分幅tif的时候需要修改一些代码达到分幅切片的目的。

首先添加一个函数

def combine(former_tile, dstile): for band in range(dstile.RasterCount): array = numpy.max([former_tile.GetRasterBand(band + 1).ReadAsArray(), dstile.GetRasterBand(band + 1).ReadAsArray()], axis=0) dstile.GetRasterBand(band + 1).WriteArray(array)

下一步找到

CreateCopy 关键词 入下图修改。

if os.path.exists(tilefilename): combine(gdal.Open(tilefilename), dstile) out_driver.CreateCopy(tilefilename, dstile, strict=0)

 

python gdal2tiles.py -v -z 1-13 --resampling=average G:\land\fromglc10v01_0_8.tif G:\tiles