#!/bin/bash

# Wikipedia sagt uebrigens, Dortmund sei bei 32U 393506 5708090
# die Schrittweite da unten sind _meter_

schrittweite=1000
unten=5690000
oben=5720000
links=380000
rechts=410000

cat <<EOF
<?xml version="1.0"?>
<gpx version="1.0" creator="Viking -- http://viking.sf.net/"
xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
xmlns="http://www.topografix.com/GPX/1/0"
xsi:schemaLocation="http://www.topografix.com/GPX/1/0
http://www.topografix.com/GPX/1/0/gpx.xsd">
<trk>
  <name>UTMGitter</name>
EOF


y=$unten
while [ $y -le $oben ] ;do
 echo "<trkseg>"
 x=$links
 while [ $x -le  $rechts ] ;do
  echo "$x $y" | cs2cs +proj=utm +zone=32 +to +proj=lonlat -f "%2f" | \
   cut -d" " -f1 | sed "s/^/<trkpt lon=\"/g;s/\t/\" lat=\"/g;;s/$/\">/g;"
  echo "</trkpt>"
  x=$(($x + $schrittweite))
 done
 echo "</trkseg>"
 y=$(($y + $schrittweite))
done

x=$links
while [ $x -le $rechts ] ;do
 echo "<trkseg>"
 y=$unten
 while [ $y -le  $oben ] ;do
  echo "$x $y" | cs2cs +proj=utm +zone=32 +to +proj=lonlat -f "%2f" | \
   cut -d" " -f1 | sed "s/^/<trkpt lon=\"/g;s/\t/\" lat=\"/g;;s/$/\">/g;"
  echo "</trkpt>"
  y=$(($y + $schrittweite))
 done
 echo "</trkseg>"
 x=$(($x + $schrittweite))
done

echo "</trk>"
echo "</gpx>"