Calculation of long term average crop indicators
From Agri4castWiki
The calculation of long term averages is done for each crop indicator and for each individual geographical feature on a 10 daily base.
For the following basic parameters average values are calculated:
- Potential above ground biomass
- Water limited above ground biomass
- Potential storage organs biomass
- Water limited storage organs biomass
- Potential Leaf Area Index
- Water limited Leaf Area Index
- Total water requirement
- Total water consumption
- Relative soil moisture
- Crop development stage
Besides the average also minimum, maximum, standard deviation and count are calculated (grid and nuts resolution only).
The data are stored in separate tables:
Two versions of each table exist; one for yields based on observed weather and one for yields based on ECMWF forcasted weather (HIS model).
| pseudo code |
|---|
select
crop_no,
ddmm, -- day converted to mmdd format
grid_no,
smu_no,
average(indicator_1),
average(indicator_2),
...
from CROP_YIELD
where year between start year and end year -- day converted to yyyy format
group by crop_no, ddmm, grid_no, smu_no -- day converted to mmdd format
select
crop_no,
mmdd, -- day converted to mmdd format
grid_no,
average(indicator_1),
average(indicator_2),
...
min(indicator_1),
min(indicator_2),
...
max(indicator_1),
max(indicator_2),
...
stddev(indicator_1),
stddev(indicator_2),
...
count(indicator_1),
count(indicator_2),
...
DECADE -- day converted to dekad in year (1-36)
from GRID_YIELD
where year between start year and end year -- day converted to yyyy format
group by crop_no, ddmm, grid_no -- day converted to mmdd format
LONG_TERM_AVERAGE_NUTS_YIELD: select crop_no,
mmdd, -- day converted to mmdd format
nuts_code,
average(indicator_1),
average(indicator_2),
...
min(indicator_1),
min(indicator_2),
...
max(indicator_1),
max(indicator_2),
...
stddev(indicator_1),
stddev(indicator_2),
...
count(indicator_1),
count(indicator_2),
...
DECADE -- day converted to dekad in year (1-36)
from NUTS_YIELD where year between start year and end year -- day converted to yyyy format group by crop_no, ddmm, nuts_code -- day converted to mmdd format |
