|
|
|
Agriculture &
Soil |
A primary study on crop production prediction using global vegetation index
Prediction Model Derivation
Mathematically, if these is an arrays [x1, y1], [x2, y2],
[xn. Yn], the regression line is y= ax+b, then the original y and regressive value y' have the relationship: y1 + y2+
..+yn = y' 1+2' 2+
.=y' n. Using this theory, the study area was divided into several subareas, and the crop production of study area could be estimated by predicting the crop production of every subarea and add them up if only the regression line for cop production can be decided.
- GVI and Crop Production
Using the average value of GVI from July to September, the study area were classified into 4 categories such as forest, agricultural area, grass land, water or desert by using the threshold classification method. For those agricultural areas of every province, the GVI volume of different time duration were calculated, which
can be expressed as follows : GVIVLMK: SIGMA (AREA)
(SIGMA
(TIME) PIXGVI) * PIXAREA
.(1) Where, GVIGLM: GVI volume of agricultural
area of certain time duration. PIXGVI: the monthly maximum value of GVI of
a certain agricultural pixel. PIXAREA; the area of
the same pixels. GVI data were usually considered as to represent the greenness of ground vegetation cover, therefore the GVI volume should be able to represent the amount go green biomass of ground vegetation cover, and the crop production should be related to the amount of green biomass in agricultural area. Since the original data had low resolution and the classification were carried out without ground
truth data, the result of classification can show the potential agricultural area but it can't represent the accurate cultivated area. In order to get better result, the GVI volume of agricultural are calculated former was calibrated using the area ratio between classified agricultural area and
real cultivated area. GVIVLM (MIDIF): Calibrated
GVI Volume AREAR: Cultivated
area AREAC: Classified agricultural area
The relationship between the calibrated GVI volume of different time duration and crop production of six province in the study area were studies and it was found that the crop production was highly related to the GVI volume of different time duration with all the correlation coefficient more than 90%. Fig. 1 shoes the relationships between the GVI volume from April to August and the crop production of that year.
- Weather Affection
From Fig. 1 it can be seen that the crop production has direct ratio with the calibrated GVI volume but the slop and segment of the linear equation for every year are different. Fig. 2 and Fig. 3 show the annual change of these tow values. If the straight line of GVI and crop production for certain year can be decided, the crop production of that year can be obtained. The weather affection was considered as the most important factors to decide the linear equation.
Average temperature, rainfall and aridity of different time duration were calculatied on province level. In order to reflect the range of affection of weather data, weight method was used. The weight of a province was decided by dividing the cultivating area of all six provinces. Table 1 shows the weight values can be used for prediction even without the cultivating area data of that year. The average temperature, rainfull and aridity of every province were multiplied by the weight of that province and added up to be the average rainfall and temperature of the study area of that year, and the relationship between these weighted weather factors and the coefficient of the crop-GVI linear function were studied. As the results, it has been fount that the weighted average temperature from June to August highly related both to the slope and segment with correlation coefficient 94% and 80% ( see Fig. 4 and Fig. 5 ). Therefore, It's possible to decide the crop-GVI equation by observing temperature from June to August.
- Prediction Processing
Using the results obtained up to now, the crop production prediction of the study ara can be carried out as follows:
Step 1: Collecting weather data from June to August and using the weight value in Table1, calculated the average temperature of this time duration, and decide the crop- GVI equation according to following equations:
SLOPE = 0.03623 * t - 0.06988------------------------(3)
SEGMENT = -127.6 * t + 2701.6------------------------(4)
where,
t : weighted temperature
The crop-GVI equation is:
Crop = SLOPE * GVIVLM + SEGMENT------------------(5)
Step 2: Calculate the GVI volume using equation (1) from June to August for every province and using the ratio of cultivating area and agricultural area, calibrate the calculated GVI volume using equation. (2)
Step 3: Add up the caliberated GVI volume of every province and introduce it to equation (5) to get the crop production prediction value of study area.
Conclusions
Fig. 6 shows the real and estimated value of crop production of study area. The average accuracy is approximately 95%.
The results obtained upto now show the possibility of using NOAA GVI data and weather data for large scale crop production prediction. Study area can be divided into several subareas and by deciding the regression line for crop production to GVI volume, the crop production of every subarea can be obtained. The summation of these predicted value forms the prediction value of the study area. The regression line can be determined by weather factors. The study should be improved in following aspects 1) the possibility of earlier prediction should be studies 2) Find out better correlationships which was used to decide the GVI crop regression line.
References
- Satellite Data Division, USA, Global Vegetarian
Index User's guide
- Xuemei Bai and Shunji Murai, (1989), "Regional
Environment Evaluation of Changjiang and Huang - Huai River Catchment
Area in China using global vegetation index and weather data" the 10th
Asian conference on Remote Sensing, Nov, 23-29 1989, Kualalumpur,
Malaysia.
- Shunji Murai and Xuemei Bai, (1990),
Habitability anaysis in China Using Global Vegetation Index, The 2nd GIS
workshop, Beijing, Aug. 6-8 1990, Beijing China.
| | | | | |
|
