1. Introduction
The spectral reflectance of the crop canopy can provide information about crop growth. Using the spectral reflectance of the crop canopy can effectively realize dynamic monitoring of crop growth characteristics and nutritional status. Moreover, nitrogen content is often the most important information for indicating crop nutritional status [
1,
2,
3,
4,
5]. The non-destructive detection of crop nitrogen content is an important application of spectral information detection technology in the field of agriculture. By using this technology, combined with the diagnosis model of crop nitrogen nutrition, we can quickly obtain the nutrition information of crop growth inside a specific field size, which is conducive to the generation of the decision-making prescription of the nitrogen application rate and to the realization of variable rate fertilization. It is practically significant for the promotion of precision agriculture technology [
6,
7,
8,
9,
10].
According to the amount of detection wavebands, the current reflectance spectrometers can be divided into super-spectral, hyper-spectral, multi-spectral and specific-spectral types of spectrometers [
11,
12]. The first three kinds of reflectance spectrometers are usually more complicated and expensive, and are mainly used in aerial remote sensing or laboratory research to explore the quantitative relationship between crop growth and vegetation indices (Vis) with specific wavebands. Both super-spectral and hyper-spectral spectrometers have much higher spectral resolution than multi-spectral spectrometers and can produce continuous spectral measurement results. Super-spectral spectrometers have the highest spectral resolution, which can generate hundreds or even thousands of detection wavebands in the measurement wavelength range (generally 350–2500 nm), and hyper-spectral spectrometers usually include dozens or, at most, hundreds of detection wavebands, while the multi-spectral spectrometers only possess a few or a dozen. Super-spectral spectrometers mostly adopt the active scanning mode as the imaging mode and the Michelson double beam interference light splitting mode for light splitting. Correspondingly, hyper-spectral spectrometers usually use the optical imaging mode and optical grating light splitting mode, while multi-spectral spectrometers adopt optical imaging technology and linear gradient optical filter for light splitting [
13,
14,
15]. The specific-spectral spectrometers are for special purposes based on growth diagnosis models in which two detection wavebands in the visible and near-infrared wavelength range are involved. Moreover, the optical filters are usually used to ensure the bandwidth of the detection wavebands [
16,
17,
18,
19,
20,
21]. Due to the characteristics of portable carrying, low cost and high suitability, specific-spectral spectrometers are widely used in the real-time diagnosis of crop growth and can be divided into passive and active types according to the property of the light source.
The earlier versions (in the 1990s) of specific-spectral spectrometers (such as GreenSeeker by NTech Industries Inc., and N-sensor by Hydro Agri Gmbh) are all belong to the category of passive light source instruments. This type of spectrometer uses sunlight as the detection light source. The biggest limitation is that the measurement results are easily affected by the variations in intensity of sunlight, sun elevation angle, atmospheric conditions, etc. In particular, they cannot be used in low illuminance environments (such as cloudy days, early morning or the late afternoon on sunny days) [
22,
23]. Bajwa et al. used a portable dual-band reflectance spectrometer (Analytical Spectral Devices, developed by Boulder Co.) to measure the spectral reflectance of a rice canopy. They found that it was appropriate to use the device at 10:00 a.m.–15:00 p.m. in clear weather to decrease the influence of the change in the sun’s elevation angle on the measurement results [
1]. Peteinatos’s review also mentioned that, with weed monitoring sensors, the measured results were easily affected by the change of ambient light environments because of changes in the cloud layer, the atmospheric variances and sun elevation angle [
24]. Yang et al. carried out a reliability analysis of Cropsense (developed by Beijing Research Center for Information Technology in Agriculture, China) and SRS_NDVI (a commercial product of Decagon Co., USA). The experimental results showed that the suitable measurement condition was related to the reliability of the spectrometers. If the intensity of sunlight was greater than 5.2 klux, the influence of light intensity and sun elevation angle was not significant, and the reliability of the spectrometers decreased obviously under other conditions [
25].
In order to overcome the limitation in application of the spectrometers based on passive light sources, spectrometers based on active light sources began to appear in the late 1990s. This type of spectrometer is characterized by its own light source with specific wavebands. The light source is illuminated in a high frequency interval mode. The purpose is to make the frequency of the spectrometer response signal caused by the active light obviously higher than that caused by the ambient light. This makes the spectrometers equally easy to operate under all ambient lighting conditions.
There were already some commercial products of reflectance spectrometers based on active light sources since 2000, such as GreenSeeker (from Trimble Navigation) [
26], Crop Circle (from Holland Scientific Inc.) [
27], CropSpec (from Topcon Precision Agriculture Inc.) [
28] and N-sensor ALS serial products (from Yara Inc.) [
29]. The working principle of these reflectance spectrometers developed by different companies is basically similar, but they have their own application characteristics and suitable application environments because of differences in detection wavebands, the range of measurement heights, and the arrangement of the light sources (vertical or oblique) [
17,
21]. GreenSeeker adopts an Light-Emitting Diode (LED) light source with center wavelengths of 660 and 770 nm. The radiation surface of the light source is elliptical, and the effective measurement height range is 60–160 cm. The Crop Circle ACS470 adopts a multi-color LED light source whose central wavelengths can be adjusted by changing optical filters. The multi-channel detectors can sense the reflected light of the target simultaneously. The N-Sensors ALS uses a flashing hernia lamp as the light source, including the left and right detectors to measure the crop canopy on both sides of the vehicle. It can reach a large size of detection area because of very high measurement height. The central wavelengths of Cropspec are 735 and 808 nm, which is from a pulsed laser light source. The light source can be installed 2–4 m high with a Field of View (FOV) of 45–55° to form a large detection area. Generally, handheld spectrometers (such as Greenseeker and CropCircle) tend to be lightweight, with LEDs as light sources, and are mainly used to measure small areas (usually in the range of 0.1–0.5 m
2) [
30]. Reflectance spectrometers for large area detection (such as Cropspec and N-sensors ALS) usually use high-power light sources (modulated laser or hernia lamp for instance) with high measurement heights (over 2 m) to increase the detection area (some are even greater than 10 m
2).
Manufacturers of spectrometers were reluctant to provide technical material because of commercial benefits. At the same time, few researchers conducted research on the general scientific or engineering issues of such spectrometers. Spectral diagnosis technology depends on the quantitative relationship model between crop growth and specific VIs. In recent years, the increasing application of hyperspectral spectrometers had greatly promoted the establishment of this relationship model. Because of the particularity of the detection wavebands used in the relationship models which were established by different researchers, it is difficult to purchase a ready-made spectrometer to carry out the relevant research. This makes the scientific research outcomes difficult to apply in real production.
In this paper, we developed a reflectance spectrometer based on an active light source with central wavelengths of 730 and 810 nm and provided a method for designing this kind of spectrometer. The VIs in the form of ratio (such as Normalized Difference Vegetation Index (NDVI) and Ratio Vegetation Index (RVI)) can be measured independent of the measurement height and the ambient light environment. The performance tests were carried out in the field to verify the feasibility and reliability of the method.
2. Materials and Methods
In this section, starting from the problems with the conventional calculation method of VIs, a new calculation method for VIs based on an active light source is proposed (an active light source is a kind of artificial light source which is equipped as a part of the spectrometer, and usually produces modulated light in two or more special wavebands), and the measurement conditions that the calculation method should meet are analyzed. Finally, a reflectance spectrometer based on this new method is introduced.
2.1. Method for Calculating Vegetation Index
The Normalized Difference Vegetation Index (NDVI) and the Ratio Vegetation Index (RVI) are the most commonly used VIs in the form of ratios, and they reflect the difference in vegetation coverage and growth conditions well, being especially suitable for detecting vegetation with vigorous growth and high coverage. But the conventional calculation method of VIs has many inconveniences and difficulties in practical applications.
2.1.1. Conventional Method for Calculating Vegetation Index
The conventional calculation formula of NDVI and RVI are shown in Equations (1) and (2), respectively.
where
Rλ is the reflectance of the specific waveband whose central wavelength is λ,
λ1 and
λ2 are central wavelengths of near-infrared and visible light bands for measurement, respectively.
The measurement of the reflectance Rλ in Equations (1) and (2) requires a specific measurement height in the case of active light, which can be expressed as the ratio of canopy to whiteboard-reflected radiation in a specific waveband.
The formula for calculating reflectance can be described as Equation (3).
where
LWλ-H is the response value of the spectrometer to the whiteboard-reflected radiation of the waveband with central wavelength
λ at the measurement height
H,
LCλ-H is the response value of the spectrometer to the canopy-reflected radiation of the waveband with central wavelength
λ at the measurement height
H.
From Equations (1) and (2), it can be known that the calculation of VIs, such as NDVI and RVI, need to obtain the reflectance of the detection light wavebands. Equation (3) shows that the reflection value of the standard whiteboard at a specific height is required as the reference value for the calculation of reflectance. Therefore, the spectrometer is required to be fixed at a certain height in actual use, which makes the actual operation very inconvenient.
2.1.2. Method for Calculating Height-Independent Vegetation Indices
To express the calculation method of VIs in a new form, the standard whiteboard response ratio coefficient
CW is introduced into the equations and defined as Equation (4):
CW can be measured by experiments. Moreover, the NDVI and RVI can be rephrased as the sum of Equations (5) and (6) with C
W.
For instance, the corresponding derivation steps of Equation (5) can be expressed as follows:
If the value of CW keeps constant independent of the measuring height, the calculation method, which is expressed by Equations (5) and (6), is fundamentally different from that expressed by Equations (1) and (2). Equations (5) and (6) indicate that NDVI and RVI can be tested independent of the measurement height, which greatly facilitates the practical application of the spectrometers since there is no need to fix the measurement height.
The test results of the reflectance sensor (developed by METER Group, Inc. in Pullman, WA, USA) showed that the ratio of red to NIR spectral irradiance α can be used as a rough approximation at midday under clear sky conditions. But this conclusion was not extensively tested since any fluctuations in α may occur with changes in atmospheric conditions or across large variations in sun elevation angle [
31]. The meaning of ratio α is similar to
CW expressed in Equation (4), except that α is the radiation ratio for a natural light environment, while
CW is that of an active light environment.
CW has more stable characteristics than α because the value of
CW is independent of the natural light environment.
To keep the CW constant with the measurement height changes, two criteria need to be met for designing spectrometers. (i) The two different spectral measurement wavebands of active light source form the same or largely the same irradiation area. (ii) The two photosensitive detectors for different spectral reflectance need to be installed within close proximity to each other, which can ensure basically the same incident angle for each reflecting point to the two photosensitive detectors.
When the functional design of the spectrometer meets the above two criteria, the signal intensity detected by the photosensitive detectors is directly proportional to the intensity of the reflected radiation because of the same FOV. Moreover, the ratio of the radiation intensity of the two spectral wavebands emitted by the active light source is unchanged, which ensures that CW maintains a constant value.
2.2. Design Scheme of Reflection Spectrum Measuring Device
In this section, a design scheme for an active light source and optical detector system is proposed to meet the design criteria proposed in
Section 2.1, and the impacts of relevant design parameters on measurement results are analyzed theoretically.
2.2.1. Layout of Active Light Source
The active light source mainly consists of an array of LEDs with specific central wavelengths and a cylindrical lens. To make the active light source with two different detection wavebands produce a uniform illumination area and make the spectrometer meet the first item of the design criteria, the LEDs with different wavebands in the array are arranged interlaced with each other in the light source layout shown in
Figure 1. The central optical axis of the LED is inclined in a vertical, downward direction to expand the illumination area. The emitting light of the LEDs is reshaped by the cylindrical lens and forms an approximate rectangular illumination area.
2.2.2. Parameter Design Requirements for Detecting Light Path System
To detect the intensity of reflected light with different wavebands, two optical detection paths are designed. The distance between the central axes of the two optical detection paths is represented by
b, and the measurement height is expressed by
H. A schematic diagram of the optical detection paths is shown in
Figure 2.
Because the intensity of the reflected light entering the optical path system from each reflecting point is proportional to the cosine value of the incident angle
θ, the relative deviation of the reflected light intensity between the two optical detection paths
E can be expressed by Equation (7).
The position
x_max where
E reaches its maximum can be calculated with Equation (8) which is obtained from the derivative of Equation (7). Moreover, the maximum value of
E, expressed by
E_max, can be calculated with Equation (9):
where
. From Equation (9),
E can be controlled within the expected range by selecting appropriate values for parameters
b and
H. The design process of determining parameter
k is actually the process of meeting the design criterion (ii).
2.3. Reflectance Spectrometer
In this section, we will introduce a self-developed spectrometer which is a specific-spectral type of spectrometer based on an active light source with 730 and 810 nm as the central wavelengths. The spectrometer only extracts the response information caused by the active light source to calculate the VIs, and the response information from ambient light sources is effectively eliminated.
2.3.1. Working Principle of the Spectrometer
The working principle of the dual-band reflectance spectrometer is shown in
Figure 3. The LEDs are illuminated in a high frequency interval mode. In this mode of operation, LEDs work for every 15 ms, during which they are illuminated at a frequency of 2000 Hz, and then stop working for 85 ms. This working mode can ensure the LEDs hold a stable central wavelength, because continuous illumination results in the obvious increase in the LEDs’ junction temperature and a red shift in the emission spectrum of LEDs [
32,
33]. The rectangular illumination area is formed after being reshaped by the cylindrical lens. The photocell is used as a photosensitive detector. The reflected light, including both the active light and the sunlight, is focused by convex lens and filtered by optical filters to form narrow band light with a specific bandwidth and central wavelengths, which is finally received by the photocell. So, the response signal of the photocell is a superposition result of active light and sunlight response signals; the amplitude of the photocell response signal is the sum of the amplitude of the response signal caused by active light and sunlight independently.
Figure 4 shows the two-minute-long sampling results of the photocell response signal when the spectrometer was fixed in a spot of a wheat field in a sunny day (there were clouds and occasionally breezes during the test.). The amplitude of the response signal caused by active light was very small compared with that caused by the sunlight. So, the amplitude changing trend of the photocell response signal was determined by the response signal caused by sunlight (or ambient light) which changed relatively slowly and was simultaneously superimposed together with the square wave signal at 2000 Hz, which was caused by the active light. Due to the obvious difference in frequency between the two types of response signals, the response signal caused by the active light can be easily extracted from the photocell response signal through high pass filtering and can maintain its own original amplitude, so as to eliminate the influence of sunlight on the measurement results.
The whole signal processing flow is expressed in
Figure 5. As shown in the processing procedure shown in
Figure 5, the response signal of the photocell, which is illustrated in
Figure 4, was firstly filtered by the high-pass filter to separate the response signal caused by active light. Then, the filtered signal, which was a signal with a high frequency square wave, was amplified again to extract the amplitude of the signal. The amplitude of the active light response signal after it was amplified for a second time is shown in
Figure 6. The method for how to extract the amplitude of the active light can be found in [
33]. Both the average and standard deviation of the amplitude of the response signal, shown in
Figure 6, were 0.545 and 0.013, respectively. It was indicated that the amplitude of the response signal caused by active light was stable, though the sunlight intensity kept changing during the test, and the response signal caused by sunlight was efficiently removed from the response of the photocell after being filtered. It should be pointed out that the fluctuation in response amplitude may be caused by circuit interference or disturbance of the crop canopy caused by the wind.
2.3.2. Main Working Elements and Performance Parameters of the Spectrometer
The self-developed reflectance spectrometer according to the above design principles is shown as
Figure 7. The main working elements of the spectrometer include LEDs with specific central wavelengths, a cylindrical lens, a convex lens, an optical filter and a photocell. The array of LEDs and the cylindrical lens constitute the active light source; the convex lens and optical filters are the main components of the optical detection paths.
(1) LEDs and Optical Filters
The central wavelengths of the reflectance spectrometer are 730 and 810 nm, and the full width at half maximum (FWHM) is less than 20nm. The optical filters (by Beijing Xintian Borui Photoelectric Technology Co., Ltd., Beijing, China) are covered on the surface of the photocells to further reduce the bandwidth of the reflected active light since the original FWHM of the LEDs (XL302440NRC-730 and XL302430HIRC-810, by Nanning Lvxing Lighting Electronics Co., Ltd., Nanning, China) cannot meet the design requirement. The light transmittance of two different optical filters is shown in
Figure 8, and the relative light intensity distribution of LED emission light before and after passing through the optical filters is shown in
Figure 9 (the testing data of the LEDs and filters are all provided by the manufacturers). The final central wavelengths of the 730 and 810 nm wavebands are 730.3 and 809.5 nm, respectively, and the FWHM values are 10.8 and 19.5 nm.
(2) Cylindrical Lens and Convex Lens
The size of the single cylindrical lens is 80 mm × 20 mm, and its focal length is 30 mm. The LED arrays are installed where the focal line of the cylindrical lens is located (26.14 mm away from the plane of the cylindrical lens). With the action of the cylindrical lens, the active light source forms a more concentrated irradiation area, which improves the detectability of the reflected light.
The focal length of the convex lens is 15 mm, and the distance between the photocell surface and the plane of the convex lens is 10.8 mm where the focus plane of the convex lens is located. The focusing function of the convex lens ensures that all the light reflected into the spectrometer falls into the detection area of the photocell. The two lenses are made of BK7 glass, which is the product of the Shotto Glass Company in Germany.
(3) Photocells
A type of silicon photocell (2CR40, by Shanghai Qinyue Electronic Technology Co., Ltd., Shanghai, China) with high sensitivity and wide spectral wavebands is used as the detector for reflected light radiance, its detection area is 10 mm × 20 mm, and its sensitivity curve is shown in
Figure 10, which shows its higher sensitivity at wavelengths of 730 and 810 nm (sensitivity data are provided by the manufacturer). The typical response time of the photocell is 70 μs which can match the requirement for detecting signals with a high frequency of 2000 Hz.
(4) Working Parameters of the Spectrometer
The working frequency of the LEDs set in the spectrometer is 2000 Hz. After being reshaped by the cylindrical lens, the irradiation of the active light forms an approximate rectangular area whose geometric dimension meets the following relationship: width = 0.25× height, length = height. The distance between the central axes of the two optical detection paths
b = 22 mm (see
Figure 2), and the designed effective range of measurement height of the spectrometer is 60–120 cm. When the measurement height changes, the theoretical maximum relative intensity deviation of the incident light entering the two optical detection paths is shown in
Figure 11 (the theoretical maximum relative intensity deviation can be calculated with Equation (9) in
Section 2.2.2). It can be seen from
Figure 11 that the maximum relative deviation of the incident light intensity is less than 2% within the effective range of the measurement height.
2.4. Calibration of the Spectrometer
To evaluate the stability and accuracy of the spectrometer, some standard grayscale plates (Hefei Xingyue Luminous Technology Applications Institute, Hefei, China) were used as the measurement objects to carry out the calibration tests. The reflectance measured by the spectrometer was compared with the nominal reflectance of the grayscale plates to evaluate the accuracy of the spectrometer. On other hand, the standard whiteboard response ratio coefficient CW was tested in different measurement height in indoor and outdoor environments to value the stability of the spectrometer.
2.4.1. Reflectance Measurement to Standard Grayscale Plates
Six standard grayscale plates were used as the test objects, whose nominal reflectances are 5%, 10%, 20%, 40%, 60% and 75%, respectively. Moreover, the standard whiteboard (made by the same producer as the grayscale board) with a nominal reflectance of 99%, was treated as the reference object. The size of all the grayscale plates and the whiteboard is 60 mm × 60 cm. The spectrometer was placed 70 cm above the grayscale plate during the tests. The experiments were carried out in an indoor laboratory environment under fluorescent lamps. To eliminate the reflection interference outside the grayscale plate as much as possible, a black cotton cloth was placed under the grayscale board. The reflectance of the black cotton cloth was less than 2%.
Once the response of the spectrometer to each grayscale plate and whiteboard were measured, the reflectance of each grayscale plate was calculated with Equation (3). The reflectance of each grayscale plate was measured three times during tests.
Figure 12 shows the linear regression results based on the measured reflectance (both at 730 and 810 nm) and the nominal reflectance of the grayscale plates. The linear regression results showed that the designed spectrometer had stable and accurate measurement results.
2.4.2. Determination of the Standard Whiteboard Response Ratio Coefficient CW
Ten measurement heights, at an interval of 5 cm, between 65 and 110 cm were set during indoor tests. The test process at each measurement height lasted for 30 s, during which the data update frequency was 10 Hz. The average value of those recorded data in 30 seconds was finally taken as the whiteboard response value of the spectrometer to calculate the
CW.
Figure 13 shows the reflection response of the spectrometer to the standard whiteboard and the
CW at different measurement heights.
Three days in different weather conditions (sunny, cloudy and overcast) were selected to carry out the outdoor tests. Four testing time points of 10:00, 12:00, 14:00 and 16:00 were selected to start the repeated tests each day, during which the
CW were measured at five heights (at an interval of 10 cm, from 60 to 100 cm). Each test (including the five measurement heights) lasted about 15 min. The minimum and maximum illuminance of ambient light were recorded by two digital illuminance meters (AS823, by Smart Sensor
®, Shenzhen, China) during each test. The test value of
CW at different measurement heights at each testing time point are shown as a boxplot in
Figure 14. The average illuminance at each testing time point outdoors is displayed as the bars in
Figure 14 for reference. The illuminance conditions during each test outdoors are described in
Table 1.
The mean value of CW in an indoor environment was 1.206, the standard deviation was 0.0034, and the coefficient of variation (CV) of CW was 0.26%. The corresponding results of CW outdoors under different weather conditions were 1.213, 0.0058 and 0.48%, respectively. It indicated that CW kept a high stability when the measurement height changed and kept constant whether indoors or outdoors. The results of the outdoor tests were highly consistent with those of the indoor tests. This provided a guarantee for the reliable measurement of field vegetation indices.
4. Conclusion
Some researches had shown that the spectral irradiance ratio of different wavebands can be used as a rough approximation during midday under clear sky conditions [
31]. Using this characteristic, the structure of the reflectance spectrometer can be simplified. However, due to the changes in the atmospheric conditions and the solar altitude angle, the ratio may have a large fluctuation, which results in poor reliability in practical applications. This is why the performance of the spectrometer based on a passive light source can be affected by the ambient light environment.
This paper provided a method for using an active light source to keep the ratio of reflected irradiance from a standard whiteboard between expected wavebands constant. We introduced the standard whiteboard response ratio coefficient CW to embody this characteristic. The VIs can be calculated independent of measurement height and ambient light environment based on CW. This effectively overcomes the defects of passive light source spectrometers.
Following the design criteria of the active light source proposed in this paper, the CW of the self-developed spectrometer showed a good characteristic of keeping constant at different measurement heights and under indoor or outdoor light environments. The field tests results showed that the measured value of VIs under different measurement heights maintained a high stability, and that the better the consistency of the measured objects inside the effective radiation area of the active light source, the higher the stability of the VI results. The experimental results verified the feasibility and validity of the proposed measurement method and the design criteria.
However, the current measurement method cannot eliminate the interference caused by soil-reflected irradiance, and the spectrometer is not suitable for large areas of bare soil because the VIs are sensitive to the optical properties of the soil background. For a given amount of vegetation, darker soil substrates result in higher vegetation index values [
35]. This is a problem for us and further research is necessary, both by us and by the manufacturers of commercial spectrometers, in order to solve it.