Changeset 9947 in ntrip for trunk


Ignore:
Timestamp:
Dec 19, 2022, 12:58:57 PM (16 months ago)
Author:
stuerze
Message:

docu updated

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/BNC/scripts/pppPlot.pl

    r9910 r9947  
    66#
    77# Plot metrics:
    8 #    - code and phase residuals in [m],
    9 #    - receiver clock errors in [m],
    10 #    - a-priori and correction values of tropospheric zenith delay in [m],
    11 #    - time offset between GPS time and Galileo/GLONASS/BDS time in [m],
    12 #    - ambiguities, given per satellite
    13 #    - elevations,  given per satellite
     8#    - NEU displacements w.r.t. coordinates in Crd file
     9#    - A-priori + correction values of tropospheric zenith delay in [m],
     10#    - GNSS Receiver clock errors in [m],
     11#    - Elevations, given per satellite
     12#    - Ambiguities, given per satellite
     13#    - Ionosphere Delay [m],
     14#    - Receiver Bias [m]
     15#    - Code and phase residuals in [m],
    1416#
    1517# Author  : Andrea Stuerze
    … …  
    105107    my ( $station, $file ) = Bnc::parsePPPLogfile( $file, $sampling );
    106108    my $EPOCHS       = $file->{'EPOCHS'};
    107     my $EPOCHS_OGR   = $file->{'EPOCHS_OGR'};
    108     my $EPOCHS_OGE   = $file->{'EPOCHS_OGE'};
    109     my $EPOCHS_OGC   = $file->{'EPOCHS_OGC'};
    110109    my $EPOCHS_G_CLK = $file->{'EPOCHS_G_CLK'};
    111110    my $EPOCHS_R_CLK = $file->{'EPOCHS_R_CLK'};
    … …  
    192191        my $chartTRP = newChart($station);
    193192        $chartTRP->set( output => $pngName );
    194         #$chartTRP->set( ylabel => "Tropospheric Delay[m]", yrange => [ " 2.0 ", " 2.6 " ] );
    195         $chartTRP->set( ylabel => "Tropospheric Delay[m]", yrange => [ " 1.0 ", " 3.0 " ] );
     193        $chartTRP->set( ylabel => "Tropospheric Delay [m]", yrange => [ " 1.0 ", " 3.0 " ] );
    196194
    197195        my $dataTRP = Chart::Gnuplot::DataSet->new(
    … …  
    400398
    401399        # AMBIGUITY_TYPE #print Dumper \%AMB;
    402         foreach my $key_ambType (%AMB) {    #print "$key_ambType \n";  ??????
     400        foreach my $key_ambType (%AMB) {    #print "$key_ambType \n";
    403401            foreach my $key_sys ( sort keys %{ $AMB{$key_ambType} } ) {
    404402
    … …  
    493491        }
    494492
    495         ######### RES #####################
    496         DEBUG "Plot Residuals";
    497         $page = $pdf->page();
    498         $page->mediabox('A4');
    499         $y             = $y0 + $dy;
    500         $headline      = sprintf ( "Residuals for station %s", $station );
    501         $headline_text = $page->text;
    502         $headline_text->font( $font1, 11 / pt );
    503         $headline_text->translate( 15 / mm, 280 / mm );
    504         $headline_text->text($headline);
    505 
    506         my $chartRES = newChart($station);
    507         $chartRES->set( legend => { position => "outside right" } );
    508 
    509         # RESIDUAL_TYPE   #print Dumper \%RES;
    510         foreach my $key_resType ( sort keys %RES ) {    #print "key_resType: $key_resType \n";
    511                                                         #SYSTEM
    512             foreach my $key_sys ( sort keys %{ $RES{$key_resType} } ) {
    513 
    514                 #print "key_sys: $key_sys \n"; #print Dumper $RES{$key_resType};
    515                 my @datasets;
    516                 my $pngName = sprintf ( "%s_RES_%s_%s.png", $station, $key_resType, $key_sys );
    517                 $chartRES->set( output => $pngName );
    518                 $chartRES->set( ylabel => "Residuals $key_resType [m]" );
    519 
    520                 if ( $key_resType =~ /^c/ ) {
    521                     #$chartRES->set( yrange => [ " -6.0 ", " 6.0 " ] );
    522                     $chartRES->set( yrange => [ " -10.0 ", " 10.0 " ] );
    523                 }
    524                 elsif ( $key_resType =~ /^l/ ) {
    525                    #$chartRES->set( yrange => [ " -0.06 ", " 0.06 " ] );
    526                     $chartRES->set( yrange => [ " -0.10 ", " 0.10 " ] );
    527                 }
    528 
    529                 elsif ( $key_resType =~ /^GIM/ ) {
    530                     $chartRES->set( yrange => [ " -6.0 ", " 6.0 " ] );
    531                 }
    532 
    533                 # SATELLITE
    534                 foreach my $key_sat ( sort keys %{ $RES{$key_resType}{$key_sys} } ) {
    535 
    536                     #print "$key_sat = $RES{$key_resType}{$key_sys}{$key_sat} \n";
    537                     $dataset = Chart::Gnuplot::DataSet->new(
    538                         xdata   => $RES{$key_resType}{$key_sys}{$key_sat}{EPOCH},  # array of epochs
    539                         ydata   => $RES{$key_resType}{$key_sys}{$key_sat}{DATA},   # array of residuals of one satellite
    540                         title   => "$key_sat",
    541                         timefmt => '%s',
    542 
    543                         #style   => "linespoints",
    544                         style => "linespoints",
    545                     );
    546                     push ( @datasets, $dataset );
    547                 }
    548                 $chartRES->plot2d(@datasets);
    549                 $y = $y - $dy;
    550                 if ( $y < 30 / mm ) {
    551                     $page = $pdf->page();
    552                     $page->mediabox('A4');
    553                     $y = $y0;
    554                 }
    555                 $png = $page->gfx();
    556                 LOGDIE("could not find image file: $!\n") unless -e $pngName;
    557                 $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
    558             }
    559         }
    560 
    561493        ######### ION #####################
    562494        if ( grep ( $_ eq "ALL", @plotTypes ) ) {
    … …  
    656588            }
    657589        }
     590
     591        ######### RES #####################
     592        DEBUG "Plot Residuals";
     593        $page = $pdf->page();
     594        $page->mediabox('A4');
     595        $y             = $y0 + $dy;
     596        $headline      = sprintf ( "Residuals for station %s", $station );
     597        $headline_text = $page->text;
     598        $headline_text->font( $font1, 11 / pt );
     599        $headline_text->translate( 15 / mm, 280 / mm );
     600        $headline_text->text($headline);
     601
     602        my $chartRES = newChart($station);
     603        $chartRES->set( legend => { position => "outside right" } );
     604
     605        # RESIDUAL_TYPE   #print Dumper \%RES;
     606        foreach my $key_resType ( sort keys %RES ) {    #print "key_resType: $key_resType \n";
     607                                                        #SYSTEM
     608            foreach my $key_sys ( sort keys %{ $RES{$key_resType} } ) {
     609
     610                #print "key_sys: $key_sys \n"; #print Dumper $RES{$key_resType};
     611                my @datasets;
     612                my $pngName = sprintf ( "%s_RES_%s_%s.png", $station, $key_resType, $key_sys );
     613                $chartRES->set( output => $pngName );
     614                $chartRES->set( ylabel => "Residuals $key_resType [m]" );
     615
     616                if ( $key_resType =~ /^c/ ) {
     617                    $chartRES->set( yrange => [ " -10.0 ", " 10.0 " ] );
     618                }
     619                elsif ( $key_resType =~ /^l/ ) {
     620                    $chartRES->set( yrange => [ " -0.10 ", " 0.10 " ] );
     621                }
     622
     623                elsif ( $key_resType =~ /^GIM/ ) {
     624                    $chartRES->set( yrange => [ " -6.0 ", " 6.0 " ] );
     625                }
     626
     627                # SATELLITE
     628                foreach my $key_sat ( sort keys %{ $RES{$key_resType}{$key_sys} } ) {
     629
     630                    #print "$key_sat = $RES{$key_resType}{$key_sys}{$key_sat} \n";
     631                    $dataset = Chart::Gnuplot::DataSet->new(
     632                        xdata   => $RES{$key_resType}{$key_sys}{$key_sat}{EPOCH},  # array of epochs
     633                        ydata   => $RES{$key_resType}{$key_sys}{$key_sat}{DATA},   # array of residuals of one satellite
     634                        title   => "$key_sat",
     635                        timefmt => '%s',
     636
     637                        #style   => "linespoints",
     638                        style => "linespoints",
     639                    );
     640                    push ( @datasets, $dataset );
     641                }
     642                $chartRES->plot2d(@datasets);
     643                $y = $y - $dy;
     644                if ( $y < 30 / mm ) {
     645                    $page = $pdf->page();
     646                    $page->mediabox('A4');
     647                    $y = $y0;
     648                }
     649                $png = $page->gfx();
     650                LOGDIE("could not find image file: $!\n") unless -e $pngName;
     651                $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
     652            }
     653        }
    658654    }    # end if ALL @plotTypes
    659655
    … …  
    663659    system ("rm *.png");
    664660    if (Common::amInteractiv ) {
    665 #       system ("evince $inputDir/$pdf_name");
     661#   system ("evince $inputDir/$pdf_name");
    666662    }
    667663}    # -----  next logfile  -----
    … …  
    688684
    689685OPTIONS:
    690   -p,--plotTypes    ALL or NEU (default)
    691   -l,--logFiles     comma separated list of BNC's PPP logfiles
    692   -s,--sampling     sampling interval in seconds for the logfile data (default: 1); for a daily logfile <30> seconds should be usefull
    693   -h,--help         show help contents
     686  -p, --plotTypes    ALL or NEU (default)
     687  -l, --logFiles     comma separated list of BNC's PPP logfiles
     688  -s, --sampling     sampling interval in seconds for the logfile data (default: 1); for a daily logfile <30> seconds should be usefull
     689  -h, --help         show help contents
    694690
    695691EXAMPLES:
Note: See TracChangeset for help on using the changeset viewer.