Hyperspectral reconstruction
I reconstructed hyperspectral cubes from SEM-EDS images of meteorites, where most voxels hold zero or one photon. I showed that the usual data-fidelity term is degenerate at these counts and replaced it with a hard box constraint, which makes the problem convex with no trade-off weight. A regularizer guided by the electron image keeps grain boundaries sharp. I ran it in PyTorch on a GPU, on a real chondrite cube of1024 × 1024 pixels × 2048 energy channels.
| Metric | Data | Result |
|---|---|---|
| PSNR against the ground truth | Synthetic, Poisson noise | +1.5 dB over raw counts |
| Spectral angle (SAM) | Synthetic, Poisson noise | 22.3° to 18.9° |
| Spectral dispersion between neighboring pixels | Real chondrite cube | 65% lower |