Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2018 / Article

Research Article | Open Access

Volume 2018 |Article ID 8419384 | 10 pages |

Navigation for Indoor Robot: Straight Line Movement via Navigator

Academic Editor: Simone Bianco
Received04 Jul 2018
Accepted02 Sep 2018
Published24 Sep 2018


Due to the need of zigzag overlay strategy, long-term linear motion is essential for sweep robot. However, the existing indoor sweep robot navigation algorithm has many problems; for instance, algorithm with high complexity demands high hardware performance and is incapable of working at night. To overcome those problems, in this paper, a new method for indoor robot Straight Line Movement via Navigator (SLMN) is proposed to ensure long linear motion of robot with an acceptable error threshold and realize multiroom navigation. Firstly, in a short time, robot runs a suitable distance when it is covered by navigator’s ultrasonic sensor. We can obtain a triangle with twice the distance between navigator and robot and the distance of robot motion. The forward angle of the robot can be conveniently obtained by the trigonometric functions. Comparing the robot’s current angle with expected angle, the robot could correct itself and realize the indoor linear navigation. Secondly, discovering dozens of the magnitude gaps between the distance of robot run and the distance between navigator and robot, we propose an optimized method using approximate scaling which increases efficiency by nearly 70.8%. Finally, to realize multiroom navigation, we introduce the conception of the depth-first search stack and a unique encode rule on rooms and navigators. It is proved by extensive quantitative evaluations that the proposed method realizes indoor full coverage at a lower cost than other state-of-the-art indoor vision navigation schemes, such as ORB-SLAM.

1. Introduction

Traditional low cost sweep robot used the random mode of sweeping, with small coverage area, low cleaning efficiency, and easy collision. Adopting the shape of Z full coverage algorithm [1] can effectively reduce the disadvantage of random sweeping, but robot must run linearly if the full coverage algorithm can work steadily [2, 3]. It could work normally based on an inertial measurement unit (IMU) [4, 5], but with the data from both sensor and PID correcting its direction, it is not enough to realize its indoor localization. Due to accumulated errors of sensors, wheel mechanical inaccuracy, and the slide problem, sweeping robot would drift away, so researchers start the indoor navigation system [6] research.

Rodas et al. proposed an optimal Bayesian fusion and localization algorithm based on Bluetooth positioning data [7]. However, Bluetooth localization error is still serious, deployment cost is huge, and the result is not satisfactory. Lin et al. proposed a fingerprint localization system based on Wi-Fi [8], which is a sample of signal intensity, and preconstruct a radio map, use the information of Wi-Fi cover, and read data from accelerometer reliably. However, a great quantity initial work should be finished before using it, and the effect is not very well. Fuentes-Pacheco et al. proposed vSLAM (visual Simultaneous Localization and Mapping) [913]; this scheme has lower deployment cost and high localization accuracy features, and also its visualization is prefect. However, when the light had changed, the result fluctuated greatly. With introducing a Kinect (a depth camera) [14], the cost of full system would rise. In [15], Binu proposed to use the naturally occurring magnetic fields inside the concrete buildings with simplified magnetic map [16] matching algorithm. Due to the existence of sensor accumulation error problems, the navigation’s accuracy is insufficient. Furthermore, Wolfgang Hess et al. propose a 2D SLAM system that combines scan-to-submap matching with loop closure detection and graph optimization, making its source codes open-source [17]. However, equally it needs a LIDAR, which costs much money.

Because of the characteristics of light moving straight forward, the infrared ray [18] wavelength is long and mostly would not scatter in the air. An infrared ray is sent by opening a small hole in the front of the navigator, so that it can be considered as a line in terms of the ultrasonic (see Figure 1). The infrared communication module consists of transmitter and receiver. The transmitter is located on the navigator, and the receiver is located on the sweeping robot. A low scattering infrared ray is sent by the transmitter, after being received by the sweeping robot, and it can be regarded as a baseline in the room (we usually put the navigator approximately on the middle of the room’s edge) [19].

Similarly, the ultrasonic module is divided into the transmitter and the receiver. The transmitter is located on the navigator, and the receiver is located on the sweeping robot. Measuring the time of running between the navigator and the sweeping robot and multiplying by the velocity of ultrasonic in the air, we can work out the distance between them. While meeting the infrared baseline, the sweeping robot would start to correct itself. And because the infrared ray has the reflection characteristic, when the infrared ray radiates to the ceiling, the reflection will occur, so that the full room would be covered [20].

In a word, to realize indoor linear navigation, the deployment cost based on Wi-Fi [21, 22] and Bluetooth is huge, and vSLAM cannot work normally in the evening. Therefore, when considering the cost problem with guaranteeing high coverage percentages, we propose to use a navigator (consisting of infrared and ultrasonic) to help localization and navigation. It mainly carries out the shape of Z strategy based on motor encoder and PID controller local optimum [23]. At the same time, infrared ray and ultrasonic are used to realize auxiliary localization and correct itself to make global optimum [18]. Comparing to existing algorithms, extensive quantitative evaluations shows that the proposed method realizes indoor full coverage at a lower cost than other state-of-the-art indoor vision navigation schemes, such as ORB-SLAM [11, 24] (see Figure 2).

The remainder of this paper is organized as follows. Section 2 provides an overview of the robot hardware platform. Section 3 explains the proposed localization and navigation algorithm. Finally, Section 4 provides experimental results and comparisons.

2. Robotic Platform

In order to verify the navigation algorithm, this paper designs a double navigation platform between sweeping robot and navigator. The sweeping robot adopts Cortex-M3 architecture 32 bit microcontroller stm32 as the main control chip, which is produced by STMicroelectronics. The navigator adopts super low power msp430 as the main control chip, which is produced by the United States TI. A sweeping robot consists of two wheels controlled by deceleration motors, equipped with photoelectric encoders to measure the speed of the two wheels, and uses six-axis gyroscope MPU6050 to calculate its posture [25]. Additionally, the sweeping robot and navigator are equipped with infrared transceiver module to communicate with each other at the top [20]. In flank, an infrared sensor and ultrasonic sensor are used to obtain the distance between them as shown in Figure 3.

3. Localization and Navigation of Indoor Robot

The navigator is located at the edge of the room, and it regularly sends an infrared ray signals without scattering. When the sweeping robot detects infrared ray signal, calculates the distance , and records the wheels roll over the distance , the ultrasonic wave starts to emit regular frequency. With the ultrasonic wave, we could measure again. As a result, by solving which form a vector triangle, we can obtain the real forward angle and correct the sweeping robot direction by comparison to the theoretical value as shown in Figure 4.

3.1. Ensure Running in a Straight Line

Traditionally, a sweeping robot can stably execute straight line because two wheels can keep themselves running with the same velocity. Because the electronic characteristics of two wheels may not be completely consistent, in the same PWM duty rate, the velocities of two motors are not exactly the same. If there is long-term running, the error is accumulated continually, and the sweeping robot would move forward one side. Therefore, it is necessity to adjust the velocity of the two wheels in real time so that the velocity of two wheels becomes the same, and robot can keep itself linearly running for a long time without any consideration to the other obstacles.

In this paper, we installed two 26-hole photoelectric tachometer encoders at the back of the motors. By measuring the number of pulses generated by the photoelectric encoder in unit time and before that measuring the circumference of the wheel L, the formula is given as follows.

Theoretically, after obtaining the forward speed of V, the velocity of the wheel can be adjusted by the PID control algorithm in the field of classical automatic control. However, considering the approximate errors produced on the formula transformation, we directly added the number of pulses N generated by the photoelectric encoder to the PID formula [26], instead of calculating the forward speed of V each wheel. Generally, there are two kinds of PID, as shown below.

Incremental PID control algorithm:

where is the output value of PID and is the desired value. is the photoelectric encoder current measurement value, and is the deviation value. In this paper, DC motor is used for robotic motion, by controlling the duty rate of PWM to control the velocity, so according to the incremental PID control algorithm, the recursion formula is

where is currently the value of PWM, which supplies motor; additionally, is the last value of PWM. , , are the offsets between expected velocity and real velocity by current, by last time, and by the time before last time. T is the sample time, usually a fixed value; ,, are the parameters which will be optimized. The algorithm’s pseudocode is as in Algorithm 1.

1 Do
2 ifn=1 or 2
3  calculate the offsets of right and left wheel ,
4  calculate the current PWM value of right and left wheel ,
5  substitute n=1 or 2 into equation(6), obtain , , ,
6 if n=3
7  substitute , , , into equation(6)
8 While (if work signal = true)
9 return,, ,,

Among them, are the PID output values and offsets of the left and right wheel, respectively. Using PID algorithm, the mechanical errors can be eliminated, and the robot can ensure that it is straight as much as possible.

3.2. Positioning and Navigation Algorithm
3.2.1. Angle and Position Correction

The two wheels of the sweeping robot are driven by two independent motors. Although using PID regulation control can ensure that the two wheels are at the same speed, the robot would shift when running a distance because there may exist slipping and mechanical error. A new method of angle and position correction based on infrared and directional ultrasonic sensors is presented in this paper as shown in Figure 5.

The navigator is placed in the middle of the wall on one side of the room, emitting vertical infrared ray and fan-shaped ultrasonic waves on both sides. The angle of motion of the robot is always parallel to the wall. When the robot runs from top to bottom and from right to left, and each time is in the middle of the room, the robot receives infrared ray and ultrasonic waves, determining its position by combining onboard devices. If the robot discovers deviation from the straight line, it corrects itself in time.

Basic work principle: With the infrared ray by the navigator treated as the baseline, we sent an ultrasonic wave with certain angles in particular frequency which can measure the precise position of robot. Based on two times of measurement, we can resolve robot’s deviation angle. Specific steps are as follows (before the algorithm is detailed, note that we briefly define the notation that is used throughout the paper in Table 1):


OThe position of navigatorAThe Expected position of robot
BThe real position of robot
CThe point where circle whose radius is equal to the distance of O to A cuts OB
The vector of O to AThe distance of O to A
The vector of O to BThe distance of O to B
The vector of A to BThe distance of A to B
The vector of O to CThe distance of O to C
The vector of A to CThe distance of A to C

Step 1. When the robot arrives at the baseline of infrared, after it receives an infrared ray, the distance between them can be obtained.

Step 2. After the robot leaves the baseline and runs in a short distance , measure the distance between them by ultrasonic wave again.

Step 3. Resolve the horizontal attitude angle from the relationship of geometry.

However, as the complexity of the formula is expensive and the gap between and is more than orders of magnitude, we propose a geometry algorithm to make the complexity reduction substantial (see Figure 6).

Follow as the center of a circle, the direction of as x axis forward, and the direction of perpendicular as y axis forward. Let down trajectory is y=f(x). The slope of the tangent at C is in .

The slope of AC: because of towards 0 relative to , -=o(x) [27], o(x) is the infinitely small quantity of, the slope of the tangent at C being in O is equal to the slope of AC.

In the circle, from the relationship of geometry (1,)=0, then =0.

Similarly, =0, from the relationship of geometry =+.

As the (10) and (12) are antitrigonometric functions, we can directly consider the internal formula. As we all know, multiplication is nearly 10 times slower than addition and subtraction, and division is nearly 20 times slower than addition and subtraction, when it is an 8-bit CPU. Now, we define the time of addition and subtraction as , then the multiplication is , and division is . For (10), the time of the whole equation is as follows.For equation (12), the time of the whole equation is as follows.

Then, we have the following.

Compared with (10), (12) reduces the computational complexity [28] and increases efficiency by nearly 70.8%. At the same time, it satisfies the actual precision requirement under the condition of being smaller. The relationship between relative errors and deviation angle from the original line is shown in Figure 7.

3.2.2. Navigation Implementation between Multiple Rooms

In order to help a robot to cross different rooms to work continuously, it should be improved autonomously. In this paper, we built topology model for multiple rooms based on the theory of graph optimization [29]. Meanwhile, the depth-first search (DFS) [30] stack is introduced to firstly select the most nearest and inmost room to work [31], which improves the efficiency of DFS. The algorithm’s pseudocode is as in Algorithm 2.

1 Define every room has a unique ID
2 Let S be a stack
3 if robot receive an infrared ray communication signals
4  Start work in regular path
5  if the robot receive another infrared ray communication signals
6   Then S.push()
7   When finish working this room
8    While S is not empty
9      = S.pop()
10      Enter the ID of room
11       if is not labeled as discovered:
12       Label as discovered
13      Goto line 3
14      Back to initial room, then the work is over
15 Else
16  Stop work

In this algorithm, we define a room as by the navigator which has a unique encode signal. We build a stack to store the room information so that the robot can select the most nearest and inmost room to work. When it finishes the current room, then it checks the stack status every time. Lastly, it goes back to the initial room and stops working.

4. Implementation and Experiments

We tested the SLMN algorithm in a m real indoor environment, and the dimension of the robot was m. Since it is an upgraded TF card, the robot could easily store the data from the electronic motor encoder. The real robot was shown in Figure 8.

After finishing its work run, the local trajectory of the robot was restored (as the errors are small, it is difficult to discover the trajectory of at a global view), as shown in Figure 9. The red and blue lines represent the robot’s true trajectory and ideal trajectory, respectively. Therefore, we can find out that when the robot passed an infrared reference line (m) by the navigator, the deviation trajectory would be effectively corrected.

By calculating the absolute value, a good view of this deviation could be created. The deviation that appears on will be effectively corrected at as shown in Figure 10.

Occasionally, resources such as memory and computational time are of primary concern, but most often cost and accuracy are what we want to measure. By analyzing several candidate algorithms for the navigation problem [6], we can identify a most efficient one. In this paper, we put the results of comparative analysis with other state-of-the-art algorithms such as vSLAM [11], LIDAR SLAM [17], and inertial navigation [32, 33] as shown in Table 2.

Compared Algorithms Precise and cost
Mean relative errors in (%)Battery life with

Radar SLAM
Inertial navigation

From Table 2, it is seen that the two SLAM systems have similar precision, which are the most highest. However, their work time is also the shortest in the same battery capacity, so it needs more power to keep it working. Although inertial navigation system has the lowest cost, its error is too enormous to work normally. As a result, since the development of new energy still faces several uncertain difficulties [3436], our proposed algorithm will cost less energy in satisfying practical applications. which than the other state-of-the-art algorithms.

5. Conclusion

On the one hand, this paper proposed a high efficiency method on indoor sweep robot line navigation. The two measured values between navigator and robot and the distance of robot motion can form a triangle then to solve the forward angle and correct the direction of forward speed with analytic geometry principle. Finally, we implemented the long time line navigation, and based on that, we propose multiroom navigation algorithm. After simulation and reality environment experiments, on the threshold of error, an indoor line navigation method based on the navigator which consists of the infrared and ultrasonic sensors, even in the evening, could resolve the problems on the other existing navigation.

On the other hand, we prepare to do more works on building a map based on the high efficiency motion to improve its practicability and integrity, so it can become self-charging and work continuously.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


This work was supported by National Key R&D Program of China Nos. 2018YFC0806900, 2016YFC0800606, 2016YFC0800310, and 2018YFC0407905; Natural Science Foundation of Jiangsu Province under Grants Nos. BK20150721 and BK20161469; and Primary Research & Development Plan of Jiangsu Province under Grants Nos. BE2015728, BE2016904, BE2017616, and BE2018754.


  1. L. Jia-jun, Z. Bi, and H. Yuan-lie, “Research on Path Planning Algorithm for Cleaning Robot Based on Improved Potential Field Grid Method,” Journal of Guangdong University of Technology, 2016. View at: Google Scholar
  2. J. Shang, X. Hu, F. Gu, D. Wang, and S. Yu, “Improvement schemes for indoor mobile location estimation: a survey,” Mathematical Problems in Engineering, vol. 2015, Article ID 397298, 32 pages, 2015. View at: Publisher Site | Google Scholar
  3. G. Liu, Y. Geng, and K. Pahlavan, “Effects of calibration RFID tags on performance of inertial navigation in indoor environment,” in Proceedings of the 2015 International Conference on Computing, Networking and Communications, ICNC 2015, pp. 945–949, USA, February 2015. View at: Google Scholar
  4. B.-W. Li and D.-Y. Yao, “Low-cost MEMS IMU navigation positioning method for land vehicle,” Zhongguo Guanxing Jishu Xuebao/Journal of Chinese Inertial Technology, vol. 22, no. 6, pp. 719–723, 2014. View at: Google Scholar
  5. S. Lynen, M. W. Achtelik, S. Weiss, M. Chli, and R. Siegwart, “A robust and modular multi-sensor fusion approach applied to MAV navigation,” in Proceedings of the 2013 26th IEEE/RSJ International Conference on Intelligent Robots and Systems: New Horizon, IROS 2013, pp. 3923–3929, Japan, November 2013. View at: Google Scholar
  6. Y. Zheng, G. Shen, L. Li, C. Zhao, M. Li, and F. Zhao, “Travi-Navi: self-deployable indoor navigation system,” IEEE/ACM Transactions on Networking, no. 99, pp. 1–15, 2014. View at: Google Scholar
  7. J. Rodas, C. J. Escudero, and D. I. Iglesia, “Bayesian filtering for a Bluetooth positioning system,” in Proceedings of the 2008 IEEE International Symposium on Wireless Communication Systems, ISWCS'08, pp. 618–622, Iceland, October 2008. View at: Google Scholar
  8. P. Lin, Q. Li, Q. Fan, X. Gao, and S. Hu, “A real-time location-based services system using WiFi fingerprinting algorithm for safety risk assessment of workers in tunnels,” Mathematical Problems in Engineering, vol. 2014, Article ID 371456, 10 pages, 2014. View at: Publisher Site | Google Scholar
  9. J. Fuentes-Pacheco, J. Ruiz-Ascencio, and J. M. Rendón-Mancha, “Visual simultaneous localization and mapping: a survey,” Artificial Intelligence Review, vol. 43, no. 1, pp. 55–81, 2015. View at: Publisher Site | Google Scholar
  10. S. Leutenegger, S. Lynen, M. Bosse, R. Siegwart, and P. Furgale, “Keyframe-based visual-inertial odometry using nonlinear optimization,” International Journal of Robotics Research, vol. 34, no. 3, pp. 314–334, 2015. View at: Publisher Site | Google Scholar
  11. R. Mur-Artal and J. D. Tardos, “ORB-SLAM2: An Open-Source SLAM System for Monocular, Stereo, and RGB-D Cameras,” IEEE Transactions on Robotics, vol. 33, no. 5, pp. 1255–1262, 2017. View at: Publisher Site | Google Scholar
  12. S. Yin, X. Li, H. Gao, and O. Kaynak, “Data-based techniques focused on modern industry: an overview,” IEEE Transactions on Industrial Electronics, vol. PP, no. 99, p. 1, 2014. View at: Publisher Site | Google Scholar
  13. R. A. Newcombe, S. J. Lovegrove, and A. J. Davison, “DTAM: Dense tracking and mapping in real-time,” in Proceedings of the 2011 IEEE International Conference on Computer Vision, ICCV 2011, pp. 2320–2327, Spain, November 2011. View at: Publisher Site | Google Scholar
  14. J. Smisek, M. Jancosek, and T. Pajdla, “3D with Kinect, Advances in Computer Vision Pattern Recognition,” 3D with Kinect, Advances in Computer Vision Pattern Recognition, vol. 21, no. 5, pp. 1154–1160, 2013. View at: Google Scholar
  15. P. K. Binu, R. A. Krishnan, and A. P. Kumar, “An efficient indoor location tracking and navigation system using simple magnetic map matching,” in Proceedings of the 2016 IEEE International Conference on Computational Intelligence and Computing Research (ICCIC), pp. 1–7, Chennai, December 2016. View at: Publisher Site | Google Scholar
  16. A. Sud, R. Gayle, E. Andersen, S. Guy, M. Lin, and D. Manocha, “Real-time navigation of independent agents using adaptive roadmaps,” in Proceedings of the ACM SIGGRAPH 2008 Classes, 2015. View at: Google Scholar
  17. W. Hess, D. Kohler, H. Rapp, and D. Andor, “Real-time loop closure in 2D LIDAR SLAM,” in Proceedings of the 2016 IEEE International Conference on Robotics and Automation, ICRA 2016, pp. 1271–1278, Sweden, May 2016. View at: Publisher Site | Google Scholar
  18. K. Wang, T. Song, T. Liang et al., “A dual-infrared-transmitter optical wireless based indoor user localization system with high accuracy,” in Proceedings of the 2017 Optical Fiber Communications Conference and Exhibition, OFC 2017, USA, March 2017. View at: Google Scholar
  19. J. Krejsa and S. Vechet, “Infrared beacons based localization of mobile robot,” Elektronika ir Elektrotechnika, no. 1, pp. 17–22, 2012. View at: Google Scholar
  20. J. Yun and M.-H. Song, “Detecting direction of movement using pyroelectric infrared sensors,” IEEE Sensors Journal, vol. 14, no. 5, pp. 1482–1489, 2014. View at: Publisher Site | Google Scholar
  21. D. Han, S. Jung, M. Lee, and G. Yoon, “Building a practical wi-fi-based indoor navigation system,” IEEE Pervasive Computing, vol. 13, no. 2, pp. 72–79, 2014. View at: Publisher Site | Google Scholar
  22. C. Yang and H.-R. Shao, “WiFi-based indoor positioning,” IEEE Communications Magazine, vol. 53, no. 3, pp. 150–157, 2015. View at: Publisher Site | Google Scholar
  23. P. Shah and S. Agashe, “Review of fractional PID controller,” Mechatronics, vol. 38, pp. 29–41, 2016. View at: Publisher Site | Google Scholar
  24. R. Mur-Artal, J. M. M. Montiel, and J. D. Tardos, “ORB-SLAM: a versatile and accurate monocular SLAM system,” IEEE Transactions on Robotics, vol. 31, no. 5, pp. 1147–1163, 2015. View at: Publisher Site | Google Scholar
  25. J. Simanek, M. Reinstein, and V. Kubelka, “Evaluation of the EKF-based estimation architectures for data fusion in mobile robots,” IEEE/ASME Transactions on Mechatronics, vol. 20, no. 2, pp. 985–990, 2015. View at: Publisher Site | Google Scholar
  26. K. J. Astrom, PID Controllers: Theory, Design and Tuning, Instrument Society of America Research Triangle Park Nc, 1995.
  27. A. Banner, “The Calculus Lifesaver,” 2009. View at: Google Scholar
  28. T. Hcormen, C. Eleiserson, R. Lrivest, and C. Stein, Introduction to Algorithms, vol. 39, Computer Education, 3rd edition, 2013.
  29. S. R. U. Murty, Graph Theory with Applications: North Holland, 1976.
  30. K. Mehlhorn, “Engineering DFS-based graph algorithms,” 2017, View at: Google Scholar
  31. R. Tarjan, “Depth-first search and linear graph algorithms,” in Proceedings of the 12th Annual Symposium on Switching and Automata Theory (swat 1971), pp. 114–121, East Lansing, MI, USA, October 1971. View at: Publisher Site | Google Scholar
  32. A. R. Jimenez, F. Seco, J. C. Prieto, and J. Guevara, “Indoor Pedestrian navigation using an INS/EKF framework for yaw drift reduction and a foot-mounted IMU,” in Proceedings of the 2010 7th Workshop on Positioning, Navigation and Communication, WPNC'10, pp. 135–143, Germany, March 2010. View at: Publisher Site | Google Scholar
  33. B. Barshan and H. F. Durrant-Whyte, “Inertial navigation systems for mobile robots,” IEEE Transactions on Robotics and Automation, vol. 11, no. 3, pp. 328–342, 1998. View at: Publisher Site | Google Scholar
  34. K. Wang, J. Pang, L. Li, S. Zhou, Y. Li, and T. Zhang, “Synthesis of hydrophobic carbon nanotubes/reduced graphene oxide composite films by flash light irradiation,” Frontiers of Chemical Science and Engineering, vol. 12, no. 3, pp. 376–382, 2018. View at: Publisher Site | Google Scholar
  35. W. Kai, L. Liwei, X. Wen et al., “Electrodeposition synthesis of PANI/MnO2/graphene composite materials and its electrochemical performance,” International Journal of Electrochemical Science, vol. 12, no. 9, pp. 8306–8314, 2017. View at: Publisher Site | Google Scholar
  36. W. Kai, “Synthesis of porous carbon by activation method and its electrochemical performance,” International Journal of Electrochemical Science, pp. 10766–10773, 2018. View at: Publisher Site | Google Scholar

Copyright © 2018 Chaozheng Zhu 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.

More related articles

1311 Views | 383 Downloads | 1 Citation
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.