Stochastic Systems: Modeling, Optimization, and ApplicationsView this Special Issue
Research Article | Open Access
Localization Accuracy of Range-Only Sensors with Additive and Multiplicative Noise
The localization accuracy of range-only sensors with both additive and multiplicative noises is investigated. The main contributions consist of three parts. First, the Cramer-Rao Lower Bound (CRLB) of range-only localization system is derived, based on which the area of the uncertain ellipse about estimation is proposed as the metric to evaluate the localization accuracy. In addition, the analytical relationships between target-sensor distance and localization accuracy and between noise and localization accuracy are derived. Second, the metric is utilized to evaluate the localization accuracy for three important regular patterns of sensor deployment, that is, triangle, square, and hexagon. Two aspects of localization accuracy have been examined and discussed, respectively, for three patterns by numerical methods, including (a) the geometric distribution of localization accuracy and (b) the average localization accuracy. Both theoretical and numerical results show that the multiplicative noise will influence significantly the localization accuracy. Third, on a lab-scale ultrasonic range-only sensor system the existence of additive and multiplicative noises is verified. It is also shown that the experimental localization accuracy is close to the previously analyzed theoretical accuracy.
It is well known that target localization plays an important role in military and security , which is also the foundation of other advanced tasks, for example, target tracking. With the great development of electronics in recent years, wireless sensor networks (WSN) based target localization systems have attracted increasing research interests [1–9]. Several kinds of sensors are widely used to collect the necessary information, for example, range-only sensors [2, 10], angle of arrival (AOA) sensors [10–12], time of arrival (TOA) sensors , and time difference of arrival (TDOA) sensors . In a WSN based localization system, the target localization problem is essentially an estimation problem. That is, a number of sensors are deployed to collect the position-related data of the target; then the data are sent to the leader node to estimate the position of the target.
In the context of statistical signal processing, the Cramer-Rao Lower Bound (CRLB) is regarded as the most important benchmark to assess the estimation accuracy of an unbiased estimator [13–16]. Therefore, CRLB is widely used in the localization system to evaluate the localization accuracy [1, 17–21]. It is also of great significance in sensor selection and deployment, motion coordination, and dynamic sensor collaboration [7, 11, 22]. With regard to different measurement models, CRLB can be derived in various forms. For example, CRLB of the widely used range-only localization system is a matrix of target-sensor azimuth angle in previous literature .
However, it is notable that this result is derived under the assumption that the measurements consist of only additive noise . In other words, this assumption ignores the multiplicative noise that indeed reflects the effect of target-sensor distance on measurement error. However,  argues that in a lot of occasions multiplicative noise cannot be neglected, and the measurement error will increase significantly with the increase of target-sensor distance. In fact, the experimental results in  even show that multiplicative noise is 3 to 4 times larger than additive noise. Unfortunately, very few studies take multiplicative noise into consideration when deriving CRLB and assessing localization accuracy for the range-only localization system.
In this paper, the localization accuracy problem of range-only sensors with both additive and multiplicative noises is investigated. The main contributions can be summarised as below.(1)Theoretically, CRLB of range-only localization system is deduced based on a more general measurement model with both additive and multiplicative noise. Then, a metric is proposed based on CRLB to evaluate the localization accuracy, and three important properties of this metric are derived. The properties indicate, surprisingly, that shorter target-sensor distance does not always lead to better localization accuracy, and proper distance will yield the best accuracy under a certain situation.(2)Numerically, the localization accuracy for three regular patterns of sensor deployment, that is, triangle, square, and hexagon, has been evaluated by the proposed metric. For each pattern, two aspects of localization accuracy have been examined and discussed, including (a) geometric distribution of localization accuracy and (b) average localization accuracy. To the best of our knowledge, there is no such extensive analysis about the localization accuracy for regular deployment patterns. Numerical results show that the multiplicative noise will influence significantly the localization accuracy and also provide some important guidelines for optimal sensor deployment.(3)Experimentally, the existence of multiplicative noise on a lab-scale ultrasonic range-only sensor system is verified and the corresponding parameters are computed. Then, the experimental localization accuracy is calculated based on real data, which is very close to the theoretical accuracy.
The rest of this paper is organized as follows. In Section 2, the measurement model of range-only localization system with both additive and multiplicative noise is introduced. Then, CRLB, the metric, and its properties are derived in Section 3. Section 4 presents the numerical analysis and discussion about the localization accuracy in terms of the proposed metric in regular deployment patterns. The experiment is presented in Section 5 and the conclusion is drawn in Section 6.
2. Measurement Model
Without loss of generality, we will consider the localization problem for a single stationary target in two-dimensional field in this paper. The presented results can be readily extended to three-dimensional cases and are also the foundation of target tracking tasks.
Let the localization system consist of homogeneous range-only sensors, whose positions are known a priori. Let the coordinate of the th sensor be , and let the coordinate of the target be . Then, the true distance between the target and sensor can be expressed as
Being different from most previous studies , which adopt the sensor measurement with only additive noise, a more general model with both additive and multiplicative noise  is considered here: where denotes the measurement of the th sensor and and are, respectively, multiplicative noise and additive noise in the th sensor.
Assume that and are independent white Gaussian noises, that is, , , which are also independent of those in the other sensors. From (2), the total measurement noise of sensor is . Clearly, it is still white Gaussian noise with the following mean and variance:
Therefore, the measurement also follows normal distribution: Then, consider the entire localization system that consists of sensors. Similar to (5), the measurement vector of sensors follows an -dimensional normal distribution: where , , , and .
Remark 1. From (3) and (4), it can be found that both the mean and the variance associated with sensor are closely related to target-sensor distance . In particular, the variance will grow in a rapid speed with the increase of .
For a target at , the probability density function (PDF) of measurement of sensor can be readily obtained according to (5): Since the measurements of different sensors are independent, the joint PDF of the entire localization system with sensors can be written as
3. CRLB and Metric
In this section, CRLB of range-only localization system with both additive and multiplicative noise is derived and a metric is proposed to evaluate the localization accuracy. Then, three important analytical properties about the relationships between target-sensor distance and localization accuracy and between noise and localization accuracy are derived.
3.1. Derivation of CRLB
It is well known that CRLB is the lower bound for the variance of any unbiased estimate. Thus, it is widely utilized in practice as an important benchmark to evaluate the estimation accuracy of an estimator [13–16, 23, 24]. The inverse of CRLB is the so-called Fisher Information Matrix (FIM), which stands for the amount of information contained in data.
By referring the definition of FIM in , an entry of FIM denoted as is where , , and . Therefore, the corresponding FIM can be expressed as
Alternatively, since the measurement given in (6) is Gaussian with nonzero mean, the entry of FIM can also be computed as below : Because the measurements of different sensors are independent, the complete FIM of the entire localization system is the sum of the individual FIM associated with each sensor . Then, taking (1), (3), and (4) into (11), we have
Let denote the azimuth angle of sensor relative to the target. Obviously, we have , . As a result, the complete FIM of sensors can be rewritten into the following compact form: where is a function of noise and target-sensor distance of sensor At last, the corresponding CRLB can be expressed as
Remark 2. Similar to previous studies, for example, [1, 2, 12], in this paper we will assume that all estimates of the target position are unbiased. Therefore, the derived CRLB corresponds to the lower bound that the variance of position estimate can achieve. That is, it can be regarded as the representation of localization accuracy. However, since CRLB is a matrix with very complex entries, it is difficult to use and to analyze in the evaluation of localization accuracy. Therefore, a scalar metric will be proposed in the next subsection, which can replace the complicated CRLB and approximately represent the localization accuracy.
3.2. A Novel Metric
Geometrically, CRLB defines an ellipse of estimation errors in both -axis and -axis . Let and denote the eigenvalues of CRLB, and consequently and represent the length of -axis and -axis of the ellipse . Then, the area of the ellipse is , which reflects the localization accuracy. That is, a smaller ellipse area indicates a more accurate estimation of target position. For simplicity, the constant is dropped and let be the metric of localization accuracy. Note that where is the intersection angle between the line of sensor to target and the line of sensor to target. Consequently, taking (17) and (18) into (16) yields
Remark 3. In some previous studies [10, 19], is also used as a metric, which is the so-called D-Criterion. Since it is monotonic to , it can reflect the variation trend of localization accuracy to some extent. However, it is important to notice that the D-Criterion is nonlinear to the area of ellipse of estimation error; thus, the scale to quantify the localization accuracy will be distorted. In assessing the average localized accuracy, it would be useful to have a scale which is linear to the area of ellipse. This is achieved by using as the metric.
Remark 4. From , of additive noise measurement model is It is clear that (20) is a special case of (19). For a localization system with only additive noise, multiplicative noise does not exist; that is, , . Therefore, from (14) and retreats to a more simple form that only depends on the intersection angle .
3.3. Analytical Properties of
Analytical properties of will be derived in this part, which indeed reveal very interesting results about the relationships between and the target-sensor distance, and between and the noise. In the following paragraphs, we will first derive the relationship between and and then the relationship between and and at last the relationship between and . The first property can be directly observed from (19).
Property 1. is a strictly monotonically decreasing function of .
Property 2. The relationship between and target-sensor distance depends on the mean and variance of multiplicative noise as follows.(i)If , is a strictly monotonically decreasing function of .(ii)If , is a strictly monotonically increasing function of when , and is a strictly monotonically decreasing function of when . achieves the maximum value at , where
Proof. The proof of Property 2 is straightforward. The partial derivative of with respect to can be written as follows: The denominator of is always positive. If , the numerator of is negative; then and is a strictly monotonically decreasing function of . Otherwise, if , the monotonicity of is related to . When , , thus is a strictly monotonically increasing function of ; when , , thus is a strictly monotonically decreasing function of . achieves the maximum at . The proof is completed.
Figure 1 illustrates the profiles of over with different . In accordance with the two aforementioned properties, the relationship between and can be readily conducted.
Remark 5. Due to the existence of multiplicative noise, becomes more complicated and is affected by both target-sensor distance and intersection angle. For the localization system with only additive noise, it is widely accepted only that the good angular diversity can lead to a low , namely, an accurate estimation of target position . However, for the localization system with both additive and multiplicative noises according to Property 2, both a good angular diversity and a proper distance are necessary to yield the best localization accuracy. Surprisingly, a shorter distance does not necessarily lead to a better accuracy, and will yield the best accuracy under a certain circumstance.
Property 3. (i) The relationship between and additive noise: is a strictly monotonically decreasing function of .
(ii) The relationship between and multiplicative noise:(a)if , is a strictly monotonically decreasing function of ;(b)if , is a strictly monotonically increasing function of when , and is a strictly monotonically decreasing function of when , where
Proof. The proof of Property 3 is straightforward.(i)The partial derivative of with respect to is
Then is a strictly monotonically decreasing function of .(ii)The partial derivative of with respect to can be written as follows:
In (25), is always positive. If , is negative; then and is a strictly monotonically decreasing function of . Otherwise, if , the monotonicity of is related to . When , , thus is a strictly monotonically increasing function of ; when , , thus is a strictly monotonically decreasing function of .
The proof is completed.
Remark 6. From Property 3(ii), surprisingly, we can conclude that the increase of multiplicative noise may not necessarily reduce the localization accuracy. When , with the further increase of , will also increase and the accuracy will be improved. However, two points should be noticed. First, ; therefore, is upper bounded by the constant even if is infinite. Second, in practice the target-sensor distance is much larger than ; therefore, always holds and case (b) barely exists in applications.
4. Numerical Analysis of under Different Sensor Deployment Patterns
In practice, it is usually desirable for target localization systems to deploy sensors in regular patterns. The regular deployment patterns outperform other irregular ones in many aspects, such as the convenience of deployment, a high degree of coverage, and connectivity . Nevertheless, the localization accuracy under different regular patterns is rarely discussed. To this end, the localization accuracy of range-only localization system with both additive and multiplicative noises will be studied in this section. Three most important regular patterns will be investigated, that is, triangle, square, and hexagon, which indeed are the basic elements to form a large localization system. From (13), (18), and (19), it is clear that is a complicated function of noise, intersection angle , and target-sensor distance . In addition, for a given deployment pattern the sensors are fixed; thus, both and are complicated functions of the target position . That is, becomes an extremely complicated function of noise and . As a result, theoretical analysis of with respect to the variations of noise and seems impossible. Therefore, numerical study will be conducted in this section. Due to page limitation, the numerical experiments will be carried out only for several typical parameter settings, which, however, are able to reveal the primary properties of .
4.1. Distribution of
In this subsection, the distribution of inside the three regular patterns, that is, triangle, square, and hexagon, will be investigated. Since the effect of additive noise on is relatively simple, we will mainly focus on the impact of multiplicative noise. To this end, let the variance of additive noise be constant for simplicity, that is, , and assume that both additive and multiplicative noises are zero-mean white Gaussian noise; that is, . The work in  indicates that in practice the distance between sensors is much larger than noise disturbance; therefore, let the side length of regular patterns be . As the experimental results in  also suggest that multiplicative noise is 3 to 4 times larger than additive noise, let the variance of multiplicative noise be , and the corresponding distributions of are shown in Figures 4, 5, and 6, respectively.
In Figure 4 with only additive noise, it is clear that the lowest always locates at the inner center of the three regular patterns. That is, the inner center achieves the best localization accuracy. On the other hand, the worst localization accuracy is located near the sensors for all patterns. This is because under additive noise, the target-sensor distance is irrelevant to localization accuracy, and the point with the best angular diversity, that is, the inner center, yields the best accuracy.
Figures 5 and 6 show that the distribution of under both additive and multiplicative noises is completely different from that under additive noise. For all patterns, the points with the best localization accuracy are located at the active sensors, while the points with poor accuracy are located at the inner centers and the midpoints on each side. It is clear that, due to the existence of multiplicative noise, the target-sensor distance significantly influences the distribution of . For more details, the profiles of along some representative lines will be examined next.
Without loss of generality, the square will be studied and similar results can readily be obtained for triangle and hexagon. Let us look at the side line , the diagonal line , and midline of the square in Figure 3(b), and the corresponding profiles of are shown in Figures 7, 8, and 9. It is clear that the minimal points of are driven from the central points to sides and corners with the increase of . In addition, the shape of profile will be maintained when is larger than a very small value, which also suggests that the shape of the distribution of with other should be similar to the shape of Figures 5 and 6.
4.2. Evaluation of Average
In many localization tasks, it is required to monitor a certain area by a certain number of sensors. Consequently, there arises a question of what the best regular pattern is to deploy these sensors. To solve this problem, it is important to investigate the average localization accuracy of different patterns.
Suppose that, in a given area , the total number of deployed sensors is ; then, the area that a single sensor can cover is . For regular deployment patterns, satisfies  where is the area of a regular pattern, stands for the number of sensors that a pattern contains, and denotes the number of patterns that a sensor connects. As we can see in Figure 10, , for triangle, , for square, and , for hexagon.
Let their side lengths be , , and , respectively. Then, we have where , , and denote the area of triangle, square, and hexagon, respectively. Consequently, the side length of each regular pattern can be derived as follows:
Let denote the average inside a certain regular pattern. Figure 11 illustrates the profiles of of three regular patterns under different situations, where and the -axis is and -axis is . The left-top subfigure shows the situation with only additive noise. It is clear that of each pattern remains constant despite the increase of , because depends only on the angular diversity other than the target-sensor distance. On the other hand, other subfigures confirm that increases with the increase of under both additive and multiplicative noises, which is consistent with our common sense. Figure 11 also indicates that the hexagon and the triangle always yield the lowest and highest , respectively, under all situations. That is, deploying the sensors in a hexagon pattern will lead to the best localization accuracy.
However, though the hexagon pattern obtains the lowest , it does not mean that the hexagon pattern is the best way to deploy sensors in practice. It is notable that, for the three deployment patterns, the numbers of active sensors are different, that is, , , , which leads to different cost, including energy, communication, and computation. Therefore, a new index is desired to represent the normalized average localization accuracy that a single active sensor can contribute. As the lower corresponds to the better accuracy, this index can be given as The comparison of of these regular patterns is illustrated in Figure 12. From the left-top subfigure with only additive noise, the hexagon pattern has the smallest ; that is, deploying the sensors in hexagon will yield the best normalized average accuracy. However, if both additive and multiplicative noises are considered and the variance of multiplicative noise is comparable with that of additive noise, the results are completely different. That is, the triangle pattern owns the best normalized average localization accuracy, while the hexagon pattern yields the worst.
Remark 7. In practice, the localization accuracy and other factors may be combined together for optimal sensor deployment. Different weights can be utilized to reflect the preference of designer under different situations. However, the same methodology can be applied to obtain the tradeoff between localization accuracy and other concerns.
5. Experimental Verification
A lab-scale experiment is set up to illustrate, first, the existence of additive noise and multiplicative noise and, then, the profiles of of a moving target. The ultrasonic range-only sensors are used in this experiment, which consists of three modules: core microprocessor Atmegal128L, wireless communication module CC2420, and ultrasonic sensor module URM37 V3.2.
In order to show the existence of additive and multiplicative noises and compute the corresponding means and variances, a 0.05 m × 0.05 m × 0.1 m rectangular target is placed in the front of an ultrasonic range-only sensor from 0.1 m to 2 m for every 0.1 m. At each point, 600 distance readings are recorded. Based on (3) and (4) and following the least square fitting procedure detailed in , we can readily obtain for additive noise , and for multiplicative noise , . It is clear that the variance of multiplicative noise is even larger than that of additive noise. However, the sensor is not unbiased. Figure 13 shows 300 distance readings when the true target-sensor distance is 1.2 m. We can see that the existence of multiplicative noise and additive noise will bring biases to the distance measurements.
In order to show the profiles of of a moving target, we set up a 1.2 m × 1.2 m square deployment with 4 ultrasonic range-only sensors for the ease of manipulating, as shown in Figure 3(b). The experimental testbed is shown in Figure 14. Now, we will move a round target with diameter of 0.1 m along diagonal line and midline , respectively, for every 0.1 m. At each point, 250 sets of distance readings will be recorded. For each set of readings, it consists of the distance measurements of 4 sensors, which can be used to readily estimate the target location. Consequently, the empirical variance matrix of localization error can be obtained for each point, as well as its determinant that is plotted in Figures 15 and 16. Meanwhile, also can be computed based on (16) for each point and plotted in Figures 15 and 16, with and . It is clear that the experimental profiles are close to the theoretical ones. It is also noticed that the determinant of variance matrix is constantly larger than for all points. Because is derived from CRLB that is the lower bound of any unbiased estimator, only reflects the lower bound of unbiased localization. In addition, the ultrasonic range-only sensor itself is biased, which may also increase the variance.
This paper studies the issue of localization accuracy for range-only sensors with both additive and multiplicative noises. The CRLB of the localization system and a novel metric of CRLB are derived. Then, three analytical properties of are derived, which are helpful for us to understand the underlying relationships between target-sensor distance and localization accuracy, as well as between noise and localization accuracy. In particular, two unexpected results about the relationships are found.
Then, the localization accuracy for three regular deployment patterns, including triangle, square, and hexagon, is evaluated numerically in terms of . It shows that the existence of multiplicative noise will significantly change the geometric distribution of localization accuracy. The average localization accuracy is also examined, and the results indicate that the hexagon always has the best average accuracy, no matter whether multiplicative noise exists or not. However, if the average accuracy is normalized by the number of active sensors, the results suggest that the triangle would be the best regular pattern to deploy the sensors.
At last, lab-scale experiment is carried out to illustrate the existence of additive and multiplicative noises. It is also shown that the experimental localization accuracy is close to , which nicely verifies the previous theoretical analysis.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work was supported by 973 Program under Grant 2013CB329503 and NSFC under Grant 61371159.
- B. S. Yaakov, X. Li, and T. Kirubarajan, Estimation with Application to Tracking and Navigation, Wiley, New York, NY, USA, 2001.
- X. Wang, M. Fu, and H. Zhang, “Target tracking in wireless sensor networks based on the combination of KF and MLE using distance measurements,” IEEE Transactions on Mobile Computing, vol. 11, no. 4, pp. 567–576, 2012.
- Y. Cheng, X. Wang, T. Caelli, X. Li, and B. Moran, “Optimal nonlinear estimation for localization of wireless sensor networks,” IEEE Transactions on Signal Processing, vol. 59, no. 12, pp. 5674–5685, 2011.
- L. Dong, “Cooperative localization and tracking of mobile ad hoc networks,” IEEE Transactions on Signal Processing, vol. 60, no. 7, pp. 3907–3913, 2012.
- N. Kantas, S. S. Singh, and A. Doucet, “Distributed maximum likelihood for simultaneous self-localization and tracking in sensor networks,” IEEE Transactions on Signal Processing, vol. 60, no. 10, pp. 5038–5047, 2012.
- G. Gaderer, P. Loschmidt, A. Nagy, R. Exel, and T. Sauter, “Localisation in wireless sensor networks,” IEEE Sensors, vol. 22, no. 4, pp. 54–69, 2005.
- S. Martínez and F. Bullo, “Optimal sensor placement and motion coordination for target tracking,” Automatica, vol. 42, no. 4, pp. 661–668, 2006.
- F. Zhao, J. Shin, and J. Reich, “Information-driven dynamic sensor collaboration,” IEEE Signal Processing Magazine, vol. 19, no. 2, pp. 61–72, 2002.
- H. Wang, K. Yao, and D. Estrin, “Information-theoretic approaches for sensor selection and placement in sensor networks for target localization and tracking,” Journal of Communications and Networks, vol. 7, no. 4, pp. 438–449, 2005.
- A. N. Bishop, B. Fidan, B. D. O. Anderson, K. Doğançay, and P. N. Pathirana, “Optimality analysis of sensor-target localization geometries,” Automatica, vol. 46, no. 3, pp. 479–492, 2010.
- L. M. Kaplan, “Global node selection for localization in a distributed sensor network,” IEEE Transactions on Aerospace and Electronic Systems, vol. 42, no. 1, pp. 113–135, 2006.
- B. Ristic and M. S. Arulampalam, “Tracking a manoeuvring target using angle-only measurements: algorithms and performance,” Signal Processing, vol. 83, no. 6, pp. 1223–1238, 2003.
- M. Steven Kay, Fundamentals of Statistical Signal Processing: Estimation Theory, vol. 1, 1993.
- Z. Yang, C. W. Chan, and Y. Wang, “A high-accuracy detection and estimation method of intermodulated sinusoids,” IEEE Transactions on Circuits and Systems I, vol. 58, no. 10, pp. 2477–2484, 2011.
- Z. Yang and C. W. Chan, “On estimation of intermodulated frequencies for sinusoidal signals with unknown fundamental frequencies,” IEEE Transactions on Signal Processing, vol. 57, no. 8, pp. 3279–3283, 2009.
- Z. Yang and C. W. Chan, “Simultaneous estimation of the input and output frequencies of nonlinear systems,” Automatica, vol. 44, no. 7, pp. 1822–1830, 2008.
- H. Godrich, A. P. Petropulu, and H. V. Poor, “Sensor selection in distributed multiple-radar architectures for localization: a knapsack problem formulation,” IEEE Transactions on Signal Processing, vol. 60, no. 1, pp. 247–260, 2012.
- T. Zhao and A. Nehorai, “Information-driven distributed maximum likelihood estimation based on Gauss-Newton method in wireless sensor networks,” IEEE Transactions on Signal Processing, vol. 55, no. 9, pp. 4669–4682, 2007.
- D. Uciński, “Optimal sensor location for parameter estimation of distributed processes,” International Journal of Control, vol. 73, no. 13, pp. 1235–1248, 2000.
- K. Doǧançay and H. Hmam, “Optimal angular sensor separation for AOA localization,” Signal Processing, vol. 88, no. 5, pp. 1248–1260, 2008.
- J. N. Ash and R. L. Moses, “On the relative and absolute positioning errors in self-localization systems,” IEEE Transactions on Signal Processing, vol. 56, no. 11, pp. 5668–5679, 2008.
- M. Zoghi and M. H. Kahaei, “Adaptive sensor selection in wireless sensor networks for target tracking,” IET Signal Processing, vol. 4, no. 5, pp. 530–536, 2010.
- J. M. Francos, “Cramér-Rao bound on the estimation accuracy of complex-valued homogeneous Gaussian random fields,” IEEE Transactions on Signal Processing, vol. 50, no. 3, pp. 710–724, 2002.
- Z. Ben-Haim and Y. C. Eldar, “The Cramér-Rao bound for estimating a sparse parameter vector,” IEEE Transactions on Signal Processing, vol. 58, no. 6, pp. 3384–3389, 2010.
- X. Bai, S. Kumar, D. Xuan, Z. Yun, and T. H. Lai, “Deploying wireless sensors to achieve both coverage and connectivity,” in Proceedings of the 7th ACM International Symposium on Mobile Ad Hoc Networking and Computing (MOBIHOC '06), pp. 131–142, May 2006.
Copyright © 2014 Xiufang Shi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.