#!/bin/bash
# usage:
#   pth2csv < line.pth > line.csv
#
# convert MagicMaps (Gauss-Krger DHDN Meridian 9) to CSV

if [ $# -ne 0 ]; then
  echo "usage: pth2csv < line.pth > line.csv"
  exit
fi


# DHDN / Gauss-Kruger zone 2
#<31466> +proj=tmerc +lat_0=0 +lon_0=6 +k=1.000000 +x_0=2500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs  <>
# DHDN / Gauss-Kruger zone 3
#<31467> +proj=tmerc +lat_0=0 +lon_0=9 +k=1.000000 +x_0=3500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs  <>
# DHDN / Gauss-Kruger zone 4
#<31468> +proj=tmerc +lat_0=0 +lon_0=12 +k=1.000000 +x_0=4500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs  <>
# DHDN / Gauss-Kruger zone 5
#<31469> +proj=tmerc +lat_0=0 +lon_0=15 +k=1.000000 +x_0=5500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs  <>

epsgcode=31467

grep "^[0-9]" | cs2cs +proj=tmerc +init=epsg:$epsgcode +to +proj=longlat +datum=WGS84 -f %.9f -s | ascii2csv
