Difference between revisions of "Calculation of long term average crop indicators"
From Agri4castWiki
| Line 1: | Line 1: | ||
{{Scientific}} | {{Scientific}} | ||
| − | The calculation of long term averages is done for each crop indicator and for each individual geographical feature on a 10 daily base. | + | The calculation of long term averages is done for each crop indicator and for each individual geographical feature on a 10 daily base. A dedicated software package called [[LTA_YIELD]] is developed to do all calculations. |
For the following basic parameters average values are calculated: | For the following basic parameters average values are calculated: | ||
Revision as of 11:27, 17 June 2011
The calculation of long term averages is done for each crop indicator and for each individual geographical feature on a 10 daily base. A dedicated software package called LTA_YIELD is developed to do all calculations.
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
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
|
