Artiﬁcial Neural Network-Based Seedling Phenotypic Information Acquisition of Plant Factory

: This work aims to construct an artiﬁcial neural network (ANN) ant colony algorithm (ACA)-based ﬁne recognition system for plant factory seedling phenotypes. To address the problems of complexity and high delay of the plant recognition system in plant factories, ﬁrst, multiple cameras at different positions are employed to collect images of seedlings and construct 3D images. Then, the mask region convolutional neural networks (MRCNN) algorithm is adopted to analyze plant phenotypes. Finally, the optimized ACA is employed to optimize the process timing in the plant factory, thereby constructing a plant factory seedling phenotype ﬁne identiﬁcation system via ANN combined with ACA. Moreover, the model performance is analyzed. The results show that plants have four stages of phenotypes, namely, the germination stage, seedling stage, rosette stage, and heading stage. The accuracy of the germination stage reaches 97.01%, and the required test time is 5.64 s. Additionally, the optimization accuracy of the process timing sequence of the proposed model algorithm is maintained at 90.26%, and the delay and energy consumption are stabilized at 20.17 ms and 17.71, respectively, when the data volume is 6000 Mb. However, the problem of image acquisition occlusion in the process of 3D image construction still needs further study. Therefore, the constructed ANN-ACA-based ﬁne recognition system for plant seedling phenotypes can optimize the process timing in a more real-time and lower energy consumption way and provide a reference for the integrated progression of unmanned intelligent recognition systems and complete sets of equipment for plant plants in the later stage.


Background
The spread of COVID-19 has brought huge changes to people's lives. Among these changes, an increasing number of people are beginning to pay more attention to new models of vegetable production suitable for indoor cultivation. As a new type of planting method, plant factories can provide a stable supply of vegetables for urban residents. As an agricultural facility that artificially controls the cultivation environment, the plant factory planting mode can control various factors, such as indoor light, temperature, humidity, and carbon dioxide concentration, which has greatly alleviated the phenomenon of grain yield reduction caused by the reduction in urban arable land [1][2][3]. However, at present, plants in factories are still mainly cultivated by manual intervention, and there are problems such as many production links and complicated process equipment. Therefore, the intelligent improvement of plant factories has attracted the attention of many scholars.
As a highly technically integrated agricultural production system, plant factories are characterized by informatization, mechanization, and intelligence. Mechanization is the key to enhancing efficiency and reducing the cost of plant factories, but while the efficiency In plant factories, analyzing plant phenotypes is an important way to determine plant growth conditions. However, using manual methods to analyze plant phenotypes is not only slow and error-prone but also labor-intensive. Therefore, many scholars use artificial intelligence algorithms to identify plant phenotypes and other growth characteristics. Jin et al. (2019) proposed a voxel-based CNN for corn stem and leaf classification and segmentation. The ability of this method to isolate the structural components of crop phenotypes utilizing deep learning was demonstrated [10]. Wu et al. (2019) used deep learning for fast, precise, and noninvasive measurements of optimal view planning pipelines. It was found that its flexibility and planning time was significantly better than stand-alone robots [11]. Nesteruk et al. (2021) proposed an image compression mechanism for data acquisition. The results showed that it had an obvious effect on the deviation of plant classification, disease identification, and phenology [12]. Yang et al. (2021) proposed a framework for intelligent detection and identification of plant stomata under feature weight transfer learning, and the results showed that this method was remarkably superior to existing methods [13]. Lv and Qiao (2020) constructed a deep belief network-based cognitive computing system model and applied it to the control of collaborative robots and related data testing. It was proven that the system was of high accuracy and security [14].
Through the research and analysis of the above scholars, it is found that most scholars apply artificial intelligence algorithms to extract relevant data features in plant phenotype analysis. However, under the background of plant factories, it is rare to analyze the phenotype of plant seedlings, and there are still many problems in plant factories, such as many production links and complicated process equipment. Therefore, artificial intelligence algorithms are introduced to acquire plant phenotypic information and optimize process timing, which is crucial to enhance the efficiency of plant factories.

Innovation and Organizational Structure
In summary, with plant factories being used today, the efficient cultivation of plants without damage is of great practical significance to the improvement of crop productivity. Therefore, in this work, aiming at the problems of many production links and miscellaneous process equipment in plant factories, multiple cameras at different locations are employed to collect seedling images and construct 3D images. Then, the MRCNN algorithm is innovatively introduced to analyze the plant phenotype, and the optimized ACA is adopted to optimize the process timing in plant factories so that a plant seedling phenotype fine management system based on ANN-ACA is established. Finally, its performance is demonstrated through experiments to provide a basis for the subsequent intelligent development of plant factories. The overall organizational structure of our work is as follows. Chapter 1 is the introduction, which describes the related development background of plant factories and related studies on plant phenotypes and proposes the innovation points and contributions of this work. In the second chapter, the methods are expounded. Multiple cameras at various locations are established to collect seedling images and design 3D images, and a fine management system of plant factory seedling phenotypes is established based on ANN-ACA. The third chapter is the results and discussion, where the experimental evaluation of the model performance is performed, and it is compared with different schemes for further discussion and analysis. The fourth chapter is the conclusion, which gives a brief description of the research results, as well as the limitations of this research and future prospects.

Temporal Analysis of Spatial Structure and Process Management in Plant Factories
The plant factory is composed of different functional rooms, such as a cultivation workshop, nursery room, harvesting and storage room, machinery room (nutrient solution pipe, carbon dioxide cylinder, control equipment, etc.), and management room (Figure 1). From the analysis of the system structure, the plant factory is based on the peripheral structure and various functional units. The all-weather operation and intelligent management of the plant factory are ensured through the construction of subsystems such as the nutrient solution circulation and control system, multilayer 3D hydroponic cultivation system, air conditioning, and purification system, carbon dioxide gas fertilizer release system, artificial light source system, and computer automatic control system [15][16][17][18].
Multiple cameras at various locations are established to collect seedling images and design 3D images, and a fine management system of plant factory seedling phenotypes is established based on ANN-ACA. The third chapter is the results and discussion, where the experimental evaluation of the model performance is performed, and it is compared with different schemes for further discussion and analysis. The fourth chapter is the conclusion, which gives a brief description of the research results, as well as the limitations of this research and future prospects.

Temporal Analysis of Spatial Structure and Process Management in Plant Factories
The plant factory is composed of different functional rooms, such as a cultivation workshop, nursery room, harvesting and storage room, machinery room (nutrient solution pipe, carbon dioxide cylinder, control equipment, etc.), and management room (Figure 1). From the analysis of the system structure, the plant factory is based on the peripheral structure and various functional units. The all-weather operation and intelligent management of the plant factory are ensured through the construction of subsystems such as the nutrient solution circulation and control system, multilayer 3D hydroponic cultivation system, air conditioning, and purification system, carbon dioxide gas fertilizer release system, artificial light source system, and computer automatic control system [15][16][17][18]. However, there are often problems of redundant production links and complicated process equipment in the practical adoption of plant factory models. Therefore, the key links and production factors of plant production in the plant factory are analyzed, and the phenotypic information of seedlings and the optimal timing rules of multi-equipment linkage operations in the factory environment are explored to realize the automatic production process and equipment linkage integrated control identification system, such as sowing, seedling raising, transplanting, harvesting, cleaning, packing, transporting, and stacking. It aims to provide a demonstration role for the development of the integration of intelligent identification systems and complete sets of equipment in unmanned plant factories.

ANN-Based Phenotype Analysis of Plant Seedlings
To accurately analyze the plant seedling phenotype of the plant factory, the ANN algorithm is employed. An extended form algorithm under Faster R-CNN, namely, the MRCNN algorithm, is selected. The algorithm can solve the orientation problem of the target image with semantic segmentation methods [19,20]. The MRCNN frame structure is shown in Figure 2. However, there are often problems of redundant production links and complicated process equipment in the practical adoption of plant factory models. Therefore, the key links and production factors of plant production in the plant factory are analyzed, and the phenotypic information of seedlings and the optimal timing rules of multi-equipment linkage operations in the factory environment are explored to realize the automatic production process and equipment linkage integrated control identification system, such as sowing, seedling raising, transplanting, harvesting, cleaning, packing, transporting, and stacking. It aims to provide a demonstration role for the development of the integration of intelligent identification systems and complete sets of equipment in unmanned plant factories.

ANN-Based Phenotype Analysis of Plant Seedlings
To accurately analyze the plant seedling phenotype of the plant factory, the ANN algorithm is employed. An extended form algorithm under Faster R-CNN, namely, the MRCNN algorithm, is selected. The algorithm can solve the orientation problem of the target image with semantic segmentation methods [19,20]. The MRCNN frame structure is shown in Figure 2. In Figure 2, the Faster R-CNN-based MRCNN algorithm is further integrated with the region proposal network (RPN) so that it can be well applied in complex life scenarios such as seedling phenotype recognition. ROIAIign determines the feature value of each point in the region of interest of the original image through bilinear interpolation and then performs pooling and other operations to improve the accuracy, providing a new solution to the misalignment alignment problem caused by pooling direct sampling [21]. In this algorithm, the backpropagation calculation is expressed as follows.
In Equation (1) The expression of the loss function L is Equation (5).  In Figure 2, the Faster R-CNN-based MRCNN algorithm is further integrated with the region proposal network (RPN) so that it can be well applied in complex life scenarios such as seedling phenotype recognition. ROIAIign determines the feature value of each point in the region of interest of the original image through bilinear interpolation and then performs pooling and other operations to improve the accuracy, providing a new solution to the misalignment alignment problem caused by pooling direct sampling [21]. In this algorithm, the backpropagation calculation is expressed as follows.
In Equation (1), y rj refers to the j-th coordinate point of the r-th region after pooling, i(r, j) refers to the original image point coordinates corresponding to point y rj , and ∆h and ∆w refer to the gradients of the horizontal and vertical coordinates, respectively. Each grid can be defined as a bin, and the size of each bin is expressed as w k × h k . The scoring equation for judging that the feature information on the (i, j)-th bin belongs to different categories is shown in Equation (2).
In Equation (2), z i,j,c refers to k 2 (c + 1) outputs of a score map, the coordinates of the upper left corner of original region of interest are represented by (x 0 , y 0 ), n is the pixel numbers in the bin, and Θ is the network parameters. Then, Equation (3) is obtained for each region of interest (ROI).
Average score r(Θ) votes for each ROI. The MRCNN algorithm replaces the traditional SoftMax classifier with a sigmoid classifier, and the corresponding calculation equation is shown in Equation (4).
The expression of the loss function L is Equation (5).
In Equation (5), L cls refers to the loss function in logarithmic form, and the classification corresponds to two categories: target and background. L box refers to the loss function of the bounding box in the original image, and L mask is the mean binary cross-entropy loss function, whose corresponding calculation equation is as follows. Therefore, according to Equation (6), the L mask in the ROI of the i(1 ≤ i ≤ c)-th category is only related to the i-th L mask , and the competition between different categories can be effectively avoided by corresponding to a mask for each class, thus improving the phenotype segmentation of seedlings. The loss principle of the multitask loss function is shown in Figure 3.
tion of the bounding box in the original image, and mask L is the mean binary cross-entropy loss function, whose corresponding calculation equation is as follows.
Therefore, according to Equation (6), the mask L in the ROI of the ( ) egory is only related to the i-th mask L , and the competition between different categories can be effectively avoided by corresponding to a mask for each class, thus improving the phenotype segmentation of seedlings. The loss principle of the multitask loss function is shown in Figure 3. From Figure 3, the mask feature predicted to be classified as k is obtained first regarding the original image, and then the mask area surrounded by the bounding box in the original image is mapped to the mask area feature of m×m size. Finally, the average binary cross-loss entropy of the m×m region is calculated and obtained.
In the seedling phenotype analysis, the center of the sponge conveyor belt in the phenotype acquisition area is set as the coordinate origin to further obtain more abundant and accurate information ( Figure 4). The images of the seedlings are collected by using eight cameras (acA1300-30 gm, Basler, Arensburg, Germany) at different positions, and then the growth maps of the seedlings in different directions are plotted. From Figure 3, the mask feature predicted to be classified as k is obtained first regarding the original image, and then the mask area surrounded by the bounding box in the original image is mapped to the mask area feature of m × m size. Finally, the average binary cross-loss entropy of the m × m region is calculated and obtained.
In the seedling phenotype analysis, the center of the sponge conveyor belt in the phenotype acquisition area is set as the coordinate origin to further obtain more abundant and accurate information ( Figure 4). The images of the seedlings are collected by using eight cameras (acA1300-30 gm, Basler, Arensburg, Germany) at different positions, and then the growth maps of the seedlings in different directions are plotted.
The key to the acquisition of seedling images is detecting different orientations of the seedlings, which currently mainly includes monocular and binocular vision detection algorithms [22,23]. The former represents the acquisition of seedling growth images by a camera. If point C is the camera, point P is the seedling in space, and the imaging point of the seedling in space is point P'(u,v), then the angle between the projection of the camera's vertical perspective on the ground and the xand y-axes of the horizontal ground is α, β. The mapping between the camera's horizontal perspective on the ground and the Y-axis of the ground plane is θ. Finally, Equations (7)-(9) are used to obtain α, β and θ, respectively.  The key to the acquisition of seedling images is detecting different orientations of the seedlings, which currently mainly includes monocular and binocular vision detection algorithms [22,23]. The former represents the acquisition of seedling growth images by a camera. If point C is the camera, point P is the seedling in space, and the imaging point of the seedling in space is point P' (u,v), then the angle between the projection of the camera's vertical perspective on the ground and the x-and y-axes of the horizontal ground is ,  . The mapping between the camera's horizontal perspective on the ground and the Y-axis of the ground plane is  . Finally, Equations (7)-(9) are used to obtain ,  and  , respectively.
In the above equations, the distance from the camera to the ground is d, the shortest distance mapped to the ground from the vertical angle is b y , and the longest distance is Where taking xoy as the coordinate system, y is the vertical distance between the seedling and the camera, x is the horizontal distance between the seedling and the camera, and L represents the actual distance between the seedling and the camera. In the above equations, the distance from the camera to the ground is d, the shortest distance mapped to the ground from the vertical angle is y b , and the longest distance is y 1 + y b . Values of d, y 1 + y b , x 1 are collected by measuring tools. S x and S y refer to the number of rows and columns of the images taken, respectively. The images taken are grayscale images of seedlings at different stages with a resolution of 1280 × 900 pixels. Where taking xoy as the coordinate system, y is the vertical distance between the seedling and the camera, x is the horizontal distance between the seedling and the camera, and L represents the actual distance between the seedling and the camera.
In binocular vision detection, the 3D coordinates of point P(x c , y c , z c ) are shown in Equation (13): Through monocular and binocular parallax analysis, it is revealed that the obtained images and seedling information have difficulty meeting the requirements of high real-time seedling growth. Therefore, for the acquisition of seedling phenotype images, an integral projection and region growing algorithm [24] fit for image position detection is adopted to detect the position of seedling in the first frame image. Machine learning algorithm is used for detection and analysis in real time to obtain accurate seedling phenotype information more efficiently.
In the integral projection coarse localization algorithm, I(x, y) represents the pixel gray value of seedling image at point (x, y). The vertical and horizontal integral projection function V(x) and L(y) are used for calculation, as shown in Equations (14) and (15): The horizontal and vertical integral projections are performed on these plant seedling images, and the midpoint of the obtained plant's coarse positioning position is set as seed point. Accurate seedling images are harvested by finely positioning seedlings using the region growing method. Eight cameras are used to obtain grayscale images of seedlings at different stages with a resolution of 1280 × 900 pixels. Background light is used to separate the seedling outlines from the background, and the seedling outlines generated from different cameras are merged. The resulting voxel spatial resolution is 0.25 mm/voxel, and the size is 240 × 240 × 300 voxels (X, Y, Z), corresponding to the actual 3D space of 60 × 60 × 76 mm. The 3D reconstruction speed is 20-60 ms per plant (i7-3.2 GHz processor is used); that is, the 3D point cloud model of approximately 10,000 seedlings is reconstructed every hour. In the acquisition of seedling images at different stages, it is agreed that the light condition is good, the temperature is kept between 18 and 22 degrees Celsius, and the relative humidity is kept between 40% and 70%, so that the effects of light and other environmental factors are ignored in the images of seedlings at different stages. The algorithm for obtaining 3D information of plant seedlings using multiple cameras with different orientations is shown in Algorithm 1. Output: 3D data information of plant seedlings 4 Rough location of seedlings by integral projection 5 Fine location of seedlings by regional growth method 6 Collecting plant seedling image I k , k ∈ {1, 2, · · · , 8} 7 For Image of each plant seedling I k 8 B k ← foreground background segmentation (I k ) 9 for each voxel v of 3D space V 10 V(v) ← 1 11 End for 12 End for 13 For each voxel v of 3D space V 14 For each camera k 15 i ← obtain the corresponding pixel (v, P k ) 16 If End if 20 End for 21 End for 22 End After the 3D data information of the plant is harvested through the above algorithm flow, the shape of the stem and leaf of the plant seedling, petioles, and other small structural information can be visualized clearly ( Figure 5). Thus, accurate plant seedling phenotype information can be obtained.

Optimization System Analysis of Plant Factory Seedling Phenotype Fine Management via Optimized ACA
After the identification of seedling phenotypes in plant factories, it is also critical to carefully manage plant growth according to the obtained phenotypic information. Therefore, the production operation control identification system in the plant factory is equipped and optimized to develop the "environment-quality-yield" intelligent and fine identification system. Figure 5 shows the framework of the ANN-ACA-based plant factory seedling phenotype fine identification system.
In Figure5, first, the center of the sponge conveyor belt in the phenotype acquisition area is set as the coordinate origin, eight cameras at different positions are employed to collect the images of the seedlings, and a 3D image of the plant seedlings is constructed. Second, MRCNN is used to extract features in 3D images of seedlings. Among them, the mask R-CNN algorithm can complete the task of instance segmentation. Based on the input plant seedling image, it can provide the boundary box of each instance and segment the instance within the box to generate the instance mask to understand the plant seedling phenotype. Finally, the optimized ACA is used to optimize and solve the optimal timing of equipment linkage operations. Combined with real-time monitoring parameter information such as the cultivation environment and operation equipment, the operation sequence of each piece of equipment is established, and the automatic production of logistics, equipment linkage, and comprehensive configuration software are realized. In this

Optimization System Analysis of Plant Factory Seedling Phenotype Fine Management via Optimized ACA
After the identification of seedling phenotypes in plant factories, it is also critical to carefully manage plant growth according to the obtained phenotypic information. Therefore, the production operation control identification system in the plant factory is equipped and optimized to develop the "environment-quality-yield" intelligent and fine identification system. Figure 5 shows the framework of the ANN-ACA-based plant factory seedling phenotype fine identification system.
In Figure 5, first, the center of the sponge conveyor belt in the phenotype acquisition area is set as the coordinate origin, eight cameras at different positions are employed to collect the images of the seedlings, and a 3D image of the plant seedlings is constructed. Second, MRCNN is used to extract features in 3D images of seedlings. Among them, the mask R-CNN algorithm can complete the task of instance segmentation. Based on the input plant seedling image, it can provide the boundary box of each instance and segment the instance within the box to generate the instance mask to understand the plant seedling phenotype. Finally, the optimized ACA is used to optimize and solve the optimal timing of equipment linkage operations. Combined with real-time monitoring parameter information such as the cultivation environment and operation equipment, the operation sequence of each piece of equipment is established, and the automatic production of logistics, equipment linkage, and comprehensive configuration software are realized. In this way, the equipment production operation control identification system and "environment-quality-yield" intelligent fine identification system are established.
In the process of optimizing the process timing using ACA [25], the Canny operator [26] is integrated into it, and a Canny-based optimized ACA is constructed. The steps of the algorithm are as follows. First, it acquires the edge heuristic information η. The improved ACA linearly combines the high and low thresholds to obtain the final threshold Th, as shown in Equation (16).
In Equation (16), a and b are constants, thL refers to the low threshold, and thH refers to the high threshold. Usually, thH defaults to 0, and thL defaults to 0.7. Then, the heuristic information η(i, j) of the pixel point (i, j) has a value range, as shown in Equation (17).
The point is an edge point Second, the edge tracking model of the Canny-based optimized ACA is established, and the transition probability is defined as follows according to the traditional ant colony edge detection algorithm.
In Equation (18), α, β refers to the control factors of pheromone and heuristic information in the process of moving and is set to α = 0.6, β = 0.2. τ i (t) refers to the size of the pheromone on pixel i at the current time t. S represents the neighborhood range of the current pixel, and the transition probability is calculated in the neighborhood of S. To reduce the amount of calculation, the following equation is used: In the setting of loop termination conditions, Width and Height are, respectively, the width and height of the image, and then the rate of steps that each ant is limited to walk is illustrated in Equation (20).
In Equation (20), dis_rate is a constant in the (0,1] interval. The path pheromone is further updated to obtain an updated objective function, as shown in Equation (21): In Equation (21), ρ ∈ [0, 1] is the pheromone volatilization rate, and ∆τ refers to the total amount of pheromone released by all ants at this pixel point in this cycle. The displacement matrix Q is used to cycle shift the unique seedling sample image in the target region, and the training set {Q i x|i = 0, 1, . . . , n − 1} is obtained. The replacement matrix Q is shown in Equation (22): Agriculture 2023, 13, 888 10 of 17 Canny-based optimized ACA is employed to calculate the optimal timing of the process. The specific optimization process is as follows. The target loss function L is the weighted sum of the confidence error L con f and the position error L loc .
L(x, c, l, g) = 1 N L con f (x, c) + αL loc (x, l, g) In Equation (23), c is the detection class confidence of target seedlings, l the prediction frame parameter, g the real frame parameter, and N the number of matched real frames, and α can adjust the ratio between the two types of losses, whose default value is 1. SoftMax loss is used to calculate the confidence error L con f .
In Equation (24), if the i-th predicted frame has a matching relationship with the j-th real frame and the matching category is p, x p ij equals 1; otherwise, it is 0, and confidenceĉ p j is presented in Equation (25): The calculation of smooth L1 (x) in Equation (25) is presented in Equation (26).
The process of training the classifier of this model is further analyzed, as described in Equation (27): That is, it finds the optimal ω so as to lowest the cost function.
In Equation (28), a i corresponds to the training sample coefficient, ϕ(X i ) is the function, which is utilized to map the training sample x to a high-dimensional feature space, and Equation (29) shows the relationship between x and x* in the high-dimensional feature space.
In Equation (31),â is the discrete Fourier transform of a, and k xx is the first-row vector of the kernel matrix K. Thus, the classifier training changes from finding the best ω to finding the optimal a.

Experiment and Evaluation Analysis
To evaluate the performance of the ANN-ACA-based plant factory seedling phenotype fine identification system, the system is generated by MATLAB simulation. Then, the vegetable lettuce in the plant factory is taken as the object, and each seedling is placed separately in rows and columns in the seedling cultivation area of the plant. Three-dimensional images of plants are obtained by using the system constructed in this work. The detailed construction details of 3D images are described in Section 2.2 above, and 408 images of 51 lettuce plants are finally obtained. The acquired image data are randomly distributed by 8:2, of which 80% are utilized for training and 20% for testing. These training samples are trained 500 times with 4 images each time, that is, 2000 iterations are carried out, that is, the model begins to converge, the loss of the verification sample begins to rise, and the loss of the training sample is fixed in a certain area. The PyTorch training framework, Intel Core i7 processor, 11 GB video memory, 1080 Ti GPU, and 16 GB memory computer are used for training.
The proposed model algorithm is compared with R-CNN [27], SPPNet [28], Fast R-CNN [29], Faster R-CNN [30], and MRCNN proposed by others. The phenotypes of plants are analyzed for accuracy according to four stages, including the germination stage, seedling stage, rosette stage, and heading stage, and the time required for the training and test data set is analyzed.
For the optimization evaluation of the process timing, the optimized ACA proposed is analyzed in contrast to the genetic algorithm (GA) [31], ACA [32], particle swarm optimization (PSO) [33], and methodology of Hou et al. (2022) [34].

System Performance Analysis
To reduce model losses, various optimizers and learning rates (LRs) are screened to optimize system performance. The selected optimizers are the standard gradient descent (GD), stochastic GD (SGD), batch GD (BGD), RMSprop algorithm, and Adam algorithm. For the LR, cosine annealing attenuations LR, exponential attenuations LR, fixed-step attenuations LR, multistep attenuations LR, and fixed value LR = 0.0001 are selected for comparison. The training and test results are illustrated in Figure 6. Figure 6a,b illustrate the optimization effects of various optimizers. Figure 6a is the loss trend of the training set, and Figure 6b is that of the test set. The training loss reaches the minimum when the RMSprop optimizer is used, followed by Adam, while the loss values of GD, BGD, and SGD are basically the same, but SGD reaches the convergence loss value the fastest. Although the RMSprop optimizer achieves an optimal training loss, the test loss is minimized when using the Adam optimizer. Therefore, the Adam optimizer has the best effect when constructing the model in this work. Figure 6c,d show the results of training loss and validation loss using iterative methods with different learning rates, respectively. Both training loss and test loss fluctuate and decrease. However, when the exponential decay learning rate iteration method is used, both training loss and test loss reach the lowest value faster than when other learning rate iteration methods are used, and the loss is minimal. Therefore, the exponential decay learning rate iteration method is chosen.

Comparison of Phenotype Recognition Accuracy of Seedlings in Plant Factories under Different Algorithms
Each algorithm is used to accurately identify and analyze the phenotypes of the seedlings according to the four stages: germination stage, seedling stage, rosette stage, and heading stage (Figure 7).   Figure 6a is the loss trend of the training set, and Figure 6b is that of the test set. The training loss reaches the minimum when the RMSprop optimizer is used, followed by Adam, while the loss values of GD, BGD, and SGD are basically the same, but SGD reaches the convergence loss value the fastest. Although the RMSprop optimizer achieves an optimal training loss, the test loss is minimized when using the Adam optimizer. Therefore, the Adam optimizer has the best effect when constructing the model in this work. Figure 6c,d show the results of training loss and validation loss using iterative methods with different learning rates, respectively. Both training loss and test loss fluctuate and decrease. However, when the exponential decay learning rate iteration method is used, both training loss and test loss reach the lowest value faster than when other learning rate iteration methods are used, and the loss is minimal. Therefore, the exponential decay learning rate iteration method is chosen.

Comparison of Phenotype Recognition Accuracy of Seedlings in Plant Factories under Different Algorithms
Each algorithm is used to accurately identify and analyze the phenotypes of the seedlings according to the four stages: germination stage, seedling stage, rosette stage, and heading stage (Figure 7). Accuracy, precision, recall, and F1 of this system model and other algorithms in analyzing the phenotypes of the four stages of plant germination stage, seedling stage, rosette stage, and heading stage are compared (Figure 7). Among the four stages, the phenotypic recognition accuracy of the germination stage is significantly better than that of the other stages, and it may be that the identification of plants from scratch is the most significant. From the comparison of different algorithms, the proposed model algorithm has a recognition accuracy of 97.01% in the germination stage, which is 1.55% superior to the model algorithms proposed by others. Furthermore, comparisons show that the precision, recall, and F1 of this model algorithm in the germination stage are 88.42%, 75.57%, Accuracy, precision, recall, and F1 of this system model and other algorithms in analyzing the phenotypes of the four stages of plant germination stage, seedling stage, rosette stage, and heading stage are compared (Figure 7). Among the four stages, the phenotypic recognition accuracy of the germination stage is significantly better than that of the other stages, and it may be that the identification of plants from scratch is the most significant. From the comparison of different algorithms, the proposed model algorithm has a recognition accuracy of 97.01% in the germination stage, which is 1.55% superior to the model algorithms proposed by others. Furthermore, comparisons show that the precision, recall, and F1 of this model algorithm in the germination stage are 88.42%, 75.57%, and 72.13%, respectively, which are significantly higher than those of the other algorithms. Hence, the phenotype recognition and prediction accuracy of the ANN-ACA-based plant factory seedling phenotype fine identification system is superior.
The comparison of time required for phenotype recognition of each algorithm is illustrated in Figure 8.

Process Optimization and Data Transmission Performance Analysis of Plant Factories under Different Algorithms
The comparison results of plant factory process optimization and data transmission performance under different algorithms are shown in Figures 9 and 10, respectively.   With the increase in iterations, a trend of decreasing first and then basically remaining unchanged is demonstrated both in the required training duration and test duration, that is, reaching convergence. The training duration and test duration of the proposed model are stable at 12.26 s and 5.64 s, respectively. Moreover, the model algorithm constructed in this work requires markedly less time for prediction than other models. Therefore, the establishment of the ANN-ACA-based plant factory seedling phenotype fine identification system can achieve a high phenotype recognition and prediction effect in a shorter time.

Process Optimization and Data Transmission Performance Analysis of Plant Factories under Different Algorithms
The comparison results of plant factory process optimization and data transmission performance under different algorithms are shown in Figures 9 and 10, respectively.

Process Optimization and Data Transmission Performance Analysis of Plant Factories under Different Algorithms
The comparison results of plant factory process optimization and data transmission performance under different algorithms are shown in Figures 9 and 10, respectively.    The time delay and energy consumption analysis of the data transmission of the three algorithms are expounded in Figure 10. With an increased data amount, the required energy consumption and delay both show an upward trend. The delay and energy consumption of the proposed model algorithm is stable at 20.17 ms and 17.71 ms, respectively, when the data are 6000 Mb, which is notably inferior to the delay and energy consumption of other models. The average time delay and energy consumption of each algorithm are in the order of the proposed model algorithm < Hou et al. (2022) < ACA < GA < PSO from small to large. Therefore, the constructed ANN-ACA-based plant factory seedling phenotype fine identification system can optimize the process timing in a more real-time manner with lower energy consumption.

Discussion
In this work, a plant factory seedling phenotype fine management system based on ANN-ACA is constructed. In the performance analysis of the system, ablation experiments are conducted from the perspective of optimizer selection and learning rate iteration. The system demonstrates the best performance when the Adam optimizer and exponential decay learning rate iteration are used. This is consistent with the views of Xu et al. (2021) [35]. In the analysis of seedling phenotype recognition accuracy at different stages, it is found that the phenotype recognition precision in the germination stage is significantly better than that in other stages, which might be the most significant in plant-fromscratch recognition. Moreover, the precision, recall, and F1 of the germination stage are 88.42%, 75.57%, and 72.13%, respectively, drastically superior to those of the other algorithms. However, the identification accuracy of the seedling stage, rosette stage, and pelleting stage was lower than that of the germination stage, which is consistent with the results obtained by Hati and Singh (2021) [36]. Further comparative analysis of plant process optimization and data transmission performance under different algorithms shows that the process timing optimization accuracy of the proposed model algorithm is maintained at 90.26%, and the time delay and energy consumption are stabilized at 20.17 ms and 17.71, respectively, when the data volume is 6000 Mb. This further verifies the accuracy of the proposed algorithm in data transmission rate and process timing optimization, which is similar to the study of Zhou et al. (2022) [37]. Thus, the model constructed in this work can effectively identify the phenotype of seedlings in plant factories and carry out fine management of their processes.

Conclusions
Today, the application field of artificial intelligence technology is increasingly extensive. Aiming at the problems of confused identification of plant phenotype, many production links, and complicated process equipment in the existing plant factory, the MRCNN algorithm is introduced to analyze plant phenotype, and the process sequence in the plant  The time delay and energy consumption analysis of the data transmission of the three algorithms are expounded in Figure 10. With an increased data amount, the required energy consumption and delay both show an upward trend. The delay and energy consumption of the proposed model algorithm is stable at 20.17 ms and 17.71 ms, respectively, when the data are 6000 Mb, which is notably inferior to the delay and energy consumption of other models. The average time delay and energy consumption of each algorithm are in the order of the proposed model algorithm < Hou et al. (2022) < ACA < GA < PSO from small to large. Therefore, the constructed ANN-ACA-based plant factory seedling phenotype fine identification system can optimize the process timing in a more real-time manner with lower energy consumption.

Discussion
In this work, a plant factory seedling phenotype fine management system based on ANN-ACA is constructed. In the performance analysis of the system, ablation experiments are conducted from the perspective of optimizer selection and learning rate iteration. The system demonstrates the best performance when the Adam optimizer and exponential decay learning rate iteration are used. This is consistent with the views of Xu et al. (2021) [35]. In the analysis of seedling phenotype recognition accuracy at different stages, it is found that the phenotype recognition precision in the germination stage is significantly better than that in other stages, which might be the most significant in plant-from-scratch recognition. Moreover, the precision, recall, and F1 of the germination stage are 88.42%, 75.57%, and 72.13%, respectively, drastically superior to those of the other algorithms. However, the identification accuracy of the seedling stage, rosette stage, and pelleting stage was lower than that of the germination stage, which is consistent with the results obtained by Hati and Singh (2021) [36]. Further comparative analysis of plant process optimization and data transmission performance under different algorithms shows that the process timing optimization accuracy of the proposed model algorithm is maintained at 90.26%, and the time delay and energy consumption are stabilized at 20.17 ms and 17.71, respectively, when the data volume is 6000 Mb. This further verifies the accuracy of the proposed algorithm in data transmission rate and process timing optimization, which is similar to the study of Zhou et al. (2022) [37]. Thus, the model constructed in this work can effectively identify the phenotype of seedlings in plant factories and carry out fine management of their processes.

Conclusions
Today, the application field of artificial intelligence technology is increasingly extensive. Aiming at the problems of confused identification of plant phenotype, many production links, and complicated process equipment in the existing plant factory, the MRCNN algorithm is introduced to analyze plant phenotype, and the process sequence in the plant factory is optimized through the optimized ACA, thus constructing a fine management system of plant factory seedling phenotype based on ANN-ACA. Experimental analysis suggests that this research algorithm can accurately identify the plant phenotype, and the recognition accuracy reaches 97.01% and can optimize the relationship of the plant process sequence according to the phenotype, which can provide an example for the progression of intelligent management systems and complete equipment integration in unmanned plant factories in the later stage. However, there are some shortcomings in this work, such as identifying plant phenotypic information efficiently by constructing 3D images, in which it is assumed that a point on a plant is visible in all camera views. However, in practical problems, there may be some plants that cannot be photographed by occlusion cameras because many plants share the same shooting system, so the solution to the occlusion problem can be studied later. At the same time, in later work, a separation machine can be designed for different kinds of plants, which makes the culture system of different seedlings more focused. In addition, the MRCNN algorithm and optimized ACA are introduced to build the model system, but the influence and combination of the two algorithms are not considered. Therefore, the combination of the MRCNN algorithm and ant colony optimization algorithm can be further promoted in follow-up work, and the performance of the seedling phenotype fine management system in plant factories can be further improved.