IAES International Journal of Robotics and Automation (IJRA) Vol. No. December 2025, pp. ISSN: 2722-2586. DOI: 10. 11591/ijra. ADC-LIO: A direct LiDAR-inertial odometry method based on adaptive distortion covariance Lixiao Yang. Youbing Feng Ocean College. Jiangsu University of Science and Technology. Zhenjiang. China Article Info ABSTRACT Article history: Focusing on the localization challenges for robots in dynamic navigation environments, this study proposes a direct LiDAR-inertial odometry (LIO) system named ADC-LIO, which achieves robust pose estimation and accurate map reconstruction using adaptive distortion covariance. ADC-LIO is engineered to address uncertain motion patterns in autonomous mobile robots, effectively integrating LiDAR scan undistortion within the Kalman filtering update process by embedding an iterative smoothing process and a backpropagation strategy. The ADC-LIO architecture enhances point cloud accuracy, improving the system's overall performance and robustness. addition, an adaptive covariance processing method is developed to resolve motion-induced sensing uncertainties, which calculates different covariances according to the error characteristics of the point cloud. This method enhances the constraints of high-quality point clouds, reduces the limitations on lowquality point clouds, and utilizes information more effectively. Experiments on the publicly available NTU-VIRAL dataset validate the effectiveness of ADCLIO, which improves pose estimation accuracy and reduces absolute position errors compared to other state-of-the-art methods, including FAST-LIO. Faster-LIO. FR-LIO, and Point-LIO. The proposed ADC-LIO is an appealing odometry method that delivers accurate, real-time, and reliable tracking and map-building results, posing a practical solution for robotic applications in structured indoor and GPS-denied outdoor environments. Received Jan 15, 2025 Revised Nov 7, 2025 Accepted Nov 11, 2025 Keywords: Dense mapping Kalman smoother Lidar inertial odometry Sensor fusion Simultaneous localization and This is an open access article under the CC BY-SA license. Corresponding Author: Youbing Feng Ocean College. Jiangsu University of Science and Technology Zhenjiang. China Email: yzfyb@just. INTRODUCTION Simultaneous localization and mapping (SLAM) technology is advancing rapidly in today's rapidly evolving robotic sector. As a cornerstone of autonomous robot navigation. SLAM directly enables critical robotic functionalities, including environment perception, path planning, and obstacle avoidance. Vision-based SLAM affords high location accuracy and is commonly employed in mobile robots owing to its lightweight and low cost . Ae. Nevertheless, algorithms that explore data from the visual domain, while providing high-resolution RGB data, do not capture depth directly and require substantial computational effort to recover the 3D scene necessary for motion planning. Opposing vision-based SLAM, light detection, and ranging (LiDAR) sensors afford enhanced LiDAR-odometry (LO) capabilities . and thus have become an integral part of self-driving cars . , . , unmanned aerial vehicles (UAV. , . , and other mobile Their high accuracy and ability to directly acquire 3D environmental information make them indispensable for robotic systems operating in GPS-denied or geometrically complex workspaces. The effectiveness of the LO system primarily depends on its precision, stability, and operational efficiency. Journal homepage: http://ijra. ISSN: 2722-2586 However, processing LiDAR data requires a large amount of computational resources, especially in high-speed motion or complex environments, which puts higher demands on real-time performance and To address this challenge, architectures that fuse LiDAR and inertial measurement units (IMU. , namely LiDAR-inertial odometry (LIO), have emerged . , . LIO provides a critical advantage for autonomous robots navigating in dynamic environments by synergizing LiDAR's precise geometry sensing with IMU's high-frequency motion tracking. This fusion strategy enhances the systemAos accuracy and reliability while reducing the computational burden through the use of intelligent algorithms, thereby enabling real-time processing. Various methods have been proposed to enhance registration accuracy, such as employing more precise probabilistic map representations . , applying finer LiDAR point segmentation . , . or selecting more informative points for registration . Additionally, several methods enhance accuracy by compensating for motion distortion in LiDAR data . , . FAST-LIO . and FAST-LIO2 . , as representative LIO methods, achieve highly accurate state estimation while maintaining computational efficiency through tightly coupled iterative Kalman filters. FAST-LIO2 further optimizes the data processing flow by introducing the incremental k-dimensional tree . kdTre. data structure for dynamic map updating and efficient nearest neighbor search, significantly improving the utility and flexibility of an LIO system. Faster-LIO . extends FAST-LIO2 by replacing the ikd-Tree with iVox and optimizing the code logic to achieve faster LIO processing. In . , the authors introduce a point-by-point . oint-wis. LIO framework that fuses LiDAR points during sampling without accumulating them in the frame. Removing point accumulation eliminates intra-frame motion distortion, allowing for high odometry output and build-up updates at near-point sampling rates. This enables the system to track swift motion and simultaneously support incremental updating and dynamic rebalancing, allowing odometry and map-building frequencies of up to 100 Hz despite limited computational resources. The work in . proposes a tightly coupled, iterative error state Kalman smoothing (ESKS) framework for state estimation, which achieves robust tracking and location in violent motion scenarios by smoothing the position in the degraded direction of the current subframe. This smoothing process is completed by combining the history and constraint information on future subframes. Nevertheless, current LIO systems still face several challenges. For instance, the IMU attempts to correct point cloud data distortion by integrating positional information. However, these errors are challenging to eliminate due to sensor noise, drift, and other dynamic influences during the integration process of the IMU. In robotic applications where sustained accuracy over long missions is crucial, motion distortions often persist in the point cloud, limiting further improvements in the LIO system's accuracy. To address this issue, this paper introduces an enhanced method named iterated backward propagation, which combines iterative Kalman smoothing with backpropagation to reduce distortion errors in the LiDAR This innovation focuses on the motion distortion patterns observed in wheeled robots during abrupt acceleration/deceleration and UAVs during aerial rotations. The approach iteratively refines the LiDAR scan's undistortion process, thereby improving the accuracy and robustness of the system and leading to more precise pose estimation and overall enhanced performance. Moreover, existing LIO systems typically assume uniform confidence levels for all points in the LiDAR scan when formulating the observation equations. However, due to the motion of the LiDAR sensor, different points in the scan exhibit varying residual distortion errors. Therefore, an adaptive distortion covariance calculation method is developed to weigh the points with different distortion uncertainty during the Kalman update process. This method enhances the weighting of high-quality points, reduces the influence of low-quality points, and ensures more efficient utilization of the available information, ultimately improving the system's overall performance and robustness. The proposed method's effectiveness is verified by qualitative and quantitative experiments on the publicly available NTU-VIRAL dataset . The experimental results demonstrate that, compared with FAST-LIO2, the proposed algorithm improves average accuracy and exhibits higher refinement in the build effect. The remainder of this paper is organized as follows. Section 2 introduces and analyzes the proposed Section 3 presents the experimental findings on the NTU-VIRAL dataset, and section 4 concludes this work. METHODOLOGY Framework overview The proposed ADC-LIO system is illustrated in Figure 1. It is designed to process high-frequency . IMU data and relatively low-frequency . LiDAR scans, outputting odometry information at the frequency of the LiDAR scan. ADC-LIO relies on an iterative error Kalman filter framework, where IMU data is used for Kalman prediction and point cloud data is used for the Kalman update process to maintain pose estimation consistency under dynamic robotic motion. During the iterative update procedure. Iterated Backward Propagation further enhances the point cloud undistortion effect, which is critical for robots operating in high-acceleration scenarios. Additionally, an adaptive distortion covariance calculation balances IAES Int J Rob & Autom. Vol. No. December 2025: 551-560 IAES Int J Rob & Autom ISSN: 2722-2586 the constraints of high-precision points and reduces the influence of erroneous points, thereby improving the overall performance and accuracy of the system. LiDAR . Motion Compensation State Update Local Map . Kd-Tre. IMUs . Forward Propagation Iterated Backward Propagation Odometry Figure 1. System overview of ADC-LIO Forward propagation Let the state variables of the system be [ ya ycIya , ya ycyycn , ya ycya , ycayca , ycayui , ya yc. , where ya ycyycn , ya ycIya , and ya ycya denote the IMU position, rotation, and speed in the global reference frame, ya yci is the gravity vector in the global reference frame, ycayca and ycayui are the IMU biases, modeled as a random walk process. The forward transfer of the IMU is computed based on the discrete IMU motion model. The first IMU frame is defined as the global reference frame in this process. The kinematic model for forward propagation is computed using the kinematic model below, which propagates the IMU's state from the initial frame through subsequent time steps. ya ycIya = ya ycIya . uiyco Oe ycayui Oe ycuyui ]O , ya ycyya = ya ycya ya. yc = ya ycIya . ca Oe ycayca Oe ycuyca ) ya yci ya yco ycayui = ycuycayui , yca. yca = ycuycayca ya . yci = 0, yaycIya = 0, yaycyya = 0 where ycayco and yuiyco are IMU measurements, ycuyca and ycuyui represents the noise measurement of ycayco and yuiyco , respectively, and . ^ is the skew-symmetric cross product matrix of vector yca OO ycI3 . Iterative backward propagation Unlike methods such as FAST-LIO2, which compensate for distortion with a single forward and backward propagation, this paper introduces an iterative backward propagation scheme integrated with Kalman smoothing. This approach effectively combines point cloud distortion removal with state updating, enhancing both processes. Specifically, smoothing the IMU position after each state update provides more accurate LiDAR coordinates for the next backward propagation iteration, further reducing residual distortions. Figure 2 illustrates the iterative process of the proposed scheme, which progressively refines the system's performance through backward propagation, state updating, and Kalman smoothing. These modules are analyzed below. Figure 2. Iterative backward propagation, state updating, and Kalman smoothing. Smoothing the IMU position after each state update enhances the accuracy of the LiDAR coordinates for the next backward propagation iteration, thereby reducing residual distortions ADC-LIO: A direct LiDAR-inertial odometry method based on adaptive distortion covariance (Lixiao Yan. A ISSN: 2722-2586 Motion compensation The LiDAR scan comprises multiple points sampled at different instants, often at higher frequencies than the IMU. This results in motion distortion of the points within each LiDAR frame due to varying sampling times. To overcome this concern, this study introduces an interpolation-based approach, where each LiDAR point is mapped to its corresponding IMU position at the sampling moment. The interpolated points are then remapped into a standard coordinate system at a unified time, thereby mitigating distortion within the LiDAR scan frame. Specifically, the interpolation method utilizes IMU trajectory data to adjust the position of each LiDAR point, employing both spatial and temporal interpolation to align the points to a standardized This process effectively reduces motion aberrations caused by the temporal discrepancies between the LiDAR and IMU data. State update After the LiDAR scan data has undergone undistortion, the undistorted point cloud is used to formulate the observation residuals for system state optimization. Further details on deriving the observation residual function are outlined in section 2. 4, and the system state update, which relies on . , is presented in section 2. By incorporating the constraints from the undistorted point cloud data, this process adjusts the position and orientation estimates, thereby improving the accuracy of the state estimation. A key aspect of this procedure is transforming the observation information from the LiDAR scan into constraints on the position and attitude, thereby effectively integrating data from multiple sensors. This integration enhances the optimization process by leveraging complementary information and refining the system's performance. Backward Kalman smoothing After the state update, the IMU position estimates obtained during the forward propagation are refined using Kalman smoothing. This process corrects the current position estimates and utilizes past and future observation data, providing a more accurate and smooth result. The corrected IMU positions, refined through Kalman smoothing, are then used in the subsequent backward propagation step to enhance the point cloud undistortion further. The advantage of Kalman smoothing over traditional filtering is that it optimizes the system state by considering a broader range of observations, including future data, thereby improving the accuracy of the current state and enhancing overall system stability. This mutual enhancementAiwhere the improved estimates from the smoothing phase enhance the following undistortion process, which improves the next updateAicreates a positive feedback loop, leading to progressively more accurate results . = yce. cuyco , ycyco ) . ycEyco 1 = . a yayco yuuy. ycEyco . a yayco yuuy. ycN . ayco yuuy. ycEyco . ayco yuuy. ycN . Oe ]Oe1 yayco = ycEyco . a yayco yuuy. ycN . cEyco 1 Oe ] ycuycoyc = ycuyco yayco . Oe ycuyco 1 Oe ]ya ycN ycEycoyc = ycEyco yayco . cEyco 1 Oe ycEyco 1 For the detailed forms of yce. cuyco , ycyco ), yayco , yayco and ycEyco , the readers can refer to . , ycuycoyc and ycEycoyc represent the outcomes of the backward smoothing process, which is initiated from the final sub-frame at time ycN , with ycuycNyc = ycuycN and ycEycNyc = ycEycN . A comprehensive description of the Kalman smoother is provided in . When applying the Kalman smoothing process to the error state, certain adjustments to the backward smoothing stage are required, as presented in . Adaptive covariance-based observation models Measurement model: After the backward propagation, all LiDAR points are considered to be sampled simultaneously at the end of the scan. Let yco be the LiDAR scan index and { yaycEyc , yc = 1,UIUIUI, yc. the points in the yco -th scan, expressed in the local LiDAR coordinate frame ya at the scan end-time. Each measured point yaycEyc is affected by measurement noise yaycuyc comprising the ranging errors and beam-direction yciyc Eliminating this noise yields the true point position in the local LiDAR coordinate frame yaycEyc . ya yciyc ycEyc = yaycEyc yaycuyc IAES Int J Rob & Autom. Vol. No. December 2025: 551-560 IAES Int J Rob & Autom ISSN: 2722-2586 Once projected to the global coordinate frame using the respective LiDAR pose ya ycNyayco = ( ya ycIyayco , ya ycEyayco ) and extrinsic yaycNya , this true point should coincide with a local planar surface in the map, i. 0 = ya ycycycN ( ya ycNyaya yaycNya ( yaycyyc yaycuyc ) Oe ya ycyc ) . where ya ycyc denotes the plane's normal vector and ya ycyc represents a point located on that plane in Figure 3. Both ya ycNyayco and yaycNyayco are components of the state vector ycUyco . Hence, the measurement for the yc -th point can be simplified from . into the compact expression: 0 = h j ( xk . L p j L n j ) . which implicitly defines the measurement model for the state vector ycUyco . LiDAR IMUs Figure 3. Observation models Unlike FAST-LIO2, which assumes that all points in the point cloud have the same confidence level, this paper considers that intense sensor motion and significant de-distortion time differences often degrade point cloud accuracy. Therefore, this paper proposes an adaptive error covariance computation scheme based on the intensity of IMU motion and the de-alignment time difference, which dynamically adjusts the error weights in the point cloud. The proposed strategy enhances the overall accuracy of the LIO system, as detailed in the process presented below. First, the motion intensity of the sensor is characterized by the acceleration and angular velocity of the IMU. Let the acceleration paradigm of the IMU at moment ycycn be . cycn )|, and the angular velocity paradigm be . cycn )|. Then, the motion intensity can be expressed as: M i = A a ( ti ) A A ( ti ) . where yu and yu are weighting coefficients regulating the relative influence of acceleration and angular velocity on the error. The more violent the sensor's motion, the greater the motion intensity ycAycn and the corresponding point cloud error. Second, consider the effect of the de-distortion time difference. For point ycyycn , whose acquisition moment is ycycn and the current moment is ycycaycycyc , the time difference yuuyc = ycycaycycyc Oe ycycn reflects the lag degree of the de-distortion process. As the time difference increases, the error in the de-distortion process gradually accumulates, expressed as: Ti = AAEt ADC-LIO: A direct LiDAR-inertial odometry method based on adaptive distortion covariance (Lixiao Yan. A ISSN: 2722-2586 where yu is the time lag error coefficient, which measures the contribution of the time lag to the error. Based on the above two factors. Mi and Ti, the covariance matrix Oc ycn for each point Pi in the point cloud can be expressed as: Oc ycn = yua02 . cAycn ycNycn )ya . where yua02 denotes the underlying noise variance and ya is the unit matrix, which indicates that the error characteristics are consistent in all directions. This adaptive covariance model enables the LIO system to dynamically adjust the confidence level of different points based on the error characteristics of the point cloud Specifically, point clouds with more significant errors . , points with high motion intensity or considerable time differenc. will receive lower confidence levels, thus weakening their impact on the system state estimation. In contrast, high-quality point clouds will increase their contribution to the system state estimation by enhancing the constraints. This adaptive covariance method effectively enhances the robustness and accuracy of the LIO system in complex dynamic environments, particularly in cases where sensor motion is violent, and the quality of point cloud data is uneven, allowing for the more effective utilization of valuable information and mitigating the negative impact of low-quality data on system performance. Kalman update The point-to-plane residuals formulated in . are employed as measurement inputs for the iterative Kalman state update: ya = . a ycN ycIOe1 ya ycEycoOe. Oe1 yaycN ycIOe1 ycuycoya 1 = ycuycoya ya . Oe Ea. cuycoya )) Oe . a Oe yay. cuycoya Oe ycuyco0 ) . where ya denotes the Jacobian matrix of the measurement model, reflecting the sensitivity of the residuals to state variations, and ycI represents the measurement noise covariance, which models the uncertainty in the observed measurements. The updated nominal state ycuycoya 1 is iteratively updated by correcting the predicted state ycuyco0 using the Kalman gain ya and the residual between the measurement model Ea. cuycoya ) and the actual During each iteration, the nominal state was corrected, and the error state was reset. When the iterations converged, the state covariance was updated using . ycEyco = . a Oe yay. ycEyco . IMU data were then used for nominal state prediction via a discrete-time propagation model: ycuyco 1 = yce. cuyco , ycyco ) . The linearized discrete-time error-state transition in Eq. was applied to propagate the error-state covariance forward, as expressed in . yuycuyco 1 = . a yayco yuuy. ayco yuuy. ycEyco 1 = . a yaya yuuy. ycEyco . a yayco yuuy. ycN . ayco yuuy. ycEyco . ayco yuuy. ycN . The outcome of this forward propagation was then applied for motion compensation and registration. RESULTS AND DISCUSSION To validate the performance of our proposed method, it is compared with Fast-LIO2 . Faster-LIO . Point-LIO . , and FR-LIO . on nine sequences of the publicly available NTU-VIRAL dataset . Figures 4 to 5 depict the mapping and simulation results of ADC-LIO in different dataset sequences. NTU-VIRAL data integrates diverse sensor configuration, including two 3D LIDARs, dual time-synchronized cameras, multiple IMUs, and four sets of ultra-wideband (UWB) nodes, to form an all-encompassing sensing system for the complex requirements of autopiloted vehicles. Through IAES Int J Rob & Autom. Vol. No. December 2025: 551-560 IAES Int J Rob & Autom ISSN: 2722-2586 flight testing in various indoor and outdoor environments. NTU VIRAL has captured high-quality multimodal data. The technical highlight of this dataset is its unique and detailed sensor integration scheme, which mimics the sensor configuration of an autonomous vehicle but exhibits more complex challenges and characteristics in an airborne application. Table 1 presents the experimental findings of the developed ADC-LIO and current state-of-the-art (SOTA) LiDAR-inertial odometry systems on the nine sequences, i. , e_01 to _03, nya_01 to _03, and sbs_01 to _03, of the NTU-VIRAL dataset. The performance is evaluated based on the Absolute Position Errors (APE), indicating that ADC-LIO achieves the slightest absolute error across all sequences. This improvement is due to our system considering the impact of motion distortion on pose estimation through adaptive covariance, which enhances the influence of points with low distortion and high accuracy during the pose estimation process, thereby increasing the overall system precision. Overall. ADC-LIO achieves an APE of approximately 50%, ranking as the second-best performing method, which demonstrates that the proposed method demonstrates superior performance compared to existing SOTA approaches. Table 1. Absolute position errors in NTU-VIRAL dataset with different algorithms Algorithm Fast-LIO2 . Faster-LIO . Point-LIO . FR-LIO . ADC-LIO e_01 e_02 e_03 nya_01 Sequence nya_02 nya_03 sbs_01 sbs_02 sbs_03 Figure 4 visualizes the corresponding experimental results of ADC-LIO. The first row illustrates the estimated trajectories represented by blue curves alongside the dense 3D point clouds generated by ADC-LIO for the e_01, nya_01, and sbs_01 sequences. The point clouds are color-coded according to the intensity information of each point, highlighting our methodAos capability to preserve fine details. The high precision of the pose estimation results in minimal ghosting effects in the dense point clouds, thereby effectively capturing the intricate structures of the surrounding environment. Additionally, the second row depicts the estimated trajectories alongside the APE metric, where the gray dashed lines represent the ground truth. The trajectories are color-coded to reflect varying magnitudes of error, providing a clear visual representation of the algorithmAos accuracy relative to the true trajectory. e_01 nya_01 sbs_01 Groundtruth LIO Figure 4. ADC-LIO mapping results under NTU-VIRAL sequences Figure 5 illustrates the localization and mapping results of sequence e_01. Specifically. Figure 5. explicitly illustrates the proposed algorithm's localization and mapping results for the e_01 Figures 5. visualize the high-quality reconstruction results for non-structural, complex objects, such as trees, while Figures 5. demonstrate the reconstruction results for structured objects, such as buildings. These high-quality reconstructions are directly derived from the iterative backpropagation method proposed in this paper, significantly enhancing the quality of point cloud ADC-LIO: A direct LiDAR-inertial odometry method based on adaptive distortion covariance (Lixiao Yan. A ISSN: 2722-2586 The adaptive distortion covariance also enhances the systemAos localization accuracy, resulting in more accurate point cloud reconstruction with reduced ghosting artifacts. Figure 5. Localization and mapping results of sequence e_01: . overall 3D reconstruction, . vegetation mapping, . mixed vegetationAestructure region, . building fayade reconstruction, and . structured object details CONCLUSION This paper proposes ADC-LIO, an accurate and robust LiDAR inertial odometry framework that utilizes an iterative error Kalman framework and an iterative backpropagation strategy to effectively integrate the point cloud de-distortion processing with the Kalman filter update process. The proposed method directly addresses the motion distortion challenges of mobile robots during high-speed navigation and aggressive In addition. ADC-LIO incorporates an adaptive covariance processing method to further enhance the constraints of high-quality point clouds. ADC-LIO is evaluated on nine sequences from the publicly available NTU-VIRAL dataset, which involves both indoor environments with significant rotational speeds and outdoor environments. The tests demonstrated that ADC-LIO produces accurate and reliable map-building results, affording an improved pose estimation accuracy compared to other SOTA approaches, including FAST-LIO. Faster-LIO. Point-LIO, and FR-LIO. In forthcoming studies, we intend to conduct a more in-depth analysis of variable covariance, aiming to improve the performance of LIO in challenging robotic scenarios such as long-term autonomous exploration and multi-robot collaborative mapping. FUNDING INFORMATION Authors state no funding involved. AUTHOR CONTRIBUTIONS STATEMENT This journal uses the Contributor Roles Taxonomy (CRediT) to recognize individual author contributions, reduce authorship disputes, and facilitate collaboration. Name of Author Lixiao Yang Youbing Feng C : Conceptualization M : Methodology So : Software Va : Validation Fo : Formal analysis ue ue ue ue ue ue ue I : Investigation R : Resources D : Data Curation O : Writing - Original Draft E : Writing - Review & Editing IAES Int J Rob & Autom. Vol. No. December 2025: 551-560 Vi : Visualization Su : Supervision P : Project administration Fu : Funding acquisition IAES Int J Rob & Autom ISSN: 2722-2586 CONFLICT OF INTEREST STATEMENT Authors state no conflict of interest. DATA AVAILABILITY Derived data supporting the findings of this study are available from the corresponding author [Feng Y] on request. REFERENCES