#!/bin/bash

REG="-R-22.75/-21.25/63.5/64"
PROJ="-Jl-22/64/63/64/1:100000"

# from Cartesian to geographical
# getting geographical coordinates of the Cartesian origin
echo From Cartesian to geographic
echo [0 0] "(origin of the Cartesian system)"
echo "0 0" | mapproject -I -C -F ${REG} ${PROJ}
# getting geographical coordinates of the point x=10km, y=-25km
echo [10000 -25000] "(point [10 km -25 km])"
echo "10000 -25000" | mapproject -I -C -F ${REG} ${PROJ}


# from geographical to Cartesian
# getting the Cartesian origin
echo From geographic to Cartesian
echo [-22 64] "(origin of the Cartesian system)"
echo "-22 64" | mapproject -C -F ${REG} ${PROJ}
# getting the Cartesian coordinates of the point -21W 63N
echo [-21 63]
echo "-21 63" | mapproject -C -F ${REG} ${PROJ}