Improved sampling scheme for LiDAR in Lissajous scanning mode | Microsystems & Nanoengineering

2022-08-20 08:58:44 By : Ms. Julie Zhu

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Microsystems & Nanoengineering volume  8, Article number: 64 (2022 ) Cite this article

MEMS light detection and ranging (LiDAR) is becoming an indispensable sensor in vehicle environment sensing systems due to its low cost and high performance. The beam scanning trajectory, sampling scheme and gridding are the key technologies of MEMS LiDAR imaging. In Lissajous scanning mode, this paper improves the sampling scheme, through which a denser Cartesian grid of point cloud data at the same scanning frequency can be obtained. By summarizing the rules of the Cartesian grid, a general sampling scheme independent of the beam scanning trajectory patterns is proposed. Simulation and experiment results show that compared with the existing sampling scheme, the resolution and the number of points per frame are both increased by 2 times with the same hardware configuration and scanning frequencies for a MEMS scanning mirror (MEMS-SM). This is beneficial for improving the point cloud imaging performance of MEMS LiDAR.

3D imaging LiDAR consists of a laser, scanning mechanism and detector. At present, the scanning mechanism driven by motors has severely restricted the realization of miniaturization and high performance of LiDAR systems, which makes it difficult to meet the needs of aerospace1,2, vehicle3,4, and surveying and mapping5 applications. With the development of MEMS technology, especially the commercialization of MEMS-SM, the National Institute of Standards and Technology (NIST) proposed the next generation LiDAR technology architecture based on MEMS-SM in 20046. The US Army Research Laboratory7,8,9,10,11,12,13,14 configured LiDAR operating parameters to form images of 256 (h) × 128 (v) pixels over a 15 × 7.5° field of view (FOV)15. Yeungnam University16 developed LiDAR to measure a range image with a resolution of 848 (h) × 480 (v) point locations at the FOV of 43.2 × 24.3°, and the angular resolution is 0.0509° in simulation17. Technical University of Denmark18, KTH Royal Institute of Technology19 have long been committed to the application of MEMS-SM in 3D imaging LiDAR. However, due to the special scanning mode of MEMS LiDAR and the limitation of angular resolution and the number of points per frame, the output points of current LiDAR cannot meet the needs of lane detection20, object detection21,22, or superresolution23.

The nonresonant-driven MEMS-SM supports three scanning modes24: point-to-point scanning mode on both axes with the laser beam stopping at each angle, then stepping to the next angle, resonant scanning mode on the X-axis and quasistatic on the Y-axis, with scanning mode on both axes. All three scanning forms are nonlinear. The spectrum of a quasistatic-driving mode, such as a triangle wave, not only consists of its fundamental frequency but also contains all of its odd harmonics, which is not suitable for high-speed scanning25,26. Moreover, spiral and rosette trajectory patterns support scanning mode with both axes27. However, the Lissajous trajectory has the special geometric characteristics necessary for the rapid reconstruction of nonrectilinear Cartesian k-space trajectories with constant sampling time intervals28. The Lissajous scheme is a time-efficient sampling scheme and is currently deployed in various fields. Its application ranges from atomic force microscopy (AFM)29, microscopy30, electroencephalography reconstruction31, frequency-modulated gyroscopes32, fast eddy current testing33, augmented reality (AR)34, and magnetic resonance imaging (MRI)35 to MEMS imaging LiDAR36.

The horizontal resolution of point cloud sampling of motor-driven mechanical LiDAR is guaranteed by the rotation of the motor, and the vertical resolution is guaranteed by the installation position of the laser, so the uniform space sampling point cloud data can be obtained by equal time interval ranging. The LiDAR point cloud sampling method based on scanning mechanisms, such as galvanometers37, pyramid polygon mirrors38, and double optical wedges39 has not been specially designed. While the new LiDAR based on MEMS-SM, or the imaging equipment using MEMS-SM as the movable part, uses Lissajous scanning during imaging mode, the scan-driving function is a sine curve with two orthogonal axes, so it is difficult to realize grid sampling in time and space. Previous studies40 have found two irrelevant parameters that determine the trajectory and proposed a design method of Lissajous trajectories. Jones41 and Hoste42 analyzed the theoretical feasibility of Lissajous trajectory gridding with knots from the mathematical point of view. Likes43 proposed that the Lissajous trajectory was original, and its sampling density was studied by Hardy et al.44. Moriguchi28 demonstrated that the Lissajous curve has the feasibility of Cartesian gridding of sampling points with constant sampling time intervals, and when \(f_x = 11\,{\rm {Hz}}\) , \(f_y = 9\,{\rm {Hz}}\) , and the sampling frequency of the point cloud is 198 Hz, the point cloud satisfies the Cartesian grid. In ref. 45, nonrectilinear Cartesian grid sampling was used to capture the statistical dependence of the two remote-sensing images for registration.

To date, the inherent link among the trajectory, the sampling scheme, and the associated complexity of the remeshing process has been investigated to only a limited extent. To improve the utilization of MEMS-SM, the angular resolution and the number of points per frame of LiDAR are improved at the same hardware configuration4,46. In this paper, the sampling scheme of MEMS LiDAR in Lissajous scanning mode is studied. First, the corresponding sampling scheme is redesigned in a more compact trajectory40. Second, through the theoretical analysis of the sampling scheme satisfying the Cartesian grid, a general sampling scheme of the Lissajous trajectory pattern is proposed.

The Lissajous scanning trajectory is driven by two single-tone harmonic waveforms of constant frequency and amplitude, which are usually generated by a 2D MEMS-SM47. A schematic diagram is shown in Fig. 1.

Lissajous trajectory driven by a MEMS-SM

According to the imaging principle of LiDAR, the angular resolution of the image can be calculated, and a schematic diagram is shown in Fig. 2.

The X-axis of angular resolution (ARx), \({\rm {AR}}_x = \arctan \frac{{{\rm {d}}x}}{R}\) , where \({\rm {d}}x\) is the X-axis coordinate difference of two adjacent points (point A and point B) in space, and R is the distance between the MEMS-SM and the target. In the nonrectilinear trajectory, two adjacent points in space are usually not adjacent in time

The sample points are all on the scanning trajectory, so first, the characteristics of the Lissajous trajectory are analyzed, which is composed of two orthogonal cosine curves:

where X and Y are the horizontal and vertical coordinates of scanning points and where \(A_x\) and \(A_y\) denote the scanning amplitude of the X-axis and Y-axis directions, respectively. Here, t is time, and \(f_x,f_y,\varphi _x,\varphi _y\) are the biaxial scanning frequencies and phases of the X-axis and Y-axis directions, respectively. If \(f_x,f_y\) are both integers and a greatest common divisor (GCD) \(f_0\) exists, then the following Eq. (2) holds:

We assume that point (\(x_0,y_0\) ) is the knot of the Lissajous trajectory, which appears at \(t_0\) .

Bringing Eq. (3) into Eq. (1), we can obtain \(y_l\) in Eq. (4):

We define the phase parameter k

Therefore, \(y_l\) can be simplified to Eq. (6)

Of all the possible values of l, two unequal \(l_1\) and \(l_2\) ensure \(y_{l_1} = y_{l_2}\) . Thus, we can obtain Eq. (7)

where i is an integer and \(\arccos \left( {\frac{{x_0}}{{A_{{{\mathrm{x}}}}}}} \right) \in [0,\pi ]\) is satisfied, the range of i is \([0,n_y - 1]\) , and bringing Eq. (7) to Eq. (3), the time when the knots appear is shown in Eq. (8).

where i is an integer. Therefore, we can obtain the time interval Δ between two adjacent points in time.

A higher sampling rate implies more points and higher resolution. However, in the case of the Cartesian grid, once the trajectory is determined, the matching sampling rate is also determined, and the sampling rate cannot be changed at will. In a Lissajous trajectory period (\(1/f_0{\kern 1pt}\) ), assuming that \(f_{\rm {s}}\) is the sampling frequency, the sampling start time is \(t_1\) and the end time is \(t_2\) , the sampling scheme in ref. 28 can be summarized as Eq. (10), and the diagram of the sampling point cloud is shown in Fig. 3a.

The green line is the Lissajous scanning trajectory, and the blue dot is the spatial position of the point. The trajectories and sampling schemes are shown below. a \(k = 0,f_s = 1/\Delta ,t_1 = 0\) , b \(k = 0,f_s = 2/\Delta ,t_1 = 0\) , c \(k = 2,f_s = 2/\Delta ,t_1 = 0\) , and d \(k = 2,f_s = 2/\Delta ,t_1 = \Delta /4\)

With the trajectory of ref. 28, we increase the sampling frequency by only twice to obtain the point cloud shown in Fig. 3b. The point cloud does not satisfy the Cartesian grid. According to the trajectory design method proposed in ref. 40, when the scanning frequency ratio remains unchanged and the phase parameter k = 2 is set, the angular resolution of the Lissajous trajectory decreases. At this time, according to the sampling scheme of ref. 28, the point cloud shown in Fig. 3c is obtained. We find that the point cloud coincides with the knots of the trajectory and does not meet the Cartesian grid. Therefore, a new sampling scheme is proposed in this paper, that is, when the scanning frequency ratio is unchanged and the phase parameter k = 2. The parameters of the proposed sampling scheme are shown in Eq. (11), and the diagram of the sampling point cloud is shown in Fig. 3d.

Through the comparison between Fig. 3, we know that a perfect sampling scheme in the same scanning frequencies has three characteristics: one is the phase parameter k = 240. Two is that the sampling interval is equal to the space interval of the trajectory knots. The last is that the starting sampling time is in the middle of the 1st and 2nd knot.

The period of k is 848, and the trajectories symmetry of the first half period and the second half period. Therefore, this paper takes \(k = 0,1,2,3\) as an example. The trajectories are shown in Figs. 3a, 4a, 3d, and 4b. We find that the point cloud is exactly the same when k = 1 and k = 3, even if they have different trajectories.

a, b When the scanning frequencies of the two axes of the scanning mirror are \(f_{x} = 11\,{\bf{Hz}},{{f}}_{{y}} = 9\,{\bf{Hz}}\) , respectively, the trajectories and sampling points under different k parametric conditions. The green line is Lissajous scanning trajectory, and the blue dot is the spatial position of the sampling point. The different parametric conditions are a \(k = 1,f_s = 1/\Delta ,t_1 = 0\) and b \(k = 3,f_s = 1/\Delta ,t_1 = 0\)

By taking the four kinds of trajectories (Figs. 3a, 4a, 3d, 4b) as an example, the improved sampling scheme proposed in this paper is used to make a difference on the grid of sampling points, and the variance of these data is used as the basis for selecting trajectories. Ref. 26 is named the scanning density. As an additional finding, the adjacent knots of the trajectories in Fig. 3a and d are identical in time and space. However, the adjacent knots of trajectories in Fig. 4a and b are different in time and space.

It can be seen from Fig. 3a that when k = 0, the point cloud can be divided into 11 × 9 grids, so the difference of the X-axis is 8 points and that of the Y-axis is 10 points. Similarly, when k = 2, there are 17 points on the X-axis and 21 points on the Y-axis in Fig. 5. When k = 1 or k = 3, the difference in the X-axis is exactly the same as with k = 0; moreover, the Y-axis is exactly the same as with k = 2. The smaller the difference is, the denser the point cloud and the smaller the angular resolution.

When the scanning frequencies of the two axes of the scanning mirror are \(f_x = 11\,{\bf{Hz}},{{f}}_{{y}} = 9\,{\bf{Hz}}\) , two-axis grid difference of the different trajectories (k = 0 or k = 2) with the improved sampling scheme are shown. The X-axis represents the number of points, and the Y-axis represents the grid difference. a X-axis grid difference of two adjacent points in space. b Y-axis grid difference of two adjacent points in space

Therefore, at the same scanning frequencies of the MEMS-SM (\(f_x = 11\ {\rm{Hz}},\,f_y = 9\ {\rm{Hz}}\) ), the number of points per frame and the angular resolution of the two axes are as shown in Table 1.

In Table 1, \({\rm {AR}}_x\) is the angular resolution in the X-axis; \({\rm {AR}}_y\) is the angular resolution in the Y-axis, and its calculation method is shown in ref. 49; PPS represents the effective points per second. We can see clearly that the algorithm performance in this paper increases 2 times when the number of sampling points is 396, and the resolution of both axes is increased by 2 times.

To verify the trajectory and sampling scheme in Fig. 3d, a MEMS LiDAR prototype is designed for imaging experiments, and its structure principle is shown in Fig. 6.

a Schematic diagram of MEMS LiDAR architecture. b Photo of the experimental platform of the MEMS LiDAR. c Photo of the laser and the receiver. The experimental equipment uses NI chassis (PXI-1082) to load computer module. The pulse transmitting module (PXI-6584) generates pulses to modulate the laser (LDM635.150.A350) amplitude. Through a beam splitter, part of the light is received by PIN (FPD610-FS-VIS), and part of the light is reflected by MEMS-SM and irradiates the target surface. The diffuse reflection light is received by APD (APD430A). The analog signals generated by PIN and APD are synchronously collected by PXI-5164 and used to calculate the time of fight and then measure the distance. The signal of Lissajous scanning trajectory is generated by signal generator PXI-5413 and controller of MEMS-SM driving module. The angle measurement is completed by PXI-6221. So far, the angle and distance measurement synchronization signal are generated by Field Programmable Gate Array (FPGA) in PXI-5164

In the experiment, the targets are the Second Gate model of Tsinghua University, which is made of light-colored wood with a reflectivity of approximately 0.8 and a size of 300 × 65 × 293 mm, and the J-20 fighter model, which is made of dark camouflage engineering plastic with baking paint printing. Part of the surface is similar to specular reflection, with a reflectivity of approximately 0.2 and a size of 210 × 14 × 13 mm, including transparent plastic of the cockpit, with a size of approximately 10 × 5 mm. For example, the parameters of the Lissajous scanning trajectory are:\(n_x:n_y = 53:51\) , the FOV of the MEMS-SM is 10 × 10°, and after amplification by the transmitting antenna, the optical field of view FOV = 30 × 30°, \(f_0 = 10{\kern 1pt} {\kern 1pt} \rm{Hz}\) . The results are shown in Fig. 7.

The scanning trajectories in ref. 28 and designed in this paper are used to scan and image the target models respectively. The photos of the different trajectories and targets are shown. The targets are the models of the Tsinghua University Second Gate and the J-20 fighter. Three sampling schemes are used for imaging. One sampling scheme in ref. 28 is based on the scan trajectory designed in ref. 28. The other two sampling schemes are based on the scanning trajectory designed in this paper. a Targets and two scanning trajectories. b Point cloud obtained by different sampling schemes. c 3D front view images obtained by different sampling schemes

From the point cloud data, the outlines of the Second Gate model and J-20 model can be clearly distinguished. The scanning frequencies are \(f_x = 530{\kern 1pt} {\kern 1pt} {\rm {Hz}}\) and \(f_y = 510{\kern 1pt} {\kern 1pt} {\rm {Hz}}\) . The number of points per frame obtained by the sampling scheme in ref. 28 is 5406, and the angular resolution is \({\rm {AR}}_x = 0.5563^\circ\) and \({\rm {AR}}_y = 0.6119^\circ\) . The number of point clouds per frame obtained by the sampling scheme by knots is 5510, and the angular resolution is \({\rm {AR}}_x = 0.2781^\circ\) and \({\rm {AR}}_y = 0.3060^\circ\) . The number of points per frame obtained by the sampling scheme proposed in this paper is 10,812, and the angular resolution is \({\rm {AR}}_x = 0.2781^\circ\) and \({\rm {AR}}_y = 0.3060^\circ\) .

At the same hardware configuration, achieving higher angular resolution and more points per frame is the main work of this paper. Therefore, this paper first designs a better sampling scheme in the specific Lissajous trajectory and then summarizes the general sampling scheme to solve the coupling problem of Lissajous trajectory, sampling scheme and gridding. In the same MEMS-SM and the same scanning frequencies, the angular resolution and the number of points per frame are increased by 2 times (simulation in Table 1, experiment in Fig. 7b), and the effect of 3D reconstruction is better (Fig. 7c). A MEMS LiDAR prototype demonstration system is built to test the effectiveness of the improved method. However, we find that the sampling frequency we set in Fig. 3a is\(2f_0n_xn_y = 198{\kern 1pt} {\kern 1pt} {\rm {Hz}}\) , but only 99 sampling points are found in Fig. 3a. This is because the other 99 sampling points completely coincide with the 99 sampling points; that is, the image frame frequency is 2 Hz. Even so, the conclusion of Table 1 is still true.

For the imaging of a Lissajous trajectory system, such as with MEMS LiDAR, the sampling scheme of the point cloud not only is related to the trajectory but also needs to match the Cartesian grid. Therefore, it is difficult to utilize a universal sampling scheme. This paper establishes a functional relationship between the time and space parameters of point clouds so that the points sampled in equal time coincide with the Cartesian grid. Compared with the existing methods, the sampling scheme in this paper is more general, and the theoretical resolution performance of point cloud images is better.

Qi, P., Zhao, X. & Palacios, R. Autonomous landing control of highly flexible aircraft based on lidar preview in the presence of wind turbulence. IEEE Trans. Aerosp. Electron. Syst. 55, 2543–2555 (2019).

Lim, H.-C. et al. Evaluation of a Geiger-mode imaging flash lidar in the approach phase for autonomous safe landing on the Moon. Adv. Space Res. 63, 1122–1132 (2019).

Xie, S., Yang, D., Jiang, K. & Zhong, Y. Pixels and 3-D points alignment method for the fusion of camera and LiDAR data. IEEE Trans. Instrum. Meas. 68, 3661–3676 (2019).

Zhao, X., Sun, P., Xu, Z., Min, H. & Yu, H. Fusion of 3D LIDAR and camera data for object detection in autonomous vehicle applications. IEEE Sens. J. 20, 4901–4913 (2020).

Börcs, A. & Benedek, C. Extraction of vehicle groups in airborne lidar point clouds with two-level point processes. IEEE Trans. Geosci. Remote Sens. 53, 1475–1489 (2015).

Stone, W. C., Juberts, M., Dagalakis, N., Stone, J. & Gorman, J. Performance analysis of next-generation LADAR for manufacturing, construction, and mobility. NIST IR 7117 https://nvlpubs.nist.gov/nistpubs/Legacy/IR/nistir7117.pdf, https://doi.org/10.6028/NIST.IR.7117 (2004).

Stann, B. L. et al. Low-cost compact ladar sensor for ground robots. In Laser Radar Technology and Applications XIV. Vol. 7323, 73230X (Proceedings of the SPIE - The International Society for Optical Engineering, 2009).

Stann, B. L. et al. MEMS-scanned ladar sensor for small ground robots. In Laser Radar Technology and Applications XV (eds. Turner, M. D. & Kamerman, G. W.) Vol. 7684, 76841E (Proceedings of the SPIE-Conference on Laser Radar Technology and Applications XV, 2010).

Stann, B. L. et al. Brassboard development of a MEMS-scanned ladar sensor for small ground robots. In Laser Radar Technology and Applications XVI (eds. Turner, M. D. & Kamerman, G. W.) Vol. 037, 80371G (Proceedings of SPIE-Conference on Laser Radar Technology and Applications XVI, 2011).

Moss, R. et al. Low-cost compact MEMS scanning ladar system for robotic applications. In Laser Radar Technology and Applications XVII (eds. Turner, M. D. & Kamerman, G. W.) Vol. 8379, 837903 (Proceedings of SPIE-Conference on Laser Radar Technology and Applications XVII, 2012).

Stann, B. L. et al. Integration and demonstration of MEMS-scanned LADAR for robotic navigation. In Unmanned Systems Technology XVI (eds. Karlsen, R. E., Gage, D. W., Shoemaker, C. M. & Gerhart, G. R.) Vol. 9084, 90840J (Proceedings of SPIE-Conference on Unmanned Systems Technology XVI, 2014).

Stann, B. L., Dammann, J. F. & Giza, M. M. Progress on MEMS-scanned ladar. In Laser Radar Technology and Applications XXI (eds. Turner, M. D. & Kamerman, G. W.) Vol. 9832, 98320L (Proceedings of SPIE-Conference on Laser Radar Technology and Applications XXI, 2016).

Stann, B. L., Dammann, J. F., Giza, M. M. & Ruff, W. C. MEMS-scanned ladar for small unmanned air vehicles. In Laser Radar Technology and Applications, Vol. XXIII (eds Turner, M. D. & Kamerman, G. W.) 12 (SPIE, 2018).

Stann, B. L., Dammann, J. F., Giza, M. M., Redman, B. C. & Ruff, W. C. RF Coherent Detection on Top of Direct Detection Lidar. Encyclopedia of Modern Optics (Second Edition), Vol. 5, 2018, pp. 1–13.

Liu, J. J., Stann, B. L., Klett, K. K., Cho, P. S. & Pellegrino, P. M. Mid and long-wave infrared free-space optical communication. In Laser Communication and Propagation through the Atmosphere and Oceans VIII (eds. van Eijk, A. M., Hammel, S. & Bos, J. P.) Vol. 11133, UNSP1113302 (Proceedings of SPIE-Conference on Laser Communication and Propagation through the Atmosphere and Oceans VIII, 2019).

Kim, G., Eom, J., Park, S. & Park, Y. Occurrence and characteristics of mutual interference between LIDAR scanners. In Photon Counting Applications 2015 (eds. Prochazka, I., Sobolewski, R. & James, R. B.) Vol. 9504, 95040K (SPIE Conference on Photon Counting Applications, 2015).

Kim, G., Eom, J. & Park, Y. A hybrid 3D LIDAR imager based on pixel-by-pixel scanning and DS-OCDMA. In Smart Photonic and Optoelectronic Integrated Circuits XVIII (eds. He, S., Lee, E.-H. & Eldada, L. A.) Vol. 9751, 975119 (Proceedings of SPIE-Conference on Smart Photonic and Optoelectronic Integrated Circuits XVIII, 2016).

Hu, Q., Pedersen, C. & Rodrigo, P. J. Eye-safe diode laser Doppler lidar with a MEMS beam-scanner. Opt. Express 24, 1934 (2016).

Errando-Herranz, C., Le Thomas, N. & Gylfason, K. B. Low-power optical beam steering by microelectromechanical waveguide gratings. Opt. Lett. 44, 855 (2019).

Wu, J., Xu, H. & Zhao, J. Automatic lane identification using the roadside LiDAR sensors. IEEE Intell. Transp.Syst. Mag. 12, 25–34 (2020).

Weng, X. & Kitani, K. Monocular 3D Object Detection with Pseudo-LiDAR Point Cloud. In IEEE International Conference on Computer Vision Workshops (eds. 2019 IEEE/CVF ICCVW) 857–866 (IEEE/CVF International Conference on Computer Vision, 2019).

Qian, R. et al. End-to-End Pseudo-LiDAR for Image-Based 3D Object Detection. 2020 IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR). p. 5880–9, (2020).

Shan, T., Wang, J., Chen, F., Szenher, P. & Englot, B. Simulation-based lidar super-resolution for ground vehicles. Robotics and Autonomous Systems. Vol. 134, 103647 (2020).

Kasturi, A., Milanovic, V., Atwood, B. H. & Yang, J. UAV-borne lidar with MEMS mirror-based scanning capability. In SPIE Defense+Security (2016).

Yong, Y. K., Bazaei, A. & Moheimani, S. O. R. Video-rate Lissajous-scan atomic force microscopy. IEEE Trans. Nanotechnol. 13, 85–93 (2014).

Hwang, K., Seo, Y.-H. & Jeong, K.-H. High resolution and high frame rate Lissajous scanning using MEMS fiber scanner. In 2016 International Conference on Optical MEMS and Nanophotonics (OMN) (IEEE, 2016).

Brand, O. & Pourkamali, S. Electrothermal excitation of resonant MEMS. In Resonant MEMS 173–201 (John Wiley & Sons, Ltd, 2015).

Moriguchi, H., Wendt, M. & Duerk, J. L. Applying the uniform resampling (URS) algorithm to a Lissajous trajectory: fast image reconstruction with optimal gridding. Magn. Reson. Med. 44, 766–781 (2000).

Wu, J.-W., Lin, Y.-T., Lo, Y.-T., Liu, W.-C. & Fu, L.-C. Lissajous hierarchical local scanning to increase the speed of atomic force microscopy IEEE Trans. Nanotechnol. 14, 810–819 (2015).

Loewke, N. O. et al. Software-based phase control, video-rate imaging, and real-time mosaicing with a Lissajous-scanned confocal microscope. IEEE Trans. Med. Imaging 39, 1127–1137 (2020).

Karacor, D., Nazlibilek, S., Sazli, M. H. & Akarsu, E. S. Discrete Lissajous figures and applications. IEEE Trans. Instrum. Meas. 63, 2963–2972 (2014).

Leoncini, M. et al. Fully integrated, 406 μA, 5 °/hr, full digital output Lissajous frequency-modulated gyroscope. IEEE Trans. Ind. Electron. 66, 7386–7396 (2019).

D’Angelo, G., Laracca, M., Rampone, S. & Betta, G. Fast eddy current testing defect classification using Lissajous figures. IEEE Trans. Instrum. Meas. 67, 821–830 (2018).

Urey, H. et al. MEMS scanners and emerging 3D and interactive Augmented Reality display applications. In 2013 Transducers & Eurosensors XXVII: The 17th International Conference on Solid-State Sensors, Actuators and Microsystems (TRANSDUCERS & EUROSENSORS XXVII) 2485–2488 (2013).

Feng, H., Gu, H., Silbersweig, D., Stern, E. & Yang, Y. Single-shot MR imaging using trapezoidal-gradient-based Lissajous trajectories. IEEE Trans. Med. Imaging 22, 925–932 (2003).

Bauweraerts, P. & Meyers, J. Reconstruction of turbulent flow fields from lidar measurements using large-eddy simulation. Journal of Fluid Mechanics 906, (2021).

Hegna, T., Pettersson, H., Laundal, K. M. & Grujic, K. 3D laser scanner system based on a galvanometer scan head for high temperature applications. In Optical Measurement Systems for Industrial Inspection VII (eds. Lehmann, P. H., Osten, W. & Gastinger, K.) Vol. 8082, 80823Z (Proceedings of SPIE-Conference on Optical Measurement Systems for Industrial Inspection VII, 2011).

Kim, J.-D., Jung, J.-K., Jeon, B.-C. & Cho, C.-D. Wide band laser heat treatment using pyramid polygon mirror. Opt. Lasers Eng. 35, 285–297 (2001).

Stevenson, G., Verdun, H. R., Stern, P. H. & Koechner, W. Testing the helicopter obstacle avoidance system. In Applied Laser Radar Technology II (ed. Kamerman, G. W.) Vol. 2472 (Proceedings of SPIE-Symposium on OE/Aerospace Sensing and Dual Use Photonics, 1995).

Wang, J., Zhang, G. & You, Z. Design rules for dense and rapid Lissajous scanning. Microsyst. Nanoeng. 6, 1–7 (2020).

Jones, V. & Przytycki, J. Lissajous knots and billiard knots. Banach Cent. Publ. 42, 145–163 (1998).

Hoste, J. & Zirbel, L. Lissajous knots and knots with Lissajous projections. Preprint at arXiv:math/0605632v1, 87–106 (2006).

Likes, R. S. Moving gradient zeugmatography. No. US 4307343 (1981).

Hardy, C. J., Cline, H. E. & Bottomley, P. A. Correcting for nonuniform k-space sampling in two-dimensional NMR selective excitation. J. Magn. Reson. (1969) 87, 639–645 (1990).

Song, Z. L., Li, S. & George, T. F. Remote sensing image registration approach based on a retrofitted SIFT algorithm and Lissajous-curve trajectories. Opt. Express 18, 513–522 (2010).

You, Y. et al. Pseudo-LiDAR++: Accurate Depth for 3D Object Detection in Autonomous Driving. In ICLR (2020).

Koh, K. H., Kobayashi, T. & Lee, C. Investigation of piezoelectric driven MEMS mirrors based on single and double S-shaped PZT actuator for 2-D scanning applications. Sens. Actuators A: Phys. 184, 149–159 (2012).

Li, X. et al. Miniature laser projection display technique based on Lissajous scanning. Acta Opt. Sin. https://doi.org/10.3788/AOS201434.0612005 (2014).

Pearson, G. N., Ridley, K. D. & Willetts, D. V. Chirp-pulse-compression three-dimensional lidar imager with fiber optics. Appl. Opt. 44, 257 (2005).

This work is supported by the National Key Research and Development Program of China under grant No. 2016YFB0500902 and the Beijing Innovation Center for Future Chips, Tsinghua University.

These authors contributed equally: Junya Wang, Gaofei Zhang

School of Mechanical Science and Engineering, Huazhong University of Science and Technology, 430074, Wuhan, China

Junya Wang & Zheng You

Department of Precision Instrument, Tsinghua University, 10084, Beijing, China

Gaofei Zhang & Zheng You

State Key Laboratory of Precision Measurement Technology and Instruments, Tsinghua University, 10084, Beijing, China

Gaofei Zhang & Zheng You

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

The project was organized and coordinated by Z.Y.

The authors declare no competing interests.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Wang, J., Zhang, G. & You, Z. Improved sampling scheme for LiDAR in Lissajous scanning mode. Microsyst Nanoeng 8, 64 (2022). https://doi.org/10.1038/s41378-022-00397-9

DOI: https://doi.org/10.1038/s41378-022-00397-9

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Microsystems & Nanoengineering (Microsyst Nanoeng) ISSN 2055-7434 (online)