---
title: Weighted Covariance Intersection (WCI)
url: https://www.emergentmind.com/topics/weighted-covariance-intersection-wci
type: topic
---

# Weighted Covariance Intersection (WCI)

Weighted Covariance Intersection (WCI) is the weighted form of covariance intersection, a fusion method for combining multiple Gaussian estimates when the cross-correlation between them is unknown or only partially known. In the standard formulation, for estimates \((x_i,P_i)\), WCI computes
\[
P_{\text{WCI}}^{-1}=\sum_{i=1}^{N}\omega_i P_i^{-1},\qquad
x_{\text{WCI}}=P_{\text{WCI}}\left(\sum_{i=1}^{N}\omega_i P_i^{-1}x_i\right),
\]
with \(\omega_i\ge 0\) and \(\sum_i \omega_i=1\). The method is designed so that the fused covariance remains a valid upper bound for all cross-correlations consistent with the available information, which makes WCI a standard conservative rule in distributed estimation and sensor fusion [2603.20402].

## 1. Classical formulation and consistency

In the CI literature, WCI usually denotes the case in which the fused covariance is explicitly parameterized by weights and those weights are chosen, often heuristically and sometimes optimally, according to a scalar uncertainty criterion. For two full-state estimates, the classical formula is
\[
P_{\text{CI}}^{-1}=\omega P_1^{-1}+(1-\omega)P_2^{-1},\qquad
x_{\text{CI}}=P_{\text{CI}}\big(\omega P_1^{-1}x_1+(1-\omega)P_2^{-1}x_2\big),
\]
with \(\omega\in[0,1]\). For \(N\) estimates, the weights lie on the simplex \(\Delta^N=\{\omega\in\mathbb{R}^N:1^\top\omega=1,\ \omega_i\ge 0\}\) [2603.20402].

The central property of CI and WCI is consistency under unknown cross-correlation. In the positive definite, or Loewner, order, if each local estimate is consistent, then the CI-fused estimate remains consistent without requiring any cross-covariance terms. In information form, this is commonly written as
\[
\Sigma_{\hat z'}^{-1}=\sum_{i=1}^N c_i \Sigma_i^{-1},\qquad
\Sigma_{\hat z'}^{-1}\mathbb{E}[\hat z']=\sum_{i=1}^N c_i \Sigma_i^{-1}\mathbb{E}[\hat z_i],
\]
with convex coefficients \(c_i\ge 0\), \(\sum_i c_i=1\). This information-space convex combination is the basic weighted mechanism underlying WCI [2108.08789].

The weight vector determines which conservative bound within the CI family is selected. Different \(\omega\) correspond to different worst-case fused covariances, and the usual design task is to choose \(\omega\) to minimize a scalar performance index \(J(P)\). Typical choices stated in the literature are \(J(P)=\operatorname{tr}(P)\), \(J(P)=\log\det(P)\), and \(J(P)=\lambda_{\max}(P)\), under the monotonicity assumption \(X\succ Y\Rightarrow J(X)>J(Y)\) [2603.20402].

## 2. Weight selection as an optimization problem

A standard interpretation of WCI is that it is optimal only within a specific parametric family of conservative bounds. In that sense, a “family-optimal” solution is optimal within the CI family rather than among all possible consistent fusion rules. Recent work formalizes this view through overlapping covariance intersection (OCI), which embeds standard CI and split covariance intersection (SCI) in a single semidefinite-programming framework [2603.20402].

For standard CI, the family-optimal weight selection problem can be written as the semidefinite program
\[
\begin{aligned}
&\underset{B\in S_+^n,\ \boldsymbol{\omega}\in\Delta^N}{\arg\min}\ J(B)\\
&\text{s.t.}\quad
\begin{bmatrix}
B & I\\
I & \sum_{i=1}^N \omega_i H_i^\top X_i^{-1}H_i
\end{bmatrix}\succeq 0,
\end{aligned}
\]
with fusion gain
\[
K^\star
=
\left(\sum_{i=1}^N \omega_i^\star H_i^\top X_i^{-1}H_i\right)^{-1}
\left[
\omega_1^\star H_1^\top X_1^{-1}\ \cdots\ \omega_N^\star H_N^\top X_N^{-1}
\right].
\]
When specialized to multiple full-state estimates, this formulation recovers the state-of-the-art family-optimal CI solutions previously reported for standard CI [2603.20402].

This optimization view clarifies two points. First, WCI is not merely a heuristic convex combination; it is a robust design against all admissible cross-correlation structures. Second, once the criterion \(J\) is fixed, weight selection can be posed as a convex SDP, rather than as an ad hoc line search. This has become particularly important in large-scale distributed estimation and cooperative localization, where off-the-shelf SDP solvers can be used for real-time or near real-time implementation [2603.20402].

## 3. Information-theoretic and structured generalizations

In Gaussian single-object estimation, WCI is also the single-object manifestation of a broader information-theoretic construction. In multiobject random finite set fusion, generalized covariance intersection (GCI) is the multiobject counterpart of CI/WCI and is exactly the weighted geometric mean of local densities. It solves a weighted information gain minimization problem,
\[
\overline f_{\rm GCI}=\arg\min_f \sum_i \omega_i D_{\rm KL}(f\|f^i),
\]
and in the single-object Gaussian case reduces to the standard CI/WCI formulas above [1903.04239].

The same literature also identifies a dual construction based on weighted information loss,
\[
\overline f_{\rm MIL}=\arg\min_f \sum_i \omega_i D_{\rm KL}(f^i\|f),
\]
whose unconstrained solution is the linear opinion pool, or arithmetic pooling. This dual rule is not standard WCI, but it is useful because it shows that WCI belongs to the weighted information-gain, or geometric-pooling, branch of conservative fusion, whereas arithmetic pooling belongs to a different design objective [1903.04239].

A separate generalization concerns structure in the state itself. Standard WCI is “monolithic” in the sense that it uses one global weight for the entire state. When probabilistic independence structure is available, tighter conservative bounds can be obtained by using multiple, non-monolithic weights. In one formulation, independent sub-states are fused with block-specific weights, and a general robust optimization over the unknown cross-covariance set yields tighter bounds than monolithic CI. On a simple problem, the non-monolithic CI solution and the general optimization scheme converge to the same solution; in large-scale target tracking, the non-monolithic algorithm yields a tighter bound and a more accurate estimate than original monolithic CI [2307.10594].

## 4. Sequential, split, and recursive forms

Several important variants of WCI address recursive fusion and implementation efficiency. In split covariance intersection filtering, each local covariance is decomposed into dependent and independent parts,
\[
P_i=P_{id}+P_{ii},
\]
and the effective covariances are defined as
\[
P_1(w)=P_{1d}/w+P_{1i},\qquad
P_2(w)=P_{2d}/(1-w)+P_{2i},
\]
with fused covariance
\[
P(w)=\bigl(P_1(w)^{-1}+P_2(w)^{-1}\bigr)^{-1}.
\]
The associated \(w\)-optimization problem
\[
w^\star=\arg\min_{w\in[0,1]}\det P(w)
\]
has been proved convex via convexity of \(\log\det P(w)\), which is a practically significant result because it guarantees that any local minimum is global and makes one-dimensional optimization straightforward [2101.10159].

A second line of work concerns sequential CI. Enhanced sequential covariance intersection (ESCI) generalizes both batch CI and one-at-a-time SCI by allowing multiple estimate pairs to be fused at each sequential step. The recursion is
\[
\left\{
\begin{array}{l}
x^e_i
=
w^e_i P^e_i (P^e_{i-1})^{-1}x^e_{i-1}
+
\sum_{j=b_{i-1}+1}^{b_i} w^e_{i,j} P^e_i P^{-1}_{r_j}x_{r_j},\\[4pt]
P^e_i
=
\Big(
w^e_i (P^e_{i-1})^{-1}
+
\sum_{j=b_{i-1}+1}^{b_i} w^e_{i,j}P^{-1}_{r_j}
\Big)^{-1},
\end{array}
\right.
\]
and it remains unbiased and consistent [2110.06598].

ESCI becomes especially notable when paired with an analytic weighting criterion. If the local importance of estimate pair \(\{x_i,P_i\}\) is represented by a positive scalar \(f(\{x_i,P_i\})\), then choosing
\[
w^e_i=\frac{\sum_{k=1}^{b_{i-1}} f(\{x_{r_k},P_{r_k}\})}{\sum_{k=1}^{b_i} f(\{x_{r_k},P_{r_k}\})},\qquad
w^e_{i,j}=\frac{f(\{x_{r_j},P_{r_j}\})}{\sum_{k=1}^{b_i} f(\{x_{r_k},P_{r_k}\})}
\]
makes the final fusion result independent of the fusion order and of the number of estimates fused in each sequential step. Examples given for \(f\) include \(1/\operatorname{tr}(P_i)\), \(1/\det(P_i)\), \(\operatorname{tr}(P_i^{-1})\), \(\det(P_i^{-1})\), and \(1/\operatorname{tr}(DP_i)\) with diagonal \(D\) for component-prioritized fusion [2110.06598].

## 5. Applications in distributed localization and multitarget fusion

WCI is especially prominent in cooperative localization. In distributed visual-inertial cooperative localization, CI is used so that each robot estimates only its own state and autocovariance while compensating for unknown correlations between robots. In one implementation, the CI prior covariance is taken as
\[
P_{\text{CI},k}:=\operatorname{Diag}\left(\frac{1}{\omega_i}P_{ii,k},\frac{1}{\omega_1}P_{11,k},\dots,\frac{1}{\omega_L}P_{LL,k}\right)\ge P_k,
\]
and fixed weights such as \(\omega_o=0.001\), \(\omega_o=0.005\), \(\omega_i=0.99\), and \(\omega_i=0.995\) were used in different update modes, with the authors noting that these weights could be found by minimizing the trace or determinant of the posterior covariance but were empirically tuned for consistent performance [2103.12770].

A related multirobot localization architecture applies CI only in the communication update while keeping standard EKF prediction and observation updates locally. In information form, robot \(i\) fuses incoming estimates by
\[
\Sigma_{i,t^+}^{-1}=\sum_{j\in C_{i,t}} c_j I^j_{i,t}+c_i \Sigma_{i,t}^{-1},
\]
\[
\bar s_{i,t^+}=\Sigma_{i,t^+}\left(\sum_{j\in C_{i,t}} c_j \bar e^j_{i,t}+c_i\Sigma_{i,t}^{-1}\bar s_{i,t}\right),
\]
with convex coefficients \(c_j\). Under a weak connectivity condition on the union of observation links in the super-neighborhood of robot \(i\), the local covariance is bounded. Empirically, this architecture is less affected by blocked or partially unavailable communication than several local-state alternatives [2108.08789].

A more specialized 3D range-based distributed cooperative localization formulation introduces an explicit weighting matrix into CI. Instead of minimizing \(\operatorname{tr}(P)\), the method selects \(\omega\) by
\[
\omega^\star=\arg\min_{\omega\in[0,1]}\operatorname{tr}(W_{i,k}P_{i,k}(\omega)),
\]
where \(W_{i,k}=S^\top S\) is derived from INS error propagation over a future horizon \(T_a\). This design addresses scale imbalance, correlation mismatch, and lack of task relevance in classical CI criteria. In the reported 3D simulations, agent 1 position RMSE was \(0.190\) m for WCI versus \(0.223\) m for CI-trace and \(0.458\) m for CI-det, while distributed operation remained more robust and scalable than centralized fusion [2508.12207].

In multitarget tracking, WCI ideas reappear in random finite set form. One paper develops adaptive, information-driven GCI for labeled multi-Bernoulli filtering, in which each Bernoulli component receives sensor-specific adaptive weights computed from a Cauchy–Schwarz divergence between prior and posterior Bernoulli components. The resulting weighted GCI is explicitly described as a multi-object, RFS-based Weighted Covariance Intersection whose weights are label-dependent and time-varying [1702.08641].

## 6. Conservatism, criticisms, and alternatives

The principal criticism of WCI is conservatism. Because it is designed to remain valid for all admissible unknown cross-correlations, its covariance bound can be substantially larger than necessary. A game-theoretic alternative formulates fusion as a minimax problem over the fusion matrix \(K\) and the unknown cross-covariance \(\Sigma_{xy}\),
\[
\min_K \sup_{\Sigma_{xy}} \operatorname{tr}(\Sigma_{zz}(K,\Sigma_{xy}))
\quad
\text{s.t. }
\begin{bmatrix}
\Sigma_{xx} & \Sigma_{xy}\\
\Sigma_{xy}^T & \Sigma_{yy}
\end{bmatrix}\succeq 0,
\]
and produces a trace-consistent estimator that is less conservative than CI in the reported examples, although it relaxes full matrix consistency to trace-consistency [1610.01045].

A Bayesian alternative assumes a prior distribution on the full covariance matrix, derives the conditional distribution of the unknown off-diagonal blocks given the known diagonal blocks, and uses Monte Carlo integration to approximate the MMSE fusion rule. In the two-node case the conditional distribution of the cross-covariance is an inverted matrix variate \(t\)-distribution. The reported simulations show that this Bayesian method works better than the popular covariance intersection method, but it replaces CI’s worst-case guarantee with a prior-dependent MMSE construction [1303.2414].

From a set-membership viewpoint, the Convex Combination Ellipsoid (CCE) method has also been proposed as a less conservative covariance estimate than CI. CCE uses the same scalar weight \(\alpha\in[0,1]\) as CI and the same fused center, but scales the CI base covariance by a factor \(k=1-d^2\). The paper emphasizes two properties that CI lacks in general: CCE does not introduce additional uncertainty that was not already present in the prior estimates, and its ellipsoid is contained in the union of the prior ellipsoids while still containing their intersection [2304.07996].

These alternatives do not eliminate the role of WCI. Rather, they clarify its place: WCI remains the canonical conservative fusion rule when unknown cross-correlations must be handled without modeling assumptions, whereas structured, Bayesian, game-theoretic, or set-membership alternatives are used when tighter bounds are sought and additional assumptions or different optimality criteria are acceptable [1610.01045].

Source: https://www.emergentmind.com/topics/weighted-covariance-intersection-wci