system_settings/output_optic_input_0.9.8h-PollyXT_Cyprus_532_210429_vf_4_GHK_0.9.8h_Py3.7.dat

Wed, 07 Jul 2021 15:42:30 +0200

author
Volker Freudenthaler
date
Wed, 07 Jul 2021 15:42:30 +0200
changeset 49
7c14c9a085ba
permissions
-rw-r--r--

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

mercurial