#!/bin/bash
# usage:
#   kal2geotiffw map.kal newmap.tif "proj projection info"
# example
#   kal2geotiffw de0805ak.kal neu.tif "+proj=utm +zone=32 +datum=WGS84"
# projection info is optional
# uses KAL-generated WorldFile information for GeoTIFF calibration info

if [ $# -ne 2 ]; then
  if [ $# -ne 3 ]; then
    echo "kal2geotiffw map.kal newmap.tif \"proj projection info\""
    exit
  fi
fi

hash gdal_translate

kalfile=$1
outfile=$2
projection=$3

mapfile=`kal2mapname < "$kalfile"`
mapfilebase=${mapfile%.*}
worldfile=$mapfilebase.wld

if [ ! -f "$worldfile" ]; then
  kal2world < "$kalfile" > "$worldfile"
fi

if [ $# -eq 2 ]; then
 gdal_translate "$mapfile" -co WORLDFILE=YES "$outfile"
else
 gdal_translate "$mapfile" -a_srs "$projection" -co WORLDFILE=YES "$outfile"
fi
