Commit 4ce95d7c27d54b8130918cf6d675beff4e8568ac

Authored by Olivier Hagolle
1 parent 492dc262
Exists in master and in 1 other branch develop

put back the comments and update according to CESBIO's GIPP

SENTINEL2_MUSCATE/S2A_TEST_GIP_L2COMM_L_ALLSITES_00001_20190626_21000101.EEF
... ... @@ -41,216 +41,388 @@
41 41 </Earth_Explorer_Header>
42 42 <Data_Block>
43 43 <Common_Parameters>
  44 + <!-- No data value -->
44 45 <No_Data>-10000</No_Data>
  46 + <!-- Maximum percentage of NoData pixels for the product to be considered as valid -->
45 47 <Max_No_Data_Percentage>99</Max_No_Data_Percentage>
  48 + <!-- Maximum percentage of cloudy pixels for the product to be considered as valid (in percentage) -->
46 49 <Max_Cloud_Percentage>90</Max_Cloud_Percentage>
  50 + <!-- Type of interpolator used in the LUT manipulation -->
47 51 <LUT_Interpolation>LINEAR</LUT_Interpolation>
  52 + <!-- Minimum value of the cosine of incidence angle (to avoid to normalise reflectances with negative values or close to infinite values) -->
  53 + <!-- This parameter is used in Snow Masking and Slope Correction -->
48 54 <Slope_Min_Cos_I>0.20</Slope_Min_Cos_I>
  55 + <!-- L2/L3 specifics products configuration values -->
49 56 <VAP_Configuration_Values>
  57 + <!--VAP quantification value -->
50 58 <VAP_Quantification_Value>0.050</VAP_Quantification_Value>
  59 + <!--VAP Nodata value (in [0;255]) -->
51 60 <VAP_No_Data_Value>0</VAP_No_Data_Value>
52 61 </VAP_Configuration_Values>
53 62 <AOT_Configuration_Values>
  63 + <!--AOT quantification value -->
54 64 <AOT_Quantification_Value>0.005</AOT_Quantification_Value>
  65 + <!--AOT Nodata value (in [0;255]) -->
55 66 <AOT_No_Data_Value>0</AOT_No_Data_Value>
56 67 </AOT_Configuration_Values>
57 68 <Cal_Adjust_Option>false</Cal_Adjust_Option>
58 69 <Cal_Adjust_Factor>1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0</Cal_Adjust_Factor>
  70 + <!-- Constant model to use in case of absent/invalid CAMS -->
59 71 <Constant_Model>CONTINEN</Constant_Model>
60 72 </Common_Parameters>
61 73 <Bands_Definition>
62 74 <Thematic_Definition>
  75 + <!-- Blue band code identifier -->
63 76 <Blue_Band_Code>B1</Blue_Band_Code>
  77 + <!-- Green band code identifier -->
64 78 <Green_Band_Code>B3</Green_Band_Code>
  79 + <!-- Red band code identifier -->
65 80 <Red_Band_Code>B4</Red_Band_Code>
  81 + <!-- Near Infra Red band code identifier -->
66 82 <NIR_Band_Code>B8A</NIR_Band_Code>
  83 + <!-- SWIR band code identifier -->
67 84 <SWIR_Band_Code>B11</SWIR_Band_Code>
  85 + <!--Band code identifier that corresponds to the 2.2µm spectral band -->
68 86 <Band_2200_Code>B12</Band_2200_Code>
  87 + <!-- Band code identifier used for the red channel of the quick look -->
69 88 <Quicklook_Red_Band_Code>B4</Quicklook_Red_Band_Code>
  89 + <!-- Band code identifier used for the green channel of the quick look -->
70 90 <Quicklook_Green_Band_Code>B3</Quicklook_Green_Band_Code>
  91 + <!-- Band code identifier used for the blue channel of the quick look -->
71 92 <Quicklook_Blue_Band_Code>B2</Quicklook_Blue_Band_Code>
72 93 </Thematic_Definition>
73 94 </Bands_Definition>
74 95 <Quicklook_Refl_Variation>
  96 + <!-- Minimum reflectance value in the red channel takes in care for the quick look rescaling -->
75 97 <Min_Ref_Red_Band>0.000</Min_Ref_Red_Band>
  98 + <!-- Maximum reflectance value in the red channel takes in care for the quick look rescaling -->
76 99 <Max_Ref_Red_Band>0.250</Max_Ref_Red_Band>
  100 + <!-- Minimum reflectance value in the greeb channel takes in care for the quick look rescaling -->
77 101 <Min_Ref_Green_Band>0.000</Min_Ref_Green_Band>
  102 + <!-- Maximum reflectance value in the green channel takes in care for the quick look rescaling -->
78 103 <Max_Ref_Green_Band>0.250</Max_Ref_Green_Band>
  104 + <!-- Minimum reflectance value in the blue channel takes in care for the quick look rescaling -->
79 105 <Min_Ref_Blue_Band>0.000</Min_Ref_Blue_Band>
  106 + <!-- Maximum reflectance value in the blue channel takes in care for the quick look rescaling -->
80 107 <Max_Ref_Blue_Band>0.250</Max_Ref_Blue_Band>
81 108 </Quicklook_Refl_Variation>
82 109 <Subsampling>
83   - <Saturation_Threshold>1.3</Saturation_Threshold>
84   - <Saturation_Threshold_Sub>0.900</Saturation_Threshold_Sub>
  110 + <!-- Threshold to consider a L2 resolution pixel as saturated -->
  111 + <Saturation_Threshold>1.3/Saturation_Threshold>
  112 + <!-- Threshold to consider a L2 coarse resolution image as saturated -->
  113 + <Saturation_Threshold_Sub>1.3</Saturation_Threshold_Sub>
85 114 </Subsampling>
86 115 <CAMS_Data_Preparation>
  116 + <!-- Boolean to activate CAMS -->
87 117 <Use_Cams_Data>false</Use_Cams_Data>
  118 + <!-- List of model to use in cams -->
88 119 <List_Of_Models>SULPHATE DUST SEASALT ORGANICM BLACKCAR</List_Of_Models>
  120 + <!-- RH sampling for LUT computation -->
89 121 <RH_Sampling>30 50 70 80 85 90 95</RH_Sampling>
  122 + <!-- Proportion of optical thickness -->
90 123 <limAOT>0.95</limAOT>
  124 + <!-- Number of hours defining the temporal window for CAMS Data validity -->
91 125 <Time_Windows_CAMS>24</Time_Windows_CAMS>
92 126 </CAMS_Data_Preparation>
93 127 <DTM_Processing>
  128 + <!--Spectral bands used in the DTM processing -->
94 129 <DTM_Viewing_Direction_Band_Code>B4</DTM_Viewing_Direction_Band_Code>
  130 + <!-- Euclidean distance beyong which a pixel is declared as hidden (in meter) -->
95 131 <Distance_Threshold>240.00</Distance_Threshold>
96 132 </DTM_Processing>
97 133 <Reduced_LUT_Computation>
  134 + <!--switch for the directional correction for plugin supporting it -->
98 135 <Directional_Correction_Option>true</Directional_Correction_Option>
  136 + <!-- Roughness coefficients, one by band (13 for SENTINEL2) -->
99 137 <Roughness_Coefficients>0.481 0.481 0.444 0.340 0.340 0.418 0.496 0.496 0.496 0.496 0.400 0.336 0.240</Roughness_Coefficients>
  138 + <!-- Volume coefficients, one by band (13 for SENTINEL2) -->
100 139 <Volume_Coefficients>0.102 0.102 0.136 0.134 0.134 0.121 0.107 0.107 0.107 0.107 0.120 0.132 0.145</Volume_Coefficients>
  140 + <!-- Azimuth and Zenith reference angle for which the reflectances after directional correction are expressed -->
101 141 <Azimuth_Reference>0.00</Azimuth_Reference>
102 142 <Zenith_Reference>0.00</Zenith_Reference>
103 143 </Reduced_LUT_Computation>
104 144 <Atmospheric_Absorption_Correction>
  145 + <!-- Option to use the default constant water vapour amount (if false use the GIP_L2WATV) -->
105 146 <Use_Default_Constant_Water_Amount>false</Use_Default_Constant_Water_Amount>
  147 + <!-- Water vapor band code identifier that corresponds to the 910nm or 940nm spectral band -->
106 148 <Water_Vapor_Band_Code>B9</Water_Vapor_Band_Code>
  149 + <!-- Water vapor reference band code identifier that corresponds to the 865nm spectral band -->
107 150 <Water_Vapor_Reference_Band_Code>B8A</Water_Vapor_Reference_Band_Code>
  151 + <!-- Water vapor reflectance threshold for the ref band(below is flagged as no_data) -->
108 152 <Water_Vapor_Reflectance_Threshold_Ref>0.1</Water_Vapor_Reflectance_Threshold_Ref>
  153 + <!-- Water vapor max value threshold -->
109 154 <Water_Vapor_Reflectance_Threshold>0.1</Water_Vapor_Reflectance_Threshold>
110 155 <Water_Vapor_Max_Threshold>7.0</Water_Vapor_Max_Threshold>
  156 + <!-- Initial width of the window used to fill VAP gaps : window_filling_width_init -->
111 157 <Water_Vapor_Init_Window_Radius>2</Water_Vapor_Init_Window_Radius>
  158 + <!-- Maximum width of the window used to fill VAP gaps -->
112 159 <Water_Vapor_Max_Window_Radius>16</Water_Vapor_Max_Window_Radius>
  160 + <!-- Structuring element for the binary mask dilatation -->
113 161 <Water_Vapor_SE_Valid_Radius>2</Water_Vapor_SE_Valid_Radius>
  162 + <!-- Ozone amount value at the center of the image (used if Use_Ozone_Amount set to true) unit: cm.atm = jacobson
  163 + note: 1 jacobson = 1000 dobson. 1 kg[O3]/m2 = 46698 dobson -->
114 164 <Ozone_Amount_Default_Value>0.300</Ozone_Amount_Default_Value>
  165 + <!-- Water amount value (used if Use_Water_Amount set to true) -->
115 166 <Water_Amount_Default_Value>2</Water_Amount_Default_Value>
116 167 </Atmospheric_Absorption_Correction>
117 168 <Rayleigh_Correction>
  169 + <!-- Default value of AOT for the Rayleigh correction -->
118 170 <Rayleigh_Correction_Default_AOT>0.05</Rayleigh_Correction_Default_AOT>
119 171 </Rayleigh_Correction>
120 172 <Geometric_Flags>
  173 + <!-- Maximum angular distance to sunglint under which a viewing direction is declared close to sunglint (in degree) -->
121 174 <Sun_Glint_Threshold unit="deg">30</Sun_Glint_Threshold>
  175 + <!-- Maximum angular distance to hotspot under which a viewing direction is declared close to hotspot (in degree) -->
122 176 <Hot_Spot_Threshold unit="deg">5</Hot_Spot_Threshold>
123 177 </Geometric_Flags>
124 178 <Cirrus_Mask>
  179 + <!-- Band code identifier that corresponds to the 1.38µm spectral band -->
125 180 <Cirrus_Band_Code>B10</Cirrus_Band_Code>
  181 + <!-- Minimum percentage of cirrus clouds within the image to raise the cirrus flag -->
126 182 <Min_Percent_Mask_Cirrus>100</Min_Percent_Mask_Cirrus>
  183 + <!-- Offset for reflectance threshold to detect a cirrus -->
127 184 <Cirrus_Mask_Threshold_Offset>0.007</Cirrus_Mask_Threshold_Offset>
  185 + <!-- Gain for reflectance threshold to detect a cirrus -->
128 186 <Cirrus_Mask_Threshold_Gain>0.000012</Cirrus_Mask_Threshold_Gain>
  187 + <!-- Cirrus altitude reference - Used in ReflectanceThresholdImage -->
129 188 <Cirrus_Alt_Ref>2000</Cirrus_Alt_Ref>
130 189 </Cirrus_Mask>
131 190 <Snow_Masking>
  191 + <!-- Structuring element for dilatation of the snow mask -->
132 192 <SE_Snow_Radius>1</SE_Snow_Radius>
  193 + <!-- Threshold on surface reflectance in the blue band -->
133 194 <NDSI_Threshold>0.4</NDSI_Threshold>
  195 + <!-- Threshold on the resd band surface reflectance absolute value -->
134 196 <Red_Reflectance_Threshold>0.12</Red_Reflectance_Threshold>
  197 + <!-- Threshold on the SWIR1 band surface reflectance absolute value -->
135 198 <SWIR1_Reflectance_Threshold>0.16</SWIR1_Reflectance_Threshold>
  199 + <!-- Boolean that indicates if the water mask has to be used -->
136 200 <Use_Water_Mask>true</Use_Water_Mask>
137 201 </Snow_Masking>
138 202 <Cloud_Masking>
  203 + <!-- Option to refine the cloud altitude instead of trust the stereoscopic altitude -->
139 204 <Refinement_Option>true</Refinement_Option>
  205 + <!-- Ratio of pixels above which the number of pixels in a cloud which have a good correlation is considered significant (in percentage) -->
140 206 <Discarded_Pixels_Ratio_Min>20</Discarded_Pixels_Ratio_Min>
  207 + <!-- Band used for absolute thresholding of reflectance over water zones (SWIR when available, NIR otherwise) -->
141 208 <Reflectance_Water_Band_Code>B12</Reflectance_Water_Band_Code>
  209 + <!-- Band code identifier used for the correlation to avoid false cloud detection -->
142 210 <Correl_Band_Code>B4</Correl_Band_Code>
  211 + <!-- Band code identifier used for the shadow detection -->
143 212 <Shadow_Band_Code>B4</Shadow_Band_Code>
  213 + <!-- Spectral band code identifier used as reference for viewing directions -->
144 214 <CLD_Viewing_Direction_Band_Code>B3</CLD_Viewing_Direction_Band_Code>
  215 + <!-- Threshold on surface reflectance in the blue band -->
145 216 <Cloud_Blue_Reflectance_Threshold>0.24</Cloud_Blue_Reflectance_Threshold>
  217 + <!-- Threshold on surface reflectance in the red band -->
146 218 <Cloud_Red_Reflectance_Threshold>0.15</Cloud_Red_Reflectance_Threshold>
  219 + <!-- Red Swir Ratio max threshold -->
147 220 <Nir_Red_Ratio_Max_Threshold>2</Nir_Red_Ratio_Max_Threshold>
  221 + <!-- Red Swir Ratio min threshold -->
148 222 <Nir_Red_Ratio_Min_Threshold>0.8</Nir_Red_Ratio_Min_Threshold>
  223 + <!-- Nir Swir Ratio min threshold -->
149 224 <Nir_Swir_Ratio_Min_Threshold>1</Nir_Swir_Ratio_Min_Threshold>
  225 + <!-- Threshold of correlation between images of date D and D-1 -->
150 226 <Correl_Threshold>80</Correl_Threshold>
  227 + <!-- Radius of the input masks (water and snow) dilatation -->
151 228 <SE_Input_Masks_Radius>1</SE_Input_Masks_Radius>
  229 + <!-- Threshold on surface reflectance in the Refl_water_band band for pixels declared water (if sunglint_flag is False) -->
152 230 <Water_Reflectance_Threshold>0.04</Water_Reflectance_Threshold>
  231 + <!-- Threshold on surface reflectance in the Refl_water_band band for pixels declared water (if sunglint_flag is True) -->
153 232 <Water_Sunglint_Reflectance_Threshold>0.16</Water_Sunglint_Reflectance_Threshold>
154   - <Min_Threshold_Var_Blue>0.016</Min_Threshold_Var_Blue>
  233 + <!-- Minimum Threshold on surface reflectance variation with time in the blue band -->
  234 + <Min_Threshold_Var_Blue>0.014</Min_Threshold_Var_Blue>
  235 + <!-- Maximum Threshold on surface reflectance variation with time in the blue band -->
155 236 <Max_Threshold_Var_Blue>0.060</Max_Threshold_Var_Blue>
  237 + <!-- Difference between the threshold used for the initial variation detection step and the one used for the extension step -->
156 238 <Delta_Threshold>0.003</Delta_Threshold>
  239 + <!-- Time delay after which Refl_Variation_Threshold is doubled (in days) -->
157 240 <Cloud_Forgetting_Duration>45</Cloud_Forgetting_Duration>
  241 + <!-- Default altitude attributed to cloud pixels, if stereoscopy has not detected any cloud (in meter) -->
158 242 <Default_Altitude>500</Default_Altitude>
  243 + <!-- Maximum cloud altitude, to find the zone where clouds outside the image could cast shadows (in meter) -->
159 244 <Max_Cloud_Altitude>8000</Max_Cloud_Altitude>
  245 + <!-- Maximum number of images in the stack -->
160 246 <Number_Of_Stack_Images>10</Number_Of_Stack_Images>
  247 + <!-- Parameter for comparison of relative variations in the blue and in the NIR -->
161 248 <K_NIR_Blue>1.20</K_NIR_Blue>
  249 + <!-- Parameter for comparison of variations in the blue and in the red -->
162 250 <K_Red_Blue>1.50</K_Red_Blue>
  251 + <!-- Radius of the correlation window -->
163 252 <N_Correl>3</N_Correl>
  253 + <!-- Minimum number of pixels to take correlation into account -->
164 254 <Min_Pixel_Correl>15</Min_Pixel_Correl>
  255 + <!-- Minimum altitude excursion around the altitude obtained via stereoscopy, to refine the altitude value (in meter) -->
165 256 <Delta_H_Min>0</Delta_H_Min>
  257 + <!-- Maximum altitude excursion around the altitude obtained via stereoscopy, to refine the altitude value (in meter) -->
166 258 <Delta_H_Max>8000</Delta_H_Max>
  259 + <!-- Step altitude excursion, to refine the altitude value (in meter) -->
167 260 <Delta_H_Step>250</Delta_H_Step>
  261 + <!-- Shadow/Cloud ratio used for clouds having less than shadowCloudRatioPixNb pixels
  262 + This ratio allows to calculate the maximum number of pixels for a given shadow (relatively to a cloud label) -->
168 263 <Shadow_Cloud_Ratio_Max>1.</Shadow_Cloud_Ratio_Max>
  264 + <!-- Shadow/Cloud ratio used for clouds having more pixels than shadowCloudRatioPixNb pixels
  265 + This ratio allows to calculate the maximum number of pixels for a given shadow (relatively to a cloud label) -->
169 266 <Shadow_Cloud_Ratio_Min>0.5</Shadow_Cloud_Ratio_Min>
  267 + <!-- Threshold on the number of pixels to use either shadowCloudRatioMax or shadowCloudRatioMin -->
170 268 <Shadow_Cloud_Ratio_Pix_Nb>10</Shadow_Cloud_Ratio_Pix_Nb>
  269 + <!-- Maximum value for the darkening threshold used to determine shadows zones (by projection) -->
171 270 <Proj_Shad_Darkening_Max>0.65</Proj_Shad_Darkening_Max>
  271 + <!-- Minimum percentile value used to determine the threshold for projected shadows -->
172 272 <Darkening_Min_Percentile>10</Darkening_Min_Percentile>
  273 + <!-- Multiplicative coefficient of the cloud cover percentage used for the determination of darkening threshold -->
173 274 <Cloud_Cover_Mult_Factor>1.1</Cloud_Cover_Mult_Factor>
  275 + <!-- Maximum value for the darkening threshold used to determine shadows zones (in the “radiometric” zone) -->
174 276 <Shad_Var_Darkening_Max>0.2</Shad_Var_Darkening_Max>
  277 + <!-- Percentile used to determine the darkening threshold for the shadows (located in the “radiometric zone”) -->
175 278 <Shad_Var_Percentile>30</Shad_Var_Percentile>
  279 + <!-- Minimum ratio of pixels (at the shadow label level) that have to be already declared as shadow (shadvar) in the previous cloud mask in order to discard the detection -->
176 280 <Shad_Var_Composite_Label_Ratio>0.8</Shad_Var_Composite_Label_Ratio>
  281 + <!-- Radius of the structuring element for cloud dilation -->
177 282 <SE_Cloud_Radius>1</SE_Cloud_Radius>
178 283 </Cloud_Masking>
179 284 <Water_Masking>
  285 + <!-- Number of possible water masks necessary to create the water mask of date D, ie. number of masks available in the PWA and TWA masks -->
180 286 <Water_Mask_Number>10</Water_Mask_Number>
  287 + <!-- NDVI threshold below which a pixel is probably covered by water -->
181 288 <Water_NDVI_Threshold>0.010</Water_NDVI_Threshold>
  289 + <!-- The reflectance of pixels with low enough NDVI has to be below this threshold to be regarded as possibly covered by water -->
182 290 <Water_Red_Reflectance_Threshold>0.05</Water_Red_Reflectance_Threshold>
  291 + <!-- No water can be detected for pixels for which the DEM slope is above this threshold -->
183 292 <Water_Slope_Threshold>0.15</Water_Slope_Threshold>
  293 + <!-- Reflectance threshold used if a water pixel is observed in sunglint conditions (where its reflectance will vary quickly with time) -->
184 294 <Water_Reflectance_Variation_Threshold>0.040</Water_Reflectance_Variation_Threshold>
  295 + <!-- Maximum reflectance for which the variation threshold is used -->
185 296 <Maximum_Sunglint_Reflectance>0.120</Maximum_Sunglint_Reflectance>
  297 + <!-- Minimum percent of tested times to declare a pixel as water -->
186 298 <Min_Percent_Possible_Water>50</Min_Percent_Possible_Water>
  299 + <!-- Standard deviation of the Gaussian filter applied to the SRTM Water mask (in pixel) -->
187 300 <SRTM_Sigma_Smoothing>10</SRTM_Sigma_Smoothing>
  301 + <!-- If a pixel is declared as water by SRTM, it brings a certain percentage of the amount necessary to declare a pixel a water -->
188 302 <SRTM_Percent>51</SRTM_Percent>
189 303 </Water_Masking>
190 304 <Rain_Flag>
  305 + <!-- Band code identifier used for the water masking -->
191 306 <Water_Band_Code>B12</Water_Band_Code>
  307 + <!-- Reflectance decrease threshold -->
192 308 <Reflectance_Decrease_Threshold>0.045</Reflectance_Decrease_Threshold>
193 309 </Rain_Flag>
194 310 <AOT_Estimation>
  311 + <!-- AOT Method : MULTITEMPORAL (MT), MULTISPECTRAL (MS), SPECTROTEMPORAL (ST)-->
195 312 <AOT_Method>SPECTROTEMPORAL</AOT_Method>
196 313 <AOT_KPondCAMS>0.2</AOT_KPondCAMS>
197 314 <AOT_HeightScale>2000</AOT_HeightScale>
  315 + <!-- Multi spectral method -->
198 316 <Multi_Spectral_Method>
  317 + <!-- AOT Band code identifier used to estimate MS AOT (blue and red bands) -->
199 318 <MS_AOT_Band_Code>B1 B4</MS_AOT_Band_Code>
  319 + <!-- Red Band code identifier used to estimate MS AOT -->
200 320 <MS_AOT_Red_Band_Code>B4</MS_AOT_Red_Band_Code>
  321 + <!-- First band for NDVI computing -->
201 322 <NDVI_AOT_Band1_Code>B8A</NDVI_AOT_Band1_Code>
  323 + <!-- Second band for NDVI computing -->
202 324 <NDVI_AOT_Band2_Code>B4</NDVI_AOT_Band2_Code>
  325 + <!-- Thresholds used to select pixels in the MS AOT method -->
203 326 <NDVI_Threshold>0.3</NDVI_Threshold>
204 327 <Red_TOA_Threshold>0.2</Red_TOA_Threshold>
205   - <MS_Slope>0.5</MS_Slope>
  328 + <!-- Slope of the MS linear relationship between the red and the blue reflectance. -->
  329 + <MS_Slope>0.52</MS_Slope>
  330 + <!-- Y intercept of the MS linear relationship between the red and the blue reflectance. -->
206 331 <MS_YIntercept>0.00</MS_YIntercept>
  332 + <!-- Slope for the multi spectral weight that depends on NDVI -->
207 333 <NDVI_Slope>0.7</NDVI_Slope>
  334 + <!-- Y Intercept for the multi spectral weight that depends on NDVI -->
208 335 <NDVI_YIntercept>0.3</NDVI_YIntercept>
  336 + <!-- Weight for the MS inversion equation -->
209 337 <W_MS>2.</W_MS>
210 338 </Multi_Spectral_Method>
  339 + <!-- Multi temporal method -->
211 340 <Multi_Temporal_Method>
  341 + <!-- Band code identifier used to estimate MT AOT -->
212 342 <MT_AOT_Band_Code>B1</MT_AOT_Band_Code>
  343 + <!-- Band used for the relative variation of reflectance between date D and D-1 (spec: band_var)-->
213 344 <Var_Band_Code>B8A</Var_Band_Code>
  345 + <!-- Maximum relative variation in the "band var" channel -->
214 346 <Relative_Var_Max>0.075</Relative_Var_Max>
  347 + <!-- Weight for the equation 2. This weight depends on the differences between the TOA reflectances of day D and D1 -->
215 348 <W_MT_equation2>2.</W_MT_equation2>
  349 + <!-- Channel for calculating the stdev within the neighborhood -->
216 350 <Stdev_Band_Code>B12</Stdev_Band_Code>
  351 + <!-- Minimum standard deviation of reflectances within neighborhood to perform the inversion -->
217 352 <Stdev_Min>0.02</Stdev_Min>
  353 + <!-- Usual standard deviation of TOA reflectances due to landscape change or instrumental noise -->
  354 + <!-- It is used to normalise the weight of the first term of the cost function -->
218 355 <Std_Noise>0.001</Std_Noise>
  356 + <!-- Initialization value for AOT(D-1) in vector [AOT(D-1), AOT(D)] -->
219 357 <AOT_initial_Dm1>0.1</AOT_initial_Dm1>
  358 + <!-- First AOT value used to determine TOA reflectance thresholds -->
220 359 <First_AOT>0.</First_AOT>
  360 + <!-- Second AOT value used to determine TOA reflectance thresholds -->
221 361 <Second_AOT>0.4</Second_AOT>
  362 + <!-- Minimum TOA reflectance value used to determine TOA reflectance thresholds -->
222 363 <TOA_Reflectance_Min>0.</TOA_Reflectance_Min>
  364 + <!-- Maximum TOA reflectance value used to determine TOA reflectance thresholds -->
223 365 <TOA_Reflectance_Max>1.</TOA_Reflectance_Max>
  366 + <!-- TOA reflectance step value used to determine TOA reflectance thresholds -->
224 367 <TOA_Reflectance_Step>0.01</TOA_Reflectance_Step>
  368 + <!-- Difference threshold between surface reflectances obtained for distant AOT values -->
225 369 <Min_Difference_Thresholds_Calculation>0.01</Min_Difference_Thresholds_Calculation>
226 370 </Multi_Temporal_Method>
  371 + <!-- Spectro temporal method -->
227 372 <Spectro_Temporal_Method>
  373 + <!-- Numerator coefficient used to compute the weight of the MT equations in the ST method -->
228 374 <A_Coefficient>1200</A_Coefficient>
  375 + <!-- Denominator coefficient used to compute the weight of the MT equations in the ST method -->
229 376 <B_Coefficient>800</B_Coefficient>
230   - <MT_Weight_Threshold>19</MT_Weight_Threshold>
  377 + <!-- Threshold ont he equation 1 weight coefficient in the ST method -->
  378 + <MT_Weight_Threshold>10</MT_Weight_Threshold>
231 379 </Spectro_Temporal_Method>
  380 + <!-- Maximum age for reflectances in the composite image. If the age of the pixel is greater than 60 (TBC) days, The pixel is not used to estimate AOT -->
232 381 <AOT_Max_Age>60</AOT_Max_Age>
  382 + <!-- Used AOT value if AOT cannot be estimated -->
233 383 <Default_AOT>0.10</Default_AOT>
  384 + <!-- Number of L2 coarse resolution pixel between each estimate AOT -->
  385 + <!-- AOT_resolution= L2_coarse_resolution * Step_AOT -->
234 386 <Step_AOT>3</Step_AOT>
  387 + <!-- Minimum sampling interval within a neighborhood -->
235 388 <Sampling_Interval_Min>2</Sampling_Interval_Min>
  389 + <!-- Maximum sampling interval within a neighborhood -->
236 390 <Sampling_Interval_Max>8</Sampling_Interval_Max>
237   - <N_Pix_Min>10</N_Pix_Min>
  391 + <!-- Minimum number of valid pixels within a neighborhood to perform AOT inversion -->
  392 + <N_Pix_Min>19</N_Pix_Min>
  393 + <!-- Value to increment the distance within neighbours if the standard deviation is not sufficient -->
238 394 <Increment_N_Distance>2</Increment_N_Distance>
  395 + <!-- Band code identifier used to estimate the maximum AOT with the Dark pixel method -->
239 396 <Dark_Band_Code>B1</Dark_Band_Code>
  397 + <!-- Structuring element for the opening AOT retrieved by the inversion radius -->
240 398 <SE_Valid_AOT_Radius>1</SE_Valid_AOT_Radius>
  399 + <!-- SE Water and snow masks dilatation radius -->
241 400 <SE_Mask_Dilatation_Radius>1</SE_Mask_Dilatation_Radius>
  401 + <!-- Neighborhood size : Half_Window_Size parameter in spec -->
242 402 <Half_Window_Size>3</Half_Window_Size>
  403 + <!-- Weight for the AOT maximal value equation -->
243 404 <W_dark>0.5</W_dark>
  405 + <!-- Weight for the equations to control the AOT to stay between the bounds of the LUT -->
244 406 <W_bounds>1000000000</W_bounds>
  407 + <!-- Minimum value for the AOT. AOT_min must be greater that the first index of the list of AOT -->
245 408 <AOT_min>0.01</AOT_min>
  409 + <!-- Levenberg Marquardt optimization function value tolerance -->
246 410 <F_Tolerance>0.00001</F_Tolerance>
  411 + <!-- Levenberg Marquardt optimization gradient tolerance -->
247 412 <G_Tolerance>0.00001</G_Tolerance>
  413 + <!-- Levenberg Marquardt optimization search space tolerance -->
248 414 <X_Tolerance>0.00001</X_Tolerance>
  415 + <!-- Levenberg Marquardt optimization step -->
249 416 <Epsilon_Function>0.00001</Epsilon_Function>
  417 + <!-- Levenberg Marquardt optimization maximum number of iteration -->
250 418 <Max_Iterations>500</Max_Iterations>
  419 + <!-- Initialization value for AOT(D) in vector [AOT(D-1), AOT(D)] -->
251 420 <AOT_initial_D>0.1</AOT_initial_D>
  421 + <!-- Initial width of the window used to fill AOT gaps : window_filling_width_init -->
252 422 <Init_Window_Radius>2</Init_Window_Radius>
  423 + <!-- Maximum width of the window used to fill AOT gaps -->
253 424 <Max_Window_Radius>16</Max_Window_Radius>
  425 + <!-- Std value for gap filling smoothing -->
254 426 <Sigma_Smoothing_AOT>2</Sigma_Smoothing_AOT>
255 427 </AOT_Estimation>
256 428 <Cirrus_Correction>
... ... @@ -263,22 +435,33 @@
263 435 <Cirrus_Max_Gamma_Value>1</Cirrus_Max_Gamma_Value>
264 436 <Cirrus_SWIR_Gamma_Factor>2</Cirrus_SWIR_Gamma_Factor>
265 437 <Cirrus_Threshold_Factor>1.</Cirrus_Threshold_Factor>
  438 + <!-- Cirrus maximum reflectance -->
266 439 <Cirrus_Max_Reflectance>0.07</Cirrus_Max_Reflectance>
267 440 </Cirrus_Correction>
268 441 <Composite>
  442 + <!-- Threshold on AOT. In case of a high AOT, the composite image is not updated with the new surface reflectance -->
269 443 <Max_AOT>0.500</Max_AOT>
  444 + <!-- Std max under which the std of the AOT image is not taking into account in the calculation of the high AOT threshold -->
270 445 <Minimum_High_AOT_Std>0.1</Minimum_High_AOT_Std>
  446 + <!-- Maximum age for reflectance in the composite image -->
271 447 <Max_Age>90</Max_Age>
  448 + <!-- Structuring element radius for dilatation of AOT_labels -->
272 449 <SE_AOT_Radius>1</SE_AOT_Radius>
  450 + <!-- Coefficient for the threshold to detect high AOT -->
273 451 <Coef_Sigma_AOT>1</Coef_Sigma_AOT>
  452 + <!-- Index of the digit that corresponds to the all cloud masks in the multi channel image cloud mask -->
274 453 <All_Clouds_Band>1</All_Clouds_Band>
275 454 </Composite>
276 455 <Environment_Correction>
  456 + <!-- option to apply the environment correction -->
277 457 <Env_Corr_Option>true</Env_Corr_Option>
  458 + <!-- Convolution radius -->
278 459 <Env_Corr_Radius>5</Env_Corr_Radius>
279 460 </Environment_Correction>
280 461 <Slope_Correction>
  462 + <!-- Minimum value of the cosine of exitence angle (to avoid a negative BRDF factor) -->
281 463 <Min_Cos_E>0.10</Min_Cos_E>
  464 + <!-- Minimum the cosine ratio value -->
282 465 <Min_Cos_Ratio>0.20</Min_Cos_Ratio>
283 466 </Slope_Correction>
284 467 </Data_Block>
... ...
SENTINEL2_MUSCATE/S2B_TEST_GIP_L2COMM_L_ALLSITES_00001_20190626_21000101.EEF
... ... @@ -41,216 +41,388 @@
41 41 </Earth_Explorer_Header>
42 42 <Data_Block>
43 43 <Common_Parameters>
  44 + <!-- No data value -->
44 45 <No_Data>-10000</No_Data>
  46 + <!-- Maximum percentage of NoData pixels for the product to be considered as valid -->
45 47 <Max_No_Data_Percentage>99</Max_No_Data_Percentage>
  48 + <!-- Maximum percentage of cloudy pixels for the product to be considered as valid (in percentage) -->
46 49 <Max_Cloud_Percentage>90</Max_Cloud_Percentage>
  50 + <!-- Type of interpolator used in the LUT manipulation -->
47 51 <LUT_Interpolation>LINEAR</LUT_Interpolation>
  52 + <!-- Minimum value of the cosine of incidence angle (to avoid to normalise reflectances with negative values or close to infinite values) -->
  53 + <!-- This parameter is used in Snow Masking and Slope Correction -->
48 54 <Slope_Min_Cos_I>0.20</Slope_Min_Cos_I>
  55 + <!-- L2/L3 specifics products configuration values -->
49 56 <VAP_Configuration_Values>
  57 + <!--VAP quantification value -->
50 58 <VAP_Quantification_Value>0.050</VAP_Quantification_Value>
  59 + <!--VAP Nodata value (in [0;255]) -->
51 60 <VAP_No_Data_Value>0</VAP_No_Data_Value>
52 61 </VAP_Configuration_Values>
53 62 <AOT_Configuration_Values>
  63 + <!--AOT quantification value -->
54 64 <AOT_Quantification_Value>0.005</AOT_Quantification_Value>
  65 + <!--AOT Nodata value (in [0;255]) -->
55 66 <AOT_No_Data_Value>0</AOT_No_Data_Value>
56 67 </AOT_Configuration_Values>
57 68 <Cal_Adjust_Option>false</Cal_Adjust_Option>
58 69 <Cal_Adjust_Factor>1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0</Cal_Adjust_Factor>
  70 + <!-- Constant model to use in case of absent/invalid CAMS -->
59 71 <Constant_Model>CONTINEN</Constant_Model>
60 72 </Common_Parameters>
61 73 <Bands_Definition>
62 74 <Thematic_Definition>
  75 + <!-- Blue band code identifier -->
63 76 <Blue_Band_Code>B1</Blue_Band_Code>
  77 + <!-- Green band code identifier -->
64 78 <Green_Band_Code>B3</Green_Band_Code>
  79 + <!-- Red band code identifier -->
65 80 <Red_Band_Code>B4</Red_Band_Code>
  81 + <!-- Near Infra Red band code identifier -->
66 82 <NIR_Band_Code>B8A</NIR_Band_Code>
  83 + <!-- SWIR band code identifier -->
67 84 <SWIR_Band_Code>B11</SWIR_Band_Code>
  85 + <!--Band code identifier that corresponds to the 2.2µm spectral band -->
68 86 <Band_2200_Code>B12</Band_2200_Code>
  87 + <!-- Band code identifier used for the red channel of the quick look -->
69 88 <Quicklook_Red_Band_Code>B4</Quicklook_Red_Band_Code>
  89 + <!-- Band code identifier used for the green channel of the quick look -->
70 90 <Quicklook_Green_Band_Code>B3</Quicklook_Green_Band_Code>
  91 + <!-- Band code identifier used for the blue channel of the quick look -->
71 92 <Quicklook_Blue_Band_Code>B2</Quicklook_Blue_Band_Code>
72 93 </Thematic_Definition>
73 94 </Bands_Definition>
74 95 <Quicklook_Refl_Variation>
  96 + <!-- Minimum reflectance value in the red channel takes in care for the quick look rescaling -->
75 97 <Min_Ref_Red_Band>0.000</Min_Ref_Red_Band>
  98 + <!-- Maximum reflectance value in the red channel takes in care for the quick look rescaling -->
76 99 <Max_Ref_Red_Band>0.250</Max_Ref_Red_Band>
  100 + <!-- Minimum reflectance value in the greeb channel takes in care for the quick look rescaling -->
77 101 <Min_Ref_Green_Band>0.000</Min_Ref_Green_Band>
  102 + <!-- Maximum reflectance value in the green channel takes in care for the quick look rescaling -->
78 103 <Max_Ref_Green_Band>0.250</Max_Ref_Green_Band>
  104 + <!-- Minimum reflectance value in the blue channel takes in care for the quick look rescaling -->
79 105 <Min_Ref_Blue_Band>0.000</Min_Ref_Blue_Band>
  106 + <!-- Maximum reflectance value in the blue channel takes in care for the quick look rescaling -->
80 107 <Max_Ref_Blue_Band>0.250</Max_Ref_Blue_Band>
81 108 </Quicklook_Refl_Variation>
82 109 <Subsampling>
83   - <Saturation_Threshold>1.3</Saturation_Threshold>
84   - <Saturation_Threshold_Sub>0.900</Saturation_Threshold_Sub>
  110 + <!-- Threshold to consider a L2 resolution pixel as saturated -->
  111 + <Saturation_Threshold>1.3/Saturation_Threshold>
  112 + <!-- Threshold to consider a L2 coarse resolution image as saturated -->
  113 + <Saturation_Threshold_Sub>1.3</Saturation_Threshold_Sub>
85 114 </Subsampling>
86 115 <CAMS_Data_Preparation>
  116 + <!-- Boolean to activate CAMS -->
87 117 <Use_Cams_Data>false</Use_Cams_Data>
  118 + <!-- List of model to use in cams -->
88 119 <List_Of_Models>SULPHATE DUST SEASALT ORGANICM BLACKCAR</List_Of_Models>
  120 + <!-- RH sampling for LUT computation -->
89 121 <RH_Sampling>30 50 70 80 85 90 95</RH_Sampling>
  122 + <!-- Proportion of optical thickness -->
90 123 <limAOT>0.95</limAOT>
  124 + <!-- Number of hours defining the temporal window for CAMS Data validity -->
91 125 <Time_Windows_CAMS>24</Time_Windows_CAMS>
92 126 </CAMS_Data_Preparation>
93 127 <DTM_Processing>
  128 + <!--Spectral bands used in the DTM processing -->
94 129 <DTM_Viewing_Direction_Band_Code>B4</DTM_Viewing_Direction_Band_Code>
  130 + <!-- Euclidean distance beyong which a pixel is declared as hidden (in meter) -->
95 131 <Distance_Threshold>240.00</Distance_Threshold>
96 132 </DTM_Processing>
97 133 <Reduced_LUT_Computation>
  134 + <!--switch for the directional correction for plugin supporting it -->
98 135 <Directional_Correction_Option>true</Directional_Correction_Option>
  136 + <!-- Roughness coefficients, one by band (13 for SENTINEL2) -->
99 137 <Roughness_Coefficients>0.481 0.481 0.444 0.340 0.340 0.418 0.496 0.496 0.496 0.496 0.400 0.336 0.240</Roughness_Coefficients>
  138 + <!-- Volume coefficients, one by band (13 for SENTINEL2) -->
100 139 <Volume_Coefficients>0.102 0.102 0.136 0.134 0.134 0.121 0.107 0.107 0.107 0.107 0.120 0.132 0.145</Volume_Coefficients>
  140 + <!-- Azimuth and Zenith reference angle for which the reflectances after directional correction are expressed -->
101 141 <Azimuth_Reference>0.00</Azimuth_Reference>
102 142 <Zenith_Reference>0.00</Zenith_Reference>
103 143 </Reduced_LUT_Computation>
104 144 <Atmospheric_Absorption_Correction>
  145 + <!-- Option to use the default constant water vapour amount (if false use the GIP_L2WATV) -->
105 146 <Use_Default_Constant_Water_Amount>false</Use_Default_Constant_Water_Amount>
  147 + <!-- Water vapor band code identifier that corresponds to the 910nm or 940nm spectral band -->
106 148 <Water_Vapor_Band_Code>B9</Water_Vapor_Band_Code>
  149 + <!-- Water vapor reference band code identifier that corresponds to the 865nm spectral band -->
107 150 <Water_Vapor_Reference_Band_Code>B8A</Water_Vapor_Reference_Band_Code>
  151 + <!-- Water vapor reflectance threshold for the ref band(below is flagged as no_data) -->
108 152 <Water_Vapor_Reflectance_Threshold_Ref>0.1</Water_Vapor_Reflectance_Threshold_Ref>
  153 + <!-- Water vapor max value threshold -->
109 154 <Water_Vapor_Reflectance_Threshold>0.1</Water_Vapor_Reflectance_Threshold>
110 155 <Water_Vapor_Max_Threshold>7.0</Water_Vapor_Max_Threshold>
  156 + <!-- Initial width of the window used to fill VAP gaps : window_filling_width_init -->
111 157 <Water_Vapor_Init_Window_Radius>2</Water_Vapor_Init_Window_Radius>
  158 + <!-- Maximum width of the window used to fill VAP gaps -->
112 159 <Water_Vapor_Max_Window_Radius>16</Water_Vapor_Max_Window_Radius>
  160 + <!-- Structuring element for the binary mask dilatation -->
113 161 <Water_Vapor_SE_Valid_Radius>2</Water_Vapor_SE_Valid_Radius>
  162 + <!-- Ozone amount value at the center of the image (used if Use_Ozone_Amount set to true) unit: cm.atm = jacobson
  163 + note: 1 jacobson = 1000 dobson. 1 kg[O3]/m2 = 46698 dobson -->
114 164 <Ozone_Amount_Default_Value>0.300</Ozone_Amount_Default_Value>
  165 + <!-- Water amount value (used if Use_Water_Amount set to true) -->
115 166 <Water_Amount_Default_Value>2</Water_Amount_Default_Value>
116 167 </Atmospheric_Absorption_Correction>
117 168 <Rayleigh_Correction>
  169 + <!-- Default value of AOT for the Rayleigh correction -->
118 170 <Rayleigh_Correction_Default_AOT>0.05</Rayleigh_Correction_Default_AOT>
119 171 </Rayleigh_Correction>
120 172 <Geometric_Flags>
  173 + <!-- Maximum angular distance to sunglint under which a viewing direction is declared close to sunglint (in degree) -->
121 174 <Sun_Glint_Threshold unit="deg">30</Sun_Glint_Threshold>
  175 + <!-- Maximum angular distance to hotspot under which a viewing direction is declared close to hotspot (in degree) -->
122 176 <Hot_Spot_Threshold unit="deg">5</Hot_Spot_Threshold>
123 177 </Geometric_Flags>
124 178 <Cirrus_Mask>
  179 + <!-- Band code identifier that corresponds to the 1.38µm spectral band -->
125 180 <Cirrus_Band_Code>B10</Cirrus_Band_Code>
  181 + <!-- Minimum percentage of cirrus clouds within the image to raise the cirrus flag -->
126 182 <Min_Percent_Mask_Cirrus>100</Min_Percent_Mask_Cirrus>
  183 + <!-- Offset for reflectance threshold to detect a cirrus -->
127 184 <Cirrus_Mask_Threshold_Offset>0.007</Cirrus_Mask_Threshold_Offset>
  185 + <!-- Gain for reflectance threshold to detect a cirrus -->
128 186 <Cirrus_Mask_Threshold_Gain>0.000012</Cirrus_Mask_Threshold_Gain>
  187 + <!-- Cirrus altitude reference - Used in ReflectanceThresholdImage -->
129 188 <Cirrus_Alt_Ref>2000</Cirrus_Alt_Ref>
130 189 </Cirrus_Mask>
131 190 <Snow_Masking>
  191 + <!-- Structuring element for dilatation of the snow mask -->
132 192 <SE_Snow_Radius>1</SE_Snow_Radius>
  193 + <!-- Threshold on surface reflectance in the blue band -->
133 194 <NDSI_Threshold>0.4</NDSI_Threshold>
  195 + <!-- Threshold on the resd band surface reflectance absolute value -->
134 196 <Red_Reflectance_Threshold>0.12</Red_Reflectance_Threshold>
  197 + <!-- Threshold on the SWIR1 band surface reflectance absolute value -->
135 198 <SWIR1_Reflectance_Threshold>0.16</SWIR1_Reflectance_Threshold>
  199 + <!-- Boolean that indicates if the water mask has to be used -->
136 200 <Use_Water_Mask>true</Use_Water_Mask>
137 201 </Snow_Masking>
138 202 <Cloud_Masking>
  203 + <!-- Option to refine the cloud altitude instead of trust the stereoscopic altitude -->
139 204 <Refinement_Option>true</Refinement_Option>
  205 + <!-- Ratio of pixels above which the number of pixels in a cloud which have a good correlation is considered significant (in percentage) -->
140 206 <Discarded_Pixels_Ratio_Min>20</Discarded_Pixels_Ratio_Min>
  207 + <!-- Band used for absolute thresholding of reflectance over water zones (SWIR when available, NIR otherwise) -->
141 208 <Reflectance_Water_Band_Code>B12</Reflectance_Water_Band_Code>
  209 + <!-- Band code identifier used for the correlation to avoid false cloud detection -->
142 210 <Correl_Band_Code>B4</Correl_Band_Code>
  211 + <!-- Band code identifier used for the shadow detection -->
143 212 <Shadow_Band_Code>B4</Shadow_Band_Code>
  213 + <!-- Spectral band code identifier used as reference for viewing directions -->
144 214 <CLD_Viewing_Direction_Band_Code>B3</CLD_Viewing_Direction_Band_Code>
  215 + <!-- Threshold on surface reflectance in the blue band -->
145 216 <Cloud_Blue_Reflectance_Threshold>0.24</Cloud_Blue_Reflectance_Threshold>
  217 + <!-- Threshold on surface reflectance in the red band -->
146 218 <Cloud_Red_Reflectance_Threshold>0.15</Cloud_Red_Reflectance_Threshold>
  219 + <!-- Red Swir Ratio max threshold -->
147 220 <Nir_Red_Ratio_Max_Threshold>2</Nir_Red_Ratio_Max_Threshold>
  221 + <!-- Red Swir Ratio min threshold -->
148 222 <Nir_Red_Ratio_Min_Threshold>0.8</Nir_Red_Ratio_Min_Threshold>
  223 + <!-- Nir Swir Ratio min threshold -->
149 224 <Nir_Swir_Ratio_Min_Threshold>1</Nir_Swir_Ratio_Min_Threshold>
  225 + <!-- Threshold of correlation between images of date D and D-1 -->
150 226 <Correl_Threshold>80</Correl_Threshold>
  227 + <!-- Radius of the input masks (water and snow) dilatation -->
151 228 <SE_Input_Masks_Radius>1</SE_Input_Masks_Radius>
  229 + <!-- Threshold on surface reflectance in the Refl_water_band band for pixels declared water (if sunglint_flag is False) -->
152 230 <Water_Reflectance_Threshold>0.04</Water_Reflectance_Threshold>
  231 + <!-- Threshold on surface reflectance in the Refl_water_band band for pixels declared water (if sunglint_flag is True) -->
153 232 <Water_Sunglint_Reflectance_Threshold>0.16</Water_Sunglint_Reflectance_Threshold>
154   - <Min_Threshold_Var_Blue>0.016</Min_Threshold_Var_Blue>
  233 + <!-- Minimum Threshold on surface reflectance variation with time in the blue band -->
  234 + <Min_Threshold_Var_Blue>0.014</Min_Threshold_Var_Blue>
  235 + <!-- Maximum Threshold on surface reflectance variation with time in the blue band -->
155 236 <Max_Threshold_Var_Blue>0.060</Max_Threshold_Var_Blue>
  237 + <!-- Difference between the threshold used for the initial variation detection step and the one used for the extension step -->
156 238 <Delta_Threshold>0.003</Delta_Threshold>
  239 + <!-- Time delay after which Refl_Variation_Threshold is doubled (in days) -->
157 240 <Cloud_Forgetting_Duration>45</Cloud_Forgetting_Duration>
  241 + <!-- Default altitude attributed to cloud pixels, if stereoscopy has not detected any cloud (in meter) -->
158 242 <Default_Altitude>500</Default_Altitude>
  243 + <!-- Maximum cloud altitude, to find the zone where clouds outside the image could cast shadows (in meter) -->
159 244 <Max_Cloud_Altitude>8000</Max_Cloud_Altitude>
  245 + <!-- Maximum number of images in the stack -->
160 246 <Number_Of_Stack_Images>10</Number_Of_Stack_Images>
  247 + <!-- Parameter for comparison of relative variations in the blue and in the NIR -->
161 248 <K_NIR_Blue>1.20</K_NIR_Blue>
  249 + <!-- Parameter for comparison of variations in the blue and in the red -->
162 250 <K_Red_Blue>1.50</K_Red_Blue>
  251 + <!-- Radius of the correlation window -->
163 252 <N_Correl>3</N_Correl>
  253 + <!-- Minimum number of pixels to take correlation into account -->
164 254 <Min_Pixel_Correl>15</Min_Pixel_Correl>
  255 + <!-- Minimum altitude excursion around the altitude obtained via stereoscopy, to refine the altitude value (in meter) -->
165 256 <Delta_H_Min>0</Delta_H_Min>
  257 + <!-- Maximum altitude excursion around the altitude obtained via stereoscopy, to refine the altitude value (in meter) -->
166 258 <Delta_H_Max>8000</Delta_H_Max>
  259 + <!-- Step altitude excursion, to refine the altitude value (in meter) -->
167 260 <Delta_H_Step>250</Delta_H_Step>
  261 + <!-- Shadow/Cloud ratio used for clouds having less than shadowCloudRatioPixNb pixels
  262 + This ratio allows to calculate the maximum number of pixels for a given shadow (relatively to a cloud label) -->
168 263 <Shadow_Cloud_Ratio_Max>1.</Shadow_Cloud_Ratio_Max>
  264 + <!-- Shadow/Cloud ratio used for clouds having more pixels than shadowCloudRatioPixNb pixels
  265 + This ratio allows to calculate the maximum number of pixels for a given shadow (relatively to a cloud label) -->
169 266 <Shadow_Cloud_Ratio_Min>0.5</Shadow_Cloud_Ratio_Min>
  267 + <!-- Threshold on the number of pixels to use either shadowCloudRatioMax or shadowCloudRatioMin -->
170 268 <Shadow_Cloud_Ratio_Pix_Nb>10</Shadow_Cloud_Ratio_Pix_Nb>
  269 + <!-- Maximum value for the darkening threshold used to determine shadows zones (by projection) -->
171 270 <Proj_Shad_Darkening_Max>0.65</Proj_Shad_Darkening_Max>
  271 + <!-- Minimum percentile value used to determine the threshold for projected shadows -->
172 272 <Darkening_Min_Percentile>10</Darkening_Min_Percentile>
  273 + <!-- Multiplicative coefficient of the cloud cover percentage used for the determination of darkening threshold -->
173 274 <Cloud_Cover_Mult_Factor>1.1</Cloud_Cover_Mult_Factor>
  275 + <!-- Maximum value for the darkening threshold used to determine shadows zones (in the “radiometric” zone) -->
174 276 <Shad_Var_Darkening_Max>0.2</Shad_Var_Darkening_Max>
  277 + <!-- Percentile used to determine the darkening threshold for the shadows (located in the “radiometric zone”) -->
175 278 <Shad_Var_Percentile>30</Shad_Var_Percentile>
  279 + <!-- Minimum ratio of pixels (at the shadow label level) that have to be already declared as shadow (shadvar) in the previous cloud mask in order to discard the detection -->
176 280 <Shad_Var_Composite_Label_Ratio>0.8</Shad_Var_Composite_Label_Ratio>
  281 + <!-- Radius of the structuring element for cloud dilation -->
177 282 <SE_Cloud_Radius>1</SE_Cloud_Radius>
178 283 </Cloud_Masking>
179 284 <Water_Masking>
  285 + <!-- Number of possible water masks necessary to create the water mask of date D, ie. number of masks available in the PWA and TWA masks -->
180 286 <Water_Mask_Number>10</Water_Mask_Number>
  287 + <!-- NDVI threshold below which a pixel is probably covered by water -->
181 288 <Water_NDVI_Threshold>0.010</Water_NDVI_Threshold>
  289 + <!-- The reflectance of pixels with low enough NDVI has to be below this threshold to be regarded as possibly covered by water -->
182 290 <Water_Red_Reflectance_Threshold>0.05</Water_Red_Reflectance_Threshold>
  291 + <!-- No water can be detected for pixels for which the DEM slope is above this threshold -->
183 292 <Water_Slope_Threshold>0.15</Water_Slope_Threshold>
  293 + <!-- Reflectance threshold used if a water pixel is observed in sunglint conditions (where its reflectance will vary quickly with time) -->
184 294 <Water_Reflectance_Variation_Threshold>0.040</Water_Reflectance_Variation_Threshold>
  295 + <!-- Maximum reflectance for which the variation threshold is used -->
185 296 <Maximum_Sunglint_Reflectance>0.120</Maximum_Sunglint_Reflectance>
  297 + <!-- Minimum percent of tested times to declare a pixel as water -->
186 298 <Min_Percent_Possible_Water>50</Min_Percent_Possible_Water>
  299 + <!-- Standard deviation of the Gaussian filter applied to the SRTM Water mask (in pixel) -->
187 300 <SRTM_Sigma_Smoothing>10</SRTM_Sigma_Smoothing>
  301 + <!-- If a pixel is declared as water by SRTM, it brings a certain percentage of the amount necessary to declare a pixel a water -->
188 302 <SRTM_Percent>51</SRTM_Percent>
189 303 </Water_Masking>
190 304 <Rain_Flag>
  305 + <!-- Band code identifier used for the water masking -->
191 306 <Water_Band_Code>B12</Water_Band_Code>
  307 + <!-- Reflectance decrease threshold -->
192 308 <Reflectance_Decrease_Threshold>0.045</Reflectance_Decrease_Threshold>
193 309 </Rain_Flag>
194 310 <AOT_Estimation>
  311 + <!-- AOT Method : MULTITEMPORAL (MT), MULTISPECTRAL (MS), SPECTROTEMPORAL (ST)-->
195 312 <AOT_Method>SPECTROTEMPORAL</AOT_Method>
196 313 <AOT_KPondCAMS>0.2</AOT_KPondCAMS>
197 314 <AOT_HeightScale>2000</AOT_HeightScale>
  315 + <!-- Multi spectral method -->
198 316 <Multi_Spectral_Method>
  317 + <!-- AOT Band code identifier used to estimate MS AOT (blue and red bands) -->
199 318 <MS_AOT_Band_Code>B1 B4</MS_AOT_Band_Code>
  319 + <!-- Red Band code identifier used to estimate MS AOT -->
200 320 <MS_AOT_Red_Band_Code>B4</MS_AOT_Red_Band_Code>
  321 + <!-- First band for NDVI computing -->
201 322 <NDVI_AOT_Band1_Code>B8A</NDVI_AOT_Band1_Code>
  323 + <!-- Second band for NDVI computing -->
202 324 <NDVI_AOT_Band2_Code>B4</NDVI_AOT_Band2_Code>
  325 + <!-- Thresholds used to select pixels in the MS AOT method -->
203 326 <NDVI_Threshold>0.3</NDVI_Threshold>
204 327 <Red_TOA_Threshold>0.2</Red_TOA_Threshold>
205   - <MS_Slope>0.5</MS_Slope>
  328 + <!-- Slope of the MS linear relationship between the red and the blue reflectance. -->
  329 + <MS_Slope>0.52</MS_Slope>
  330 + <!-- Y intercept of the MS linear relationship between the red and the blue reflectance. -->
206 331 <MS_YIntercept>0.00</MS_YIntercept>
  332 + <!-- Slope for the multi spectral weight that depends on NDVI -->
207 333 <NDVI_Slope>0.7</NDVI_Slope>
  334 + <!-- Y Intercept for the multi spectral weight that depends on NDVI -->
208 335 <NDVI_YIntercept>0.3</NDVI_YIntercept>
  336 + <!-- Weight for the MS inversion equation -->
209 337 <W_MS>2.</W_MS>
210 338 </Multi_Spectral_Method>
  339 + <!-- Multi temporal method -->
211 340 <Multi_Temporal_Method>
  341 + <!-- Band code identifier used to estimate MT AOT -->
212 342 <MT_AOT_Band_Code>B1</MT_AOT_Band_Code>
  343 + <!-- Band used for the relative variation of reflectance between date D and D-1 (spec: band_var)-->
213 344 <Var_Band_Code>B8A</Var_Band_Code>
  345 + <!-- Maximum relative variation in the "band var" channel -->
214 346 <Relative_Var_Max>0.075</Relative_Var_Max>
  347 + <!-- Weight for the equation 2. This weight depends on the differences between the TOA reflectances of day D and D1 -->
215 348 <W_MT_equation2>2.</W_MT_equation2>
  349 + <!-- Channel for calculating the stdev within the neighborhood -->
216 350 <Stdev_Band_Code>B12</Stdev_Band_Code>
  351 + <!-- Minimum standard deviation of reflectances within neighborhood to perform the inversion -->
217 352 <Stdev_Min>0.02</Stdev_Min>
  353 + <!-- Usual standard deviation of TOA reflectances due to landscape change or instrumental noise -->
  354 + <!-- It is used to normalise the weight of the first term of the cost function -->
218 355 <Std_Noise>0.001</Std_Noise>
  356 + <!-- Initialization value for AOT(D-1) in vector [AOT(D-1), AOT(D)] -->
219 357 <AOT_initial_Dm1>0.1</AOT_initial_Dm1>
  358 + <!-- First AOT value used to determine TOA reflectance thresholds -->
220 359 <First_AOT>0.</First_AOT>
  360 + <!-- Second AOT value used to determine TOA reflectance thresholds -->
221 361 <Second_AOT>0.4</Second_AOT>
  362 + <!-- Minimum TOA reflectance value used to determine TOA reflectance thresholds -->
222 363 <TOA_Reflectance_Min>0.</TOA_Reflectance_Min>
  364 + <!-- Maximum TOA reflectance value used to determine TOA reflectance thresholds -->
223 365 <TOA_Reflectance_Max>1.</TOA_Reflectance_Max>
  366 + <!-- TOA reflectance step value used to determine TOA reflectance thresholds -->
224 367 <TOA_Reflectance_Step>0.01</TOA_Reflectance_Step>
  368 + <!-- Difference threshold between surface reflectances obtained for distant AOT values -->
225 369 <Min_Difference_Thresholds_Calculation>0.01</Min_Difference_Thresholds_Calculation>
226 370 </Multi_Temporal_Method>
  371 + <!-- Spectro temporal method -->
227 372 <Spectro_Temporal_Method>
  373 + <!-- Numerator coefficient used to compute the weight of the MT equations in the ST method -->
228 374 <A_Coefficient>1200</A_Coefficient>
  375 + <!-- Denominator coefficient used to compute the weight of the MT equations in the ST method -->
229 376 <B_Coefficient>800</B_Coefficient>
230   - <MT_Weight_Threshold>19</MT_Weight_Threshold>
  377 + <!-- Threshold ont he equation 1 weight coefficient in the ST method -->
  378 + <MT_Weight_Threshold>10</MT_Weight_Threshold>
231 379 </Spectro_Temporal_Method>
  380 + <!-- Maximum age for reflectances in the composite image. If the age of the pixel is greater than 60 (TBC) days, The pixel is not used to estimate AOT -->
232 381 <AOT_Max_Age>60</AOT_Max_Age>
  382 + <!-- Used AOT value if AOT cannot be estimated -->
233 383 <Default_AOT>0.10</Default_AOT>
  384 + <!-- Number of L2 coarse resolution pixel between each estimate AOT -->
  385 + <!-- AOT_resolution= L2_coarse_resolution * Step_AOT -->
234 386 <Step_AOT>3</Step_AOT>
  387 + <!-- Minimum sampling interval within a neighborhood -->
235 388 <Sampling_Interval_Min>2</Sampling_Interval_Min>
  389 + <!-- Maximum sampling interval within a neighborhood -->
236 390 <Sampling_Interval_Max>8</Sampling_Interval_Max>
237   - <N_Pix_Min>10</N_Pix_Min>
  391 + <!-- Minimum number of valid pixels within a neighborhood to perform AOT inversion -->
  392 + <N_Pix_Min>19</N_Pix_Min>
  393 + <!-- Value to increment the distance within neighbours if the standard deviation is not sufficient -->
238 394 <Increment_N_Distance>2</Increment_N_Distance>
  395 + <!-- Band code identifier used to estimate the maximum AOT with the Dark pixel method -->
239 396 <Dark_Band_Code>B1</Dark_Band_Code>
  397 + <!-- Structuring element for the opening AOT retrieved by the inversion radius -->
240 398 <SE_Valid_AOT_Radius>1</SE_Valid_AOT_Radius>
  399 + <!-- SE Water and snow masks dilatation radius -->
241 400 <SE_Mask_Dilatation_Radius>1</SE_Mask_Dilatation_Radius>
  401 + <!-- Neighborhood size : Half_Window_Size parameter in spec -->
242 402 <Half_Window_Size>3</Half_Window_Size>
  403 + <!-- Weight for the AOT maximal value equation -->
243 404 <W_dark>0.5</W_dark>
  405 + <!-- Weight for the equations to control the AOT to stay between the bounds of the LUT -->
244 406 <W_bounds>1000000000</W_bounds>
  407 + <!-- Minimum value for the AOT. AOT_min must be greater that the first index of the list of AOT -->
245 408 <AOT_min>0.01</AOT_min>
  409 + <!-- Levenberg Marquardt optimization function value tolerance -->
246 410 <F_Tolerance>0.00001</F_Tolerance>
  411 + <!-- Levenberg Marquardt optimization gradient tolerance -->
247 412 <G_Tolerance>0.00001</G_Tolerance>
  413 + <!-- Levenberg Marquardt optimization search space tolerance -->
248 414 <X_Tolerance>0.00001</X_Tolerance>
  415 + <!-- Levenberg Marquardt optimization step -->
249 416 <Epsilon_Function>0.00001</Epsilon_Function>
  417 + <!-- Levenberg Marquardt optimization maximum number of iteration -->
250 418 <Max_Iterations>500</Max_Iterations>
  419 + <!-- Initialization value for AOT(D) in vector [AOT(D-1), AOT(D)] -->
251 420 <AOT_initial_D>0.1</AOT_initial_D>
  421 + <!-- Initial width of the window used to fill AOT gaps : window_filling_width_init -->
252 422 <Init_Window_Radius>2</Init_Window_Radius>
  423 + <!-- Maximum width of the window used to fill AOT gaps -->
253 424 <Max_Window_Radius>16</Max_Window_Radius>
  425 + <!-- Std value for gap filling smoothing -->
254 426 <Sigma_Smoothing_AOT>2</Sigma_Smoothing_AOT>
255 427 </AOT_Estimation>
256 428 <Cirrus_Correction>
... ... @@ -263,22 +435,33 @@
263 435 <Cirrus_Max_Gamma_Value>1</Cirrus_Max_Gamma_Value>
264 436 <Cirrus_SWIR_Gamma_Factor>2</Cirrus_SWIR_Gamma_Factor>
265 437 <Cirrus_Threshold_Factor>1.</Cirrus_Threshold_Factor>
  438 + <!-- Cirrus maximum reflectance -->
266 439 <Cirrus_Max_Reflectance>0.07</Cirrus_Max_Reflectance>
267 440 </Cirrus_Correction>
268 441 <Composite>
  442 + <!-- Threshold on AOT. In case of a high AOT, the composite image is not updated with the new surface reflectance -->
269 443 <Max_AOT>0.500</Max_AOT>
  444 + <!-- Std max under which the std of the AOT image is not taking into account in the calculation of the high AOT threshold -->
270 445 <Minimum_High_AOT_Std>0.1</Minimum_High_AOT_Std>
  446 + <!-- Maximum age for reflectance in the composite image -->
271 447 <Max_Age>90</Max_Age>
  448 + <!-- Structuring element radius for dilatation of AOT_labels -->
272 449 <SE_AOT_Radius>1</SE_AOT_Radius>
  450 + <!-- Coefficient for the threshold to detect high AOT -->
273 451 <Coef_Sigma_AOT>1</Coef_Sigma_AOT>
  452 + <!-- Index of the digit that corresponds to the all cloud masks in the multi channel image cloud mask -->
274 453 <All_Clouds_Band>1</All_Clouds_Band>
275 454 </Composite>
276 455 <Environment_Correction>
  456 + <!-- option to apply the environment correction -->
277 457 <Env_Corr_Option>true</Env_Corr_Option>
  458 + <!-- Convolution radius -->
278 459 <Env_Corr_Radius>5</Env_Corr_Radius>
279 460 </Environment_Correction>
280 461 <Slope_Correction>
  462 + <!-- Minimum value of the cosine of exitence angle (to avoid a negative BRDF factor) -->
281 463 <Min_Cos_E>0.10</Min_Cos_E>
  464 + <!-- Minimum the cosine ratio value -->
282 465 <Min_Cos_Ratio>0.20</Min_Cos_Ratio>
283 466 </Slope_Correction>
284 467 </Data_Block>
... ...
SENTINEL2_MUSCATE_CAMS/S2A_TEST_GIP_L2COMM_L_ALLSITES_00001_20190626_21000101.EEF
... ... @@ -41,216 +41,388 @@
41 41 </Earth_Explorer_Header>
42 42 <Data_Block>
43 43 <Common_Parameters>
  44 + <!-- No data value -->
44 45 <No_Data>-10000</No_Data>
  46 + <!-- Maximum percentage of NoData pixels for the product to be considered as valid -->
45 47 <Max_No_Data_Percentage>99</Max_No_Data_Percentage>
  48 + <!-- Maximum percentage of cloudy pixels for the product to be considered as valid (in percentage) -->
46 49 <Max_Cloud_Percentage>90</Max_Cloud_Percentage>
  50 + <!-- Type of interpolator used in the LUT manipulation -->
47 51 <LUT_Interpolation>LINEAR</LUT_Interpolation>
  52 + <!-- Minimum value of the cosine of incidence angle (to avoid to normalise reflectances with negative values or close to infinite values) -->
  53 + <!-- This parameter is used in Snow Masking and Slope Correction -->
48 54 <Slope_Min_Cos_I>0.20</Slope_Min_Cos_I>
  55 + <!-- L2/L3 specifics products configuration values -->
49 56 <VAP_Configuration_Values>
  57 + <!--VAP quantification value -->
50 58 <VAP_Quantification_Value>0.050</VAP_Quantification_Value>
  59 + <!--VAP Nodata value (in [0;255]) -->
51 60 <VAP_No_Data_Value>0</VAP_No_Data_Value>
52 61 </VAP_Configuration_Values>
53 62 <AOT_Configuration_Values>
  63 + <!--AOT quantification value -->
54 64 <AOT_Quantification_Value>0.005</AOT_Quantification_Value>
  65 + <!--AOT Nodata value (in [0;255]) -->
55 66 <AOT_No_Data_Value>0</AOT_No_Data_Value>
56 67 </AOT_Configuration_Values>
57 68 <Cal_Adjust_Option>false</Cal_Adjust_Option>
58 69 <Cal_Adjust_Factor>1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0</Cal_Adjust_Factor>
  70 + <!-- Constant model to use in case of absent/invalid CAMS -->
59 71 <Constant_Model>CONTINEN</Constant_Model>
60 72 </Common_Parameters>
61 73 <Bands_Definition>
62 74 <Thematic_Definition>
  75 + <!-- Blue band code identifier -->
63 76 <Blue_Band_Code>B1</Blue_Band_Code>
  77 + <!-- Green band code identifier -->
64 78 <Green_Band_Code>B3</Green_Band_Code>
  79 + <!-- Red band code identifier -->
65 80 <Red_Band_Code>B4</Red_Band_Code>
  81 + <!-- Near Infra Red band code identifier -->
66 82 <NIR_Band_Code>B8A</NIR_Band_Code>
  83 + <!-- SWIR band code identifier -->
67 84 <SWIR_Band_Code>B11</SWIR_Band_Code>
  85 + <!--Band code identifier that corresponds to the 2.2µm spectral band -->
68 86 <Band_2200_Code>B12</Band_2200_Code>
  87 + <!-- Band code identifier used for the red channel of the quick look -->
69 88 <Quicklook_Red_Band_Code>B4</Quicklook_Red_Band_Code>
  89 + <!-- Band code identifier used for the green channel of the quick look -->
70 90 <Quicklook_Green_Band_Code>B3</Quicklook_Green_Band_Code>
  91 + <!-- Band code identifier used for the blue channel of the quick look -->
71 92 <Quicklook_Blue_Band_Code>B2</Quicklook_Blue_Band_Code>
72 93 </Thematic_Definition>
73 94 </Bands_Definition>
74 95 <Quicklook_Refl_Variation>
  96 + <!-- Minimum reflectance value in the red channel takes in care for the quick look rescaling -->
75 97 <Min_Ref_Red_Band>0.000</Min_Ref_Red_Band>
  98 + <!-- Maximum reflectance value in the red channel takes in care for the quick look rescaling -->
76 99 <Max_Ref_Red_Band>0.250</Max_Ref_Red_Band>
  100 + <!-- Minimum reflectance value in the greeb channel takes in care for the quick look rescaling -->
77 101 <Min_Ref_Green_Band>0.000</Min_Ref_Green_Band>
  102 + <!-- Maximum reflectance value in the green channel takes in care for the quick look rescaling -->
78 103 <Max_Ref_Green_Band>0.250</Max_Ref_Green_Band>
  104 + <!-- Minimum reflectance value in the blue channel takes in care for the quick look rescaling -->
79 105 <Min_Ref_Blue_Band>0.000</Min_Ref_Blue_Band>
  106 + <!-- Maximum reflectance value in the blue channel takes in care for the quick look rescaling -->
80 107 <Max_Ref_Blue_Band>0.250</Max_Ref_Blue_Band>
81 108 </Quicklook_Refl_Variation>
82 109 <Subsampling>
83   - <Saturation_Threshold>1.3</Saturation_Threshold>
84   - <Saturation_Threshold_Sub>0.900</Saturation_Threshold_Sub>
  110 + <!-- Threshold to consider a L2 resolution pixel as saturated -->
  111 + <Saturation_Threshold>1.3/Saturation_Threshold>
  112 + <!-- Threshold to consider a L2 coarse resolution image as saturated -->
  113 + <Saturation_Threshold_Sub>1.3</Saturation_Threshold_Sub>
85 114 </Subsampling>
86 115 <CAMS_Data_Preparation>
  116 + <!-- Boolean to activate CAMS -->
87 117 <Use_Cams_Data>true</Use_Cams_Data>
  118 + <!-- List of model to use in cams -->
88 119 <List_Of_Models>SULPHATE DUST SEASALT ORGANICM BLACKCAR</List_Of_Models>
  120 + <!-- RH sampling for LUT computation -->
89 121 <RH_Sampling>30 50 70 80 85 90 95</RH_Sampling>
  122 + <!-- Proportion of optical thickness -->
90 123 <limAOT>0.95</limAOT>
  124 + <!-- Number of hours defining the temporal window for CAMS Data validity -->
91 125 <Time_Windows_CAMS>24</Time_Windows_CAMS>
92 126 </CAMS_Data_Preparation>
93 127 <DTM_Processing>
  128 + <!--Spectral bands used in the DTM processing -->
94 129 <DTM_Viewing_Direction_Band_Code>B4</DTM_Viewing_Direction_Band_Code>
  130 + <!-- Euclidean distance beyong which a pixel is declared as hidden (in meter) -->
95 131 <Distance_Threshold>240.00</Distance_Threshold>
96 132 </DTM_Processing>
97 133 <Reduced_LUT_Computation>
  134 + <!--switch for the directional correction for plugin supporting it -->
98 135 <Directional_Correction_Option>true</Directional_Correction_Option>
  136 + <!-- Roughness coefficients, one by band (13 for SENTINEL2) -->
99 137 <Roughness_Coefficients>0.481 0.481 0.444 0.340 0.340 0.418 0.496 0.496 0.496 0.496 0.400 0.336 0.240</Roughness_Coefficients>
  138 + <!-- Volume coefficients, one by band (13 for SENTINEL2) -->
100 139 <Volume_Coefficients>0.102 0.102 0.136 0.134 0.134 0.121 0.107 0.107 0.107 0.107 0.120 0.132 0.145</Volume_Coefficients>
  140 + <!-- Azimuth and Zenith reference angle for which the reflectances after directional correction are expressed -->
101 141 <Azimuth_Reference>0.00</Azimuth_Reference>
102 142 <Zenith_Reference>0.00</Zenith_Reference>
103 143 </Reduced_LUT_Computation>
104 144 <Atmospheric_Absorption_Correction>
  145 + <!-- Option to use the default constant water vapour amount (if false use the GIP_L2WATV) -->
105 146 <Use_Default_Constant_Water_Amount>false</Use_Default_Constant_Water_Amount>
  147 + <!-- Water vapor band code identifier that corresponds to the 910nm or 940nm spectral band -->
106 148 <Water_Vapor_Band_Code>B9</Water_Vapor_Band_Code>
  149 + <!-- Water vapor reference band code identifier that corresponds to the 865nm spectral band -->
107 150 <Water_Vapor_Reference_Band_Code>B8A</Water_Vapor_Reference_Band_Code>
  151 + <!-- Water vapor reflectance threshold for the ref band(below is flagged as no_data) -->
108 152 <Water_Vapor_Refl