Skip to content

Commit 6caffec

Browse files
committed
Adding new script from the Contest
1 parent 64f0583 commit 6caffec

File tree

6 files changed

+102
-1
lines changed

6 files changed

+102
-1
lines changed

README.md

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -181,12 +181,13 @@ Dedicated to supplying data for [Copernicus services](http://www.esa.int/Our_Act
181181
- [Monthly composite](sentinel-2/monthly_composite)
182182
- [Water Bodies Mapping - WBM](sentinel-2/water_bodies_mapping-wbm)
183183

184-
#### Other available scripts
184+
#### Other available scripts and indices
185185
- [Selective Enhancement based on Indices](sentinel-2/selective_enhancement_based_on_indices)
186186
- [Homage to Mondrian](sentinel-2/homage_to_mondrian) - artistic script
187187
- [Index visualisation](sentinel-2/index_visualization) - universal script for visualisation of indices
188188
- [NDVI on L2A Vegetation and natural Colours](sentinel-2/ndvi-on-vegetation-natural_colours)
189189
- [PUCK](sentinel-2/puck) - Perceptually-Uniform Color Map Kit
190+
- [MixLAI](sentinel-2/mixlai) - Mix Leaf Area Index
190191

191192
#### Scripts including machine learning techniques (eo-learn)
192193

sentinel-2/mixlai/README.md

Lines changed: 93 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,93 @@
1+
# Mix Leaf Area Index (MixLAI)
2+
3+
<a href="#" id='togglescript'>Show</a> script or [download](script.js){:target="_blank"} it.
4+
<div id='script_view' style="display:none">
5+
{% highlight javascript %}
6+
{% include_relative script.js %}
7+
{% endhighlight %}
8+
</div>
9+
10+
## Evaluate and visualize
11+
- [EO Browser](https://apps.sentinel-hub.com/eo-browser/?zoom=12&lat=34.59422&lng=46.58478&themeId=DEFAULT-THEME&datasetId=S2L2A&fromTime=2020-09-20T00%3A00%3A00.000Z&toTime=2020-09-20T23%3A59%3A59.999Z&visualizationUrl=https%3A%2F%2Fservices.sentinel-hub.com%2Fogc%2Fwms%2Fbd86bcc0-f318-402b-a145-015f85b9427e&evalscript=dmFyIE1peExBSSA9ICgoQjA4LyhCMDQrQjExKSkrKEIwOC8oQjA0K0IxMikpKS8yLjA7IC8vIGNhbGN1bGF0ZSB0aGUgaW5kZXgKcmV0dXJuIFtNaXhMQUldCmlmIChNaXhMQUkgPCAwLjUpIHsKcmV0dXJuIFswLjMsIDAuMiwgMC43XSAKfQppZiAoTWl4TEFJIDwgMC43KSB7CnJldHVybiBbMC4yLCAwLjYsIDAuM10gCn0KaWYgKE1peExBSSA8IDEuNSkgewpyZXR1cm4gWzAuNSwgMC44LCAwLjJdIAp9CmVsc2UgewpyZXR1cm4gWzEsIDAuNCwgMF0gCn07Cg%3D%3D#custom-script){:target="_blank"}
12+
13+
14+
## General description of the script
15+
16+
Globally, improving remotely sensed characterization of biophysical properties of vegetation is of paramount importance for a variety of applications. Changes of Plant Leaf Area and Its Relationships with Soil Factors in the Process of Grassland Desertification. In particular, the non-destructive estimation of the Leaf Area Index (LAI) from earth observation data has been topical for decades. LAI is defined as the one-sided green leaf area per unit ground area. This script represents the leaf area index using the Sentinel 2 satellite. There are examples of this equation, but this equation uses the value of several different bands to calculate the leaf index value in a simple way.
17+
18+
To determine the validity of the index, 14 Exam (point) MCD15A3H product LAI data from around the world were used. The average correlation obtained is 0.594 which is a satisfactory and reassuring result. This data has been used in Iran and using another experimental equation that is the result of a research project under my management, the yields of crops such as dry and irrigated wheat, dry and irrigated barley, and rice have been well predicted.
19+
20+
## Details of the script
21+
22+
Leaf area index is one of the most important parameters in determining plant yields. On the other hand, agricultural scientists calculate this index using special hardware. But because of the global reach of the index in this method, the basis of validation is the data of MODIS product: MCD15A3H (500 M, 4 DAY).
23+
24+
Random data analysis has been selected from all over the world. Product data is obtained from the following address: https://lpdaacsvc.cr.usgs.gov/appeears/explore.
25+
26+
Pearson correlation and mean error methods were used to check the validity of the data. The average error is obtained with this method:
27+
28+
MAE = (ABS(MODIS - S2)/C); MAE: mean absolute error; C: count image (Picture of the same day) in the year: Sep 2019-Sep2020.
29+
30+
exam point of world:
31+
32+
long lat
33+
34+
-3.787585 40.089773
35+
36+
-4.991665 33.647475
37+
38+
75.28038 24.862526
39+
40+
24.106951 -31.276369
41+
42+
-82.253509 34.582171
43+
44+
-101.92016 22.943836
45+
46+
-73.503856 3.385921
47+
48+
149.130662 -30.375819
49+
50+
22.50151 45.943571
51+
52+
11.464104 43.512098
53+
54+
11.664784 51.26062
55+
56+
1.229418 44.669143
57+
58+
-122.048149 44.590588
59+
60+
-67.949633 -5.25761
61+
62+
## Author of the script
63+
64+
Hamid Rahimi
65+
66+
## Description of representative images
67+
68+
Exam Result
69+
70+
![Exam result](fig/exam_result.jpg)
71+
72+
Export Chart
73+
74+
![Export chart](fig/export_chart.jpg)
75+
76+
Export Image
77+
78+
![Export image](fig/export_image.jpg)
79+
80+
## References
81+
82+
[1] Gitelson, Anatoly A., ViÒa, AndrÈs, Arkebauer, Timothy J., Rundquist, Donald C., Keydan, Galina, Leavitt, Bryan, 2003. Remote Estimation of Leaf Area Index and Green Leaf Biomass in Maize Canopies.
83+
84+
[2] Lymburner, L., Beggs, P.J., Jacobson, C.R., 2000. Estimation of canopy-average surface-specific leaf area using Landsat TM data. Photogramm. Eng. Remote Sens. 66, 183ñ191.
85+
86+
[3] Mikeladze G, Gavashelishvili A, Akobia I, Metreveli V (2020). Estimation of forest cover change using Sentinel-2 multi-spectral imagery in Georgia (the Caucasus). iForest 13: 329-335. - DOI: 10.3832/ifor3386-013.
87+
88+
[4] Roya Mourad, Hadi Jaafar, Martha Anderson, Feng Gao. (2020) Assessment of Leaf Area Index Models Using Harmonized Landsat and Sentinel-2 Surface Reflectance Data over a Semi-Arid Irrigated Landscape. Remote Sensing 12:19 pages 3121.
89+
90+
## Credits
91+
92+
The present script is the result of combining and simplifying the previous methods. But the initial idea of the script was based on the following basic method:
93+
Lymburner, L., Beggs, P.J., Jacobson, C.R., 2000. Estimation of canopy-average surface-specific leaf area using Landsat TM data. Photogramm. Eng. Remote Sens. 66, 183ñ191.

sentinel-2/mixlai/fig/exam_result.jpg

80.8 KB
Loading
169 KB
Loading
163 KB
Loading

sentinel-2/mixlai/script.js

Lines changed: 7 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,7 @@
1+
// Mix Leaf Area Index with Sentinel 2 Bands//
2+
3+
//Creat by: Dr.Hamid Rahimi//
4+
5+
var MixLAI = ((B08/(B04+B11))+(B08/(B04+B12)))/2.0; // calculate the index
6+
7+
return [MixLAI];

0 commit comments

Comments
 (0)