3 #=======================================================================
5 # File ID: 082106ec-f924-11dd-b757-0001805bf4b1
6 # Converts between various GPS formats
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 #=======================================================================
16 use Time
::Local qw
{ timegm_nocheck
};
19 push(@INC, "$ENV{'HOME'}/bin/src/gpstools");
33 # Initial values for command line arguments {{{
38 'double-y-scale' => 0,
44 'output-format' => "gpsml",
50 'save-to-file' => "\n", # \n = undefined, it’s banned in filenames anyway.
53 'strip-whitespace' => 0,
63 $progname =~ s/^.*\/(.*?)$/$1/;
64 our $VERSION = "0.00";
66 Getopt
::Long
::Configure
("bundling");
68 # Command line options {{{
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'},
95 ) || die("$progname: Option error. Use -h for help.\n");
99 my $PAUSE_LIMIT = 2 * 60; # Antall sekunder mellom to punkter det må til før en move legges inn.
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.
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);
113 'gpsbabel' => '/usr/local/bin/gpsbabel',
118 if ($Opt{'output-format'} =~ /^(gpx|pgtab)$/) {
119 $Opt{'require'} .= "p";
122 'ele' => ($Opt{'require'} =~ /e/) ?
1 : 0,
123 'position' => ($Opt{'require'} =~ /p/) ?
1 : 0,
124 'time' => ($Opt{'require'} =~ /t/) ?
1 : 0,
126 $Opt{'require'} =~ /[^ept]/
127 && die("$0: Unknown flag in --require (-r) value\n");
129 $Opt{'debug'} && ($Debug = 1);
130 $Opt{'help'} && usage
(0);
131 if ($Opt{'version'}) {
136 if ($Opt{'pos1'} =~ /^($DIGIT+),($DIGIT+)$/) {
140 if ($Opt{'pos2'} =~ /^($DIGIT+),($DIGIT+)$/) {
155 if ($Opt{'epoch'} && $Opt{'short-date'}) {
156 die("$progname: Cannot mix the --epoch (-e) and --short-date (-s) options\n");
159 if ($Opt{'inside'} && $Opt{'outside'}) {
160 die("$progname: Cannot mix the --inside and --outside options\n");
163 # To avoid printing out extra "/> at the start of svg output:
164 my $svg_start_thing = "";
168 if (defined($Opt{'round'})) {
169 my $R = $Opt{'round'};
170 $R =~ s/([a-z]+)=(\d+)/($Round{$1}=$2, "")/eg;
173 length($Opt{'undefined'}) && ($Udef = $Opt{'undefined'});
175 $Opt{'save-to-file'} eq "\n" && print_header
(*STDOUT
);
182 my $from_stdin = scalar(@ARGV) ?
0 : 1;
184 $from_stdin && push(@ARGV, "-");
186 for $curr_file (@ARGV) {
187 # Scan through stdin or specified files and send every GPS entry to
190 print(STDERR
"$progname: Opening \"$curr_file\" for read\n") if $Opt{'verbose'};
191 if (open(CurrFP
, "<$curr_file")) {
196 'year' => '', 'month' => '', 'day' => '',
197 'hour' => '', 'min' => '', 'sec' => '',
200 'lat' => '', 'lon' => '',
207 $Opt{'epoch'} && ($Dat{'date-format'} = "epoch");
208 $Opt{'short-date'} && ($Dat{'date-format'} = "short");
210 if ($Opt{'save-to-file'} ne "\n") {
211 push(@first_lines, $_);
213 s/^# error // && ($Dat{'error'} = "error");
214 s/^# ?// && ($Dat{'error'} = "desc");
216 if (m
#^<(e?tp)\b(.*?)>(.*?)</(e?tp)>\s*$#) {
217 # gpsml — The main storage format {{{
218 my ($Elem, $Props, $Data) =
220 my $err_str = ($Props =~ /\berr="(.*?)"/) ?
$1 : "error";
221 $Elem eq "etp" && ($Dat{'error'} = $err_str);
223 $Data =~ m
#<time>(.*?)</time># && ($Time = $1);
225 (\d\d\d\d
)-?
(\d\d
)-?
(\d\d
)[T
](\d\d
):?
(\d\d
):?
([\d\
.]+?
)Z
227 ($Dat{'year'}, $Dat{'month'}, $Dat{'day'},
228 $Dat{'hour'}, $Dat{'min'}, $Dat{'sec'}) =
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);
239 } elsif (m
#^<break\b.*?/>#) {
241 } elsif (m
#^<(title|pause)\b.*?>(.*?)</(title|pause)>#) {
245 } elsif (m
#^<desc\b.*?>(.*$)#s) {
246 $Dat{'what'} = "desc";
248 until ($Txt =~ m
#</desc>#s) {
251 $Txt =~ s
#^(.*)(</desc>.*$)#$1#s;
256 $xml_data .= join("", <CurrFP
>);
257 if (!length($Opt{'output-format'})) {
258 $Opt{'output-format'} = "gpx";
259 print_header
(*STDOUT
);
261 read_xmlfile
($xml_data);
265 } elsif (m
#^(\d+)\t($DIGIT+)\t($DIGIT+)\t($DIGIT)#) {
266 # CSV format, epoch style {{{
267 my ($ep_time, $lon_val, $lat_val, $Alt) =
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 Ⅱ
279 (\d\d\d\d
)-?
(\d\d
)-?
(\d\d
)[T\
](\d\d
):?
(\d\d
):?
(\d\d
)Z?
\t
280 ($DIGIT+)\t($DIGIT+)\t($DIGIT)
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'}) =
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 {{{
296 # N60.41630 E5.31675\t
297 # 09.02.2006 20:24:37 (UTC)\t
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) =
309 # Nødløsning for å unngå at variabler
311 "\t\t\t\t\t\t\t\t\t\t"
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",
326 $Legtime_hour, $Legtime_min, $Legtime_sec,
327 $Legspeed, $Legspeed_unit,
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.
343 defined($Dat{'year'})
344 && $Dat{'year'} =~ /\d/
345 && $Dat{'year'} < 1900
346 ) && ($Dat{'year'} += 2000);
349 } elsif (/^Track\t(.*?)\t/) {
350 $Dat{'title'} = txt_to_xml
($1);
351 $Dat{'what'} = "title";
357 (\d\d
)/(\d\d)/(\d\d\d\d
)\
(\d\d
):(\d\d
):(\d\d
)\t
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) =
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));
382 1\ (\S+)\ (\S+)\ (\S+)\ (\S+)\x20
383 (\d\d)/(\d\d)/(\d\d\d\d)\ (\d\d):(\d\d):(\d\d)
386 # 1 60.3938222 5.3238754 17.3 0 09/01/2002 14:18:23 {{{
387 ($Dat{'lat
'}, $Dat{'lon
'}, $Dat{'speed
'},
389 $Dat{'month
'}, $Dat{'day
'}, $Dat{'year
'},
390 $Dat{'hour
'}, $Dat{'min
'}, $Dat{'sec
'}) =
398 # @020721221336N6048353E00701826S015-00001E4859N1673U0000 {{{
408 (\d\d) # Latitude degree
409 (\d\d) # Latitude minute
410 (\d\d\d) # Latitude minute decimals
412 (\d\d\d) # Longitude degree
413 (\d\d) # Longitude minute
414 (\d\d\d) # Longitude minute degree
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,
432 my $ep_time = timegm_nocheck(
433 $Dat{'sec
'}, $Dat{'min
'}, $Dat{'hour
'},
434 $Dat{'day
'}, $Dat{'month
'}-1, $Dat{'year
'}
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(
444 my $tmp_lat = sprintf("%${Flat}f",
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;
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
'}) =
460 $Dat{'error
'} = "nosignal";
463 } elsif (/^xmaplog /) {
466 ($Opt{'output
-format
'} eq "csv")
467 && ($Opt{'save
-to
-file
'} eq "\n")
469 } elsif (/^Pause: /) {
470 # NOP, is here to cope with old files I’ve lying around.
471 } elsif ($Dat{'error
'} eq "desc") {
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";
482 $Opt{'verbose
'} && warn("Line $.: Unknown: \"$_\"\n");
487 warn("$progname: $curr_file: Cannot open file for read: $!\n");
492 print_footer(*STDOUT);
498 my $Txt = join("", @_);
499 $Txt =~ s/<!--(.*?)-->//gs;
500 $Txt =~ s#(<gpx\b.*?>.*?</gpx>)#print_gpx($1)#gse;
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 ";
514 if ($Opt{'output
-format
'} =~ /^(pgwtab|pgwupd)$/) {
518 <wpt\b(.*?)>(.*?)</wpt>
523 my ($Lat, $Lon, $Name, $Ele, $Type, $Time, $Cmt, $Desc, $Src, $Sym) =
524 ('\N
', '\N
', '\N
', '\N
', '\N
', '\N
', '\N
', '\N
', '\N
', '\N
');
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)));
547 if (length($Opt{'round
'})) {
548 if (defined($Round{'lat
'}) && length($Lat)) {
549 ($Lat = 1.0 * sprintf("%.$Round{'lat
'}f", $Lat));
551 if (defined($Round{'lon
'}) && length($Lon)) {
552 ($Lon = 1.0 * sprintf("%.$Round{'lon
'}f", $Lon));
554 if (defined($Round{'ele
'}) && $Ele ne '\N
') {
555 ($Ele = 1.0 * sprintf("%.$Round{'ele
'}f", $Ele));
559 if ($Opt{'output
-format
'} eq "pgwtab") {
573 } elsif ($Opt{'output
-format
'} eq "pgwupd") {
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;",
591 <trk\b(.*?)>(.*?)</trk>
597 <name\b(.*?)>(.*?)</name>
600 $tmp_dat{'title
'} = $2;
601 $tmp_dat{'what
'} = "title";
602 $tmp_dat{'error
'} = "";
603 print_entry(%tmp_dat);
608 <trkseg\b(.*?)>(.*?)</trkseg>
614 <trkpt\b(.*?)>(.*?)</trkpt>
617 my ($attr_trkpt, $el_trkpt) =
620 'year
' => '', 'month
' => '', 'day
' => '',
621 'hour
' => '', 'min
' => '', 'sec
' => '',
624 'lat
' => '', 'lon
' => '',
630 ($attr_trkpt =~ /\blon="(.*?)"/) && ($Dat{'lon
'} = $1);
631 ($attr_trkpt =~ /\blat="(.*?)"/) && ($Dat{'lat
'} = $1);
632 ($el_trkpt =~ m#<ele\b.*?>(.*?)</ele>#) && ($Dat{'ele
'} = $1);
636 <time>(\d\d\d\d)-?(\d\d)-?(\d\d)T
637 (\d\d):?(\d\d):?([\d\.]+)Z</time>
640 ($Dat{'year
'}, $Dat{'month
'}, $Dat{'day
'},
641 $Dat{'hour
'}, $Dat{'min
'}, $Dat{'sec
'}) =
642 ($1, $2, $3, $4, $5, $6);
658 local *OutFP = shift;
659 if ($Opt{'output
-format
'} eq "gpsml") {
661 "<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n",
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") {
670 qq{<?xml version="1.0" encoding="UTF-8" standalone="no"?>\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},
680 qq{$Spc$Spc$Spc$Spc<trkseg>\n},
682 } elsif ($Opt{'output
-format
'} eq "ps") {
683 print(OutFP ps_header(532, 6034, 533, 6040));
685 } elsif ($Opt{'output
-format
'} eq "svg") {
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",
701 local *OutFP = shift;
702 if ($Opt{'output
-format
'} eq "gpsml") {
707 } elsif ($Opt{'output
-format
'} eq "gpx") {
709 "$Spc$Spc$Spc$Spc</trkseg>\n",
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");
718 } elsif ($Opt{'output
-format
'} eq "ps") {
724 } elsif ($Opt{'output
-format
'} eq "svg") {
725 print(OutFP "\"/>\n</svg>\n");
731 # Print a GPS entry with time, latitude, longitude and elevation in
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;
742 if (!$Req{'position
'} && $Opt{'output
-format
'} eq "gpsml") {
743 $print_pos = (length($Dat{'lat
'}) || length($Dat{'lon
'})) ? 1 : 0;
745 $print_pos = (length($Dat{'lat
'}) && length($Dat{'lon
'})) ? 1 : 0;
748 $Dat{'lat
'} = $Dat{'lon
'} = "";
750 my $print_ele = length($Dat{'ele
'}) ? 1 : 0;
751 my $print_desc = length($Dat{'desc
'}) ? 1 : 0;
753 # D("print_entry(\"" . join("\", \"", @_) . "\");");
756 if (length($Opt{'round
'})) {
757 for my $Tmp (qw{ lat lon ele }) {
758 if (defined($Round{$Tmp}) && length($Dat{$Tmp})) {
760 ($Dat{$Tmp} = 1.0 * sprintf("%.$Round{$Tmp}f", $Dat{$Tmp}));
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);
775 $ep_time = timegm_nocheck(
776 $Dat{'sec
'}, $Dat{'min
'}, $Dat{'hour
'},
777 $Dat{'day
'}, $Dat{'month
'} - 1, $Dat{'year
'}
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;
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
'})) {
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")
803 # FIXME: Make --fix work with gpx.
804 if ($Opt{'fix
'} && ($Opt{'output
-format
'} !~ /^gpx$/)) {
805 $Dat{'error
'} = "chrono";
807 } elsif ($last_time == $ep_time && !length($Dat{'error
'})) {
809 "%s: $curr_file: \"%sZ\": Duplicated time\n",
810 $progname, sec_to_string($last_time, "T")
812 # FIXME: Make --fix work with gpx.
813 if ($Opt{'fix
'} && ($Opt{'output
-format
'} !~ /^gpx$/)) {
814 $Dat{'error
'} = "duptime";
828 if ($Opt{'save
-to
-file
'} ne "\n") {
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;
836 for (my $a = 1; (-e $file_name) && ($a < 1000); $a++) {
837 $file_name = "$base_name.dup_$a";
840 die("$progname: $base_name: File already exists, and ran " .
841 "out of attempts to create unique file name\n");
843 if ($Opt{'verbose
'}) {
844 warn("$progname: $base_name: File already exists, using " .
845 "unique name \"$file_name\" instead\n");
848 if (open(ToFP, ">", $file_name)) {
857 ) || die("$progname: $file_name: Cannot write to file: $!\n");
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));
867 print("$progname: Saved unconverted data to \"$file_name\"\n");
871 die("$progname: $file_name: Cannot create file: $!\n");
879 if ($Dat{'what
'} eq "tp") {
881 if ($Opt{'require'}) {
882 $Req{'time'} && !$print_time && return;
883 $Req{'position
'} && !$print_pos && return;
884 $Req{'ele
'} && !$print_ele && return;
887 if ($Opt{'inside
'} || $Opt{'outside
'}) {
889 ($Dat{'lat
'} < $lat1) ||
890 ($Dat{'lat
'} > $lat2) ||
891 ($Dat{'lon
'} < $lon1) ||
892 ($Dat{'lon
'} > $lon2)
894 $Opt{'inside
'} && return;
896 $Opt{'outside
'} && return;
900 if ($Opt{'output
-format
'} eq "ps") {
907 && ($Dat{'lon
'} eq $last_lon)
908 && ($Dat{'lat
'} eq $last_lat)
910 if ($Opt{'output
-format
'} eq 'gpsml
') {
911 $Dat{'error
'} = "dup";
920 $Opt{'create
-breaks
'}
921 && $ep_time-$last_time > $PAUSE_LIMIT
924 $pause_len = $ep_time-$last_time;
925 # D("pause_len set to '$pause_len'");
928 $Line .= pause_entry($pause_len, $ep_time, $last_time);
933 # Valid data was found, send to stdout {{{
934 unless ($first_time) {
935 $first_time = $ep_time;
937 $Line .= gen_entry($print_pos, $pause_len, $print_time, $ep_time, $print_ele, %Dat);
941 if (!$last_time && $Opt{'output
-format
'} eq "ps") {
942 $Line .= "$Dat{'lon
'} $Dat{'lat
'} m\n";
947 if ($Opt{'output
-format
'} eq "gpsml") {
948 $Line = "<break/>\n$Line";
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";
961 $print_time && ($last_time = $ep_time);
963 $last_lon = $Dat{'lon
'};
964 $last_lat = $Dat{'lat
'};
966 $last_line = $data_line;
967 $svg_start_thing = "\"/>\n";
972 # Generate trackpoint entry, calls trackpoint() {{{
973 my ($print_pos, $pause_len, $print_time, $ep_time, $print_ele, %Dat) = @_;
975 if ($Opt{'double
-y
-scale
'} && length($Dat{'lat
'})) {
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",
988 } elsif ($Opt{'output
-format
'} eq "pgtab") {
989 if ($Dat{'what
'} eq "tp" && !length($Dat{'error
'})) {
990 $Dat{'format
'} = "pgtab";
991 $Line .= trackpoint(%Dat);
993 } elsif ($Opt{'output
-format
'} eq "xgraph") {
994 if ($print_pos && !length($Dat{'error
'})) {
995 $Dat{'format
'} = "xgraph";
996 $Line .= trackpoint(%Dat);
998 } elsif($Opt{'output
-format
'} eq "gpstrans") {
999 if ($print_pos && !length($Dat{'error
'})) {
1000 $Dat{'format
'} = "gpstrans";
1001 $Line .= trackpoint(%Dat);
1003 } elsif($Opt{'output
-format
'} eq "gpx") {
1004 if ($Dat{'what
'} eq "tp") {
1005 $Dat{'format
'} = "gpx";
1006 $Line .= trackpoint(%Dat);
1008 } elsif ($Opt{'output
-format
'} eq "clean") {
1009 if ($Dat{'what
'} eq "tp" && !length($Dat{'error
'})) {
1010 $Dat{'format
'} = "clean";
1011 $Line .= trackpoint(%Dat);
1013 } elsif ($Opt{'output
-format
'} eq "ps") {
1016 ? "f\n$Dat{'lon
'} $Dat{'lat
'} m\n"
1017 : "$Dat{'lon
'} $Dat{'lat
'} l\n"
1019 } elsif ($Opt{'output
-format
'} eq "svg") {
1021 ($last_lon == 1000) || $pause_len
1023 "$svg_start_thing<path\n",
1024 " stroke=\"blue\"\n",
1025 " stroke-width=\"0.001\"\n",
1028 "M $Dat{'lon
'} $Dat{'lat
'}\n")
1029 : "L $Dat{'lon
'} $Dat{'lat
'}\n"
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";
1036 } elsif ($Opt{'output
-format
'} eq "csv") {
1038 if (!length($Dat{'error
'})) {
1039 $Dat{'format
'} = "csv";
1044 : $Opt{'short
-date
'}
1045 ? "$Dat{'year
'}$Dat{'month
'}$Dat{'day
'}T" .
1046 "$Dat{'hour
'}$Dat{'min
'}$Dat{'sec
'}Z"
1047 : "$Dat{'year
'}-$Dat{'month
'}-$Dat{'day
'}T" .
1048 "$Dat{'hour
'}:$Dat{'min
'}:$Dat{'sec
'}Z"
1052 $print_ele ? $Dat{'ele
'} : "", # Elevation
1057 } elsif ($Opt{'output
-format
'} eq "pgwtab") {
1058 # FIXME: NOP at the moment.
1060 die("$progname: \"$Opt{'output
-format
'}\": " .
1061 "Unknown output format\n");
1069 my ($pause_len, $ep_time, $last_time) = @_;
1072 if ($Opt{'output
-format
'} eq "gpsml") {
1073 $Line .= sprintf("<pause>%s</pause>\n",
1074 sec_to_readable($ep_time-$last_time));
1075 } elsif ($Opt{'output
-format
'} eq "clean") {
1076 $pause_len && ($Line .= "\n");
1077 } elsif ($Opt{'output
-format
'} eq "csv") {
1078 $Line .= sprintf("# Pause: %s\n# move\n",
1079 sec_to_readable($ep_time-$last_time));
1080 } elsif ($Opt{'output
-format
'} eq "xgraph") {
1081 $pause_len && ($Line .= "move ");
1089 # Send a Postscript header to stdout {{{
1090 my ($bl_lon, $bl_lat, $br_lon, $br_lat) = @_;
1091 my $Date = sec_to_string(time);
1093 "%!PS-Adobe-3.0 EPSF-3.0\n",
1094 "%%Creator: gpst\n",
1096 "%%CreationDate: $Date\n",
1097 "%%BoundingBox: $bl_lon $bl_lat $br_lon $br_lat\n",
1098 "%%DocumentData: Clean7Bit\n",
1101 "/bd { bind def } bind def\n",
1102 "/incompound false def\n",
1103 "/m { moveto } bd\n",
1104 "/l { lineto } bd\n",
1105 "/c { curveto } bd\n",
1106 "/F { incompound not {fill} if } bd\n",
1107 "/f { closepath F } bd\n",
1108 "/S { stroke } bd\n",
1109 "/*u { /incompound true def } bd\n",
1110 "/*U { /incompound false def f} bd\n",
1111 "/k { setcmykcolor } bd\n",
1121 # Print program version {{{
1122 print("$progname v$VERSION\n");
1127 # Send the help message to stdout {{{
1130 if ($Opt{'verbose
'}) {
1136 Converts between various GPS formats.
1138 Usage: $progname [options] [file [files [...]]]
1139 $progname -S [file [files [...]]]
1140 $progname -u [file [files [...]]]
1145 Check for broken chronology, warn about entries with an old
1148 Skip duplicated coordinates.
1150 Use seconds since 1970-01-01 00:00:00 GMT as date format.
1152 Comment out entries which is obviously wrong. Use together with
1153 --chronology to fix those kind of errors. Does not work with GPX
1156 Used by the pgwupd format. Specifies from which date waypoints
1157 should be updated. No checks for valid date format here, let
1158 PostgreSQL take care of that. All variants it understands can be
1163 Print only trackpoints inside a rectangle specified by --pos1 and
1166 Use x as undefined value. Default: "$Udef".
1167 -o, --output-format x
1168 Use output format x:
1183 Print only trackpoints outside a rectangle specified by --pos1 and
1187 Specifies one corner where x is in "lat,lon" format (decimal
1188 degrees, negative for west or south) of area rectangle used by the
1189 --inside and --outside options.
1191 Specify requirements for trackpoints to be written. x is a string
1192 with the following flags:
1194 Print only waypoints which have an elevation.
1196 Print only waypoints which have a position.
1198 Print only waypoints which have a timestamp.
1199 -R, --round x=y[,x2=y2[...]]
1200 Round trackpoint element x to y decimals. Example:
1201 --round lat=4,lon=5,ele=1
1203 Use short date format.
1204 -S, --save-to-file x
1205 Save the unconverted data to a file with a filename starting with
1206 the timestamp of the first trackpoint. The parameter string x is
1207 added at the end of the filename. For the time being this option
1208 will ignore all other options. Note: If several files are specified
1209 on the command line, all data will be saved into only one file. This
1210 behaviour may change in the future.
1212 Create breaks in track between points with a difference more than
1213 $PAUSE_LIMIT seconds.
1214 -T x, --time-shift x
1215 Move time of trackpoint x seconds forwards or backwards. x can be a
1216 positive or negative integer.
1218 Increase level of verbosity. Can be repeated.
1220 Print version information.
1221 -w, --strip-whitespace
1222 Strip all unnecessary whitespace.
1223 -y, --double-y-scale
1224 Double Y scale (latitude) to get it right in gnuplot.
1226 Print debugging messages.
1234 # Print a status message to stderr based on verbosity level {{{
1235 my ($verbose_level, $Txt) = @_;
1237 if ($Opt{'verbose'} >= $verbose_level) {
1238 print(STDERR
"$progname: $Txt\n");
1246 # Copyleft © Øyvind A. Holm <sunny@sunbase.org>
1248 This program is free software
: you can redistribute it
and/or modify it
1249 under the terms of the GNU General Public License as published by the
1250 Free Software Foundation
, either version
3 of the License
, or (at your
1251 option
) any later version
.
1253 This program is distributed
in the hope that it will be useful
, but
1254 WITHOUT ANY WARRANTY
; without even the implied warranty of
1255 MERCHANTABILITY
or FITNESS FOR A PARTICULAR PURPOSE
.
1256 See the GNU General Public License
for more details
.
1258 You should have received a copy of the GNU General Public License along
1260 If
not, see L
<http
://www
.gnu
.org
/licenses/>.
1263 # vim: set fenc=UTF-8 ft=perl fdm=marker ts=4 sw=4 sts=4 et fo+=w :