Index: trunk/BNC/src/bnchelp.html =================================================================== --- trunk/BNC/src/bnchelp.html (revision 10091) +++ trunk/BNC/src/bnchelp.html (revision 10092) @@ -2390,62 +2390,83 @@ 


-! SP3 File 1: Input/CNE1SSRRTS_20222410000_01D_05S_ORB.SP3
-! SP3 File 2: Input/BKG1SSRRTS_20222410000_01D_05S_ORB.SP3
+! SP3 File 1: /data/stuerze/tools/BNC/Example_Configs/Input/CNE1MGXRTS_20222410000_01D_05S_ORB.SP3
+! SP3 File 2: /data/stuerze/tools/BNC/Example_Configs/Input/BKG1MGXRTS_20222410000_01D_05S_ORB.SP3
+!
+! Clock residuals and orbit differences in [m]
+! ----------------------------------------------------------------------------
 !
 !  Epoch                PRN  radial   along   out        clk    clkRed   iPRN
 ! ----------------------------------------------------------------------------
-..
-2022-08-29_08:01:30.000 E02 -0.0023  0.1159 -0.0047     0.0187  0.0210    60
-2022-08-29_08:01:30.000 E03 -0.0567 -0.0411  0.0146    -0.0093  0.0474    61
-2022-08-29_08:01:30.000 E04  0.0695 -0.0406  0.0880    -0.0219 -0.0914    62
-2022-08-29_08:01:30.000 E05 -0.1042  0.1161 -0.0503    -0.0101  0.0942    63
-2022-08-29_08:01:30.000 E07  0.0354  0.0384 -0.0346     0.0134 -0.0219    65
-2022-08-29_08:01:30.000 E08 -0.0312  0.0312 -0.0481     0.0130  0.0442    66
-2022-08-29_08:01:30.000 E09 -0.0548  0.1941  0.0263    -0.0113  0.0436    67
-2022-08-29_08:01:30.000 E12  0.0402 -0.0023  0.0012    -0.0051 -0.0454    70
-2022-08-29_08:01:30.000 E13  0.0420 -0.0608 -0.0525    -0.0006 -0.0427    71
-2022-08-29_08:01:30.000 E15  0.0478  0.0348 -0.0353    -0.0122 -0.0600    73
-2022-08-29_08:01:30.000 E21  0.0639  0.0033  0.0292    -0.0167 -0.0806    79
-2022-08-29_08:01:30.000 E24  0.0900 -0.1113 -0.0385     0.0047 -0.0853    82
-2022-08-29_08:01:30.000 E25  0.0106  0.0340 -0.0564     0.0023 -0.0083    83
-2022-08-29_08:01:30.000 E26  0.0598  0.0506 -0.0983    -0.0011 -0.0609    84
-2022-08-29_08:01:30.000 E27  0.0763 -0.0972 -0.0480     0.0041 -0.0722    85
-2022-08-29_08:01:30.000 E30  0.0385  0.0266  0.0028     0.0092 -0.0293    88
-2022-08-29_08:01:30.000 E31  0.1017 -0.1834 -0.0211     0.0021 -0.0997    89
-2022-08-29_08:01:30.000 E33  0.0827  0.0385 -0.0741    -0.0060 -0.0887    91
-2022-08-29_08:01:30.000 E36  0.0099 -0.0769  0.0384     0.0269  0.0170    94
-
+2022-08-29_08:01:30.000 G01  0.0297 -0.0526 -0.0198    -0.0005 -0.0301     1
+2022-08-29_08:01:30.000 G02  0.0225 -0.0179  0.0055     0.0054 -0.0171     2
+2022-08-29_08:01:30.000 G03  0.0210 -0.0387 -0.0153     0.0025 -0.0185     3
+2022-08-29_08:01:30.000 G06  0.0650  0.0599 -0.0260    -0.0074 -0.0724     6
+2022-08-29_08:01:30.000 G07  0.0112  0.0774  0.0331    -0.0021 -0.0133     7
+2022-08-29_08:01:30.000 G08  0.0241  0.0725  0.0009    -0.0008 -0.0248     8
+2022-08-29_08:01:30.000 G09 -0.0145  0.0669  0.0069     0.0051  0.0196     9
+2022-08-29_08:01:30.000 G11  0.0375  0.0650 -0.0080    -0.0094 -0.0469    11
+2022-08-29_08:01:30.000 G12  0.0634  0.0229 -0.0003     0.0078 -0.0556    12
+2022-08-29_08:01:30.000 G13  0.0163  0.0551  0.0904     0.0010 -0.0153    13
+2022-08-29_08:01:30.000 G14  0.0365  0.0646 -0.0398     0.0034 -0.0332    14
+2022-08-29_08:01:30.000 G15  0.0175  0.0892  0.0512    -0.0032 -0.0207    15
+2022-08-29_08:01:30.000 G16  0.0453  0.0408  0.0346    -0.0029 -0.0483    16
+2022-08-29_08:01:30.000 G17  0.0381 -0.0064 -0.0701     0.0020 -0.0362    17
+2022-08-29_08:01:30.000 G18 -0.0886  0.0031  0.0721    -0.0055  0.0831    18
+2022-08-29_08:01:30.000 G19  0.0310 -0.0190 -0.0264    -0.0000 -0.0311    19
+2022-08-29_08:01:30.000 G20  0.0199  0.0308  0.0476     0.0106 -0.0093    20
+2022-08-29_08:01:30.000 G21  0.0163  0.0079 -0.0464     0.0047 -0.0116    21
+2022-08-29_08:01:30.000 G22  0.0471  0.0173 -0.0315    -0.0023 -0.0494    22
+2022-08-29_08:01:30.000 G23 -0.0354 -0.0042  0.0247    -0.0035  0.0319    23
+2022-08-29_08:01:30.000 G24  0.0110  0.0030 -0.0323     0.0112  0.0002    24
+2022-08-29_08:01:30.000 G25  0.0037  0.0712 -0.0175    -0.0019 -0.0056    25
+2022-08-29_08:01:30.000 G26  0.0421 -0.0002  0.0511    -0.0009 -0.0430    26
+2022-08-29_08:01:30.000 G27  0.0164  0.1203  0.0347    -0.0041 -0.0205    27
+2022-08-29_08:01:30.000 G29 -0.0093  0.0141  0.0570     0.0095  0.0188    29
+2022-08-29_08:01:30.000 G30  0.0199  0.0805  0.0529    -0.0074 -0.0273    30
+2022-08-29_08:01:30.000 G31  0.0336 -0.0079 -0.0281    -0.0132 -0.0468    31
+2022-08-29_08:01:30.000 G32 -0.0185 -0.0200 -0.0103     0.0019  0.0205    32
 ..
 !
-! RMS[m]
+! Summary
+! -----------------------------------------------------------------------------------------------------------------
 !
-!    PRN  radial   along   out     nOrb    clk   clkRed   nClk    Offset
-! ----------------------------------------------------------------------
-!    E02  0.0072  0.1160  0.0267    165  0.0110  0.0158    165   -0.0345
-!    E03  0.0610  0.0387  0.0151    165  0.0083  0.0617    165   -0.0376
-!    E04  0.0830  0.0619  0.0819    165  0.0112  0.0828    165    0.0649
-!    E05  0.0966  0.1193  0.0418    165  0.0058  0.0965    165   -0.0073
-!    E07  0.0296  0.0309  0.0389    165  0.0102  0.0296    165   -0.0904
-!    E08  0.0371  0.0193  0.0424    165  0.0080  0.0372    165   -0.1017
-!    E09  0.0436  0.1744  0.0422    165  0.0044  0.0433    165   -0.0300
-!    E12  0.0490  0.0033  0.0107    165  0.0042  0.0489    165   -0.1174
-!    E13  0.0336  0.0393  0.0661    165  0.0046  0.0341    165   -0.0683
-!    E15  0.0548  0.0273  0.0375    165  0.0064  0.0547    165    0.0536
-!    E21  0.0667  0.0116  0.0301    165  0.0079  0.0671    165   -0.0384
-!    E24  0.0913  0.1101  0.0332    165  0.0050  0.0913    165    0.0946
-!    E25  0.0099  0.0412  0.0628    165  0.0033  0.0101    165    0.0797
-!    E26  0.0621  0.0516  0.0946    165  0.0033  0.0622    165    0.0300
-!    E27  0.0734  0.0841  0.0522    165  0.0052  0.0733    165    0.1015
-!    E30  0.0382  0.0380  0.0072    165  0.0046  0.0384    165    0.0682
-!    E31  0.0937  0.1700  0.0234    165  0.0089  0.0943    165    0.1110
-!    E33  0.0926  0.0340  0.0654    165  0.0032  0.0926    165    0.0010
-!    E36  0.0059  0.0630  0.0456    165  0.0089  0.0073    165   -0.0790
+!     PRN   radialRMS   alongRMS    outRMS     3DRMS      nOrb    clkRMS    clkRedRMS   clkRedSig   nClk    Offset
+!              [mm]       [mm]        [mm]      [mm]      [-]      [ns]        [ns]        [ns]      [-]     [ns]
+! -----------------------------------------------------------------------------------------------------------------
 ..
+!     G01       49.4       90.1        4.0      102.8     165       0.01       0.16       0.16       165     -2.71
+!     G02       39.8       34.2        6.2       52.8     165       0.01       0.13       0.13       165      4.13
+!     G03       27.8       57.7       12.8       65.3     165       0.01       0.09       0.09       165     -2.18
+!     G06       61.8       50.7       22.2       82.9     165       0.01       0.21       0.21       165     -2.70
+!     G07       22.2       80.3       31.7       89.2     165       0.01       0.07       0.07       165     -1.26
+!     G08       18.8      102.1       12.7      104.6     165       0.01       0.06       0.06       165      1.17
+!     G09        9.8       87.2        5.6       87.9     165       0.01       0.03       0.03       165      0.21
+!     G11       54.8       65.5        6.9       85.6     165       0.03       0.18       0.18       165     -1.93
+!     G12       64.4       48.2        7.3       80.8     165       0.02       0.21       0.21       165     -1.45
+!     G13       17.1       45.9       88.8      101.4     165       0.01       0.06       0.06       165     -0.67
+!     G14       34.9       79.4       32.6       92.7     165       0.01       0.12       0.12       165     -2.07
+!     G15       33.4       38.1       45.0       67.8     165       0.01       0.11       0.11       165     -2.19
+!     G16       50.0       67.9       51.2       98.7     165       0.01       0.17       0.17       165      3.27
+!     G17       44.1       12.0       73.8       86.8     165       0.01       0.15       0.15       165     -1.67
+!     G18       78.7       15.0       65.4      103.4     165       0.01       0.26       0.26       165     -2.37
+!     G19       29.1       12.7       38.2       49.7     165       0.01       0.10       0.10       165      6.27
+!     G20       19.0       25.2       33.4       45.9     165       0.02       0.06       0.06       165      5.58
+!     G21        8.6       13.3       44.5       47.3     165       0.01       0.03       0.03       165      5.40
+!     G22       53.5       18.3       29.0       63.5     165       0.01       0.18       0.18       165      1.54
+!     G23       29.8       35.8       16.1       49.2     165       0.01       0.10       0.10       165     -1.84
+!     G24       13.4       20.2       28.4       37.3     165       0.02       0.04       0.04       165     -2.67
+!     G25        9.2       88.4       28.0       93.2     165       0.02       0.04       0.04       165      1.68
+!     G26       45.2       33.4       60.5       82.5     165       0.02       0.15       0.15       165      0.36
+!     G27       22.4      128.4       39.9      136.3     165       0.01       0.07       0.07       165      0.59
+!     G29       17.2       42.3       59.7       75.1     165       0.02       0.06       0.06       165     -2.56
+!     G30       19.9       91.7       61.2      112.0     165       0.01       0.07       0.07       165      2.15
+!     G31       30.3       20.1       39.4       53.6     165       0.04       0.11       0.11       165     -1.97
+!     G32       16.7       28.9        8.2       34.4     165       0.01       0.06       0.06       165     -2.11
 !
-!  Total  0.0552  0.0877  0.0690  10395  0.0075  0.0551  10395
+!   Total       55.2       87.7       69.0      124.5   10395       0.02       0.18       0.18     10395

-The first part of this output uses the following abbreviations: +The first part of this output following string 'following string' uses the following abbreviations:

@@ -2461,17 +2482,21 @@

-The second part following string 'RMS' provides a summary of the comparison using the following +The second part following string 'Summary' provides a summary of the comparison using the following abbreviations:
+! PRN radialRMS alongRMS outRMS 3DRMS nOrb clkRMS clkRedRMS clkRedSig nClk Offset +! [mm] [mm] [mm] [mm] [-] [ns] [ns] [ns] [-] [ns]

- - - + + + + - - + + + - +
'PRN'  Satellite specification
'radial'  RMS of radial component of orbit coordinate differences [m]
'along'  RMS of along track component of orbit coordinate differences [m]
'out'  RMS of out-of-plane component of orbit coordinate differences [m]
'radialRMS'  RMS of the radial component of orbit coordinate differences [mm]
'alongRMS'  RMS of the along track component of orbit coordinate differences [mm]
'outRMS'  RMS of the out-of-plane component of orbit coordinate differences [mm]
'3DRMS'  3D RMS of the orbit coordinate differences [mm]
'nOrb'  Number of epochs used in in orbit comparison
'clk'  RMS of clock differences [m]
'clkRed'  RMS of clock differences after reduction of radial orbit differences [m]
'clkRMS'  RMS of clock differences [ns]
'clkRedRMS'  RMS of The clock differences after reduction of radial orbit differences [ns]
'clkRedSig'  Sigma of the clock differences after reduction of radial orbit differences [ns]
'nClk'  Number of epochs use in clock comparisons
'Offset'  Clock offset [m]
'Offset'  Clock offset [ns]

Index: trunk/BNC/src/orbComp/sp3Comp.cpp =================================================================== --- trunk/BNC/src/orbComp/sp3Comp.cpp (revision 10091) +++ trunk/BNC/src/orbComp/sp3Comp.cpp (revision 10092) @@ -217,13 +217,32 @@ // Compute clock residuals // ----------------------- + const string all = "ZZZ"; for (unsigned ie = 0; ie < epochs.size(); ie++) { - map& dc = epochs[ie]->_dc; + map& dc = epochs[ie]->_dc; + map& dcRed = epochs[ie]->_dcRed; + const map& dr = epochs[ie]->_dr; for (map::iterator it = dc.begin(); it != dc.end(); it++) { const t_prn& prn = it->first; if (satIndex(clkSats, prn) != -1) { int index = epochs.size() + satIndex(clkSats, prn); - dc[prn] = it->second - xx[ie] - xx[index]; + dc[prn] = it->second - xx[ie] - xx[index]; stat[prn.toString()]._offset = xx[index]; - } + if (dr.find(prn) != dr.end()){ + dcRed[prn] = dc[prn] - dr.find(prn)->second[0]; // clock minus radial component + stat[prn.toString()]._dcRedMean += dcRed[prn]; + stat[all ]._dcRedMean += dcRed[prn]; + stat[prn.toString()]._nc += 1; + stat[all ]._nc += 1; + } + } + } + } + + // Compute Clock Mean + // ------------------ + for (map::iterator it = stat.begin(); it != stat.end(); it++) { + t_stat& stat = it->second; + if (stat._nc > 0) { + stat._dcRedMean = stat._dcRedMean / stat._nc; } } @@ -342,9 +361,11 @@ } - // Print Residuals - // --------------- + // Print epoch-wise Clock Residuals + // -------------------------------- const string all = "ZZZ"; out.setf(ios::fixed); + out << "!\n! Clock residuals and orbit differences in [m]\n" + "! ----------------------------------------------------------------------------\n"; out << "!\n! Epoch PRN radial along out clk clkRed iPRN" "\n! ----------------------------------------------------------------------------\n"; @@ -353,4 +374,5 @@ const map& dr = epochs[ii]->_dr; const map& dc = epochs[ii]->_dc; + const map& dcRed = epochs[ii]->_dcRed; for (map::const_iterator it = dr.begin(); it != dr.end(); it++) { const t_prn& prn = it->first; @@ -362,18 +384,20 @@ << setw(7) << setprecision(4) << rao[2] << " "; stat[prn.toString()]._rao += SP(rao, rao); // Schur product + stat[all ]._rao += SP(rao, rao); stat[prn.toString()]._nr += 1; - stat[all]._rao += SP(rao, rao); - stat[all]._nr += 1; - if (dc.find(prn) != dc.end()) { + stat[all ]._nr += 1; + if (dc.find(prn) != dc.end() && dcRed.find(prn) != dc.end()) { double clkRes = dc.find(prn)->second; - double clkResRed = clkRes - it->second[0]; // clock minus radial component + double clkResRed = dcRed.find(prn)->second; out << setw(7) << setprecision(4) << clkRes << ' ' << setw(7) << setprecision(4) << clkResRed; - stat[prn.toString()]._dc += clkRes * clkRes; - stat[prn.toString()]._dcRed += clkResRed * clkResRed; - stat[prn.toString()]._nc += 1; - stat[all]._dc += clkRes * clkRes; - stat[all]._dcRed += clkResRed * clkResRed; - stat[all]._nc += 1; + stat[prn.toString()]._dcRMS += clkRes * clkRes; + stat[all ]._dcRMS += clkRes * clkRes; + stat[prn.toString()]._dcRedRMS += clkResRed * clkResRed; + stat[all ]._dcRedRMS += clkResRed * clkResRed; + stat[prn.toString()]._dcRedSig += (clkResRed - stat[prn.toString()]._dcRedMean) * + (clkResRed - stat[prn.toString()]._dcRedMean); + stat[all ]._dcRedSig += (clkResRed - stat[all ]._dcRedMean) * + (clkResRed - stat[all ]._dcRedMean); } else { @@ -388,7 +412,9 @@ // Print Summary // ------------- - out << "!\n! RMS[m]\n"; - out << "!\n! PRN radial along out nOrb clk clkRed nClk Offset" - "\n! ----------------------------------------------------------------------\n"; + out << "!\n! Summary"; + out << "\n! -----------------------------------------------------------------------------------------------------------------\n"; + out << "!\n! PRN radialRMS alongRMS outRMS 3DRMS nOrb clkRMS clkRedRMS clkRedSig nClk Offset " + "\n! [mm] [mm] [mm] [mm] [-] [ns] [ns] [ns] [-] [ns] " + "\n! -----------------------------------------------------------------------------------------------------------------\n"; for (map::iterator it = stat.begin(); it != stat.end(); it++) { const string& prn = it->first; @@ -398,22 +424,28 @@ stat._rao[1] = sqrt(stat._rao[1] / stat._nr); stat._rao[2] = sqrt(stat._rao[2] / stat._nr); + stat._rao3DRMS = stat._rao.NormFrobenius(); if (prn == all) { - out << "!\n! Total "; + out << "!\n! Total "; } else { - out << "! " << prn << ' '; - } - out << setw(7) << setprecision(4) << stat._rao[0] << ' ' - << setw(7) << setprecision(4) << stat._rao[1] << ' ' - << setw(7) << setprecision(4) << stat._rao[2] << ' ' - << setw(6) << stat._nr << " "; + out << "! " << prn << ' '; + } + // orbit values in millimeters + out << setw(10) << setprecision(1) << stat._rao[0] * 1e3 << ' ' + << setw(10) << setprecision(1) << stat._rao[1] * 1e3 << ' ' + << setw(10) << setprecision(1) << stat._rao[2] * 1e3 << ' ' + << setw(10) << setprecision(1) << stat._rao3DRMS * 1e3 << ' ' + << setw( 7) << stat._nr << " "; if (stat._nc > 0) { - stat._dc = sqrt(stat._dc / stat._nc); - stat._dcRed = sqrt(stat._dcRed / stat._nc); - out << setw(7) << setprecision(4) << stat._dc << ' ' - << setw(7) << setprecision(4) << stat._dcRed << ' ' - << setw(6) << stat._nc << " "; + stat._dcRMS = sqrt(stat._dcRMS / stat._nc); + stat._dcRedRMS = sqrt(stat._dcRedRMS / stat._nc); + stat._dcRedSig = sqrt(stat._dcRedSig / stat._nc); + // clock values in nano seconds + out << setw(10) << setprecision(2) << stat._dcRMS / t_CST::c * 1e9 << ' ' + << setw(10) << setprecision(2) << stat._dcRedRMS / t_CST::c * 1e9 << ' ' + << setw(10) << setprecision(2) << stat._dcRedSig / t_CST::c * 1e9 << ' ' + << setw( 9) << stat._nc << " "; if (prn != all) { - out << setw(9) << setprecision(4) << stat._offset; + out << setw(9) << setprecision(2) << stat._offset / t_CST::c * 1e9; } } Index: trunk/BNC/src/orbComp/sp3Comp.h =================================================================== --- trunk/BNC/src/orbComp/sp3Comp.h (revision 10091) +++ trunk/BNC/src/orbComp/sp3Comp.h (revision 10092) @@ -38,5 +38,5 @@ class t_sp3Comp : public QThread { Q_OBJECT - + public: t_sp3Comp(QObject* parent); @@ -47,10 +47,10 @@ signals: void finished(); - + public slots: public: virtual void run(); - + private: class t_epoch { @@ -60,4 +60,5 @@ std::map _xyz; std::map _dc; + std::map _dcRed; }; @@ -65,15 +66,20 @@ public: t_stat() { - _rao.ReSize(3); - _rao = 0.0; - _dc = 0.0; - _dcRed = 0.0; - _offset = 0.0; - _nr = 0; - _nc = 0; + _rao.ReSize(3); _rao = 0.0; + _rao3DRMS = 0.0; + _dcRMS = 0.0; + _dcRedRMS = 0.0; + _dcRedMean = 0.0; + _dcRedSig = 0.0; + _offset = 0.0; + _nr = 0; + _nc = 0; } ColumnVector _rao; - double _dc; - double _dcRed; + double _rao3DRMS; + double _dcRMS; + double _dcRedRMS; + double _dcRedMean; + double _dcRedSig; double _offset; int _nr;