Easy Leaf Area Automated Digital Image Analysis for Rapid and Accurate Measurement of Leaf Area1

Abstract

Accurate measurement of leaf area index (LAI) is important for agricultural analysis such as the estimation of crop yield, which makes its measurement work important. There are mainly two ways to obtain LAI: ground station measurement and remote sensing satellite monitoring. Recently, reliable progress has been made in long-term automatic LAI observation using wireless sensor network (WSN) technology under certain conditions. We developed and designed an LAI measurement system (LAIS) based on a wireless sensor network to select and improve the appropriate algorithm according to the image collected by the sensor, to get a more realistic leaf area index. The corn LAI was continuously observed from May 30 to July 16, 2015. Research on hardware has been published, this paper focuses on improved system algorithm and data verification. By improving the finite length average algorithm, the data validation results are as follows: (1) The slope of the fitting line between LAIS measurement data and the real value is 0.944, and the root means square error (RMSE) is 0.264 (absolute error ~ 0–0.6), which has high consistency with the real value. (2) The measurement error of LAIS is less than LAI2000, although the result of our measurement method will be higher than the actual value, it is due to the influence of weeds on the ground. (3) LAIS data can be used to support the retrieval of remote sensing products. We find a suitable application situation of our LAIS system data, and get our application value as ground monitoring data by the verification with remote sensing product data, which supports its application and promotion in similar research in the future.

Introduction

The leaf area index (LAI) is defined as half of all the leaf areas per unit surface area¹. It can provide quantitative information of structural description for vegetation canopy structure and energy exchange. As an important parameter of vegetation community growth analysis and ecosystem monitoring, LAI has been widely used in agriculture, forestry, ecology, biology, and other fields^2,3,4,5,6. Therefore, the acquisition of LAI has become essential at home and abroad, especially when researching the above fields. According to the space coverage, we divide the acquisition of LAI into two types: remote sensing satellite image retrieval and ground station measurement.

Satellite sensing provides products over a very large spatial scale, which is a huge advantage to do regional or global research. Song et al. used remotely sensed leaf area index (LAI) data to monitor spatial–temporal changes of winter wheat phenology in response to climate warming across the North China Plain, and they found high spatial heterogeneity of winter wheat phenology in pixel scale across the whole area, which could not be detected in previous site-based studies⁵. Reygadas et al. analyzed long time series of leaf area index data collected by the Moderate Resolution Imaging Spectroradiometer (MODIS) to determine the trend of forest degradation in Mexico⁷.

Although the satellite remote sensing inversion technology of LAI is becoming more and more mature and widely used, it is still indispensable to measure LAI data on the ground. On the one hand, it can provide more detailed and specific vegetation information on the local scale, on the other hand, the ground measurement data of LAI can be used as the validation data for its remote sensing inversion to improve the accuracy of LAI remote sensing products^8,9.

The LAI ground measurement methods include the direct measurement method and the indirect measurement method^10,11. The commonly used direct method is the LAI leaf length and width (LAILLW) method, which uses the sampling method to remove leaves from plants, measure single leaf area one by one and finally convert it into leaf area index¹². When the sample is suitably representative, the results of direct measurements are more accurate, and direct measurements are thus often regarded as the real leaf area index. However, the method inevitably causes damage to plants and it is with the time-consuming and labor-intensive character¹³. Therefore, indirect methods can be introduced. It reduces the damage caused to the vegetation and the environment, and it is also faster than the direct measurement method. According to the measurement, the principle used in the indirect method is usually divided into the radiometric-based method and the image-based method¹². The LAI2000 and LAI2200 are a series of instruments for analyzing plant canopy parameters and they are common devices for indirect measurement of LAI¹⁴. Optical sensors are used to measure the transmitted light above and below the canopy, and the leaf area index is automatically derived from a specific calculation model¹⁵. Although LAI2000 measurement is accurate and quick, it is not compatible with long-term monitoring of the spatial and temporal dynamics of leaf area development^12,16.

The development and application of wireless sensor network (WSN) technology provide a possibility for automatic continuous observation of vegetation structure parameters. Up to now, many researchers have done studies on the design and product validation of LAI measurement systems based on WSN. Due to the differences in system design and research areas, the conclusions of these verification works are also different. Qu et al. designed an automatic crop structure parameters measurement system¹⁷ based on WSN in which the parameters were calculated by measuring the canopy transmittance at different solar height angles. However, the data reliability of the system needs to be further tested, because the obtained data is only compared with the LAI2000, rather than the true value acquired via the direct measurement. Guo et al. created an empirical model¹⁸ that contains the LAI parameter through a large number of measurements in apple orchards, but some determination coefficients are less than 0.9 with low accuracy. Assessment of FAPAR from ESA Sentinel-2 using WSN found considerable discrepancies between image data and measurements especially in forest ecosystems¹⁹. These models measure fewer parameters and are only compared with the LAI value measured by the indirect method. WSN, as the data for rapid ground monitoring, also has certain reference results for the verification of remote sensing satellite inversion results. Therefore, the model designed in this paper not only collects the leaf area index, but also collects the radiative transfer rate, and uses a more accurate algorithm to calculate the value of LAI.

In addition, the Internet of things has also made some new progress in the field of plant monitoring. Sun²⁰ established an intelligent spacing selection model of Internet of things communication nodes under energy-saving constraints. Their hardware model improves the performance of the Internet of things system. Kundu et al.²¹ proposed Deep learning and IoT-based solutions which makes the model more suitable for automating disease detection and proved that the proposed model is effective in providing a low-cost and handy tool for farmers to improve crop yield and product quality. Dhaka et al.²² presented a survey of the existing literature in applying deep Convolutional Neural Networks to predict plant diseases from leaf images. They provide a platform for scholars in the field of applying deep learning techniques for the identification and classification of plant leaf diseases.The findings of these studies provide some ideas for our research.

As a result, a set of automatic, continuous, and low-cost leaf area index observation systems (LAIS), which was composed of the hardware structure of the sensor and software part of the image processing have been designed by our team, and the research results have been published in Sensors journal in 2015²³. Our previous research mainly focused on the sensor layout, data acquisition, transmission, and storage. In data processing, we simply adapted the finite length averaging method²⁴ to estimate LAI from images. However, the finite length averaging method is based on the assumption of Poisson leaf distribution, for crop canopy which has a non-Poisson structure and high vegetation coverage²⁵, it may cause underestimation of LAI up to 25%²⁶. So, as an improvement for the finite length averaging method, our recent research proposed a new formula to estimate LAI from a gap image. The hardware sensor system is also updated with a new technique to acquire the image of better quality. A series of experiments were carried out in Huailai experimental station, Hebei Province. The purpose of this experiment includes not only testing the performance of LAIS equipment and the processing algorithm of the image acquired by LAIS but also validating the remote sensing LAI products based on the ground LAI observation data. Moreover, the work carried out in Huailai experimental station also includes the ground measurement of radiation flux, surface temperature, albedo, and other parameters²⁷, as well as the verification of remote sensing products.

This paper mainly focuses on the verification between the data obtained by the upgraded system and the direct measurement methods and satellite remote sensing products, so as to get the applicability of WSN in different situations. The content of this study includes the acquisition of different types of data in the early stage and the accurate calculation of leaf area index–image processing algorithm. Besides, the verification and comparison between different data.By comparing with LAI2000 measurement and remote sensing LAI product, we try to find a suitable situation for the application of our LAIS system and demonstrate its potential for verification of remote sensing product.

Materials and methods

Study area

With the advanced observational techniques, abundant data accumulation, and ability to carry on multi-scale experiments, the Huailai Remote Sensing Station and around (for short Huailai Station), located in Huailai, Hebei province, China (40.349°N, 115.785°E), becomes one of the ideal study areas for the observation and validation of the LAI²⁷. The Huailai Station is mainly covered by corn and some weeds. So, we mainly use LAIS to monitor the growth cycle of corn (in April 2015, we submitted an application for plant collection permission to Huailai Remote Sensing Station and obtained approval.)

Huailai WSN vegetation monitoring system includes 6 sets of monitoring equipment, and its distribution is shown in Fig. 1 as follows, in which red dot represents LAIS Node, purple frame represents MODIS pixel, red frame represents observation area. The observation system is designed for the application of remote sensing pixel scale authenticity tests. The observation scale is a 1 km MODIS pixel on the pixel scale, and the actual coverage area is 2 km * 1.5 km. The six sets of equipment cover the core area of the test station and the surrounding typical growth plot, which is a good representative of the 1 km pixel scale.

Figure 1

Equipment distribution of WSN vegetation monitoring network in Huailai (red dot represents LAIS Node; purple frame represents the footprint of a MODIS pixel.

Full size image

Each piece of equipment consists of two cameras which were only one camera with two different angles in previous work²³ set up at a height of 2.5–4 m above the ground (Fig. 2), one for vertical downward observation and the other for inclined observation, which can take canopy photos regularly every day at its fixed position. The observation system obtained the photos of the corn canopy from May to August, but the corn did not grow in August. Therefore, in this study, we selected the photos taken by the vertical observation camera of the corn sample plot in the experimental station from May to July 2015.

Figure 2

The design of the LAIS node.

Full size image

Related work—data acquisition

Data collection using LAIS

The data collection complies with the plant guidelines statement: "LAI-2000 Plant Canopy Analyzer Instrution Manual" (Supplementary Information 2) (https://www.licor.com/env/, Last visit time: 21 October 2021). Existing facilities such as the high poles and the wireless sensor network in the experimental station have proved convenient for the installation of the LAI measurement system. LAIS uses the GEO001 digital serial camera that is suitable for a variety of embedded image acquisition modes. The specification of the camera includes: the total field of view is 120°, the maximum image size is 2176 * 1920 (approximately 5 million pixels), mounted at a height of 3 m, the spatial resolution at ground level is about 3 mm. The acquired image is simultaneously stored in a flash card in two formats: the JPEG format merits in less file size thus suitable for quick wireless transfer; the RAW format, which is the user data in our analysis, contains 3 channel binary image in 10 bits bit-depth. Compared to our previous work, an important new feature of this camera is the programmable cut-off filter. As we know, unlike scientific sensor which has the precise spectral response to each band, the digital camera is cheap and can only acquire the so-called RGB image. Usual digital cameras have one NIR cut-off filter to exclude the near-infrared light. The GEO001 camera, which was a commercial camera produced by Zhongshan Yunteng Photographic Equipment Co., Ltd, has two cut-off filters: one is the NIR cut-off filter, another is a blue cut-off filter. Switching on the NIR cut-off filter results in an ordinary color image as in a usual household digital camera. While the blue cut-off filter is switched on and NIR cut-off filter is switched off, near-infrared light is allowed to reach the detector array and blue light is blocked, resulting in false-color images as in Fig. 3b. Adding near-infrared light can increase illumination in the shadow area, and blocking blue light can alleviate the disturbance of sun glint, so, switching to a blue cut-off filter helps to improve the image quality when the direct sunlight is strong such as around noon time.

Figure 3

Three images on July 2 of site 1: (a) and (c) are true-color images obtained at 05:31 a.m. and 6:32 p.m., and (b) is a false-color image when the blue filter is removed at 1:28 p.m.

Full size image

To acquire an image in the best illumination condition and avoid the influence of rain or other unsuitable weather, the image acquisition device based on WSN was set up to acquire images three times per day: 5:30 a.m., 1:30 p.m., and 6:30 p.m. According to our experience, when the canopy is open (sparse vegetation), usually images acquired at 6:30 p.m. are the best for classification because the direct sunlight is weak; when the canopy is closed (dense vegetation), the illumination on the soil background is very poor in all time, and classification is difficult. So, the camera is programmed to switch to a blue cut-off filter when acquiring images at 1:30 p.m., while the images acquired at other times were with NIR cut-off filter, resulting in true color images, as shown in Fig. 3.

LAILLW data and LAI2000 data

To evaluate the accuracy of the improved finite length averaging method proposed in this study, a field experiment was carried out to measure LAI by manual sampling (Supplementary Information 3,4). A field sampling scheme covering the corn growing season (late May to early July) was designed (Supplementary Information 1). The LAI of corn in the experimental area was measured by the quadrat harvesting method, and the validation data of LAI of corn in each growth period were obtained. Considering the rapid growth of the corn, the sampling experiment period was set as 1 week, but due to the actual work in summer and the influence of rainfall, six effective measurements were carried out in the field experiment: May 30, June 7, June 13, June 20, July 4 and July 16.

The LAILLW method, which is also known as the shape factor method, involves outdoor and indoor measurements. The formulas are:

$${\text{L}} = {\text{S}}*{\text{N}}$$

(1)

$${\text{f}} = {{\text{S}} /{\left( {\sum\limits_{i = 1}^m {{\text{len}}*{\text{wid}}} } \right)}}$$

(2)

where L represents the leaf area index, S refers to the area of a single plant, and N refers to the number of plants in a unit area. The shape factor ƒ is the ratio of the S to the value multiplied by the length and width of all leaves in the plant.

To reduce measurement errors, 10 plants were selected in the sample, and the length and width of each leaf on each corn were recorded with a ruler. To obtain the shape factor, representative corn plants were cut next to the sample (not in the image coverage area) and the true area of each leaf was obtained by software, and the shape factor was derived from this²³. Through the length and width of 10 strains measured in the field, and the shape factor obtained, the total leaf area of 10 corns can be calculated, and the average leaf area of one plant is finally obtained. The LAI value under the LAILLW method is obtained.

Using the difference between the solar radiation values of the upper and lower canopies, the LAI2000 canopy analyzer can obtain LAI and set up a corresponding point folder to save the measured data for subsequent collation. 10 measurement points were selected for each site, and the average value was the final result for each site. To reduce the effects of the solar altitude angle on measurement accuracy, the experiments were repeated every two hours.

To make it easier to record the date of data acquisition, the data were summarized in the order day of the year (DOY). For example, 30 May 2015 is the 150th day in the year and its DOY is 150. The DOY information of data acquisition using the LAILLW method and LAI2000 is specifically shown in Table 1.

Table 1 The DOY information of data acquisition using the LAILLW and LAI2000.

Full size table

MODIS LAI data

MODIS leaf area index data was downloaded from the United States Geological Survey (https://modis.gsfc.nasa.gov/data/dataprod/mod15.php), named MCD15A2Hv006. It is an 8-day composite dataset with a 500-m pixel size. The algorithm chooses the best pixel available from all the acquisitions of both MODIS sensors located on NASA's Terra and Aqua satellites from within the 8 days.

In the comparison of MODIS LAI data, as the pixel of the satellite product is in 500 m resolution, it is not recommended to directly compare single node LAIS measurement with the MODIS LAI product because of the scale mismatch. Though complicated upscaling approaches have been discussed and implemented in Huailai station for other parameters²⁸, it is not the purpose of this study So, we simply averaged the LAI in all the LAIS nodes to compare to the average MODIS LAI product in the 3 * 3 nearest pixels (1.5 km * 1.5 km), referred to as MODIS LAI_Mean in a later context, which approximately covers the area of all LAIS nodes. Time matching was carried out by selecting the date of the MODIS product closest to the date of the handheld LAI2000 measurement. The following Table 2 is obtained by taking 3 * 3 pixels closest to the LAIS Nodes.

Table 2 MODIS leaf area index of 3 * 3 pixels around Huailai experimental station.

Full size table

Improved LAIS methods

In previous work, we have deployed sensors and cameras, and also have an automatic image processing and preliminary method of calculating LAI²³. Figure 4 is a flow chart of our work. The previous articles focused on hardware and system implementation but did not pay much attention to performance. On this basis, we upgrade the image classification method and LAI calculation method, which will be explained in detail below.

Figure 4

Flow chart of leaf area index measurement system based on WSN.

Full size image

Image preprocessing and classification methods

Because of weather-related factors such as water vapor and dust or inaccurate exposure, a small number of the photographs are not clear. Besides, some of the image data cannot be decoded because of unstable communications and other factors. Therefore, it is necessary to check and select the photographs that meet the processing requirements before binary image processing. Currently, the selection process is carried out by human visual inspection based on the following principles: (1) when the canopy is open (sparse vegetation), the image at 6:30 p.m. is preferred, when the vegetation the canopy is closed (sparse vegetation), the image at 1:30 p.m. is preferred; (2) if the preferred image is not clear, other clear image acquired on the same day should be used; if all the images are not clear, then this day is marked as a failure.

If we decided to use the image acquired at 1:30 p.m. It is also necessary to convert it from a false-color image to a true-color-like image (as shown in Fig. 3b) in which the leaves are shown in green color. The conversion is carried out by multiplying the vector of DN (digital number) of 3 bands with a coefficient matrix which is provided by the camera manufacturer. Another preprocessing is to choose the near nadir-view area of the image for further processing. As the off-nadir-view area of the image is subject to large geometric distortion as well as saturation of fraction of vegetation cover (FVC), they are not used in this study. The images are clipped to an ROI (region of interest) of about 2 * 2 square meters in ground area, with a maximum view zenith angle less than 30°.

The study of the color spatial distributions of the crop images is helpful for the classification of the images and extraction of the image information. The color of the image pixel is the most direct and effective element that can be used to describe the image²⁹. Because the red–green–blue (RGB) color space has the characteristic of a clear and convenient expression of information. When corn leaves are small, the crops in the fields are sparse, and most of them are soil background in the images. The soil in a lower hue is similar to the corn in terms of R and B components, while it has an overlap with the corn in G components when soil is in a higher hue. This makes it difficult to classify sparse corn scenes only by RGB space, so it is necessary to consider the characteristics of hue, luminosity, and saturation (HLS) spatial components.

Statistical analysis showed that the component values of the crop leave in the RGB color space were in the ranges of G > R and G > B while the corresponding values for the soil follow the law that B < G < R³⁰. In the HLS color space, the H component has a specific distribution law³¹. Therefore, the distributions of R, G, B, and H components were used as an important criterion for crop image classification³².

In this study, two threshold classification methods were proposed, as shown in Table 3 The default thresholds (t1, t2, t3, t4, and t5) were recommended by the system when operating each method.

Table 3 The relationship of two threshold classification methods.

Full size table

The condition "t1 < H < t2" corresponds to the typical range of H component of green leaves in HLS color space, and "S > t3" is for the occasion when the junction of green leaves and soil in the H component is not very high. The intersection of two conditions can accurately distinguish the general green leaves. "G > t4" can distinguish green leaves with high illumination and "R < t5 and B < t5" can obtain green leaves in the case of dry soil background. The three conditions can get higher classification accuracy. In the experiments, the crop canopy images that were obtained when the leaves are small and the corn is sparse in the field were suitable for use with method 1 and the first four corresponding default thresholds were 80, 160, 18, and 240. When "t5 = − 1" the soil is moist and when "t5 = 40" the soil is a little dry.

When the leaves were larger and they occluded each other, changes in the chlorophyll content and enhancement of the sun-light resulted in a low contrast between the leaf and the background in the true-color images when the growth period reached two months into the panicle stage, so it was no longer appropriate to use classification method 1. The default thresholds in method 2 were 5, 5, − 1, 80, and − 1. False-color images are acquired in raw format. First, need to convert the RAW format to BMP format. Then replace the order of the image bands of G and B in professional software such as ENVI.

In the classification of crop canopy images, only the sub-region in the center of the photograph was used because of the geometric distortion of the camera photographs, and the corresponding area in each binary image had dimensions of approximately 2 × 2 m. Figure 5 shows the original and binary images of site 1 at 1:30 p.m. and 6:30 p.m. on July 1. By visually comparing the original image with the binary image (Fig. 5), it can be seen that the classification result of the false-color image is better when the leaves are large and occluded, the contrast between the lower leaves and the soil is low. This is also another significant finding in this validation experiment and it's of practical significance for the application of the proposed system to farmland monitoring.

Figure 5

Original images and binary images of site 1 on July 1: (a) the original image at 1:30 p.m.; (b) preprocessing result of image (a); (c) the binary image of (a); (d) original image at 6:30 p.m.; (e) binary image of (d).

Full size image

During the application, it is necessary to select the appropriate image as well as classification methods by the practical situation such as the crop type, the crop growth stage, and the image quality. We have mentioned in "Data collection using LAIS" section that when the canopy is sparse, we prefer to use the image at around 6:30 p.m.; if the illumination on soil is too poor as the result of the dense canopy, we prefer to use the false-color image at around 1:30 p.m., and the above color space transformation is required. Sometimes the weather is bad and the preferred images are not of good quality, we need to choose one of the three images at 5:30 am 1:30 p.m., and 6:30 p.m. according to the situation. Currently, this choice needs human visual inspection, but it is a quick and simple operation, and won't waste much time as in the field measurement. Although most of the image processing is automated, human supervision is still essential, as the image quality is important to distinguish leaves from soil background, which is a prerequisite for calculating LAI. Figure 6 shows the binary images for the different dates in site 1.

Figure 6

The time series of extracted leaf-cover from the digital images in site 1: (a) May 30; (b) June 7; (c) June 13; (d) July 4; (e) July 16; (f) August 1.

Full size image

The improved finite length averaging method for LAI estimation

The classification discriminates pixels of green leaf from soil background in the LAIS acquired image. So, the fractional vegetation coverage (FVC) and gap probability, which equals 1-FVC, can be derived by dividing the number of leaf pixels by the total number of pixels in the ROI, assuming the change of view zenith angle within the ROI can be neglected. Theoretically, LAI can be related to FVC through the widely used Beer-Lambert law.

However, the original Beer-Lambert law only applies to uniformly distributed and infinitesimal leaves. For the real canopy, leaf angle distribution and clumping index (CI) should be considered^1,24. The finite length averaging method was proposed in 1986 to simultaneously estimate LAI and CI from field measurement of canopy gaps. It also applies to gap data generated from canopy photography and still is the recommended method to estimate LAI and CI for crop canopy up to now. The formula can be summarized as:

$$L = - \frac{\cos (\theta )}{{mG(\theta )}}\sum\limits_{i = 1}^{m} {\ln \left( {P_{i} (\theta )} \right)}$$

(3)

$$\Omega = \frac{{m\ln \left( {\frac{1}{m}\sum\nolimits_{i = 1}^{m} {P_{i} (\theta )} } \right)}}{{\sum\nolimits_{i = 1}^{m} {\ln \left( {P_{i} (\theta )} \right)} }}$$

(4)

where L denotes LAI and \(\Omega\) denotes CI; \({P}_{i}(\theta )\) is the gap probability when observation zenith angle is \(\theta\), which is 0 in this case; \(G(\theta )\) is related to leaf angle distribution and is normalized projection area of leaf in the observation direction; the footnote \(i\) denotes the \({i}{th}\) sample line/sample rectangle, and \({\text{m}}\) is the total number of sample lines/sample squares. When dealing with gap data generated from canopy photography, usually the sample square is adopted instead of the sample line.

However, Lang and Xiang²⁴ also pointed out two limitations with this approach: (1) If the gap probability is equal to 0, then the estimation is meaningless as infinity occurs. This is often the case of a high LAI. (2) the method assumes that the leaf size should be sufficiently small relative to the side of the sample square. If the sample square is too small, then the precondition cannot be satisfied, and usually, an overestimation of LAI occurs.

As the camera of LAIS is fixed on the top of a pole, the acquired sample image is of limited area. Then we face the problem of either using a small sample rectangle or the total number of sample rectangles is insufficient. To solve the problem of a small sample square, we proposed empirical formula to replace the log function in characterizing the relationship between gap probability and LAI in the sample square based on computer simulations³³. The new formulas correct the shortcomings of over-estimation and instability of log function when the canopy is dense and the side length of the sample square is short. The revised formula is:

$$\mathrm{L}=\frac{1}{m}\sum_{i=1}^{m}\frac{cos{\theta }_{i}}{G\left({\theta }_{i}\right)}f({P}_{i}\left({\theta }_{i}\right),D)$$

(5)

where D denotes the equivalent leaf length, which is defined as the square root of the average area of the single leaf; and \(f\) is the proposed empirical formula, in the form of:

$$f\left(P,D\right)=\left(1-P\right){P}^{1/{A}_{1}}+\frac{\left(1-{P}^{1/{A}_{2}}\right){(1-P)}^{1/{A}_{3}}}{1/{A}_{4}+{A}_{5}{P}^{{A}_{6}}}$$

(6)

$${\mathrm{A}}_{i}={a}_{i1}+{a}_{i2}(W/D)+{a}_{i3}\mathrm{log}(W/D) i=\mathrm{1,2}\dots ,6$$

(7)

where \({a}_{ij}\) (\(\mathrm{i}=1,...,6; j=\mathrm{1,2},3\)) are empirical coefficients with values in Table 4. W is the side length of the sample squire. Our study also found that the optimal setting for the side length of the sample square is about 3 times of equivalent leaf length in most cases of crop or grass scenes.

Table 4 The coefficients in \(f\left(P,D\right)\) estimation formula for fusiform leaf, using sample square.

Full size table

Results

Verification and comparison of LAIS data and ground monitoring data

The LAI that was extracted using the improved algorithm and the fractional vegetation coverage (FVC) are shown in Fig. 7. The FVC refers to the percentage of the vertical projection area of vegetation in the total sample area. After one and a half months of corn planting, LAI growth accelerated (see Fig. 7) which is by the typical growth cycle of corn. The WSN-based system increases the time–frequency of data acquisition and detailed analysis of the growth of the crop can be achieved. A small fluctuation was observed in the LAI curve, which was a normal phenomenon. Excessive illumination and a lack of water in the leaves are the main causes of this phenomenon. Both FVC and LAI are unitless ratio values, so are plotted in one picture (Fig. 7).

Figure 7

Change curve of LAI and FVC in site 1.

Full size image

The LAI data of the area under study can be acquired using a combination of images classification and the improved algorithm. The data obtained by the LAILLW method were considered as true values. A comparison of values obtained using the LAIS, the LAI 2000, and the LAILLW methods is shown in Fig. 8. The results showed that the LAIS data and the LAI2000 data were both in good agreement with the true values using the LAILLW method. The slope of the fitting line of the LAIS is closer to 1 (slope = 0.944), with a very small offset, and the root-mean-square error (RMSE) is smaller (RMSE = 0.264) than the LAI2000. When the LAI is large, there is a significant underestimation of the LAI2000 measurements (see Fig. 8b).

Figure 8

Scatter distribution of the LAI estimation: (a) the LAIS; (b) the LAI2000.

Full size image

Table 5 shows the acquired LAI data using the LAILLW method and the absolute errors between the other two methods (LAI2000 method and the LAIS method) and the true value. As can be seen from the table, in the early stage of corn growth, the LAIS measurement error is smaller than the LAI2000. The LAI2000 is to calculate the LAI based on the difference in radiance values above and below the canopy, and the LAI2000 is not very sensitive to radiation measurements when the corn leaves are small. Therefore, the measurement of LAI2000 in the early stage of corn growth is not very accurate, and the LAIS method is to obtain LAI by image processing, which can avoid such defects. The main problem of the LAIS method at this stage is that the measured value is higher than the true value. The main reason for this phenomenon is that there are weeds in the acquired image, so the area of the extracted crop is overestimated, and the estimated LAI value is larger than the true value. This factor is not considered when using the LAILLW method for field measurements. When the LAI value is greater than 1, there is a significant underestimation of the LAI2000 measurement, while the LAIS estimate is relatively stable.

Table 5 The LAILLW data and the absolute errors between LAI2000 and LAIS data.

Full size table

Validation of remote sensing product with LAIS data

We plotted the ground measurement data together with the large-scale MODIS LAI data in Fig. 9. We can see from the figure that the change curves of LAIS and LAI2000 almost coincide with the real leaf area index (LAI_LAILLW), while the MODIS data is higher than our real value. It is because the remote sensing pixel value is the average of the region, and the growth trend of corn is not the same in different plots of the region, and even there are other crops. So, the verification of remote sensing products is not easy to do. The more observation points on the ground, the better. WSN just provides this possibility. So, our work in the future is to make the LAI observation of WSN more convenient.

Figure 9

Variation curves of different LAI products in different observation periods.

Full size image

Discussion

In this paper, based on the conventional LAI measurement method, a leaf area index (LAI) sensor which can automatically and continuously be monitored in real-time is developed, and two threshold algorithms are designed according to the growth stage of corn in different periods. Qu,y et al. also realized the automatic acquisition of leaf area index with the help of WSN system¹⁷, but they did not collect the image of the leaf, and the validation data is also the data of LAI2000 and did not use LAILLW (the true value of the LAI), so the results still need to be further verified. The measurement principle of our system is to obtain the image of the leaf itself, which is closer to the real value than the measurement principle of LAI2000, and the result also proves that this is the case. Similarly, Bauer, Jan et al. presented a modification of commercial PAR sensors³⁴, which is also compared with the measurement data of LAI2200, but no image is used to obtain the leaf area. However, in our experiments, the leaf area index was calculated by image classification, which is more suitable for the small leaf of corn. The crop canopy images that were obtained when the leaves are small and the corn is sparse in the field were suitable for use with method 1 ("Image preprocessing and classification methods" Method1) .When the leaves were larger and they occluded each other, changes in the chlorophyll content and enhancement of the sun-light resulted in a low contrast between the leaf and the background in the true color images when the growth period reached two months into the panicle stage, it was no longer appropriate to use classification method 1, so we use method 2 ("Image preprocessing and classification methods" Method2).The improvement of image processing algorithm makes us get more accurate leaf area and character factor.

However, it is a disadvantage that human operation is needed to select the best quality image for classification. In the future, it is expected that this visual inspection process can be modeled to achieve full automation for data processing. Compared with the traditional field measurement methods, it is fast and accurate. The system based on a wireless sensor network not only saves time and labor but also increases the time and frequency of data acquisition, which can realize the detailed analysis of crop growth. At the same time, our verification results show that the measurement error of LAIS is less than LAI2000 in the early growth stage of corn. Although the result of our measurement method will be higher than the real value, it is due to the influence of weeds on the ground. When the LAI value is greater than 1, the measured value of LAI2000 is seriously underestimated, while the estimated value of LAIS is relatively stable. We also used LAI2200 in the experiment, but because LAI2000 was used twice, to keep the data consistent, we use LAI2000 to name them. Their measurement principle and accuracy are the same. Perhaps the new LAI instrument has been improved, but as a system that can automatically obtain crop images and monitor crop growth status in real-time, our measurement system and method have great advantages, because it not only obtains more detailed growth information but also saves manpower and material resources. What's more, our results are reliable, especially at the early stage of corn growth.

Another important issue is the bidirectional reflectance effect of vegetation canopy, which has a subtle impact on LAI estimation from satellite image^35,36 and ground-based measurement as well. Part of the bidirectional reflectance effect comes from the solar illumination angle. The LAIS image taken at 1:30 p.m. is of a small solar zenith angle, the direct sunlight can better penetrate dense vegetation canopy to illuminate soil background; on the other hand, strong direct light may cause glare in leaf and become a disturbance in the classification process. So, when the vegetation is sparse, we prefer to use the LAIS image taken at 6:30 p.m. when the sky diffuse light dominates. Anyway, as the LAI estimation is based on the binary images as in Fig. 6, the effect of illumination angle will not affect the final result so long as the classification is successful. Another part of the bidirectional reflectance effect comes views angle. If the camera takes image off-nadir-view, then a view angle correction must be applied to the gap probability \({P}_{i}(\theta )\), and leaf angle distribution, which is represented by \(G(\theta ),\) should be considered. Our strategy in this study is to fix the camera in nadir view, and use only the part of the image with a view angle less than 30°. However, the off-nadir-view image can also be used in LAI estimation. The optimal view angle and its directionality in LAI estimation should be discussed in our future research.

Conclusions

The work in this paper improved the measurement of LAI to an extent and the data of the previously developed system is compared and verified. The system verification was performed by comparing the results with those from LAILLW measurements. A comparative analysis of measurement results from the LAIS and commonly used instrument LAI2000 was also performed. At the same time, using the data obtained by the LAIS to validate the remote sensing products of the leaf area index, it is concluded that the wireless sensor network monitoring system can provide test values for large-scale spatial data. Our data results are also true and reliable and can be applied to LAI inversion in the future.

This experiment improves the finite length average method, which improves the accuracy of the calculated LAI. The slope of the fitting line of the LAIS is 0.944, which displays a highly consistent with the true LAI that was acquired using the LAILLW method. The RMSE is 0.264 of the LAIS, smaller than that of LAI2000 (RMSE = 0.373). The absolute error is 0–0.6 and the relative error of most LAIS data is less than 30%. The classification results for the false-color image containing near-infrared information are better than those that are obtained from the true color image when the blades of corn are large and form the occlusion. The overestimation of LAI in the system is largely caused by the existence of grass. The strong sunlight is the main factor, reducing the discrimination between the soil and vegetation, and thus producing an underestimation when the leaves are in small size.

The outcome that the LAIS showed a better performance than LAI2000 supports its use to the application and popularization on farmland. With the verification of MODIS LAI, we conclude that the leaf area index measurement system based on WSN can provide test values for large-scale spatial remote sensing data. Of course, our system also needs to be improved. The further accuracy of the image recognition algorithm will have an important impact on the accuracy of the leaf area index.

References

1. Chen, J. M. et al. Leaf area index of boreal forests: Theory, techniques, and measurements. J. Geophys. Res. Atmos. 102(D24), 29429–29443 (1997).

   ADS  Article  Google Scholar

2. Fang, H. et al. Continuous estimation of canopy leaf area index (LAI) and clumping index over broadleaf crop fields: An investigation of the PASTIS-57 instrument and smartphone applications. Agric. For. Meteorol. 253–254, 48–61 (2018).

   ADS  Article  Google Scholar

3. Jin, H. et al. Intercomparison and validation of MODIS and GLASS leaf area index (LAI) products over mountain areas: A case study in southwestern China. Int. J. Appl. Earth Obs. Geoinf. 55, 52–67 (2017).

   ADS  Article  Google Scholar

4. Melnikova, I. et al. Estimation of leaf area index in a mountain forest of central Japan with a 30-m spatial resolution based on landsat operational land imager imagery: An application of a simple model for seasonal monitoring. Remote Sens. 10(2), 179 (2018).

   ADS  Article  Google Scholar

5. Song, Y. et al. Using MODIS LAI data to monitor spatio-temporal changes of winter wheat phenology in response to climate warming. Remote Sens. 12(5), 786 (2020).

   ADS  Article  Google Scholar

6. Zhao, X. et al. Retrieval and application of leaf area index over China using HJ-1 data. Geomat. Nat. Haz. Risk 8(2), 478–495 (2016).

   Article  Google Scholar

7. Reygadas, Y., Jensen, J. L. & Moisen, G. G. Forest degradation assessment based on trend analysis of MODIS-leaf area index: A case study in Mexico. Remote Sens. 11(21), 2503 (2019).

   ADS  Article  Google Scholar

8. Fang, H. et al. Validation of global moderate resolution leaf area index (LAI) products over croplands in northeastern China. Remote Sens. Environ. 233, 111377 (2019).

   ADS  Article  Google Scholar

9. Xu, B. et al. An integrated method for validating long-term leaf area index products using global networks of site-based measurements. Remote Sens. Environ. 209, 134–151 (2018).

   ADS  Article  Google Scholar

10. Olivas, P. C. et al. Comparison of direct and indirect methods for assessing leaf area index across a tropical rain forest landscape. Agric. For. Meteorol. 177, 110–116 (2013).

    ADS  Article  Google Scholar

11. Černý, J., Haninec, P. & Pokorný, R. Leaf area index estimated by direct, semi-direct, and indirect methods in European beech and sycamore maple stands. J. For. Res. 31(3), 827–836 (2018).

    Article  Google Scholar

12. Yan, G. et al. Review of indirect optical measurements of leaf area index: Recent advances, challenges, and perspectives. Agric. For. Meteorol. 265, 390–411 (2019).

    ADS  Article  Google Scholar

13. Jiapaer, G. et al. Comparison of non-destructive LAI determination methods and optimization of sampling schemes in an open Populus euphratica ecosystem. Urban For. Urban Green. 26, 114–123 (2017).

    Article  Google Scholar

14. Vincent, G. et al. Mapping plant area index of tropical evergreen forest by airborne laser scanning. A cross-validation study using LAI2200 optical sensor. Remote Sens. Environ. 198, 254–266 (2017).

    ADS  Article  Google Scholar

15. Sumnall, M. et al. Assessing the transferability of statistical predictive models for leaf area index between two airborne discrete return LiDAR sensor designs within multiple intensely managed Loblolly pine forest locations in the south-eastern USA. Remote Sens. Environ. 176, 308–319 (2016).

    ADS  Article  Google Scholar

16. Bauer, J. et al. Processing and filtering of leaf area index time series assessed by in-situ wireless sensor networks. Comput. Electron. Agric. 165, 104867 (2019).

    Article  Google Scholar

17. Qu, Y. et al. Design and experiment of crop structural parameters automatic measurement system. Trans. Chin. Soc. Agric. Eng. 28(2), 160–165 (2012).

    Google Scholar

18. Guo, X. M. et al. A model with leaf area index and apple size parameters for 2.4 GHz radio propagation in apple orchards. Precis. Agric. 16(2), 180–200 (2014).

    Article  Google Scholar

19. Putzenlechner, B. et al. Validation of Sentinel-2 fAPAR products using ground observations across three forest ecosystems. Remote Sens. Environ. 232, 111310 (2019).

    ADS  Article  Google Scholar

20. Sun, J.-S., Zhu, T. & Wozniak, M. Intelligent spacing selection model under energy-saving constraints for the selection of communication nodes in the Internet of Things. Mob. Netw. Appl. (2021).

21. Kundu, N. et al. IoT and interpretable machine learning based framework for disease prediction in pearl millet. Sensors 21(16), 5386 (2021).

    ADS  Article  Google Scholar

22. Dhaka, V. S. et al. A survey of deep convolutional neural networks applied for prediction of plant leaf diseases. Sensors 21(14), 4749 (2021).

    ADS  Article  Google Scholar

23. Li, X. et al. The design and implementation of the leaf area index sensor. Sensors 15(3), 6250–6269 (2015).

    ADS  Article  Google Scholar

24. Lang, A. R. G. & Xiang, Y. Q. Estimation of leaf-area index from transmission of direct sunlight in discontinuous canopies. Agric. For. Meteorol. 37(3), 229–243 (1986).

    ADS  Article  Google Scholar

25. Chen, J. M. & Black, T. A. Foliage area and architecture of plant canopies from sunfleck size distributions. Agric. For. Meteorol. 60(3–4), 249–266 (1992).

    ADS  Article  Google Scholar

26. Hu, R. et al. Indirect measurement of leaf area index on the basis of path length distribution. Remote Sens. Environ. 155, 239–247 (2014).

    ADS  Article  Google Scholar

27. Dou, B. et al. Wireless sensor network of typical land surface parameters and its preliminary applications for coarse-resolution remote sensing pixel. Int. J. Distrib. Sensor Netw. 12(4), 9639021 (2016).

    Article  Google Scholar

28. Wu, X. et al. Coarse scale in situ albedo observations over heterogeneous snow-free land surfaces and validation strategy: A case of MODIS albedo products preliminary validation over northern China. Remote Sens. Environ. 184, 25–39 (2016).

    ADS  Article  Google Scholar

29. Su, J. et al. Machine learning-based crop drought mapping system by UAV remote sensing RGB imagery. Unmanned Syst. 8(1), 71–83 (2020).

    Article  Google Scholar

30. Qu, Y., Liu, S. & Xia, J. Image processing methods to determine the fractional vegetation cover of winter wheat using digital camera. Arid Land Geogr. 33(6), 997–1003 (2010).

    Google Scholar

31. Yang, W. et al. Greenness identification based on HSV decision tree. Inf. Process. Agric. 2(3–4), 149–160 (2015).

    Google Scholar

32. Guijarro, M. et al. Automatic segmentation of relevant textures in agricultural images. Comput. Electron. Agric. 75(1), 75–83 (2011).

    Article  Google Scholar

33. Liu, Q. et al. Improvement of limited length average method in field measurement of leaf area index. Chin. Bull. Bot. 53(5), 671–685 (2018) (in Chinese).

    Google Scholar

34. Bauer, J. et al. On the potential of Wireless Sensor Networks for the in-situ assessment of crop leaf area index. Comput. Electron. Agric. 128, 149–159 (2016).

    Article  Google Scholar

35. Zhang, S. et al. Development of a new BRDF-resistant vegetation index for improving the estimation of leaf area index. Remote Sens. 8(11), 947 (2016).

    ADS  Article  Google Scholar

36. Qi, J. et al. Leaf area index estimates using remotely sensed data and BRDF models in a semiarid region. Remote Sens. Environ. 73(1), 18–30 (2000).

    ADS  Article  Google Scholar

Download references

Acknowledgements

We would also like to thank the Huailai Science Experiment Station in the Aerospace Information Research Institute.

Funding

This work was supported by National key research and development plan subject "Watershed Non-point Source Pollution Prevention and Control Technology and Application Demonstration Project" (2021YFC3201500) and the National Natural Science Foundation of China (No. 41476161), the WSN application extended of the open fund of the state laboratory of remote sensing science (No. OFSLRSS201626) and the Fundamental Research Funds for the Central Universities (No. 312231103).

Author information

Authors and Affiliations

1. State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences, No. 8, Da Yang Fang, An Wai, Chao Yang District, Beijing, 100012, China

   Rongjin Yang

2. State Key Laboratory of Remote Sensing Science, College of Global Change and Earth System Science, Beijing Normal University, No.19, Xinjiekou Wai Street, Haidian District, Beijing, 100875, China

   Lu Liu, Qiang Liu, Xiuhong Li, Xuejie Hao, Yushuang Ma & Qiao Song

3. Higher Institute of Computer Modeling and Their Applications, Clermont Auvergne University, Clermont Auvergne, Clermont-Ferrand, France

   Lizeyan Yin

Contributions

R.Y., X.L. and Q.L. wrote the main manuscript text; L.L., Y.L. and Y.M. prepared figures. X.H., Q.S. and Y.M. prepared tables. All authors reviewed the manuscript.

Corresponding author

Correspondence to Xiuhong Li.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

About this article

Cite this article

Yang, R., Liu, L., Liu, Q. et al. Validation of leaf area index measurement system based on wireless sensor network. Sci Rep 12, 4668 (2022). https://doi.org/10.1038/s41598-022-08373-z

Download citation

• Received: 28 November 2021

• Accepted: 28 February 2022

• Published: 18 March 2022

• DOI : https://doi.org/10.1038/s41598-022-08373-z

Further reading

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

englishwaakeen.blogspot.com

Source: https://www.nature.com/articles/s41598-022-08373-z

Share this post