Changeset 10383 in ntrip


Ignore:
Timestamp:
Mar 11, 2024, 12:52:02 PM (5 weeks ago)
Author:
stuerze
Message:

minor changes

Location:
trunk/BNC/scripts
Files:
2 edited

Legend:

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

    r10259 r10383  
    493493    my $epochSec     = 0;
    494494    my $epochDiff    = 0;
    495     my ( @hlp, @N, @E, @U, %SATNUM, @TRPs, @CLKs,   @OFF_GLOs, @OFF_GALs, @OFF_BDSs);
    496     my ( @EPOCHs, @EPOCHs_CLK, @EPOCHs_OFF_GLO, @EPOCHs_OFF_GAL, @EPOCHs_OFF_BDS );
     495    my ( @hlp, @N, @E, @U, %SATNUM, @TRP, @CLK, @OFF_GLO, @OFF_GAL, @OFF_BDS);
     496    my ( @EPOCH, @EPOCH_CLK, @EPOCH_OFF_GLO, @EPOCH_OFF_GAL, @EPOCH_OFF_BDS );
    497497    my ( %AMB,        %RES,            %ELE,            %ION, %BIA );
    498498    my ( $station,    $lki,            $sys,            $sat, $amb );
     
    530530
    531531        if ( $ln =~ /\bdN\b/ ) {
    532             push ( @EPOCHs, $epochSec );    # besser $epo ?
     532            push ( @EPOCH, $epochSec );    # besser $epo ?
    533533            $old_epochSec = $epochSec;
    534534
    535 #2015-08-27_13:59:50.000 DIEP1 X = 3842152.9054 +- 0.0242 Y = 563402.0331 +- 0.0176 Z = 5042888.5182 +- 0.0319 dN = 0.0130 +- 0.0193 dE = -0.0032 +- 0.0178 dU = -0.0248 +- 0.0349
     535            #2015-08-27_13:59:50.000 DIEP1 X = 3842152.9054 +- 0.0242 Y = 563402.0331 +- 0.0176 Z = 5042888.5182 +- 0.0319 dN = 0.0130 +- 0.0193 dE = -0.0032 +- 0.0178 dU = -0.0248 +- 0.0349
    536536            $station = $hlp[1];
    537537
     
    585585        }
    586586        elsif ( $ln =~ /\bREC_CLK\b/ ) {
    587             push ( @EPOCHs_CLK, $epochSec );
    588             push ( @CLKs,       $hlp[2] + $hlp[3] );
     587            push ( @EPOCH_CLK, $epochSec );
     588            push ( @CLK,       $hlp[2] + $hlp[3] );
    589589        }
    590590        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] );
     591            push ( @EPOCH_OFF_GLO, $epochSec );
     592            push ( @OFF_GLO,       $hlp[2] + $hlp[3] );
    593593        }
    594594        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] );
     595            push ( @EPOCH_OFF_GAL, $epochSec );
     596            push ( @OFF_GAL,       $hlp[2] + $hlp[3] );
    597597        }
    598598        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] );
     599            push ( @EPOCH_OFF_BDS, $epochSec );
     600            push ( @OFF_BDS,       $hlp[2] + $hlp[3] );
    601601        }
    602602        elsif ( $ln =~ /\bSATNUM\b/ ) {     # 2015-09... SATNUM G  8
     
    604604        }
    605605        elsif ( $ln =~ /\bTRP\b/ ) {        # 2015-08... TRP  2.3803 +0.1009 +- 0.0324
    606             push ( @TRPs, $hlp[2] + $hlp[3] );
     606            push ( @TRP, $hlp[2] + $hlp[3] );
    607607        }
    608608    }    # -----  next line  -----
     
    611611    close $fh;
    612612
    613     my $nof_epochs = scalar @EPOCHs;
     613    my $nof_epochs = scalar @EPOCH;
    614614    DEBUG(   "$station: epochs:$nof_epochs, North displac.: "
    615615           . scalar @N
     
    618618           . ", Up displac.: "
    619619           . scalar @U
    620            . ", TRPs:"
    621            . scalar @TRPs
     620           . ", TRP:"
     621           . scalar @TRP
    622622           );
    623     if ( $nof_epochs != scalar @N )                                { LOGDIE "number of epochs and residuals not equal\n" }
    624     if ( $nof_epochs != scalar @TRPs )                             { LOGDIE "number of epochs and TRPs not equal\n" }
    625     if ( @CLKs     && scalar @EPOCHs_CLK != scalar @CLKs )         { LOGDIE "number of epochs and CLKs not equal\n" }
    626     if ( @OFF_GLOs && scalar @EPOCHs_OFF_GLO != scalar @OFF_GLOs ) { LOGDIE "number of epochs and OFF_GLOs not equal\n" }
    627     if ( @OFF_GALs && scalar @EPOCHs_OFF_GAL != scalar @OFF_GALs ) { LOGDIE "number of epochs and OFF_GALs not equal\n" }
    628     if ( @OFF_BDSs && scalar @EPOCHs_OFF_BDS != scalar @OFF_BDSs ) { LOGDIE "number of epochs and OFF_BDSs not equal\n" }
     623    if ( $nof_epochs != scalar @N )                              { LOGDIE "number of epochs and residuals not equal\n" }
     624    if ( $nof_epochs != scalar @TRP )                            { LOGDIE "number of epochs and TRP not equal\n" }
     625    if ( @CLK     && scalar @EPOCH_CLK     != scalar @CLK )      { LOGDIE "number of epochs and CLK not equal\n" }
     626    if ( @OFF_GLO && scalar @EPOCH_OFF_GLO != scalar @OFF_GLO )  { LOGDIE "number of epochs and OFF_GLO not equal\n" }
     627    if ( @OFF_GAL && scalar @EPOCH_OFF_GAL != scalar @OFF_GAL )  { LOGDIE "number of epochs and OFF_GAL not equal\n" }
     628    if ( @OFF_BDS && scalar @EPOCH_OFF_BDS != scalar @OFF_BDS )  { LOGDIE "number of epochs and OFF_BDS not equal\n" }
    629629
    630630    if ( !$station ) { WARN "could not grep stationname from file: $file\n" }
    631631
    632632    my %data = (
    633                  EPOCHS     => \@EPOCHs,
    634                  N          => \@N,
    635                  E          => \@E,
    636                  U          => \@U,
    637                  SATNUM     => \%SATNUM,
    638                  TRPs       => \@TRPs,
    639                  CLKs      => \@CLKs,
    640                  OFF_GLOs  => \@OFF_GLOs,
    641                  OFF_GALs  => \@OFF_GALs,
    642                  OFF_BDSs  => \@OFF_BDSs,
    643                  RES        => \%RES,
    644                  AMB        => \%AMB,
    645                  ELE        => \%ELE,
    646                  ION        => \%ION,
    647                  BIA        => \%BIA,
     633                 EPOCH   => \@EPOCH,
     634                 N       => \@N,
     635                 E       => \@E,
     636                 U       => \@U,
     637                 SATNUM  => \%SATNUM,
     638                 TRP     => \@TRP,
     639                 CLK     => \@CLK,
     640                 OFF_GLO => \@OFF_GLO,
     641                 OFF_GAL => \@OFF_GAL,
     642                 OFF_BDS => \@OFF_BDS,
     643                 RES     => \%RES,
     644                 AMB     => \%AMB,
     645                 ELE     => \%ELE,
     646                 ION     => \%ION,
     647                 BIA     => \%BIA,
    648648    );
    649649
  • trunk/BNC/scripts/pppPlot.pl

    r10259 r10383  
    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 $EPOCH       = $file->{'EPOCH'};
     109    my $EPOCH_OFF_GLO = $file->{'EPOCH_OFF_GLO'};
     110    my $EPOCH_OFF_GAL = $file->{'EPOCH_OFF_GAL'};
     111    my $EPOCH_OFF_BDS = $file->{'EPOCH_OFF_BDS'};
     112    my $EPOCH_CLK     = $file->{'EPOCH_CLK'};
    113113    my %AMB            = %{ $file->{'AMB'} };
    114114    my %RES            = %{ $file->{'RES'} };
     
    128128    $u = pdl( $file->{'U'} );
    129129    ( $mean, $prms, $median, $min, $max, $adev, $rms_u ) = stats($u);
    130     $trp = pdl( $file->{'TRPs'} );
     130    $trp = pdl( $file->{'TRP'} );
    131131    ( $mean, $prms, $median, $min, $max, $adev, $rms_trp ) = stats($trp);
    132132    $str_rms_n   = sprintf ( " %.2f ", $rms_n );
     
    158158
    159159    my $dataN = Chart::Gnuplot::DataSet->new(
    160                                               xdata   => $EPOCHS,
     160                                              xdata   => $EPOCH,
    161161                                              ydata   => $file->{'N'},
    162162                                              title   => "Displacements N, RMS + -$str_rms_n m",
     
    165165    );
    166166    my $dataE = Chart::Gnuplot::DataSet->new(
    167                                               xdata   => $EPOCHS,
     167                                              xdata   => $EPOCH,
    168168                                              ydata   => $file->{'E'},
    169169                                              title   => "Displacements E, RMS + -$str_rms_e m",
     
    172172    );
    173173    my $dataU = Chart::Gnuplot::DataSet->new(
    174                                               xdata   => $EPOCHS,
     174                                              xdata   => $EPOCH,
    175175                                              ydata   => $file->{'U'},
    176176                                              title   => "Displacements U, RMS + -$str_rms_u m",
     
    187187    ######### TRP #####################
    188188    if ( grep ( $_ eq "ALL", @plotTypes ) ) {
    189         DEBUG "Plot TRPs";
     189        DEBUG "Plot TRP";
    190190        my $pngName  = sprintf ( "%s_TRP.png", $station );
    191191        my $chartTRP = newChart($station);
     
    194194
    195195        my $dataTRP = Chart::Gnuplot::DataSet->new(
    196                                                     xdata   => $EPOCHS,
    197                                                     ydata   => $file->{'TRPs'},
     196                                                    xdata   => $EPOCH,
     197                                                    ydata   => $file->{'TRP'},
    198198                                                    title   => "Tropospheric Delay, RMS + -$str_rms_trp m",
    199199                                                    timefmt => '%s',
     
    213213
    214214        ######### CLK #####################
    215         if ( scalar @{ $file->{'CLKs'} } < 1 ) {
    216             DEBUG "No CLKs found";
     215        if ( scalar @{ $file->{'CLK'} } < 1 ) {
     216            DEBUG "No CLK found";
    217217        }
    218218        else {
    219             DEBUG "Plot CLKs";
     219            DEBUG "Plot CLK";
    220220            $pngName = sprintf ( "%s_CLK.png", $station );
    221221            my $chartCLK = newChart($station);
     
    224224
    225225            $dataset = Chart::Gnuplot::DataSet->new(
    226                                                      xdata   => $EPOCHS,
    227                                                      ydata   => $file->{'CLKs'},
     226                                                     xdata   => $EPOCH,
     227                                                     ydata   => $file->{'CLK'},
    228228                                                     title   => "Receiver clock",
    229229                                                     timefmt => '%s',
     
    246246        ######### OFF_GLO #####################
    247247        if ( grep ( $_ eq "ALL", @plotTypes ) ) {
    248             if ( scalar @{ $file->{'OFF_GLOs'} } < 1 ) {
    249                 DEBUG "No OFF_GLOs found";
     248            if ( scalar @{ $file->{'OFF_GLO'} } < 1 ) {
     249                DEBUG "No OFF_GLO found";
    250250            }
    251251            else {
    252                 DEBUG "Plot OFF_GLOs";
     252                DEBUG "Plot OFF_GLO";
    253253                my $pngName  = sprintf ( "%s_OFF_R.png", $station );
    254254                my $chartOFF_GLO = newChart($station);
     
    257257
    258258                $dataset = Chart::Gnuplot::DataSet->new(
    259                                                          xdata   => $EPOCHS_OFF_GLO,
    260                                                          ydata   => $file->{'OFF_GLOs'},
     259                                                         xdata   => $EPOCH_OFF_GLO,
     260                                                         ydata   => $file->{'OFF_GLO'},
    261261                                                         title   => "Receiver Offset GLONASS",
    262262                                                         timefmt => '%s',
     
    279279        ######### OFF_GAL #####################
    280280        if ( grep ( $_ eq "ALL", @plotTypes ) ) {
    281             if ( scalar @{ $file->{'OFF_GALs'} } < 1 ) {
    282                 DEBUG "No OFF_GALs found";
     281            if ( scalar @{ $file->{'OFF_GAL'} } < 1 ) {
     282                DEBUG "No OFF_GAL found";
    283283            }
    284284            else {
    285                 DEBUG "Plot OFF_GALs";
     285                DEBUG "Plot OFF_GAL";
    286286                my $pngName  = sprintf ( "%s_OFF_E.png", $station );
    287287                my $chartOFF_GAL = newChart($station);
     
    290290
    291291                $dataset = Chart::Gnuplot::DataSet->new(
    292                                                          xdata   => $EPOCHS_OFF_GAL,
    293                                                          ydata   => $file->{'OFF_GALs'},
     292                                                         xdata   => $EPOCH_OFF_GAL,
     293                                                         ydata   => $file->{'OFF_GAL'},
    294294                                                         title   => "Receiver Offset Galileo",
    295295                                                         timefmt => '%s',
     
    312312        ######### OFF_BDS #####################
    313313        if ( grep ( $_ eq "ALL", @plotTypes ) ) {
    314             if ( scalar @{ $file->{'OFF_BDSs'} } < 1 ) {
    315                 DEBUG "No OFF_BDSs found";
     314            if ( scalar @{ $file->{'OFF_BDS'} } < 1 ) {
     315                DEBUG "No OFF_BDS found";
    316316            }
    317317            else {
    318                 DEBUG "Plot OFF_BDSs";
     318                DEBUG "Plot OFF_BDS";
    319319                my $pngName  = sprintf ( "%s_OFF_C.png", $station );
    320320                my $chartOFF_BDS = newChart($station);
     
    323323
    324324                $dataset = Chart::Gnuplot::DataSet->new(
    325                                                          xdata   => $EPOCHS_OFF_BDS,
    326                                                          ydata   => $file->{'OFF_BDSs'},
     325                                                         xdata   => $EPOCH_OFF_BDS,
     326                                                         ydata   => $file->{'OFF_BDS'},
    327327                                                         title   => "Receiver Offset Beidou",
    328328                                                         timefmt => '%s',
Note: See TracChangeset for help on using the changeset viewer.