Changeset 9910 in ntrip for trunk/BNC/scripts/Bnc.pm
- Timestamp:
- Nov 24, 2022, 10:29:49 PM (17 months ago)
- File:
-
- 1 edited
-
trunk/BNC/scripts/Bnc.pm (modified) (10 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
Note:
See TracChangeset
for help on using the changeset viewer.
