Changeset 9910 in ntrip
- Timestamp:
- Nov 24, 2022, 10:29:49 PM (23 months ago)
- Location:
- trunk/BNC/scripts
- Files:
-
- 2 edited
-
Bnc.pm (modified) (10 diffs)
-
pppPlot.pl (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/scripts/Bnc.pm
r9597 r9910 488 488 $startSec = time () - $goBackSecs; 489 489 } 490 my $isFreePPP = 1; # Free or closed PPP version491 490 my $epo; 492 491 my $old_epochSec = 0; … … 494 493 my $epochDiff = 0; 495 494 my ( 496 @hlp, @EPOCHs, @N, @E, @U, %SATNUM, @TRPs, @CLKs,497 @G_CLKs, @R_CLKs, @E_CLKs, @C_CLKs, @OGRs, @OGEs, @OGCs, @OFFGLOs495 @hlp, @EPOCHs, @N, @E, @U, %SATNUM, @TRPs, 496 @G_CLKs, @R_CLKs, @E_CLKs, @C_CLKs, 498 497 ); 499 my ( @EPOCHs_OGE, @EPOCHs_OGR, @EPOCHs_OGC );500 498 my ( @EPOCHs_G_CLK, @EPOCHs_R_CLK, @EPOCHs_E_CLK, @EPOCHs_C_CLK ); 501 499 my ( %AMB, %RES, %ELE, %ION, %BIA ); … … 516 514 517 515 # PPP of Epoch 2015-08-27_14:00:15.000 518 if ( $ln =~ /PPP of Epoch (\d{4}-\d{2}-\d{2}_\d{2}:\d{2}:\d{2})\.\d+/ ) { # closed PPP 519 $isFreePPP = 0; 520 $epo = $1; #print "$epo\n"; 521 } 522 elsif ( $ln =~ /Point Positioning of Epoch (\d{4}-\d{2}-\d{2}_\d{2}:\d{2}:\d{2})\.\d+/ ) { # free PPP 523 $isFreePPP = 1; 524 $epo = $1; 516 if ( $ln =~ /PPP of Epoch (\d{4}-\d{2}-\d{2}_\d{2}:\d{2}:\d{2})\.\d+/ ) { 517 $epo = $1; #print "$epo\n"; 525 518 } 526 519 else { ERROR "strange line: \"$ln\""; next } 527 520 528 #my $date = sprintf ( "%s %s", split ( /_/, $epo ) );529 521 my $tp = Time::Piece->strptime( $epo, '%Y-%m-%d_%H:%M:%S' ); 530 522 $epochSec = $tp->epoch(); … … 539 531 @hlp = split ( /\s+/, $ln ); 540 532 541 if ( $ln =~ /\bOFFGLO\b/ ) { # ... OFFGLO 8.417 +- 28.479 542 push ( @OFFGLOs, $hlp[2] ); 543 } 544 545 elsif ( $ln =~ /\bdN\b/ ) { 533 if ( $ln =~ /\bdN\b/ ) { 546 534 push ( @EPOCHs, $epochSec ); # besser $epo ? 547 535 $old_epochSec = $epochSec; 548 536 549 #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.0349537 #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 550 538 $station = $hlp[1]; 551 539 … … 559 547 } 560 548 elsif ( ( $ln =~ /\bAMB\b/ ) && ( $ln !~ /RESET/ ) ) { 561 if ($isFreePPP) { # 2015-10... AMB G05 -6.754 +- 0.086 nEpo = 633 562 $sat = $hlp[2]; 563 $sys = substr ( $sat, 0, 1 ); 564 $amb = $hlp[4]; 565 push @{ $AMB{$sys}{$sat}{EPOCH} }, $epochSec; 566 push @{ $AMB{$sys}{$sat}{DATA} }, $amb; 567 push @{ $AMB{$sys}{$sat}{NUMEPO} }, $hlp[9]; 568 } 569 else { # 2015-08... AMB lIF G04 253.0000 -8.9924 +- 1.7825 el = 22.03 epo = 86 570 $lki = $hlp[2]; 571 $sat = $hlp[3]; 572 $sys = substr ( $sat, 0, 1 ); 573 $amb = $hlp[4] + $hlp[5]; 574 push @{ $AMB{$lki}{$sys}{$sat}{EPOCH} }, $epochSec; 575 push @{ $AMB{$lki}{$sys}{$sat}{DATA} }, $amb; 576 push @{ $AMB{$lki}{$sys}{$sat}{NUMEPO} }, $hlp[13]; 577 push @{ $ELE{$sys}{$sat}{EPOCH} }, $epochSec; 578 push @{ $ELE{$sys}{$sat}{DATA} }, $hlp[10]; 579 } 549 # 2015-08... AMB lIF G04 253.0000 -8.9924 +- 1.7825 el = 22.03 epo = 86 550 $lki = $hlp[2]; 551 $sat = $hlp[3]; 552 $sys = substr ( $sat, 0, 1 ); 553 $amb = $hlp[4] + $hlp[5]; 554 push @{ $AMB{$lki}{$sys}{$sat}{EPOCH} }, $epochSec; 555 push @{ $AMB{$lki}{$sys}{$sat}{DATA} }, $amb; 556 push @{ $AMB{$lki}{$sys}{$sat}{NUMEPO} }, $hlp[13]; 557 push @{ $ELE{$sys}{$sat}{EPOCH} }, $epochSec; 558 push @{ $ELE{$sys}{$sat}{DATA} }, $hlp[10]; 580 559 } 581 560 elsif ( $ln =~ /\bRES\b/ && $ln !~ /Neglected/ ) { 582 if ($isFreePPP) { # 2015-10... RES R08 L3 -0.0069 583 $sat = $hlp[2]; 584 $lki = $hlp[3]; 585 } 586 else { # 2015-08... RES lIF G30 -0.0076 587 $sat = $hlp[3]; 588 $lki = $hlp[2]; 589 } 561 # 2015-08... RES lIF G30 -0.0076 562 $sat = $hlp[3]; 563 $lki = $hlp[2]; 590 564 $sys = substr ( $sat, 0, 1 ); 591 565 … … 610 584 push @{ $BIA{$lki}{$sys}{DATA} }, $hlp[4] + $hlp[5]; 611 585 } 612 elsif ( $ln =~ /\bCLK\b/ ) {613 if ($isFreePPP) { push ( @CLKs, $hlp[2] ) }614 else { push ( @CLKs, $hlp[2] + $hlp[3] ) }615 }616 617 586 # REC_CLK in BNC 2.13 618 #elsif ( $ln =~ /\bREC_CLK\b/ ) {619 elsif ( $ln =~ /\bREC_CLK\s{8}/ ) {620 push ( @CLKs, $hlp[2] + $hlp[3] );621 }622 587 elsif ( $ln =~ /\bREC_CLK G\b/ ) { 623 588 push ( @EPOCHs_G_CLK, $epochSec ); … … 635 600 push ( @EPOCHs_C_CLK, $epochSec ); 636 601 push ( @C_CLKs, $hlp[3] + $hlp[4] ); 637 }638 elsif ( $ln =~ /\bOGR\b/ ) { # 2015-08... OGR 52.6806 -3.8042 +- 9.0077639 push ( @EPOCHs_OGR, $epochSec );640 push ( @OGRs, $hlp[2] + $hlp[3] ); # only https so far641 }642 elsif ( $ln =~ /\bOGE\b/ ) { # 2015-08... OGE 52.6806 -3.8042 +- 9.0077643 push ( @EPOCHs_OGE, $epochSec );644 push ( @OGEs, $hlp[2] + $hlp[3] ); # only https so far645 }646 elsif ( $ln =~ /\bOGC\b/ ) { # 2015-08... OGC 52.6806 -3.8042 +- 9.0077647 push ( @EPOCHs_OGC, $epochSec );648 push ( @OGCs, $hlp[2] + $hlp[3] ); # only https so far649 602 } 650 603 elsif ( $ln =~ /\bSATNUM\b/ ) { # 2015-09... SATNUM G 8 … … 668 621 . scalar @U 669 622 . ", TRPs:" 670 . scalar @TRPs 671 . ", CLKs:" 672 . scalar @CLKs 673 . ", OFFGLOs:" 674 . scalar @OFFGLOs ); 623 . scalar @TRPs ); 675 624 if ( $nof_epochs != scalar @N ) { LOGDIE "number of epochs and residuals not equal\n" } 676 625 if ( $nof_epochs != scalar @TRPs ) { LOGDIE "number of epochs and TRPs not equal\n" } 677 if ( @CLKs && $nof_epochs != scalar @CLKs ) { LOGDIE "number of epochs and CLKs not equal\n" }678 626 if ( @G_CLKs && scalar @EPOCHs_G_CLK != scalar @G_CLKs ) { LOGDIE "number of epochs and G_CLKs not equal\n" } 679 627 if ( @R_CLKs && scalar @EPOCHs_R_CLK != scalar @R_CLKs ) { LOGDIE "number of epochs and R_CLKs not equal\n" } 680 628 if ( @E_CLKs && scalar @EPOCHs_E_CLK != scalar @E_CLKs ) { LOGDIE "number of epochs and E_CLKs not equal\n" } 681 629 if ( @C_CLKs && scalar @EPOCHs_C_CLK != scalar @C_CLKs ) { LOGDIE "number of epochs and C_CLKs not equal\n" } 682 if ( @OGRs && scalar @EPOCHs_OGR != scalar @OGRs ) { LOGDIE "number of epochs and OGRs not equal\n" }683 if ( @OGEs && scalar @EPOCHs_OGE != scalar @OGEs ) { LOGDIE "number of epochs and OGEs not equal\n" }684 if ( @OGCs && scalar @EPOCHs_OGC != scalar @OGCs ) { LOGDIE "number of epochs and OGCs not equal\n" }685 630 686 631 my %data = ( … … 691 636 SATNUM => \%SATNUM, 692 637 TRPs => \@TRPs, 693 CLKs => \@CLKs,694 638 G_CLKs => \@G_CLKs, 695 639 R_CLKs => \@R_CLKs, 696 640 E_CLKs => \@E_CLKs, 697 641 C_CLKs => \@C_CLKs, 698 OGRs => \@OGRs,699 OGEs => \@OGEs,700 OGCs => \@OGCs,701 OFFGLOs => \@OFFGLOs,702 642 RES => \%RES, 703 643 AMB => \%AMB, … … 707 647 ); 708 648 709 return ( $station, \%data, $isFreePPP);649 return ( $station, \%data, 0 ); 710 650 } 711 651 -
trunk/BNC/scripts/pppPlot.pl
r9836 r9910 192 192 my $chartTRP = newChart($station); 193 193 $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 " ] ); 195 196 196 197 my $dataTRP = Chart::Gnuplot::DataSet->new( … … 213 214 $png->image( $pdf->image_png($pngName), $x, $y, $width, $height ); 214 215 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 }245 216 ######### GPS CLK ##################### 246 217 if ( scalar @{ $file->{'G_CLKs'} } < 1 ) { … … 363 334 LOGDIE("could not find image file: $!\n") unless -e $pngName; 364 335 $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 }464 336 } 465 337 … … 647 519 648 520 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 " ] ); 650 523 } 651 524 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 " ] ); 653 527 } 654 528 … … 787 661 $pdf->end(); 788 662 789 #system ("rm *.png");663 system ("rm *.png"); 790 664 if (Common::amInteractiv ) { 791 system ("evince $inputDir/$pdf_name");665 # system ("evince $inputDir/$pdf_name"); 792 666 } 793 667 } # ----- next logfile -----
Note:
See TracChangeset
for help on using the changeset viewer.
