Wed, 07 Jul 2021 15:42:30 +0200
2nd order polynomial fits to
- K(LDRCal) and
- LDR-uncertainty max. and min depending on LDRtrue
Volker@49 | 1 | From C:\Projects\Earlinet-Asos\Lidarsystems\calc_lidar_correction_parameters |
Volker@49 | 2 | Running GHK_0.9.8h_Py3.7.py |
Volker@49 | 3 | Reading input file optic_input_0.9.8h-PollyXT_Cyprus_532_210429_vf_4.py |
Volker@49 | 4 | for Lidar system: li , PollyXT Cyprus 210429 532 approx 4 |
Volker@49 | 5 | --- Input parameters: value ±error / ±steps ---------------------- |
Volker@49 | 6 | Laser: Qin = 0.9672± 0.0100/ 1 |
Volker@49 | 7 | Vin = 0.0000± 0.0000/ 0 |
Volker@49 | 8 | Rotation alpha = 91.6500± 0.2400/ 1 |
Volker@49 | 9 | => DOP 0.9672 (degree of polarisation) |
Volker@49 | 10 | Optic: Diatt., Tunpol, Retard., Rotation (deg) |
Volker@49 | 11 | Emitter 0.0000 ± 0.0200 / 0, 1.0000, 0±180/ 0, 0.0000± 1.0000/ 0 |
Volker@49 | 12 | Receiver 0.0000 ± 0.0200 / 0, 1.0000, 0±180/ 0, 0.0000± 0.5000/ 0 |
Volker@49 | 13 | Calibrator 0.999280± 0.000210/ 1, 0.4000, 0± 0/ 0, 0.0000± 0.1000/ 0 |
Volker@49 | 14 | Pol.-filter ------ |
Volker@49 | 15 | ERT, RotT : 1.0000± 0.0000/ 0, 0.0000± 1.0000/ 1 |
Volker@49 | 16 | ERR, RotR : 1.0000± 0.0000/ 0, 0.0000± 1.0000/ 0 |
Volker@49 | 17 | PBS ------ |
Volker@49 | 18 | TP,TS : 0.5000± 0.0100/ 1, 0.5000± 0.0100/ 1 |
Volker@49 | 19 | RP,RS : 0.0020± 0.0010/ 1, 0.9980± 0.0010/ 1 |
Volker@49 | 20 | DT,TT,DR,TR,Y : 0.0000, 0.5000, -0.9960, 0.5000, -1 |
Volker@49 | 21 | Combined PBS + Pol.-filter ------ |
Volker@49 | 22 | DT,TT,DR,TR : 0.0000, 0.5000, -0.9960, 0.5000 |
Volker@49 | 23 | LDRCal during calibration in calibration range: 0.110± 0.100/ 1 |
Volker@49 | 24 | --- Additional ND filter attenuation (transmission) during the calibration --- |
Volker@49 | 25 | TCalT,TCalR : 1.0000± 0.0100/ 0, 1.0000± 0.0001/ 0 |
Volker@49 | 26 | |
Volker@49 | 27 | Linear polarizer calibrator is located before the receiver |
Volker@49 | 28 | Rotation error epsilon is considered also for normal measurements = False |
Volker@49 | 29 | The PBS incidence plane is perpendicular to the reference plane |
Volker@49 | 30 | The laser polarisation in the reference plane is finally reflected => the parallel laser polarisation is detected in the transmitted channel. |
Volker@49 | 31 | RS_RP_depend_on_TS_TP = False |
Volker@49 | 32 | |
Volker@49 | 33 | =========================================================================================================== |
Volker@49 | 34 | GR , GT , HR , HT , K(0.110), K(0.004), K(0.05) , K(0.1) , K(0.2) , K(0.3) , K(0.45) |
Volker@49 | 35 | 1.00000, 1.00000,-0.96173, 0.00000, 0.97068 , 0.96369, 0.96690, 0.97008, 0.97564, 0.98033, 0.98615 |
Volker@49 | 36 | =========================================================================================================== |
Volker@49 | 37 | |
Volker@49 | 38 | 2nd order polyfit: K(LDRcal) = a + (b * LDRcal) + (c * LDRcal^2) |
Volker@49 | 39 | (LDRcal is the volume linear depolarisation ratio in the range used for the polarisation calibration) |
Volker@49 | 40 | a = 0.96352 |
Volker@49 | 41 | b = 0.06904 |
Volker@49 | 42 | c = -0.04188 |
Volker@49 | 43 | |
Volker@49 | 44 | Errors from neglecting GHK corrections and/or calibration: |
Volker@49 | 45 | LDRtrue, LDRunCorr,1/LDRunCorr, LDRsimx, 1/LDRsimx, LDRCorr |
Volker@49 | 46 | 0.00400, 0.04732, 21.13427, 0.04593, 21.77255, 0.00400 |
Volker@49 | 47 | 0.02000, 0.07828, 12.77531, 0.07598, 13.16114, 0.02000 |
Volker@49 | 48 | 0.10000, 0.21956, 4.55449, 0.21313, 4.69204, 0.10000 |
Volker@49 | 49 | 0.30000, 0.49670, 2.01327, 0.48214, 2.07407, 0.30000 |
Volker@49 | 50 | 0.45000, 0.65439, 1.52814, 0.63520, 1.57430, 0.45000 |
Volker@49 | 51 | LDRsimx = LDR of the nominal system directly from measured signals without calibration and GHK-corrections |
Volker@49 | 52 | LDRunCorr = LDR of the nominal system directly from measured signals with calibration but without GHK-corrections; electronic amplifications = 1 assumed |
Volker@49 | 53 | LDRCorr = LDR calibrated and GHK-corrected |
Volker@49 | 54 | |
Volker@49 | 55 | Errors from signal noise: |
Volker@49 | 56 | Signal counts: NI, NCalT, NCalR, NILfac, nNCal, nNI, stdev(NI)/NI = 100000000, 100000000, 100000000, 1, 0, 0, 0.00010 |
Volker@49 | 57 | |
Volker@49 | 58 | |
Volker@49 | 59 | Lidar ID: PollyXT Cyprus 210429 532 approx 4 |
Volker@49 | 60 | |
Volker@49 | 61 | Correction K(LDRCal in calibration range) |
Volker@49 | 62 | LDRCal, K(LDRCal) |
Volker@49 | 63 | 0.004, 0.96369 |
Volker@49 | 64 | 0.050, 0.96690 |
Volker@49 | 65 | 0.100, 0.97008 |
Volker@49 | 66 | 0.200, 0.97564 |
Volker@49 | 67 | 0.300, 0.98033 |
Volker@49 | 68 | 0.450, 0.98615 |
Volker@49 | 69 | |
Volker@49 | 70 | minimum and maximum values of the distributions of possibly measured LDR for different LDRtrue |
Volker@49 | 71 | LDRtrue , LDRmin, LDRmax |
Volker@49 | 72 | 0.0040,-0.0029, 0.0114 |
Volker@49 | 73 | 0.0200, 0.0127, 0.0279 |
Volker@49 | 74 | 0.1000, 0.0906, 0.1102 |
Volker@49 | 75 | 0.3000, 0.2863, 0.3151 |
Volker@49 | 76 | 0.4500, 0.4340, 0.4679 |
Volker@49 | 77 | |
Volker@49 | 78 | LDRCal is the volume linear depolarisation ratio in the range used for the polarisation calibration. K(LDRCal) is the variability of K depending on LDRCal. |
Volker@49 | 79 | If LDRCal is not known, its uncertainty has to be included in the input file for the error calcualtion below. |
Volker@49 | 80 | LDR_min-true(LDR_true) and LDR_max-true(LDR_true) are the difference between the lowest and highest possible retrieved LDR, respectively, and the true LDR depending on the true LDR. They represent the negative and positive uncertainty of the retrieved LDR resulting from all known uncertainties of the optical input parameters. |
Volker@49 | 81 | ----------------------------------------------- |
Volker@49 | 82 | 2nd order polyfit = a + (b * LDR) + (c * LDR^2) |
Volker@49 | 83 | fit-parameter, K(LDRCal), LDR_min-true(LDRtrue), LDR_max-true(LDRtrue) |
Volker@49 | 84 | a, +0.96352, -0.00674, +0.00727 |
Volker@49 | 85 | b, +0.06904, -0.02805, +0.03080 |
Volker@49 | 86 | c, -0.04188, +0.01659, -0.01587 |
Volker@49 | 87 | |
Volker@49 | 88 | ################################################################################## |
Volker@49 | 89 | Further statistical analyses of the uncertainties, their distribution and sources: |
Volker@49 | 90 | |
Volker@49 | 91 | LDRtrue , mean , median, max-mean, min-mean, std, excess_kurtosis, skewness |
Volker@49 | 92 | 0.00400, 0.00404, 0.00402, 0.00739,-0.00685, 0.00431, -1.37505, 0.01060 |
Volker@49 | 93 | 0.02000, 0.02004, 0.02003, 0.00788,-0.00729, 0.00432, -1.35168, 0.01361 |
Volker@49 | 94 | 0.10000, 0.10006, 0.10001, 0.01019,-0.00938, 0.00447, -1.11574, 0.02519 |
Volker@49 | 95 | 0.30000, 0.30013, 0.30010, 0.01508,-0.01366, 0.00523, -0.47128, 0.06082 |
Volker@49 | 96 | 0.45000, 0.45018, 0.45000, 0.01792,-0.01600, 0.00607, -0.34137, 0.12718 |
Volker@49 | 97 | |
Volker@49 | 98 | minimum and maximum values of the distributions of signal ratios and calibration factors for different LDRtrue |
Volker@49 | 99 | LDRtrue , LDRsim, (max-min)/2, relerr |
Volker@49 | 100 | 0.00400, 0.04596, 0.01325, 0.28826 |
Volker@49 | 101 | 0.02000, 0.07602, 0.01353, 0.17791 |
Volker@49 | 102 | 0.10000, 0.21319, 0.01478, 0.06935 |
Volker@49 | 103 | 0.30000, 0.48224, 0.01725, 0.03578 |
Volker@49 | 104 | 0.45000, 0.63531, 0.01866, 0.02937 |
Volker@49 | 105 | LDRtrue , Etax , (max-min)/2, relerr |
Volker@49 | 106 | 0.00400, 0.97066, 0.03207, 0.03304 |
Volker@49 | 107 | 0.02000, 0.97066, 0.03207, 0.03304 |
Volker@49 | 108 | 0.10000, 0.97066, 0.03207, 0.03304 |
Volker@49 | 109 | 0.30000, 0.97066, 0.03207, 0.03304 |
Volker@49 | 110 | 0.45000, 0.97066, 0.03207, 0.03304 |
Volker@49 | 111 | LDRtrue , Etapx , (max-min)/2, relerr |
Volker@49 | 112 | 0.00400, 0.97203, 0.03152, 0.03243 |
Volker@49 | 113 | 0.02000, 0.97203, 0.03152, 0.03243 |
Volker@49 | 114 | 0.10000, 0.97203, 0.03152, 0.03243 |
Volker@49 | 115 | 0.30000, 0.97203, 0.03152, 0.03243 |
Volker@49 | 116 | 0.45000, 0.97203, 0.03152, 0.03243 |
Volker@49 | 117 | LDRtrue , Etamx , (max-min)/2, relerr |
Volker@49 | 118 | 0.00400, 0.96930, 0.03303, 0.03408 |
Volker@49 | 119 | 0.02000, 0.96930, 0.03303, 0.03408 |
Volker@49 | 120 | 0.10000, 0.96930, 0.03303, 0.03408 |
Volker@49 | 121 | 0.30000, 0.96930, 0.03303, 0.03408 |
Volker@49 | 122 | 0.45000, 0.96930, 0.03303, 0.03408 |
Volker@49 | 123 | |
Volker@49 | 124 | LDRCorr: contribution to the total error (per mill) |
Volker@49 | 125 | of individual parameter errors of combined parameter errors |
Volker@49 | 126 | at LDRtrue 0.004 0.020 0.100 0.300 0.450 0.004 0.020 0.100 0.300 0.450 |
Volker@49 | 127 | |
Volker@49 | 128 | Qin 725 680 519 318 230 21 19 13 4 0 |
Volker@49 | 129 | RotL 34 32 25 15 12 1 1 1 1 1 |
Volker@49 | 130 | DiC 15 24 62 133 185 7 8 15 28 37 |
Volker@49 | 131 | TP 31 49 104 145 138 9 9 7 3 0 |
Volker@49 | 132 | TS 31 49 104 145 138 9 9 7 3 0 |
Volker@49 | 133 | RP 141 132 101 62 45 4 4 3 1 0 |
Volker@49 | 134 | RS 0 0 0 0 0 0 0 0 0 0 |
Volker@49 | 135 | RotaT 0 0 0 0 0 0 0 0 0 0 |
Volker@49 | 136 | LDRCal 21 33 84 182 252 10 12 19 34 44 |
Volker@49 | 137 | Sum 999 999 999 999 1000 61 62 64 74 83 |
Volker@49 | 138 | |
Volker@49 | 139 | LDRsim: contribution to the total error (per mill) |
Volker@49 | 140 | of individual parameter errors of combined parameter errors |
Volker@49 | 141 | at LDRtrue 0.004 0.020 0.100 0.300 0.450 0.004 0.020 0.100 0.300 0.450 |
Volker@49 | 142 | |
Volker@49 | 143 | Qin 745 706 550 310 202 15 14 11 6 4 |
Volker@49 | 144 | RotL 35 33 26 14 9 0 0 0 0 0 |
Volker@49 | 145 | DiC 0 0 0 0 0 0 0 0 0 0 |
Volker@49 | 146 | TP 68 108 258 425 465 18 17 13 9 9 |
Volker@49 | 147 | TS 2 4 30 134 216 1 1 4 8 9 |
Volker@49 | 148 | RP 148 143 121 88 73 2 2 2 1 1 |
Volker@49 | 149 | RS 3 5 14 28 34 1 1 1 1 1 |
Volker@49 | 150 | RotaT 0 0 0 0 0 0 0 0 0 0 |
Volker@49 | 151 | LDRCal 0 0 0 0 0 0 0 0 0 0 |
Volker@49 | 152 | Sum 1000 1000 1000 1000 1000 37 35 30 26 25 |
Volker@49 | 153 | |
Volker@49 | 154 | EtaX, D90: contribution to the total error (per mill) |
Volker@49 | 155 | of individual parameter errors of combined parameter errors |
Volker@49 | 156 | at LDRtrue 0.004 0.020 0.100 0.300 0.450 0.004 0.020 0.100 0.300 0.450 |
Volker@49 | 157 | |
Volker@49 | 158 | Qin 10 10 10 10 10 3 3 3 3 3 |
Volker@49 | 159 | RotL 0 0 0 0 0 0 0 0 0 0 |
Volker@49 | 160 | DiC 136 136 136 136 136 24 24 24 24 24 |
Volker@49 | 161 | TP 312 312 312 312 312 7 7 7 7 7 |
Volker@49 | 162 | TS 294 294 294 294 294 13 13 13 13 13 |
Volker@49 | 163 | RP 32 32 32 32 32 0 0 0 0 0 |
Volker@49 | 164 | RS 30 30 30 30 30 1 1 1 1 1 |
Volker@49 | 165 | RotaT 0 0 0 0 0 0 0 0 0 0 |
Volker@49 | 166 | LDRCal 186 186 186 186 186 29 29 29 29 29 |
Volker@49 | 167 | Sum 1000 1000 1000 1000 1000 78 78 78 78 78 |
Volker@49 | 168 | |
Volker@49 | 169 | Etapx, +45°: contribution to the total error (per mill) |
Volker@49 | 170 | of individual parameter errors of combined parameter errors |
Volker@49 | 171 | at LDRtrue 0.004 0.020 0.100 0.300 0.450 0.004 0.020 0.100 0.300 0.450 |
Volker@49 | 172 | |
Volker@49 | 173 | Qin 9 9 9 9 9 3 3 3 3 3 |
Volker@49 | 174 | RotL 6 6 6 6 6 3 3 3 3 3 |
Volker@49 | 175 | DiC 132 132 132 132 132 23 23 23 23 23 |
Volker@49 | 176 | TP 317 317 317 317 317 7 7 7 7 7 |
Volker@49 | 177 | TS 300 300 300 300 300 13 13 13 13 13 |
Volker@49 | 178 | RP 33 33 33 33 33 0 0 0 0 0 |
Volker@49 | 179 | RS 31 31 31 31 31 1 1 1 1 1 |
Volker@49 | 180 | RotaT 0 0 0 0 0 0 0 0 0 0 |
Volker@49 | 181 | LDRCal 172 172 172 172 172 28 28 28 28 28 |
Volker@49 | 182 | Sum 1000 1000 1000 1000 1000 79 79 79 79 79 |
Volker@49 | 183 | |
Volker@49 | 184 | Etamx -45°: contribution to the total error (per mill) |
Volker@49 | 185 | of individual parameter errors of combined parameter errors |
Volker@49 | 186 | at LDRtrue 0.004 0.020 0.100 0.300 0.450 0.004 0.020 0.100 0.300 0.450 |
Volker@49 | 187 | |
Volker@49 | 188 | Qin 10 10 10 10 10 4 4 4 4 4 |
Volker@49 | 189 | RotL 6 6 6 6 6 4 4 4 4 4 |
Volker@49 | 190 | DiC 139 139 139 139 139 27 27 27 27 27 |
Volker@49 | 191 | TP 303 303 303 303 303 7 7 7 7 7 |
Volker@49 | 192 | TS 284 284 284 284 284 13 13 13 13 13 |
Volker@49 | 193 | RP 31 31 31 31 31 0 0 0 0 0 |
Volker@49 | 194 | RS 29 29 29 29 29 1 1 1 1 1 |
Volker@49 | 195 | RotaT 0 0 0 0 0 0 0 0 0 0 |
Volker@49 | 196 | LDRCal 197 197 197 197 197 33 33 33 33 33 |
Volker@49 | 197 | Sum 999 999 999 999 999 89 89 89 89 89 |