gpst

   1 #!/usr/bin/perl -w
   2
   3 #=======================================================================
   4 # gpst
   5 # File ID: 082106ec-f924-11dd-b757-0001805bf4b1
   6 # Converts between various GPS formats
   7 #
   8 # Character set: UTF-8
   9 # ©opyleft 2002– Øyvind A. Holm <sunny@sunbase.org>
  10 # License: GNU General Public License version 3 or later, see end of
  11 # file for legal stuff.
  12 #=======================================================================
  13
  14 use strict;
  15 use Getopt::Long;
  16 use Time::Local qw { timegm_nocheck };
  17
  18 BEGIN {
  19     push(@INC, "$ENV{'HOME'}/bin/src/gpstools");
  20 }
  21
  22 use GPST;
  23 use GPSTdate;
  24 use GPSTdebug;
  25 use GPSTgeo;
  26 use GPSTxml;
  27
  28 $| = 1;
  29
  30 our $Debug = 0;
  31
  32 our %Opt = (
  33     # Initial values for command line arguments {{{
  34
  35     'chronology' => 0,
  36     'create-breaks' => 0,
  37     'debug' => 0,
  38     'double-y-scale' => 0,
  39     'epoch' => 0,
  40     'fix' => 0,
  41     'from-date' => "",
  42     'help' => 0,
  43     'inside' => 0,
  44     'output-format' => "gpsml",
  45     'outside' => 0,
  46     'pos1' => "",
  47     'pos2' => "",
  48     'require' => "",
  49     'round' => "",
  50     'save-to-file' => "\n", # \n =  undefined, it’s banned in filenames anyway.
  51     'short-date' => 0,
  52     'skip-dups' => 0,
  53     'strip-whitespace' => 0,
  54     'time-shift' => 0,
  55     'undefined' => "",
  56     'verbose' => 0,
  57     'version' => 0,
  58
  59     # }}}
  60 );
  61
  62 our $progname = $0;
  63 $progname =~ s/^.*\/(.*?)$/$1/;
  64 our $VERSION = "0.00";
  65
  66 Getopt::Long::Configure("bundling");
  67 GetOptions(
  68     # Command line options {{{
  69
  70     "chronology" => \$Opt{'chronology'},
  71     "create-breaks|t" => \$Opt{'create-breaks'},
  72     "debug" => \$Opt{'debug'},
  73     "double-y-scale|y" => \$Opt{'double-y-scale'},
  74     "epoch|e" => \$Opt{'epoch'},
  75     "fix" => \$Opt{'fix'},
  76     "from-date=s" => \$Opt{'from-date'},
  77     "help|h" => \$Opt{'help'},
  78     "inside" => \$Opt{'inside'},
  79     "output-format|o=s" => \$Opt{'output-format'},
  80     "outside" => \$Opt{'outside'},
  81     "pos1=s" => \$Opt{'pos1'},
  82     "pos2=s" => \$Opt{'pos2'},
  83     "require|r=s" => \$Opt{'require'},
  84     "round|R=s" => \$Opt{'round'},
  85     "save-to-file|S=s" => \$Opt{'save-to-file'},
  86     "short-date|s" => \$Opt{'short-date'},
  87     "skip-dups|d" => \$Opt{'skip-dups'},
  88     "strip-whitespace|w" => \$Opt{'strip-whitespace'},
  89     "time-shift|T=i" => \$Opt{'time-shift'},
  90     "undefined|n=s" => \$Opt{'undefined'},
  91     "verbose|v+" => \$Opt{'verbose'},
  92     "version" => \$Opt{'version'},
  93
  94     # }}}
  95 ) || die("$progname: Option error. Use -h for help.\n");
  96
  97 my %Dat;
  98
  99 my $PAUSE_LIMIT = 2 * 60; # Antall sekunder mellom to punkter det må til før en move legges inn.
 100 my $Udef = "?";
 101 my $DIGIT = '[0-9\.\-\+]'; # Used in regexps
 102 $GPST::Spc = $Opt{'strip-whitespace'} ? "" : " ";
 103 my $Spc = $GPST::Spc; # FIXME
 104 my $found_move = 0; # Settes til 1 hvis en /^# move$/ blir funnet.
 105 my $first_time = 0;
 106 my $last_time = 0;
 107 my ($last_lon, $last_lat, $last_line) =
 108    (     1000,      1000,         ""); # Vi kan jo teoretisk sett være i Greenwich eller på ekvator
 109 my ($lat1, $lon1, $lat2, $lon2) =
 110    (-1000, -1000,  1000,  1000);
 111
 112 our %Cmd = (
 113     'gpsbabel' => '/usr/local/bin/gpsbabel',
 114 );
 115
 116 my %Poscount = ();
 117
 118 if ($Opt{'output-format'} =~ /^(gpx|pgtab)$/) {
 119     $Opt{'require'} .= "p";
 120 }
 121 my %Req = (
 122     'ele' => ($Opt{'require'} =~ /e/) ? 1 : 0,
 123     'position' => ($Opt{'require'} =~ /p/) ? 1 : 0,
 124     'time' => ($Opt{'require'} =~ /t/) ? 1 : 0,
 125 );
 126 $Opt{'require'} =~ /[^ept]/
 127     && die("$0: Unknown flag in --require (-r) value\n");
 128
 129 $Opt{'debug'} && ($Debug = 1);
 130 $Opt{'help'} && usage(0);
 131 if ($Opt{'version'}) {
 132     print_version();
 133     exit(0);
 134 }
 135
 136 if ($Opt{'pos1'} =~ /^($DIGIT+),($DIGIT+)$/) {
 137     $lat1 = $1;
 138     $lon1 = $2;
 139 }
 140 if ($Opt{'pos2'} =~ /^($DIGIT+),($DIGIT+)$/) {
 141     $lat2 = $1;
 142     $lon2 = $2;
 143 }
 144 if ($lat1 > $lat2) {
 145     my $Tmp = $lat1;
 146     $lat1 = $lat2;
 147     $lat2 = $Tmp;
 148 }
 149 if ($lon1 > $lon2) {
 150     my $Tmp = $lon1;
 151     $lon1 = $lon2;
 152     $lon2 = $Tmp;
 153 }
 154
 155 if ($Opt{'epoch'} && $Opt{'short-date'}) {
 156     die("$progname: Cannot mix the --epoch (-e) and --short-date (-s) options\n");
 157 }
 158
 159 if ($Opt{'inside'} && $Opt{'outside'}) {
 160     die("$progname: Cannot mix the --inside and --outside options\n");
 161 }
 162
 163 # To avoid printing out extra "/> at the start of svg output:
 164 my $svg_start_thing = "";
 165
 166 my %Round = ();
 167
 168 if (defined($Opt{'round'})) {
 169     my $R = $Opt{'round'};
 170     $R =~ s/([a-z]+)=(\d+)/($Round{$1}=$2, "")/eg;
 171 }
 172
 173 length($Opt{'undefined'}) && ($Udef = $Opt{'undefined'});
 174
 175 $Opt{'save-to-file'} eq "\n" && print_header(*STDOUT);
 176
 177 my @first_lines;
 178 my $xml_data;
 179 my $data_line = "";
 180 our $curr_file = "";
 181
 182 my $from_stdin = scalar(@ARGV) ? 0 : 1;
 183
 184 $from_stdin && push(@ARGV, "-");
 185
 186 for $curr_file (@ARGV) {
 187     # Scan through stdin or specified files and send every GPS entry to
 188     # print_entry()
 189     # {{{
 190     print(STDERR "$progname: Opening \"$curr_file\" for read\n") if $Opt{'verbose'};
 191     if (open(CurrFP, "<$curr_file")) {
 192         # {{{
 193         while (<CurrFP>) {
 194             $data_line = $_;
 195             %Dat = (
 196                 'year' => '', 'month' => '', 'day' => '',
 197                 'hour' => '', 'min' => '', 'sec' => '',
 198                 'epoch' => '',
 199                 'date-format' => '',
 200                 'lat' => '', 'lon' => '',
 201                 'ele' => '',
 202                 'desc' => '',
 203                 'error' => "",
 204                 'what' => 'tp',
 205             );
 206
 207             $Opt{'epoch'} && ($Dat{'date-format'} = "epoch");
 208             $Opt{'short-date'} && ($Dat{'date-format'} = "short");
 209
 210             if ($Opt{'save-to-file'} ne "\n") {
 211                 push(@first_lines, $_);
 212             }
 213             s/^# error // && ($Dat{'error'} = "error");
 214             s/^# ?// && ($Dat{'error'} = "desc");
 215             $xml_data = "";
 216             if (m#^<(e?tp)\b(.*?)>(.*?)</(e?tp)>\s*$#) {
 217                 # gpsml — The main storage format {{{
 218                 my ($Elem, $Props, $Data) =
 219                    (   $1,     $2,    $3);
 220                 my $err_str = ($Props =~ /\berr="(.*?)"/) ? $1 : "error";
 221                 $Elem eq "etp" && ($Dat{'error'} = $err_str);
 222                 my $Time = "";
 223                 $Data =~ m#<time>(.*?)</time># && ($Time = $1);
 224                 $Time =~ s{
 225                     (\d\d\d\d)-?(\d\d)-?(\d\d)[T ](\d\d):?(\d\d):?([\d\.]+?)Z
 226                 } {
 227                     ($Dat{'year'}, $Dat{'month'}, $Dat{'day'},
 228                      $Dat{'hour'},   $Dat{'min'}, $Dat{'sec'}) =
 229                     (          $1,            $2,          $3,
 230                                $4,            $5,          $6);
 231                     "";
 232                 }ex;
 233                 $Data =~ m#<lat>($DIGIT*?)</lat># && ($Dat{'lat'} = $1);
 234                 $Data =~ m#<lon>($DIGIT*?)</lon># && ($Dat{'lon'} = $1);
 235                 $Data =~ m#<ele>($DIGIT*?)</ele># && ($Dat{'ele'} = $1);
 236                 $Data =~ m#<desc>(.*?)</desc># && ($Dat{'desc'} = $1);
 237                 print_entry(%Dat);
 238                 # }}}
 239             } elsif (m#^<break\b.*?/>#) {
 240                 $found_move = 1;
 241             } elsif (m#^<(title|pause)\b.*?>(.*?)</(title|pause)>#) {
 242                 $Dat{'what'} = $1;
 243                 $Dat{$1} = $2;
 244                 print_entry(%Dat);
 245             } elsif (m#^<desc\b.*?>(.*$)#s) {
 246                 $Dat{'what'} = "desc";
 247                 my $Txt = $1;
 248                 until ($Txt =~ m#</desc>#s) {
 249                     $Txt .= <CurrFP>;
 250                 }
 251                 $Txt =~ s#^(.*)(</desc>.*$)#$1#s;
 252                 $Dat{'desc'} = $Txt;
 253                 print_entry(%Dat);
 254             } elsif (/<gpx\b/) {
 255                 $xml_data = $_;
 256                 $xml_data .= join("", <CurrFP>);
 257                 if (!length($Opt{'output-format'})) {
 258                     $Opt{'output-format'} = "gpx";
 259                     print_header(*STDOUT);
 260                 }
 261                 read_xmlfile($xml_data);
 262                 last;
 263             } elsif (/^move$/) {
 264                 $found_move = 1;
 265             } elsif (m#^(\d+)\t($DIGIT+)\t($DIGIT+)\t($DIGIT)#) {
 266                 # CSV format, epoch style {{{
 267                 my ($ep_time, $lon_val, $lat_val, $Alt) =
 268                    (      $1,       $2,       $3,   $4);
 269                 $Dat{'epoch'} = $ep_time;
 270                 ($Dat{'sec'}, $Dat{'min'}, $Dat{'hour'},
 271                  $Dat{'day'}, $Dat{'month'}, $Dat{'year'},
 272                  $Dat{'wday'}, $Dat{'yday'}) = gmtime($ep_time);
 273                 $Dat{'month'}++; # Urgh Ⅰ
 274                 $Dat{'year'} += 1900; # Urgh Ⅱ
 275                 print_entry(%Dat);
 276                 # }}}
 277             } elsif (
 278                 m{^
 279                     (\d\d\d\d)-?(\d\d)-?(\d\d)[T\ ](\d\d):?(\d\d):?(\d\d)Z?\t
 280                     ($DIGIT+)\t($DIGIT+)\t($DIGIT)
 281                 }x
 282             ) {
 283                 # CSV format, human-readable date format {{{
 284                 ($Dat{'year'}, $Dat{'month'}, $Dat{'day'},
 285                  $Dat{'hour'}, $Dat{'min'}, $Dat{'sec'},
 286                  $Dat{'lon'}, $Dat{'lat'}, $Dat{'ele'}) =
 287                 ($1, $2, $3,
 288                  $4, $5, $6,
 289                  $7, $8, $9);
 290                 print_entry(%Dat);
 291                 # }}}
 292             } elsif (/^Trackpoint\t/) {
 293                 # Trackpoint\tN60.41630 E5.31675\t09.02.2006 20:24:37 (UTC)\t13.6 m\t\t93.9 m\t00:00:06\t56 kph\t123° true {{{
 294                 #
 295                 # Trackpoint\t
 296                 # N60.41630 E5.31675\t
 297                 # 09.02.2006 20:24:37 (UTC)\t
 298                 # 13.6 m\t
 299                 # \t
 300                 # 93.9 m\t
 301                 # 00:00:06\t
 302                 # 56 kph\t
 303                 # 123° true
 304                 my $Orig = $_;
 305                 $Orig =~ s/[\r\n]+$//;
 306                 my ($Marker_f, $Position_f, $Time_f, $Alt_f, $Depth_f,
 307                     $Leglength_f, $Legtime_f, $Legspeed_f, $Legcourse_f) =
 308                        split(/\t/, $Orig .
 309                                    # Nødløsning for å unngå at variabler
 310                                    # blir udefinert.
 311                                    "\t\t\t\t\t\t\t\t\t\t"
 312                        );
 313                 # D(join("",
 314                 #     "Position_f=\"$Position_f\" \x7B\x7B\x7B\n",
 315                 #     "Time_f=\"$Time_f\"\n",
 316                 #     "Alt_f=\"$Alt_f\"\n",
 317                 #     "Depth_f=\"$Depth_f\"\n",
 318                 #     "Leglength_f=\"$Leglength_f\"\n",
 319                 #     "Legtime_f=\"$Legtime_f\"\n",
 320                 #     "Legspeed_f=\"$Legspeed_f\"\n",
 321                 #     "Legcourse_f=\"$Legcourse_f\" \x7D\x7D\x7D\n",
 322                 # ));
 323                 my ($NS, $WE,
 324                     $Alt_unit,
 325                     $Leglength,
 326                     $Legtime_hour, $Legtime_min, $Legtime_sec,
 327                     $Legspeed, $Legspeed_unit,
 328                     $Legcourse
 329                 ) = ("", "", "", "", "", "", "", "", "", "", "", "", "",
 330                      "", "", "", "", "", "", "", "", "", "");
 331                 ($Position_f =~ /^(N|S)([\d\.]+) (W|E)([\d\.]+)/) &&
 332                     ($NS = $1, $Dat{'lat'} = $2, $WE = $3, $Dat{'lon'} = $4);
 333                 ($Time_f =~ /^(\d+)\.(\d+)\.(\d+) (\d+):(\d+):(\d+) \((.+?)\)/) &&
 334                     ($Dat{'day'} = $1, $Dat{'month'} = $2, $Dat{'year'} = $3,
 335                      $Dat{'hour'} = $4, $Dat{'min'} = $5, $Dat{'sec'} = $6);
 336                 ($Alt_f =~ /^($DIGIT+) (.*?)/) &&
 337                     ($Dat{'ele'} = $1, $Alt_unit = $2);
 338                 # D("ele = \"$Dat{'ele'}\"");
 339                 ($NS eq "S") && ($Dat{'lat'} = 0-$Dat{'lat'});
 340                 ($WE eq "W") && ($Dat{'lon'} = 0-$Dat{'lon'});
 341                 # MapSource in win xp writes YYYY, but YY in win98se.
 342                 (
 343                     defined($Dat{'year'})
 344                         && $Dat{'year'} =~ /\d/
 345                         && $Dat{'year'} < 1900
 346                 ) && ($Dat{'year'} += 2000);
 347                 print_entry(%Dat);
 348                 # }}}
 349             } elsif (/^Track\t(.*?)\t/) {
 350                 $Dat{'title'} = txt_to_xml($1);
 351                 $Dat{'what'} = "title";
 352                 $found_move = 1;
 353                 print_entry(%Dat);
 354             } elsif (
 355                 m{^
 356                     T\t
 357                     (\d\d)/(\d\d)/(\d\d\d\d)\ (\d\d):(\d\d):(\d\d)\t
 358                     (.+)\xB0(.+)'(.+)"\t
 359                     (.+)\xB0(.+)'(.+)"
 360                 }x
 361             ) {
 362                 # T 09/01/2002 11:51:26 60°23'36.3" 5°19'35.9" {{{
 363                 my ($lat_d, $lat_m, $lat_s, $lon_d, $lon_m, $lon_s);
 364                 ($Dat{'month'}, $Dat{'day'}, $Dat{'year'},
 365                  $Dat{'hour'}, $Dat{'min'}, $Dat{'sec'},
 366                  $lat_d, $lat_m, $lat_s,
 367                  $lon_d, $lon_m, $lon_s) =
 368                 ($1, $2, $3,
 369                  $4, $5, $6,
 370                  $7, $8, $9,
 371                  $10, $11, $12);
 372                 my $Flat = defined($Round{'lat'}) ? ".$Round{'lat'}" : "";
 373                 my $Flon = defined($Round{'lon'}) ? ".$Round{'lon'}" : "";
 374                 $Dat{'lat'} = sprintf("%${Flat}f",
 375                                       1.0*($lat_d+($lat_m/60)+($lat_s/3600)));
 376                 $Dat{'lon'} = sprintf("%${Flon}f",
 377                                       1.0*$lon_d+($lon_m/60)+($lon_s/3600));
 378                 print_entry(%Dat);
 379                 # }}}
 380             } elsif (
 381                 m{^
 382                     1\ (\S+)\ (\S+)\ (\S+)\ (\S+)\x20
 383                     (\d\d)/(\d\d)/(\d\d\d\d)\ (\d\d):(\d\d):(\d\d)
 384                 }x
 385             ) {
 386                 # 1 60.3938222 5.3238754 17.3 0 09/01/2002 14:18:23 {{{
 387                 ($Dat{'lat'}, $Dat{'lon'}, $Dat{'speed'},
 388                  $Dat{'unkn'},
 389                  $Dat{'month'}, $Dat{'day'}, $Dat{'year'},
 390                  $Dat{'hour'}, $Dat{'min'}, $Dat{'sec'}) =
 391                 ($1, $2, $3,
 392                  $4,
 393                  $5, $6, $7,
 394                  $8, $9, $10);
 395                 print_entry(%Dat);
 396                 # }}}
 397             } elsif (/^
 398                 # @020721221336N6048353E00701826S015-00001E4859N1673U0000 {{{
 399                     # Regexp {{{
 400                     (@) # @
 401                     (\d\d) # Year
 402                     (\d\d) # Month
 403                     (\d\d) # Day
 404                     (\d\d) # Hours
 405                     (\d\d) # Minutes
 406                     (\d\d) # Seconds
 407                     ([NS]) # N|S
 408                     (\d\d) # Latitude degree
 409                     (\d\d) # Latitude minute
 410                     (\d\d\d) # Latitude minute decimals
 411                     ([EW]) # E|W
 412                     (\d\d\d) # Longitude degree
 413                     (\d\d) # Longitude minute
 414                     (\d\d\d) # Longitude minute degree
 415                     (....) # Accurancy
 416                     (......) # Elevation
 417                     (...............)
 418                     # }}}
 419                     /x) {
 420                 my ($NS, $EW, $lat_deg, $lat_degmin, $lat_mindec, $lon_deg,
 421                     $lon_degmin, $lon_mindec);
 422                 ($Dat{'year'}, $Dat{'month'}, $Dat{'day'}, $Dat{'hour'},
 423                  $Dat{'min'}, $Dat{'sec'}, $NS, $lat_deg,
 424                  $lat_degmin, $lat_mindec, $EW,
 425                  $lon_deg, $lon_degmin, $lon_mindec,
 426                  $Dat{'accur'}, $Dat{'ele'}, $Dat{'unknown'}) =
 427                 ($2+2000, $3, $4, $5,
 428                  $6, $7, $8, $9,
 429                  $10, $11, $12,
 430                  $13, $14, $15,
 431                  $16, $17, $18);
 432                 my $ep_time = timegm_nocheck(
 433                     $Dat{'sec'}, $Dat{'min'}, $Dat{'hour'},
 434                     $Dat{'day'}, $Dat{'month'}-1, $Dat{'year'}
 435                 );
 436                 $last_time = $ep_time;
 437                 my $Flat = defined($Round{'lat'}) ? ".$Round{'lat'}" : "";
 438                 my $Flon = defined($Round{'lon'}) ? ".$Round{'lon'}" : "";
 439                 my $tmp_lon = sprintf(
 440                                       "%${Flon}f",
 441                                       $lon_deg +
 442                                       $lon_degmin/60 +
 443                                       $lon_mindec/60000);
 444                 my $tmp_lat = sprintf("%${Flat}f",
 445                                       $lat_deg +
 446                                       $lat_degmin/60 +
 447                                       $lat_mindec/60000);
 448                 ($NS eq "S") && ($tmp_lat = 0-$tmp_lat);
 449                 ($EW eq "W") && ($tmp_lon = 0-$tmp_lon);
 450                 $Dat{'lat'} = $tmp_lat;
 451                 $Dat{'lon'} = $tmp_lon;
 452                 print_entry(%Dat);
 453                 # }}}
 454             } elsif (/^(@)(\d\d)(\d\d)(\d\d)(\d\d)(\d\d)(\d\d)(_{42})/) {
 455                 # @020721221336__________________________________________ {{{
 456                 ($Dat{'year'}, $Dat{'month'}, $Dat{'day'},
 457                  $Dat{'hour'}, $Dat{'min'}, $Dat{'sec'}, $Dat{'rest'}) =
 458                 ($2+2000, $3, $4,
 459                  $5, $6, $7, $8);
 460                 $Dat{'error'} = "nosignal";
 461                 print_entry(%Dat);
 462                 # }}}
 463             } elsif (/^xmaplog /) {
 464                 # NOP
 465             } elsif (/^$/) {
 466                 ($Opt{'output-format'} eq "csv")
 467                     && ($Opt{'save-to-file'} eq "\n")
 468                     && print("\n");
 469             } elsif (/^Pause: /) {
 470                 # NOP, is here to cope with old files I’ve lying around.
 471             } elsif ($Dat{'error'} eq "desc") {
 472                 my $Comment = $_;
 473                 if (defined($Comment)) {
 474                     $Comment =~ s/^\s*(.*?)\s*$/$1/;
 475                     if ($Opt{'output-format'} eq "gpsml") {
 476                         $Dat{'desc'} = txt_to_xml($Comment);
 477                         $Dat{'what'} = "desc";
 478                         print_entry(%Dat);
 479                     }
 480                 }
 481             } else {
 482                 $Opt{'verbose'} && warn("Line $.: Unknown: \"$_\"\n");
 483             }
 484         }
 485         # }}}
 486     } else {
 487         warn("$progname: $curr_file: Cannot open file for read: $!\n");
 488     }
 489     # }}}
 490 }
 491
 492 print_footer(*STDOUT);
 493
 494 exit(0);
 495
 496 sub read_xmlfile {
 497     # {{{
 498     my $Txt = join("", @_);
 499     $Txt =~ s/<!--(.*?)-->//gs;
 500     $Txt =~ s#(<gpx\b.*?>.*?</gpx>)#print_gpx($1)#gse;
 501     # }}}
 502 }
 503
 504 sub print_gpx {
 505     # {{{
 506     my $Orig = shift;
 507     my $Str = $Orig;
 508     # D("print_xml_gps(\"$Orig\")\n");
 509     $Str =~ s/<!--(.*?)-->//gs;
 510     my $fromdate_str = "";
 511     if ($Opt{'from-date'}) {
 512         $fromdate_str = "date >= '$Opt{'from-date'}' AND ";
 513     }
 514     if ($Opt{'output-format'} =~ /^(pgwtab|pgwupd)$/) {
 515         # {{{
 516         $Str =~
 517         s{
 518             <wpt\b(.*?)>(.*?)</wpt>
 519         }
 520         {
 521             my $attr_wpt = $1;
 522             my $el_wpt = $2;
 523             my ($Lat, $Lon, $Name, $Ele, $Type, $Time, $Cmt, $Desc, $Src, $Sym) =
 524                ('\N', '\N',  '\N', '\N',  '\N',  '\N', '\N',  '\N', '\N', '\N');
 525
 526             $attr_wpt =~ /.*lat="($DIGIT+?)"/s &&
 527                 ($Lat = postgresql_copy_safe($1));
 528             $attr_wpt =~ /.*lon="($DIGIT+?)"/s &&
 529                 ($Lon = postgresql_copy_safe($1));
 530             $el_wpt =~ /.*<name\b(.*?)>(.*?)<\/name>/s &&
 531                 ($Name = postgresql_copy_safe(xml_to_txt($2)));
 532             $el_wpt =~ /.*<ele\b(.*?)>(.*?)<\/ele>/s &&
 533                 ($Ele = postgresql_copy_safe(xml_to_txt($2)));
 534             $el_wpt =~ /.*<type\b(.*?)>(.*?)<\/type>/s &&
 535                 ($Type = postgresql_copy_safe(xml_to_txt($2)));
 536             $el_wpt =~ /.*<time\b(.*?)>(.*?)<\/time>/s &&
 537                 ($Time = postgresql_copy_safe(xml_to_txt($2)));
 538             $el_wpt =~ /.*<cmt\b(.*?)>(.*?)<\/cmt>/s &&
 539                 ($Cmt = postgresql_copy_safe(xml_to_txt($2)));
 540             $el_wpt =~ /.*<desc\b(.*?)>(.*?)<\/desc>/s &&
 541                 ($Desc = postgresql_copy_safe(xml_to_txt($2)));
 542             $el_wpt =~ /.*<src\b(.*?)>(.*?)<\/src>/s &&
 543                 ($Src = postgresql_copy_safe(xml_to_txt($2)));
 544             $el_wpt =~ /.*<sym\b(.*?)>(.*?)<\/sym>/s &&
 545                 ($Sym = postgresql_copy_safe(xml_to_txt($2)));
 546
 547             if (length($Opt{'round'})) {
 548                 if (defined($Round{'lat'}) && length($Lat)) {
 549                     ($Lat = 1.0 * sprintf("%.$Round{'lat'}f", $Lat));
 550                 }
 551                 if (defined($Round{'lon'}) && length($Lon)) {
 552                     ($Lon = 1.0 * sprintf("%.$Round{'lon'}f", $Lon));
 553                 }
 554                 if (defined($Round{'ele'}) && $Ele ne '\N') {
 555                     ($Ele = 1.0 * sprintf("%.$Round{'ele'}f", $Ele));
 556                 }
 557             }
 558
 559             if ($Opt{'output-format'} eq "pgwtab") {
 560                 print(
 561                     join("\t",
 562                         "($Lat,$Lon)",
 563                         $Name,
 564                         $Ele,
 565                         $Type,
 566                         $Time,
 567                         $Cmt,
 568                         $Desc,
 569                         $Src,
 570                         $Sym
 571                     ) . "\n"
 572                 );
 573             } elsif ($Opt{'output-format'} eq "pgwupd") {
 574                 $Name =~ s/'/''/gs;
 575                 print(join("\n",
 576                     "BEGIN;",
 577                     "$Spc${Spc}UPDATE logg SET name = clname(coor) " .
 578                         "WHERE $fromdate_str(point($Lat,$Lon) <-> coor) < 0.05;",
 579                     "$Spc${Spc}UPDATE logg SET dist = cldist(coor) " .
 580                         "WHERE $fromdate_str(point($Lat,$Lon) <-> coor) < 0.05;",
 581                     "COMMIT;"
 582                 ) . "\n");
 583             }
 584             "";
 585         }gsex;
 586         # }}}
 587     } else {
 588         # {{{
 589         $Str =~
 590         s{
 591             <trk\b(.*?)>(.*?)</trk>
 592         }
 593         {
 594             my $el_trk = $2;
 595             $el_trk =~
 596             s{
 597                 <name\b(.*?)>(.*?)</name>
 598             }{
 599                 my %tmp_dat = ();
 600                 $tmp_dat{'title'} = $2;
 601                 $tmp_dat{'what'} = "title";
 602                 $tmp_dat{'error'} = "";
 603                 print_entry(%tmp_dat);
 604                 "";
 605             }sex;
 606             $el_trk =~
 607             s{
 608                 <trkseg\b(.*?)>(.*?)</trkseg>
 609             }
 610             {
 611                 my $el_trkseg = $2;
 612                 $el_trkseg =~
 613                 s{
 614                     <trkpt\b(.*?)>(.*?)</trkpt>
 615                 }
 616                 {
 617                     my ($attr_trkpt, $el_trkpt) =
 618                        (         $1,        $2);
 619                     %Dat = (
 620                         'year' => '', 'month' => '', 'day' => '',
 621                         'hour' => '', 'min' => '', 'sec' => '',
 622                         'epoch' => '',
 623                         'date-format' => '',
 624                         'lat' => '', 'lon' => '',
 625                         'ele' => '',
 626                         'desc' => '',
 627                         'error' => "",
 628                         'what' => 'tp',
 629                     );
 630                     ($attr_trkpt =~ /\blon="(.*?)"/) && ($Dat{'lon'} = $1);
 631                     ($attr_trkpt =~ /\blat="(.*?)"/) && ($Dat{'lat'} = $1);
 632                     ($el_trkpt   =~ m#<ele\b.*?>(.*?)</ele>#) && ($Dat{'ele'} = $1);
 633                     if (
 634                         $el_trkpt =~
 635                         m{
 636                             <time>(\d\d\d\d)-?(\d\d)-?(\d\d)T
 637                             (\d\d):?(\d\d):?([\d\.]+)Z</time>
 638                         }x
 639                     ) {
 640                         ($Dat{'year'}, $Dat{'month'}, $Dat{'day'},
 641                          $Dat{'hour'}, $Dat{'min'}, $Dat{'sec'}) =
 642                         ($1, $2, $3, $4, $5, $6);
 643                     }
 644                     print_entry(%Dat);
 645                     "";
 646                 }gsex;
 647                 $found_move = 1;
 648             }gsex;
 649             $found_move = 1;
 650         }gsex;
 651         # }}}
 652     }
 653     # }}}
 654 }
 655
 656 sub print_header {
 657     # {{{
 658     local *OutFP = shift;
 659     if ($Opt{'output-format'} eq "gpsml") {
 660         print(OutFP join("",
 661             "<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n",
 662             "<gpsml>\n",
 663             "<track>\n",
 664         ));
 665     } elsif ($Opt{'output-format'} eq "gpstrans") {
 666         print(OutFP "Format: DMS  UTC Offset:   0.00 hrs  " .
 667                     "Datum[100]: WGS 84\n");
 668     } elsif ($Opt{'output-format'} eq "gpx") {
 669         print(OutFP join("",
 670             qq{<?xml version="1.0" encoding="UTF-8" standalone="no"?>\n},
 671             qq{<gpx\n},
 672             qq{$Spc${Spc}version="1.1"\n},
 673             qq{$Spc${Spc}creator="gpst - http://sunny256.github.com/gpstools/"\n},
 674             qq{$Spc${Spc}xmlns="http://www.topografix.com/GPX/1/1"\n},
 675             qq{$Spc${Spc}xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"\n},
 676             qq{$Spc${Spc}xsi:schemaLocation="http://www.topografix.com/GPX/1/1 },
 677             qq{http://www.topografix.com/GPX/1/1/gpx.xsd"\n},
 678             qq{>\n},
 679             qq{$Spc$Spc<trk>\n},
 680             qq{$Spc$Spc$Spc$Spc<trkseg>\n},
 681         ));
 682     } elsif ($Opt{'output-format'} eq "ps") {
 683         print(OutFP ps_header(532, 6034, 533, 6040));
 684         print(OutFP "*u\n");
 685     } elsif ($Opt{'output-format'} eq "svg") {
 686         print(OutFP join("",
 687             "<?xml version=\"1.0\" standalone=\"no\"?>\n",
 688             "<!DOCTYPE svg PUBLIC \"-//W3C//DTD SVG 1.1//EN\"\n",
 689             "$Spc$Spc\"http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd\">\n",
 690             "<svg height=\"1000\" width=\"1000\" viewBox=\"23 70 2 2\"\n",
 691             "$Spc${Spc}xmlns=\"http://www.w3.org/2000/svg\" version=\"1.1\">\n",
 692             "$Spc$Spc<title></title>\n",
 693             "$Spc$Spc<desc></desc>\n",
 694         ));
 695     }
 696     # }}}
 697 }
 698
 699 sub print_footer {
 700     # Print footer {{{
 701     local *OutFP = shift;
 702     if ($Opt{'output-format'} eq "gpsml") {
 703         print(OutFP join("",
 704             "</track>\n",
 705             "</gpsml>\n",
 706         ));
 707     } elsif ($Opt{'output-format'} eq "gpx") {
 708         print(OutFP join("",
 709             "$Spc$Spc$Spc$Spc</trkseg>\n",
 710             "$Spc$Spc</trk>\n",
 711             "</gpx>\n",
 712         ));
 713     } elsif ($Opt{'output-format'} eq "poscount") {
 714         while (my ($l_name, $l_val) = each %Poscount) {
 715             $l_name =~ /^(.+?),(.+?)$/
 716                 && print(OutFP "$1\t$2\t$l_val\n");
 717         }
 718     } elsif ($Opt{'output-format'} eq "ps") {
 719         print(OutFP join("",
 720             "*U\n",
 721             "%%Trailer\n",
 722             "%%EOF\n",
 723         ));
 724     } elsif ($Opt{'output-format'} eq "svg") {
 725         print(OutFP "\"/>\n</svg>\n");
 726     }
 727     # }}}
 728 }
 729
 730 sub print_entry {
 731     # Print a GPS entry with time, latitude, longitude and elevation in
 732     # various formats
 733     # {{{
 734     my %Dat = @_;
 735     defined($Dat{'desc'}) || ($Dat{'desc'} = "");
 736     defined($Dat{'ele'}) || ($Dat{'ele'} = "");
 737     defined($Dat{'lat'}) || ($Dat{'lat'} = "");
 738     defined($Dat{'lon'}) || ($Dat{'lon'} = "");
 739     defined($Dat{'year'}) || ($Dat{'year'} = "");
 740     my $print_time = length($Dat{'year'}) ? 1 : 0;
 741     my $print_pos;
 742     if (!$Req{'position'} && $Opt{'output-format'} eq "gpsml") {
 743         $print_pos = (length($Dat{'lat'}) || length($Dat{'lon'})) ? 1 : 0;
 744     } else {
 745         $print_pos = (length($Dat{'lat'}) && length($Dat{'lon'})) ? 1 : 0;
 746     }
 747     if (!$print_pos) {
 748         $Dat{'lat'} = $Dat{'lon'} = "";
 749     }
 750     my $print_ele = length($Dat{'ele'}) ? 1 : 0;
 751     my $print_desc = length($Dat{'desc'}) ? 1 : 0;
 752     my $Line = "";
 753     # D("print_entry(\"" . join("\", \"", @_) . "\");");
 754     my $ep_time;
 755
 756     if (length($Opt{'round'})) {
 757         for my $Tmp (qw{ lat lon ele }) {
 758             if (defined($Round{$Tmp}) && length($Dat{$Tmp})) {
 759                 # D("Tmp = '$Tmp'");
 760                 ($Dat{$Tmp} = num_expand(1.0 * sprintf("%.$Round{$Tmp}f", $Dat{$Tmp})));
 761             }
 762         }
 763     }
 764
 765     if ($Opt{'output-format'} eq "poscount") {
 766         if (!length($Dat{'error'})) {
 767             my $Name = "$Dat{'lon'},$Dat{'lat'}";
 768             defined($Poscount{$Name}) || ($Poscount{$Name} = 0);
 769             $Poscount{$Name}++;
 770         }
 771         return;
 772     }
 773
 774     if ($print_time) {
 775         $ep_time = timegm_nocheck(
 776             $Dat{'sec'}, $Dat{'min'}, $Dat{'hour'},
 777             $Dat{'day'}, $Dat{'month'} - 1, $Dat{'year'}
 778         );
 779         if ($Opt{'time-shift'}) {
 780             # D("ep_time før: '$ep_time'");
 781             $ep_time += $Opt{'time-shift'};
 782             # D("ep_time etter: '$ep_time'");
 783             ($Dat{'sec'}, $Dat{'min'},$Dat{'hour'}, $Dat{'day'},
 784              $Dat{'month'}, $Dat{'year'}) = gmtime($ep_time);
 785             $Dat{'year'} += 1900;
 786             $Dat{'month'}++;
 787         }
 788         $Dat{'epoch'} = $ep_time;
 789         $Dat{'year'}  = sprintf("%04u", $Dat{'year'});
 790         $Dat{'month'} = sprintf("%02u", $Dat{'month'});
 791         $Dat{'day'}   = sprintf("%02u", $Dat{'day'});
 792         $Dat{'hour'}  = sprintf("%02u", $Dat{'hour'});
 793         $Dat{'min'}   = sprintf("%02u", $Dat{'min'});
 794         $Dat{'sec'}   = sprintf("%02u", $Dat{'sec'});
 795         if ($Opt{'chronology'}) {
 796             if ($last_time > $ep_time && !length($Dat{'error'})) {
 797                 warn(sprintf(
 798                     "%s: $curr_file: \"%sZ\": Next date is %s in the past (%sZ)\n",
 799                     $progname, sec_to_string($last_time, "T"),
 800                     sec_to_readable($last_time-$ep_time),
 801                     sec_to_string($ep_time, "T")
 802                 ));
 803                 # FIXME: Make --fix work with gpx.
 804                 if ($Opt{'fix'} && ($Opt{'output-format'} !~ /^gpx$/)) {
 805                     $Dat{'error'} = "chrono";
 806                 }
 807             } elsif ($last_time == $ep_time && !length($Dat{'error'})) {
 808                 warn(sprintf(
 809                     "%s: $curr_file: \"%sZ\": Duplicated time\n",
 810                     $progname, sec_to_string($last_time, "T")
 811                 ));
 812                 # FIXME: Make --fix work with gpx.
 813                 if ($Opt{'fix'} && ($Opt{'output-format'} !~ /^gpx$/)) {
 814                     $Dat{'error'} = "duptime";
 815                 }
 816             }
 817         }
 818     } else {
 819         $ep_time = 0;
 820         $Dat{'year'} = 0;
 821         $Dat{'month'} = 0;
 822         $Dat{'day'} = 0;
 823         $Dat{'hour'} = 0;
 824         $Dat{'min'} = 0;
 825         $Dat{'sec'} = 0;
 826     }
 827
 828     if ($Opt{'save-to-file'} ne "\n") {
 829         # {{{
 830         $print_time || return;
 831         my $base_name = "$Dat{'year'}$Dat{'month'}$Dat{'day'}T" .
 832                         "$Dat{'hour'}$Dat{'min'}$Dat{'sec'}Z" .
 833                         "$Opt{'save-to-file'}";
 834         my $file_name = $base_name;
 835         if (-e $file_name) {
 836             for (my $a = 1; (-e $file_name) && ($a < 1000); $a++) {
 837                 $file_name = "$base_name.dup_$a";
 838             }
 839             if (-e $file_name) {
 840                 die("$progname: $base_name: File already exists, and ran " .
 841                     "out of attempts to create unique file name\n");
 842             }
 843             if ($Opt{'verbose'}) {
 844                 warn("$progname: $base_name: File already exists, using " .
 845                      "unique name \"$file_name\" instead\n");
 846             }
 847         }
 848         if (open(ToFP, ">", $file_name)) {
 849             print_header(*ToFP);
 850             print(ToFP (
 851                 $from_stdin
 852                     ? @first_lines
 853                     : ()),
 854                 (length($xml_data)
 855                     ? $xml_data
 856                     : <>)
 857             ) || die("$progname: $file_name: Cannot write to file: $!\n");
 858             print_footer(*ToFP);
 859             close(ToFP);
 860             if ($Opt{'output-format'} eq "gpsml") {
 861                 printf("<include>%s</include>\n",
 862                        txt_to_xml($file_name));
 863             } elsif ($Opt{'output-format'} eq "gpx") {
 864                 printf("<!-- Saved unconverted data to \"%s\" -->\n",
 865                        txt_to_xml($file_name));
 866             } else {
 867                 print("$progname: Saved unconverted data to \"$file_name\"\n");
 868             }
 869             exit 0;
 870         } else {
 871             die("$progname: $file_name: Cannot create file: $!\n");
 872         }
 873         # }}}
 874     }
 875
 876     my $pause_len = 0;
 877     my $do_print = 1;
 878
 879     if ($Dat{'what'} eq "tp") {
 880         # {{{
 881         if ($Opt{'require'}) {
 882             $Req{'time'} && !$print_time && return;
 883             $Req{'position'} && !$print_pos && return;
 884             $Req{'ele'} && !$print_ele && return;
 885         }
 886
 887         if ($Opt{'inside'} || $Opt{'outside'}) {
 888             if (
 889                 ($Dat{'lat'} < $lat1) ||
 890                 ($Dat{'lat'} > $lat2) ||
 891                 ($Dat{'lon'} < $lon1) ||
 892                 ($Dat{'lon'} > $lon2)
 893             ) {
 894                 $Opt{'inside'} && return;
 895             } else {
 896                 $Opt{'outside'} && return;
 897             }
 898         }
 899
 900         if ($Opt{'output-format'} eq "ps") {
 901             $Dat{'lon'} *= 100;
 902             $Dat{'lat'} *= 100;
 903         }
 904
 905         if (
 906             $Opt{'skip-dups'}
 907             && ($Dat{'lon'} eq $last_lon)
 908             && ($Dat{'lat'} eq $last_lat)
 909         ) {
 910             if ($Opt{'output-format'} eq 'gpsml') {
 911                 $Dat{'error'} = "dup";
 912             } else {
 913                 $do_print = 0;
 914             }
 915         } else {
 916                 $do_print = 1;
 917         }
 918
 919         if (
 920             $Opt{'create-breaks'}
 921             && $ep_time-$last_time > $PAUSE_LIMIT
 922             && $last_time
 923         ) {
 924             $pause_len = $ep_time-$last_time;
 925             # D("pause_len set to '$pause_len'");
 926         }
 927
 928         $Line .= pause_entry($pause_len, $ep_time, $last_time);
 929         # }}}
 930     }
 931
 932     if ($do_print) {
 933         # Valid data was found, send to stdout {{{
 934         unless ($first_time) {
 935             $first_time = $ep_time;
 936         }
 937         $Line .= gen_entry($print_pos, $pause_len, $print_time, $ep_time, $print_ele, %Dat);
 938         # }}}
 939     }
 940
 941     if (!$last_time && $Opt{'output-format'} eq "ps") {
 942         $Line .= "$Dat{'lon'} $Dat{'lat'} m\n";
 943     }
 944
 945     if ($do_print) {
 946         if ($found_move) {
 947             if ($Opt{'output-format'} eq "gpsml") {
 948                 $Line = "<break/>\n$Line";
 949             }
 950             (!$pause_len && ($Opt{'output-format'} eq "xgraph"))
 951                 && ($Line .= "move $Line");
 952             ($Opt{'output-format'} eq "clean") && ($Line .= "\n");
 953             if ($Opt{'output-format'} eq "gpx") {
 954                 $Line .= "$Spc$Spc$Spc$Spc</trkseg>\n" .
 955                          "$Spc$Spc$Spc$Spc<trkseg>\n";
 956             }
 957             $found_move = 0;
 958         }
 959         print($Line);
 960     }
 961     $print_time && ($last_time = $ep_time);
 962     if ($print_pos) {
 963         $last_lon  = $Dat{'lon'};
 964         $last_lat  = $Dat{'lat'};
 965     }
 966     $last_line = $data_line;
 967     $svg_start_thing = "\"/>\n";
 968     # }}}
 969 }
 970
 971 sub gen_entry {
 972     # Generate trackpoint entry, calls trackpoint() {{{
 973     my ($print_pos, $pause_len, $print_time, $ep_time, $print_ele, %Dat) = @_;
 974     my $Line = "";
 975     if ($Opt{'double-y-scale'} && length($Dat{'lat'})) {
 976         $Dat{'lat'} *= 2;
 977     }
 978     if ($Opt{'output-format'} eq "gpsml") {
 979         if ($Dat{'what'} eq "tp") {
 980             $Dat{'format'} = "gpsml";
 981             $Line .= trackpoint(%Dat);
 982         } elsif ($Dat{'what'} =~ /^(pause|desc|title)$/) {
 983             $Line .= sprintf("<%s>%s</%s>\n",
 984                              $1,
 985                              $Dat{$1},
 986                              $1);
 987         }
 988     } elsif ($Opt{'output-format'} eq "pgtab") {
 989         if ($Dat{'what'} eq "tp" && !length($Dat{'error'})) {
 990             $Dat{'format'} = "pgtab";
 991             $Line .= trackpoint(%Dat);
 992         }
 993     } elsif ($Opt{'output-format'} eq "xgraph") {
 994         if ($print_pos && !length($Dat{'error'})) {
 995             $Dat{'format'} = "xgraph";
 996             $Line .= trackpoint(%Dat);
 997         }
 998     } elsif($Opt{'output-format'} eq "gpstrans") {
 999         if ($print_pos && !length($Dat{'error'})) {
1000             $Dat{'format'} = "gpstrans";
1001             $Line .= trackpoint(%Dat);
1002         }
1003     } elsif($Opt{'output-format'} eq "gpx") {
1004         if ($Dat{'what'} eq "tp") {
1005             $Dat{'format'} = "gpx";
1006             $Line .= trackpoint(%Dat);
1007         }
1008     } elsif ($Opt{'output-format'} eq "clean") {
1009         if ($Dat{'what'} eq "tp" && !length($Dat{'error'})) {
1010             $Dat{'format'} = "clean";
1011             $Line .= trackpoint(%Dat);
1012         }
1013     } elsif ($Opt{'output-format'} eq "ps") {
1014         $Line .= (
1015             $pause_len
1016                 ? "f\n$Dat{'lon'} $Dat{'lat'} m\n"
1017                 : "$Dat{'lon'} $Dat{'lat'} l\n"
1018         );
1019     } elsif ($Opt{'output-format'} eq "svg") {
1020         $Line .= (
1021             ($last_lon == 1000) || $pause_len
1022                 ? join("",
1023                     "$svg_start_thing<path\n",
1024                     "  stroke=\"blue\"\n",
1025                     "  stroke-width=\"0.001\"\n",
1026                     "  fill=\"none\"\n",
1027                     "  d=\"\n",
1028                     "M $Dat{'lon'} $Dat{'lat'}\n")
1029                 : "L $Dat{'lon'} $Dat{'lat'}\n"
1030         );
1031     } elsif ($Opt{'output-format'} eq "ygraph") {
1032         if (!length($Dat{'error'})) {
1033             my $Time = $print_time ? ($ep_time - $first_time) * 1 : 0;
1034             $Line .= "\"Time = $Time.0\n$Dat{'lon'} $Dat{'lat'}\n\n";
1035         }
1036     } elsif ($Opt{'output-format'} eq "csv") {
1037         # {{{
1038         if (!length($Dat{'error'})) {
1039             $Dat{'format'} = "csv";
1040             $Dat{'lat'} = num_expand($Dat{'lat'});
1041             $Dat{'lon'} = num_expand($Dat{'lon'});
1042             $Dat{'ele'} = num_expand($Dat{'ele'});
1043             $Line .= join("\t",
1044                 $print_time
1045                     ? $Opt{'epoch'}
1046                         ? $ep_time
1047                         : $Opt{'short-date'}
1048                             ? "$Dat{'year'}$Dat{'month'}$Dat{'day'}T" .
1049                               "$Dat{'hour'}$Dat{'min'}$Dat{'sec'}Z"
1050                             : "$Dat{'year'}-$Dat{'month'}-$Dat{'day'}T" .
1051                               "$Dat{'hour'}:$Dat{'min'}:$Dat{'sec'}Z"
1052                     : "",
1053                 $Dat{'lon'},
1054                 $Dat{'lat'},
1055                 $print_ele ? $Dat{'ele'} : "", # Elevation
1056                 "\n"
1057             );
1058         }
1059         # }}}
1060     } elsif ($Opt{'output-format'} eq "pgwtab") {
1061         # FIXME: NOP at the moment.
1062     } else {
1063         die("$progname: \"$Opt{'output-format'}\": " .
1064             "Unknown output format\n");
1065     }
1066     return($Line);
1067     # }}}
1068 } # gen_entry()
1069
1070 sub pause_entry {
1071     # {{{
1072     my ($pause_len, $ep_time, $last_time) = @_;
1073     my $Line = "";
1074     if ($pause_len) {
1075         if ($Opt{'output-format'} eq "gpsml") {
1076             $Line .= sprintf("<pause>%s</pause>\n",
1077                              sec_to_readable($ep_time-$last_time));
1078         } elsif ($Opt{'output-format'} eq "clean") {
1079             $pause_len && ($Line .= "\n");
1080         } elsif ($Opt{'output-format'} eq "csv") {
1081             $Line .= sprintf("# Pause: %s\n# move\n",
1082                              sec_to_readable($ep_time-$last_time));
1083         } elsif ($Opt{'output-format'} eq "xgraph") {
1084             $pause_len && ($Line .= "move ");
1085         }
1086     }
1087     return($Line);
1088     # }}}
1089 } # pause_entry()
1090
1091 sub ps_header {
1092     # Send a Postscript header to stdout {{{
1093     my ($bl_lon, $bl_lat, $br_lon, $br_lat) = @_;
1094     my $Date = sec_to_string(time);
1095     return(join("",
1096         "%!PS-Adobe-3.0 EPSF-3.0\n",
1097         "%%Creator: gpst\n",
1098         "%%Title:\n",
1099         "%%CreationDate: $Date\n",
1100         "%%BoundingBox: $bl_lon $bl_lat $br_lon $br_lat\n",
1101         "%%DocumentData: Clean7Bit\n",
1102         "%%EndComments\n",
1103         "%%BeginProlog\n",
1104         "/bd { bind def } bind def\n",
1105         "/incompound false def\n",
1106         "/m { moveto } bd\n",
1107         "/l { lineto } bd\n",
1108         "/c { curveto } bd\n",
1109         "/F { incompound not {fill} if } bd\n",
1110         "/f { closepath F } bd\n",
1111         "/S { stroke } bd\n",
1112         "/*u { /incompound true def } bd\n",
1113         "/*U { /incompound false def f} bd\n",
1114         "/k { setcmykcolor } bd\n",
1115         "/K { k } bd\n",
1116         "%%EndProlog\n",
1117         "%%BeginSetup\n",
1118         "%%EndSetup\n",
1119     ));
1120     # }}}
1121 }
1122
1123 sub print_version {
1124     # Print program version {{{
1125     print("$progname v$VERSION\n");
1126     # }}}
1127 } # print_version()
1128
1129 sub usage {
1130     # Send the help message to stdout {{{
1131     my $Retval = shift;
1132
1133     if ($Opt{'verbose'}) {
1134         print("\n");
1135         print_version();
1136     }
1137     print(<<END);
1138
1139 Converts between various GPS formats.
1140
1141 Usage: $progname [options] [file [files [...]]]
1142        $progname -S [file [files [...]]]
1143        $progname -u [file [files [...]]]
1144
1145 Options:
1146
1147   --chronology
1148     Check for broken chronology, warn about entries with an old
1149     timestamp.
1150   -d, --skip-dups
1151     Skip duplicated coordinates.
1152   -e, --epoch
1153     Use seconds since 1970-01-01 00:00:00 GMT as date format.
1154   --fix
1155     Comment out entries which is obviously wrong. Use together with
1156     --chronology to fix those kind of errors. Does not work with GPX
1157     output yet.
1158   --from-date x
1159     Used by the pgwupd format. Specifies from which date waypoints
1160     should be updated. No checks for valid date format here, let
1161     PostgreSQL take care of that. All variants it understands can be
1162     used here.
1163   -h, --help
1164     Show this help.
1165   --inside
1166     Print only trackpoints inside a rectangle specified by --pos1 and
1167     --pos2.
1168   -n, --undefined x
1169     Use x as undefined value. Default: "$Udef".
1170   -o, --output-format x
1171     Use output format x:
1172       clean
1173       csv
1174       gpsml (Default)
1175       gpstrans
1176       gpx (Not complete)
1177       pgtab
1178       pgwtab
1179       pgwupd
1180       poscount
1181       ps (Unfinished)
1182       svg (Unfinished)
1183       xgraph
1184       ygraph
1185   --outside
1186     Print only trackpoints outside a rectangle specified by --pos1 and
1187     --pos2.
1188   --pos1 x
1189   --pos2 x
1190     Specifies one corner where x is in "lat,lon" format (decimal
1191     degrees, negative for west or south) of area rectangle used by the
1192     --inside and --outside options.
1193   -r, --require x
1194     Specify requirements for trackpoints to be written. x is a string
1195     with the following flags:
1196       e
1197         Print only waypoints which have an elevation.
1198       p
1199         Print only waypoints which have a position.
1200       t
1201         Print only waypoints which have a timestamp.
1202   -R, --round x=y[,x2=y2[...]]
1203     Round trackpoint element x to y decimals. Example:
1204       --round lat=4,lon=5,ele=1
1205   -s, --short-date
1206     Use short date format.
1207   -S, --save-to-file x
1208     Save the unconverted data to a file with a filename starting with
1209     the timestamp of the first trackpoint. The parameter string x is
1210     added at the end of the filename. For the time being this option
1211     will ignore all other options. Note: If several files are specified
1212     on the command line, all data will be saved into only one file. This
1213     behaviour may change in the future.
1214   -t, --create-breaks
1215     Create breaks in track between points with a difference more than
1216     $PAUSE_LIMIT seconds.
1217   -T x, --time-shift x
1218     Move time of trackpoint x seconds forwards or backwards. x can be a
1219     positive or negative integer.
1220   -v, --verbose
1221     Increase level of verbosity. Can be repeated.
1222   --version
1223     Print version information.
1224   -w, --strip-whitespace
1225     Strip all unnecessary whitespace.
1226   -y, --double-y-scale
1227     Double Y scale (latitude) to get it right in gnuplot.
1228   --debug
1229     Print debugging messages.
1230
1231 END
1232     exit($Retval);
1233     # }}}
1234 } # usage()
1235
1236 sub msg {
1237     # Print a status message to stderr based on verbosity level {{{
1238     my ($verbose_level, $Txt) = @_;
1239
1240     if ($Opt{'verbose'} >= $verbose_level) {
1241         print(STDERR "$progname: $Txt\n");
1242     }
1243     # }}}
1244 } # msg()
1245
1246 __END__
1247
1248 # Law talk {{{
1249 # Copyleft © Øyvind A. Holm <sunny@sunbase.org>
1250 #
1251 This program is free software: you can redistribute it and/or modify it
1252 under the terms of the GNU General Public License as published by the
1253 Free Software Foundation, either version 3 of the License, or (at your
1254 option) any later version.
1255
1256 This program is distributed in the hope that it will be useful, but
1257 WITHOUT ANY WARRANTY; without even the implied warranty of
1258 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
1259 See the GNU General Public License for more details.
1260
1261 You should have received a copy of the GNU General Public License along
1262 with this program.
1263 If not, see L<http://www.gnu.org/licenses/>.
1264 # }}}
1265
1266 # vim: set fenc=UTF-8 ft=perl fdm=marker ts=4 sw=4 sts=4 et fo+=w :