Changeset 10005 in ntrip for trunk/BNC/scripts
- Timestamp:
- Mar 16, 2023, 12:30:50 PM (21 months ago)
- Location:
- trunk/BNC/scripts
- Files:
-
- 2 edited
-
Bnc.pm (modified) (6 diffs)
-
pppPlot.pl (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/scripts/Bnc.pm
r9910 r10005 304 304 my @mpts_def = (); 305 305 foreach (@mpts) { 306 307 306 # user:passwd@igs-ip.net:2101/ASPA0 RTCM_3.0 ASM -14.33 189.28 no 1 308 307 if ( $_ =~ … … 493 492 my $epochDiff = 0; 494 493 my ( 495 @hlp, @EPOCHs, @N, @E, @U, %SATNUM, @TRPs, 496 @G_CLKs, @R_CLKs, @E_CLKs, @C_CLKs, 494 @hlp, @EPOCHs, @N, @E, @U, 495 %SATNUM, @TRPs, 496 @CLKs, @OFF_GLOs, @OFF_GALs, @OFF_BDSs, 497 497 ); 498 my ( @EPOCHs_ G_CLK, @EPOCHs_R_CLK, @EPOCHs_E_CLK, @EPOCHs_C_CLK);498 my ( @EPOCHs_CLK, @EPOCHs_OFF_GLO, @EPOCHs_OFF_GAL, @EPOCHs_OFF_BDS ); 499 499 my ( %AMB, %RES, %ELE, %ION, %BIA ); 500 500 my ( $station, $lki, $sys, $sat, $amb ); … … 584 584 push @{ $BIA{$lki}{$sys}{DATA} }, $hlp[4] + $hlp[5]; 585 585 } 586 # REC_CLK in BNC 2.13 587 elsif ( $ln =~ /\bREC_CLK G\b/ ) { 588 push ( @EPOCHs_G_CLK, $epochSec ); 589 push ( @G_CLKs, $hlp[3] + $hlp[4] ); 590 } 591 elsif ( $ln =~ /\bREC_CLK R\b/ ) { 592 push ( @EPOCHs_R_CLK, $epochSec ); 593 push ( @R_CLKs, $hlp[3] + $hlp[4] ); 594 } 595 elsif ( $ln =~ /\bREC_CLK E\b/ ) { 596 push ( @EPOCHs_E_CLK, $epochSec ); 597 push ( @E_CLKs, $hlp[3] + $hlp[4] ); 598 } 599 elsif ( $ln =~ /\bREC_CLK C\b/ ) { 600 push ( @EPOCHs_C_CLK, $epochSec ); 601 push ( @C_CLKs, $hlp[3] + $hlp[4] ); 586 elsif ( $ln =~ /\bREC_CLK\b/ ) { 587 push ( @EPOCHs_CLK, $epochSec ); 588 push ( @CLKs, $hlp[3] + $hlp[4] ); 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] ); 602 601 } 603 602 elsif ( $ln =~ /\bSATNUM\b/ ) { # 2015-09... SATNUM G 8 … … 607 606 push ( @TRPs, $hlp[2] + $hlp[3] ); 608 607 } 609 610 608 } # ----- next line ----- 611 609 … … 622 620 . ", TRPs:" 623 621 . scalar @TRPs ); 624 if ( $nof_epochs != scalar @N ) { LOGDIE "number of epochs and residuals not equal\n" }625 if ( $nof_epochs != scalar @TRPs ) { LOGDIE "number of epochs and TRPs not equal\n" }626 if ( @ G_CLKs && scalar @EPOCHs_G_CLK != scalar @G_CLKs ) { LOGDIE "number of epochs and G_CLKs not equal\n" }627 if ( @ R_CLKs && scalar @EPOCHs_R_CLK != scalar @R_CLKs ) { LOGDIE "number of epochs and R_CLKs not equal\n" }628 if ( @ E_CLKs && scalar @EPOCHs_E_CLK != scalar @E_CLKs ) { LOGDIE "number of epochs and E_CLKs not equal\n" }629 if ( @ C_CLKs && scalar @EPOCHs_C_CLK != scalar @C_CLKs ) { LOGDIE "number of epochs and C_CLKs not equal\n" }622 if ( $nof_epochs != scalar @N ) { LOGDIE "number of epochs and residuals not equal\n" } 623 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 if ( @OFF_GLOs && scalar @EPOCHs_OFF_GLO != scalar @OFF_GLOs ) { LOGDIE "number of epochs and OFF_GLOs not equal\n" } 626 if ( @OFF_GALs && scalar @EPOCHs_OFF_GAL != scalar @OFF_GALs ) { LOGDIE "number of epochs and OFF_GALs not equal\n" } 627 if ( @OFF_BDSs && scalar @EPOCHs_OFF_BDS != scalar @OFF_BDSs ) { LOGDIE "number of epochs and OFF_BDSs not equal\n" } 630 628 631 629 my %data = ( … … 636 634 SATNUM => \%SATNUM, 637 635 TRPs => \@TRPs, 638 G_CLKs => \@G_CLKs,639 R_CLKs => \@R_CLKs,640 E_CLKs => \@E_CLKs,641 C_CLKs => \@C_CLKs,636 CLKs => \@CLKs, 637 OGRs => \@OFF_GLOs, 638 OGEs => \@OFF_GALs, 639 OGCs => \@OFF_BDSs, 642 640 RES => \%RES, 643 641 AMB => \%AMB, -
trunk/BNC/scripts/pppPlot.pl
r9947 r10005 8 8 # - NEU displacements w.r.t. coordinates in Crd file 9 9 # - A-priori + correction values of tropospheric zenith delay in [m], 10 # - GNSS Receiver clock errors in [m],10 # - Receiver clock error and offsets in [m], 11 11 # - Elevations, given per satellite 12 12 # - Ambiguities, given per satellite … … 106 106 # ----------------------------------------------------------------------------- 107 107 my ( $station, $file ) = Bnc::parsePPPLogfile( $file, $sampling ); 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'};113 my %AMB = %{ $file->{'AMB'} };114 my %RES = %{ $file->{'RES'} };115 my %ELE = %{ $file->{'ELE'} };116 my %ION = %{ $file->{'ION'} };117 my %BIA = %{ $file->{'BIA'} };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'}; 113 my %AMB = %{ $file->{'AMB'} }; 114 my %RES = %{ $file->{'RES'} }; 115 my %ELE = %{ $file->{'ELE'} }; 116 my %ION = %{ $file->{'ION'} }; 117 my %BIA = %{ $file->{'BIA'} }; 118 118 119 119 # ----------------------------------------------------------------------------- … … 212 212 $png->image( $pdf->image_png($pngName), $x, $y, $width, $height ); 213 213 214 ######### GPSCLK #####################215 if ( scalar @{ $file->{' G_CLKs'} } < 1 ) {216 DEBUG "No GPSCLKs found";214 ######### CLK ##################### 215 if ( scalar @{ $file->{'CLKs'} } < 1 ) { 216 DEBUG "No CLKs found"; 217 217 } 218 218 else { 219 DEBUG "Plot GPSCLKs";220 $pngName = sprintf ( "%s_CLK _G.png", $station );219 DEBUG "Plot CLKs"; 220 $pngName = sprintf ( "%s_CLK.png", $station ); 221 221 my $chartCLK = newChart($station); 222 222 $chartCLK->set( output => $pngName ); … … 224 224 225 225 $dataset = Chart::Gnuplot::DataSet->new( 226 xdata => $EPOCHS _G_CLK,227 ydata => $file->{' G_CLKs'},228 title => " GPSReceiver clock",226 xdata => $EPOCHS, 227 ydata => $file->{'CLKs'}, 228 title => "Receiver clock", 229 229 timefmt => '%s', 230 230 style => "linespoints", … … 242 242 $png->image( $pdf->image_png($pngName), $x, $y, $width, $height ); 243 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 => "linespoints", 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("could not 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 => "linespoints", 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("could not 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 => "linespoints", 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("could not find image file: $!\n") unless -e $pngName; 333 $png->image( $pdf->image_png($pngName), $x, $y, $width, $height ); 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 } 334 343 } 335 344
Note:
See TracChangeset
for help on using the changeset viewer.
