#!/bin/bash

input=dem_eu_30m.tif
output=out/subwindow.tif

gdal_calc.py -A $input --A_band=1 --outfile=hi_byte.tif --calc="(A>>8)&255" --type=Byte
gdal_calc.py -A $input --A_band=1 --outfile=lo_byte.tif --calc="A&255" --type=Byte
gdal_merge.py -separate -o output_multi_band.tif hi_byte.tif lo_byte.tif



