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

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

mercurial