Ground data processing & production of the high resolution biophysical maps
Date : 16/02/2018
Moncada site
Ref: EOLAB_18R06
Issue : 1.00
GROUND DATA PROCESSING & PRODUCTION OF THE HIGH RESOLUTION BIOPHYSICAL MAPS
“MONCADA SITE - SPAIN” 22ND MAY, 18TH JULY AND 8TH NOVEMBER, 2017
Contributing authors: Name Fernando Camacho Beatriz Fuster David Vinué
Organization EOLAB EOLAB EOLAB
E-mail
[email protected] [email protected] [email protected]
www.eolab.es
Ground data processing & production of the high resolution biophysical maps
Date : 16/02/2018
Moncada site
Ref: EOLAB_18R06
Issue : 1.00
TABLE OF CONTENTS Table of Contents ............................................................................................................ 2 List of Figures .................................................................................................................. 4 1.
Background of the Document ................................................................................. 8 1.1.
Executive Summary ...........................................................................................................8
1.2.
Portfolio ............................................................................................................................8
1.3.
Scope and Objectives.........................................................................................................9
1.4.
Content of the Document ..................................................................................................9
2.
Introduction ......................................................................................................... 10
3.
Study area ........................................................................................................... 12 3.1.
Location ...........................................................................................................................12
3.2.
Description of The Test Site .............................................................................................13
4.
Ground measurements ......................................................................................... 14 4.1.
Material and Methods .....................................................................................................14
4.1.1 4.1.2 4.1.3
Digital Hemispheric Photographs (DHP ..................................................................................... 14 AccuparLP80-Ceptometer.......................................................................................................... 17 LI-COR LAI-2200C plant canopy analyser ................................................................................... 18
4.2.
Spatial sampling scheme .................................................................................................20
4.3.
Ground data ....................................................................................................................21
4.3.1. 4.3.2.
5.
Data processing ......................................................................................................................... 21 Content of the Ground Dataset ................................................................................................. 25
Evaluation of the sampling ................................................................................... 31 5.1.
Principles .........................................................................................................................31
5.2.
Evaluation Based On NDVI...............................................................................................31
5.3.
Evaluation Based On Convex Hull: Product Quality Flag. .................................................32
6.
Production of ground-based maps ........................................................................ 34 6.1.
Imagery ...........................................................................................................................34
6.2.
The Transfer Function......................................................................................................34
6.2.1. 6.2.2. 6.2.3.
The regression method.............................................................................................................. 34 Band combination...................................................................................................................... 35 The selected Transfer Function ................................................................................................. 37 www.eolab.es
6.3.
Ground data processing & production of the high resolution biophysical maps
Date : 16/02/2018
Moncada site
Ref: EOLAB_18R06
Issue : 1.00
The High Resolution Ground Based Maps........................................................................41
6.3.1.
Mean Values .............................................................................................................................. 43
7.
Conclusions .......................................................................................................... 44
8.
References ........................................................................................................... 45
www.eolab.es
Ground data processing & production of the high resolution biophysical maps
Date : 16/02/2018
Moncada site
Ref: EOLAB_18R06
Issue : 1.00
LIST OF FIGURES Figure 1: Location and 3x3 km2 study area of Moncada site in Spain. 12 Figure 2: RGB color composition (red-green-blue) of Level-2A Reflectance Sentinel-2A imagery over the Moncada site 5x5 km2 (22nd May, 18th July and 8th November, 2017). Red box represents the central 3x3 km2. 13 Figure 3: Accupar LP80-Ceptometer 18 Figure 4: LAI-2200C device. 19 Figure 5: LAI-2000 optical sensor with 5 zenith angles 20 Figure 6: Location of ESUs per campaigns (May in purple, July in Green and November in blue) and the center point in yellow. Red box represents the central 3x3 km2. 21 Figure 7: Examples of Digital Hemispherical Photographs acquired in Moncada (Spain) during the field campaigns in 2017. Top: First field campaign (22nd May). Middle: Second field campaign (18th July). Bottom: Third field campaign (8th November). 22 Figure 8: Results of the CAN-EYE processing carried out on a lettuce crop ESU (ESU 36) during the first third campaign (8th November, 2017). (a) DHP images. (b) Classified images. (c) Average gap fraction and (d) Clumping factor versus view zenith angle. 22 Figure 9: Intercomparison of the calculated biophysical variable LAI (left side) and LAIeff (right side) over the ESUs with different methods: CEV5.1, CEV6.1 and Miller´s formula. Moncada (Spain) during the campaign of 22nd of May, 2017). 23 Figure 10: Intercomparison of the calculated biophysical variable LAI (left side) and LAIeff (right side) over the ESUs with different methods: CEV5.1, CEV6.1 and Miller´s formula. Moncada (Spain) during the campaign of 18th of July, 2017). 23 Figure 11: Intercomparison of the calculated biophysical variable LAI (left side) and LAIeff (right side) over the ESUs with different methods: CEV5.1, CEV6.1 and Miller´s formula. Moncada (Spain) during the campaign of 8th of November, 2017). 24 Figure 12: Inter-comparison of the measured biophysical variables over the ESUs. LAI versus FAPAR, Moncada (Spain). Right: First field campaign (22 nd May). Middle: Second field campaign (18th July). Left: Third field campaign (8th November). LAIeff in blue and LAI in cyan. 24 Figure 13: Inter-comparison of the measured biophysical variables over the ESUs. FAPAR versus FCOVER, Moncada (Spain). Right: First field campaign (22 nd May). Middle: Second field campaign (18th July). Left: Third field campaign (8th November). 25 Figure 14: LAIeff measurements at field by ESU in Moncada, 2017. Top: First field campaign (22 nd May). Middle: Second field campaign (18th July). Bottom: Third field campaign (8th November). 26 nd Figure 15: LAI measurements at field by ESU in Moncada, 2017. Top: First field campaign (22 May). Middle: Second field campaign (18th July). Bottom: Third field campaign (8th November). 27 nd Figure 16: FAPAR measurements at field by ESU in Moncada, 2017. Top: First field campaign (22 May). Middle: Second field campaign (18th July). Bottom: Third field campaign (8th November). 28 Figure 17: FCOVER measurements at field by ESU in Moncada, 2017. Top: First field campaign (22nd May). Middle: Second field campaign (18th July). Bottom: Third field campaign (8th November). 28 Figure 18: Distribution of the measured biophysical variables over the ESUs, Moncada site- Spain, during the first campaign on 22nd May, 2017. 29 Figure 19: Distribution of the measured biophysical variables over the ESUs, Moncada site- Spain, during the second campaign on 18th July, 2017. 29 www.eolab.es
Ground data processing & production of the high resolution biophysical maps
Date : 16/02/2018
Moncada site
Ref: EOLAB_18R06
Issue : 1.00
Figure 20: Distribution of the measured biophysical variables over the ESUs, Moncada site- Spain, during the first campaign on 8 th November, 2017. 30 Figure 21: Comparison of NDVI distribution between ESUs (green dots) and over the whole image (blue line). Right: First field campaign (22nd May). Middle: Second field campaign (18th July). Left: Third field campaign (8th November). 31 2 Figure 22: Convex hull test over 5x5 km : clear and dark blue correspond to the pixels belonging to the ‘strict’ and ‘large’ convex hulls. Red corresponds to the pixels for which the transfer function is extrapolating. Grey corresponds to soil mask. Left: First field campaign with 4B and NDVI (22nd May). Middle: Second field campaign (18th July). Right: Third field campaign (8th November). Black box represents the central 3x3 km2. 33 Figure 23: Test of multiple regressions (TF) applied on different band combinations. Band combinations are given in abscissa (1=G, 2=RED, 3=NIR and 4=SWIR). The weighted root mean square error (RMSE) is presented in red along with the cross-validation RMSE in green. The numbers indicate the number of data used for the robust regression with a weight lower than 0.7 that could be considered as outliers (Moncada, 22 nd of May, 2017). 35 Figure 24: Test of multiple regressions (TF) applied on different band combinations. Band combinations are given in abscissa (1=G, 2=RED, 3=NIR and 4=SWIR). The weighted root mean square error (RMSE) is presented in red along with the cross-validation RMSE in green. The numbers indicate the number of data used for the robust regression with a weight lower than 0.7 that could be considered as outliers (Moncada, 18 th July, 2017). 36 Figure 25: Test of multiple regressions (TF) applied on different band combinations. Band combinations are given in abscissa (1=G, 2=RED, 3=NIR and 4=SWIR). The weighted root mean square error (RMSE) is presented in red along with the cross-validation RMSE in green. The numbers indicate the number of data used for the robust regression with a weight lower than 0.7 that could be considered as outliers (Moncada, 8th November, 2017). 37 Figure 26: LAI, LAIeff, FAPAR and FCOVER results for regression on NDVI. Full dots: Weight>0.7. Empty dots: 0
