Researches of spectral
feature and growing monitoring of rice
Cha Zhongxing Institute of Land survey and
Planning of China
Wang Yanyi Jiangsu Academy of
Agrcultural Sciences
Qu Boolin Institute of Land survey
and Planning of China
Abstract Using the
method of power regression analysis, we find correlations between the
Ratio Vegetation Index D/C (infrered to red, obtained from paddy spectral
data) and LAI or TDM are very high. Before rice heading , the correlation
coefficients (c.c.) are 0.96 and 0.93 respectively; after rice heading,
the c.c. between D/C and TDM is -0.90 (5 varieties, 173 samples). Using
the regression model obtained from one year to forecasting the LAI and TDM
of another year, the precisions are: 96.4% and 94.6% respectively before
rice handing, after rice heading it is 95.4% for forecasting TDM.
Introduction Rice is one of the major crops in China.
There are a great number of rice varieties and their growing periods and
plant types are very different. So it is necessary to create a method to
monitor the rice growth for the paddy management and yield estimation. In
order to predict rice yield, agronomists usually measuse LAI and TDM
during the rice growth period, but their methods are laborious and
time-consuming.
With the development of remote sensing technology,
it has been recognized that the crop spectral feature (such as light
absorption, transmission and reflection) are closely determined by crop
physiological characters which can express growing vigour and yield
components of rice. From this point of view, the research was carried out
the explore the relationship between rice spectral parameter and agronomic
parameters during rice growing period and to monitor dynamically the rice
growth.
Experiment Method Experiments were carried out
in Jiangsu Academy of Agricultural Sciences, Nanjing, 1987-1989,with 4
rice varieties (Yangan No.2-mediumjaponica, Xiu Shui No.4-late japonica,
Nanjing No.11-medium indica, Shanyou No.63-hybrid medium indica), 3
fertilization levels and 2 replications. The 0.2 hectare experimental
field was divided into 24 plots. Rice seedlings were transplanted on June
10th every year. Rice plant height, leaf area index(LAI), total dry
biomass (TDM), tiller developments and phonological phases were measured
every 15 days from the 7th day after transplantation.
The Exotech
100 Radiometer used in the experiment has 4 spectral bands: a 0.45-0.52um,
b:052-0.60um, C: 0.63-0.69um, D:0.76-0.90um ( the same as the first 4
bands of TM) . All spectral data were collected by a polycordes on
selected clear, windless days at 10:00-14:00 Binjing time. The Radiometer
is set perpendicularly downward, 2m above the rice canopy. Five testing
points were set in every plots. The average value of 5 tests was regarded
as the spectral reflectance of the paddy plots. The average value of 5
tests was made once for every 1 or 2 weeks during the whole rice growing
period.
Results and discussion
- Fig. 1 shows that the paddy spectral features of blue, green and red
bands are: low reflective, low transmissive and highly absorbed. The
main reason is that the pigment especially the chlorophy11) strongly
absorbs the photon. In red band the reflectance is very low because more
than 90% of incident sun light energy is absorbed for rich
photosynthesis. During the period from seedling transplantation to
heading, long with the increase of rich leaf photosynthesis, the
chlorophy11content in rice plant increases rapidly. With the arrival of
rice heading time, rice reflectance in red bands increases. After milk
stage, decreases but the rice panicle needs more nutrient, therefore the
rice reflectance in red band increases more rapildly.
Fig. 1 Paddy reflectance
character
There is a little reflectance peak in the
green band and it increases with the growth of rice. Because the inner
tissue of the rice leaf has high reflectance and high transmission in
the near infrared (NIR) band, so there is a high reflectance and low
absorption in the NIR band Fig. 2 shows that he NIR reflectance
increases with the development during the growing period of the rice.
Fig. 2 Near infrared reflectance of
rice
From seedling transplantation to heading stage,
the reflectance of NIR increases with the increase of LAI, then a
tendency of steadiness follows. During milk stage, along with large
amount of nutrient transferred into the panicle and the change of inner
tissue of the leaf, the reflectance of NIR decreases (Fig. 3)
Fig. 3 Relationship between LAI and NIR
reflectance.
Agronomist had indicated that the value
of rice biomass depends and LAI, and the former is closely related to
yield. It is, therefore, reasonable to rely on R and NIR bands in
studying the relationship between paddy spectrum and rice growth.
- Relationship between paddy spectrum and LAI or TDM
Initial
analysis of data shows that the relationship between rice reflectance of
single band and LAI or TDM is not very ideal. But by using regression
analysis form R and NIR band with 4 spectrum combinations with 6
mathmodels,, we have selected the best spectrum combination (D/C) and
the best correlation model F( Y=AXB).
Because the reflectance
before and after rice heading are different in data analysis, the rice
growing period is divided into 2 parts.
Tab. 1 shows the poer
regression coefficient between D/C and LAI or TDM is highest (0.961 and
0.928 respectively).
Using the D/C in 1988 and 1989 as one
correlated variables and the LAI or TDM in 1988 and 1989 as the other
correlated variables, we got the power regression coefficients of
individual variety and total varieties of rice. Tab. 2 shows, before
rice heading, correlation between D./C and LAI as well as TDM are
significant for either individual variety and total varieties.
Using the D/C of 4 rice varieties in 1988 or 5 rice varieties in
1989 to make regression analysis with LAI or TDM (before rice heading,
82 samples), the results are:
In 1988 |
YL = 0.301
X0.865 |
( r = 0.900) ......(1) |
YT = 8.684
X1.100 |
( r = 0.892) ......(2) |
In 1989 |
YL = 0.345
X0.800 |
( r = 0.961) ......(3) |
YT= 12.44
X0.992 |
( r = 0.928) ......(4)
| ( X=D/C, YL = LAI and
YT = TDM )
Table. 1 The comparison of 6 regression analysis between
LAI or TDM and 4 composition of bands in 1989
LAI or TDM |
Composition of bands |
Mathmode1 |
Y = AX+B |
Y = (A / X)+B |
Y = 1 / A+BX |
Y = X / A+BX |
Y=AeBX |
Y = AXB |
Correlation Coefficient |
L A I |
D/C |
0.654 |
-0.669 |
-0.648 |
0.880 |
0.772 |
0.961 |
D-C |
0.585 |
-0.716 |
-0.592 |
0.821 |
0.774 |
0.832 |
(D-C) / (D+C) |
0.665 |
-0.880 |
-0.614 |
0.868 |
0.905 |
0.889 |
C / (A+B+C) |
-0.715 |
0.868 |
0.716 |
-0.695 |
-0.912 |
-0.880 |
T D M |
D/C |
0.684 |
-0.672 |
-0.620 |
0.873 |
0.764 |
0.928 |
D-C |
0.575 |
-0.708 |
-0.612 |
0.813 |
0.767 |
0.814 |
(D-C) / (D+C) |
0.704 |
-0.873 |
-0.625 |
0.863 |
0.900 |
0.871 |
C / (A+B+C) |
-0.711 |
0.853 |
0.702 |
-0.643 |
-0.903 |
-0.840 |
Table 2. Power correlation coefficient between D/C and
TDM or LAI
Varceties |
Before rice heading in 1988 |
Before rice heading in 1989 |
before rice heading in 1988 & 1989
|
After rice heading in 1988 |
After rice heading in 1989 |
SAM.* |
LAI |
TDM |
SAM. |
LAI |
TDM |
SAM. |
LAI |
TDM |
SAM |
TDM |
SAM. |
TDM |
B.P.** |
A.P.** |
B.P. |
A.P. |
Yangen No. 2 |
23 |
0.895 |
0.876 |
23 |
0.987 |
0.945 |
40 |
0.942 |
0.915 |
30 |
-0.678 |
-0.912 |
27 |
-0.689 |
-0.933 |
Nianjing No. 11 |
18 |
0.944 |
0.918 |
19 |
0.965 |
0.944 |
37 |
0.949 |
0.926 |
12 |
-0.304 |
-0.902 |
20 |
-0.487 |
-0.900 |
Shanyou No. 63 |
18 |
0.931 |
0.903 |
20 |
0.961 |
0.931 |
38 |
0.925 |
0.901 |
23 |
-0.854 |
-0.914 |
29 |
-0.691 |
-0.930 |
Xiu Shiu No. 4 |
23 |
0.932 |
0.927 |
18 |
0.949 |
0.937 |
41 |
0.940 |
0.911 |
30 |
-0.416 |
-0.912 |
14 |
-0.738 |
-0.916 |
ZaoDan 8 |
11 |
0.964 |
0.900 |
12 |
-0.374 |
-0.831 |
Total Varceties |
82 |
0.900 |
0.892 |
91 |
0.961 |
0.928 |
173 |
0.933 |
0.903 |
95 |
-0.473 |
-0.904 |
103 |
-0.541 |
-0.901 | * SAM. : Samples; ** B.P.:
Before standardized processing A.P.: After Standardized processing
When using the total of 1988 and 1989 (before rice
heading 173 samples) to make the regression analysis with D/C and LAI or
TDM, the reaults are:
Y = 0.354
X0.794 |
(r =
0.933)......(5) |
Y = 10.982 X1.018 |
(r =
0.903)......(6) |
From
the above results of regression analysis, it is found that power
correlation between D/C and LAI or TDM before heading is not limited by
rice varieties.
After rice heading, the reflectance from paddy
canopy is the mixed one from leaves and panicles of rice. Especially
after the beginning of milk stage, because of differences in speed of
grain filling and changes of the ratio of leaf and panicle, the spectral
reflectances of different plots diverse greatly and the correlation
between D/C and TDM (single variety or total varieties) is very low. But
using following formula to make standardized process, the result is
satisfactory .
Eq. 7
(X1' is the value
after standardized processing, Xi is observed value and K is the number
of observations after rice heading) Using 103 samples (after rice
heading) to make the regression analysis between D/C and TDM, the
correlation reached to extreme significance level (Tab. 2) and the
formula as follows;
Y = 0.924
X0.741 |
(r =
-0.901)......(8) |
(X
is D/C and Y is TDM)
- Test of mathematical model
Putting the graded D/C average
value of 5 varieties in 1989 (before rice heading , 91 samples) into
formula 1 and formula 2 (from 1988) and calculating the predicted value
of LAI and TDM, the average rediction precision is 96.4% for TDM
(Tab.3).
Grading the processed D/C value in 1988 (98 samples,
after rice heading) into 5 classes, and putting it to formula 8, the
prediction precisions for TDM are 91.2% to 98.5% (compare with the
observed TDM value), the average precision is 95.4% (Tab.4)
Table.3 The comparison between predicted value and
observed value of LAI and TDM before rice heading.
D / C |
range |
< 10 |
10-20 |
20-30 |
30-40 |
< 50 |
samples |
20 |
22 |
14 |
25 |
4 |
average value |
4.56 |
19.81 |
24.29 |
34.75 |
41.57 |
LAI |
predicted value |
1.12 |
3.89 |
4.75 |
6.48 |
7.41 |
observed value |
1.16 |
3.56 |
4.94 |
6.40 |
7.32 |
precision |
96.5% |
91.1% |
96.1% |
99.4% |
98.8% |
TDM |
predicted value |
48.09 |
231.90 |
290.20 |
430.30 |
510.21 |
observed value |
52.56 |
255.68 |
298.62 |
416.14 |
522.12 |
precision |
91.1% |
90.3% |
97.1% |
96.7% |
97.7% |
Table 4. The
comparison between predicted value and observed value of TDM after
rice heading
D / C |
range |
< 0.5 |
0.5-1.0 |
1.0-1.5 |
1.5-2.0 |
< 2.0 |
S.P.V. |
samples |
6 |
10 |
18 |
56 |
8 |
TDM |
predicted |
698.89 |
783.64 |
864.06 |
988.61 |
1116.22 |
observed value |
763.22 |
802.13 |
876.89 |
1019.38 |
1040.05 |
precision |
91.5% |
97.7% |
98.5% |
96.9% |
92.9% |
- The relationship between TDM and rice yield
Making linear
regression analysis between the TDM value after rice heading and rice
yield from 16 plots in 1989, the result is
Y = 0.707X +
315.069 |
(r =
0.803)......(9) |
(X is
average TDM after heading, Y is rice yield)
If biomass
of each variety is measured according to its growing stages, perhaps the
correlation between average TDM and yield will be much more
;significant. Conclusion
- Because there is some inner relationship between the reflectance of
red band and the chlorophy11 content of rice leaf and between the
reflectance of NIR band and the LAI of rice, D/C is very successful in
reflecting the rice growth.
- Before rice heading, the power correlation (. C.) between D/C and
LAI or TDM is much closer than other kinds of combination. After rice
heading, when standardized process is made, the P.C between D/C and TDM
reaches to a very significance level.
- Differences of rice varieties can be neglected in setting the
correlation model between D/C and LAI or TDM.
- It is quite possible to use "D/C - LAI" and "D/C - TDM" models in
monitoring rice growth.
- Rice yield can be predicted indirectly by suing the data of paddy
spectrum.
References
- N. K. Patll etc. The Relationship between the Rice spectral
Reflection and Rice Yield. Int J. Remote Sensing, 6(5), 1985.
- E. W. Lemaster etc. Seasonal Inspection of Model of Wheat Suits
Spectral Albedo, Photogrammetic Engineering and Remote Sensing, 9,
1980.
|