Commit 3bfddf2309e7411f11e900642dcf5167b55b5bc4

Authored by Olivier Hagolle
1 parent 43908cef
Exists in master and in 1 other branch develop

add comments that were suppressed and update to the values tested by O.Hagolle

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