To visualize SFSM outputs, we should generate the outputs in a format suitable for web rendering.
warped = gdal.Warp(
"",
src_ds,
format="MEM",
dstSRS="EPSG:3857",
)
gdal.Translate(
output_path,
warped,
format="COG",
creationOptions=[
"COMPRESS=LZW",
"OVERVIEWS=NONE",
],
)