Changeset 10388 in ntrip for trunk/BNC/scripts

Timestamp:

Mar 12, 2024, 3:10:17 PM (11 months ago)

Author:

stuerze

Message:

changes regarding PPP

Location:

trunk/BNC/scripts

Files:

• Bnc.pm (modified) (4 diffs)
• pppPlot.pl (modified) (5 diffs)

trunk/BNC/scripts/Bnc.pm

@@ -493 +493 @@
     my $epochSec     = 0;
     my $epochDiff    = 0;
-    my ( @hlp, @N, @E, @U, %SATNUM, @TRP, @CLK, @OFF_GLO, @OFF_GAL, @OFF_BDS);
-    my ( @EPOCH, @EPOCH_CLK, @EPOCH_OFF_GLO, @EPOCH_OFF_GAL, @EPOCH_OFF_BDS );
+    my ( @hlp, @N, @E, @U, %SATNUM, @TRP, @CLK, @OFF_GPS, @OFF_GLO, @OFF_GAL, @OFF_BDS);
+    my ( @EPOCH, @EPOCH_CLK, @EPOCH_OFF_GPS, @EPOCH_OFF_GLO, @EPOCH_OFF_GAL, @EPOCH_OFF_BDS );
     my ( %AMB,        %RES,            %ELE,            %ION, %BIA );
     my ( $station,    $lki,            $sys,            $sat, $amb );
@@ -588 +588 @@
             push ( @CLK,       $hlp[2] + $hlp[3] );
         }
+        elsif ( $ln =~ /\bOFF_GPS\b/ ) {                # 2015-08... OFF_GPS 52.6806 -3.8042 +- 9.0077
+            push ( @EPOCH_OFF_GPS, $epochSec );
+            push ( @OFF_GPS,       $hlp[2] + $hlp[3] );
+        }
         elsif ( $ln =~ /\bOFF_GLO\b/ ) {                # 2015-08... OFF_GLO 52.6806 -3.8042 +- 9.0077
-#            push ( @EPOCH_OFF_GLO, $epochSec );
+            push ( @EPOCH_OFF_GLO, $epochSec );
             push ( @OFF_GLO,       $hlp[2] + $hlp[3] );
         }
         elsif ( $ln =~ /\bOFF_GAL\b/ ) {                # 2015-08... OFF_GAL 52.6806 -3.8042 +- 9.0077
-#            push ( @EPOCH_OFF_GAL, $epochSec );
+            push ( @EPOCH_OFF_GAL, $epochSec );
             push ( @OFF_GAL,       $hlp[2] + $hlp[3] );
         }
         elsif ( $ln =~ /\bOFF_BDS\b/ ) {                # 2015-08... OFF_BDS 52.6806 -3.8042 +- 9.0077
-#            push ( @EPOCH_OFF_BDS, $epochSec );
+            push ( @EPOCH_OFF_BDS, $epochSec );
             push ( @OFF_BDS,       $hlp[2] + $hlp[3] );
         }
@@ -624 +628 @@
     if ( $nof_epochs != scalar @TRP )                            { LOGDIE "number of epochs and TRP not equal\n" }
     if ( @CLK     && scalar @EPOCH_CLK     != scalar @CLK )      { LOGDIE "number of epochs and CLK not equal\n" }
-#    if ( @OFF_GLO && scalar @EPOCH_OFF_GLO != scalar @OFF_GLO )  { LOGDIE "number of epochs and OFF_GLO not equal\n" }
-#    if ( @OFF_GAL && scalar @EPOCH_OFF_GAL != scalar @OFF_GAL )  { LOGDIE "number of epochs and OFF_GAL not equal\n" }
-#    if ( @OFF_BDS && scalar @EPOCH_OFF_BDS != scalar @OFF_BDS )  { LOGDIE "number of epochs and OFF_BDS not equal\n" }
+    if ( @OFF_GPS && scalar @EPOCH_OFF_GPS != scalar @OFF_GPS )  { LOGDIE "number of epochs and OFF_GPS not equal\n" }
+    if ( @OFF_GLO && scalar @EPOCH_OFF_GLO != scalar @OFF_GLO )  { LOGDIE "number of epochs and OFF_GLO not equal\n" }
+    if ( @OFF_GAL && scalar @EPOCH_OFF_GAL != scalar @OFF_GAL )  { LOGDIE "number of epochs and OFF_GAL not equal\n" }
+    if ( @OFF_BDS && scalar @EPOCH_OFF_BDS != scalar @OFF_BDS )  { LOGDIE "number of epochs and OFF_BDS not equal\n" }
 
     if ( !$station ) { WARN "could not grep stationname from file: $file\n" }
@@ -638 +643 @@
                  TRP     => \@TRP,
                  CLK     => \@CLK,
+                 OFF_GPS => \@OFF_GPS,
                  OFF_GLO => \@OFF_GLO,
                  OFF_GAL => \@OFF_GAL,

trunk/BNC/scripts/pppPlot.pl

@@ -107 +107 @@
     my ( $station, $file ) = Bnc::parsePPPLogfile( $file, $sampling );
     my $EPOCH       = $file->{'EPOCH'};
+    my $EPOCH_OFF_GPS = $file->{'EPOCH_OFF_GPS'};
     my $EPOCH_OFF_GLO = $file->{'EPOCH_OFF_GLO'};
     my $EPOCH_OFF_GAL = $file->{'EPOCH_OFF_GAL'};
@@ -243 +244 @@
         }
 
+        ######### OFF_GPS #####################
+        if ( grep ( $_ eq "ALL", @plotTypes ) ) {
+            if ( scalar @{ $file->{'OFF_GPS'} } < 1 ) {
+                DEBUG "No OFF_GPS found";
+            }
+            else {
+                DEBUG "Plot OFF_GPS";
+                my $pngName  = sprintf ( "%s_OFF_G.png", $station );
+                my $chartOFF_GPS = newChart($station);
+                $chartOFF_GPS->set( output => $pngName );
+                $chartOFF_GPS->set( ylabel => "Offset GPS [m]" );
+
+                $dataset = Chart::Gnuplot::DataSet->new(
+                                                         xdata   => $EPOCH_OFF_GPS,
+                                                         ydata   => $file->{'OFF_GPS'},
+                                                         title   => "Receiver Offset GPS",
+                                                         timefmt => '%s',
+                                                         style   => "linespoints",
+                );
+                $chartOFF_GPS->plot2d($dataset);    #system ("display $pngName&");
+                $y = $y - $dy;
+
+                if ( $y < 30 / mm ) {
+                    $page = $pdf->page();
+                    $page->mediabox('A4');
+                    $y = $y0;
+                }
+                $png = $page->gfx();
+                LOGDIE("could not find image file: $!\n") unless -e $pngName;
+                $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
+            }
+        }
+
 
         ######### OFF_GLO #####################
@@ -257 +291 @@
 
                 $dataset = Chart::Gnuplot::DataSet->new(
-                                                         #xdata   => $EPOCH_OFF_GLO,
-                                                         xdata   => $EPOCH,
+                                                         xdata   => $EPOCH_OFF_GLO,
                                                          ydata   => $file->{'OFF_GLO'},
                                                          title   => "Receiver Offset GLONASS",
@@ -291 +324 @@
 
                 $dataset = Chart::Gnuplot::DataSet->new(
-                                                         #xdata   => $EPOCH_OFF_GAL,
-                                                         xdata   => $EPOCH,
+                                                         xdata   => $EPOCH_OFF_GAL,
                                                          ydata   => $file->{'OFF_GAL'},
                                                          title   => "Receiver Offset Galileo",
@@ -325 +357 @@
 
                 $dataset = Chart::Gnuplot::DataSet->new(
-                                                         #xdata   => $EPOCH_OFF_BDS,
-                                                         xdata   => $EPOCH,
+                                                         xdata   => $EPOCH_OFF_BDS,
                                                          ydata   => $file->{'OFF_BDS'},
                                                          title   => "Receiver Offset Beidou",