Changeset 10233 in ntrip


Ignore:
Timestamp:
Nov 1, 2023, 9:02:50 AM (7 months ago)
Author:
stuerze
Message:

minor changes

Location:
trunk/BNC/scripts
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/scripts/Bnc.pm

    r10090 r10233  
    493493    my $epochSec     = 0;
    494494    my $epochDiff    = 0;
    495     my ( @hlp, @EPOCHs, @N, @E, @U, %SATNUM, @TRPs, @CLKs, @OFF_GLOs, @OFF_GALs, @OFF_BDSs, );
    496     my ( @EPOCHs_CLK, @EPOCHs_OFF_GLO, @EPOCHs_OFF_GAL, @EPOCHs_OFF_BDS );
     495    my ( @hlp, @N, @E, @U, %SATNUM, @TRPs, @G_CLKs, @R_CLKs, @E_CLKs, @C_CLKs );
     496    my ( @EPOCHs, @EPOCHs_G_CLK, @EPOCHs_R_CLK, @EPOCHs_E_CLK, @EPOCHs_C_CLK);
    497497    my ( %AMB,        %RES,            %ELE,            %ION, %BIA );
    498498    my ( $station,    $lki,            $sys,            $sat, $amb );
     
    584584            push @{ $BIA{$lki}{$sys}{DATA} },  $hlp[4] + $hlp[5];
    585585        }
    586         elsif ( $ln =~ /\bREC_CLK\b/ ) {
    587             push ( @EPOCHs_CLK, $epochSec );
    588             push ( @CLKs,       $hlp[2] + $hlp[3] );
    589         }
    590         elsif ( $ln =~ /\bOFF_GLO\b/ ) {    # 2015-08... OFF_GLO 52.6806 -3.8042 +- 9.0077
    591             push ( @EPOCHs_OFF_GLO, $epochSec );
    592             push ( @OFF_GLOs,       $hlp[2] + $hlp[3] );
    593         }
    594         elsif ( $ln =~ /\bOFF_GAL\b/ ) {    # 2015-08... OFF_GAL 52.6806 -3.8042 +- 9.0077
    595             push ( @EPOCHs_OFF_GAL, $epochSec );
    596             push ( @OFF_GALs,       $hlp[2] + $hlp[3] );
    597         }
    598         elsif ( $ln =~ /\bOFF_BDS\b/ ) {    # 2015-08... OFF_BDS 52.6806 -3.8042 +- 9.0077
    599             push ( @EPOCHs_OFF_BDS, $epochSec );
    600             push ( @OFF_BDSs,       $hlp[2] + $hlp[3] );
     586        elsif ( $ln =~ /\bREC_CLK  G\b/ ){
     587            push ( @EPOCHs_G_CLK, $epochSec);
     588                push ( @G_CLKs, $hlp[3] + $hlp[4] );
     589        }
     590        elsif ( $ln =~ /\bREC_CLK  R\b/ ){
     591            push ( @EPOCHs_R_CLK, $epochSec);
     592                push ( @R_CLKs, $hlp[3] + $hlp[4] );
     593        }
     594        elsif ( $ln =~ /\bREC_CLK  E\b/ ){
     595            push ( @EPOCHs_E_CLK, $epochSec);
     596                push ( @E_CLKs, $hlp[3] + $hlp[4] );
     597        }
     598        elsif ( $ln =~ /\bREC_CLK  C\b/ ){
     599            push ( @EPOCHs_C_CLK, $epochSec);
     600                push ( @C_CLKs, $hlp[3] + $hlp[4] );
    601601        }
    602602        elsif ( $ln =~ /\bSATNUM\b/ ) {     # 2015-09... SATNUM G  8
     
    619619           . scalar @U
    620620           . ", TRPs:"
    621            . scalar @TRPs );
     621           . scalar @TRPs
     622           . ", GPS CLKs:"
     623           . scalar @G_CLKs
     624           . ", GLO CLKs:"
     625           . scalar @R_CLKs
     626           . ", GAL CLKs:"
     627           . scalar @E_CLKs
     628           . ", BDS CLKs:"
     629           . scalar @E_CLKs
     630           );
    622631    if ( $nof_epochs != scalar @N )                    { LOGDIE "number of epochs and residuals not equal\n" }
    623632    if ( $nof_epochs != scalar @TRPs )                 { LOGDIE "number of epochs and TRPs not equal\n" }
    624     if ( @CLKs && scalar @EPOCHs_CLK != scalar @CLKs ) { LOGDIE "number of epochs and CLKs not equal\n" }
    625 
    626     if ( @OFF_GLOs && scalar @EPOCHs_OFF_GLO != scalar @OFF_GLOs ) {
    627         LOGDIE "number of epochs and OFF_GLOs not equal\n";
    628     }
    629     if ( @OFF_GALs && scalar @EPOCHs_OFF_GAL != scalar @OFF_GALs ) {
    630         LOGDIE "number of epochs and OFF_GALs not equal\n";
    631     }
    632     if ( @OFF_BDSs && scalar @EPOCHs_OFF_BDS != scalar @OFF_BDSs ) {
    633         LOGDIE "number of epochs and OFF_BDSs not equal\n";
    634     }
     633    if ( @G_CLKs && scalar @EPOCHs_G_CLK != scalar @G_CLKs ) {LOGDIE "number of epochs and G_CLKs not equal\n";}
     634    if ( @R_CLKs && scalar @EPOCHs_R_CLK != scalar @R_CLKs ) {LOGDIE "number of epochs and R_CLKs not equal\n";}
     635    if ( @E_CLKs && scalar @EPOCHs_E_CLK != scalar @E_CLKs ) {LOGDIE "number of epochs and E_CLKs not equal\n";}
     636    if ( @C_CLKs && scalar @EPOCHs_C_CLK != scalar @C_CLKs ) {LOGDIE "number of epochs and C_CLKs not equal\n";}
    635637    if ( !$station ) { WARN "could not grep stationname from file: $file\n" }
    636638
     
    642644                 SATNUM     => \%SATNUM,
    643645                 TRPs       => \@TRPs,
    644                  CLKs       => \@CLKs,
    645                  CLK_EPOCHS => \@EPOCHs_CLK,
    646                  OFF_GLOs   => \@OFF_GLOs,
    647                  OFF_GALs   => \@OFF_GALs,
    648                  OFF_BDSs   => \@OFF_BDSs,
     646                 G_CLKs  => \@G_CLKs,
     647                 R_CLKs  => \@R_CLKs,
     648                 E_CLKs  => \@E_CLKs,
     649                 C_CLKs  => \@C_CLKs,
    649650                 RES        => \%RES,
    650651                 AMB        => \%AMB,
  • trunk/BNC/scripts/pppPlot.pl

    r10197 r10233  
    106106    # -----------------------------------------------------------------------------
    107107    my ( $station, $file ) = Bnc::parsePPPLogfile( $file, $sampling );
    108     my $EPOCHS         = $file->{'EPOCHS'};
    109     my $EPOCHS_OFF_GLO = $file->{'EPOCHS_OFF_GLO'};
    110     my $EPOCHS_OFF_GAL = $file->{'EPOCHS_OFF_GAL'};
    111     my $EPOCHS_OFF_BDS = $file->{'EPOCHS_OFF_BDS'};
    112     my $EPOCHS_CLK     = $file->{'EPOCHS_CLK'};
     108    my $EPOCHS       = $file->{'EPOCHS'};
     109    my $EPOCHS_G_CLK = $file->{'EPOCHS_G_CLK'};
     110    my $EPOCHS_R_CLK = $file->{'EPOCHS_R_CLK'};
     111    my $EPOCHS_E_CLK = $file->{'EPOCHS_E_CLK'};
     112    my $EPOCHS_C_CLK = $file->{'EPOCHS_C_CLK'};
    113113    my %AMB            = %{ $file->{'AMB'} };
    114114    my %RES            = %{ $file->{'RES'} };
     
    212212        $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
    213213
    214         ######### CLK #####################
    215         if ( scalar @{ $file->{'CLKs'} } < 1 ) {
    216             DEBUG "No CLKs found";
     214        ######### GPS CLK #####################
     215        if ( scalar @{ $file->{'G_CLKs'} } < 1 ) {
     216          DEBUG "No GPS CLKs found";
    217217        }
    218218        else {
    219             DEBUG "Plot CLKs";
    220             $pngName = sprintf ( "%s_CLK.png", $station );
    221             my $chartCLK = newChart($station);
    222             $chartCLK->set( output => $pngName );
    223             $chartCLK->set( ylabel => "CLK [m]" );
    224 
    225             $dataset = Chart::Gnuplot::DataSet->new(
    226                                                      xdata   => $EPOCHS,
    227                                                      ydata   => $file->{'CLKs'},
    228                                                      title   => "Receiver clock",
    229                                                      timefmt => '%s',
    230                                                      style   => "linespoints",
    231             );
    232             $chartCLK->plot2d($dataset);
    233             $y = $y - $dy;
    234 
    235             if ( $y < 30 / mm ) {
    236                 $page = $pdf->page();
    237                 $page->mediabox('A4');
    238                 $y = $y0;
    239             }
    240             $png = $page->gfx();
    241             LOGDIE("could not find image file: $!\n") unless -e $pngName;
    242             $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
    243         }
    244 
    245 
    246         ######### OFF_GLO #####################
    247         if ( grep ( $_ eq "ALL", @plotTypes ) ) {
    248             if ( scalar @{ $file->{'OFF_GLOs'} } < 1 ) {
    249                 DEBUG "No OFF_GLOs found";
    250             }
    251             else {
    252                 DEBUG "Plot OFF_GLOs";
    253                 my $pngName  = sprintf ( "%s_OFF_R.png", $station );
    254                 my $chartOFF_GLO = newChart($station);
    255                 $chartOFF_GLO->set( output => $pngName );
    256                 $chartOFF_GLO->set( ylabel => "Offset GLONASS [m]" );
    257 
    258                 $dataset = Chart::Gnuplot::DataSet->new(
    259                                                          xdata   => $EPOCHS_OFF_GLO,
    260                                                          ydata   => $file->{'OFF_GLOs'},
    261                                                          title   => "Receiver Offset GLONASS",
    262                                                          timefmt => '%s',
    263                                                          style   => "linespoints",
    264                 );
    265                 $chartOFF_GLO->plot2d($dataset);    #system ("display $pngName&");
    266                 $y = $y - $dy;
    267 
    268                 if ( $y < 30 / mm ) {
    269                     $page = $pdf->page();
    270                     $page->mediabox('A4');
    271                     $y = $y0;
    272                 }
    273                 $png = $page->gfx();
    274                 LOGDIE("could not find image file: $!\n") unless -e $pngName;
    275                 $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
    276             }
    277         }
    278 
    279         ######### OFF_GAL #####################
    280         if ( grep ( $_ eq "ALL", @plotTypes ) ) {
    281             if ( scalar @{ $file->{'OFF_GALs'} } < 1 ) {
    282                 DEBUG "No OFF_GALs found";
    283             }
    284             else {
    285                 DEBUG "Plot OFF_GALs";
    286                 my $pngName  = sprintf ( "%s_OFF_E.png", $station );
    287                 my $chartOFF_GAL = newChart($station);
    288                 $chartOFF_GAL->set( output => $pngName );
    289                 $chartOFF_GAL->set( ylabel => "Offset Galileo [m]" );
    290 
    291                 $dataset = Chart::Gnuplot::DataSet->new(
    292                                                          xdata   => $EPOCHS_OFF_GAL,
    293                                                          ydata   => $file->{'OFF_GALs'},
    294                                                          title   => "Receiver Offset Galileo",
    295                                                          timefmt => '%s',
    296                                                          style   => "linespoints",
    297                 );
    298                 $chartOFF_GAL->plot2d($dataset);    #system ("display $pngName&");
    299                 $y = $y - $dy;
    300 
    301                 if ( $y < 30 / mm ) {
    302                     $page = $pdf->page();
    303                     $page->mediabox('A4');
    304                     $y = $y0;
    305                 }
    306                 $png = $page->gfx();
    307                 LOGDIE("could not find image file: $!\n") unless -e $pngName;
    308                 $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
    309             }
    310         }
    311 
    312         ######### OFF_BDS #####################
    313         if ( grep ( $_ eq "ALL", @plotTypes ) ) {
    314             if ( scalar @{ $file->{'OFF_BDSs'} } < 1 ) {
    315                 DEBUG "No OFF_BDSs found";
    316             }
    317             else {
    318                 DEBUG "Plot OFF_BDSs";
    319                 my $pngName  = sprintf ( "%s_OFF_C.png", $station );
    320                 my $chartOFF_BDS = newChart($station);
    321                 $chartOFF_BDS->set( output => $pngName );
    322                 $chartOFF_BDS->set( ylabel => "Offset Beidou [m]" );
    323 
    324                 $dataset = Chart::Gnuplot::DataSet->new(
    325                                                          xdata   => $EPOCHS_OFF_BDS,
    326                                                          ydata   => $file->{'OFF_BDSs'},
    327                                                          title   => "Receiver Offset Beidou",
    328                                                          timefmt => '%s',
    329                                                          style   => "linespoints",
    330                 );
    331                 $chartOFF_BDS->plot2d($dataset);    #system ("display $pngName&");
    332                 $y = $y - $dy;
    333 
    334                 if ( $y < 30 / mm ) {
    335                     $page = $pdf->page();
    336                     $page->mediabox('A4');
    337                     $y = $y0;
    338                 }
    339                 $png = $page->gfx();
    340                 LOGDIE("could not find image file: $!\n") unless -e $pngName;
    341                 $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
    342             }
     219          DEBUG "Plot GPS CLKs";
     220          $pngName = sprintf ( "%s_CLK_G.png", $station );
     221          my $chartCLK = newChart($station);
     222          $chartCLK->set( output => $pngName );
     223          $chartCLK->set( ylabel => "CLK [m]" );
     224
     225          $dataset = Chart::Gnuplot::DataSet->new(
     226                                                   xdata   => $EPOCHS_G_CLK,
     227                                                   ydata   => $file->{'G_CLKs'},
     228                                                   title   => "GPS Receiver clock",
     229                                                   timefmt => '%s',
     230                                                   style   => "points",
     231          );
     232          $chartCLK->plot2d($dataset);
     233          $y = $y - $dy;
     234
     235          if ( $y < 30 / mm ) {
     236              $page = $pdf->page();
     237              $page->mediabox('A4');
     238              $y = $y0;
     239          }
     240          $png = $page->gfx();
     241          LOGDIE("Unable to find image file: $!\n") unless -e $pngName;
     242          $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
     243        }
     244         ######### GLONASS CLK #####################
     245        if ( scalar @{ $file->{'R_CLKs'} } < 1 ) {
     246          DEBUG "No GLONASS CLKs found";
     247        }
     248        else {
     249          DEBUG "Plot GLONASS CLKs";
     250          $pngName = sprintf ( "%s_CLK_R.png", $station );
     251          my $chartCLK = newChart($station);
     252          $chartCLK->set( output => $pngName );
     253          $chartCLK->set( ylabel => "CLK [m]" );
     254
     255          $dataset = Chart::Gnuplot::DataSet->new(
     256                                                   xdata   => $EPOCHS_R_CLK,
     257                                                   ydata   => $file->{'R_CLKs'},
     258                                                   title   => "GLONASS Receiver clock",
     259                                                   timefmt => '%s',
     260                                                   style   => "points",
     261          );
     262          $chartCLK->plot2d($dataset);
     263          $y = $y - $dy;
     264
     265          if ( $y < 30 / mm ) {
     266              $page = $pdf->page();
     267              $page->mediabox('A4');
     268              $y = $y0;
     269          }
     270          $png = $page->gfx();
     271          LOGDIE("Unable to find image file: $!\n") unless -e $pngName;
     272          $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
     273        }
     274         ######### Galileo CLK #####################
     275        if ( scalar @{ $file->{'E_CLKs'} } < 1 ) {
     276          DEBUG "No Galileo CLKs found";
     277        }
     278        else {
     279          DEBUG "Plot Galileo CLKs";
     280          $pngName = sprintf ( "%s_CLK_E.png", $station );
     281          my $chartCLK = newChart($station);
     282          $chartCLK->set( output => $pngName );
     283          $chartCLK->set( ylabel => "CLK [m]" );
     284
     285          $dataset = Chart::Gnuplot::DataSet->new(
     286                                                   xdata   => $EPOCHS_E_CLK,
     287                                                   ydata   => $file->{'E_CLKs'},
     288                                                   title   => "Galileo Receiver clock",
     289                                                   timefmt => '%s',
     290                                                   style   => "points",
     291          );
     292          $chartCLK->plot2d($dataset);
     293          $y = $y - $dy;
     294
     295          if ( $y < 30 / mm ) {
     296              $page = $pdf->page();
     297              $page->mediabox('A4');
     298              $y = $y0;
     299          }
     300          $png = $page->gfx();
     301          LOGDIE("Unable to find image file: $!\n") unless -e $pngName;
     302          $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
     303        }
     304
     305        ######### Beidou CLK #####################
     306        if ( scalar @{ $file->{'C_CLKs'} } < 1 ) {
     307          DEBUG "No BDS CLKs found";
     308        }
     309        else {
     310          DEBUG "Plot BDS CLKs";
     311          $pngName = sprintf ( "%s_CLK_C.png", $station );
     312          my $chartCLK = newChart($station);
     313          $chartCLK->set( output => $pngName );
     314          $chartCLK->set( ylabel => "CLK [m]" );
     315
     316          $dataset = Chart::Gnuplot::DataSet->new(
     317                                                   xdata   => $EPOCHS_C_CLK,
     318                                                   ydata   => $file->{'C_CLKs'},
     319                                                   title   => "BDS Receiver clock",
     320                                                   timefmt => '%s',
     321                                                   style   => "points",
     322          );
     323          $chartCLK->plot2d($dataset);
     324          $y = $y - $dy;
     325
     326          if ( $y < 30 / mm ) {
     327              $page = $pdf->page();
     328              $page->mediabox('A4');
     329              $y = $y0;
     330          }
     331          $png = $page->gfx();
     332          LOGDIE("Unable to find image file: $!\n") unless -e $pngName;
     333          $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
    343334        }
    344335
Note: See TracChangeset for help on using the changeset viewer.