gdalbuildvrt -of GTIFF -srcnodata 0 input.tif myVrt.vrt

# convert to EPSG 4326
gdalwarp -s_srs EPSG:26913 -t_srs EPSG:4326 -r bilinear -ot Byte -of GTiff -co COMPRESS=JPEG -co JPEG_QUALITY=75 bac_area_1941.img baca2.tif

# make pyramided  overview of tif file. speeds up the rio_tiler which is kind of a hunk of shit. 
gdaladdo -r average  baca1941/baca4.tif 2 4 8 16
# or
gdaladdo -r average --config COMPRESS_OVERVIEW LZW --config BIGTIFF_OVERVIEW IF_NEEDED baca1941/baca4.tif 2 4 8 16
# apply color map
gdaldem color-relief baca1941/baca3.tif color.txt  -co COMPRESS=LZW -alpha baca1941/baca4.tif
# color map for the above command
256    256 256 256 
255    255 255 255
128    128 128 128
2      2   2   2
1      0   0   0   0  
0      0   0   0   0
nv     0   0   0   0