Kalman filter for radio source power and direction of arrival estimation

Images are an important source of information for spacecraft navigation. Based on an image and a known attitude, triangulation techniques (intersection or resection) are often used for positioning and navigation. In the resection problem, an observer estimate its unknown location by using angle measurements to points at known locations (i.e., landmarks), the localization performance depending on the accuracy of the angle measurements. As a contribution to resection for spacecraft navigation, we considers the dynamic image estimation problem based on radio interferometry, i.e., image of radio source power, where the measurements are sample covariance matrices (SCMs). Considering the case where several measurements are available as well as a known dynamic linear model of image evolution, a.k.a a linear state model, the minimum mean-squared error image estimator (MMSE) is given by the Kalman filter (KF) or one of its variants. However standard Kalman-like filters are not a priori suitable for the problem at hand since the measurements (i.e., SCMs) cannot be formulated analytically as a function of state parameters to be estimated (i.e., radio source power). In fact, this lack of analytical formulation can be circumvented by a statistical linear fitting allowing the SCMs to be expressed in terms of the state. This linear fitting introduces an additive residual noise, equivalent to a measurement noise, whose covariance matrix depends on the current state, a non-standard case for a measurement model. The covariance matrix of the residual noise is derived whatever the distributions of the radio sources and of the additive noise at the samples level, unveiling the contribution of their multivariate kurtosis. The proposed method is evaluated on simulated data representative of a dynamic radio interferometry framework. The results show that the proposed method is capable of effectively tracking moving radio sources in complex scenes with theoretical guaranties when the signal multivariate kurtosis is known.