- School of Marine Science and Technology, Northwestern Polytechnical University, Xi'an, China
This paper presents a multi-objective evolutionary algorithm of bio-inspired geomagnetic navigation for Autonomous Underwater Vehicle (AUV). Inspired by the biological navigation behavior, the solution was proposed without using a priori information, simply by magnetotaxis searching. However, the existence of the geomagnetic anomalies has significant influence on the geomagnetic navigation system, which often disrupts the distribution of the geomagnetic field. An extreme value region may easily appear in abnormal regions, which makes AUV lost in the navigation phase. This paper proposes an improved bio-inspired algorithm with behavior constraints, for sake of making AUV escape from the abnormal region. First, the navigation problem is considered as the optimization problem. Second, the environmental monitoring operator is introduced, to determine whether the algorithm falls into the geomagnetic anomaly region. Then, the behavior constraint operator is employed to get out of the abnormal region. Finally, the termination condition is triggered. Compared to the state-of- the-art, the proposed approach effectively overcomes the disturbance of the geomagnetic abnormal. The simulation result demonstrates the reliability and feasibility of the proposed approach in complex environments.
Introduction
Autonomous Underwater Vehicle (AUV) has been widely used for both civilian and military applications, such as laying pipelines, ocean data collection, underwater equipment maintenance, and laying mines (Wadhams, 2012; Wynn et al., 2014; Shi et al., 2017). Most navigation systems rely on dead-reckoning that is by using inertial navigation system and velocity information from the Doppler sonar. However, the cumulative error becomes difficult to handle during the navigation, in contrast to the GPS and the acoustic transponder networks (Caiti et al., 2014).
Various techniques have been developed to eliminate the dead-reckoning errors (Hao et al., 2008; Yi et al., 2008; Shen et al., 2017). The geomagnetic navigation, which is originated from animal behaviors due to the earth's magnetic fields, has been widely used (Gould, 1984; Lohmann, 2010). It overcomes the following drawbacks: the error accumulation from the inertial navigation system, and the rapid attenuation signal from the satellite navigation system (Goldenberg, 2006). Thus, the geomagnetic field plays a key role in navigation, as it provides the position information in large scale environments (Teixeira and Pascoal, 2008). Typically, each point on the near-earth space has an unique magnetic field vector with respect to the corresponding coordinate, while the geomagnetic navigation is employed to provide a reliable navigation reference (Fu-qing, 2006; Zhou et al., 2008).
Based on the characteristics of the geomagnetism, numbers of methods have been developed (Caifa et al., 2011). The principle idea is employed based on the matching algorithms, such as MSD (Mean Square Difference), MAD (Mean Absolute Difference), ICCP (Iterative Closest Contour Point Algorithm) (Jia et al., 2012; Xie et al., 2013; Chong et al., 2014). However, traditional matching algorithms strongly depend on a priori geomagnetic map, which is quite challenging to acquire in practice.
To address this issue, the bio-inspired geomagnetic navigation method was proposed from animal behaviors (Liu et al., 2013). The most informative experimental paradigms have verified that the animals (like sea turtles and pigeons) can geomagnetic navigation to reach their goal locations (Paolo et al., 2003). Schulten concluded that animals can navigate by relying on geomagnetic sensitivity (Schulten, 1982). Winklhofer proposed the magneto-reception mechanisms to explain animals' navigation behavior (Winklhofer, 2009). Mole rats and salmons are also sensitive to the earth's magnetic field (Kimchi et al., 2004; Hays, 2013). In conclusion, the bio-inspired geomagnetic navigation provides a natural solution for navigation problem without a priori geomagnetic map.
However, the geomagnetic anomaly often influences the bio-inspired geomagnetic navigation system which is caused by large iron ores. Boström showed that many geomagnetic blind spots were existed with respect to some migrating animals (Boström et al., 2012). Kiliowska pointed out the geomagnetic anomalies have interfered the whales' navigation system and caused them to beach (Klinowska, 1986). Lohmann pointed out the geomagnetic anomalies made the migrating turtles lost into Mexico (Lohmann et al., 2007). Dennis pointed out the pigeons were also easily lost when they were released in the geomagnetic anomaly area (Dennis et al., 2007). With respect to an AUV navigation, the geomagnetic anomalies disrupt the distribution of the geomagnetic field, which makes AUV lost in the navigation phase. In our previous works, the bio-inspired navigation algorithm was easily trapped in a local minimum point, caused by the geomagnetic anomalies (Liu et al., 2014). This happens where the distribution of the multiple geomagnetic was changed to the unsmooth and discontinuous regions (“concave” or “convex”). This caused the multi-objective search to converge to a local minimum easily (see Figure 1).
This paper presents a multi-objective evolutionary algorithm to address the local minimum problem by using the behavior constraint strategy. Main contributions are summarized as follows: First, it determines whether the algorithm is falling into the abnormal area. Second, it provides the strategy to jump out of the geomagnetic anomaly area. Once the algorithm is trapped in the local minimum, the behavior constraints are utilized to jump out the abnormal area.
The rest of this paper is organized as follows: Section Problem Formulation briefly introduces both the multi-objective search and the local minimum problems, respectively. Section Multi-Objective Evolutionary Algorithm with Geomagnetic Anomaly Presents the Multi-Objective Evolutionary Algorithm in Scenarios of Geomagnetic Anomaly. Section Numerical Simulations gives the main results and theoretical analysis. Finally, the conclusion is given in Section Conclusion.
Problem Formulation
In this paper, the motion of AUV on the horizontal plane is considered in the 2D Cartesian coordinate system. This simplification is justified due to the fact that the difference of the geomagnetic in vertical is often negligible. Therefore, the kinematic equations of the motion are introduced as follows:
where ΔT is the sample period, k is the time instant, v(k) is the kinematic velocity, θ(k) is the heading of AUV in a time instant. Here, assuming AUV moves with a constant velocity, formula (1) can be written as follows:
where L is the step size, L = v(k)ΔT.
Mathematical Description of the Multi-Objective Search Problem
The geomagnetic fields include multiple geomagnetic components (Caifa et al., 2011), which can be described as follows:
where B is the set of the geomagnetic components vector, B1, B2, ⋯, Bn are defined as geomagnetic components, such as the north magnetic field Bx, the east magnetic field By, the downward magnetic field Bz, the total intensity BF, the horizontal magnetic field BH, the declination angle BD, and the declination angle BI.
The bio-inspired navigation is thus the convergence process of the geomagnetic components from the current position to the target position (Liu et al., 2014). Therefore, the navigation process can be considered as the multi-objective searching problem as follow:
where is the geomagnetic component of the target position; is the geomagnetic component of the current position; is the difference of the ith geomagnetic component between the target position and the current position; Fk is the objective function.
Considering the different magnitude of the geomagnetic components, the objective function is normalized as follows:
The errors between the current position and the target position can be assumed as the geomagnetic trend, while the searching process is terminated when the error converges to ε. This is expressed as:
where ε is a fixed value.
The bio-inspired navigation is thus considered as a posteriori searching problem, in presence of unknown geomagnetic components. Once the geomagnetic error between the current position and the target position converges to zero, the navigation process is terminated.
Geomagnetic Searching in Anomaly Environment
The geomagnetic anomaly field has significant influences on the geomagnetic navigation system, where the related factors are analyzed according to the distribution characteristics of the geomagnetic field (Talwani, 1965). Usually the geomagnetic anomaly is given by:
where B is the amplitude of the abnormal region.
Considering the geomagnetic anomaly for the bio-inspired algorithm, the difference of the ith geomagnetic component is expressed as:
Combining formulas (7) and (8), the convergence condition is given by:
where ξ is bigger than ε if the geomagnetic anomalies exists. This result shows that the algorithm is easily trapped in the local minimization.
Multi-Objective Evolutionary Algorithm with Geomagnetic Anomaly
This section gives the general idea of the proposed algorithm in geomagnetic anomaly, the experiment will be given in next section.
The Search Principle
A bio-inspired navigation method is investigated, which is based on a simple assumption that a homing animal only senses and compares the variation of the geomagnetic field to reach its home.
Evolutionary algorithm is often used to solve problems in highly complex spaces (Cliff et al., 1993; Droste et al., 2002; Vrugt and Robinson, 2007; Peng et al., 2016; Peng and Wu, 2017). Based on Darwin's evolutionary theory, ethologists have been concerned with the evolution of animals' behavior (Alerstam et al., 2003). Searching behavior can be described as active movement, in which animals attempt to find resources such as food, mates, nesting sites (Bell, 1992). Based on this, McFarland has proposed a number of mathematical analogs by using adaptive evolutionary behaviors (McFarland and Bösser, 1993).
Hence, by utilizing the evolution of animals behaviors, the evolutionary algorithm is employed for searching the geomagnetic space. In evolutionary algorithm, each possible solution is defined as the individual, whereas a set of individuals is defined as the population. The evolution method requires maintaining a population of various individuals, according to the operators such as selection, mutation, and so on. The repeated process of recombination, selection, and mutation leads the individuals to adapt to the environment, whereas the selection operator is implemented by a task-oriented evaluation function: the better the AUV performs its task, the more offspring it has.
However, the local anomalies from the geomagnetic field have significant influences on the proposed geomagnetic navigation system. The next section explains the solution to the anomaly problem.
The Search Algorithm with Geomagnetic Anomaly
This section discusses how to deal with the local minimum problem within the geomagnetic anomaly areas. Our solution has two phases: the first phase is to determine whether the algorithm falls into the geomagnetic anomaly area, and the second is to determine whether the algorithm arrives at the target location.
The first phase refers to the environmental monitoring operator Ψ(k), which can be expressed as:
where T1 is the period of time, E(·) is described as mean values.
The termination condition in second phase is determined as:
where ρ is a small value zero, φ is the scope of the desired target location. The steps for the multi-objective evolutionary algorithm which can address the challenge of the geomagnetic anomaly are given as follows:
Step 1: Population Initialization. Randomly generate N individuals in the population space Q (Q = {θj|j = 1, 2, ⋯, N }), θ is the set of samples and given by:
where θ is the sampling interval.
Step 2: Individual Selection. Randomly select the sample θj, and the probability of each individual is given by:
when the current number of iterations is k, the probability of the selected individual θi can be obtained as
where δ is a symbolic function.
Step 3: Population Updating Rules. The termination is calculated if the selected individual accomplishes a successful search. Thus the individual updating rules can be divided into two parts: for one that moves toward the target direction (), the selected individual heading is reserved in the population. And for the other that moves away from the target direction (), the selected individual heading is replaced.
Step 4: Environmental Monitoring Operator. The monitoring factor Ψ(k) is used to determine whether the algorithm falls into the geomagnetic anomaly area (Ψ(k) ≥ μ, μ is the threshold value).
Step 5: Motion Constraint Operator. Once the search algorithm is trapped into the local minimum area, the maximum probability θi is preserved. This can be expressed as:
When formula (17) is satisfied, the algorithm jumps out of the local minimum area.
Step 6: Termination condition. The navigation is terminated if formula (11) is satisfied. The workflow of the proposed evolutionary algorithm is shown in Figure 2.
Numerical Simulations
To verify the effectiveness of the proposed multi-objective evolutionary algorithm, numerical simulations are performed.
Simulation Setup
The Word Magnetic Model (WMM2010) is used to provide the real time geomagnetic data (Maus et al., 2012). In simulation, a rectangular area is selected from north latitude 30° and east longitude 100° (N.30° and E.100°) to north latitude 45° and east longitude 135° (N.45° and E.135°). Considering the non-relevance of the geomagnetic components, only three geomagnetic components are used, which are the north magnetic field Bx, the east magnetic field By, and the total intensity BF, respectively.
Meanwhile, it is assumed that the geomagnetic navigation would encounter a geomagnetic anomaly region. Here, the multi-mode function is utilized to construct the abnormal geomagnetic environment (see in Figure 3), where the intensity of the highest abnormal geomagnetic field value is −5,000nT.
In Figure 3, the starting position and the target position are depicted by using the three geomagnetic components at B0 = [31, 464nT, −2, 036nT, 52, 508nT] and Bt = [26, 290nT, −4, 229nT, 54, 254nT], respectively. The circle, “◦”, stands for the starting position, and the star, “*,” stands for the target position. The related parameters are shown in Table 1.
Simulation Results
To evaluate the performance of the proposed approach, two cases are evaluated: one occurs without considering geomagnetic anomaly and the other one with geomagnetic anomaly.
Figure 4 illustrates the navigation trajectory without the interference of geomagnetic anomalies, whereas Figure 5 illustrates the moving trajectory with the interference of geomagnetic anomalies. As shown in Figure 5, the previous algorithm is easily trapped into the local minimization in geomagnetic anomaly areas, in which individuals move toward different directions (Liu et al., 2014). It is observed that the navigation performance strongly depends on the geomagnetic anomalies.
It is concluded that the original searching algorithm is not capable for the anomalies of geomagnetic fields, in which AUV moves toward different directions and ultimately fails.
Figure 6 illustrates the proposed approach, where AUV successfully overcomes the influence of the geomagnetic anomalies. During the navigation phase, the individual directions are randomly changed due to the probability selected in the population. Zoomed figure on the left corner shows the two motion constraints, which are performed using the monitoring factor Ψ(k). Two time periods “c” and “d” are within the geomagnetic anomaly area, while “a” and “b” represent the variation anomalies of the geomagnetic components, respectively. By monitoring the changes of the geomagnetic components, the proposed approach effectively detects the abnormal regions. Then, the statistic characteristics of the convergence state is utilized to the behavior constraints, for the purpose of jumping out the abnormal region.
Convergence Performance of the Search Algorithm
The convergence curves of three geomagnetic components are shown in Figure 7. It is observed that the convergence curves present violent shakings in three geomagnetic components.
To better evaluate the performance in presence of geomagnetic anomalies, the convergence properties of three geomagnetic components are provided. As depicted in Figure 8, there is a divergence trend for the three geomagnetic components in periods of “a” and “b.” The convergence properties of the geomagnetic components are improved after using behavior constraints, demonstrating that the normalized error are able to converge to a stable state. It is observed that the improved algorithm could converge the curves of the geomagnetic components to zero in finite steps.
To explore the advantages of the proposed algorithm, the traveling time employed by the AUV in both scenarios is also compared in Figure 9. It shows that the proposed approach could converge to a stable state in scenarios of geomagnetic anomalies, and the navigation process is terminated in 540 steps.
Discussion
In this paper, the multi-objective evolutionary algorithm has been proposed to address the AUV navigation problem. As a result of geomagnetic anomaly, the distribution of the geomagnetic components is often disrupted over the searching space. Therefore, the environmental monitoring and the behavior constraints are used to tackle the local optimal problem. The time complexity of the proposed algorithm is O(k3) in the geomagnetic anomalies.
For the geomagnetic navigation, the optimal path between the starting position and the destination is quite challenging without a priori geomagnetic map. Also, the AUV navigation in oceans strongly influences by the geomagnetic anomaly effects. However, the proposed approach indicates that the searching strategy performs optimally for AUV navigation in the geomagnetic anomaly areas.
Thus, the multi-objective searching algorithm with the environmental monitoring and the behavior constraints is proposed to ensure the success and safety of AUV navigation.
Conclusion
This paper presents a novel strategy for bio-inspired geomagnetic navigation in presence of geomagnetic anomalies. Inspired by the biological navigation, in our previous work we proposed an evolutionary schema which helps the AUV to reach the destination without geomagnetic map. However, it suffers the disturbance from geomagnetic anomalies and often converges to the local minimum point. To tackle the problem, this paper proposes an improved navigation model by introducing constraints strategy to make AUV escape from the abnormal regions. Simulation results show that the proposed model effectively overcomes the disturbance of the geomagnetic abnormal for AUV navigation.
Author Contributions
HL wrote the paper and performed the experiments. ML provided some ideas to improve and perfect the paper. FZ reviewed and edited the manuscript. All authors read and approved the manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (NSFC) under Grants 51379176.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors are also grateful for the constructive suggestions from reviewers that significantly enhance the presentation of this paper.
References
Alerstam, T., Hedenström, A., and Åkesson, S. (2003). Long-distance migration: evolution and determinants. Oikos 103, 247–260. doi: 10.1034/j.1600-0706.2003.12559.x
Bell, W. J. (1992). Searching behaviour: the behaviour ecology of finding resources. J. Anim. Ecol. 61, 370–371.
Boström, J. E., Åkesson, S., and Alerstam, T. (2012). Where on earth can animals use a geomagnetic bicoordinate map for navigation? Ecography 35, 1039–1047. doi: 10.1111/j.1600-0587.2012.07507.x
Caifa, G., Anliang, L., Hong, C., and Huabo, Y. (2011). Algorithm for geomagnetic navigation and its validity evaluation. Int. Conf. Comput. Sci. Autom. Eng. 1, 573–577. doi: 10.1109/CSAE.2011.5953286
Caiti, A., Di Corato, F., Fenucci, D., Ridolfi, A., Pugi, L., Monni, N., et al. (2014). Experimental results with a mixed USBL/LBL system for AUV navigation. Underwater Commun. Netw. 2014, 1–4. doi: 10.1109/UComms.2014.7017129
Chong, K., Zhang, X., and Fan, L. (2014). Geomagnetic-matching aided navigation based on ICCP algorithm. J. Basic Sci. Eng. 22, 598–605. doi: 10.3969/j.issn.1005-0930.2014.03.018
Cliff, D., Husbands, P., and Harvey, I. (1993). Explorations in evolutionary robotics. Adaptive Behav. 2, 73–110. doi: 10.1177/105971239300200104
Dennis, T. E., Rayner, M. J., and Walker, M. M. (2007). Evidence that pigeons orient to geomagnetic intensity during homing. Proc. R. Soc. Lond. B Biol. Sci. 274, 1153–1158. doi: 10.1098/rspb.2007.3768
Droste, S., Jansen, T., and Wegener, I. (2002). On the analysis of the (1+1) evolutionary algorithm. Theor. Comp. Sci. 276, 51–81. doi: 10.1016/S0304-3975(01)00182-7
Fu-qing, P. (2006). Geomagnetic model and geomagnetic navigation. Hydrogr. Surv. Chart. 2, 73–75. doi: 10.3969/j.issn.1671-3044.2006.02.021
Goldenberg, F. (2006). Geomagnetic navigation beyond the magnetic compass. Proc. IEEE/ION Plans 2001, 684–694. doi: 10.1109/PLANS.2006.1650662
Gould, J. L. (1984). Magnetic field sensitivity in animals. Annu. Rev. Physiol. 46, 585–598. doi: 10.1146/annurev.ph.46.030184.003101
Hao, Y. L., Zhao, Y. F., and Jun-Feng, H. U. (2008). Preliminary analysis on the application of geomagnetic field matching in underwater vehicle navigation. Prog. Geophys. 23, 594–598. doi: 10.6038/pg20080601
Hays, G. C. (2013). Animal navigation: salmon track magnetic variation. Curr. Biol. 23, R144–R145. doi: 10.1016/j.cub.2013.01.025
Jia, L., Wang, Y. G., Shan, B., and Yan, T. (2012). Simulation of single element geomagnetic matching navigation based on intensified MAD. Mod. Def. Technol. 1, 90–94. doi: 10.3969/j.issn.1009-086x.2012.01.019
Kimchi, T., Etienne, A. S., and Terkel, J. (2004). A subterranean mammal uses the magnetic compass for path integration. Proc. Natl. Acad. Sci. U.S.A. 101, 1105–1109. doi: 10.1073/pnas.0307560100
Klinowska, M. (1986). Cetacean live stranding dates relate to geomagnetic disturbances. Aquat. Mammals 11, 109–119.
Liu, M., Li, H., and Liu, K. (2014). Geomagnetic navigation of AUV without a priori magnetic map[C]//OCEANS 2014-TAIPEI. IEEE 2014, 1–5. doi: 10.1109/OCEANS-TAIPEI.2014.6964431
Liu, M., Liu, K., Yang, P., and Lei, X. (2013). Bio-inspired navigation based on geomagnetic robotics and biomimetics (ROBIO). IEEE International Conference on IEEE 2013 (Shenzhen), 2339–2344.
Lohmann, K. J. (2010). Animal behaviour: magnetic-field perception. Nature 464, 1140–1142. doi: 10.1038/4641140a
Lohmann, K. J., Lohmann, C. M. F., and Putman, N. F. (2007). Magnetic maps in animals: nature's GPS. J. Exp. Biol. 210, 3697–3705. doi: 10.1242/jeb.001313
Maus, S., Macmillan, S., McLean, S., and Hamilton, B. (2012). The US/UK World Magnetic Model for 2010-2015. Available online at: https://www.ngdc.noaa.gov/geomag/WMM/DoDWMM.shtml
Paolo, L., Hays, G. C., and Floriano, P. (2003). A review of long-distance movements by marine turtles, and the possible role of ocean currents. Oikos 103, 293–302. doi: 10.1034/j.1600-0706.2003.12123.x
Peng, X., and Wu, Y. (2017). Large-scale cooperative co-evolution using nichingbased multi-modal optimization and adaptive fast clustering. Swarm Evol. Comput. doi: 10.1016/j.swevo.2017.03.001. [Epub ahead of print].
Peng, X., Liu, K., and Jin, Y. (2016). A dynamic optimization approach to the design of cooperative co-evolutionary algorithms. Knowledge-Based Syst. 109, 174–186. doi: 10.1016/j.knosys.2016.07.001
Schulten, K. (1982). “Magnetic field effects in chemistry and biology,” in Festkörperprobleme 22, eds K. Schulten, H. Staerk, and A. Weller (Berlin Heidelberg: Springer), 61–83.
Shen, C., Shi, Y., and Buckham, B. (2017). Integrated path planning and tracking control of an AUV: a unified receding horizon optimization approach. IEEE/ASME Trans. Mech. 22, 1163–1173. doi: 10.1109/TMECH.2016.2612689
Shi, Y., Shen, C., Fang, H., and Li, H. (2017). Advanced control in marine mechatronic systems: a survey. IEEE/ASME Trans. Mech. 22, 1121–1131. doi: 10.1109/tmech.2017.2660528
Talwani, M. (1965). Computation with the help of a digital computer of magnetic anomalies caused by bodies of arbitrary shape. Geophysics 30, 797–817. doi: 10.1190/1.1439654
Teixeira, F., and Pascoal, A. (2008). “Geophysical navigation of autonomous underwater vehicles using geomagnetic information,” in Proceedings of NGCUV'08-IFAC Workshop on Navigation, Guidance and Control of Underwater Vehicles (Ireland).
Vrugt, J. A., and Robinson, B. A. (2007). Improved evolutionary optimization from genetically adaptive multimethod search. Proc. Natl. Acad. Sci. U.S.A. 104, 708–711. doi: 10.1073/pnas.0610471104
Wadhams, P. (2012). The use of autonomous underwater vehicles to map the variability of under-ice topography. Ocean Dyn. 62, 439–447. doi: 10.1007/s10236-011-0509-1
Winklhofer, M. (2009). The physics of geomagnetic-field transduction in animals. IEEE Trans. Magn. 45, 5259–5265. doi: 10.1109/TMAG.2009.2017940
Wynn, R. B., Huvenne, V. A. I., Le Bas, T. P., Murton, B. J., Connelly, D. P., Bett, B. J., et al. (2014). Autonomous underwater vehicles (AUVs): their past, present and future contributions to the advancement of marine geosciences. Mar. Geol. 352, 451–468. doi: 10.1016/j.margeo.2014.03.012
Xie, W., Qu, Z., and Li, Q. (2013). A fast algorithm of the geomagnetic correlation matching based on MSD. Adv. Intell. Syst. Res. 45, 59–62. doi: 10.2991/case-13.2013.15
Yi, L., Lei, Y., and Tong, Q. (2008). “Underwater geomagnetic navigation based on ICP algorithm,” in IEEE International Conference on Robotics and Biomimetics (Sanya), 2115–2120.
Keywords: geomagnetic navigation, bio-inspired navigation, geomagnetic anomaly, local optimal, evolutionary algorithm
Citation: Li H, Liu M and Zhang F (2017) Geomagnetic Navigation of Autonomous Underwater Vehicle Based on Multi-objective Evolutionary Algorithm. Front. Neurorobot. 11:34. doi: 10.3389/fnbot.2017.00034
Received: 12 October 2016; Accepted: 23 June 2017;
Published: 12 July 2017.
Edited by:
Qiang Huang, Beijing Institute of Technology, ChinaReviewed by:
Yang Shi, University of Victoria, CanadaTara Julia Hamilton, Western Sydney University, Australia
Copyright © 2017 Li, Liu and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Feihu Zhang, ZmVpaHUuemhhbmdAbndwdS5lZHUuY24=