Authors: Behnood Rasti, Bikram Koirala, Paul Scheunders
Organizations: Faculty of Electrical Engineering and Computer Science, Technische Universität Berlin, 10623 Berlin, Germany · Berlin Institute for the Foundations of Learning and Data (BIFOLD), 10623 Berlin, Germany · Imec-VisionLab, Department of Physics, University of Antwerp, Belgium
This paper proposes a semisupervised geometric unmixing approach called minimum simplex semisupervised unmixing (MiSiSUn). The geometry of the data was incorporated for the first time into library-based unmixing using a simplex-volume-flavored penalty based on an archetypal analysis-type linear model. The experimental results were performed on two simulated datasets considering different levels of mixing ratios and spatial instruction at varying input noise. MiSiSUn considerably outperforms state-of-the-art semisupervised unmixing methods. The improvements vary from 1 dB to over 3 dB in different scenarios. The proposed method was also applied to a real dataset where visual interpretation is close to the geological map. MiSiSUn was implemented using PyTorch, which is open-source and available at https://github.com/BehnoodRasti/MiSiSUn. Moreover, we provide a dedicated Python package for Semisupervised Unmixing, which is open-source and includes all the methods used in the experiments for the sake of reproducibility.
TABLE I: Hyperparameters used in the experiments in different scenarios.
Fig. 1: Sim1 dataset: a) Band number 70 (1050 nm) b) Endmembers.
Fig. 2: Endmembers used to simulate the Sim2 dataset.
Fig. 3: Abundance SRE in dB obtained by the selected state-of-the-art methods on the Sim1 dataset for different noise levels.
Fig. 4: Estimated abundances by the selected state-of-the-art methods on the Sim1 dataset for a given input SNR (20 dB).
Fig. 5: Abundance SRE in dB obtained by the selected state-of-the-art methods on the Sim2 dataset at a given input SNR (30 dB) for different pixel purity levels.
Fig. 6: Abundance maps obtained by the different methods on Cuprite, along with the geological reference map and the produced GT abundance maps from this map.
Fig. 7: Endmembers obtained by the different low-rank methods to Cuprite, compared with the USGS spectral library. The ground-truth endmembers are depicted in red, the estimated ones in black.
SAD (in degree)
SUnAA
FaSUn
SUns
MiSiSUn
Alunite
7.68
1.48
3.68
5.97
Chalcedony
6.09
3.43
3.06
7.33
Kaolinite
2.46
2.67
16.59
2.78
Montmorillonite
4.83
6.55
3.91
6.53
Sphene
8.26
16.23
10.92
3.45
TABLE II: SAD between estimated and library endmembers for the Cuprite dataset. The Best Performances are Shown in Bold.
# Pixels
# Bands
# Endmembers
# Atoms
SUnSAL
S2WSU
CLSUnSAL
MUA_SLIC
SUnCNN
SUnAA
FaSUn
SUnS
MiSiSUn
Sim1
11025
224
6
498
41.6
107.4
55.8
37.1
68.0
189.5
126.4
129.5
122.0
Sim2
10000
224
6
498
20.5
84.7
56.9
38.9
60.0
159.0
111.3
115.3
113.2
Cuprite
47750
188
12
498
47.7
321.4
175.6
116.9
565.8
195.5
134.4
134.0
138.7
Sim1_extend
99225
224
6
498
339.6
884.2
341.8
325.2
477.4
1631.4
148.9
176.6
152.5
TABLE III: Processing time. The best results are in boldface, and the second best are underlined.
Semantic segmentation and hyperspectral unmixing are two central problems in spectral image analysis. The former assigns each pixel a discrete label corresponding to its material class, whereas the latter estimates pure material spectra, called endmembers, and, for each pixel, a vector representing material abundances in the observed scene. Despite their complementarity, these two problems are usually addressed independently. This paper aims to bridge these two lines of work by formally showing that, under the linear mixing model, pixel classification by dominant materials induces polyhedral-cone regions in the spectral space. We leverage this fundamental property to propose a direct segmentation-to-unmixing pipeline that performs blind hyperspectral unmixing from any semantic segmentation by constructing a polyhedral-cone partition of the space that best fits the labeled pixels. Signed distances from pixels to the estimated regions are then computed, linearly transformed via a change of basis in the distance space, and projected onto the probability simplex, yielding an initial abundance estimate. This estimate is used to extract endmembers and recover final abundances via matrix pseudo-inversion. Because the segmentation method can be freely chosen, the user gains explicit control over the unmixing process, while the rest of the pipeline remains essentially deterministic and lightweight. Beyond improving interpretability, experiments on three real datasets demonstrate the effectiveness of the proposed approach when associated with appropriate clustering algorithms, and show consistent improvements over recent deep and non-deep state-of-the-art methods. The code is available at: https://github.com/antoine-bottenmuller/polyhedral-unmixing
Antoine Bottenmuller, Etienne Decencière, Petr Dokládal
CMM, PSL, STIM · Mines Paris, PSL University, Centre for Mathematical Morphology (CMM), 77300 Fontainebleau, France
Signal separation problems are common in science. A prominent example of this occurs during the use of diffraction or spectroscopy to identify the individual components of a mixture by measuring it. In the simplest case, the measured signal is a linear combination of basis patterns corresponding to the constituent parts. The unmixing problem is to infer all or some of these basis patterns and abundances of components from measurements of distinct mixtures. One of the core challenges of this task is the variation of the basis from mixture to mixture due to noise and the exact physics of the measurement process. This is usually addressed with tailored model-based and parametric methods that are then limited in use to specific application domains by the nature of the assumptions made. We propose a novel geometric approach to unmixing problems which views the generation of data during measurement through a metric space lens, thereby shifting the focus from parametrised models to a general relationship between basis transformations and the corresponding geometry. We take advantage of the optimal transport distances to capture commonly occurring basis variations, and use minimisation of in-class variance of candidate solutions to drive the optimisation. We pay special attention to the one-dimensional case due to its practical importance and availability of efficient distance and transport map routines. The effectiveness of our approach is demonstrated on a range of unmixing tasks using random Gaussian mixture models, simulated powder X-ray diffraction, and laboratory hyperspectral imaging datasets.
Dumitru Mirauta, Vladimir V. Gusev, Michael W. Gaultois +2
The Linear Mixing Model (LMM) dominates spectral unmixing for its simplicity, but fails under multiple scattering; existing nonlinear models compensate by applying a fixed regime uniformly across entire scenes. We propose Physics-Guided Regime Unmixing (PGRU), which estimates a pixel-wise scalar ξi∈[0,1] from observable physical features to activate nonlinear mixing only where justified. Residuals from the Generalized Bilinear Model (GBM), the Post-Nonlinear Mixing Model (PPNM), and Hapke are combined via learned attention, yielding interpretable regime maps. Experiments on Samson, Jasper Ridge, and Urban show consistent improvements over baselines, with physical coherence ρ>0.90.
Paula Pacheco, Pablo Granitto, Juan B. Cabral
GVT-CONAE, Centro Espacial Teófilo Tabanera, Argentina · CONICET, Argentina · FAMAF-UNC, Córdoba, Argentina +1