Changeset 10233 in ntrip
- Timestamp:
- Nov 1, 2023, 9:02:50 AM (6 months ago)
- Location:
- trunk/BNC/scripts
- Files:
-
- 2 edited
-
Bnc.pm (modified) (4 diffs)
-
pppPlot.pl (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/scripts/Bnc.pm
r10090 r10233 493 493 my $epochSec = 0; 494 494 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); 497 497 my ( %AMB, %RES, %ELE, %ION, %BIA ); 498 498 my ( $station, $lki, $sys, $sat, $amb ); … … 584 584 push @{ $BIA{$lki}{$sys}{DATA} }, $hlp[4] + $hlp[5]; 585 585 } 586 elsif ( $ln =~ /\bREC_CLK \b/ ){587 push ( @EPOCHs_ CLK, $epochSec);588 push ( @CLKs, $hlp[2] + $hlp[3] );589 } 590 elsif ( $ln =~ /\b OFF_GLO\b/ ) { # 2015-08... OFF_GLO 52.6806 -3.8042 +- 9.0077591 push ( @EPOCHs_ OFF_GLO, $epochSec);592 push ( @OFF_GLOs, $hlp[2] + $hlp[3] );593 } 594 elsif ( $ln =~ /\b OFF_GAL\b/ ) { # 2015-08... OFF_GAL 52.6806 -3.8042 +- 9.0077595 push ( @EPOCHs_ OFF_GAL, $epochSec);596 push ( @OFF_GALs, $hlp[2] + $hlp[3] );597 } 598 elsif ( $ln =~ /\b OFF_BDS\b/ ) { # 2015-08... OFF_BDS 52.6806 -3.8042 +- 9.0077599 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] ); 601 601 } 602 602 elsif ( $ln =~ /\bSATNUM\b/ ) { # 2015-09... SATNUM G 8 … … 619 619 . scalar @U 620 620 . ", 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 ); 622 631 if ( $nof_epochs != scalar @N ) { LOGDIE "number of epochs and residuals not equal\n" } 623 632 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";} 635 637 if ( !$station ) { WARN "could not grep stationname from file: $file\n" } 636 638 … … 642 644 SATNUM => \%SATNUM, 643 645 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, 649 650 RES => \%RES, 650 651 AMB => \%AMB, -
trunk/BNC/scripts/pppPlot.pl
r10197 r10233 106 106 # ----------------------------------------------------------------------------- 107 107 my ( $station, $file ) = Bnc::parsePPPLogfile( $file, $sampling ); 108 my $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_C LK = $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'}; 113 113 my %AMB = %{ $file->{'AMB'} }; 114 114 my %RES = %{ $file->{'RES'} }; … … 212 212 $png->image( $pdf->image_png($pngName), $x, $y, $width, $height ); 213 213 214 ######### CLK #####################215 if ( scalar @{ $file->{' CLKs'} } < 1 ) {216 DEBUG "NoCLKs found";214 ######### GPS CLK ##################### 215 if ( scalar @{ $file->{'G_CLKs'} } < 1 ) { 216 DEBUG "No GPS CLKs found"; 217 217 } 218 218 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 ); 343 334 } 344 335
Note:
See TracChangeset
for help on using the changeset viewer.
