#!/bin/bash

touch -d "4 days ago" /tmp/4dayspeaks

if [ /home/www/geo/topo/topographic_isolation_viefinderpanoramas.txt -ot /tmp/4dayspeaks ] ; then
 rm -f /tmp/overpasssaettel.osm /tmp/peakquery.ql
 echo '[out:xml][timeout:5000];(node["natural"="volcano"];node["natural"="peak"]);out body;>;out meta qt;' > /tmp/peakquery.ql
 wget -O /tmp/overpassgipfel.osm --post-file=/tmp/peakquery.ql "http://overpass-api.de/api/interpreter"
 
 n=`wc -l  /tmp/overpassgipfel.osm | cut -d" " -f1`
 d=`date`
 
 if [ $n -gt 400000 ] ; then
  (
   echo "# Dominanz von Gipfeln nach Hoehendaten (Max. 100km)"
   echo "#"
   echo '# erzeugt mit "hole (node["natural"="peak"];node["natural"="volcano"];)'
   echo "#              von overpass-api.de und schiebe das Ergebnis durch osmconvert und isolation -f hoehendaten.tiff\""
   echo "#"
   echo "#"
   echo "# Programm dazu: https://github.com/der-stefan/OpenTopoMap/blob/master/mapnik/tools/isolation.c"
   echo "# Dokumentation: https://wiki.openstreetmap.org/wiki/User:Maxbe/Dominanz_von_Gipfeln"
   echo "#"
   echo "# Hoehendaten: 3-Bogensekunden-DEM von http://viewfinderpanoramas.org/ mit freundlicher Genehmigung von Jonathan de Ferranti"
   echo "# ID und Koordinaten von OSM Mitwirkenden http://www.openstreetmap.org (Lizenz: ODBL)"
   echo "#"
   echo "# Stand: $d"
   echo "#"
   echo "# ID;LON;LAT;Dominanz"
   osmconvert /tmp/overpassgipfel.osm --out-csv --csv-separator=";" --csv="@id @lon @lat ele" | /home/www/bin/isolation -r 100000 -f /hoehen/viewfinder/viewfinderDEM.tiff | sort -n 
  ) > /tmp/peakisolation.txt
  chmod 0644 /tmp/peakisolation.txt
  cp /tmp/peakisolation.txt /home/www/geo/topo/topographic_isolation_viefinderpanoramas.txt
  chmod 0644 /home/www/geo/topo/topographic_isolation_viefinderpanoramas.txt 
  (
   echo "# Vergleich von ele und Satellitendaten"
   echo "#"
   echo "# Hoehendaten: 3-Bogensekunden-DEM von http://viewfinderpanoramas.org/ mit freundlicher Genehmigung von Jonathan de Ferranti"
   echo "# ID und Koordinaten von OSM Mitwirkenden http://www.openstreetmap.org (Lizenz: ODBL)"
   echo "#"
   echo "# Stand: $d"
   echo "#"
   echo "# ID;LON;LAT;OSM-ele;SRTM-ele"
   osmconvert /tmp/overpassgipfel.osm --out-csv --csv-separator=";" --csv="@id @lon @lat ele" | /home/www/bin/peakele -f /hoehen/viewfinder/viewfinderDEM.tiff | sort -n 
  ) > /tmp/peakele.txt
  chmod 0644 /tmp/peakele.txt
  cp /tmp/peakele.txt /home/www/geo/topo/peakele_viewfinderpanoramas.txt
  chmod 0644 /home/www/geo/topo/peakele_viewfinderpanoramas.txt
 fi
fi