Hi,I have to retile about 100,000 images of size about 2667x2667 pixels for geoserver. The source images' nodata area were badly damaged due to compression so I decided to use VRT pixel function max to make the result tiles seamless.However, the output was quite weird -- Killed appeared and the scr
Thank you, the following code worked perfectly:
import os
import osgeo.gdal
name, ext = os.path.splitext(os.path.basename(in_path))
srcDS = osgeo.gdal.OpenEx(in_path)
ds = osgeo.gdal.VectorTranslate(out_path, srcDS, SQLStatement="SELECT
ST_Union(geometry) FROM " + name , SQLDialect='sqlite')
Che
