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(my $curr_fp, "<$curr_file")) {
 192         # {{{
 193         while (<$curr_fp>) {
 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 .= <$curr_fp>;
 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("", <$curr_fp>);
 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     my $out_fp = shift;
 659     if ($Opt{'output-format'} eq "gpsml") {
 660         print($out_fp 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($out_fp "Format: DMS  UTC Offset:   0.00 hrs  " .
 667                     "Datum[100]: WGS 84\n");
 668     } elsif ($Opt{'output-format'} eq "gpx") {
 669         print($out_fp 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($out_fp ps_header(532, 6034, 533, 6040));
 684         print($out_fp "*u\n");
 685     } elsif ($Opt{'output-format'} eq "svg") {
 686         print($out_fp 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     my $out_fp = shift;
 702     if ($Opt{'output-format'} eq "gpsml") {
 703         print($out_fp join("",
 704             "</track>\n",
 705             "</gpsml>\n",
 706         ));
 707     } elsif ($Opt{'output-format'} eq "gpx") {
 708         print($out_fp 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($out_fp "$1\t$2\t$l_val\n");
 717         }
 718     } elsif ($Opt{'output-format'} eq "ps") {
 719         print($out_fp join("",
 720             "*U\n",
 721             "%%Trailer\n",
 722             "%%EOF\n",
 723         ));
 724     } elsif ($Opt{'output-format'} eq "svg") {
 725         print($out_fp "\"/>\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         # FIXME: Sort output in some way
 767         if (!length($Dat{'error'})) {
 768             $Dat{'lat'} = num_expand($Dat{'lat'});
 769             $Dat{'lon'} = num_expand($Dat{'lon'});
 770             my $Name = "$Dat{'lon'},$Dat{'lat'}";
 771             defined($Poscount{$Name}) || ($Poscount{$Name} = 0);
 772             $Poscount{$Name}++;
 773         }
 774         return;
 775     }
 776
 777     if ($print_time) {
 778         $ep_time = timegm_nocheck(
 779             $Dat{'sec'}, $Dat{'min'}, $Dat{'hour'},
 780             $Dat{'day'}, $Dat{'month'} - 1, $Dat{'year'}
 781         );
 782         if ($Opt{'time-shift'}) {
 783             # D("ep_time før: '$ep_time'");
 784             $ep_time += $Opt{'time-shift'};
 785             # D("ep_time etter: '$ep_time'");
 786             ($Dat{'sec'}, $Dat{'min'},$Dat{'hour'}, $Dat{'day'},
 787              $Dat{'month'}, $Dat{'year'}) = gmtime($ep_time);
 788             $Dat{'year'} += 1900;
 789             $Dat{'month'}++;
 790         }
 791         $Dat{'epoch'} = $ep_time;
 792         $Dat{'year'}  = sprintf("%04u", $Dat{'year'});
 793         $Dat{'month'} = sprintf("%02u", $Dat{'month'});
 794         $Dat{'day'}   = sprintf("%02u", $Dat{'day'});
 795         $Dat{'hour'}  = sprintf("%02u", $Dat{'hour'});
 796         $Dat{'min'}   = sprintf("%02u", $Dat{'min'});
 797         $Dat{'sec'}   = sprintf("%02u", $Dat{'sec'});
 798         if ($Opt{'chronology'}) {
 799             if ($last_time > $ep_time && !length($Dat{'error'})) {
 800                 warn(sprintf(
 801                     "%s: $curr_file: \"%sZ\": Next date is %s in the past (%sZ)\n",
 802                     $progname, sec_to_string($last_time, "T"),
 803                     sec_to_readable($last_time-$ep_time),
 804                     sec_to_string($ep_time, "T")
 805                 ));
 806                 # FIXME: Make --fix work with gpx.
 807                 if ($Opt{'fix'} && ($Opt{'output-format'} !~ /^gpx$/)) {
 808                     $Dat{'error'} = "chrono";
 809                 }
 810             } elsif ($last_time == $ep_time && !length($Dat{'error'})) {
 811                 warn(sprintf(
 812                     "%s: $curr_file: \"%sZ\": Duplicated time\n",
 813                     $progname, sec_to_string($last_time, "T")
 814                 ));
 815                 # FIXME: Make --fix work with gpx.
 816                 if ($Opt{'fix'} && ($Opt{'output-format'} !~ /^gpx$/)) {
 817                     $Dat{'error'} = "duptime";
 818                 }
 819             }
 820         }
 821     } else {
 822         $ep_time = 0;
 823         $Dat{'year'} = 0;
 824         $Dat{'month'} = 0;
 825         $Dat{'day'} = 0;
 826         $Dat{'hour'} = 0;
 827         $Dat{'min'} = 0;
 828         $Dat{'sec'} = 0;
 829     }
 830
 831     if ($Opt{'save-to-file'} ne "\n") {
 832         # {{{
 833         $print_time || return;
 834         my $base_name = "$Dat{'year'}$Dat{'month'}$Dat{'day'}T" .
 835                         "$Dat{'hour'}$Dat{'min'}$Dat{'sec'}Z" .
 836                         "$Opt{'save-to-file'}";
 837         my $file_name = $base_name;
 838         if (-e $file_name) {
 839             for (my $a = 1; (-e $file_name) && ($a < 1000); $a++) {
 840                 $file_name = "$base_name.dup_$a";
 841             }
 842             if (-e $file_name) {
 843                 die("$progname: $base_name: File already exists, and ran " .
 844                     "out of attempts to create unique file name\n");
 845             }
 846             if ($Opt{'verbose'}) {
 847                 warn("$progname: $base_name: File already exists, using " .
 848                      "unique name \"$file_name\" instead\n");
 849             }
 850         }
 851         if (open(my $to_fp, ">", $file_name)) {
 852             print_header(*$to_fp);
 853             print($to_fp (
 854                 $from_stdin
 855                     ? @first_lines
 856                     : ()),
 857                 (length($xml_data)
 858                     ? $xml_data
 859                     : <>)
 860             ) || die("$progname: $file_name: Cannot write to file: $!\n");
 861             print_footer(*$to_fp);
 862             close($to_fp);
 863             if ($Opt{'output-format'} eq "gpsml") {
 864                 printf("<include>%s</include>\n",
 865                        txt_to_xml($file_name));
 866             } elsif ($Opt{'output-format'} eq "gpx") {
 867                 printf("<!-- Saved unconverted data to \"%s\" -->\n",
 868                        txt_to_xml($file_name));
 869             } else {
 870                 print("$progname: Saved unconverted data to \"$file_name\"\n");
 871             }
 872             exit 0;
 873         } else {
 874             die("$progname: $file_name: Cannot create file: $!\n");
 875         }
 876         # }}}
 877     }
 878
 879     my $pause_len = 0;
 880     my $do_print = 1;
 881
 882     if ($Dat{'what'} eq "tp") {
 883         # {{{
 884         if ($Opt{'require'}) {
 885             $Req{'time'} && !$print_time && return;
 886             $Req{'position'} && !$print_pos && return;
 887             $Req{'ele'} && !$print_ele && return;
 888         }
 889
 890         if ($Opt{'inside'} || $Opt{'outside'}) {
 891             if (
 892                 ($Dat{'lat'} < $lat1) ||
 893                 ($Dat{'lat'} > $lat2) ||
 894                 ($Dat{'lon'} < $lon1) ||
 895                 ($Dat{'lon'} > $lon2)
 896             ) {
 897                 $Opt{'inside'} && return;
 898             } else {
 899                 $Opt{'outside'} && return;
 900             }
 901         }
 902
 903         if ($Opt{'output-format'} eq "ps") {
 904             $Dat{'lon'} *= 100;
 905             $Dat{'lat'} *= 100;
 906         }
 907
 908         if (
 909             $Opt{'skip-dups'}
 910             && ($Dat{'lon'} eq $last_lon)
 911             && ($Dat{'lat'} eq $last_lat)
 912         ) {
 913             if ($Opt{'output-format'} eq 'gpsml') {
 914                 $Dat{'error'} = "dup";
 915             } else {
 916                 $do_print = 0;
 917             }
 918         } else {
 919                 $do_print = 1;
 920         }
 921
 922         if (
 923             $Opt{'create-breaks'}
 924             && $ep_time-$last_time > $PAUSE_LIMIT
 925             && $last_time
 926         ) {
 927             $pause_len = $ep_time-$last_time;
 928             # D("pause_len set to '$pause_len'");
 929         }
 930
 931         $Line .= pause_entry($pause_len, $ep_time, $last_time);
 932         # }}}
 933     }
 934
 935     if ($do_print) {
 936         # Valid data was found, send to stdout {{{
 937         unless ($first_time) {
 938             $first_time = $ep_time;
 939         }
 940         $Line .= gen_entry($print_pos, $pause_len, $print_time, $ep_time, $print_ele, %Dat);
 941         # }}}
 942     }
 943
 944     if (!$last_time && $Opt{'output-format'} eq "ps") {
 945         $Line .= "$Dat{'lon'} $Dat{'lat'} m\n";
 946     }
 947
 948     if ($do_print) {
 949         if ($found_move) {
 950             if ($Opt{'output-format'} eq "gpsml") {
 951                 $Line = "<break/>\n$Line";
 952             }
 953             (!$pause_len && ($Opt{'output-format'} eq "xgraph"))
 954                 && ($Line .= "move $Line");
 955             ($Opt{'output-format'} eq "clean") && ($Line .= "\n");
 956             if ($Opt{'output-format'} eq "gpx") {
 957                 $Line .= "$Spc$Spc$Spc$Spc</trkseg>\n" .
 958                          "$Spc$Spc$Spc$Spc<trkseg>\n";
 959             }
 960             $found_move = 0;
 961         }
 962         print($Line);
 963     }
 964     $print_time && ($last_time = $ep_time);
 965     if ($print_pos) {
 966         $last_lon  = $Dat{'lon'};
 967         $last_lat  = $Dat{'lat'};
 968     }
 969     $last_line = $data_line;
 970     $svg_start_thing = "\"/>\n";
 971     # }}}
 972 }
 973
 974 sub gen_entry {
 975     # Generate trackpoint entry, calls trackpoint() {{{
 976     my ($print_pos, $pause_len, $print_time, $ep_time, $print_ele, %Dat) = @_;
 977     my $Line = "";
 978     if ($Opt{'double-y-scale'} && length($Dat{'lat'})) {
 979         $Dat{'lat'} *= 2;
 980     }
 981     if ($Opt{'output-format'} eq "gpsml") {
 982         if ($Dat{'what'} eq "tp") {
 983             $Dat{'format'} = "gpsml";
 984             $Line .= trackpoint(%Dat);
 985         } elsif ($Dat{'what'} =~ /^(pause|desc|title)$/) {
 986             $Line .= sprintf("<%s>%s</%s>\n",
 987                              $1,
 988                              $Dat{$1},
 989                              $1);
 990         }
 991     } elsif ($Opt{'output-format'} eq "pgtab") {
 992         if ($Dat{'what'} eq "tp" && !length($Dat{'error'})) {
 993             $Dat{'format'} = "pgtab";
 994             $Line .= trackpoint(%Dat);
 995         }
 996     } elsif ($Opt{'output-format'} eq "xgraph") {
 997         if ($print_pos && !length($Dat{'error'})) {
 998             $Dat{'format'} = "xgraph";
 999             $Line .= trackpoint(%Dat);
1000         }
1001     } elsif($Opt{'output-format'} eq "gpstrans") {
1002         if ($print_pos && !length($Dat{'error'})) {
1003             $Dat{'format'} = "gpstrans";
1004             $Line .= trackpoint(%Dat);
1005         }
1006     } elsif($Opt{'output-format'} eq "gpx") {
1007         if ($Dat{'what'} eq "tp") {
1008             $Dat{'format'} = "gpx";
1009             $Line .= trackpoint(%Dat);
1010         }
1011     } elsif ($Opt{'output-format'} eq "clean") {
1012         if ($Dat{'what'} eq "tp" && !length($Dat{'error'})) {
1013             $Dat{'format'} = "clean";
1014             $Line .= trackpoint(%Dat);
1015         }
1016     } elsif ($Opt{'output-format'} eq "ps") {
1017         $Line .= (
1018             $pause_len
1019                 ? "f\n$Dat{'lon'} $Dat{'lat'} m\n"
1020                 : "$Dat{'lon'} $Dat{'lat'} l\n"
1021         );
1022     } elsif ($Opt{'output-format'} eq "svg") {
1023         $Line .= (
1024             ($last_lon == 1000) || $pause_len
1025                 ? join("",
1026                     "$svg_start_thing<path\n",
1027                     "  stroke=\"blue\"\n",
1028                     "  stroke-width=\"0.001\"\n",
1029                     "  fill=\"none\"\n",
1030                     "  d=\"\n",
1031                     "M $Dat{'lon'} $Dat{'lat'}\n")
1032                 : "L $Dat{'lon'} $Dat{'lat'}\n"
1033         );
1034     } elsif ($Opt{'output-format'} eq "ygraph") {
1035         if (!length($Dat{'error'})) {
1036             $Dat{'lat'} = num_expand($Dat{'lat'});
1037             $Dat{'lon'} = num_expand($Dat{'lon'});
1038             $Dat{'ele'} = num_expand($Dat{'ele'});
1039             my $Time = $print_time ? ($ep_time - $first_time) * 1 : 0;
1040             $Line .= "\"Time = $Time.0\n$Dat{'lon'} $Dat{'lat'}\n\n";
1041         }
1042     } elsif ($Opt{'output-format'} eq "csv") {
1043         # {{{
1044         if (!length($Dat{'error'})) {
1045             $Dat{'format'} = "csv";
1046             $Dat{'lat'} = num_expand($Dat{'lat'});
1047             $Dat{'lon'} = num_expand($Dat{'lon'});
1048             $Dat{'ele'} = num_expand($Dat{'ele'});
1049             $Line .= join("\t",
1050                 $print_time
1051                     ? $Opt{'epoch'}
1052                         ? $ep_time
1053                         : $Opt{'short-date'}
1054                             ? "$Dat{'year'}$Dat{'month'}$Dat{'day'}T" .
1055                               "$Dat{'hour'}$Dat{'min'}$Dat{'sec'}Z"
1056                             : "$Dat{'year'}-$Dat{'month'}-$Dat{'day'}T" .
1057                               "$Dat{'hour'}:$Dat{'min'}:$Dat{'sec'}Z"
1058                     : "",
1059                 $Dat{'lon'},
1060                 $Dat{'lat'},
1061                 $print_ele ? $Dat{'ele'} : "", # Elevation
1062                 "\n"
1063             );
1064         }
1065         # }}}
1066     } elsif ($Opt{'output-format'} eq "pgwtab") {
1067         # FIXME: NOP at the moment.
1068     } else {
1069         die("$progname: \"$Opt{'output-format'}\": " .
1070             "Unknown output format\n");
1071     }
1072     return($Line);
1073     # }}}
1074 } # gen_entry()
1075
1076 sub pause_entry {
1077     # {{{
1078     my ($pause_len, $ep_time, $last_time) = @_;
1079     my $Line = "";
1080     if ($pause_len) {
1081         if ($Opt{'output-format'} eq "gpsml") {
1082             $Line .= sprintf("<pause>%s</pause>\n",
1083                              sec_to_readable($ep_time-$last_time));
1084         } elsif ($Opt{'output-format'} eq "clean") {
1085             $pause_len && ($Line .= "\n");
1086         } elsif ($Opt{'output-format'} eq "csv") {
1087             $Line .= sprintf("# Pause: %s\n# move\n",
1088                              sec_to_readable($ep_time-$last_time));
1089         } elsif ($Opt{'output-format'} eq "xgraph") {
1090             $pause_len && ($Line .= "move ");
1091         }
1092     }
1093     return($Line);
1094     # }}}
1095 } # pause_entry()
1096
1097 sub ps_header {
1098     # Send a Postscript header to stdout {{{
1099     my ($bl_lon, $bl_lat, $br_lon, $br_lat) = @_;
1100     my $Date = sec_to_string(time);
1101     return(join("",
1102         "%!PS-Adobe-3.0 EPSF-3.0\n",
1103         "%%Creator: gpst\n",
1104         "%%Title:\n",
1105         "%%CreationDate: $Date\n",
1106         "%%BoundingBox: $bl_lon $bl_lat $br_lon $br_lat\n",
1107         "%%DocumentData: Clean7Bit\n",
1108         "%%EndComments\n",
1109         "%%BeginProlog\n",
1110         "/bd { bind def } bind def\n",
1111         "/incompound false def\n",
1112         "/m { moveto } bd\n",
1113         "/l { lineto } bd\n",
1114         "/c { curveto } bd\n",
1115         "/F { incompound not {fill} if } bd\n",
1116         "/f { closepath F } bd\n",
1117         "/S { stroke } bd\n",
1118         "/*u { /incompound true def } bd\n",
1119         "/*U { /incompound false def f} bd\n",
1120         "/k { setcmykcolor } bd\n",
1121         "/K { k } bd\n",
1122         "%%EndProlog\n",
1123         "%%BeginSetup\n",
1124         "%%EndSetup\n",
1125     ));
1126     # }}}
1127 }
1128
1129 sub print_version {
1130     # Print program version {{{
1131     print("$progname v$VERSION\n");
1132     # }}}
1133 } # print_version()
1134
1135 sub usage {
1136     # Send the help message to stdout {{{
1137     my $Retval = shift;
1138
1139     if ($Opt{'verbose'}) {
1140         print("\n");
1141         print_version();
1142     }
1143     print(<<END);
1144
1145 Converts between various GPS formats.
1146
1147 Usage: $progname [options] [file [files [...]]]
1148        $progname -S [file [files [...]]]
1149        $progname -u [file [files [...]]]
1150
1151 Options:
1152
1153   --chronology
1154     Check for broken chronology, warn about entries with an old
1155     timestamp.
1156   -d, --skip-dups
1157     Skip duplicated coordinates.
1158   -e, --epoch
1159     Use seconds since 1970-01-01 00:00:00 GMT as date format.
1160   --fix
1161     Comment out entries which is obviously wrong. Use together with
1162     --chronology to fix those kind of errors. Does not work with GPX
1163     output yet.
1164   --from-date x
1165     Used by the pgwupd format. Specifies from which date waypoints
1166     should be updated. No checks for valid date format here, let
1167     PostgreSQL take care of that. All variants it understands can be
1168     used here.
1169   -h, --help
1170     Show this help.
1171   --inside
1172     Print only trackpoints inside a rectangle specified by --pos1 and
1173     --pos2.
1174   -n, --undefined x
1175     Use x as undefined value. Default: "$Udef".
1176   -o, --output-format x
1177     Use output format x:
1178       clean
1179       csv
1180       gpsml (Default)
1181       gpstrans
1182       gpx (Not complete)
1183       pgtab
1184       pgwtab
1185       pgwupd
1186       poscount
1187       ps (Unfinished)
1188       svg (Unfinished)
1189       xgraph
1190       ygraph
1191   --outside
1192     Print only trackpoints outside a rectangle specified by --pos1 and
1193     --pos2.
1194   --pos1 x
1195   --pos2 x
1196     Specifies one corner where x is in "lat,lon" format (decimal
1197     degrees, negative for west or south) of area rectangle used by the
1198     --inside and --outside options.
1199   -r, --require x
1200     Specify requirements for trackpoints to be written. x is a string
1201     with the following flags:
1202       e
1203         Print only waypoints which have an elevation.
1204       p
1205         Print only waypoints which have a position.
1206       t
1207         Print only waypoints which have a timestamp.
1208   -R, --round x=y[,x2=y2[...]]
1209     Round trackpoint element x to y decimals. Example:
1210       --round lat=4,lon=5,ele=1
1211   -s, --short-date
1212     Use short date format.
1213   -S, --save-to-file x
1214     Save the unconverted data to a file with a filename starting with
1215     the timestamp of the first trackpoint. The parameter string x is
1216     added at the end of the filename. For the time being this option
1217     will ignore all other options. Note: If several files are specified
1218     on the command line, all data will be saved into only one file. This
1219     behaviour may change in the future.
1220   -t, --create-breaks
1221     Create breaks in track between points with a difference more than
1222     $PAUSE_LIMIT seconds.
1223   -T x, --time-shift x
1224     Move time of trackpoint x seconds forwards or backwards. x can be a
1225     positive or negative integer.
1226   -v, --verbose
1227     Increase level of verbosity. Can be repeated.
1228   --version
1229     Print version information.
1230   -w, --strip-whitespace
1231     Strip all unnecessary whitespace.
1232   -y, --double-y-scale
1233     Double Y scale (latitude) to get it right in gnuplot.
1234   --debug
1235     Print debugging messages.
1236
1237 END
1238     exit($Retval);
1239     # }}}
1240 } # usage()
1241
1242 sub msg {
1243     # Print a status message to stderr based on verbosity level {{{
1244     my ($verbose_level, $Txt) = @_;
1245
1246     if ($Opt{'verbose'} >= $verbose_level) {
1247         print(STDERR "$progname: $Txt\n");
1248     }
1249     # }}}
1250 } # msg()
1251
1252 __END__
1253
1254 # Law talk {{{
1255 # Copyleft © Øyvind A. Holm <sunny@sunbase.org>
1256 #
1257 This program is free software: you can redistribute it and/or modify it
1258 under the terms of the GNU General Public License as published by the
1259 Free Software Foundation, either version 3 of the License, or (at your
1260 option) any later version.
1261
1262 This program is distributed in the hope that it will be useful, but
1263 WITHOUT ANY WARRANTY; without even the implied warranty of
1264 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
1265 See the GNU General Public License for more details.
1266
1267 You should have received a copy of the GNU General Public License along
1268 with this program.
1269 If not, see L<http://www.gnu.org/licenses/>.
1270 # }}}
1271
1272 # vim: set fenc=UTF-8 ft=perl fdm=marker ts=4 sw=4 sts=4 et fo+=w :