Literature DB >> 35130303

Algorithm for wireless sensor networks in ginseng field in precision agriculture.

Changcheng Li1,2,3, Deyun Chen1, Chengjun Xie4, You Tang2,3.   

Abstract

In the research on energy-efficient networking methods for precision agriculture, a hot topic is the energy issue of sensing nodes for individual wireless sensor networks. The sensing nodes of the wireless sensor network should be enabled to provide better services with limited energy to support wide-range and multi-scenario acquisition and transmission of three-dimensional crop information. Further, the life cycle of the sensing nodes should be maximized under limited energy. The transmission direction and node power consumption are considered, and the forward and high-energy nodes are selected as the preferred cluster heads or data-forwarding nodes. Taking the cropland cultivation of ginseng as the background, we put forward a particle swarm optimization-based networking algorithm for wireless sensor networks with excellent performance. This algorithm can be used for precision agriculture and achieve optimal equipment configuration in a network under limited energy, while ensuring reliable communication in the network. The node scale is configured as 50 to 300 nodes in the range of 500 × 500 m2, and simulated testing is conducted with the LEACH, BCDCP, and ECHERP routing protocols. Compared with the existing LEACH, BCDCP, and ECHERP routing protocols, the proposed networking method can achieve the network lifetime prolongation and mitigate the decreased degree and decreasing trend of the distance between the sensing nodes and center nodes of the sensor network, which results in a longer network life cycle and stronger environment suitability. It is an effective method that improves the sensing node lifetime for a wireless sensor network applied to cropland cultivation of ginseng.

Entities:  

Mesh:

Year:  2022        PMID: 35130303      PMCID: PMC8820603          DOI: 10.1371/journal.pone.0263401

Source DB:  PubMed          Journal:  PLoS One        ISSN: 1932-6203            Impact factor:   3.240


Introduction

Precision agriculture, a modern agricultural management strategy and operational technology system based on spatial information management and mutation analysis, is of great significance for improving the efficiency of utilization of agricultural resources and ensuring sustainable agricultural development. In a precision agriculture system, numerous monitoring data need to be collected from the ginseng cultivation fields, including temperature °C, humidity %RH, ground temperature °C, ground humidity %RH, illumination Lux, CO₂ ppm, wind direction velocity m/s, and precipitation mm. The monitoring environment is diverse and complex. Besides, changes in channel conditions resulting from factors like occlusion are also inevitable during the ginseng growth. Hence, in precision agriculture systems, numerous scholars have focused on improving the energy utilization and network life-cycle in a complex communication context under limited energy conditions of sensor equipment. In broad terms, precision agriculture refers to an agricultural endeavor that aims to deliver management at an adequate time and location with the appropriate potency. It practices regulated inputs by partitioning the croplands into management blocks or zones rather than offering uniform management. Given the sustained or elevated yield at lower inputs, precision agriculture should mitigate environmental impacts and improve income for the peasants in theory [1]. To achieve precision agriculture, it is necessary to understand the spatial changes in crop and edaphic parameters within croplands [2]. Data communication which serves as an intermediate link between data collection and platform application in precision agriculture, plays a vital role in the stable operation of the system [3]. However, since the sampling points of agricultural information are highly scattered, their geographical distribution is uneven, and power supply is difficult, the construction of wired networks is difficult. Thus, the wireless sensor network has become a basic communication technology in precision agriculture due to its low construction cost, small size, convenient deployment, and flexible network [4, 5]. Wireless Sensor Networks (WSNs), capable of supporting diverse sectors, including manufacturing, education, agriculture, and environmental surveillance, have permanently changed all lives. However, the wireless presence disallows the replacement of the sensor node cells in case the deployment site is remote or unsupervised. To this end, a few studies have attempted to prolong the life expectancy of the nodes. Although cluster-based routing has greatly helped solve this problem, the cluster head (CH) selection can still be bettered through critical variable introduction [6]. However, to improve the flexibility of the network deployment, most devices in the wireless sensor networks use batteries as their main energy source [7]. B. Chatterjee et al. To minimize the network energy expenditure, a co-optimization strategy (real-time) has been formulated, achieving total life-cycle maximization for the batteries. The communication–computation energy tradeoffs for a meshed WSN structure are described by Chatterjee et al. [8], which is arranged throughout a college campus2400-acre in size to monitor the water nitrate content and temperature humidity for smart agricultural purposes via multiple sensors. Deng et al. develop, implement, and elucidate a hardware platform suitable for nodes of a self-powered WSN and mainly aims to prolong the service life of WSN nodes configured in the field context by offering a system for capturing hybrid energy. Apart from minimizing power expenditure through optimization of sensors, microcontrollers, radio frequency (RF) transceivers, etc., the energy sources also need to be taken into consideration, rather than habitual practices of cell replacements or recharges [9]. Exploiting the TENG-based direct sensory transmission (TDST), F. Wen et al. put forward a short-range self-powered (SS) WSN that is free of cells. Wen et al. [10] improved the output of TENG signals and strengthened their frequency through instantaneous short-term discharge by manipulating either a mechanical switch or a diode/switch combo. In this way, signals can be transmitted directly, where no external power supplier or extra wireless module is needed. For energy-harvesting underwater optical WSNs (EH-UOWSNs), N. Saeed et al. develops a localization strategy that is based on received signal strength (RSS), which though encounters great difficulty of predicting ranges due to the optical noise sources in optical devices and the seawater-induced disruption of channels. For viable nodes modeled by utilizing the features of optical communication channels (underwater), their RSSs are determined to achieve the localization of network [11]. In addition to satisfying the basic collection tasks, the device uses a large part of its energy for transmitting information to other devices. Therefore, in the agricultural environment, the topics of improving the energy utilization efficiency of the sensor devices in the communication process and improving the network life cycle have been rigorously researched by many scholars. In view of the complex and changeable wireless environment in the precision agriculture systems, this paper proposes a novel energy-efficient networking algorithm by dynamically adjusting the network architecture and the signal transmitting power of each node based on the particle swarm optimization (PSO) algorithm, which fully considers the residual energy variation of nodes during communication. Compared to other algorithms, the present method extends the WSN life-cycle greatly. The remainder of this work is arranged as follows: Section II introduces a precision agriculture system with ginseng cultivation as an example. Section III describes the currently available networking algorithms for the nodes of wireless sensor networks in precision agriculture systems. Section IV describes the node model of the wireless sensor network and introduces the networking algorithm corresponding to the energy-saving wireless sensor nodes based on power control. Results of the validation performed on the algorithm have been presented in Section V by comparing them with the corresponding results obtained using the existing protocols. Section VI provides a summary of this study.

The precision agriculture system for ginseng

As a crucial ancient Chinese herbal material, Panax ginseng exerts prominent functions on a variety of diseases [12]. As an important Chinese medicinal material, ginseng has the functions of recovering meridians and astringing consciousness loss, replenishing the spleen and kidney, and calming the spirit with marked effects for relieving nerves and improving the immune system [13]. The main ginseng-producing areas are distributed in Jilin, Liaoning, and Heilongjiang in China and South Korea, and North Korea. As a perennial herb, Panax ginseng prefers a cold and sunny climate, and ginseng growth requires strict illumination conditions. In addition, in terms of its growth environment, the relative water content of the soil needs to be above 70–80%, and the pH value required is preferably 6.0–6.5. In order to ensure the sustainability of ginseng cultivation, China has undergone a major transformation from “cutting forest to grow ginseng” to “cultivating ginseng in crop fields”. Ensuring that the trace elements are available in sufficient quantities in the soil and the growth environment in the farmland meets the requirements for ginseng cultivation are the key aspects that impact the emergence rate and yield of ginseng. As the Internet of Things (IoT) and novel sensing technologies are evolving persistently, a good solution has been provided for ginseng cultivation in crop fields. As shown in Fig 1, the precision agriculture system consists of a sensing layer, a transport layer, as well as a platform for data application. The sensing layer consists of a series of sensing devices and mainly implements the collection of key information such as the soil moisture, its pH value, and the light intensity during growing ginseng in the crop fields. The transport layer is the main transmission path for connecting the collected data to the application platform with data forwarding and intermediate data processing functions. The application platform reprocesses the data according to actual needs and provides a series of targeted and professional functions of data application such as environmental monitoring systems and expert systems in agriculture.
Fig 1

Structure of a precision agriculture system.

In the precision agriculture system for ginseng cultivation in a farm field, the influence of the cristate leaf on the transmission of wireless signals should not be ignored during its growing process due to the diversity of the collection parameters such as light, soil moisture, and pH value, and the deployment of sensors at different heights. Furthermore, the limited energy available to the devices in the wireless sensor network demands a set of efficient transmission strategies. In order to ensure the long-term durability of the devices under different pressures, extensive scholarly efforts have been devoted to the energy conservation in these devices from the aspects of hardware design, network routing, and data with compressive sensing. This study presents an energy-saving networking algorithm for a precision agriculture system that combines the characteristics of the wireless sensor devices and the advantages of the particle swarm optimization technique.

Related works

In the research studies conducted on the routing algorithms of wireless sensor networks based on precision agriculture systems, based on the concept of a mobile node, Prasad R et al. used miniature aerial vehicles as the receiver of regional data [14]. This reduces the energy consumption caused by multi-hop data transmission and ensures reliable network communication with low power consumption. However, periodic small-scale collection systems, such as ginseng growing, undoubtedly increase the equipment cost. In order to achieve balanced energy consumption by the network and to avoid the “premature death” of nodes caused by the “hot spot” problems, a scheme called energy-efficient load-balanced clustering (EELBC) has been formulated, which upgrades the network energy efficiency through number equilibration of cluster members [15]. Considering the shortcomings of the low-power adaptive clustering hierarchy (LEACH) protocol in the random cluster head mechanism. Kandris et al. have modeled the network as a linear system [16] and, starting from the current and future residual energy of the nodes, have proposed the equalized cluster head election routing protocol (ECHERP) that effectively improves the overall life cycle of the network. A differential evolution (DE)-based clustering strategy has been proposed [17] for WSNs, which inhibits the rapid channel (high-load) death to extend the network life-cycle. However, the wireless presence disallows the replacement of the sensor node cells in case the deployment site is remote or unsupervised. To this end, a few studies have attempted to prolong the life expectancy of the nodes. Although cluster-based routing has greatly helped solve this problem, the cluster head (CH) selection can still be upgraded through the critical variable introduction. Moreover, the main emphasis has been laid on the choice of CH or the inter-node transmission of data. For network efficiency optimization, the potential of meta-heuristic strategies can be bright [18]. Furthermore, the layout of sensors has been explored under restricted conditions of sensor and App communication ranges [19]. In Double-tiered WSNs, the algorithms are superior to the established clustering approaches, including Agglomerative Clustering, Minimum Energy Routing, PSO, Divisive Clustering, as well as Relay Node Placement and its Improved approach regarding mean performance [20]. On surface networks, the paths of varying homotopy types are calculated by developing an approach that is free of GPS and tessellation. According to the virtual planar coordinate calculations of nodes, a data package is sent by a source node to its destination in a greedy manner. Upon failure of an existing path, other greedy path options (varying homotopy types) are available for the source node to achieve the package delivery. The proposed distributed algorithms are scalable to the size of surface networks and their genus number. With the proposed design, a novel scheme of multichain routing can be generated for data package transmission in the energy-balanced WSNs. The multichain routing technique can compute the optimal inter-node distance of transmission through mathematical modeling. Regarding the design standard for optimal distance, the farthest sending node within a communication scope and the direction to sink are chosen so that the multichain architectures can be established with precise direction, distance, and multipath. Shorter chain routing allows cutting unwarranted energy consumption, thereby achieving life-cycle extension for the global network [21]. Initially, a framework for smart agriculture is designed based on SWIPT. The next step involves the exploration of an energy performance-optimizing strategy for green communication. In this process, the pairing of subcarriers is optimized together with their power allocation. A communication procedure comprises two stages. To tackle the optimization problem proposed, a Lagrangian dual formulation is utilized to put forward an efficient algorithm of iterative optimization [22]. A composite routing measurement method has been proposed for combining energy sensing, reliability awareness, robustness awareness, and resource awareness, which increases the adaptability of the system to different scenarios [23]. Azharuddin et al. [24] adopted distributed clustering and routing ideas from the perspective of fault tolerance to solve the problem of network reorganization after the failure of cluster head nodes. Under a proposed model, a brand-new problem of sensor clustering is formulated, where the sensors are divided into varying clusters, and a sensor enabling transmission quantity maximization of redundant suppressed data is chosen for each cluster [25]. Li et al. [26] develop an effective routing scheme called DMARL for UOWSNs, which exploits multiagent reinforcement learning. Initially, the network is created as a distributed multiagent architecture through modeling. The design of the routing scheme takes into account the link quality and residual energy to better accommodate the dynamic environment and help better extend the network lifespan. For faster convergence of the algorithm for reinforcement learning, two additional optimization measures are developed as well. For the distributed system, they offer a reward mechanism on this basis. Cao et al. [27] put forward enhanced versions of distributed parallel algorithms after checking the prior knowledge-based and the stochastic grouping and an additional optimizer. Their algorithms are superior to counterparts in terms of performance, which exhibit remarkably shorter time expenditure than the serial algorithms owing to the distributed parallelism nature. Hence, within a considerably restricted time, a higher performance can be attained by their algorithms. An optimization problem is formulated by Baek et al. [28] to maximize the post-data transmission minimum residual sensor energy in UAV routing (energy-saving scheme), which depends on data collection and is constrained by the travel distance of UAV. Focusing on the life-cycle maximization of a WSN, Zhang et al. [29] develop a sensor configuration scheme called EDAGD, which comprises several algorithms. Its performance is superior to the classic random deployment algorithms. Künzel et al. [30] propose a Q-learning reliable routing strategy, which has a weighting agent for achieving weight adjustment of a most-advanced graph-routing algorithm. The sets of weights are denoted by the agent state, and in the process of network operation, the weights are altered by actions. Upon shortening of mean network latency or extension of a network life expectancy, the agent is rewarded accordingly. Liu et al. [31] develop TPE-FTED, a fault-tolerant algorithm for event region detection, where the faulty nodes are recognized based on the extraction of trajectory patterns. By learning the probabilistic model online, the sensing value distributions in varying sensing states are elucidated by individual nodes. A trajectory can be generated with a particular group of probabilistic models, which is indicative of the occurrence of a special event. TPE-FTED performs the pattern matching and examines the spatiotemporal restraints by utilizing the implicit information from the trajectories, so that the faulty node declaration can be recognized. For efficiency enhancement of network having QoS parameters (multi-constrained), the problem is modeled as a multi constrained optimal path issue and preserved by developing an algorithm based on distributed learning automaton (DLA). Exploiting the pros of DLA, the proposed method finds the minimum number of nodes to conform to the QoS requirements. A few QoS routing constraints are considered during path selection, such as end-to-end delay and reliability. Algorithm outperforms existing advanced algorithms regarding energy efficiency and end-to-end delay [32]. Given above, the use of UAVs for data acquisition in dense WSNs is proposed, where the compressive data gathering (CDG)based on projections serves as a new solution technique. Data en-route are gathered from a big sensor node set by CDG to choose the projection nodes working as cluster heads (CHs), so that the number of required transmissions can be lowered to remarkably save energy and prolong the network life-cycle [33]. By fusing analysis of edge data with that of cloud data, a brand-new method for automated detection of anomalies in heterogeneous WSNs is developed. The edge data analysis exploits an entirely unsupervised algorithm of artificial neural network, while with the cloud data analysis, the multi-parameterized edit distance algorithm is utilized [34]. However, none of the abovementioned algorithms take the transmit power tunability of wireless sensor devices into consideration, resulting in insufficient utilization of the device performance during network construction. In this study, taking the characteristics of power control and wireless sensor networks into consideration, the network has been established from three aspects, namely, device performance, energy consumption, and location information. When wireless communication is inadequate, it is necessary to increase the transmission power of the high-energy devices to reduce the unnecessary data forwarding nodes. On the other hand, when wireless communication is adequate, it is necessary to decrease the transmission power of the neighboring devices for reducing the overall increased energy consumption in the network caused by the data loss or retransmission resulting from the interference of the adjacent communication channels.

The proposed algorithm

This section describes the algorithm proposed herein for an energy-efficient WSN that is applied to ginseng cultivation in the precision agriculture system. The algorithm incorporates power control into the network construction by adjusting the communication radius of each sensor node to achieve balanced energy consumption in the network. In addition, the proposed algorithm preferentially selects the forward and high-energy nodes as the cluster heads or the data forwarding nodes, so that the transmission orientation and energy expenditure of the nodes can be taken into account.

Abstraction of the network model

Fig 2 shows a schematic of the network structure of the precision agriculture system based on the clustering model. Devices in the network can be divided into three types: central node, cluster head node, and cluster member node. The central nodes, having a fixed power supply and strong information processing capabilities, provide data interaction between the wireless sensor network and the internet. The same hardware structure is employed for the cluster head nodes and the cluster member nodes in the network with the battery power statically distributed in a certain area. The cluster member node is mainly responsible for responding to the collection instruction and sending the response information to the corresponding cluster head node. In addition to the abovementioned functions, the cluster head node is also responsible for the functions of member management and data forwarding within the cluster.
Fig 2

Network diagram of the precision agriculture system.

Energy-based model for the network

For energy expenditure analysis of the wireless sensor nodes in the WSN, the widely acknowledged energy expenditure model for wireless communication is adopted in this study [16]. Where Etx and Erx represent the energy expenditures from transmission and reception of data, respectively, l is the length of the transmitted data, Eelec denotes the energy expenditure by the node for each bit of forwarded/received data, and d denotes the distance between sending and receiving nodes. In case the distance is below the threshold d, we use a free space model; otherwise, we use a multipath model. ε and ε are the power magnifications in the different models. d can be calculated using the following formula: Because E depends on the hardware, such as a digital coder, modulator, and filter, the default value of this system is fixed in this study.

Description of the proposed algorithm

In order to facilitate research and analysis in this study, a wireless sensor network has been abstracted with a two-way communication function as an undirected connectivity graph, G = (V, E), where V represents a set of all wireless sensor devices, and E is a possible set of communication paths between all devices. The transmit power can be altered for all wireless sensor devices by changing the communication radius. For a unified analysis of energy consumption in the network, the network life cycle is defined as the ratio of the number of transfers from the commencement of data transmission to the “death” of the first network node. In order to quantitatively represent the transmit power of the device, the “power adjustment factor” for each wireless sensor device is set according to the following expression: And its value is in the range of [0.5, 1.5]. With the introduction of the “power adjustment coefficient”, the wireless network construction becomes a non-deterministic polynomial problem. In order to achieve rapid optimization of the network structure, particle swarm optimization has been introduced as the main solution algorithm. Fig 3 illustrates the flowchart of our algorithm.
Fig 3

Flowchart of the proposed algorithm.

Phase I: Initialization of the proposed algorithm

Suppose the number of particles in the PSO algorithm is m, the maximum number of cycles is p, and the network optimization index is f = max(Life), where Life is the value of the life cycle of the network. The position update formula of each particle is given by Where W is the “power adjustment factor” matrix, W = {w |i = 1,…,N}, N is the number of nodes in the sensor network, λ is the inertia coefficient, C1, and C2 are the self-learning factors, and W and W represent the individual optimal value and the group optimal value, respectively.

Phase II: Network construction under particle condition

(1) Cluster head selection phase Step 1: The central node sends a request for route creation to the nodes within the communication range in a broadcast form, and each node obtains the information of the minimum hops, h, between the central node and itself. Step 2: The communication radius, R, of each node is calculated according to the “power adjustment factor” matrix, and the degree of discretization, DisRatio, among a node and its adjacent nodes is obtained as follows: where dis(j, CHi)denotes the communication distance between node i and its adjacent node j, and N denotes the number of neighboring nodes. The smaller the index, the higher is the concentration of the neighboring nodes and the cluster heads. Step 3: Within the current communication radius, the maximum number of data transfers, Count (i), is estimated considering the node i as the cluster head node: Where E(i) refers to the residual energy of the current node, N denotes the number of adjacent nodes, η is the data fusion coefficient, and Ris the estimated communication distance when data is sent. Step 4: The “fitness factor”, fitness(i), is calculated considering the node i as the cluster head: whereα and β are the proportional coefficients of each item, and the system assumes their values to be –1 and1, respectively. Step 5: Among the neighboring nodes, the one with the largest “fitness factor” is voted as the cluster head node by all the nodes in the network. (2) Networking phase Step 6: The selection coefficient, C(j, i), of each cluster head node from the neighboring nodes is calculated using the cluster member node j, and the cluster head node with the larger coefficient is selected for networking. where E0(i) refers to the residual energy of cluster head node i, Etx(j,i) represents the communication loss between the nodes, and h(j)min/h(i)min guarantees that the cluster head, which is closer to the central node is selected for networking. φ and γ are the proportional coefficients of each item, and the system assumes their values to be –1 and1, respectively. Step 7: The aggregate information of the adjacent cluster heads obtained from the cluster head node (with m as an example) is used for calculating the selection coefficient CH(m, n) of the candidate forwarding node, and the cluster head node with maximum coefficient is chosen as the data forwarding node. where count(n) is the maximum number of data transmissions that can be performed after the node n is selected as the cluster head, h(m)/h(n) guarantees that the cluster head, which is closer to the central node, is preferentially selected to forward data, θ is the proportionality factor for this item, and its value has been taken to be 1 in this study. Step 1: From the abovementioned procedure, the life cycle values of each network are calculated under the condition of m different particles, along with a record of the individual optimal value, W, and the group optimal value, W, in the objective function as well as their corresponding “power adjustment factor” matrices, with the number of cycles increased by one. Step 2: The particle position is updated using the particle update formula, and the steps in Phase II are repeated until the program termination conditions are met.

Phase III: Particle update phase

In the clustered network structure, the functional characteristics of the cluster head node determine whether its energy consumption is much larger than the cluster member nodes. In this study, the cluster head rotation mechanism has also been used for optimizing the network. The rotation threshold of the cluster head is set to q. When 1/5 of the network life cycle is greater than q, either 1/5 of the life cycle of the nodes in the network or q instances of data transmissions are used as the trigger conditions for the network reconstruction.

Performance analysis of the algorithm

Description of the experimental parameters

For verification analysis of our algorithm, the MATLAB is used for conducting the simulations in the range of 500 × 500 m2 and a comparison has been made between the results obtained using the proposed algorithm and the typical LEACH [35], BCDCP [36], and ECHERP [37] algorithms. During the simulations, each wireless sensor node in the WSN is identical regarding initial energy and performance of wireless communication and is randomly fixed in the entire area. The physical structure of the precision agriculture system network is as follows: The central node is located outside the measurement area equipped with an independent power supply. The physical position of wireless sensor equipment in the WSN cannot be changed, but it has the function of transmitting power adjustment using multi-hop as the method of transmission. The parameter settings in the experiment are shown in Table 1.
Table 1

Parameter settings.

ParameterValue
No. of nodes50–300
C1, C20.5
m 10
p 30
E elec 50 nJ/bit
λ 0.5
As shown in Table 1, 50 to 300 wireless sensors are deployed in the network at random; and their physical position cannot be changed, but such equipment has the function of transmitting power adjustment; and power consumption of the node for transmitting/receiving 1bit data Eelec = 50nJ/bit. In the proposed algorithm, Particle Swarm Optimization (PSO) is taken as the basis of algorithm optimization; and in the particle update rules, the number of particles m = 10, selected self-learning factors C1 and C2 = 0.5, inertial coefficient λ = 0.5, and the maximum number of optimizations of the particle p = 30.

Performance analysis of the networking algorithm

Effect of the network size

The number of network nodes, which is an important factor that affects the communication distance between the network nodes, has an important effect on the network life cycle. For networking performance verification of our algorithm over varying sizes of the network, the LEACH, BCDCP, and ECHERP algorithms were used in the simulations with 50–300 nodes selected in an area of 500 × 500 m2. The simulation results are shown in Fig 4. With the gradual addition of the network nodes, the network life cycle for different algorithms is observed to increase gradually, and the proposed algorithm exhibits superior performance in improving the network life cycle. The algorithm comparison on network life cycle under different numbers of nodes is shown in Table 2.
Fig 4

Network life cycle as a function of the number of nodes deployed.

Table 2

Algorithm comparison on network lifecycle under different numbers of nodes.

Number of nodesAlgorithm50100150200250300
LEACH4718301380189021982258
BCDCP4889601603226024602680
ECHERP49311801900230125972732
Proposed algorithm50213022301258727302886

Effect of the position of the central node

In the precision agriculture system, the location of the host determines whether the central node should be deployed outside the monitoring range of the wireless sensor network or not. As an important part of the sensor network, the distance of communication between the central node and the sensing equipment greatly affects the size of the network life cycle. To confirm how this factor affects the network life cycle, a distance of 50 ~ 170 m from the edge of the monitoring range was selected for the simulation and the results thus obtained are shown in Fig 5. As the distance between the device and the position of the central node increases gradually, the network life cycle is observed to decrease gradually. When the position is greater than 130 m, the downward trend becomes steeper. The role of the distance between sensing and central nodes appears crucial. The algorithm proposed in this work exhibits a lower decline rate and downward trend than the compared algorithms. The algorithm comparison on network life cycle under different distances of the central node is shown in Table 3.
Fig 5

Network life cycle as a function of the distance of the device from the position of the central node.

Table 3

Algorithm comparison on network life cycle under different distances of central node.

Distance of central nodeAlgorithm5075100125150175
LEACH18961800168714091001690
BCDCP23862205209816891201834
ECHERP25002402238520451480960
Proposed algorithm259825232443230315131000

Analysis of transmission frequency of the network node

The analysis of the data transmission process of all nodes in the network requires verification of the balance in the energy consumption during data transmission by the network. In this study, 200 sensor nodes were selected in the range of 500 × 500 m2 in the simulation, and the central node was assumed to be 100 m away from the edge of the monitoring range. The failure time of the network node for different algorithms is shown in Fig 6. Compared to the other algorithms, the large-scale failure of the network nodes in the algorithm proposed in this study and the breakdown in its communication occur later with the increasing frequency of data collection. The algorithm comparison on network life cycle under different distances between the central node and other nodes is shown in Table 4.
Fig 6

Failure time of all nodes in the network corresponding to different algorithms.

Table 4

Algorithm comparison on network life cycle under difference distances between central node and other nodes.

Number of other nodesAlgorithm150100500
LEACH4102554070548579
BCDCP6940760285699509
ECHERP73408402923010301
Proposed algorithm75808501947410503

Conclusion

In order to enable limited energy usage by the wireless sensors in the precision agriculture system, an algorithm for high-efficiency wireless sensor networks for a ginseng field has been proposed in this study. Exploiting sensor nodes’ adjustable transmission power, the particle swarm optimization technique has been used to achieve the communication range adjustment of sensing devices in the network. From the aspect of reliable communication, the optimal deployment of the limited energy of the devices in the network has been achieved. The simulation results show that compared to the existing LEACH, BCDCP, and ECHERP algorithms, the proposed algorithm has a higher life cycle for different network sizes with good performance at different deployment locations of the central node. Therefore, the environmental growth variables and indices of farm fields in Jilin, Liaoning, and Heilongjiang in China, and the remote mountainous areas in South Korea and North Korea, have been monitored in real-time using the proposed algorithm. The optimal growth environment of the ginseng fields has been improved using the sensing nodes of the wireless sensor network, which has a significant application value in increasing the yield of ginseng and reducing labor costs. 26 Nov 2021
PONE-D-21-28152
Algorithm for Wireless Sensor Networks in Ginseng Field in Precision Agriculture
PLOS ONE Dear Dr. Chen, Thank you for submitting your manuscript to PLOS ONE. After careful consideration, we feel that it has merit but does not fully meet PLOS ONE’s publication criteria as it currently stands. Therefore, we invite you to submit a revised version of the manuscript that addresses the points raised during the review process. Please submit your revised manuscript by Jan 10 2022 11:59PM. If you will need more time than this to complete your revisions, please reply to this message or contact the journal office at plosone@plos.org. When you're ready to submit your revision, log on to https://www.editorialmanager.com/pone/ and select the 'Submissions Needing Revision' folder to locate your manuscript file. Please include the following items when submitting your revised manuscript:
If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter. Guidelines for resubmitting your figure files are available below the reviewer comments at the end of this letter. A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). You should upload this letter as a separate file labeled 'Response to Reviewers'. A marked-up copy of your manuscript that highlights changes made to the original version. You should upload this as a separate file labeled 'Revised Manuscript with Track Changes'. An unmarked version of your revised paper without tracked changes. You should upload this as a separate file labeled 'Manuscript'. If applicable, we recommend that you deposit your laboratory protocols in protocols.io to enhance the reproducibility of your results. Protocols.io assigns your protocol its own identifier (DOI) so that it can be cited independently in the future. For instructions see: https://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols. Additionally, PLOS ONE offers an option for publishing peer-reviewed Lab Protocol articles, which describe protocols hosted on protocols.io. Read more information on sharing protocols at https://plos.org/protocols?utm_medium=editorial-email&utm_source=authorletters&utm_campaign=protocols. We look forward to receiving your revised manuscript. Kind regards, Chi-Tsun Cheng, Ph.D., M.Sc., B.Eng. Academic Editor PLOS ONE Journal Requirements: When submitting your revision, we need you to address these additional requirements. 1. Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at https://journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and https://journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf 2. Thank you for stating the following in the Acknowledgments/ Financial Disclosure Section of your manuscript: This work was supported by the National Natural Science Foundation of China (Project No. 60972127, 61072111 and 60672156) and Key Project of Jilin Provincial Science and Technology Department (Project No. 20100503) and the Project for Science and Tech-nology Center and Science and Technology Service Platform (Project No. 20180623004TC). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We note that you have provided additional information within the Acknowledgements Section that is not currently declared in your Funding Statement. Please note that funding information should not appear in the Acknowledgments section or other areas of your manuscript. We will only publish funding information present in the Funding Statement section of the online submission form. Please remove any funding-related text from the manuscript and let us know how you would like to update your Funding Statement. Currently, your Funding Statement reads as follows: This work was supported by the National Natural Science Foundation of China (Project No. 60972127, 61072111 and 60672156) and Key Project of Jilin Provincial Science and Technology Department (Project No. 20100503) and the Project for Science and Tech-nology Center and Science and Technology Service Platform (Project No. 20180623004TC). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Please include your amended statements within your cover letter; we will change the online submission form on your behalf. 3. PLOS requires an ORCID iD for the corresponding author in Editorial Manager on papers submitted after December 6th, 2016. Please ensure that you have an ORCID iD and that it is validated in Editorial Manager. To do this, go to ‘Update my Information’ (in the upper left-hand corner of the main menu), and click on the Fetch/Validate link next to the ORCID field. This will take you to the ORCID site and allow you to create a new iD or authenticate a pre-existing iD in Editorial Manager. Please see the following video for instructions on linking an ORCID iD to your Editorial Manager account: https://www.youtube.com/watch?v=_xcclfuvtxQ 4. We note you have included a table to which you do not refer in the text of your manuscript. Please ensure that you refer to Table 1 in your text; if accepted, production will need this reference to link the reader to the Table. 5. Please review your reference list to ensure that it is complete and correct. If you have cited papers that have been retracted, please include the rationale for doing so in the manuscript text, or remove these references and replace them with relevant current references. Any changes to the reference list should be mentioned in the rebuttal letter that accompanies your revised manuscript. If you need to cite a retracted article, indicate the article’s retracted status in the References list and also include a citation and full reference for the retraction notice. Additional Editor Comments: Further elaboration on the experiment setup is required. The rationale of the selected evaluation criteria has not be explained clearly. The connection between the results (i.e. figures) and the intrepretations in the analysis section should be highlighted and extended further. [Note: HTML markup is below. Please do not edit.] Reviewers' comments: Reviewer's Responses to Questions Comments to the Author 1. Is the manuscript technically sound, and do the data support the conclusions? The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented. Reviewer #1: Yes Reviewer #2: Yes ********** 2. Has the statistical analysis been performed appropriately and rigorously? Reviewer #1: Yes Reviewer #2: N/A ********** 3. Have the authors made all data underlying the findings in their manuscript fully available? The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified. Reviewer #1: Yes Reviewer #2: Yes ********** 4. Is the manuscript presented in an intelligible fashion and written in standard English? PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here. Reviewer #1: Yes Reviewer #2: Yes ********** 5. Review Comments to the Author Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters) Reviewer #1: The paper proposed an algorithm for high-efficiency wireless sensor networks for a ginseng field explained well. But parameters settings and description of the experimental parameters need more explanation. May be accepted for the publication. Reviewer #2: The paper presents Algorithm for Wireless Sensor Networks in Ginseng Field in Precision Agriculture. The overall architecture and main component are describes. Since statistical analysis should be performed rigorously, you can add a table contains many aspects comparing existing algorithms to the proposed algorithm and specify the method of testing that had been used in the study as well. ********** 6. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files. If you choose “no”, your identity will remain anonymous but your review may still be made public. Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy. Reviewer #1: No Reviewer #2: No [NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files.] While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com/. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email PLOS at figures@plos.org. Please note that Supporting Information files do not need this step.
16 Dec 2021 List of Responses Dear Editors and Reviewers: Thank you for your letter and for the reviewers’ comments concerning our manuscript entitled “Algorithm for Wireless Sensor Networks in Ginseng Field in Precision Agriculture”, (Manuscript Number: PONE-D-21-28152). Those comments are all valuable and very helpful for revising and improving our paper, as well as the important guiding significance to our researches. We have studied comments carefully and have made corrections which we hope meet with approval. Revised portions are marked in red in the paper. The main corrections in the paper and responses to the reviewer’s comments are as follows: Part one: Journal additional requirements revision 1. Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at https://journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and https://journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf Response: As per your advice and PLOS ONE's style requirements, the following changes are made: the first word of both the title and the subtitle, as well as all the proper nouns and generic names are modified. The postal addresses of the corresponding authors are deleted, and their initials are placed in brackets after their e-mail addresses. The fonts of the level-one headings, level-two headings and level-three headings are changed into 18pt bold, 16pt bold and 14pt bold respectively. The formats of level-four headings 3.3.2.1 Cluster head selection phase and 3.3.2.2 Networking phase are changed into (1) Cluster head selection phase and (2) Networking phase. The word "Figure" used to refer to the figures in the paper is abbreviated as "Fig". The revisions are marked in red. 2. Thank you for stating the following in the Acknowledgments/ Financial Disclosure Section of your manuscript: This work was supported by the National Natural Science Foundation of China (Project No. 60972127, 61072111 and 60672156) and Key Project of Jilin Provincial Science and Technology Department (Project No. 20100503) and the Project for Science and Tech-nology Center and Science and Technology Service Platform (Project No. 20180623004TC). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We note that you have provided additional information within the Acknowledgements Section that is not currently declared in your Funding Statement. Please note that funding information should not appear in the Acknowledgments section or other areas of your manuscript. We will only publish funding information present in the Funding Statement section of the online submission form. Please remove any funding-related text from the manuscript and let us know how you would like to update your Funding Statement. Currently, your Funding Statement reads as follows: This work was supported by the National Natural Science Foundation of China (Project No. 60972127, 61072111 and 60672156) and Key Project of Jilin Provincial Science and Technology Department (Project No. 20100503) and the Project for Science and Tech-nology Center and Science and Technology Service Platform (Project No. 20180623004TC). The funders had no role instudy design, data collection and analysis, decision to publish, or preparation of the manuscript. Please include your amended statements within your cover letter; we will change the online submission form on your behalf. Response: As per your advice and PLOS ONE's relevant requirements, the funding-related text is removed from our manuscript. The funding information declared in the Funding Statement section of our online submission form is correct, and "None" is presented in the Acknowledgements section. 3. PLOS requires an ORCID iD for the corresponding author in Editorial Manager on papers submitted after December 6th, 2016. Please ensure that you have an ORCID iD and that it is validated in Editorial Manager. To do this, go to ‘Update my Information’ (in the upper left-hand corner of the main menu), and click on the Fetch/Validate link next to the ORCID field. This will take you to the ORCID site and allow you to create a new iD or authenticate a pre-existing iD in Editorial Manager. Please see the following video for instructions on linking an ORCID iD to your Editorial Manager account: https://www.youtube.com/watch?v=_xcclfuvtxQ Response: As per your advice, an ORCID iD has been obtained, i.e. 0000-0003-2081-9999, and it has been validated in Editorial Manager. 4. We note you have included a table to which you do not refer in the text of your manuscript. Please ensure that you refer to Table 1 in your text; if accepted, production will need this reference to link the reader to the Table. Response: As per your advice, the table number "Table I" is changed to "Table 1", and a description is added to refer to Table 1 in the text of our manuscript, and the sentence "The parameter settings in the experiment are shown in Table 1." is added. The corrections are marked in red. 5. Please review your reference list to ensure that it is complete and correct. If you have cited papers that have been retracted, please include the rationale for doing so in the manuscript text, or remove these references and replace them with relevant current references. Any changes to the reference list should be mentioned in the rebuttal letter that accompanies your revised manuscript. If you need to cite a retracted article, indicate the article’s retracted status in the References list and also include a citation and full reference for the retraction notice. As per your advice, the whole reference list has been checked, and reference No.13 is changed from "Chinese Pharmacopoeia. Part 12015; 8" to "State Pharmacopoeia Commission. Chinese Pharmacopoeia. Beijing: China Medical Science and Technology Press; 2015. PP. 260-261". Part two, Responses to the reviewer's comments: 1. Additional Editor Comments: Further elaboration on the experiment setup is required. The rationale of the selected evaluation criteria has not be explained clearly. The connection between the results (i.e. figures) and the intrepretations in the analysis section should be highlighted and extended further. Response: As per your advice, further elaboration on the experiment setup is provided, Table 2, Table 3 and Table 4 are added to explain the rationale of the evaluation criteria, and algorithm comparison is highlighted and extended further with respect to the experiment results in different situations. 2. Reviewers #1: The paper proposed an algorithm for high-efficiency wireless sensor networks for a ginseng field explained well. But parameters settings and description of the experimental parameters need more explanation. May be accepted for the publication. Response: As per the reviewer's advice, a detailed description of the experimental parameters is added, which is placed after Table 1 and marked in red, as follows: As shown in Table 1, 50 to 300 wireless sensors are deployed in the network at random; and their physical position cannot be changed, but such equipment has the function of transmitting power adjustment; and power consumption of the node for transmitting/receiving 1bit data Eelec=50nJ/bit. In the proposed algorithm, Particle Swarm Optimization (PSO) is taken as the basis of algorithm optimization; and in the particle update rules, number of particles m=10, selected self-learning factors C1 and C2=0.5, inertial coefficient λ=0.5, and maximum number of optimizations of the particle p=30. 3. Reviewer #2: The paper presents Algorithm for Wireless Sensor Networks in Ginseng Field in Precision Agriculture. The overall architecture and main component are describes. Since statistical analysis should be performed rigorously, you can add a table contains many aspects comparing existing algorithms to the proposed algorithm and specify the method of testing that had been used in the study as well. Response: As per the reviewer's advice, Table 2, Table 3 and Table 4 are added, in which algorithm comparison is provided on network life cycle under different numbers of nodes, on network life cycle under difference distances of central node, and on network life cycle under different distances between central node and other nodes, as follows. In the algorithm comparison, the existing algorithms are contrasted with the proposed algorithm, and the test method is indicated. The revisions are marked in red. Table 2 Algorithm comparison on network lifecycle under different numbers of nodes Number of nodes Algorithm 50 100 150 200 250 300 LEACH 471 830 1380 1890 2198 2258 BCDCP 488 960 1603 2260 2460 2680 ECHERP 493 1180 1900 2301 2597 2732 Proposed algorithm 502 1302 2301 2587 2730 2886 Table 3 Algorithm comparison on network life cycle under different distances of central node Distance of central node Algorithm 50 75 100 125 150 175 LEACH 1896 1800 1687 1409 1001 690 BCDCP 2386 2205 2098 1689 1201 834 ECHERP 2500 2402 2385 2045 1480 960 Proposed algorithm 2598 2523 2443 2303 1513 1000 Table 4 Algorithm comparison on network life cycle under difference distances between central node and other nodes Number of other nodes Algorithm 150 100 50 0 LEACH 4102 5540 7054 8579 BCDCP 6940 7602 8569 9509 ECHERP 7340 8402 9230 10301 Proposed algorithm 7580 8501 9474 10503 We tried our best to improve the manuscript and made some changes to the manuscript. These changes will not influence the content and framework of the paper. And here we did not list the changes but marked in red in the revised paper. We appreciate for Editors/Reviewers/ warm work earnestly, and hope that the correction will meet with approval. Once again, thank you very much for your comments and suggestions. Submitted filename: renamed_766a8.docx Click here for additional data file. 19 Jan 2022 Algorithm for wireless sensor networks in ginseng field in precision agriculture PONE-D-21-28152R1 Dear Dr. Chen, We’re pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it meets all outstanding technical requirements. Within one week, you’ll receive an e-mail detailing the required amendments. When these have been addressed, you’ll receive a formal acceptance letter and your manuscript will be scheduled for publication. An invoice for payment will follow shortly after the formal acceptance. To ensure an efficient process, please log into Editorial Manager at http://www.editorialmanager.com/pone/, click the 'Update My Information' link at the top of the page, and double check that your user information is up-to-date. If you have any billing related questions, please contact our Author Billing department directly at authorbilling@plos.org. If your institution or institutions have a press office, please notify them about your upcoming paper to help maximize its impact. If they’ll be preparing press materials, please inform our press team as soon as possible -- no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org. Kind regards, Chi-Tsun Cheng, Ph.D., M.Sc., B.Eng. Academic Editor PLOS ONE Additional Editor Comments (optional): All issues have been resolved in the current version. The paper is recommended to be accepted. Reviewers' comments: Reviewer's Responses to Questions Comments to the Author 1. If the authors have adequately addressed your comments raised in a previous round of review and you feel that this manuscript is now acceptable for publication, you may indicate that here to bypass the “Comments to the Author” section, enter your conflict of interest statement in the “Confidential to Editor” section, and submit your "Accept" recommendation. Reviewer #1: All comments have been addressed Reviewer #2: All comments have been addressed ********** 2. Is the manuscript technically sound, and do the data support the conclusions? The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented. Reviewer #1: Yes Reviewer #2: Yes ********** 3. Has the statistical analysis been performed appropriately and rigorously? Reviewer #1: Yes Reviewer #2: Yes ********** 4. Have the authors made all data underlying the findings in their manuscript fully available? The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified. Reviewer #1: No Reviewer #2: Yes ********** 5. Is the manuscript presented in an intelligible fashion and written in standard English? PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here. Reviewer #1: Yes Reviewer #2: Yes ********** 6. Review Comments to the Author Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters) Reviewer #1: (No Response) Reviewer #2: (No Response) ********** 7. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files. If you choose “no”, your identity will remain anonymous but your review may still be made public. Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy. Reviewer #1: No Reviewer #2: No 27 Jan 2022 PONE-D-21-28152R1 Algorithm for wireless sensor networks in ginseng field in precision agriculture Dear Dr. Chen: I'm pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now with our production department. If your institution or institutions have a press office, please let them know about your upcoming paper now to help maximize its impact. If they'll be preparing press materials, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information please contact onepress@plos.org. If we can help with anything else, please email us at plosone@plos.org. Thank you for submitting your work to PLOS ONE and supporting open access. Kind regards, PLOS ONE Editorial Office Staff on behalf of Dr. Chi-Tsun Cheng Academic Editor PLOS ONE
  5 in total

Review 1.  Perspectives for Remote Sensing with Unmanned Aerial Vehicles in Precision Agriculture.

Authors:  Wouter H Maes; Kathy Steppe
Journal:  Trends Plant Sci       Date:  2018-12-15       Impact factor: 18.313

2.  An ion-selective crown ether covalently grafted onto chemically exfoliated MoS2 as a biological fluid sensor.

Authors:  Anastasios Stergiou; Christina Stangel; Ruben Canton-Vitoria; Ryo Kitaura; Nikos Tagmatarchis
Journal:  Nanoscale       Date:  2021-05-20       Impact factor: 7.790

Review 3.  Insights into the antitumor mechanism of ginsenosides Rg3.

Authors:  Zongyu Liu; Tongjun Liu; Wei Li; Jiannan Li; Cuizhu Wang; Kai Zhang
Journal:  Mol Biol Rep       Date:  2021-03-04       Impact factor: 2.316

Review 4.  Energy-Efficient Wireless Sensor Networks for Precision Agriculture: A Review.

Authors:  Haider Mahmood Jawad; Rosdiadee Nordin; Sadik Kamel Gharghan; Aqeel Mahmood Jawad; Mahamod Ismail
Journal:  Sensors (Basel)       Date:  2017-08-03       Impact factor: 3.576

5.  A Scalable Context-Aware Objective Function (SCAOF) of Routing Protocol for Agricultural Low-Power and Lossy Networks (RPAL).

Authors:  Yibo Chen; Jean-Pierre Chanet; Kun-Mean Hou; Hongling Shi; Gil de Sousa
Journal:  Sensors (Basel)       Date:  2015-08-10       Impact factor: 3.576

  5 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.