#!/bin/bash
# usage:
#   jgw2kal-utm mapname.jpg utmzone
#
# convert JPG UTM WorldFile calibration file to Glopus KAL

if [ $# -ne 2 ]; then
  echo 'Syntax: jgw2kal-utm mapname.jpg utmzone'
  exit
fi

name=$1
base=${name%.*}
utmzone=$2

image=$name
jgw=$base.jgw
kal=$base.kal

tempcoo=$base.coordinates

# get image size
sizex=`gdalinfo "$image" | grep "^Size is " | sed "s/Size is \(.*\), \(.*\)/\1/g"`
sizey=`gdalinfo "$image" | grep "^Size is " | sed "s/Size is \(.*\), \(.*\)/\2/g"`

world2pt $sizex $sizey < $jgw | cut -s -f 1-2 -d " " | \
invproj +proj=utm +zone=$utmzone +to +proj=longlat +datum=WGS84 -f %.9f > $tempcoo

world2pt $sizex $sizey < $jgw | cut -s -f 3-4 -d " " | \
paste $tempcoo - | asc2kal "$image" $sizex $sizey > $kal

rm $tempcoo
