src/gpx.cc

   1 /* gpx.cc
   2  * Export from Aven as GPX.
   3  */
   4 /* Copyright (C) 2012 Olaf Kähler
   5  * Copyright (C) 2012,2013,2014,2015,2016 Olly Betts
   6  *
   7  * This program is free software; you can redistribute it and/or modify
   8  * it under the terms of the GNU General Public License as published by
   9  * the Free Software Foundation; either version 2 of the License, or
  10  * (at your option) any later version.
  11  *
  12  * This program is distributed in the hope that it will be useful,
  13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15  * GNU General Public License for more details.
  16  *
  17  * You should have received a copy of the GNU General Public License
  18  * along with this program; if not, write to the Free Software
  19  * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301 USA
  20  */
  21
  22 #ifdef HAVE_CONFIG_H
  23 # include <config.h>
  24 #endif
  25
  26 #include "gpx.h"
  27
  28 #include "export.h" // For LABELS, etc
  29
  30 #include <stdio.h>
  31 #include <string>
  32 #include <time.h>
  33 #include <math.h>
  34
  35 #include "useful.h"
  36 #include <proj_api.h>
  37
  38 #include "aven.h"
  39 #include "message.h"
  40
  41 using namespace std;
  42
  43 #define WGS84_DATUM_STRING "+proj=longlat +ellps=WGS84 +datum=WGS84"
  44
  45 static void
  46 html_escape(FILE *fh, const char *s)
  47 {
  48     while (*s) {
  49         switch (*s) {
  50             case '<':
  51                 fputs("&lt;", fh);
  52                 break;
  53             case '>':
  54                 fputs("&gt;", fh);
  55                 break;
  56             case '&':
  57                 fputs("&amp;", fh);
  58                 break;
  59             default:
  60                 PUTC(*s, fh);
  61         }
  62         ++s;
  63     }
  64 }
  65
  66 GPX::GPX(const char * input_datum)
  67     : pj_input(NULL), pj_output(NULL), in_trkseg(false), trk_name(NULL)
  68 {
  69     if (!(pj_input = pj_init_plus(input_datum))) {
  70         wxString m = wmsg(/*Failed to initialise input coordinate system “%s”*/287);
  71         m = wxString::Format(m.c_str(), input_datum);
  72         throw m;
  73     }
  74     if (!(pj_output = pj_init_plus(WGS84_DATUM_STRING))) {
  75         wxString m = wmsg(/*Failed to initialise output coordinate system “%s”*/288);
  76         m = wxString::Format(m.c_str(), WGS84_DATUM_STRING);
  77         throw m;
  78     }
  79 }
  80
  81 GPX::~GPX()
  82 {
  83     if (pj_input)
  84         pj_free(pj_input);
  85     if (pj_output)
  86         pj_free(pj_output);
  87     free((void*)trk_name);
  88 }
  89
  90 const int *
  91 GPX::passes() const
  92 {
  93     static const int default_passes[] = { LABELS|ENTS|FIXES|EXPORTS, LEGS|SURF, 0 };
  94     return default_passes;
  95 }
  96
  97 /* Initialise GPX routines. */
  98 void GPX::header(const char * title, const char *, time_t datestamp_numeric,
  99                  double, double, double, double, double, double)
 100 {
 101     fputs(
 102 "<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n"
 103 "<gpx version=\"1.0\" creator=\"" PACKAGE_STRING " (aven) - https://survex.com/\""
 104 " xmlns:xsi=\"http://www.w3.org/2001/XMLSchema-instance\""
 105 " xmlns=\"http://www.topografix.com/GPX/1/0\""
 106 " xsi:schemaLocation=\"http://www.topografix.com/GPX/1/0"
 107 " http://www.topografix.com/GPX/1/0/gpx.xsd\">\n", fh);
 108     if (title) {
 109         fputs("<name>", fh);
 110         html_escape(fh, title);
 111         fputs("</name>\n", fh);
 112         trk_name = strdup(title);
 113     }
 114     if (datestamp_numeric != time_t(-1)) {
 115         struct tm * tm = gmtime(&datestamp_numeric);
 116         if (tm) {
 117             char buf[32];
 118             if (strftime(buf, sizeof(buf), "%Y-%m-%dT%H:%M:%SZ", tm)) {
 119                 fputs("<time>", fh);
 120                 fputs(buf, fh);
 121                 fputs("</time>\n", fh);
 122             }
 123         }
 124     }
 125     // FIXME: optional in GPX, but perhaps useful:
 126     // <bounds minlat="..." minlon="..." maxlat="..." maxlon="..." />
 127     // NB Not necessarily the same as the bounds in survex coords translated
 128     // to WGS84 lat+long...
 129 }
 130
 131 void
 132 GPX::line(const img_point *p1, const img_point *p, unsigned /*flags*/, bool fPendingMove)
 133 {
 134     if (fPendingMove) {
 135         if (in_trkseg) {
 136             fputs("</trkseg><trkseg>\n", fh);
 137         } else {
 138             fputs("<trk>", fh);
 139             if (trk_name) {
 140                 fputs("<name>", fh);
 141                 html_escape(fh, trk_name);
 142                 fputs("</name>", fh);
 143             }
 144             fputs("<trkseg>\n", fh);
 145             in_trkseg = true;
 146         }
 147         double X = p1->x, Y = p1->y, Z = p1->z;
 148         pj_transform(pj_input, pj_output, 1, 1, &X, &Y, &Z);
 149         X = deg(X);
 150         Y = deg(Y);
 151         // %.8f is at worst just over 1mm.
 152         fprintf(fh, "<trkpt lon=\"%.8f\" lat=\"%.8f\"><ele>%.2f</ele></trkpt>\n", X, Y, Z);
 153     }
 154     double X = p->x, Y = p->y, Z = p->z;
 155     pj_transform(pj_input, pj_output, 1, 1, &X, &Y, &Z);
 156     X = deg(X);
 157     Y = deg(Y);
 158     // %.8f is at worst just over 1mm.
 159     fprintf(fh, "<trkpt lon=\"%.8f\" lat=\"%.8f\"><ele>%.2f</ele></trkpt>\n", X, Y, Z);
 160 }
 161
 162 void
 163 GPX::label(const img_point *p, const char *s, bool /*fSurface*/, int type)
 164 {
 165     double X = p->x, Y = p->y, Z = p->z;
 166     pj_transform(pj_input, pj_output, 1, 1, &X, &Y, &Z);
 167     X = deg(X);
 168     Y = deg(Y);
 169     // %.8f is at worst just over 1mm.
 170     fprintf(fh, "<wpt lon=\"%.8f\" lat=\"%.8f\"><ele>%.2f</ele><name>", X, Y, Z);
 171     html_escape(fh, s);
 172     fputs("</name>", fh);
 173     // Add a "pin" symbol with colour matching what aven shows.
 174     switch (type) {
 175         case FIXES:
 176             fputs("<sym>Pin, Red</sym>", fh);
 177             break;
 178         case EXPORTS:
 179             fputs("<sym>Pin, Blue</sym>", fh);
 180             break;
 181         case ENTS:
 182             fputs("<sym>Pin, Green</sym>", fh);
 183             break;
 184     }
 185     fputs("</wpt>\n", fh);
 186 }
 187
 188 void
 189 GPX::footer()
 190 {
 191     if (in_trkseg)
 192         fputs("</trkseg></trk>\n", fh);
 193     fputs("</gpx>\n", fh);
 194 }