Changeset 9910 in ntrip for trunk/BNC/scripts/pppPlot.pl


Ignore:
Timestamp:
Nov 24, 2022, 10:29:49 PM (17 months ago)
Author:
stuerze
Message:

ppp plot script updated

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/scripts/pppPlot.pl

    r9836 r9910  
    192192        my $chartTRP = newChart($station);
    193193        $chartTRP->set( output => $pngName );
    194         $chartTRP->set( ylabel => "Tropospheric Delay[m]", yrange => [ " 2.0 ", " 2.6 " ] );
     194        #$chartTRP->set( ylabel => "Tropospheric Delay[m]", yrange => [ " 2.0 ", " 2.6 " ] );
     195        $chartTRP->set( ylabel => "Tropospheric Delay[m]", yrange => [ " 1.0 ", " 3.0 " ] );
    195196
    196197        my $dataTRP = Chart::Gnuplot::DataSet->new(
     
    213214        $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
    214215
    215         ######### CLK #####################
    216         if ( scalar @{ $file->{'CLKs'} } < 1 ) {
    217             DEBUG "No CLKs found";
    218         }
    219         else {
    220             DEBUG "Plot CLKs";
    221             $pngName = sprintf ( "%s_CLK.png", $station );
    222             my $chartCLK = newChart($station);
    223             $chartCLK->set( output => $pngName );
    224             $chartCLK->set( ylabel => "CLK [m]" );
    225 
    226             $dataset = Chart::Gnuplot::DataSet->new(
    227                                                      xdata   => $EPOCHS,
    228                                                      ydata   => $file->{'CLKs'},
    229                                                      title   => "Receiver clock",
    230                                                      timefmt => '%s',
    231                                                      style   => "linespoints",
    232             );
    233             $chartCLK->plot2d($dataset);
    234             $y = $y - $dy;
    235 
    236             if ( $y < 30 / mm ) {
    237                 $page = $pdf->page();
    238                 $page->mediabox('A4');
    239                 $y = $y0;
    240             }
    241             $png = $page->gfx();
    242             LOGDIE("could not find image file: $!\n") unless -e $pngName;
    243             $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
    244         }
    245216        ######### GPS CLK #####################
    246217        if ( scalar @{ $file->{'G_CLKs'} } < 1 ) {
     
    363334            LOGDIE("could not find image file: $!\n") unless -e $pngName;
    364335            $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
    365         }
    366 
    367         ######### OGR #####################
    368         if ( grep ( $_ eq "ALL", @plotTypes ) ) {
    369             if ( scalar @{ $file->{'OGRs'} } < 1 ) {
    370                 DEBUG "No OGRs found";
    371             }
    372             else {
    373                 DEBUG "Plot OGRs";
    374                 my $pngName  = sprintf ( "%s_OGR.png", $station );
    375                 my $chartOGR = newChart($station);
    376                 $chartOGR->set( output => $pngName );
    377                 $chartOGR->set( ylabel => "OGR [m]" );
    378 
    379                 $dataset = Chart::Gnuplot::DataSet->new(
    380                                                          xdata   => $EPOCHS_OGR,
    381                                                          ydata   => $file->{'OGRs'},
    382                                                          title   => "Offset GPS - GLONASS",
    383                                                          timefmt => '%s',
    384                                                          style   => "linespoints",
    385                 );
    386                 $chartOGR->plot2d($dataset);    #system ("display $pngName&");
    387                 $y = $y - $dy;
    388 
    389                 if ( $y < 30 / mm ) {
    390                     $page = $pdf->page();
    391                     $page->mediabox('A4');
    392                     $y = $y0;
    393                 }
    394                 $png = $page->gfx();
    395                 LOGDIE("could not find image file: $!\n") unless -e $pngName;
    396                 $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
    397             }
    398         }
    399 
    400         ######### OGE #####################
    401         if ( grep ( $_ eq "ALL", @plotTypes ) ) {
    402             if ( scalar @{ $file->{'OGEs'} } < 1 ) {
    403                 DEBUG "No OGEs found";
    404             }
    405             else {
    406                 DEBUG "Plot OGEs";
    407                 my $pngName  = sprintf ( "%s_OGE.png", $station );
    408                 my $chartOGE = newChart($station);
    409                 $chartOGE->set( output => $pngName );
    410                 $chartOGE->set( ylabel => "OGE [m]" );
    411 
    412                 $dataset = Chart::Gnuplot::DataSet->new(
    413                                                          xdata   => $EPOCHS_OGE,
    414                                                          ydata   => $file->{'OGEs'},
    415                                                          title   => "Offset GPS - Galileo",
    416                                                          timefmt => '%s',
    417                                                          style   => "linespoints",
    418                 );
    419                 $chartOGE->plot2d($dataset);    #system ("display $pngName&");
    420                 $y = $y - $dy;
    421 
    422                 if ( $y < 30 / mm ) {
    423                     $page = $pdf->page();
    424                     $page->mediabox('A4');
    425                     $y = $y0;
    426                 }
    427                 $png = $page->gfx();
    428                 LOGDIE("could not find image file: $!\n") unless -e $pngName;
    429                 $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
    430             }
    431         }
    432 
    433         ######### OGC #####################
    434         if ( grep ( $_ eq "ALL", @plotTypes ) ) {
    435             if ( scalar @{ $file->{'OGCs'} } < 1 ) {
    436                 DEBUG "No OGCs found";
    437             }
    438             else {
    439                 DEBUG "Plot OGCs";
    440                 my $pngName  = sprintf ( "%s_OGC.png", $station );
    441                 my $chartOGC = newChart($station);
    442                 $chartOGC->set( output => $pngName );
    443                 $chartOGC->set( ylabel => "OGC [m]" );
    444 
    445                 $dataset = Chart::Gnuplot::DataSet->new(
    446                                                          xdata   => $EPOCHS_OGC,
    447                                                          ydata   => $file->{'OGCs'},
    448                                                          title   => "Offset GPS - BDS",
    449                                                          timefmt => '%s',
    450                                                          style   => "linespoints",
    451                 );
    452                 $chartOGC->plot2d($dataset);    #system ("display $pngName&");
    453                 $y = $y - $dy;
    454 
    455                 if ( $y < 30 / mm ) {
    456                     $page = $pdf->page();
    457                     $page->mediabox('A4');
    458                     $y = $y0;
    459                 }
    460                 $png = $page->gfx();
    461                 LOGDIE("could not find image file: $!\n") unless -e $pngName;
    462                 $png->image( $pdf->image_png($pngName), $x, $y, $width, $height );
    463             }
    464336        }
    465337
     
    647519
    648520                if ( $key_resType =~ /^c/ ) {
    649                     $chartRES->set( yrange => [ " -6.0 ", " 6.0 " ] );
     521                    #$chartRES->set( yrange => [ " -6.0 ", " 6.0 " ] );
     522                    $chartRES->set( yrange => [ " -10.0 ", " 10.0 " ] );
    650523                }
    651524                elsif ( $key_resType =~ /^l/ ) {
    652                     $chartRES->set( yrange => [ " -0.06 ", " 0.06 " ] );
     525                   #$chartRES->set( yrange => [ " -0.06 ", " 0.06 " ] );
     526                    $chartRES->set( yrange => [ " -0.10 ", " 0.10 " ] );
    653527                }
    654528
     
    787661    $pdf->end();
    788662
    789     #system ("rm *.png");
     663    system ("rm *.png");
    790664    if (Common::amInteractiv ) {
    791         system ("evince $inputDir/$pdf_name");
     665#       system ("evince $inputDir/$pdf_name");
    792666    }
    793667}    # -----  next logfile  -----
Note: See TracChangeset for help on using the changeset viewer.