Stationary subspace analysis for spatial data
Abstract
Stationary subspace analysis (SSA) is a blind source separation framework that decomposes linearly mixed multivariate data into stationary and nonstationary components. We extend SSA to spatially indexed data by introducing spatial stationary subspace analysis (spSSA), which explicitly accounts for spatial dependence. We propose three estimation procedures for the unmixing matrix based on first- and second-order spatial statistics. Each procedure targets a different type of nonstationarity and can be formulated as the solution to a generalized eigenvalue problem. To address situations where multiple forms of nonstationarity are present simultaneously, we combine the three procedures using approximate joint diagonalization. Simulation studies demonstrate that this combined approach yields superior separation performance. When the dimension of the nonstationary subspace is known, the proposed methods reliably recover the latent stationary and nonstationary components. However, determining this dimension remains a fundamental challenge in SSA, for which no generally accepted solution currently exists. Building on our estimation procedures, we propose a novel data augmentation approach to estimate the dimension of the nonstationary subspace and demonstrate its effectiveness through simulation studies. The proposed methodology is easily transferable to time series settings, making it of broader methodological interest.
Versions
Related
- Subspace concentration condition and Bogomolov--Gieseker inequality
- A proof of the Braun-Etzion-Vardy bound for binary subspace codes
- Discrete homotopy of graphs and k-parabolic subspace arrangements
- Bias-Corrected Subspace Intersection: Minimax-Optimal Shared Subspace Estimation in Multi-View Data
- A complex associated to a family of ideals and applications to subspace arrangements
BibTeX
@misc{saarela2026stationary,
title = {{Stationary subspace analysis for spatial data}},
author = {Perttu Saarela and Klaus Nordhausen and Jaakko Pere and Anne M. Ruiz},
howpublished = {arXiv},
year = {2026},
month = may,
doi = {10.48550/arxiv.2605.19861},
eprint = {2605.19861},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2605.19861},
}
RIS
TY - UNPB
TI - Stationary subspace analysis for spatial data
AU - Perttu Saarela
AU - Klaus Nordhausen
AU - Jaakko Pere
AU - Anne M. Ruiz
PY - 2026
DA - 2026/05/19
PB - arXiv
DO - 10.48550/arxiv.2605.19861
UR - https://arxiv.org/abs/2605.19861
AB - Stationary subspace analysis (SSA) is a blind source separation framework that decomposes linearly mixed multivariate data into stationary and nonstationary components. We extend SSA to spatially indexed data by introducing spatial stationary subspace analysis (spSSA), which explicitly accounts for spatial dependence. We propose three estimation procedures for the unmixing matrix based on first- and second-order spatial statistics. Each procedure targets a different type of nonstationarity and can be formulated as the solution to a generalized eigenvalue problem. To address situations where multiple forms of nonstationarity are present simultaneously, we combine the three procedures using approximate joint diagonalization. Simulation studies demonstrate that this combined approach yields superior separation performance. When the dimension of the nonstationary subspace is known, the proposed methods reliably recover the latent stationary and nonstationary components. However, determining this dimension remains a fundamental challenge in SSA, for which no generally accepted solution currently exists. Building on our estimation procedures, we propose a novel data augmentation approach to estimate the dimension of the nonstationary subspace and demonstrate its effectiveness through simulation studies. The proposed methodology is easily transferable to time series settings, making it of broader methodological interest.
ER -