Next Article in Journal
A Cross-Layer, Anomaly-Based IDS for WSN and MANET
Next Article in Special Issue
Measuring the Human Ultra-Weak Photon Emission Distribution Using an Electron-Multiplying, Charge-Coupled Device as a Sensor
Previous Article in Journal
Determining Forest Duff Water Content Using a Low-Cost Standing Wave Ratio Sensor
Previous Article in Special Issue
Surface Plasmon Resonance and Bending Loss-Based U-Shaped Plastic Optical Fiber Biosensors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optical Method for Estimating the Chlorophyll Contents in Plant Leaves

by
Madaín Pérez-Patricio
1,
Jorge Luis Camas-Anzueto
1,*,
Avisaí Sanchez-Alegría
1,
Abiel Aguilar-González
3,
Federico Gutiérrez-Miceli
1,
Elías Escobar-Gómez
1,
Yvon Voisin
2,
Carlos Rios-Rojas
1 and
Ruben Grajales-Coutiño
1
1
Tecnológico Nacional de México, 29000 Tuxtla Gutiérrez, Chiapas, Mexico
2
Le2i FRE2005, CNRS, Arts et Métiers, Univ. Bourgogne Franche-Comté, av. des Plaines de l’Yonne, BP16, 89010 Auxerre CEDEX, France
3
Institut Pascal, Université Clermont Auvergne, 63178 Clermont Ferrand, France
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(2), 650; https://doi.org/10.3390/s18020650
Submission received: 3 January 2018 / Revised: 16 February 2018 / Accepted: 18 February 2018 / Published: 22 February 2018
(This article belongs to the Special Issue Optical Biochemical Sensor Systems and Applications)

Abstract

:
This work introduces a new vision-based approach for estimating chlorophyll contents in a plant leaf using reflectance and transmittance as base parameters. Images of the top and underside of the leaf are captured. To estimate the base parameters (reflectance/transmittance), a novel optical arrangement is proposed. The chlorophyll content is then estimated by using linear regression where the inputs are the reflectance and transmittance of the leaf. Performance of the proposed method for chlorophyll content estimation was compared with a spectrophotometer and a Soil Plant Analysis Development (SPAD) meter. Chlorophyll content estimation was realized for Lactuca sativa L., Azadirachta indica, Canavalia ensiforme, and Lycopersicon esculentum. Experimental results showed that—in terms of accuracy and processing speed—the proposed algorithm outperformed many of the previous vision-based approach methods that have used SPAD as a reference device. On the other hand, the accuracy reached is 91% for crops such as Azadirachta indica, where the chlorophyll value was obtained using the spectrophotometer. Additionally, it was possible to achieve an estimation of the chlorophyll content in the leaf every 200 ms with a low-cost camera and a simple optical arrangement. This non-destructive method increased accuracy in the chlorophyll content estimation by using an optical arrangement that yielded both the reflectance and transmittance information, while the required hardware is cheap.

1. Introduction

Due to escalating population growth, food crop demand has increased. To obtain higher-yielding crops, several food producers constantly monitor the crop´s health since it is well-known that stress conditions affect photosynthetic activity, which is closely related to yield. There are variables related to crop health such as chlorophyll fluorescence, chlorophyll, calcium, and nitrogen contents. One variable often used and which has obtained satisfactory results is the chlorophyll content in the leaves. This is due to the high correlation found between chlorophyll content and the health of plants [1]. Several methods for estimating chlorophyll content can be found in the current literature and are based on the transmittance or reflectance of the leaf; nevertheless, chloroplast arrangement in the cells is modified by the intensity, color, and duration of the incident light, which produces variations in the values obtained with measurement devices. Therefore, the analyzed plant must be placed under light-controlled conditions before being measured in order to define the chloroplast arrangement [2] as this yields more accurate results.
Devices such as the SPAD 502 (Soil Plant Analysis Development) of Minolta have been used to estimate chlorophyll contents [3,4,5,6,7]. Its main advantage is that it can provide fast estimations with high accuracy. In general, the SPAD measures the leaf absorbance between 650 and 940 nm [2]. Based on the measured absorbance, a SPAD numerical value proportional to the chlorophyll content is then computed. The main disadvantage of the SPAD system is that it only estimates the absorbance at one point of the leaf under analysis, calculating the chlorophyll content only within a small spatial location on the leaf. To solve this problem, iterative measurements at different spatial locations must be performed. Then, the average value of all the measured points corresponds to the approximated chlorophyll content in the plant leaf being analyzed.
A vision-based approach is frequently used for chlorophyll content estimation. In this context, vision-based techniques estimate the chlorophyll content using information obtained via image processing algorithms, some of them capturing images using airborne cameras [8,9]. These techniques have the advantage of monitoring large land areas; however, their cost is high, and there is no control for ambient lighting changes that might decrease their accuracy. There are other approaches that capture images at different wavelengths. In particular, multispectral and hyperspectral images taken in a range of 300–1200 nm have been used to analyze the health of food crops [10,11,12,13,14,15,16,17,18,19]. In general, multispectral/hyperspectral-based approaches can obtain relatively high accuracy and real-time processing; unfortunately, they have a high cost and large size.
Other trends that aim to reduce cost and processing time have been previously proposed by researchers [20,21,22,23,24,25] who have used images with visible spectrum (300–700 nm) information and a single RGB (Red-Green-Blue) color space camera to determine the correlation between chlorophyll and nitrogen content. Several of these approaches have been tested in controlled conditions, and are not well-suited for field conditions. In [26], a robotic system was used to transport the camera along the crop, and five images of the leaves were obtained to reduce noise. The mean of the RGB values was used as a reference to determine lettuce calcium deficiency via statistical computations.
In several cases, multiple images are required making the procedures not suitable for real-time processing. For example, Wiwart et al. [27] required three images of successive nodes for chlorophyll estimation, while Pagola et al. [4] required four images. Finally, high accuracy and relatively low processing time were achieved using neural networks as proposed by Odabas et al. and Vesali et al. [28,29], where the leaf being analyzed was placed in front of a smartphone camera and the color transmittance was used as the base for the chlorophyll estimation.
The objective of this work was to propose an accurate, low-cost, and real-time approach for chlorophyll estimation. In this paper, reflectance and transmittance in three broad bands (R,G,B) were proposed as the base parameters. There are two contributions attributed to this work: the first is a novel mathematical formulation that uses reflectance and transmittance for estimating chlorophyll content that delivers robust and accurate chlorophyll estimations for hydroponic food crops; the second is a novel optical arrangement capable of capturing base parameters (reflectance/transmittance). The proposed approach, which reached high accuracy with low cost, could be a promising solution for the current industrial plantation procedures. Experimental results showed that the proposed algorithm reached over R2 = 0.97 for a Lactuca sativa L. hydroponic crop and achieved a fast estimation of the chlorophyll contents in the plant leaf every 200 ms.
The rest of this manuscript is organized as follows: Section 2 presents the proposed approach in detail; the experimental results and a comparison to previous works are detailed in Section 3; and finally, Section 4 concludes this manuscript.

2. Materials and Methods

To validate the proposed optical configuration, a data set was prepared consisting of different plants from the same crop, each with different health states and different chlorophyll content. A Lactuca sativa L. hydroponic crop was studied because of its rapid growth. On this crop, a nutritional solution was applied to obtain different health levels for the plants within the crop. A solution able to provide the daily nutritional macronutrient requirements is typically composed of 4.0810 mg of nitrogen, 0.5531 mg of phosphorus, and 0.4881 mg of potassium. Twenty-seven nutritional solutions were prepared, with the concentration of macronutrients varying among them. For the first one, 100% of the daily requirements of each nutrient were applied: that is, 4.0810 mg of nitrogen, 0.5531 mg of phosphorus, and 0.4884 of potassium. For the second and third solutions, nitrogen was reduced to 2.0405 mg and 0.0 mg—corresponding to 50% and 0% of the daily requirement, respectively—while the concentrations of phosphorus and potassium were maintained at 0.5531 mg and 0.4884 mg, respectively. This procedure was repeated for each component of the solution. Each nutritional solution was applied to three different plants, and the experiment was replicated six times. The hydroponic system was placed in a greenhouse with a controlled temperature. The temperature was set to 26 °C during the day and 18 °C at night. For the first 10 days, only water was applied. From Day 11 to 60, 450 mL of solution was applied every third day. Finally, leaf samples were taken after 50 days, when the plant had reached maturity and deterioration due to age was minimal. Chlorophyll content estimation was performed using a SPAD meter (SPAD 502) that computes the optical density difference at two wavelengths: 650 nm and 940 nm. The measurement area of the device was 2 mm × 3 mm. Eight measurements were taken from different zones of the leaf and the mean of measurements was retained as the chlorophyll content value.
For the second experiment, fresh leaves of Lactuca sativa L., Azadirachta indica L., Lycopersicum sculentum L., and Canavalia ensiforme L. were collected. Thirteen leaves grouped in three classes according to their greenness, from yellowish to dark green, were used. In a mortar, 0.5 g of fresh plant material was placed and macerated. The whole process was carried out in a place with low luminosity. From the macerate, 4 mL of 99% acetone was mixed with 2 mL ethanol (2:1 v/v), placed in 10 mL tubes, and mixed for 1 min stirring, ensuring complete contact of the plant material. These were then left to stand for 30 min in the freezer in the dark, and centrifuged for 10 min at 2000 rpm. They were covered with aluminum foil and 5 mL of acetone/ethanol (2:1 v/v) was added and stirred for 1 min. Absorbance readings were performed at wavelengths of 663 nm and 645 nm. The control was acetone/ethanol (2:1 v/v).
The obtained values were substituted in the following formulas, described in [30], for the estimation of photosynthetic pigments.
Chlorophyll a (mg/g) = (12.7 * A663) − (2.59 * A645)
Chlorophyll b (mg/g) = (22.9 * A645) − (4.7 * A663)
Chlorophyll total (mg/g) = (8.2 * A663) + (20.2 * A645)
where A663 and A645 are the absorbance measured from 663 nm and 645 nm, respectively. The spectrophotometer was adjusted to zero using the acetone/ethanol mixture.

2.1. Acquisition Technique and Base Parameters

Previous work has demonstrated that there are two crucial issues in the chlorophyll estimation procedures. The first is related to the acquisition technique used to capture base parameters. As reported in [2], several acquisition parameters such as the type of light emitted, intensity, or duration from the source and the leaf side (adaxial or abaxial) being illuminated can affect the reflectance and transmittance data. In general, acquisition techniques sensitive to noise yield poor results in terms of accuracy. In addition, environmental conditions such as rain or dust can affect their performance [31]. The second issue relates to the base parameters used in the chlorophyll estimation process. Robust base parameters should deliver accurate results. Nevertheless, previous vision-based approaches have demonstrated that it is difficult to obtain robust parameters for chlorophyll content estimation [32].
For the acquisition technique, previous vision-based approaches have reached low-cost solutions, while compromising accuracy. This is because accuracy is highly related to the environmental conditions during the image acquisition procedures. Therefore, computer vision systems using natural illumination tend to deliver poor performance in terms of accuracy. In this work, we proposed that one alternative to decrease noise (induced by the natural illumination) could be through the use of artificial light sources with known parameters of the light spectrum, thus improving the performance of the chlorophyll estimation. As a result, a novel portable device suitable for chlorophyll estimations was proposed, which uses controlled illumination conditions that help reduce the noise. Furthermore, chlorophyll estimation is performed within a closed environment, thus eliminating environmental perturbations. This means that the presented device is capable of provide health indicators with low noise sensitivity and without environmental perturbations. Prior to the image acquisition stage, all the samples were placed in darkness to define the chloroplast arrangement and to reduce the error in the chlorophyll content estimation.
For the base parameters, previous works have proposed several indices that relate chlorophyll content to the color components in the leaf being analyzed [5,19]. In this work, two different base parameters dependent on the value of each channel of color reflected or transmitted in the RGB space were proposed. Both parameters were obtained simultaneously by using a relatively simple optical configuration. As a result, it was possible to obtain accurate and economic chlorophyll content estimations.

2.2. Image Acquisition

To capture the base parameters, an optical arrangement was proposed, as seen in Figure 1. A color camera (model DBK 31AU03) was used to provide a 1024 × 768 Bayer pattern image. In a Bayer image, only one of the three color components is obtained. The remaining components are computed by interpolation. A Sony ICX204AK sensor shot the image in the RAW format in the ‘grbg’ Bayer pattern with 8 bits per pixel.
The camera sensor was a Sony ICX204 with dimensions of 5.8 × 4.92 mm2, and the lens used had an 8 mm focal length. The image acquisition device used a MCWHD2 led lamp (800 mW) as the light source. The nominal wavelength spectrum was 400–700 nm, and the maximum irradiance measured at a distance of 200 mm was about 121 µmol m−2s−1. In addition, the optical arrangement used a 50 × 70 mm2 glass that held the leaf to be analyzed. Finally, at the bottom, a 70 mm square mirror reflected the leaf reflectance information. Black-matte plastic pieces were used to control the light flow inside the device and reduce the entry of external light, which would otherwise induce measurement errors. The camera was placed at a 39° horizontal inclination (to observe simultaneously both the upper and lower part of the leaf), with a 100 mm vertical distance and a horizontal gap of 20 mm with respect to the glass in the center. The viewing angle—which determines the longitude of the visual field as a function of the distance between the camera and the glass where the leaf is placed—is calculated by using Equation (4), where α is the vertical or horizontal viewing angle, s is the vertical or horizontal dimension, and f is the focal distance. Using this configuration, both the glass and the mirror could be observed. As a result, the bisector d (Figure 1) represents the necessary distance between the vision field and the camera, and is calculated by Equation (5), where d is the bisector and l is the visual field longitude. The camera was focused on the visual field line, and the diaphragm opening was adjusted to f  =  f/8. Finally, to acquire the images, a small part (about 2 cm) of the lettuce leaves were inserted in the optical system and the images were acquired in RAW format, as shown in Figure 2.
α = 2 × tan 1 ( s / ( 2 × f )
d = l / ( 2 × tan ( α 2 ) )
As can be observed in Figure 2, the obtained image shows both the adaxial leaf side in the left part of the image and the abaxial leaf side in the right part of the image. The left part of the image produces the transmittance information, while the right part yields the reflectance information. The shape of both leaf sides is not similar due to the perspective distortion. A non-flat leaf must be placed between two glasses to help introduce the leaf in the device.

2.3. Image Processing

The image processing algorithm is outlined as follows.
  • Convert the Bayer image to an RGB image I(x,y,3).
  • Calibrate the system using the Macbeth Color Checker table.
  • Compute the binary image (Ib(x,y)) from the green color using active contours.
  • Compute the reflectance information by using the right part of the Ib’(x,y) image, and the transmittance information by using the left part of the Ib’(x,y).
Each one of the steps was developed as follows: the camera produced a Bayer image where only one color channel was provided for each pixel, and the other two color channels were obtained by interpolation. In this case, the input Bayer image (Figure 2) was processed by the Homogeneous edge direct algorithm [33] to obtain an RGB image, as shown in Figure 3. The color calibration of the camera was performed by using a linear model and the Macbeth Color Checker table as described in [34].
To separate the leaf from the background a method based on active contours has been used [35]. In Figure 4, the results of the binarization are shown. The blue line specifies the initial state of the active contour and the boundaries of the object are indicated in red and yellow color.
The right part of Figure 3 contains the reflectance information (Ibr(x, y)) and the left part contains the transmittance (Ibt(x, y)). The proposed method computed the reflectance (Rc) and transmittance (Tc) as shown in Equations (6) and (7), where c  =  [R, G, B], I(x, y, c) is the input RGB image (Figure 3), and Ibr(x, y), Ibt(x, y) are the binary images corresponding to the reflectance and transmittance, respectively; as seen in Figure 4b.
R c = ( 1 c ( I ( x , y , c ) ×   I b r ( x , y ) ) ) / n
T c = ( 1 c ( I ( x , y , c ) ×   I b t ( x , y ) ) ) / n  

2.4. Chlorophyll Content Estimation by Linear Regression

An association between the base parameters (Rc and Tc) and chlorophyll content was obtained by linear regression, providing high performance in terms of accuracy and processing time. In general, the main contribution of this work was the mathematical formulation, which used reflectance/transmittance as the base parameters, as well as basic image processing algorithms and a simple linear regression method. As a result, this formulation should provide accurate measurements with low cost and a compact system design.

3. Results

MATLAB 2015b was used to perform all image processing steps and the linear regression. Figure 5 and Figure 6 present the chlorophyll contents in SPAD values in relation with the average values for transmittance and reflectance values for each of the R, G, and B channels.
On one hand, as can be seen in Figure 5, reflectance decreased in all channels as the amount of chlorophyll increased. This was because higher chlorophyll contents tends to absorb more energy in the form of light. On the other hand, Figure 6 shows that higher chlorophyll contents reduced the transmittance value for all the channels. This was due to the higher chlorophyll contents which reduced the light that passed through the leaf.
To build a chlorophyll estimation framework, any statistical regression model has to fulfill the linear dependence of chlorophyll with respect to the reflectance/transmittance indicators. In this work, linear regression was used because of its mathematical simplicity, which involves high processing speed with low computational requirements. It was tested with different combinations of the base parameters proposed in this work. It was demonstrated that all the proposed parameters delivered high accuracy in terms of chlorophyll estimation; in particular, the combination of all of them provided accuracy superior to most of the previous vision-based approaches, and were similar to the SPAD-based measurements. Table 1 shows the results of using the simple linear regression models for each value showed in Figure 5 and Figure 6. The resulting expression for estimating the SPAD values is
SPAD = 40.87 − 14.59 × Tr
The maximum R2  =  0.94 was for the red channel transmittance Tr, and the smallest standard deviation obtained SD  =  1.19 SPAD was also for Tr. Even if the normalized root mean square error (NRMSE) value for Tr was 0.2817, the best variable for estimating the chlorophyll content was the transmittance in the red channel.
In Table 2, the results of applying the multiple linear regression models (using two variables) for each value presented in Figure 5 and Figure 6 are shown. The resulting expression for estimating the SPAD values is
SPAD = 40.34 + 4.14 × Rr − 17.65 × Tr
The best R2  =  0.97 was for the combination of reflectance and transmittance in the red channel (Rr, Tr). In addition, this combination had the best standard deviation SD  =  0.83 SPAD. However, the best NRMSE  =  0.25 was for reflectance and transmittance in the green channel (Rb, Tb), which had an R2  =  0.90, and a SD close to the unit. Since multiple linear regression requires more computational power than simple linear regression, it was possible to conclude that a single variable model was one of the best options for chlorophyll estimations in terms of simplicity. Simple linear regression using transmittance as a health indicator combined both relatively high accuracy (R2 = 0.93) and a very simple mathematical formulation.
Thirteen plants were also used to evaluate the proposal where the chlorophyll content was determined using a spectrophotometer, as shown in Table 3. The mean and the standard deviation of the chlorophyll content, the R2, and the NRMSE of the chlorophyll estimation using linear regression are shown.
Relationships between the chlorophyll content and transmittance values were determined by linear regression yielding Equations (10)–(18). Tr has been used because it has been demonstrated previously that it yields the most reliable performance to chlorophyll content estimation.
Chlorophyll a (Canavalia ensiforme leaves) = 51.64 − 161.71 × Tr
Chlorophyll b (Canavalia ensiforme leaves) = 19.63 − 49.51 × Tr
Total chlorophyll (Canavalia ensiforme leaves) = 69.97 − 203.01 × Tr
Chlorophyll a (Azadirachta indica leaves) = 27.21 − 49.03 × Tr
Chlorophyll b (Azadirachta indica leaves) = 10.76 − 20.31 × Tr
Total chlorophyll (Azadirachta indica leaves) = 37.97 − 69.35 × Tr
Chlorophyll a (Lycopersicon esculentum leaves) = 114.41 − 862.14 × Tr
Chlorophyll b (Lycopersicon esculentum leaves) = 37.95 − 281.80 × Tr
Total chlorophyll (Lycopersicon esculentum leaves) = 152.31 − 1143.50 × Tr
Like other methods based on computer vision, the equation used for the estimation of chlorophyll for each crop is different. On the other hand, the values of R2 obtained for Azadirachta indica and Lycopersicon esculentum leaves showed that this proposal was adequate for chlorophyll content estimation (a, b, and total). Nevertheless, the lowest values of R2 obtained for the Canavalia ensiforme leaves showed that non-flat leaves introduced errors into the chlorophyll content estimation. The performance of the chlorophyll content estimation can be improved by using a more robust method to select pixels used in the estimation process, or by implementing a method to increase flatness of the leaf.

3.1. Processing Speed and System Size

In Table 4, the processing speeds for the different versions of the proposed algorithm are presented. Clearly, the increase in the number of variables used in the regression model also increased the processing time. Nevertheless, the proposed algorithm reached real-time processing. Regarding the system size, the device used was as illustrated in Figure 1, which was a portable and small device (100 mm × 200 mm). Currently, the device must be connected to a computer that hosts the image processing algorithms and the chlorophyll content estimation algorithm. Due to the simplicity of the proposed approach, however, it can be implemented in either a small processor embedded in the optical device, or in a smart camera. The developed device can then maintain the same size, and be a promising solution for portable chlorophyll content estimation.

3.2. Comparison with Previous Work

Let F (Rr, Tr) be the best performance in terms of accuracy. In Equation (9), the corresponding adjusted model is presented, where F is the estimated chlorophyll value while Rr and Tr are the base parameters in the mathematical formulation. In Table 5, the accuracy comparisons with respect to previous vision-based approaches are shown. This method has been compared to previous works that have used SPAD measurements as reference values. In the case of multiple crops in the same approach, we considered the best R2 value. As observed, the proposed measurement approach outperformed most of the previous vision-based approaches because previous formulations used base parameters that had low stability under illumination changes. For this work, given the wide availability of low-cost RGB cameras, the proposed solution could be a promising chlorophyll estimation framework with a similar accuracy and size as SPAD, but with a lower cost and lower processing time (near 200 ms). As a work in progress, we will apply the proposed approach to other hydroponic crops (Raphanus sativus, Coriandrum sativum, and Brassica oleracea).

4. Conclusions

In this work, a novel chlorophyll estimation framework was proposed based on an optical arrangement capable of capturing base parameters (reflectance/transmittance). It was demonstrated that reflectance/transmittance delivered robust and discriminant values for the chlorophyll content estimations. The proposed optical arrangement captured the reflectance and transmittance changes in the leaf, analyzing them simultaneously and with low cost. The experimental results were encouraging: they demonstrated that a multiple linear regression algorithm using two variables achieved 97% accuracy for the hydroponic crop analyzed using a chlorophyll estimation method such as SPAD, while reaching estimated values close to the real values obtained with the spectrophotometer. In addition, the proposed approach delivered fast measurements with low cost and allowed a compact system design. To determine if other crop types present similar patterns, the mathematical model can be applied to several different food crops. In any case, the acquisition device and structure of the learning algorithm have to remain the same, therefore, only a color comparison of the crop being analyzed and an adjustment of the multiple regression models are needed.
For the developed device, the initial formulation was made by using the full area of the sample. However, the experimental results demonstrated that a similar performance could be obtained with small regions of the leaf being analyzed. Thus, the measurement device size could be smaller than in the presented version. As a work in progress, a smart camera is currently being developed. The chlorophyll values could then be displayed on a compact device, which would be useful for farmers and the current industrial plantation procedures.

Acknowledgments

This work was supported by grant No. 5660.15-P from the Tecnológico Nacional de México.

Author Contributions

Madaín Pérez-Patricio, Federico Gutiérrez-Miceli, Yvon Voisin and Jorge Luis Camas Anzueto conceived and designed the experiments; Avisaí Sánchez-Alegría and Abiel Aguilar-González performed the experiments; Carlos Rios-Rojas, Elías Escobar-Gómez and Ruben Grajales-Coutiño analyzed the data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Muñoz-Huerta, F.R.; Guevara-Gonzalez, G.R.; Contreras-Medina, M.L.; Torres-Pacheco, I.; Prado-Olivarez, J.; Ocampo-Velazquez, R.V. A Review of Methods for Sensing the Nitrogen Status in Plants: Advantages, Disadvantages and Recent Advances. Sensors 2013, 13, 10823–10843. [Google Scholar] [CrossRef] [PubMed]
  2. Nauš, J.; Prokopová, J.; Řebíček, J.; Špundová, Y.M. SPAD chlorophyll meter reading can be pronouncedly affected by chloroplast movement. Photosynth. Res. 2010, 105, 265–271. [Google Scholar] [CrossRef] [PubMed]
  3. Chang, S.X.; Robison, D.J. Nondestructive and rapid estimation of hardwood foliar nitrogen status using the SPAD-502 chlorophyll meter. For. Ecol. Manag. 2003, 181, 331–338. [Google Scholar] [CrossRef]
  4. Pagola, M.; Ortiz, R.; Irigoyen, I.; Bustince, H.; Barrenechea, E.; Aparicio-Tejo, P.; Lamsfus, C.; Lasa, B. New method to assess barley nitrogen nutrition status based on image colour analysis: Comparison with SPAD-502. Comput. Electron. Agric. 2009, 65, 213–218. [Google Scholar] [CrossRef]
  5. Hu, H.; Zhang, G.; Zheng, K. Modeling Leaf Image, Chlorophyll Fluorescence, Reflectance from SPAD Readings. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 4368–4373. [Google Scholar] [CrossRef]
  6. Fernandez-Jaramillo, A.A.; de Jesus Romero-Troncoso, R.; Duarte-Galvan, C.; Torres-Pacheco, I.; Guevara-Gonzalez, R.G.; Contreras-Medina, L.M.; Herrera-Ruiz, J.; Millan-Almaraz, J.R. FPGA-based chlorophyll fluorescence measurement system with arbitrary light stimulation waveform using direct digital synthesis. Measurement 2015, 75, 12–22. [Google Scholar] [CrossRef]
  7. Wan, W.; Hua, D.; Le, J.; He, T.; Yan, Z.; Zhou, C. Study of laser-induced chlorophyll fluorescence lifetime measurement and its correction. Measurement 2015, 60, 64–70. [Google Scholar] [CrossRef]
  8. Saberioon, M.M.; Amin, M.S.M.; Anuar, A.R.; Gholizadeh, A.; Wayayok, A.; Khairunniza-Bejo, S. Assessment of rice leaf chlorophyll content using visible bands at different growth stages at both the leaf and canopy scale. Int. J. Appl. Earth Obs. Geoinform. 2014, 32, 35–45. [Google Scholar] [CrossRef]
  9. Delegido, J.; van Wittenberghe, S.; Verrelst, J.; Ortiz, V.; Veroustraete, F.; Valcke, R.; Samson, R.; Rivera, J.P.; Tenjo, C.; Moreno, J. Chlorophyll content mapping of urban vegetation in the city of Valencia based on the hyperspectral NAOC index. Ecol. Indic. 2014, 40, 34–42. [Google Scholar] [CrossRef]
  10. Noh, H.; Zhang, Q. Shadow effect on multi-spectral image for detection of nitrogen deficiency in corn. Comput. Electron. Agric. 2012, 83, 52–57. [Google Scholar] [CrossRef]
  11. Zhang, X.; Liu, F.; He, Y.; Gong, X. Detecting macronutrients content and distribution in oilseed rape leaves based on hyperspectral imaging. Biosyst. Eng. 2013, 115, 56–65. [Google Scholar] [CrossRef]
  12. Li, G.; Wan, S.; Zhou, J.; Yang, Z.; Qin, P. Leaf chlorophyll fluorescence, hyperspectral reflectance, pigments content, malondialdehyde and proline accumulation responses of castor bean (Ricinus. communis L.) seedlings to salt stress levels. Ind. Crops Prod. 2010, 31, 13–19. [Google Scholar] [CrossRef]
  13. Zhang, Y.; Chen, J.M.; Miller, J.R.; Noland, T.L. Leaf chlorophyll content retrieval from airborne hyperspectral remote sensing imagery. Remote Sens. Environ. 2008, 112, 3234–3247. [Google Scholar] [CrossRef]
  14. Diezma, B.; Lleó, L.; Roger, J.M.; Herrero-Langreo, A.; Lunadei, L.; Ruiz-Altisent, M. Examination of the quality of spinach leaves using hyperspectral imaging. Postharvest Biol. Technol. 2013, 85, 8–17. [Google Scholar] [CrossRef]
  15. Wu, C.; Niu, Z.; Tang, Q.; Huang, W. Estimating chlorophyll content from hyperspectral vegetation indices: Modeling and validation. Agric. For. Meteorol. 2008, 148, 1230–1241. [Google Scholar] [CrossRef]
  16. Darvishzadeh, R.; Skidmore, A.; Schlerf, M.; Atzberger, C.; Corsi, F.; Cho, M. LAI and chlorophyll estimation for a heterogeneous grassland using hyperspectral measurements. ISPRS J. Photogramm. Remote Sens. 2008, 63, 409–426. [Google Scholar] [CrossRef]
  17. Elarab, M.; Ticlavilca, A.M.; Torres-Rua, A.F.; Maslova, I.; McKee, M. Estimating chlorophyll with thermal and broadband multispectral high resolution imagery from an unmanned aerial system using relevance vector machines for precision agriculture. Int. J. Appl. Earth Obs. Geoinform. 2015, 43, 32–42. [Google Scholar] [CrossRef]
  18. Carmona, F.; Rivas, R.; Fonnegra, D.C. Vegetation Index to estimate chlorophyll content from multispectral remote sensing data. Eur. J. Remote Sens. 2015, 48, 319–326. [Google Scholar] [CrossRef]
  19. Tewari, V.K.; Arudra, A.K.; Kumar, S.P.; Pandey, V.; Chandel, N.S. Estimation of plant nitrogen content using digital image processing. Agric. Eng. Int. CIGR J. 2013, 15, 78–86. [Google Scholar]
  20. Fernández-Pacheco, D.G.; Escarabajal-Henarejos, D.; Ruiz-Canales, A.; Conesa, J.; Molina-Martínez, J.M. A digital image-processing-based method for determining the crop coefficient of lettuce crops in the southeast of Spain. Biosyst. Eng. 2014, 117, 23–34. [Google Scholar] [CrossRef]
  21. Juan, W.; Changzhou, W.; Jinqiang, G.; Yongwen, L. A Method Based on Digital Image Analysis for Estimating Crop Canopy Parameters. In Proceedings of the 2011 Fourth International Conference on Intelligent Computation Technology and Automation, Shenzhen, China, 28–29 March 2011. [Google Scholar]
  22. Wang, Y.; Wang, D.; Zhang, G.; Wang, J. Estimating nitrogen status of rice using the image segmentation of G-R thresholding method. Field Crops Res. 2013, 149, 33–39. [Google Scholar] [CrossRef]
  23. Moghaddam, P.A.; Derafshi, M.H.; Shirzad, V. Estimation of single leaf chlorophyll content in sugar beet using machine vision. Turk. J. Agric. For. 2011, 35, 563–568. [Google Scholar]
  24. Ma, L.; Fang, J.; Chen, Y.; Gong, S. Color Analysis of Leaf Images of Deficiencies and Excess Nitrogen Content in Soybean Leaves. In Proceedings of the 2010 International Conference on E-Product E-Service and E-Entertainment, Henan, China, 7–9 November 2010. [Google Scholar]
  25. Lee, K.; Lee, B. Estimation of rice growth and nitrogen nutrition status using color digital camera image analysis. Eur. J. Agron. 2013, 48, 57–65. [Google Scholar] [CrossRef]
  26. Story, D.; Kacira, M.; Kubota, C.; Akoglu, A.; An, L. Lettuce calcium deficiency detection with machine vision computed plant features in controlled environments. Comput. Electron. Agric. 2010, 74, 238–243. [Google Scholar] [CrossRef]
  27. Wiwart, M.; Fordoński, G.; Żuk-Gołaszewska, K.; Suchowilska, E. Early diagnostics of macronutrient deficiencies in three legume species by color image analysis. Comput. Electron. Agric. 2009, 65, 125–132. [Google Scholar] [CrossRef]
  28. Odabas, M.S.; Simsek, H.; Lee, C.W.; İseri, İ. Multilayer Perceptron Neural Network Approach to Estimate Chlorophyll Concentration Index of Lettuce (Lactuca sativa L.). Commun. Soil Sci. Plant Anal. 2017, 48, 162–169. [Google Scholar] [CrossRef]
  29. Vesali, F.; Omid, M.; Kaleita, A.; Mobli, H. Development of an Android App to Estimate Chlorophyll Content of Corn Leaves Based on Contact Imaging. Comput. Electron. Agric. 2015, 116, 211–220. [Google Scholar] [CrossRef]
  30. Gu, D.-D.; Wang, W.-Z.; Hu, J.-D.; Zhang, X.-M.; Wang, J.-B.; Wang, B.-S. Nondestructive Determination of Total Chlorophyll Content in Maize Using Three-Wavelength Diffuse Reflectance. J. Appl. Spectrosc. 2016, 83, 541–547. [Google Scholar] [CrossRef]
  31. Hu, H.; Zhang, J.; Sun, X.; Zhang, X. Estimation of leaf chlorophyll content of rice using image color analysis. Can. J. Remote Sens. 2013, 39, 185–190. [Google Scholar] [CrossRef]
  32. Dey, A.K.; Sharma, M.; Meshram, M.R. An Analysis of Leaf Chlorophyll Measurement Method Using Chlorophyll Meter and Image Processing Technique. Procedia Comput. Sci. 2016, 85, 286–292. [Google Scholar] [CrossRef]
  33. Langseth, R.; Gaddam, V.R.; Stensland, H.K.; Griwodz, C.; Halvorsen, P. An Evaluation of Debayering Algorithms on GPU for Real-Time Panoramic Video Recording. In Proceedings of the 2014 IEEE International Symposium on Multimedia, Taichung, Taiwan, 10–12 December 2014. [Google Scholar]
  34. Heikkinen, V.; Lenz, R.; Jetsu, T.; Parkkinen, J.; Hauta-Kasari, M.; Jääskeläinen, T. Evaluation and unification of some methods for estimating reflectance spectra from RGB images. JOSA A 2008, 25, 2444–2458. [Google Scholar] [CrossRef] [PubMed]
  35. Chan, T.F.; Vese, L.A. Active contours without edges. IEEE Trans. Image Process. 2001, 10, 266–277. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Optical system for the image acquisition: (a) Schematic view, (b) optical system test.
Figure 1. Optical system for the image acquisition: (a) Schematic view, (b) optical system test.
Sensors 18 00650 g001
Figure 2. Image of the adaxial and abaxial leaf side in Bayer format.
Figure 2. Image of the adaxial and abaxial leaf side in Bayer format.
Sensors 18 00650 g002
Figure 3. Image in RGB format.
Figure 3. Image in RGB format.
Sensors 18 00650 g003
Figure 4. Leaf and background, separation process.
Figure 4. Leaf and background, separation process.
Sensors 18 00650 g004
Figure 5. Reflectance according to the chlorophyll content in Soil Plant Analysis Development (SPAD) values.
Figure 5. Reflectance according to the chlorophyll content in Soil Plant Analysis Development (SPAD) values.
Sensors 18 00650 g005
Figure 6. Transmittance according to the chlorophyll content in SPAD values.
Figure 6. Transmittance according to the chlorophyll content in SPAD values.
Sensors 18 00650 g006
Table 1. Regression models (simple linear regression).
Table 1. Regression models (simple linear regression).
Independent VariablesR2Standard Deviation (SPAD)NRMSE
Rr0.781.690.05
Rg0.802.600.07
Rb0.762.270.05
Tr0.941.190.28
Tg0.911.350.19
Tb0.921.300.26
Table 2. Multiple linear regression models (using two variables).
Table 2. Multiple linear regression models (using two variables).
Independent VariablesR2Standard Deviation (SPAD)NRMSE
Rr, Tr0.970.830.36
Rg, Tg0.960.940.43
Rb, Tb0.921.340.25
Table 3. Chlorophyll a, b and total content in Canavalia ensiforme, Azadirachta indica, and Lycopersicon esculentum leaves.
Table 3. Chlorophyll a, b and total content in Canavalia ensiforme, Azadirachta indica, and Lycopersicon esculentum leaves.
Canavalia ensiforme LeavesAzadirachta indica LeavesLycopersicon esculentum Leaves
Mean, stdR2NRMSEMean, stdR2NRMSEMean, stdR2NRMSE
Tr (%)0.16 ± 0.07 0.21 ± 0.09 0.29 ± 0.05
Chlorophyll a (µg/mL)26.07 ± 14.020.730.0116.47 ± 4.340.910.0221.97 ± 3.370.960.02
Chlorophyll b (µg/mL)11.80 ± 4.630.630.036.31 ± 1.730.980.308.12 ± 0.840.990.03
Chlorophyll b}total37.86 ± 18.490.660.1522.77 ± 6.060.940.0230.08 ± 3.530.970.03
Table 4. Processing speed of the proposed algorithm.
Table 4. Processing speed of the proposed algorithm.
CaseProcessing Time (ms)
F (Tr)150
F (Rg, Tr)188
Table 5. Vision-based approaches for chlorophyll content estimation.
Table 5. Vision-based approaches for chlorophyll content estimation.
ApproachAccuracy (R2)
H. Noh and Q. Zhang (2012), Whole area0.86
H. Noh and Q. Zhang (2012), Bright area0.87
H. Noh and Q. Zhang (2012), Corn area0.85
Tewari et al. (2013)0.94
Hao Hu et al. (2014), Green Value0.74
Hao Hu et al. (2014), Red Value0.75
Pagola et al. (2009), IpcaM40.92
Pagola et al. (2009), IpcaM20.92
Moghaddam et al. (2011), MLPN0.94
Moghaddam et al. (2011), R, B (regression)0.88
Kawashima et al. (1998), NORMALIZED ‘r’0.79
Kawashima et al. (1998), NORMALIZED ‘g’0.76
This work, F (Rr, Tr)0.97

Share and Cite

MDPI and ACS Style

Pérez-Patricio, M.; Camas-Anzueto, J.L.; Sanchez-Alegría, A.; Aguilar-González, A.; Gutiérrez-Miceli, F.; Escobar-Gómez, E.; Voisin, Y.; Rios-Rojas, C.; Grajales-Coutiño, R. Optical Method for Estimating the Chlorophyll Contents in Plant Leaves. Sensors 2018, 18, 650. https://doi.org/10.3390/s18020650

AMA Style

Pérez-Patricio M, Camas-Anzueto JL, Sanchez-Alegría A, Aguilar-González A, Gutiérrez-Miceli F, Escobar-Gómez E, Voisin Y, Rios-Rojas C, Grajales-Coutiño R. Optical Method for Estimating the Chlorophyll Contents in Plant Leaves. Sensors. 2018; 18(2):650. https://doi.org/10.3390/s18020650

Chicago/Turabian Style

Pérez-Patricio, Madaín, Jorge Luis Camas-Anzueto, Avisaí Sanchez-Alegría, Abiel Aguilar-González, Federico Gutiérrez-Miceli, Elías Escobar-Gómez, Yvon Voisin, Carlos Rios-Rojas, and Ruben Grajales-Coutiño. 2018. "Optical Method for Estimating the Chlorophyll Contents in Plant Leaves" Sensors 18, no. 2: 650. https://doi.org/10.3390/s18020650

APA Style

Pérez-Patricio, M., Camas-Anzueto, J. L., Sanchez-Alegría, A., Aguilar-González, A., Gutiérrez-Miceli, F., Escobar-Gómez, E., Voisin, Y., Rios-Rojas, C., & Grajales-Coutiño, R. (2018). Optical Method for Estimating the Chlorophyll Contents in Plant Leaves. Sensors, 18(2), 650. https://doi.org/10.3390/s18020650

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop