Make proper GeoTiff COG with overviews through GDAL

make_cog_from_tiff.sh

#!/bin/bash

# https://trac.osgeo.org/gdal/wiki/CloudOptimizedGeoTIFF

input=$1
output=$2
oversample="nearest"  # or nearest, average

gdal_translate $1 interm_step1.tif -co TILED=YES -co COMPRESS=DEFLATE

gdaladdo -r $oversample interm_step1.tif 2 4 8 16 32

gdal_translate interm_step1.tif $2 -co TILED=YES -co COMPRESS=LZW -co COPY_SRC_OVERVIEWS=YES \
      -co BIGTIFF=IF_SAFER  -co BLOCKXSIZE=512 -co BLOCKYSIZE=512 --config GDAL_TIFF_OVR_BLOCKSIZE 512

rm interm_step1.tif