#!/bin/sh

#echo executed ${0} {$1} ${2} ${3} ${4}
#echo path $PATH

SIZE=800 # szerokosc i wysokosc obrazka
DPI=120 # rozdzielczosc w DITA - nie ruszac
SCALE=500 # domyslna skala
BORDER=5 # liczba metrow, o ktora nalezy powiekszyc zasieg z kazdej strony wzgledem danych wejsciowych
BORDER_CUT="5"

IMAGE_FILE="${PWD}/image" # nazwa plikow obrazow wraz ze sciezka - sam prefix (zostana dodane oznaczenia wspolrzednych kafelkow oraz dodane rozszerzenie)
XML_FILE="${PWD}/images.xml" # nazwa wynikowego pliku XML wraz ze sciezka
FONT_DIR="${PWD}/fonts"

MAP_SCHEMA_DIR="./map_schema"
MAP_SCHEMA_HEAD="${MAP_SCHEMA_DIR}/head.map.inc"
MAP_SCHEMA_TAIL="${MAP_SCHEMA_DIR}/tail.map.inc"
MAP_SCHEMA_GML="${MAP_SCHEMA_DIR}/gml_layer.map.inc"
MAP_SCHEMA_WMS="${MAP_SCHEMA_DIR}/wms_layer.map.inc"


error() {
	echo '<?xml version="1.0" encoding="UTF-8"?>'
	echo "<Error>${1}</Error>"
	exit 1
}
usage() {
	cat >&2 << EOT
Usage: ${0} -g|--gml FILE [OPTION]...

 FILE - GML input file

Options:
 -s, --scale SCALE		Scale, defult 500 (means 1:500)
 -b, --border METERS		Size of border in meters that will be added
				  around GML data (default 5)
 -cb, --cut-border METERS	Size of additional border in meters
				  while splitting image (default 5)
 -w, --wms WMS_LAYER[s]		URLs and layer names (as following arguments)
				  of external WMS layers (URL|LAYER_NAME)

Example:

 ${0} -g gml/projekt_3625.gml -s 500 -w "https://biuro.biall-net.pl/wms/rastry|Przywidz_Jesionowa" \\
    "https://biuro.biall-net.pl/wms/podklady|MDCP_Przywidz_osiedle_jesionowa_lin" \\
    "https://biuro.biall-net.pl/wms/podklady|MDCP_Przywidz_osiedle_jesionowa_textx"
EOT
	exit 1
}

declare -a WMS=()
while [ "${1}" ]; do
	case ${1} in
		"-s"|"--scale")
			PARAM="scale"
			;;
		"-b"|"--border")
			PARAM="border"
			;;
		"-cb"|"--cut-border")
			PARAM="cut_border"
			;;
		"-g"|"--gml")
			PARAM="gml"
			;;
		"-w"|"--wms")
			PARAM="wms"
			;;
		*)
			case ${PARAM} in
				"scale")
					SCALE=${1}
					;;
				"border")
					BORDER=${1}
					;;
				"cut_border")
					CUT_BORDER={1}
					;;
				"gml")
					FILE=${1}
					;;
				"wms")
					WMS+=("${1}")
					;;
				*)
					usage
					;;
			esac
	esac
	shift
done

if [ -z "${FILE}" ]; then
	usage
fi

if [ ! -e ${FILE} ]; then
	error "File ${FILE} not found."
fi


TMP=`mktemp -d`

size=`echo "${SIZE}/${DPI}*2.54*${SCALE}/100" | bc -l`
size_cut=`echo "${size} - ${BORDER_CUT}" | bc -l`

ogr2ogr -t_srs EPSG:2177 ${TMP}/2177.shp ${FILE} 2>/dev/null
extent=(`ogrinfo -al -so ${TMP}/2177.shp | grep Extent | sed 's/[()]/|/g;s/, /|/g' | awk -F '|' '{print $2,$3,$5,$6}'`)
xmin=${extent[0]}
ymin=${extent[1]}
xmax=${extent[2]}
ymax=${extent[3]}
echo "Project extent: ${xmin} ${ymin} ${xmax} ${ymax}"

xmin=`echo "${xmin} - 5" | bc -l`
ymin=`echo "${ymin} - 5" | bc -l`
xmax=`echo "${xmax} + 5" | bc -l`
ymax=`echo "${ymax} + 5" | bc -l`
echo "New extent after 5m border added: ${xmin} ${ymin} ${xmax} ${ymax}"

lx=`echo "(${xmax} - ${xmin} - ${BORDER_CUT}) / ${size_cut} + ((${xmax} - ${xmin} - ${BORDER_CUT}) % ${size_cut} > 0)" | bc`
ly=`echo "(${ymax} - ${ymin} - ${BORDER_CUT}) / ${size_cut} + ((${ymax} - ${ymin} - ${BORDER_CUT}) % ${size_cut} > 0)" | bc`
echo "Images count (width x height): ${lx}x${ly}"

xmin=`echo "${xmin} - (${lx} * ${size_cut} + ${BORDER_CUT} - (${xmax} - ${xmin})) / 2" | bc -l`
ymin=`echo "${ymin} - (${ly} * ${size_cut} + ${BORDER_CUT} - (${ymax} - ${ymin})) / 2" | bc -l`
xmax=`echo "${xmin} + ${lx} * ${size_cut} + ${BORDER_CUT}" | bc -l`
ymax=`echo "${ymin} + ${ly} * ${size_cut} + ${BORDER_CUT}" | bc -l`
echo "New extent after centering project features: ${xmin} ${ymin} ${xmax} ${ymax}"

MAP="${TMP}/map.map"
EXTENT="${xmin} ${ymin} ${xmax} ${ymax}"

sed "s/__EXTENT__/${EXTENT}/g;s|__FONT_DIR__|${FONT_DIR}|g" ${MAP_SCHEMA_HEAD} > ${MAP}

if [ ${#WMS[@]} -gt 0 ]; then
	for i in `seq 0 $[${#WMS[@]}-1]`; do
		URL=`echo ${WMS[i]} | cut -d '|' -f 1`
		LAYER_NAME=`echo ${WMS[i]} | cut -d '|' -f 2`
		sed "s|__URL__|${URL}|g;s|__EXTENT__|${EXTENT}|g;s|__LAYER_NAME__|${LAYER_NAME}|g" ${MAP_SCHEMA_WMS} >> ${MAP}
	done
fi

sed "s|__GML_FILE__|${FILE}|g" ${MAP_SCHEMA_GML} >> ${MAP}
cat ${MAP_SCHEMA_TAIL} >> ${MAP}

echo '<?xml version="1.0" encoding="UTF-8"?>
<images>' > ${XML_FILE}
for y in `seq 1 ${ly}`; do
	for x in `seq 1 ${lx}`; do
		printf "Generating ${x}x${y} "
		lxmin=`echo "${xmin} + ${size_cut} * (${x} -1)" | bc -l`
		lymin=`echo "${ymin} + ${size_cut} * (${y} -1)" | bc -l`
		lxmax=`echo "${lxmin} + ${size}" | bc -l`
		lymax=`echo "${lymin} + ${size}" | bc -l`
		printf "(BBOX: ${lxmin} ${lymin} ${lxmax} ${lymax})..."
		file="${IMAGE_FILE}-${y}-${x}.png"
		count=`ogrinfo -so -al -spat ${lxmin} ${lymin} ${lxmax} ${lymax} ${TMP}/2177.shp | grep "Feature Count:" | sed 's/.*: //g'`
		if [ ${count} -gt 0 ]; then
			echo "<image x=\"${x}\" y=\"${y}\">" >> ${XML_FILE}
			echo "  <xmin>${lxmin}</xmin>" >> ${XML_FILE}
			echo "  <ymin>${lymin}</ymin>" >> ${XML_FILE}
			echo "  <xmax>${lxmax}</xmax>" >> ${XML_FILE}
			echo "  <ymax>${lymax}</ymax>" >> ${XML_FILE}
			echo "  <file>${file}</file>" >> ${XML_FILE}
			echo "</image>" >> ${XML_FILE}
			shp2img -m "${MAP}" -o ${file} -e ${lxmin} ${lymin} ${lxmax} ${lymax}
			echo " done."
		else
			echo " no features found - skipping."
		fi
	done
done
echo '</images>' >> ${XML_FILE}
echo "===== XML (${XML_FILE}) ====="
cat ${XML_FILE}

rm -rf ${TMP}