Task- and Setup-Agnostic Robot Localization and State Estimation with Factor Graphs

IEEE Transactions on Robotics (T-RO), 2026
1Robotic Systems Lab, ETH Zürich  ·  2Max Planck Institute for Intelligent Systems  ·  3NASA Jet Propulsion Laboratory, California Institute of Technology
Holistic Fusion illustration: robots, measurements, and a factor graph
Holistic Fusion (HF) fuses IMU, absolute, landmark, and local measurements from an arbitrary number of sensors and reference frames in one factor graph.
Overview

Video

A short overview of the approach and the real-world deployments on the ANYmal quadruped, the RACER off-road vehicle, and the HEAP walking excavator.

Paper

Abstract

Seamless operation of mobile robots in challenging environments requires low-latency local motion estimation and accurate global localization. While most sensor-fusion approaches are designed for specific scenarios, this work introduces a flexible open-source solution for task- and setup-agnostic multimodal sensor fusion distinguished by its generality and usability. Holistic Fusion formulates sensor fusion as a combined estimation problem of i) the local and global robot state and ii) a (theoretically unlimited) number of dynamic variables, including automatic alignment of reference frames; this formulation fits countless real-world applications without conceptual modifications, offering a comprehensive solution beyond hard-coded/task-specific approaches. The proposed factor-graph formulation enables direct fusion of an arbitrary number of absolute, local, and landmark measurements expressed with respect to different frames by explicitly including them as states in the optimization and modeling their evolution as random walks. Moreover, local smoothness and consistency receive particular attention to prevent estimation jumps. Holistic Fusion enables low-latency and smooth online state estimation on typical robot hardware while simultaneously providing low-drift global localization at the IMU measurement rate. The efficacy of this released framework is demonstrated in five real-world scenarios on three robotic platforms with distinct task requirements, highlighting the advantages of fusing multiple absolute measurement types.

Approach

What is Holistic Fusion?

HF acts as the central state-estimation module of a robot. Beyond the robot state, it estimates the relationship between all reference frames of the fused measurements and external modules, and it feeds high-rate beliefs back to them.

Role and information streams of Holistic Fusion
Intended role of HF: a central fusion module that also estimates the (relative) transformations between the reference frames of measurements and external modules, e.g., a SLAM system, and propagates high- or low-rate beliefs back to them.

Objectives

  • O1 Real-time robot state at IMU rate: locally smooth in the odometry frame, globally accurate in the world frame.
  • O2 Flexible fusion of any number and type of delayed, out-of-order measurements without prior engineering.
  • O3 Automatic alignment of all reference frames, i.e., a synchronized localization manager.
  • O4 Online and offline extrinsic calibration between sensor frames.
  • O5 Easy use for new setups: adding a new measurement type should be simple.
Illustrative Holistic Fusion scenario on ANYmal
An illustrative scenario: IMU, GNSS, LiDAR-registration poses, and local velocity are fused directly without preprocessing. Global and non-global absolute measurements drift against each other; HF explicitly estimates the shift between their reference frames, aligns them at local keyframes, and models the evolution as a random walk. The graph on top shows the resulting factor graph with states and factors in matching colors.

Key Ideas

Compared to existing sensor-fusion formulations, the state vector of HF is not fixed. Dynamic context variables are allocated as needed, and part of the logic that is usually hard-coded is offloaded to the optimization.

1

Holistic state variables

Next to the IMU-rate robot navigation states, HF dynamically creates reference-frame alignment states, landmark states, and global (e.g., calibration) states and optimizes all of them jointly.

2

Drift modeled as a random walk

A single rigid (Umeyama) alignment cannot capture the drift between frames. HF models the evolution of each reference-frame transformation as a multivariate discrete-time random walk on SE(3), added as a zero-mean between factor.

3

Local keyframe alignment

Alignment variables are parametrized around local keyframes instead of the global origin. This removes the growing rotational lever arm and makes automatic alignment numerically stable over long missions.

4

Smooth, non-jumping odometry

New information (e.g., returning GNSS) may jump the world-frame belief. HF integrates the body-frame velocity from the smoother window to provide a smooth odometry-frame estimate for control and local mapping.

5

Online: high-rate, out-of-order

IMU-based propagation runs at full rate while a fixed-lag iSAM2 smoother optimizes asynchronously. Dense states allow delayed and out-of-order measurements to be inserted without rewiring the graph.

6

Offline batch and pseudo ground truth

The offline graph mirrors the online one and is initialized from the online belief, converging in few iterations. It yields IMU-rate optimized trajectories usable as post-mission pseudo ground truth.

Factor graph structure of Holistic Fusion
Structural overview of the HF factor graph. IMU measurements drive the creation of the navigation states; all other states are created dynamically based on the provided measurements: a) reference-frame alignment states, b) landmark states, c) global states, d) example measurement types depicted as factors. Every new measurement type follows the same four-step template (Steps A to D), so Jacobians are derived automatically via GTSAM expression factors.

Supported measurement types

Four categories cover most practical robotic state-estimation problems, including non-gravity-aligned and drifting measurements expressed in arbitrary sensor frames with respect to arbitrary reference frames.

  • IMU measurements (core sensor, defines the state-creation rate).
  • Absolute measurements w.r.t. a reference frame, e.g., GNSS positions or SLAM poses in a map frame.
  • Landmark measurements, e.g., foot contact points or feature locations.
  • Local and relative measurements, e.g., wheel or RADAR velocities and odometry deltas.
Smooth odometry generation by integrating the body-frame state
Motivation and illustration of the smooth odometry generation by integrating the body-frame state estimated in the current smoother window.
Open Source

The Framework

HF is released as a documented C++ framework with main dependencies on GTSAM and Eigen. The core is independent of any middleware; ROS 1 and ROS 2 wrappers and a set of robot examples are provided.

Software structure of Holistic Fusion
High-level overview of the software structure and examples. HF Core is fully generic and templated, HF Interface provides concrete implementations, and HF ROS and ROS 2 are the corresponding middleware wrappers.

Core library

holistic_fusion: middleware-independent C++ library. New measurement types are added by implementing a small interface; Jacobians come from GTSAM expression factors.

ROS 1 & ROS 2

holistic_fusion_ros and holistic_fusion_ros2 provide callbacks, logging, visualization, TF handling, and message advertisement.

Examples

Ready-to-run estimators for several platforms, evaluated in the paper or used in teaching and ground-truth generation.

  • ANYmal quadruped (IMU, GNSS, leg odometry, LiDAR poses)
  • HEAP walking excavator (IMU, dual GNSS, LiDAR poses)
  • SuperMegaBot (ROS 1 and ROS 2)
  • Leica total station and GNSS ground-truth generation
  • All examples

Get started with the installation guide, the system overview, and the parameter reference.

Evaluation

Experiments

HF is the default localization and state-estimation solution on all evaluated platforms. Its suitability is demonstrated in five real-world scenarios on three robots with different sensor suites, motions, and environments.

Evaluated robotic platforms
Evaluated platforms: three ANYmal datasets (autonomous hiking, parkour, and indoor locomotion with motion capture), the RACER off-road vehicle, and the HEAP walking excavator. The SuperMegaBot and the ground-truth generation setup are available as open-source examples.
PlatformSensors (all: IMU, LiDAR)MotionsEnvironments
ANYmal
quadrupedal robot
Leg odometry, single GNSS antenna Dynamic, vertical (climbing), high-acceleration stomping, multi-contact, leg slip, long distances Indoors and outdoors, mixed
RACER
off-road vehicle
Three LiDARs, RADAR, GNSS, single wheel encoder Highly dynamic, wheel slip, unpaved ground, few geometric features Outdoors
HEAP
walking excavator
Two GNSS antennas Multi-hour operation, tracking and control in the world frame Outdoors, covered by building structures
Autonomous forest hike with ANYmal
Fully autonomous hiking experiment with ANYmal in a forest. A: online and offline trajectories with highlighted LiDAR drift and GNSS loss. B: the aligned LiDAR trajectory at three keyframes, locally aligning perfectly at each. C: evolution of the map frame w.r.t. the world frame, estimated as a random walk.
Elevation map in odometry frame vs. world frame
Local elevation map expressed in the odometry frame vs. the world frame on the forest hike. The return of GNSS causes an update jump in the world frame that corrupts the map, while the odometry estimate stays smooth and remains suitable for local mapping and navigation.

Highlights

  • Fusing multiple absolute measurement types (e.g., GNSS and LiDAR poses) improves over any single one.
  • Reference-frame drift is estimated well online and offline; without random-walk modeling, alignment errors grow over time.
  • Local keyframe alignment is essential: alignment around the global origin degrades far from it.
  • HF Odom is markedly smoother than the world estimate and than purely local estimators, including during slip.
  • Offline optimization of hour-long missions with hundreds of thousands of variables takes about a minute.
Data

Datasets

The recorded data used for the ANYmal and HEAP examples is publicly available. Combined with the provided replay launch files, the estimators can be run on the recordings out of the box.

ANYmal

Quadrupedal robot recordings with IMU, GNSS, leg odometry, and LiDAR, corresponding to the anymal_estimator_graph example.

HEAP

Walking-excavator recordings with IMU, dual GNSS, and LiDAR, corresponding to the excavator_holistic_graph example.

Download datasets (Google Drive)
Reference

Citation

If you find this work or the code useful, please consider citing the following publications.

Holistic Fusion (T-RO 2026)

@article{nubert2026holistic,
  title     = {Holistic Fusion: Task- and Setup-Agnostic Robot Localization and State Estimation with Factor Graphs},
  author    = {Nubert, Julian and Tuna, Turcan and Frey, Jonas and Cadena, Cesar and Kuchenbecker, Katherine J. and Khattak, Shehryar and Hutter, Marco},
  journal   = {IEEE Transactions on Robotics},
  year      = {2026},
  publisher = {IEEE},
  doi       = {10.1109/TRO.2026.3714645}
}

Graph-Based Multi-Sensor Fusion (ICRA 2022)

@inproceedings{nubert2022graph,
  title        = {Graph-based Multi-sensor Fusion for Consistent Localization of Autonomous Construction Robots},
  author       = {Nubert, Julian and Khattak, Shehryar and Hutter, Marco},
  booktitle    = {IEEE International Conference on Robotics and Automation (ICRA)},
  year         = {2022},
  organization = {IEEE}
}
People

Authors

Robotic Systems Lab, ETH Zürich · Max Planck Institute for Intelligent Systems · NASA Jet Propulsion Laboratory

Julian Nubert
ETH Zürich · MPI-IS · JPL
Turcan Tuna
ETH Zürich
Jonas Frey
ETH Zürich · JPL
Cesar Cadena
ETH Zürich
Katherine J. Kuchenbecker
MPI-IS
Shehryar Khattak
JPL
Marco Hutter
ETH Zürich

Corresponding author: Julian Nubert, nubertj@ethz.ch

Thanks

Acknowledgments

The authors thank their colleagues at ETH Zürich and NASA JPL for help conducting the robot experiments and evaluations and for using GMSF and HF on their robots. Special thanks go to Takahiro Miki and the ANYmal Hike team at RSL; Nikita Rudin and David Hoeller for the ANYmal Parkour experiments; Patrick Spieler for running the deployments on the JPL RACER vehicle; the entire HEAP team at RSL and Gravis Robotics; Thomas Mantel and the teaching assistants of the ETH Robotic Summer School for their help on the SuperMegaBots; and Mayank Mittal for help generating the renderings.

This work is supported in part by the Sony Research Grant 2023, the EU Horizon 2020 programme grant agreements No. 852044, 101016970, and 101070405, EU Horizon 2021 programme grant agreement No. 101070596, the ETH Zurich Research Grant No. 21-1 ETH-27, the Swiss National Science Foundation (SNSF) through project No. 227617, the NCCR Digital Fabrication, and the Max Planck ETH Center for Learning Systems. This research was partially conducted at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (80NM0018D0004). This work was partially supported by the Defense Advanced Research Projects Agency (DARPA).