Research Article  Open Access
Tingting Qian, Xiuguo Zheng, Juan Yang, Yeying Xu, Yan Wang, Qiang Zhou, Shenglian Lu, Xiaotao Ding, "Optimal Utilization of Light Energy in Semiclosed Greenhouse Using ThreeDimensional Cucumber Model", Scientific Programming, vol. 2020, Article ID 8855063, 12 pages, 2020. https://doi.org/10.1155/2020/8855063
Optimal Utilization of Light Energy in Semiclosed Greenhouse Using ThreeDimensional Cucumber Model
Abstract
In the east of China, low temperature and light energy in winter are the main factors for the decline in cucumber yield, as well as in greenhouses without supplementary light. Optimal utilization of light energy is critical to increase cucumber yield. In this study, experimental measurements were conducted in two scenarios, April to May (AprMay) and November to December (NovDec) 2015, respectively, to analyze leaf development, dry matter accumulation, and yield of cultivated cucumber. Statistical analysis showed that leaves grew in NovDec had larger leaf area and lower dry matter than leaves grew in AprMay. This revealed that the dry matter accumulation rate per unit area was lower in winter. To be precise, the yield 0.174 kg/m^{2} per day in NovDec was 35.3% lower than the yield in AprMay. Environmental monitoring data showed that there was no significant difference in the average temperature between two scenarios, but the light intensity in NovDec was only 2/3 of that in AprMay. Threedimensional (3D) cucumber canopy models were used in this study to quantify the effects of weak light on dry matter production in NovDec. Three 3D canopies of cucumber were reconstructed with 20, 25, and 30 leaves per plant, respectively, by using a parametric modeling method. Light interception of three canopies from 8:00 to 15:00 on 4 November 2015 was simulated by using the radiositygraphic combined model (RGM) with an hourly time step. CO_{2} assimilation per unit area was calculated using the FvCB photosynthetic model. As a result, the effects of light intensity and CO_{2} concentration on the photosynthetic rate were considered. The leaf photosynthesis simulation result showed that during the daytime in winter, the RuBP regenerationlimited assimilation A_{j} was always less than the Rubiscolimited assimilation A_{c}. This means that the limiting factor affecting the photosynthesis rate in winter was rather light intensity. As the CO_{2} concentration in the greenhouse was utmost in the morning, increasing the light intensity and therefrom increasing the canopy light interception at this time will be highly beneficial to increase the yield production. Through a comparative analysis of photosynthetic characteristics in these three virtual 3D canopies, the 25leaf canopy was the bestperformed canopy structure in photosynthetic production in winter. This study provides an insight into the light deficiency for yield production in winter and a solution to make optimal use of light in the greenhouse.
1. Introduction
Cucumber is one of the most popular crops grown in facility agriculture. It requires relatively high light intensity and temperature. Low temperature and deficient light energy in winter in the east and north of China are the main factors for the decline of cucumber yield. With the promotion and application of semiclosed greenhouses in China, the excess solar energy in summer was collected and stored, and this energy should be reused to heat the greenhouse in winter [1]. As a result, the impact of low temperature in winter on yield has been greatly improved [2]. Nevertheless, in greenhouses without supplementary light, low light intensity remains an important cause affecting winter production [1]. For most greenhouse vegetables, 1% light increment results in 0.7–1% increase in harvestable product [3]. Therefore, optimal utilization of light energy should be critical for increasing the yield of cucumber in winter.
Optimization and control of the light environment in greenhouse have been widely studied. In order to increase the incident light intensity, supplementary light sources are often added on top of the canopy during the growing period as early as the 1980s [4]. For most greenhouse fruiting vegetables, such as cucumber and tomato, vegetative growth is concentrated on the top and middle of the canopy while reproductive growth occurs at the middle and bottom. In recent years, lighting strategies including light source type and distribution, spectral composition, and lighting duration have been widely studied. Hovi et al. [5] have shown that interlighting by using HPS (highpressure sodium) can increase by 10% of cucumber fruit yield. However, these kinds of light sources are very difficult to be used as interlighting in narrow row arrangements, as they produce high heat. Hao et al. [6] studied the effect of HPS and LED (lightemitting diodes) light source and light position (on the top or inside of the canopy) on minicucumber growth rate and fruit yield. The result shows that the plants grown with LED interlighting achieved higher energy use efficiency than that with the top HPS. And, once canopy closure achieves, increasing total canopy light interception, especially increasing middle and lower light interception, is more conducive to yield [7].
In cucumber canopy, light interception is directly affected by canopy structure, which depends on plant growth and cultivation management. Accurate measurement of total light interception in the canopy has always been a great challenge. In previous studies, Beer–Lambert’s law was widely used to calculate canopy light interception, but this law assumes that the canopy is a uniform medium, and ignored the influence of leaf distribution and inclination angle on light transmission [8]. For horticultural crops planted in rows, the canopy structure exhibits distinct heterogeneity. Since light distributes significantly different between and within crop rows, the orientation of the leaf directly determines its light interception [9]. The 3D structural model, which delivers a structure close to the real canopy, can accurately describe the distribution and location of organs in 3D space. Combing with ray casting, radiosity, and Monte Carlo raytracing approaches, local light conditions and light interception can be simulated at the leaf level and upscale to the canopy level [10]. Thus, the accuracy of the light interception estimation can be greatly improved.
Recent studies have shown that the 3D structural model is extremely useful for light interceptionrelated calculation, for example, evaluating the effect of leaf area distribution and leaf angles on canopy light interception for “ideotype” selection [11, 12] and optimizing canopy structure pruning management [13]. In many FSPM (functional and structural plant models) research studies, 3D structural models were always combined with photosynthetic models and dry matter distribution models at the aim to design crops with higher yield potential and to enhance resource use efficiency [14–16].
In the present study, we combined 3D cucumber canopy model, raytracing method, and photosynthesis model together to analyze light conditions inside the canopy. By virtual experiments, we discussed a solution in greenhouse management to increase the photosynthetic production of cucumber crops in winter.
2. Materials and Methods
2.1. Field Experiments and Measurement on Plants
Field experiments were conducted in May and November 2015, respectively, in a semiclosed Venlotype greenhouse at Shanghai Chongming Island (lat. 31°34′N, long. 121°41′E). The detailed information of the experimental greenhouse can be referenced in Ding et al. [2]. Minicucumbers (Cucumis Sativus L. cv. Deltastar) were planted on rockwool slabs along two parallel lines and each plant is located with a growing wire at a stem density of 2.8 stems/m^{2}. Air temperature, relative humidity, and CO_{2} concentration in the greenhouse were recorded automatically at 5 min intervals using one climate sensor (Priva, De Lier, the Netherlands) per greenhouse, placed at the top of the canopy.
2.1.1. Measurement 1: Organ Developing Measurement
Two measurements were conducted in 2015 from 15^{th} April to 14^{th} May and from 4^{th} November to 3 December, and for convenience, AprMay and NovDec will be used to refer the two experimental scenarios, respectively, hereinafter. Every day in the experiment, five newly appeared phytomers, including an internode with the directly attached leaves and petioles, were selected to be marked as the initial point of organ expansion. In total, 150 phytomers corresponding to 450 organs were marked in the thirty consecutive days. Morphological measurements were performed on all marked phytomers, respectively, on May 16, 2015, and December 3, 2015. The measurement contents include leaf length (LL), leaf width (LW), petiole length (PL), petiole diameter (PD), internode length (IL), and internode diameter (ID). Leaf length was measured from the lamina tip to the intersection end of petiole along the lamina midrib. Leaf width was measured from tip to tip between the widest lamina lobes, as referenced in Cho et al. [17]. The diameter of the petiole and internode was measured in the middle of the organ. The manual measurement error is ±1 mm. After the measurement, the leaves and petioles were destructively dried and weighted.
2.1.2. Measurement 2: Canopy Structure Measurement
Canopy structure measurements were conducted on 4 November 2015. A total of 8 cucumber plants from 2 (rows) (plants) adjacent plants were selected, and all organs of the whole plant were measured manually. The measurement contents include leaf length, leaf width, petiole length, petiole diameter, internode length, and internode diameter. After the measurement, organs were segmented, dried, and weighed one by one.
2.2. 3D Reconstruction of the Canopy
2.2.1. Organ Expansion Modeling
Based on measurement 1, the leastsquares method was used to fit statistical models for organ development. We used logistics function for leaf length, leaf width, petiole diameter, and internode diameter and Boltzmann function for petiole length, internode length, leaf dry weight (LDW), and petiole dry weight (PDW). The parameters of descriptive models for organ development are listed in Tables 1 and 2. The sum of the thermal time (GDD) is calculated independently: GDD = , where T is the temperature measured by the sensor, T_{base} is the base temperature, 10°C for minicucumber, and 288 is an invariant for the number of measurements per day. The duration is from the initial date till the end of measurement and valid temperature ranging from 10 to 40°C is applied for cucumber.


2.2.2. Parametric Construction of 3D Cucumber Model
Parameterized modeling method was used to construct the 3D cucumber canopies. The shape of each organ was determined by morphological parameters, which were generated by using the equations described above. Organ morphology was composed of skeletons and triangular mesh. Bsplines curve was used to represent these skeletons in the computer. The control points on Bspline curve can determine the size of the organs (Figure 1), and the spatial distance of the control points was determined by the organ morphological parameters, which in turn determined the topology of the plant [18, 19].
2.2.3. 3D Canopy Construction for Cucumber
Three 3D canopies of cucumber were reconstructed with 20, 25, and 30 leaves per plant, respectively, by using the parameterized modeling method. Each virtual canopy was consisted of sixteen plants (4 rows4 plants). The interplant distance in one line was 40 cm, the line spacing was 70 cm (inside) and 110 cm (path), respectively (Figure 2).
(a)
(b)
(c)
2.3. Light Interception Modeling
Light interception of three canopies from 8:00 to 15:00 on 4 November 2015 was simulated by using the radiositygraphic combined model (RGM) with an hourly time step. Distribution of 3D diffuse photosynthetically active radiation (DPAR) in the canopy was simulated discretely by dividing the sky hemisphere into a grid of fine cells that allowed for the anisotropic distribution of DPAR over the sky hemisphere. Diffuse radiation was calculated using the method of Wen et al. [10]. It was assumed that a leaf absorbs 87%, transmits 7%, and reflects 6% of the incident photosynthetically active radiation [15]. Total PAR and proportion of diffuse light intensity were measured using SunScan Canopy Analysis System (Delta Company, UK).
2.4. Leaf Photosynthesis Modeling
FvCB model (Farquhar–von Caemmerer–Berry biochemical photosynthesis model) was applied in leaf photosynthesis modeling. According to Farquhar et al. [20], assimilation rate A is the minimum of the Rubiscolimited (A_{c}, μ mol·CO_{2}·m^{−2}·s^{−1}) or RuBP regenerationlimited (A_{j}, μ mol·CO_{2}·m^{−2}·s^{−1}) photosynthesis rate:where C_{c} (chloroplastic CO_{2} concentration, μ mol·mol^{−1}) is calculated by using the following equation:
When substituting equation (3) into equations (1) and (2), we get the quadratic equation of variable A_{c} and A_{j}, and the solutions are as follows:
J (rate of electron transport, μ mol·m^{−2}s^{−1}) was calculated by the following equation:where is the maximum rate of Rubisco carboxylation (130.84 μ mol·CO_{2}·m^{−2}·s^{−1}); Γ is the CO_{2} compensation point in the absence of dark respiration (43.02 μ mol·mol^{−1}); K_{c} (404 μ mol·mol^{−1}) and K_{o} (278 mmol·mol^{−1}) are Michaelis–Menten constants of Rubisco for CO_{2} and O_{2}; O is the mole fraction of O_{2} at the site of carboxylation (210 mmol·mol^{−1}); R_{d} is the respiration rate (1.08 μ mol·CO_{2}·m^{−2}·s^{−1}); is mesophyll conductance (0.8 mol·m^{−2}·s^{−1}); C_{i} is ambient CO_{2} concentration measured by CO_{2} sensor; J_{max} is the maximum electron transport rate (125.66 μ mol·m^{−2}·s^{−1}); K is a constant describing the conversion efficiency of intercepted PAR (I) to J (0.425 mol·mol^{−1} photon); and θ is a constant convexity factor describing the response of J to I (0.7) [15].
2.5. Model Calibration and Validation
The structural parameters measured in measurement 2 were used for model validation. The rootmeansquared deviation (RMSD), bias, and accuracy (%) were determined following the approach described by Kahlen and Stützel [21]:where x_{i} and y_{i} are the simulated and measured values, respectively. n is the number of simulated or measured values.
3. Results
3.1. Organ Development Analysis
The growth of cucumber organs conforms to the sigmoid growth curve. The maximum asymptote (y = max), the maximum growth rate point (x_{0}, y_{0}), and the slope () at the inflection point can be used as characteristic values describing its growth process (Figure 3). The curve characteristic values based on the results of measurement 1 are shown in Table 3 and Figure 4. The of the leaf area expansion curves in AprMay was 4.58, while the maximum expansion rate of leaf area occurs when the GDD reaches 92.75 Cd (Figure 4(a)). The maximum expansion rate of the petiole and internode length in AprMay was 0.21 and 0.30, while the maximum expansion rate of the petiole appears at 67 Cd, and the maximum expansion rate of the internode appears at 41 Cd (Figure 4(c)). It indicated that internode reaches the maximum expansion rate before leaf and petiole reaching the maximum expansion rate, and the elongation of the internode is nearly completed. The dry matter accumulation of the petiole was less than that of the leaf, and the maximum values were 0.68 g and 1.91 g, respectively (Figure 4(b)). The dry matter of the leaf was 3 times the weight of the petiole.

(a)
(b)
(c)
(d)
The leaf expansion curve of cucumber in NovDec was similar to that of AprMay, and the leaf length expansion rate was lower than that of leaf width. Different from AprMay, the maximum leaf expansion rate in NovDec occurred when GDD reached 76 Cd, earlier than AprMay. The maximum expansion rates of petiole length and internode length were both 0.21. The maximum expansion rate of the petiole appeared at 52 Cd, while the maximum expansion rate of the internode appeared at 27 Cd. In NovDec, the maximum dry matter accumulation of petiole was 0.44 g and 1.38 g, and the maximum dry matter accumulation in winter leaf was significantly smaller than that in AprMay (Table 3).
3.2. Analysis of Differences in Organ Development in Different Scenarios
The rate of leaf area expansion in AprMay was less than that in NovDec, while the dry weight of leaves was greater in AprMay than in NovDec (Figure 4). That is, the leaf dry matter per unit area in AprMay was greater than that in NovDec. From the changing trend of the dry matter accumulation curve, there was unobvious difference between the dry matter accumulation in the early AprMay and NovDec, but the rapid growth period of the dry matter accumulation curve in AprMay was obviously longer than that in NovDec. Therefore, the accumulation ratio of leaf dry weight per unit area was greater in AprMay than in NovDec (Figure 5). The elongation of the petiole and internode had no significant difference in the two scenarios. Although the diameter change of the petiole and internode was slightly different during the expansion period, the final diameter was still relatively close, with no obvious difference.
3.3. Methodological Accuracy of Organ Expansion Models
The data obtained in experiment 2 was used for model verification. The relationship between the measured and simulated values of leaf length, leaf width, petiole length, and internode length is shown in Figure 6, while the bias, RMSD, and accuracy of the model simplifications are shown in Table 4. The simulated value was close to the measured value, and the accuracy of the four simulations was higher than 0.99.
(a)
(b)
(c)
(d)

3.4. Analysis of Leaf Light Interception and Photosynthetic Characteristics in 3D Canopies
Comparing the temperature and light environment conditions between AprMay and NovDec, the temperature difference was small, but the light difference was large. Since light directly participates in photosynthetic reactions during the accumulation of dry matter, there was a significant difference in yield in AprMay and NovDec (Figure 7). The yield was 0.174 kg/m^{2} per day in NovDec and it was 35.3% lower than the yield in AprMay. To further analyze the canopy photosynthesis characteristics, it is necessary to combine the canopy radiation interception model and photosynthesis model.
In this study, light interception of three canopies from 8:00 to 15:00 on 4 November 2015 was simulated by using the radiositygraphic combined model. Combining the 3D structure model and the radiositygraphic combined model, the instantaneous light interception of each cucumber leaf in the canopy was calculated. The light interception data of the leaves of 4 plants inside the canopy were selected for analysis. The blocking effect of the surroundings was fully considered. The average light interception of a single leaf is shown in Figure 8. The light interception of the middle and lower canopy leaves was significantly lower than that of the upper part, and the light of the top leaves was relatively small due to the small leaf area. The lower leaves in the 20leaf canopy have lower light interception than the other two canopies. The difference in light interception between the three upper canopy leaves was not obvious.
FvCB photosynthesis model was used in leaf assimilation rate modeling considering carbon dioxide assimilation mechanisms. The light interception of each blade calculated from the light interception simulation is input into the FvCB model as the incident light intensity (I). The FvCB model calculates the assimilation rate of Rubiscolimited (A_{c}, μ mol·CO_{2}·m^{−2}·s^{−1}) and RuBP regenerationlimited (A_{j}, μ mol·CO_{2}·m^{−2}·s^{−1}), respectively. The data showed that during the daytime in NovDec, the A_{j} value was always less than the A_{c} value. It means that the limiting factor affecting the photosynthesis rate in NovDec was light intensity. The high concentration of CO_{2} in the morning can also promote photosynthesis (Figure 9). In the morning and afternoon, light intensity was only about half of the noon, while the photosynthesis rate was about 85% of noon. High CO_{2} concentration can make up part of the effect of low light intensity on the photosynthetic rate in NovDec. There was no significant difference in the photosynthetic rate of the upper leaves of the three canopies, but affected by the light interception of the leaves, the photosynthetic rate of the lower leaves of the 20leaf canopy was significantly lower than that of the other two canopies.
4. Discussion
4.1. The Contribution of 3D Model in Simulation Research
In the light interception calculation, the consistency of the 3D model and the true structure of the canopy is a crucial factor that determines the reliability of the simulation results. At the field level, many research studies were conducted on 3D point cloud acquisition of field crop canopy based on UAV (unmanned aerial vehicle) and LiDAR [22, 23]. Based on these 3D point clouds, plant height, ground cover, and aboveground biomass were obtained [24, 25]. But these methodologies presented an underestimation of biomass and leaf area index and cannot extract the architecture parameters of individual plants from the canopy. In order to obtain more detailed 3D point clouds at the individual plant level, 3D laser scanner, 3D digitizing, and multiview imaging were used in plant and organ 3D point cloud acquisition [26, 27]. Using these methods, many specific morphological parameters of various crops have been achieved, such as leaf azimuthal angle, leaf length, and leaf area of maize [28]. Since the last few years, many 3D point cloud acquisition methods were used in plant 3D structural model construction at different scales, such as LIDAR, 3D laser scanner, 3D digitizing, and multiview imaging [29]. However, there are still many challenges in building a sufficiently accurate 3D model because plant morphologies are diverse for different species, and there is none uniform solution for all kinds of plants [28]. The parametric modeling method applied in this study is also one of the alternative methods for constructing threedimensional plants. Based on the descriptive model of organ expansion, a plant population with any number stands can be constructed. However, since the organ expansion model is a regression model, its adaptability is limited, and the simulation of the threedimensional structure model is also constrained by the model. Therefore, we separately constructed the organ expansion curve equations in different scenarios (AprMay and NovDec) to construct the canopy growth accordingly.
4.2. Light Energy on Cucumber Development
The effect of light on plant growth and development not only directly affects the photosynthetic reaction of leaves, but also affects plant morphogenesis. In most of the plant growth models, leaf area is an important parameter directly related to biomass accumulation and structural development [30]. However, the results of this study showed that leaf area is not completely positively related to dry matter. Cucumber grown in winter has larger leaf area but less dry matter, resulting in a smaller specific leaf area (SLA). Therefore, growth is more directly related to specific leaf area than to leaf area [31].
4.3. Photosynthesis Modeling
The photosynthetic capacity of leaves is affected not only by the amount of light intercepted, but also by the age of the leaves and the ambient temperature. In particular, cucumbers belong to C3 type plants, which have stronger photorespiration and growth respiration. They consume more assimilates than C4 plants to maintain growth and protect photosynthetic organs. In many research studies, various forms of equations were used to describe the response of FvCB model parameters to temperature, such as temperature sensitivity (Q10), polynomials, exponential functions, and normal distribution functions [32]. The change in temperature dependence of photosynthesis is caused by several factors such as C_{i}, J_{max}, and J_{max} to ratio, which were also reported to change with growth temperature. The response of each factor to temperature seems to differ among species, and there is no report yet in cucumber. In the research of Zhou et al. [33], it is found that leaf aging significantly affected A_{n}, , R_{dark}, , J_{max}, R_{dark}/A_{n}, and , but leaf aging had little effect on the ratio, suggesting coordination between and processes regardless of leaf age. Further studies are needed to reconfirm the leaf aging effect on the temperature sensitivity of R_{dark} and to elucidate the mechanisms involved.
4.4. Optimal Utilization of Light Energy
In a semiclosed greenhouse, the CO_{2} concentration in the morning is significantly higher than at noon and afternoon (Figure 9). The main reason is that the respiration of plants at night accumulates a lot of CO_{2}. From the photosynthesis simulation results, the high concentration of CO_{2} in the morning can obviously promote the leaf photosynthesis rate (Figure 9). At this stage, light intensity is the main factor limiting the assimilation rate. By noon, the light intensity increased, but the CO_{2} concentration decreased significantly and became the main factor limiting the rate of assimilation. Therefore, increasing light intensity in the morning or supplementing CO_{2} at noon can increase photosynthetic output. In the results of leaf light interception and photosynthesis rate analysis (Figures 8 and 10), the light interception and photosynthetic rate of the lower leaves in the 20leaf canopy were significantly lower than those in the other two canopies. The main reason is that, in the parametric modeling process, the blade azimuth model (φ_{L}) we choose will make the blade distribution ratio in each azimuth range change with the LAI. As the LAI increases, more leaves will turn orientation towards rows with better lighting conditions to intercept more light. It can be seen that when the light conditions in winter are weak, maintaining a reasonable LAI is more conducive to canopy photosynthesis. From the results of this study, the 25leaf canopy is the best canopy structure for light interception and photosynthesis. However, considering that the light interception of leaves under 20 leaves is weak and the leaf function is weak due to the effects of leaf senescence, the number of leaves in the canopy during winter should be controlled within 25 leaves.
5. Conclusions
Low temperature and low light energy in the winter of east China are important reasons for the decline in cucumber yield, especially in greenhouses without supplementary light. This study provides an important management basis for the optimal use of light in the greenhouse in winter. Using 3D modeling technology, the virtual structure model of cucumber canopy was established, and combined light interception modeling and photosynthesis modeling, the light environment in the canopy and the photosynthetic performance of leaves were quantitatively studied to obtain cucumber management strategies in winter. By using the FvCB model, Rubiscolimited (A_{c}) and RuBP regenerationlimited (A_{j}) assimilation rate of each leaf in three canopies with different leaves per plant were simulated, respectively. Based on the simulated result, we analyzed the effects of light and CO_{2} changes on the photosynthesis at different time periods and different leaf positions. Through a comparative analysis of photosynthetic characteristics in the virtual canopy, we have come up with strategies to improve the canopy photosynthetic capacity, and a suitable canopy structure that is conducive to winter canopy photosynthetic yields.
Data Availability
All data used during the study are available from the corresponding author by request (email: qiantingting@saas.sh.cn).
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Authors’ Contributions
Tingting Qian and Xiuguo Zheng contributed equally to this work.
Acknowledgments
The authors would like to thank Dr. Feng Lu for her support during the preparation of this manuscript. This work was supported by Shanghai Agriculture Applied Technology Development Program of China (no. G2015060402), National Natural Science Foundation of China (no. 61762013), and Shanghai Science and Technology Commission Project of China (no. 19DZ2281300).
References
 T. Boulard, J.C. Roy, J.B. Pouillard, H. Fatnassi, and A. Grisey, “Modelling of micrometeorology, canopy transpiration and photosynthesis in a closed greenhouse using computational fluid dynamics,” Biosystems Engineering, vol. 158, pp. 110–133, 2017. View at: Publisher Site  Google Scholar
 X. Ding, Y. Jiang, D. Hui et al., “Model simulation of cucumber yield and microclimate analysis in a semiclosed greenhouse in China,” HortScience, vol. 54, no. 3, pp. 547–554, 2019. View at: Publisher Site  Google Scholar
 L. Marcelis, A. Broekhuijsen, E. Meinen, E. Nijs, and M. Raaphorst, “Quantification of the growth response to light quantity of greenhouse grown crops,” in Proceedings of the International Symposium on Artificial Lighting in Horticulture, pp. 97–104, Lillehammer, Norway, June 2005. View at: Google Scholar
 L. Marcelis, F. Maas, and E. Heuvelink, “The latest developments in the lighting technologies in Dutch horticulture,” in Proceedings of the International ISHS Symposium on Artificial Lighting, pp. 35–42, Québec City, Canada, November 2000. View at: Google Scholar
 T. Hovi, J. Näkkilä, and R. Tahvonen, “Interlighting improves production of yearround cucumber,” Scientia Horticulturae, vol. 102, no. 3, pp. 283–294, 2004. View at: Publisher Site  Google Scholar
 X. Hao, X. Guo, X. Chen, and S. Khosla, “Interlighting in minicucumbers: interactions with overhead lighting and plant density,” Acta Horticulturae, vol. 1107, no. 1107, pp. 291–296, 2015. View at: Publisher Site  Google Scholar
 X. Hao, X. Guo, J. Lanoue et al., “A review on smart application of supplemental lighting in greenhouse fruiting vegetable production,” Acta Horticulturae, no. 1227, pp. 499–506, 2017. View at: Publisher Site  Google Scholar
 J. H. Kim, J. W. Lee, T. I. Ahn, J. H. Shin, K. S. Park, and J. E. Son, “Sweet pepper (Capsicum annuum l.) canopy photosynthesis modeling using 3D plant architecture and light raytracing,” Frontiers in Plant Science, vol. 7, p. 1321, 2016. View at: Publisher Site  Google Scholar
 T.T. Qian, S.L. Lu, C.J. Zhao, X.Y. Guo, W.L. Wen, and J.J. Du, “Heterogeneity analysis of cucumber canopy in the solar greenhouse,” Journal of Integrative Agriculture, vol. 13, no. 12, pp. 2645–2655, 2014. View at: Publisher Site  Google Scholar
 W. Wen, X. Guo, B. Li et al., “Estimating canopy gap fraction and diffuse light interception in 3D maize canopy using hierarchical hemispheres,” Agricultural and Forest Meteorology, vol. 276277, Article ID 107594, 2019. View at: Publisher Site  Google Scholar
 A. J. Burgess, R. Retkute, T. Herman, and E. H. Murchie, “Exploring relationships between canopy architecture, light distribution, and photosynthesis in contrasting rice genotypes using 3D canopy reconstruction,” Frontiers in Plant Science, vol. 8, p. 734, 2017. View at: Publisher Site  Google Scholar
 R. P. A. Perez, C. Fournier, L. CabreraBosquet et al., “Changes in the vertical distribution of leaf area enhanced light interception efficiency in maize over generations of selection,” Plant, Cell & Environment, vol. 42, no. 7, pp. 2105–2119, 2019. View at: Publisher Site  Google Scholar
 L. Tang, D. Yin, and C. Chen, “Optimal design of plant canopy based on light interception: a case study with loquat,” Frontiers in Plant Science, vol. 10, p. 364, 2019. View at: Publisher Site  Google Scholar
 G. BuckSorlin, P. H. B. de Visser, M. Henke et al., “Towards a functionalstructural plant model of cutrose: simulation of light environment, light absorption, photosynthesis and interference with the plant structure,” Annals of Botany, vol. 108, no. 6, pp. 1121–1134, 2011. View at: Publisher Site  Google Scholar
 T.W. Chen, M. Henke, P. H. B. de Visser et al., “What is the most prominent factor limiting photosynthesis in different layers of a greenhouse cucumber canopy?” Annals of Botany, vol. 114, no. 4, pp. 677–688, 2014. View at: Publisher Site  Google Scholar
 R. P. A. Perez, J. Dauzat, B. Pallas et al., “Designing oil palm architectural ideotypes for optimal light interception and carbon assimilation through a sensitivity analysis of leaf traits,” Annals of Botany, vol. 121, no. 5, pp. 909–926, 2017. View at: Publisher Site  Google Scholar
 Y. Y. Cho, S. Oh, M. M. Oh, and J. E. Son, “Estimation of individual leaf area, fresh weight, and dry weight of hydroponically grown cucumbers (Cucumis sativus l.) using leaf length, width, and SPAD value,” Scientia Horticulturae, vol. 111, no. 4, pp. 330–334, 2007. View at: Publisher Site  Google Scholar
 T. Qian, X. Zheng, X. Guo, W. Wen, J. Yang, and S. Lu, “Influence of temperature and light gradient on leaf arrangement and geometry in cucumber canopies: structural phenotyping analysis and modelling,” Information Processing in Agriculture, vol. 6, no. 2, pp. 224–232, 2019. View at: Publisher Site  Google Scholar
 S. Lu, C. Zhao, X. Guo, W. Wen, and B. Xiao, “A new paradigm for fast interactive design of crops,” Intelligent Automation and Soft Computing, vol. 16, pp. 1147–1155, 2010. View at: Google Scholar
 G. D. Farquhar, S. von Caemmerer, and J. A. Berry, “A biochemical model of photosynthetic CO_{2} assimilation in leaves of C3 species,” Planta, vol. 149, no. 1, pp. 78–90, 1980. View at: Publisher Site  Google Scholar
 K. Kahlen and H. Stützel, “Modelling photomodulated internode elongation in growing glasshouse cucumber canopies,” New Phytologist, vol. 190, no. 3, pp. 697–708, 2011. View at: Publisher Site  Google Scholar
 L. Lei, C. Qiu, Z. Li et al., “Effect of leaf occlusion on leaf area index inversion of maize using UAVLiDAR data,” Remote Sensing, vol. 11, no. 9, p. 1067, 2019. View at: Publisher Site  Google Scholar
 R. P. A. Perez, E. Costes, F. Théveny, S. Griffon, J.P. Caliman, and J. Dauzat, “3D plant model assessed by terrestrial LiDAR and hemispherical photographs: a useful tool for comparing light interception among oil palm progenies,” Agricultural and Forest Meteorology, vol. 249, pp. 250–263, 2018. View at: Publisher Site  Google Scholar
 J. A. JimenezBerni, D. M. Deery, P. RozasLarraondo et al., “High throughput determination of plant height, ground cover, and aboveground biomass in wheat with LiDAR,” Frontiers in Plant Science, vol. 9, p. 237, 2018. View at: Publisher Site  Google Scholar
 S. Madec, F. Baret, B. de Solan et al., “Highthroughput phenotyping of plant height: comparing unmanned aerial vehicles and ground LiDAR estimates,” Frontiers in Plant Science, vol. 8, p. 2002, 2017. View at: Publisher Site  Google Scholar
 B. Elnashef, S. Filin, and R. N. Lati, “Tensorbased classification and segmentation of threedimensional point clouds for organlevel plant phenotyping and growth analysis,” Computers and Electronics in Agriculture, vol. 156, pp. 51–61, 2019. View at: Publisher Site  Google Scholar
 S. Wu, W. Wen, B. Xiao et al., “An accurate skeleton extraction approach from 3D point clouds of maize plants,” Frontiers in Plant Science, vol. 10, p. 248, 2019. View at: Publisher Site  Google Scholar
 Y. Wang, W. Wen, S. Wu et al., “Maize plant phenotyping: comparing 3D laser scanning, multiview stereo reconstruction, and 3D digitizing estimates,” Remote Sensing, vol. 11, no. 1, p. 63, 2018. View at: Publisher Site  Google Scholar
 M. VázquezArellano, H. Griepentrog, D. Reiser, and D. Paraforos, “3D imaging systems for agricultural applicationsa review,” Sensors, vol. 16, no. 5, p. 618, 2016. View at: Publisher Site  Google Scholar
 L. F. M. Marcelis, E. Heuvelink, and J. Goudriaan, “Modelling biomass production and yield of horticultural crops: a review,” Scientia Horticulturae, vol. 74, no. 12, pp. 83–111, 1998. View at: Publisher Site  Google Scholar
 S. M. Weraduwage, J. Chen, F. C. Anozie, A. Morales, S. E. Weise, and T. D. Sharkey, “The relationship between leaf area growth and biomass accumulation in Arabidopsis thaliana,” Frontiers in Plant Science, vol. 6, 2015. View at: Publisher Site  Google Scholar
 K. Hikosaka, K. Ishikawa, A. Borjigidai, O. Muller, and Y. Onoda, “Temperature acclimation of photosynthesis: mechanisms involved in the changes in temperature dependence of photosynthetic rate,” Journal of Experimental Botany, vol. 57, no. 2, pp. 291–302, 2005. View at: Publisher Site  Google Scholar
 H. Zhou, M. Xu, H. Pan, and X. Yu, “Leafage effects on temperature responses of photosynthesis and respiration of an alpine oak, Quercus aquifolioides, in southwestern China,” Tree Physiology, vol. 35, no. 11, pp. 1236–1248, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Tingting Qian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.