Next Article in Journal
The Role of Emission Sources and Atmospheric Sink in the Seasonal Cycle of CH4 and δ13-CH4: Analysis Based on the Atmospheric Chemistry Transport Model TM5
Next Article in Special Issue
Winter Potato Water Footprint Response to Climate Change in Egypt
Previous Article in Journal
High Level of Ammonium Nitrogen Increases Net Ecosystem Productivity in a Quercus liaotungensis Forest in Northern China
Previous Article in Special Issue
Quantifying Agricultural Drought Severity for Spring Wheat Based on Response of Leaf Photosynthetic Features to Progressive Soil Drying
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Response of Population Canopy Color Gradation Skewed Distribution Parameters of the RGB Model to Micrometeorology Environment in Begonia Fimbristipula Hance

1
Jiangsu Meteorological Bureau, Nanjing 210008, China
2
Longyan Company of Fujian Provincial Tobacco Corporation, Longyan 364000, China
3
QinGengRen Modern Agricultural Science and Technology Development (Huai’an) Co., Ltd., Huai’an 223001, China
4
Jiangsu Key Laboratory for the Research and Utilization of Plant Resources, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China
5
Department of Chemistry, Faculty of Science, Carleton University, Ottawa, ON K1S 5B6, Canada
6
Jiangsu Provincial Key Laboratory of Agricultural Meteorology, College of Applied Meteorology, Nanjing University of Information Science and Technology, Nanjing 210044, China
*
Author to whom correspondence should be addressed.
Submission received: 26 April 2022 / Revised: 18 May 2022 / Accepted: 19 May 2022 / Published: 30 May 2022
(This article belongs to the Special Issue Agrometeorology)

Abstract

:
The high quality and efficient production of greenhouse vegetation depend on micrometeorology environmental adjusting such as system warming and illumination supplement. In order to improve the quantity, quality, and efficiency of greenhouse vegetation, it is necessary to figure out the relationship between the crop growth conditions and environmental meteorological factors, which could give constructive suggestions for precise control of the greenhouse environment and reduce the running costs. The parameters from the color information of the plant canopy reflect the internal physiological conditions, thus, the RGB model has been widely used in the color analysis of digital pictures of leaves. We take photographs of Begonia Fimbristipula Hance (BFH) growing in the greenhouse at a fixed time every day and measure the meteorological factors. The results showed that the color scale for the single leaf, single plant, and the populated canopy of the BFH photographs all have skewed cumulative distribution histograms. The color gradation skewness-distribution (CGSD) parameters of the RGB model were increased from 4 to 20 after the skewness analysis, which greatly expanded the canopy leaf color information and could simultaneously describe the depth and distribution characteristics of the canopy color. The 20 CGSD parameters were sensitive to the micrometeorology factors, especially to the radiation and temperature accumulation. The multiple regression models of mean, median, mode, and kurtosis parameters to microclimate factors were established, and the spatial models of skewness parameters were optimized. The models can well explain the response of canopy color to microclimate factors and can be used to monitor the variation of plant canopy color under different micrometeorology.

1. Introduction

Meteorological conditions are crucial factors that determine crop growth and leaf development [1,2]. Facility agriculture is an efficient crop production through the precise regulating of microclimate within the greenhouse [3,4,5]. Clarifying the relationships between plant physiological and ecological indicators and meteorological conditions are the basic method for quantitative assessment of agro-meteorological conditions, assessing the agro-meteorological disaster losses, predicting the yield, and applying the intelligent micrometeorology control of agricultural production. Thus, it is the right path to find the determining parameters that have a close relationship with meteorological factors and reflect plant growth status sensitively.
Digital image technology is well-developed and high-resolution camera equipment is widely used, with digital imaging technology having the obvious advantages of high resolution and low cost in the research of plant phenotypes [6,7]. A digital color image is not only convenient to obtain [8], but also contains abundant information about plant morphology, structure, and color gradation [9,10,11], which can reflect the internal status of plants. The RGB color model is the popularly used color representation for digital images [12]. The parameters of the traditional RGB model are mainly the mean values of each color channel and their combined values, which can approximately describe the color depth of the leaves and estimate the chlorophyll content of the leaves [13,14,15,16,17,18]. Moreover, the RGB model reflects the soil moisture and nutrient levels [19,20,21,22]. However, the mean values of image parameters are not comprehensive enough in the description of color distribution [23], as a result, the RGB model is generally used in the studies of water, fertility, and disease which may drastically change the leaf color [8,24,25,26]. The RGB model is also an effective method to reflect the relationship between the image information and the variation of environmental factors.
The quantitative analysis of mega-high quality digital photographs could explain the genotype-by-environment interactions in plant phenotypes [27]. The variation of meteorological factors affects the physiological and biochemical characters of plants slightly, such as influencing the openness of the stomata [28], changing the arrangement of the chloroplast [29], and increasing or decreasing the photosynthesis rate of the leaves [30,31]. All these affections change the plant’s color under a visible spectrum [32]. Therefore, the image color information parameters are used as representative references to reflect the meteorological environment of the growing plant. Skewed analysis expanded the RGB model parameters from the mean to the mean, median, mode, skewness, and kurtosis, with the number of parameters from 4 to 20 [33]. The color gradation skewness-distribution (CGSD) parameters describe leaf color information more accurately and comprehensively than the traditional mean parameters [33], including whether the color gradation distribution of a single plant varies the same as the leaf, and the field scale determines its sensitivity in response to the environmental meteorological factors. The CGSD for all scales varying consistently helps to construct an effective model between the external meteorological factors and morphological parameters of the plant. The quantitative evaluation of agro-meteorological conditions for a plant can be achieved by combining canopy color information indexes corresponding to different morphological parameters of plants, to provide information that supports precise management of plant production.
Begonia Fimbristipula Hance (BFH) is a vegetation type widely used for medicine and food. BFH has abundant anthocyanins, polyphenols, and flavonoids and a strong antioxidant capacity [34,35,36]. BFH is cultivated in facilities that have a warm and humid environment to ensure annual cultivation [37]. After planting BFH in a glass greenhouse and collecting the meteorological factors and canopy photographs, the relationship between changes in canopy color and the meteorological factors in the facility environment was analyzed. The purpose of this research is: (1) to quantify the mathematical relationship between CGSD parameters and meteorological factors, and increase the accuracy of quantifying the changes in population canopy color with meteorological changes; (2) to construct population canopy color–meteorological response models by using the relationship of CGSD parameters and meteorological factors, including temperature, humidity, and solar radiation, and; (3) to provide new methods from the RGB model to monitor the patterns of canopy color variations due to meteorological environment changes. The results may construct the crop growth model based on the CGSD parameters of the RGB model, and provide a new approach for intelligent control of facility agricultural production.

2. Materials and Methods

2.1. Plant Material and Growth Conditions

The glass greenhouse for the experiment is in Guli Modern Agriculture Demonstration Park, Jiangning District, Nanjing, Jiangsu, China (24°57′ N, 116°30′ E). The BFH were potted by the rooting of cuttings and planted with a total nutrient matrix which was in the proportion of N: P2O5: K2O in a 100 g matrix at 20:10:20, trace elements included Cu, Fe, Zn, Mn, B, and Mo ≥ 1%, and pH 5.85. Moreover, the matrix contained an active accelerator and sustained release agent of poly-fertilizer. The BFH experiment had 2 replicates, and each replicate contained 24 potted plants. The first growing period of BFH was from 1 December 2018 to 17 February 2019, and after cutting off the 10 cm top length on 18 February, BFH continued growing in the second growth period from 19 February 2019 to 11 April 2019.

2.2. Meteorological Data Acquisition

Meteorological data were obtained from the agricultural meteorological observation station (DZZ4, Jiangsu Provincial Radio Science Research Institute Co., LTD., Wuxi, China) in the facility. The observed meteorological elements included daily meteorological elements, such as the daily mean temperature (Tdm), daily mean relative humidity (RHdm), daily mean ground temperature (TGdm), daily mean soil temperature at 10 cm (TSdm−10c), daily mean vapor pressure (VPdm), daily mean dew point temperature (TDdm), daily total global radiation (GRdt), and daily total photosynthetically active radiation (PARdt). The accumulated temperature (AT), accumulated global radiation (AGR), and accumulated photosynthetically active radiation (APAR) were calculated from Tdm, GRdt, and PARdt.
AT = T dm
AGR = GR dt
APAR = PAR dt

2.3. Canopy Image Collection

Canopy color images of BFH were collected in the facility by a monitoring camera (DH-SD-65F630U-HN-Q, Zhejiang Dahua Technology Co., Ltd., Hangzhou, China) with an image resolution of 1920 × 1080, and the installation height was 280 cm (Figure 1). The camera with a vertical lens was set with a fixed focal length shooting and automatic white balance and adopted a fixed time shooting mode to take a photograph at set times of 9:03 a.m. every day. The 117 images from 3 groups without direct sunlight were selected for analysis, including 102 images in the first growth period (1 December 2018 to 17 February 2019) of 2 replicates with 51 pictures each, and 15 images in the second growth period (19 February 2019 to 11 April 2019).

2.4. Cutting and Denoising of the Image

Adobe Photoshop V8.0 software (San Jose, CA, USA) and MATLAB2016R software (referred to as MATLAB [Math Works, Natick, MA, USA]) were primarily used to cut and denoise the BFH original image [38,39] (Figure 2a).
  • Adobe Photoshop CS software (San Jose, CA, USA) was primarily used to intercept the range of 600 × 600 in the lower-left corner of the image, and the processed image was saved in a JPG image format (Figure 2b).
  • The rgb2hsv function of MATLAB was used to convert RGB images into HSV images. Double cycle operation was used to set the H-value of the image background to 0 (i.e., black), while the H value of the plant remains unchanged. The hsv2rgb function was used to convert the processed HSV image into an RGB image (Figure 2c).
  • The double cycle operation of MATLAB was used again to filter the color threshold value of the image processed in the previous step. The color opacity of the black part of the image was adjusted to 0 (that is, completely transparent), and the color image of the target leaf or canopy with the transparent background was saved as a PNG image mode (Figure 2d).

2.5. Information Collection of the RGB Image

Transformation of double-precision arrays of the image. MATLAB was used to extract and analyze BFH RGB images. After reading color images by using the imread function, image (:, :, 1), image (:, :, 2), image (:, :, 3), and the rgb2gray function were respectively used to read every pixel gradation of red, green and blue channels, as well as gray-level images. The full circulation algorithm was used to retrieve and record the non-black pixel index codes of these pixels, which were combined with an array of leaf color gradation. Then the double function was used to transform it into double precision arrays again.
Establishment of the canopy color gradation skewness-distribution (CGSD) parameters table. The mean, median, mode, std, skewness and kurtosis functions were used to acquire the mean, median, mode, standard Deviation, skewness, and kurtosis of the double-precision arrays of red, green, and blue channels, as well as gray-level images [33]. The CGSD parameters were obtained, including RMean, RMedian, RMode, RSkewness, RKurtosis, GMean, GMedian, GMode, GSkewness, GKurtosis, BMean, BMedian, BMode, BSkewness, BKurtosis, YMean, YMedian, YMode, YSkewness, and YKurtosis (Table 1). Finally, the CGSD parameters tables of color gradation distribution of BFH RGB images were formed.
Array distribution normality testing. The lillietest and jbtest functions were used to conduct the Lilliefors and Jarque-Bera tests of normal distribution for the color gradation distribution of red, green, and blue channels, as well as gray-level images of BFH RGB images.

2.6. Prediction Model Construction and Goodness-of-Fit Detection

Correlation Analysis of 20 CGSD Parameters to Microclimate Factors. The Cor package of R was used to analyze the relationship between 20 CGSD parameters of population canopy in RGB images of BFH and the corresponding microclimate factors, including daily meteorological factors (Tdm, RHdm, TGdm, TSdm−10c, VPdm, TDdm, GRdt, and PARdt) and photothermal accumulation factors (AT, AGR, and APAR), which were collected for significant examination with a double tail inspection.
Linear prediction models establishment. By using R, the linear prediction models (Y1-Y20) of 20 CGSD parameters of population canopy in RGB images of BFH were established by a regression approach based on the least-square method, with meteorological factors as the independent variables. The probability of F-to-enter of the models was set to 0.050 or less, while the probability of F-to-remove was set to 0.100 or more. Decision coefficient optimization, significance test of regression model and regression coefficient, and collinearity diagnosis of independent variables of the regression model were conducted for the alternative regression models in turn [40], and the optimal regression models were finally determined.
Spatial polynomial model establishment. Two photothermal accumulation factors with the correlation coefficient of highest absolute value with CGSD parameters were selected as the independent variable, and the spatial polynomial models (Z1–Z4) of skewness of red, green, and blue channels, as well as gray-level images, were established by using polynomial fitting in the Curve Fitting Tool of R.

3. Results

3.1. Skew Analysis of the Distribution of Leaf Color Gradation of the RGB Images

After analyzing the cumulative distribution of color gradation of the 117 selected images of BFH, we found that all the red, green, and blue channels, and the gray-level for single leaf, single plant, and population canopy images showed a left-skewed distribution (Figure 3). Lilliefors and Jarque-Bera normality tests also indicated that color gradation data did not distribute normally (Supplementary Table S1). Twenty color gradient skewness-distribution (CGSD) parameters were obtained including the mean, median, mode, skewness, and kurtosis for the red, green, and blue channels, as well as the gray-level image, for each canopy image, respectively. These parameters can describe not only the depth of canopy color but also its distribution.
The main distribution of color gradation is positive skewness for RGB channels and the gray-level image near the value of 0. This increment of positive skewness as the canopy develops reflected the leaf color transit from light green to dark green (Figure 4).

3.2. Correlation Analysis Microclimate Factors in Glasshouse and Population Canopy CGSD Parameters

The relationship between 20 CGSD parameters and the microclimate factors including daily meteorological factors (Tdm, RHdm, TGdm, TSdm−10c, VPdm, TDdm, GRdt, and PARdt) and accumulative factors (AT, AGR, and APAR) are shown in Table 2. Table 2 indicates that temperature, humidity, ground temperature, soil temperature at 10 cm, and solar radiation were significantly correlated with the CGSD parameters. The accumulated temperature and solar radiation (AT, AGR, and APAR) had an extremely significant correlation with the 20 CGSD parameters, with the correlation coefficients higher than 0.8. The parameters of mean, median, and mode, which represented color depth, showed a significant positive correlation with the meteorological factors, while the parameters of skewness and kurtosis representing color uniformity showed an obviously negative correlation.

3.3. Multiple Linear Relationships of the Microclimate Factors and Population Canopy CGSD Parameters

Multiple linear relationships between population canopy CGSD parameters and microclimate factors were established by the stepwise regression method using the ordinary least square method (OLS) (Table 3). Table 3 showed that all the multivariate determination coefficients (R2) of the population canopy color-meteorological response models equation were greater than 0.7, indicating that the equation well explained the response of BFH population canopy CGSD parameters to microclimate factors. The R2 of the red and green channels, as well as the gray-level image, was 0.692–0.954, which was better than those of the blue channel with 0.674–0.791. The model between the mean of green channel and accumulated temperature, daily mean relative humidity, and the daily total photosynthetically active radiation had the highest R2 of 0.954. Moreover, all the models were significant (p < 0.01).
All the models were verified by between-group and outside-group samples (Table 4 and Figure 5, Figure 6 and Figure 7), which indicated that the population canopy color fitting models performed well in predicting the mean, median, mode, and kurtosis. The accuracy of prediction of the skewness was about 80%, while other CGSD parameters exceeded 90%, with the mean of the blue channel, as well as the gray-level image exceeding 95% (Table 4) in the same modeling group samples. Similarly, the accuracy of prediction of the mean, median, and mode of RGB channels and the gray-level image was 90%, the kurtosis was about 75%, and the skewness was 53–58% in the between-group prediction. The accuracy of prediction of other CGSD parameters was 80% with the exception of the skewness and kurtosis of the red channel, the mode, and skewness of the green channel, and the median of the gray-level image in the outside-group. Moreover, the accuracy of prediction of the kurtosis of the green and blue channels, and the gray-level image exceeded 85%.

3.4. Optimization of Canopy Color Skewness-Meteorological Models

Although the skewness of the red and green channels, and the gray-level image performed well, the between-group and outside-group prediction accuracy were generally lower than other CGSD parameters. The skewness was non-linear to microclimate factors [33]. Therefore, we fitted the polynomial Z1-Z4 that incorporates spatial bidirectional patterns for the RGB channels, and the gray-level image, by taking two accumulation factors, i.e., the accumulated temperature (x1) and the accumulated photosynthetically active radiation (x5) as the independent variables, which had the highest absolute value of the correlation coefficient.
Z1 = 0.3629 + 0.6382 × 10−3 x1 + 2.626 × 10−3 x5 + 5.312 × 10−6 x12 − 1.072 × 10−5 x1 x5
Z2 = 0.2242 + 0.5163 × 10−3 x1 + 3.11 × 10−3 x5 + 3.08 × 10−6 x12 − 7.463 × 10−6 x1 x5
Z3 = 0.3215 + 0.6477 × 10−3 x1 + 2.673 × 10−3 x5 + 3.93 × 10−6 x12 − 8.614 × 10−6 x1 x5
Z4 = 0.2786 + 0.5116 × 10−3 x1 + 3.114 × 10−3 x5 + 4.244 × 10−6 x12 − 9.285 × 10−6 x1 x5
After comparing the multiple stepwise regression and spatial bidirectional models using the skewness, the (R2) of the spatial bidirectional models was higher than the multiple stepwise regression models with lower RMSE. This indicated that the nonlinear models had better fitting results, or the skewness of the population canopy image of BFH has a close relationship to the microclimate in the glass greenhouse (Table 5).
The spatial bidirectional models for the skewness have 5.7 percentage points higher accuracy in prediction, and the number of outliers significantly decreased. Compared to the regression model, the average accuracy of the spatial bidirectional models was improved by 7.63, 5.75, and 6.25 percentage points in the modeling group, the between-group, and the outside-group, respectively. The RGB channels and gray-level image were all improved by 6.90, 10.87, 3.47, and 1.72 percentage points, respectively (Table 6).

4. Discussion

The high quality and efficient production of greenhouse plants depend on the environmental regulation focusing on the plants’ growth in the facility. The changing physiological parameters result from the variation of temperature, radiation, and humidity rapidly and instantaneously [41,42,43,44]. Ineffective control of environmental factors will restrict the improvement of glass greenhouse yield and crop quality [45]. Therefore, coupling the relationship between the plant growth state indicated by color parameters and controlling the meteorological factors will help to promote the high quality and efficient production in facility plants, realize the precise control of the facility environment, and decrease the costs.
Digital color images’ information can reflect the growing state of plants [11], which provides a convenient way of plant growth monitoring. Research on conducting qualitative and quantitative descriptions of phenotypic traits for plant appearance by using digital imaging technology has increased in smart agriculture [6]. The RGB color model has been widely used to process digital color images to study chlorophyll content and its related plant nutritional status and stress response [24,46,47]. The RGB model and leaf color parameters can describe changes in tobacco leaf color depth and homogeneity [33]. Now, we found that images of the RGB channels and the gray-level images all showed a skewed distribution at the single leaf, plant, and population canopy images scale in BFH.
Environmental changes lead to changes in the color of plants’ leaves [21,26,48,49,50]. Our work indicated that temperature, humidity, ground temperature, soil temperature at 10 cm, and solar radiation in glass greenhouses were significantly correlated with multiple CGSD parameters of canopy images in BFH. Moreover, the accumulative value of temperature and solar radiation (AT, AGR, and APAR) had an extremely significant correlation with 20 CGSD parameters, with the correlation coefficients beyond 0.8. The mean, median, and mode parameters representing color depth have a positive relationship with the accumulative parameters, while the skewness and kurtosis parameters representing the color distribution showed a significantly negative correlation with the accumulative factors. The deep canopy color of BFH images has higher accumulative factors depending on the analysis of the peak changes of the cumulative distribution diagram of canopy color and provides a new approach to quantitative plant growing conditions.
The multi-dimensionality of CGSD parameters enables them to have a basis for building relationships with complex meteorological factors. Canopy color–meteorological responding models have high prediction accuracy after calculating the image mean, median, mode, and kurtosis parameters and skewness distribution using the spatial bidirectional models. The change of accumulative factors is generally linear, while skewness is calculated from a higher-order function. Although the skewness was significantly related to microclimate factors, it cannot be simply described and fitted by linear fitting [33]. Compared to the multiple stepwise regression models, the spatial bidirectional models had better fitness and higher accuracy.
The digital color image has been applied to automatic observation in agrometeorology [51,52], and is proposed as a new technical approach that can be combined with crop growth models [53], which agriculture meteorological services will automatically and intelligently apply in future services [54,55]. The relationship between the CGSD parameters of the population canopy and the environment meteorological factors makes it possible to use the canopy color information to construct a plant growth model from digital images, which can provide new methods for quantitatively evaluating agricultural meteorological conditions, assessing agricultural meteorological disaster loss, predicting yield and quality, and supporting precise management of crop production.

5. Conclusions

Analyzing the canopy images’ color histograms skewed distribution of BFH in the glasshouse, we found that 20 population canopy CGSD parameters were sensitive to the microclimate factors, especially the radiation and temperature accumulation factors. The canopy color–meteorological responding models showed that the CGSD parameters have a close relationship with microclimate factors, which expands the potential application of the RGB model in monitoring the variations of plant or leaf color. The CGSD parameters can describe the canopy color information comprehensively and accurately and are considered a new indicator for quantitative evaluation of plant growth conditions. These parameters help construct a plant growth condition model and provide a new method for agrometeorological assessment and prediction, and, ultimately, support the precise management of crop production.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/atmos13060890/s1, Table S1: Normal test of color distribution of a BFH population canopy.

Author Contributions

P.Z. and Z.C. are the Co-first authors. The work presented here was carried out in collaboration among all authors. P.Z., Z.C., F.W., C.L. and X.J. designed the study. P.Z. and R.W. conducted the research. Z.A., T.B., and X.X. processed the data, P.Z., Z.C., F.W. and X.X. analyzed the data result. P.Z. wrote the entire manuscript. C.L. edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

We are grateful for the grants from the “333 project” research project for the high-level talent of Jiangsu Province (Grant number BRA2019348) and the National Natural Science Foundation of China (Grant number 51309132).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All datasets generated and analyzed for this study are included in the article/Supplementary Materials.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Zhou, H.L.; Zhou, G.S.; He, Q.J.; Zhou, L.; Ji, Y.H.; Zhou, M.Z. Environmental explanation of maize specific leaf area under varying water stress regimes. Environ. Exp. Bot. 2019, 171, 103932. [Google Scholar] [CrossRef]
  2. Nawaz, R.; Abbasi, N.A.; Hafiz, I.A.; Khalid, A. Impact of climate variables on growth and development of Kinnow fruit (Citrus nobilis Lour x Citrus deliciosa Tenora) grown at different ecological zones under climate change scenario. Sci. Hortic. 2019, 260, 108868. [Google Scholar] [CrossRef]
  3. Kempkes, F.; De Zwart, H.; Muñoz, P.; Montero, J.; Baptista, F.; Giuffrida, F.; Gilli, C.; Stępowska, A.; Stanghellini, C. Heating and dehumidification in production greenhouses at northern latitudes: Energy use. Acta Hortic. 2017, 445–452. [Google Scholar] [CrossRef]
  4. Amani, M.; Foroushani, S.; Sultan, M.; Bahrami, M. Comprehensive review on dehumidification strategies for agricultural greenhouse applications. Appl. Therm. Eng. 2020, 181, 115979. [Google Scholar] [CrossRef]
  5. Guo, Y.; Zhao, H.; Zhang, S.; Wang, Y.; Chow, D. Modeling and optimization of environment in agricultural greenhouses for improving cleaner and sustainable crop production. J. Clean. Prod. 2020, 285, 124843. [Google Scholar] [CrossRef]
  6. Chen, D.; Neumann, K.; Friedel, S.; Kilian, B.; Chen, M.; Altmann, T.; Klukas, C. Dissecting the Phenotypic Components of Crop Plant Growth and Drought Responses Based on High-Throughput Image Analysis. Plant Cell 2014, 26, 4636–4655. [Google Scholar] [CrossRef] [Green Version]
  7. He, J.Q.; Harrison, R.J.; Li, B. A novel 3D imaging system for strawberry phenotyping. Plant Methods 2017, 13, 1–8. [Google Scholar] [CrossRef]
  8. Sancho-Adamson, M.; Trillas, M.I.; Bort, J.; Fernandez-Gallego, J.A.; Romanyà, J. Use of RGB Vegetation Indexes in Assessing Early Effects of Verticillium Wilt of Olive in Asymptomatic Plants in High and Low Fertility Scenarios. Remote Sens. 2019, 11, 607. [Google Scholar] [CrossRef] [Green Version]
  9. Liu, Z.; Hu, H.; Yu, H.; Yang, X.; Yang, H.; Ruan, C.; Wang, Y.; Tang, J. Relationship between leaf physiologic traits and canopy color indices during the leaf expansion period in an oak forest. Ecosphere 2015, 6, art259. [Google Scholar] [CrossRef]
  10. Grosskinsky, D.K.; Syaifullah, S.J.; Roitsch, T. Integration of multi-omics techniques and physiological phenotyping within a holistic phenomics approach to study senescence in model and crop plants. J. Exp. Bot. 2017, 69, 825–844. [Google Scholar] [CrossRef] [Green Version]
  11. Vasseur, F.; Bresson, J.; Wang, G.; Schwab, R.; Weigel, D. Image-based methods for phenotyping growth dynamics and fitness components in Arabidopsis thaliana. Plant Methods 2018, 14, 1–11. [Google Scholar] [CrossRef]
  12. Barker, J.; Zhang, N.Q.; Sharon, J.; Steeves, R.; Wang, X.; Wei, Y.; Poland, J. Development of a field-based high-throughput mobile phenotyping platform. Comput. Electron. Agric. 2016, 122, 74–85. [Google Scholar] [CrossRef] [Green Version]
  13. Dey, A.K.; Sharma, M.; Meshram, M. An Analysis of Leaf Chlorophyll Measurement Method Using Chlorophyll Meter and Image Processing Technique. Procedia Comput. Sci. 2016, 85, 286–292. [Google Scholar] [CrossRef] [Green Version]
  14. Yadav, S.P.; Ibaraki, Y.; Gupta, S.D. Estimation of the chlorophyll content of micropropagated potato plants using RGB based image analysis. Plant Cell Tissue Organ. Cult. (PCTOC) 2009, 100, 183–188. [Google Scholar] [CrossRef]
  15. Adamsen, F.J.; Pinter, P.J.; Barnes, E.M.; Lamorte, R.L.; Wall, G.W.; Leavitt, S.W.; Kimball, B.A. Measuring Wheat Senescence with a Digital Camera. Crop Sci. 1999, 39, 719–724. [Google Scholar] [CrossRef]
  16. Hu, H.; Liu, H.-Q.; Zhang, H.; Zhu, J.-H.; Yao, X.-G.; Zhang, X.-B.; Zheng, K.-F. Assessment of Chlorophyll Content Based on Image Color Analysis, Comparison with SPAD-502. In Proceedings of the 2nd International Conference on Information Engineering and Computer Science (ICIECS), Wuhan, China, 25–26 December 2010; pp. 476–478. [Google Scholar] [CrossRef]
  17. Ali, M.M.; Al-Ani, A.; Eamus, D.; Tan, D.K.Y.A. New image processing based technique to determine chlorophyll in plants. Am. Eurasian J. Agric. Environ. Sci. 2012, 12, 1323–1328. [Google Scholar] [CrossRef]
  18. Han, W.T.; Sun, Y.; Xu, T.F.; Chen, X.W.; Su, K.O. Detecting maize leaf water status by using digital RGB images. Int. J. Agric. Biol. Eng. 2014, 7, 45–53. [Google Scholar] [CrossRef]
  19. Barbedo, J.G.A. Detection of nutrition deficiencies in plants using proximal images and machine learning: A review. Comput. Electron. Agric. 2019, 162, 482–492. [Google Scholar] [CrossRef]
  20. Humplík, J.F.; Lazár, D.; Husičková, A.; Spíchal, L. Automated phenotyping of plant shoots using imaging methods for analysis of plant stress responses—A review. Plant Methods 2015, 11, 1–10. [Google Scholar] [CrossRef] [Green Version]
  21. Humplík, J.F.; Lazár, D.; Fürst, T.; Husičková, A.; Hýbl, M.; Spíchal, L. Automated integrative high-throughput phenotyping of plant shoots: A case study of the cold-tolerance of pea (Pisum sativum L.). Plant Methods 2015, 11, 20. [Google Scholar] [CrossRef] [Green Version]
  22. Bresson, J.; Bieker, S.; Riester, L.; Doll, J.; Zentgraf, U. A guideline for leaf senescence analyses: From quantification to physiological and molecular investigations. J. Exp. Bot. 2017, 69, 769–786. [Google Scholar] [CrossRef]
  23. Li, L.; Zhang, Q.; Huang, D. A Review of Imaging Techniques for Plant Phenotyping. Sensors 2014, 14, 20078–20111. [Google Scholar] [CrossRef]
  24. Neilson, E.H.; Edwards, A.M.; Blomstedt, C.K.; Berger, B.; Møller, B.L.; Gleadow, R.M. Utilization of a high-throughput shoot imaging system to examine the dynamic phenotypic responses of a C4 cereal crop plant to nitrogen and water deficiency over time. J. Exp. Bot. 2015, 66, 1817–1832. [Google Scholar] [CrossRef]
  25. Bai, G.; Jenkins, S.; Yuan, W.; Graef, G.L.; Ge, Y. Field-Based Scoring of Soybean Iron Deficiency Chlorosis Using RGB Imaging and Statistical Learning. Front. Plant Sci. 2018, 9, 1002–1011. [Google Scholar] [CrossRef] [Green Version]
  26. Gracia-Romero, A.; Kefauver, S.C.; Vergara-Díaz, O.; Zaman-Allah, M.; Prasanna, B.M.; Cairns, J.; Araus, J.L. Comparative Performance of Ground vs. Aerially Assessed RGB and Multispectral Indices for Early-Growth Evaluation of Maize Performance under Phosphorus Fertilization. Front. Plant Sci. 2017, 8, 2004. [Google Scholar] [CrossRef] [Green Version]
  27. Tester, M.; Langridge, P. Breeding Technologies to Increase Crop Production in a Changing World. Science 2010, 327, 818–822. [Google Scholar] [CrossRef]
  28. Bunce, J.A. Responses of stomatal conductance to light, humidity and temperature in winter wheat and barley grown at three concentrations of carbon dioxide in the field. Global Change Biol 2000, 6, 371–382. [Google Scholar] [CrossRef]
  29. Min-Wha, J.; Mohammad, B.A.; Eun-Joo, H.; Kee-Yoeup, P. Photosynthetic pigments, morphology and leaf gas exchange during ex vitro acclimatization of micropropagated CAM Doritaenopsis plantlets under relative humidity and air temperature. Environ. Exp. Bot 2006, 55, 183–194. [Google Scholar] [CrossRef]
  30. Urban, J.; Ingwers, M.W.; McGuire, M.A.; Teskey, R.O. Increase in leaf temperature opens stomata and decouples net photosynthesis from stomatal conductance in Pinus taeda and Populus deltoides x nigra. J. Exp. Bot. 2017, 68, 1757–1767. [Google Scholar] [CrossRef] [PubMed]
  31. Wang, B.; Cai, W.; Li, J.; Wan, Y.; Li, Y.; Guo, C.; Wilkes, A.; You, S.; Qin, X.; Gao, Q.; et al. Leaf photosynthesis and stomatal conductance acclimate to elevated [CO2] and temperature thus increasing dry matter productivity in a double rice cropping system. Field Crop. Res. 2020, 248, 107735. [Google Scholar] [CrossRef]
  32. Gitelson, A.A.; Gritz, Y.; Merzlyak, M.N. Relationships between leaf chlorophyll content and spectral reflectance and algorithms for non-destructive chlorophyll assessment in higher plant leaves. J. Plant Physiol. 2003, 160, 271–282. [Google Scholar] [CrossRef]
  33. Chen, Z.; Wang, F.; Zhang, P.; Ke, C.; Zhu, Y.; Cao, W.; Jiang, H. Skewed distribution of leaf color RGB model and application of skewed parameters in leaf color description model. Plant Methods 2020, 16, 23–28. [Google Scholar] [CrossRef]
  34. Fukuoka, N.; Suzuki, T.; Minamide, K.; Hamada, T. Effect of Shading on Anthocyanin and Non-flavonoid Polyphenol Biosynthesis of Gynura bicolor Leaves in Midsummer. HortScience 2014, 49, 1148–1153. [Google Scholar] [CrossRef] [Green Version]
  35. Ren, J.; Guo, S.S.; Xin, X.L.; Chen, L. Changes in volatile constituents and phenols from Gynura bicolor DC grown in elevated CO2 and LED lighting. Sci. Hortic. 2014, 175, 243–250. [Google Scholar] [CrossRef]
  36. Shimizu, Y.; Maeda, K.; Kato, M.; Shimomura, K. Methyl jasmonate induces anthocyanin accumulation in Gynura bicolor cultured roots. Vitr. Cell Dev. Biol. Plant 2010, 46, 460–465. [Google Scholar] [CrossRef]
  37. Wu, C.-C.; Chang, Y.-P.; Wang, J.-J.; Liu, C.-H.; Wong, S.-L.; Jiang, C.-M.; Hsieh, S.-L. Dietary administration of Gynura bicolor (Roxb. Willd.) DC water extract enhances immune response and survival rate against Vibrio alginolyticus and white spot syndrome virus in white shrimp Litopeneaus vannamei. Fish Shellfish Immunol. 2014, 42, 25–33. [Google Scholar] [CrossRef]
  38. Dobrescu, A.; Scorza, L.C.T.; Tsaftaris, S.A.; McCormick, A.J. A “Do-It-Yourself” phenotyping system: Measuring growth and morphology throughout the diel cycle in rosette shaped plants. Plant Methods 2017, 13, 1–12. [Google Scholar] [CrossRef] [Green Version]
  39. Sun, Y.; Zhu, S.; Yang, X.; Weston, M.V.; Wang, K.; Shen, Z.; Xu, H.; Chen, L. Nitrogen diagnosis based on dynamic characteristics of rice leaf image. PLoS ONE 2018, 13, e0196298. [Google Scholar] [CrossRef] [Green Version]
  40. Gai, J.Y. Experimental Statistical Method; China Agriculture Press: Beijing, China, 2000; pp. 193–208. [Google Scholar]
  41. Hand, D. Effects of Atmospheric Humidity on Greenhouse Crops. Acta Hortic. 1988, 229, 143–158. [Google Scholar] [CrossRef]
  42. Shamshiri, R.R.; Jones, J.W.; Thorp, K.; Ahmad, D.; Man, H.C.; Taheri, S. Review of optimum temperature, humidity, and vapour pressure deficit for microclimate evaluation and control in greenhouse cultivation of tomato: A review. Int. Agrophys. 2018, 32, 287–302. [Google Scholar] [CrossRef]
  43. Li, G.; Tang, L.; Zhang, X.; Dong, J.; Xiao, M. Factors affecting greenhouse microclimate and its regulating techniques: A review. IOP Conf. Ser. Earth Environ. Sci. 2018, 167, 012019. [Google Scholar] [CrossRef]
  44. Singh, M.C.; Singh, J.P.; Pandey, S.K.; Cutting, N.G.; Sharma, P.; Shrivastav, V.; Sharma, P. A Review of Three Commonly Used Techniques of Controlling Greenhouse Microclimate. Int. J. Curr. Microbiol. Appl. Sci. 2018, 7, 3491–3505. [Google Scholar] [CrossRef]
  45. Zhang, S.H.; Guo, Y.; Zhao, H.J.; Wang, Y.; Chow, D.; Fang, Y. Methodologies of control strategies for improving energy efficiency in agricultural greenhouses. J. Clean. Prod. 2020, 274, 122695. [Google Scholar] [CrossRef]
  46. Zhu, J.-X.; Deng, J.-S.; Shi, Y.-Y.; Chen, Z.-L.; Han, N.; Wang, K. Diagnoses of rice nitrogen status based on characteristics of scanning leaf. Spectrosc. Spect. Anal. 2009, 29, 2171–2175. [Google Scholar] [CrossRef]
  47. Chaerle, L.; Hagenbeek, D.; Vanrobaeys, X.; Van Der Straeten, D. Early detection of nutrient and biotic stress in Phaseolus vulgaris. Int. J. Remote Sens. 2007, 28, 3479–3492. [Google Scholar] [CrossRef]
  48. Gous, P.W.; Meder, R.; Fox, G.P. Near Infrared Spectral Assessment of Stay-Green Barley Genotypes under Heat Stress. J. Near Infrared Spectrosc. 2015, 23, 145–153. [Google Scholar] [CrossRef]
  49. Cai, J.; Okamoto, M.; Atieno, J.; Sutton, T.; Li, Y.; Miklavcic, S.J. Quantifying the Onset and Progression of Plant Senescence by Color Image Analysis for High Throughput Applications. PLoS ONE 2016, 11, e0157102. [Google Scholar] [CrossRef] [Green Version]
  50. Schmalko, M.E.; Scipioni, P.G.; Ferreyra, D. Effect of Water Activity and Temperature in Color and Chlorophylls Changes in Yerba Mate Leaves. Int. J. Food Prop. 2005, 8, 313–322. [Google Scholar] [CrossRef]
  51. Shibayama, M.; Sakamoto, T.; Takada, E.; Inoue, A.; Morita, K.; Takahashi, W.; Kimura, A. Estimating Paddy Rice Leaf Area Index with Fixed Point Continuous Observation of Near Infrared Reflectance Using a Calibrated Digital Camera. Plant Prod. Sci. 2011, 14, 30–46. [Google Scholar] [CrossRef] [Green Version]
  52. Sritarapipat, T.; Rakwatin, P.; Kasetkasem, T. Automatic Rice Crop Height Measurement Using a Field Server and Digital Image Processing. Sensors 2014, 14, 900–926. [Google Scholar] [CrossRef] [Green Version]
  53. Sun, Y.Y.; Shen, S.H. Research progress in application of crop growth models. Chin. J. Agrometeorol. 2019, 40, 444–459. [Google Scholar] [CrossRef]
  54. Yu, Z.H.; Cao, Z.G.; Wu, X.; Bai, X.D.; Qin, Y.M.; Zhuo, W.; Xiao, Y.; Zhang, X.F.; Xue, H.X. Automatic image-based detection technology for two critical growth stages of maize: Emergence and three-leaf stage. Agric. For. Meteorol. 2013, 174, 65–84. [Google Scholar] [CrossRef]
  55. Li, Y.; Cao, Z.; Lu, H.; Xiao, Y.; Zhu, Y.; Cremers, A.B. In-field cotton detection via region-based semantic image segmentation. Comput. Electron. Agric. 2016, 127, 475–486. [Google Scholar] [CrossRef]
Figure 1. Schematic map of microclimate observation and digital camera in the glass greenhouse.
Figure 1. Schematic map of microclimate observation and digital camera in the glass greenhouse.
Atmosphere 13 00890 g001
Figure 2. Flow map of Cutting and Denoising of the Image.
Figure 2. Flow map of Cutting and Denoising of the Image.
Atmosphere 13 00890 g002
Figure 3. RGB model color gradation distribution of single leaf, plant canopy, and population canopy images of Begonia Fimbristipula Hance (BFH). The cumulative frequency histogram of red, green, and blue channels, as well as gray-level images, were drawn using the imhist function of MATLAB. The X-axis is the color gradation value, and the Y-axis is the cumulative frequency.
Figure 3. RGB model color gradation distribution of single leaf, plant canopy, and population canopy images of Begonia Fimbristipula Hance (BFH). The cumulative frequency histogram of red, green, and blue channels, as well as gray-level images, were drawn using the imhist function of MATLAB. The X-axis is the color gradation value, and the Y-axis is the cumulative frequency.
Atmosphere 13 00890 g003
Figure 4. RGB model color gradation distribution of population canopy images of BFH on 3 different dates. The cumulative frequency histogram of red, green, and blue channels, as well as gray-level images, were drawn using the imhist function of MATLAB. The X-axis is the color gradation value, and the Y-axis is the cumulative frequency.
Figure 4. RGB model color gradation distribution of population canopy images of BFH on 3 different dates. The cumulative frequency histogram of red, green, and blue channels, as well as gray-level images, were drawn using the imhist function of MATLAB. The X-axis is the color gradation value, and the Y-axis is the cumulative frequency.
Atmosphere 13 00890 g004
Figure 5. Prediction accuracy interval analysis of the population canopy color-meteorological response models of the modeling group. The Boxplot function of MATLAB was used to draw the accuracy percentile distribution diagram of the population canopy CGSD parameters prediction model. The blue box in the figure represents the prediction accuracy of the model distributed between 25% and 75%. The shorter the box, the better the degree of data concentration. The red line in the box represents the median distribution of prediction accuracy, and the dotted line at both ends outside the box represents the endpoint value (maximum and minimum) of the prediction data distribution. The shorter the distance between the two endpoint values, the better the degree of data concentration. The red cross outside the box represents abnormal data. The fewer abnormal data, the better the degree of data distribution concentration.
Figure 5. Prediction accuracy interval analysis of the population canopy color-meteorological response models of the modeling group. The Boxplot function of MATLAB was used to draw the accuracy percentile distribution diagram of the population canopy CGSD parameters prediction model. The blue box in the figure represents the prediction accuracy of the model distributed between 25% and 75%. The shorter the box, the better the degree of data concentration. The red line in the box represents the median distribution of prediction accuracy, and the dotted line at both ends outside the box represents the endpoint value (maximum and minimum) of the prediction data distribution. The shorter the distance between the two endpoint values, the better the degree of data concentration. The red cross outside the box represents abnormal data. The fewer abnormal data, the better the degree of data distribution concentration.
Atmosphere 13 00890 g005
Figure 6. Prediction accuracy interval analysis of the population canopy color-meteorological response models of the between-group. The Boxplot function of MATLAB was used to draw the accuracy percentile distribution diagram of the population canopy CGSD parameters prediction model. The blue box in the figure represents the prediction accuracy of the model distributed between 25% and 75%. The shorter the box, the better the degree of data concentration. The red line in the box represents the median distribution of prediction accuracy, and the dotted line at both ends outside the box represents the endpoint value (maximum and minimum) of the prediction data distribution. The shorter the distance between the two endpoint values, the better the degree of data concentration. The red cross outside the box represents abnormal data. The fewer abnormal data, the better the degree of data distribution concentration.
Figure 6. Prediction accuracy interval analysis of the population canopy color-meteorological response models of the between-group. The Boxplot function of MATLAB was used to draw the accuracy percentile distribution diagram of the population canopy CGSD parameters prediction model. The blue box in the figure represents the prediction accuracy of the model distributed between 25% and 75%. The shorter the box, the better the degree of data concentration. The red line in the box represents the median distribution of prediction accuracy, and the dotted line at both ends outside the box represents the endpoint value (maximum and minimum) of the prediction data distribution. The shorter the distance between the two endpoint values, the better the degree of data concentration. The red cross outside the box represents abnormal data. The fewer abnormal data, the better the degree of data distribution concentration.
Atmosphere 13 00890 g006
Figure 7. Prediction accuracy interval analysis of the population canopy color-meteorological response models of the outside-group. The Boxplot function of MATLAB was used to draw the accuracy percentile distribution diagram of the population canopy CGSD parameters prediction model. The blue box in the figure represents the prediction accuracy of the model distributed between 25% and 75%. The shorter the box, the better the degree of data concentration. The red line in the box represents the median distribution of prediction accuracy, and the dotted line at both ends outside the box represents the endpoint value (maximum and minimum) of the prediction data distribution. The shorter the distance between the two endpoint values, the better the degree of data concentration. The red cross outside the box represents abnormal data. The fewer abnormal data, the better the degree of data distribution concentration.
Figure 7. Prediction accuracy interval analysis of the population canopy color-meteorological response models of the outside-group. The Boxplot function of MATLAB was used to draw the accuracy percentile distribution diagram of the population canopy CGSD parameters prediction model. The blue box in the figure represents the prediction accuracy of the model distributed between 25% and 75%. The shorter the box, the better the degree of data concentration. The red line in the box represents the median distribution of prediction accuracy, and the dotted line at both ends outside the box represents the endpoint value (maximum and minimum) of the prediction data distribution. The shorter the distance between the two endpoint values, the better the degree of data concentration. The red cross outside the box represents abnormal data. The fewer abnormal data, the better the degree of data distribution concentration.
Atmosphere 13 00890 g007
Table 1. Color gradation skewness-distribution (CGSD) parameters.
Table 1. Color gradation skewness-distribution (CGSD) parameters.
Channel CGSD Parameters
MeanMedianModeSkewnessKurtosis
Red Channel (R)RMeanRMedianRModeRSkewnessRKurtosis
Green Channel (G)GMeanGMedianGModeGSkewnessGKurtosis
Blue Channel (B)BMeanBMedianBModeBSkewnessBKurtosis
Gray Image (Y)YMeanYMedianYModeYSkewnessYKurtosis
Table 2. Spearman correlation coefficient of population canopy CGSD parameters of BFH and microclimate factors.
Table 2. Spearman correlation coefficient of population canopy CGSD parameters of BFH and microclimate factors.
TdmRHdmTGdmTSdm−10cVPdmTDdmGRdtPARdtATAGRAPAR
RMean−0.2490.606 **0.278 *0.629 **−0.025−0.029−0.24−0.251−0.865 **−0.856 **−0.860 **
RMedian−0.2530.605 **0.280 *0.641 **−0.024−0.033−0.227−0.235−0.843 **−0.836 **−0.840 **
RMode−0.2570.584 **0.315 *0.690 **−0.041−0.045−0.247−0.251−0.815 **−0.813 **−0.816 **
RSkewness0.225−0.586 **−0.292 *−0.641 **0.0030.010.2140.2240.853 **0.846 **0.850 **
RKurtosis0.142−0.503 **−0.265−0.519 **−0.054−0.0470.1130.1310.792 **0.774 **0.779 **
GMean−0.282 *0.633 **0.2360.588 **−0.067−0.053−0.344 *−0.358 **−0.950 **−0.947 **−0.949 **
GMedian−0.289 *0.641 **0.2460.613 **−0.068−0.057−0.343 *−0.355 *−0.931 **−0.929 **−0.931 **
GMode−0.295 *0.607 **0.250.623 **−0.088−0.079−0.326 *−0.334 *−0.901 **−0.901 **−0.902 **
GSkewness0.231−0.598 **−0.307 *−0.668 **0.0080.010.2570.2680.866 **0.864 **0.867 **
GKurtosis0.157−0.507 **−0.261−0.526 **−0.031−0.0330.1760.1970.855 **0.844 **0.848 **
BMean−0.332 *0.625 **0.2350.632 **−0.105−0.112−0.284 *−0.293 *−0.794 **−0.789 **−0.793 **
BMedian−0.317 *0.619 **0.2520.654 **−0.088−0.097−0.272−0.278 *−0.790 **−0.787 **−0.790 **
BMode−0.2630.577 **0.321 *0.716 **−0.047−0.054−0.292 *−0.298 *−0.750 **−0.751 **−0.754 **
BSkewness0.243−0.588 **−0.294 *−0.656 **0.0230.0280.2420.2520.864 **0.860 **0.863 **
BKurtosis0.205−0.517 **−0.238−0.527 **0.0120.0170.1640.1810.821 **0.808 **0.812 **
YMean−0.281 *0.632 **0.250.610 **−0.06−0.053−0.313 *−0.327 *−0.922 **−0.917 **−0.920 **
YMedian−0.283 *0.630 **0.2580.630 **−0.059−0.056−0.304 *−0.314 *−0.899 **−0.896 **−0.898 **
YMode−0.2710.608 **0.285 *0.665 **−0.056−0.051−0.316 *−0.323 *−0.884 **−0.882 **−0.884 **
YSkewness0.229−0.593 **−0.305 *−0.661 **0.0050.010.2390.250.858 **0.854 **0.857 **
YKurtosis0.159−0.508 **−0.259−0.524 **−0.033−0.0310.1540.1730.832 **0.818 **0.822 **
Note: A Spearman correlation analysis of SPSS software was used on 20 CGSD parameters (RMean, RMedian, RMode, RSkewness, RKurtosis, GMean, GMedian, GMode, GSkewness, GKurtosis, BMean, BMedian, BMode, BSkewness, BKurtosis, YMean, YMedian, YMode, YSkewness, and YKurtosis) and the microclimate factors including daily meteorological factors (Tdm, RHdm, TGdm, TSdm−10c, VPdm, TDdm, GRdt, and PARdt) and photothermal accumulation factors (AT, AGR, and APAR) (n = 51). ** indicates significant correlation according to a two-tailed test (p < 0.01), * indicates significant correlation according to a two-tailed test (p < 0.05), the same as below.
Table 3. Population canopy color-meteorological response models of BFH and their goodness of fit.
Table 3. Population canopy color-meteorological response models of BFH and their goodness of fit.
ModelR-SquareAdjusted
R-Square
RMSEF ValueSignificance F
RMeanY1 = 11.556 − 0.017x1 + 0.524x2 + 0.549x30.8480.8392.45687.6340.000
RMedianY2 = − 28.395 − 0.015x1 + 0.629x2 + 0.703x3 + 1.871x40.8680.8572.63075.8600.000
RModeY3 = 14.284 − 0.024x5 + 2.596x40.7460.7353.66870.4580.000
RSkewnessY4 = 2.100 + 0.001x1 − 0.199x4 + 0.107x60.8030.7910.14263.9680.000
RKurtosisY5 = 5.371 − 0.010 x1 − 0.024x7 − 0.028 x20.7060.6880.38237.7030.000
GMeanY6 = 27.424 − 0.022x1 + 0.464x2 + 0.349x30.9540.9511.609327.3370.000
GMedianY7 = 34.402 − 0.036x5 + 0.366 x20.9060.9022.548230.3590.000
GModeY8 = 30.294 − 0.037x5 + 0.342 x20.8450.8393.456130.8210.000
GSkewnessY9 = 1.501 + 0.030x5 − 0.052x7 − 0.013 x20.8310.8200.13876.9000.000
GKurtosisY10 = 2.824 + 0.002x10.7320.7260.310133.5390.000
BMeanY11 = 30.430 − 0.013x1 + 0.341x20.7030.6913.13656.7870.000
BMedianY12 = −6.307 − 0.012x1 + 0.380x2 + 2.094 x40.7520.7363.47047.4830.000
BModeY13 = 4.464 − 0.020 x5 + 3.267x40.6950.6834.02854.7810.000
BSkewnessY14 = 1.716 + 0.001x1 − 0.083x40.7910.7820.14590.1640.000
BKurtosisY15 = 2.874 + 0.002x10.6740.6680.346101.3770.000
YMeanY16 = 38.927 − 0.019x1 + 0.307x20.8920.8882.240198.3720.000
YMedianY17 = 30.819 − 0.021x1 + 0.368x20.8550.8482.956140.9730.000
YModeY18 = 26.665 − 0.032 x5 + 2.144 x40.8250.8183.440113.2130.000
YSkewnessY19 = 1.818 + 0.001x1 − 0.092 x40.7840.7750.15687.3020.000
YKurtosisY20 = 2.901 + 0.002x10.6920.6860.355110.0290.000
Note: x1 is the accumulated temperature, x2 is the daily mean relative humidity, x3 is the daily total photosynthetically active radiation, x4 is the daily mean soil temperature at 10 cm, x5 is the accumulated photosynthetically active radiation, x6 is the daily minimum ground temperature, x7 is the accumulated global radiation.
Table 4. Analytical results of prediction accuracy of population canopy color-meteorological response models.
Table 4. Analytical results of prediction accuracy of population canopy color-meteorological response models.
ModelModeling GroupBetween-Group PredictionOutside-Group PredictionAverage
Prediction Sample DataEliminate Abnormal DataPredictive AccuracyPrediction Sample DataEliminate Abnormal DataPredictive AccuracyPrediction Sample DataEliminate Abnormal DataPredictive Accuracy
RMeanY151096.27%51092.89%15090.92%94.11%
RMedianY251096.02%51091.98%15084.90%92.84%
RModeY351092.89%51092.47%15087.87%92.06%
RSkewnessY451080.98%511053.21%15370.30%68.80%
RKurtosisY551090.48%51072.53%15073.56%80.49%
GMeanY651097.98%51093.55%15093.71%95.50%
GMedianY751096.29%51091.97%15080.08%92.33%
GModeY851094.14%51092.38%15077.33%91.22%
GSkewnessY951176.29%51952.40%15735.12%63.73%
GKurtosisY1051092.98%51076.64%15086.75%85.06%
BMeanY1151095.46%51092.05%15084.66%92.59%
BMedianY1251094.53%51090.57%15085.42%91.64%
BModeY1351092.66%51092.65%15084.46%91.60%
BSkewnessY1451182.21%51856.11%15382.39%72.36%
BKurtosisY1551092.17%51074.21%15084.12%83.31%
YMeanY1651096.95%51093.33%15089.60%94.43%
YMedianY1751090.26%51093.64%15067.68%88.84%
YModeY1851094.11%51093.15%15087.27%92.82%
YSkewnessY1951280.43%511058.03%15384.64%73.81%
YKurtosisY2051092.56%51076.04%15086.53%84.59%
Note: Predictive accuracy = (1 − |predictive value-measured value|/measured) × 100%.
Table 5. Comparison of multiple stepwise regression models and spatial bidirectional models for the skewness.
Table 5. Comparison of multiple stepwise regression models and spatial bidirectional models for the skewness.
ParametersModelR-SquareAdjusted
R-Square
RMSESignificant F
RSkewnessY4 0.8030.7910.1420.000
Z1Equation (4)0.8840.8740.1110.000
GSkewnessY9 0.8310.8200.1380.000
Z2Equation (5)0.9020.8930.1070.000
BSkewnessY14 0.7910.7820.1450.000
Z3Equation (6)0.9000.8910.1020.000
YSkewnessY19 0.7840.7750.1560.000
Z4Equation (7)0.8940.8840.1110.000
Table 6. Comparison of the prediction accuracy of the multiple stepwise regression models and spatial bidirectional models for the skewness.
Table 6. Comparison of the prediction accuracy of the multiple stepwise regression models and spatial bidirectional models for the skewness.
ParametersModelModeling GroupBetween-GroupOutside-GroupAverage
Prediction Sample DataEliminate Abnormal DataPredictive AccuracyPrediction Sample DataEliminate Abnormal DataPredictive AccuracyPrediction Sample DataEliminate Abnormal DataPredictive Accuracy
RSkewnessY451080.98%511053.21%15370.30%68.80%
Z151090.02%51160.56%15375.58%77.76%
GSkewnessY951176.29%51952.40%15735.12%63.73%
Z251088.31%51160.82%15373.71%74.60%
BSkewnessY1451182.21%51856.11%15382.39%72.36%
Z351090.16%51161.27%15375.61%75.83%
YSkewnessY1951280.43%511058.03%15384.64%73.81%
Z451089.05%51160.10%15382.34%75.53%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, P.; Chen, Z.; Wang, F.; Wang, R.; Bao, T.; Xie, X.; An, Z.; Jian, X.; Liu, C. Response of Population Canopy Color Gradation Skewed Distribution Parameters of the RGB Model to Micrometeorology Environment in Begonia Fimbristipula Hance. Atmosphere 2022, 13, 890. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos13060890

AMA Style

Zhang P, Chen Z, Wang F, Wang R, Bao T, Xie X, An Z, Jian X, Liu C. Response of Population Canopy Color Gradation Skewed Distribution Parameters of the RGB Model to Micrometeorology Environment in Begonia Fimbristipula Hance. Atmosphere. 2022; 13(6):890. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos13060890

Chicago/Turabian Style

Zhang, Pei, Zhengmeng Chen, Fuzheng Wang, Rong Wang, Tingting Bao, Xiaoping Xie, Ziyue An, Xinxin Jian, and Chunwei Liu. 2022. "Response of Population Canopy Color Gradation Skewed Distribution Parameters of the RGB Model to Micrometeorology Environment in Begonia Fimbristipula Hance" Atmosphere 13, no. 6: 890. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos13060890

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