pith. sign in

arxiv: 2306.08793 · v5 · pith:DKURR3URnew · submitted 2023-06-15 · ⚛️ physics.optics · physics.bio-ph

Adaptive Optical Multi-Spectral Matrix Approach for Label-free High-resolution Imaging through Complex Scattering Media

Yiwen Zhang , Minh Dinh , Zeyu Wang , Tianhao Zhang , Tianhang Chen , Chia Wei Hsu This is my paper

Pith reviewed 2026-05-24 08:43 UTC · model grok-4.3

classification ⚛️ physics.optics physics.bio-ph
keywords scattering matrix tomographyadaptive opticsimaging through scattering medialabel-free imagingvolumetric imagingnonconvex optimizationZernike regularization
0
0 comments X

The pith

Scattering Matrix Tomography inverts multi-spectral matrices with Zernike-regularized nonconvex optimization to recover object fields through strong high-order aberrations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper presents Scattering Matrix Tomography (SMT) as a method that reformulates imaging through complex scattering media into a numerical optimization problem solved from spectrally resolved scattering matrix measurements. SMT applies Zernike-mode regularization and a coarse-to-fine nonconvex strategy to undo aberrations that defeat standard adaptive optics and matrix approaches. This yields label-free volumetric imaging at depths over three transport mean free paths in colloids and a depth-to-resolution ratio above 900 in ex vivo mouse brain tissue. A reader would care because the technique operates noninvasively both inside and outside the media, extending high-resolution access into regimes where multiple scattering previously prevented recovery of fine detail.

Core claim

SMT reformulates imaging through complex media as a numerical optimization and employs Zernike-mode wavefront regularization and coarse-to-fine nonconvex optimization strategy to reverse severe aberrations, enabling noninvasive high-resolution volumetric imaging in multiple scattering regime based on spectrally-resolved matrix measurement.

What carries the argument

Scattering Matrix Tomography (SMT), which inverts the spectrally-resolved scattering matrix via Zernike-regularized nonconvex optimization to recover the true object field under high-order aberrations.

If this is right

  • SMT enables volumetric imaging inside opaque media where conventional adaptive optics and matrix methods cannot correct strong aberrations.
  • The approach achieves a depth-over-resolution ratio above 900 beneath ex vivo mouse brain tissue.
  • SMT operates noninvasively and label-free both inside and outside the scattering media.
  • It supports applications in medical imaging, biological science, device inspection, and colloidal physics.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the multi-wavelength matrix inversion generalizes, the same optimization could be tested on dynamic or thicker samples to map the practical depth limit.
  • The method's reliance on spectral diversity suggests potential for combining with other spectral channels to reduce the number of required measurements.

Load-bearing premise

The scattering matrix measured at multiple wavelengths can be inverted via Zernike-regularized nonconvex optimization to recover the true object field even when high-order aberrations dominate.

What would settle it

A controlled test in which the object field is known in advance but the SMT reconstruction fails to match it at the reported depths in mouse brain tissue or colloid while conventional methods also fail would falsify the claim.

read the original abstract

Imaging through complex scattering media is severely limited by aberrations and scattering which obscure images and reduce resolution. Confocal and temporal gatings partly filter out multiple scattering but are severely degraded by wavefront distortions. Adaptive optics restore resolution by correcting low-order aberrations and matrix-based imaging enables more complex wavefront corrections. However, they struggle to undo high-order aberrations under strong scattering, preventing imaging at greater depths. To address these challenges, we present Scattering Matrix Tomography (SMT), an approach that makes full use of the wavefront engineering capability of scattering matrix and extreme adaptive optics. SMT reformulates imaging through complex media as a numerical optimization and employs Zernike-mode wavefront regularization and coarse-to-fine nonconvex optimization strategy to reverse severe aberrations, enabling noninvasive high-resolution volumetric imaging in multiple scattering regime. Based on the spectrally-resolved matrix measurement, SMT achieves a depth-over-resolution ratio above 900 beneath $ex~vivo$ mouse brain tissue and volumetric imaging at over three transport mean free paths inside an opaque colloid, where conventional methods fail to correct strong aberrations under these challenging conditions. SMT is noninvasive, label-free, and works both inside and outside the scattering media, making it suitable for various applications, including medical imaging, biological science, device inspection, and colloidal physics.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 0 minor

Summary. The paper introduces Scattering Matrix Tomography (SMT), which reformulates imaging through scattering media as a numerical optimization problem solved via Zernike-mode wavefront regularization and a coarse-to-fine nonconvex strategy applied to spectrally-resolved scattering matrix measurements. It claims this enables noninvasive, label-free high-resolution volumetric imaging at depths where conventional adaptive optics and matrix methods fail, with reported performance of depth-over-resolution ratio >900 beneath ex vivo mouse brain tissue and volumetric imaging at >3 transport mean free paths in an opaque colloid.

Significance. If the central claims hold, SMT would constitute a meaningful advance over existing adaptive-optics and scattering-matrix techniques by addressing high-order aberrations in the multiple-scattering regime. The experimental demonstrations in biological tissue and colloidal samples, together with the noninvasive and label-free character, would be relevant to biomedical imaging and colloidal physics.

major comments (2)
  1. [Abstract (SMT reformulation and optimization strategy)] The section describing the SMT reformulation and optimization strategy (Zernike-mode regularization plus coarse-to-fine nonconvex optimization applied to the multi-spectral matrix) provides no analysis of solution uniqueness, basin of attraction, or regularization strength under dominant high-order aberrations. This is load-bearing for the claim that the recovered field corresponds to the true object rather than a data-consistent artifact at depths >3 tmfp.
  2. [Abstract (performance claims)] No quantitative ground-truth recovery metrics, error bars, or comparison against known phantoms are supplied for the >3 tmfp colloid experiment or the depth-over-resolution ratio of 900; the performance numbers are presented as direct experimental outcomes without the supporting validation data needed to substantiate the central claims.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the constructive comments, which help clarify the presentation of our optimization approach and experimental validation. We address each major comment below and outline planned revisions.

read point-by-point responses

  1. Referee: [Abstract (SMT reformulation and optimization strategy)] The section describing the SMT reformulation and optimization strategy (Zernike-mode regularization plus coarse-to-fine nonconvex optimization applied to the multi-spectral matrix) provides no analysis of solution uniqueness, basin of attraction, or regularization strength under dominant high-order aberrations. This is load-bearing for the claim that the recovered field corresponds to the true object rather than a data-consistent artifact at depths >3 tmfp.

    Authors: The coarse-to-fine strategy combined with Zernike regularization is intended to guide the optimizer toward physically plausible solutions by progressively incorporating higher-order modes while penalizing unphysical wavefronts. Empirical validation is provided through consistent recovery of known structures across spectral channels and direct comparison to conventional matrix methods that fail under the same conditions. We agree a dedicated discussion of regularization strength selection (via L-curve or cross-validation) and simulation-based exploration of the basin of attraction would strengthen the manuscript; this will be added as a new subsection in the Methods together with supplementary simulation results. revision: yes

  2. Referee: [Abstract (performance claims)] No quantitative ground-truth recovery metrics, error bars, or comparison against known phantoms are supplied for the >3 tmfp colloid experiment or the depth-over-resolution ratio of 900; the performance numbers are presented as direct experimental outcomes without the supporting validation data needed to substantiate the central claims.

    Authors: The reported depth-over-resolution ratio is obtained from the measured tissue thickness divided by the experimentally determined lateral resolution (FWHM of resolved features), and the >3 tmfp depth is calculated from independently measured scattering parameters of the colloid. Resolution is further corroborated by the ability to separate particles whose size and spacing are known a priori. We will augment the Results and Supplementary Information with error bars from repeated acquisitions, quantitative metrics (e.g., structural similarity indices on simulated phantoms), and additional phantom comparisons to make the validation more explicit. revision: yes

standing simulated objections not resolved
  • Rigorous mathematical proof of global solution uniqueness for the nonconvex optimization under arbitrary high-order aberrations (such proofs are generally intractable for this class of problems; we can only provide empirical and simulation-based support).

Circularity Check

0 steps flagged

No circularity; experimental outcomes independent of any self-referential derivation

full rationale

The paper presents SMT as an experimental method that acquires spectrally-resolved scattering matrices and applies Zernike-regularized nonconvex optimization to reconstruct object fields. Reported metrics (depth-over-resolution >900, imaging >3 tmfp) are direct experimental measurements on ex vivo brain tissue and colloids, not quantities derived from the paper's own equations or fitted parameters. No self-definitional loop, fitted-input-called-prediction, or load-bearing self-citation chain exists in the described chain; the optimization is a computational procedure whose validity is assessed by external experimental benchmarks rather than internal consistency alone. The derivation is therefore self-contained against external data.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review prevents exhaustive enumeration; the method implicitly relies on the existence of a measurable scattering matrix and the convergence of the stated optimizer, but no explicit free parameters or invented entities are named.

pith-pipeline@v0.9.0 · 5770 in / 1147 out tokens · 21418 ms · 2026-05-24T08:43:34.349081+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Taking advantage of multiple scattering for Optical Reflection Tomography

    physics.optics 2025-05 unverdicted novelty 7.0

    A novel optimization algorithm reconstructs 3D refractive index in reflection-mode microscopy by exploiting multiply-scattered waves from uncontrolled background structures, using weighted time loss, positivity constr...

  2. Label-free subcellular 3D imaging of oocytes and embryos via reflection matrix microscopy

    physics.optics 2024-10 unverdicted novelty 6.0

    A multi-spectral reflection matrix imaging platform achieves label-free subcellular 3D visualization of oocytes and blastocysts at 300 nm resolution by computationally correcting aberrations and multiple scattering.

Reference graph

Works this paper leans on

45 extracted references · 45 canonical work pages · cited by 2 Pith papers · 1 internal anchor

  1. [1]

    minimize defocus error

    VRM spectral dispersion compensation The very first step of VRM is to “minimize defocus error” through “a numerical propagation to the depth with the maximum reflectance (zmax)” [22]: Rmax(kout ∥, kin ∥,ω)≡R(kout ∥, kin ∥,ω)ei(kout z −kin z )zmax. (S43) Doing so places the reference plane of the reflection matrix at zmax. This step is only mentioned on pa...

    work page

  2. [2]

    Meanwhile, VRM uses ϕin(kin ∥,ω) and ϕout(kout ∥,ω) as free variables for every pair of ( kin ∥,ω) and every pair of ( kout ∥,ω)

    To avoid overfitting and to alleviate the curse of dimensionality in an optimization problem, SMT employs only three variables for dispersion compensation: ˜a1, ˜a2, ˜a3, each having a physical meaning (group velocity, group velocity dispersion, and asymmetric pulse distortion). Meanwhile, VRM uses ϕin(kin ∥,ω) and ϕout(kout ∥,ω) as free variables for eve...

    work page

  3. [3]

    (S20) sums over a 3D volume, which avoids overfitting

    The spatial summation of SMT dispersion compensation in Eq. (S20) sums over a 3D volume, which avoids overfitting. The spatial summation of VRM dispersion compensation in Eq. (S49) sums over a 2D slice at a single depth zmax, which can create serious artifacts. This is the reason that the volumetric VRM image in Fig. 5 of the main text shows bright slices...

    work page

  4. [4]

    (2) of the main text) used in the SMT objective function of Eq

    The SMT image ( ISMT, defined in Eq. (2) of the main text) used in the SMT objective function of Eq. (S20) employs both input spatial gating (through a summation over kin) and output spatial gating (through a summa- tion over kout). The confocal spatial gate is crucial for the optimization, as shown in Supplementary Fig. 14. Meanwhile, theIin andIout of E...

    work page

  5. [5]

    temporal truncation

    The SMT image ISMT used in the SMT dispersion compensation also employees temporal gating through a sum- mation over all frequencies ω. The temporal gate is crucial for the optimization, as shown in Supplementary Fig. 14. The VRM frequency summation in Eqs. (S50)–(S51) only sums over two frequencies, ω′and ωc, so the temporal gate is gone, which further r...

    work page

  6. [6]

    closed-loop accumulation of single-scattering

    VRM angular wavefront correction and reconstruction After estimating the spectro-angular dispersion relative to the center frequencyωc, VRM performs another iterative phase conjugation to estimate the angular wavefront distortion ϕin(kin ∥,ωc) and ϕout(kout ∥,ωc) at ωc, and then adds the same wavefront correction of−ϕin(kin ∥,ωc) and−ϕout(kout ∥,ωc) to al...

    work page

  7. [7]

    In CLASS and VRM, every ϕin(kin ∥) and every ϕout(kout ∥) is a free variable

    To avoid overfitting, SMT regularize the problem by expanding ϕin(kin ∥) and ϕout(kout ∥) in Zernike polynomials and only using the Zernike coefficients as the free variables. In CLASS and VRM, every ϕin(kin ∥) and every ϕout(kout ∥) is a free variable

    work page

  8. [8]

    VRM starts with the ϕin(kin ∥,ω) and ϕout(kout ∥,ω) from its dispersion compensation step

    SMT starts from an initial guess of no wavefront correction, ϕin(kin ∥) = ϕout(kout ∥) = 0. VRM starts with the ϕin(kin ∥,ω) and ϕout(kout ∥,ω) from its dispersion compensation step

    work page

  9. [9]

    There is no such progression in CLASS and VRM

    To avoid getting trapped in a poor local optimum, SMT first optimizes ϕin/out over the full volume and then progressively optimizes over smaller zones while using the previous ϕin/out as the initial guess. There is no such progression in CLASS and VRM. As shown in Fig. 2c–f of the main text, optimizing the wavefront inside smaller zones without regulariza...

    work page

  10. [10]

    FOV in air

    SMT also increases the number of Zernike terms when progressing to smaller zones to account for the relation between isoplanatic patch size and the smoothness of wavefront variation. This is not possible in CLASS and VRM. 30 While the VRM wavefront correction uses multiple slices, it cannot reverse the over-optimization of the intensity at the single slic...

    work page

  11. [11]

    Verrier and M

    N. Verrier and M. Atlan, Off-axis digital hologram reconstruction: some practical considerations, Appl. Opt. 50, H136 (2011)

    work page 2011

  12. [12]

    H. Yu, T. R. Hillman, W. Choi, J. O. Lee, M. S. Feld, R. R. Dasari, and Y. Park, Measuring large optical transmission matrices of disordered media, Phys. Rev. Lett. 111, 153902 (2013)

    work page 2013

  13. [13]

    C. Wang, X. Lee, T. Cui, Y. Qu, Y. Li, H. Li, and Q. Wang, Reducing the influence of direct reflection on return signal detection in a 3 D imaging lidar system by rotating the polarizing beam splitter, Appl. Opt. 55, 1559 (2016)

    work page 2016

  14. [14]

    Phantom TMX 7510, https://www.phantomhighspeed.com/products/cameras/tmx/7510

    work page

  15. [15]

    Refractive index of SCHOTT multiple purpose D263TECO, https://refractiveindex.info/?shelf=glass&book= SCHOTT-multipurpose&page=D263TECO

    work page

  16. [16]

    Daimon and A

    M. Daimon and A. Masumura, Measurement of the refractive index of distilled water from the near-infrared region to the ultraviolet region, Appl. Opt. 46, 3811 (2007)

    work page 2007

  17. [17]

    Gupta, P

    V. Gupta, P. T. Probst, F. R. Goßler, A. M. Steiner, J. Schubert, Y. Brasse, T. A. F. K¨ onig, and A. Fery, Mechanotunable 36 1470152515800.0FWHMalongx(µm)1.51.00.5 147015251580z(µm) 1.0FWHMalongz(µm)1.52.0 0.5 1470152515800.0FWHMalongy(µm)1.51.00.5 Supplementary Fig. 24. Experimental FWHM resolution of nanoparticle colloid. Each experimental data point i...

    work page 2019

  18. [18]

    M. N. Polyanskiy, Refractive index database, https://refractiveindex.info

    work page

  19. [19]

    Wojtkowski, V

    M. Wojtkowski, V. J. Srinivasan, T. H. Ko, J. G. Fujimoto, A. Kowalczyk, and J. S. Duker, Ultrahigh-resolution, high- speed, Fourier domain optical coherence tomography and methods for dispersion compensation, Opt. Express 12, 2404 (2004)

    work page 2004

  20. [20]

    S. G. Johnson, The NLopt nonlinear-optimization package, http://github.com/stevengj/nlopt

    work page

  21. [21]

    D. C. Liu and J. Nocedal, On the limited memory BFGS method for large scale optimization, Math. Program. 45, 503 (1989)

    work page 1989

  22. [22]

    Fricker, Zernike polynomials, https://www.mathworks.com/matlabcentral/fileexchange/ 7687-zernike-polynomials (2023), MATLAB Central File Exchange

    P. Fricker, Zernike polynomials, https://www.mathworks.com/matlabcentral/fileexchange/ 7687-zernike-polynomials (2023), MATLAB Central File Exchange

    work page 2023

  23. [23]

    Calzavara-Pinton, C

    P. Calzavara-Pinton, C. Longo, M. Venturini, R. Sala, and G. Pellacani, Reflectance confocal microscopy for in vivo skin imaging, Photochem. Photobiol. 84, 1421 (2008)

    work page 2008

  24. [24]

    S. Kang, S. Jeong, W. Choi, H. Ko, T. D. Yang, J. H. Joo, J.-S. Lee, Y.-S. Lim, Q.-H. Park, and W. Choi, Imaging deep within a scattering medium using collective accumulation of single-scattered waves, Nat. Photon. 9, 253 (2015)

    work page 2015

  25. [25]

    C. W. Hsu, A. Goetschy, Y. Bromberg, A. D. Stone, and H. Cao, Broadband coherent enhancement of transmission and absorption in disordered media, Phys. Rev. Lett. 115, 223901 (2015)

    work page 2015

  26. [26]

    S. G. Adie, B. W. Graf, A. Ahmad, P. S. Carney, and S. A. Boppart, Computational adaptive optics for broadband optical interferometric tomography of biological tissue, Proc. Natl. Acad. Sci. U.S.A. 109, 7175 (2012)

    work page 2012

  27. [27]

    N. D. Shemonski, F. A. South, Y.-Z. Liu, S. G. Adie, P. S. Carney, and S. A. Boppart, Computational high-resolution optical imaging of the living human retina, Nat. Photon. 9, 440–443 (2015)

    work page 2015

  28. [28]

    F. A. South, Y.-Z. Liu, P.-C. Huang, T. Kohlfarber, and S. A. Boppart, Local wavefront mapping in tissue using compu- tational adaptive optics oct, Opt. Lett. 44, 1186 (2019)

    work page 2019

  29. [29]

    S. Liu, F. Xia, X. Yang, M. Wu, L. A. Bizimana, C. Xu, and S. G. Adie, Closed-loop wavefront sensing and correction in the mouse brain with computed optical coherence microscopy, Biomed. Opt. Express 12, 4934 (2021)

    work page 2021

  30. [30]

    Kumar, W

    A. Kumar, W. Drexler, and R. A. Leitgeb, Subaperture correlation based digital adaptive optics for full field optical coherence tomography, Opt. Express 21, 10850 (2013)

    work page 2013

  31. [31]

    Hillmann, H

    D. Hillmann, H. Spahr, C. Hain, H. Sudkamp, G. Franke, C. Pf¨ affle, C. Winter, and G. H¨ uttmann, Aberration-free volumetric high-speed imaging of in vivo retina, Sci. Rep. 6, 35209 (2016). 37

    work page 2016

  32. [32]

    Lee, D.-Y

    Y.-R. Lee, D.-Y. Kim, Y. Jo, M. Kim, and W. Choi, Exploiting volumetric wave correlation for enhanced depth imaging in scattering medium, Nat. Commun. 14, 1878 (2023)

    work page 2023

  33. [33]

    S. Kang, P. Kang, S. Jeong, Y. Kwon, T. D. Yang, J. H. Hong, M. Kim, K.-D. Song, J. H. Park, J. H. Lee, et al. , High- resolution adaptive optical imaging within thick scattering media using closed-loop accumulation of single scattering, Nat. Commun. 8, 2157 (2017)

    work page 2017

  34. [34]

    S. Yoon, H. Lee, J. H. Hong, Y.-S. Lim, and W. Choi, Laser scanning reflection-matrix microscopy for aberration-free imaging through intact mouse skull, Nat. Commun. 11, 5721 (2020)

    work page 2020

  35. [35]

    Z. Wang, Y. Zhang, and C. W. Hsu, Full-wave simulations of tomographic optical imaging inside scattering media, preprint at https://arxiv.org/abs/2308.07244 (2023)

    work page arXiv 2023

  36. [36]

    Titanium dioxide (TiO2) Nanoparticles XRD Analyses: An Insight

    T. Theivasanthi and M. Alagar, Titanium dioxide (TiO2) nanoparticles XRD analyses: an insight, arXiv:1307.1091 (2013)

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  37. [37]

    Mohanty, L

    S. Mohanty, L. B. Larsen, J. Trifol, P. Szabo, H. V. R. Burri, C. Canali, M. Dufva, J. Emn´ eus, and A. Wolff, Fabrication of scalable and structured tissue engineering scaffolds using water dissolvable sacrificial 3 D printed moulds, Mater. Sci. Eng. C 55, 569 (2015)

    work page 2015

  38. [38]

    M¨ atzler,MATLAB functions for Mie scattering and absorption, Institute of Applied Physics Research Report No

    C. M¨ atzler,MATLAB functions for Mie scattering and absorption, Institute of Applied Physics Research Report No. 2002-08, University of Bern (2002)

    work page 2002

  39. [39]

    J. R. DeVore, Refractive indices of rutile and sphalerite, J. Opt. Soc. Am. 41, 416 (1951)

    work page 1951

  40. [40]

    L. F. Rojas-Ochoa, J. M. Mendez-Alcaraz, J. J. S´ aenz, P. Schurtenberger, and F. Scheffold, Photonic properties of strongly correlated colloidal liquids, Phys. Rev. Lett. 93, 073903 (2004)

    work page 2004

  41. [41]

    Carminati and J

    R. Carminati and J. C. Schotland, Principles of Scattering and Transport of Light (Cambridge University Press, 2021)

    work page 2021

  42. [42]

    Duelk and K

    M. Duelk and K. Hsu, SLEDs and swept source laser technology for OCT, in Optical Coherence Tomography: Technology and Applications, edited by W. Drexler and J. G. Fujimoto (Springer, 2015) pp. 527–561

    work page 2015

  43. [43]

    Tinevez, N

    J.-Y. Tinevez, N. Perry, J. Schindelin, G. M. Hoopes, G. D. Reynolds, E. Laplantine, S. Y. Bednarek, S. L. Shorte, and K. W. Eliceiri, Trackmate: An open and extensible platform for single-particle tracking, Methods 115, 80 (2017)

    work page 2017

  44. [44]

    Ershov, M.-S

    D. Ershov, M.-S. Phan, J. W. Pylv¨ an¨ ainen, S. U. Rigaud, L. Le Blanc, A. Charles-Orszag, J. R. W. Conway, R. F. Laine, N. H. Roy, D. Bonazzi, G. Dum´ enil, G. Jacquemet, and J.-Y. Tinevez, Trackmate 7: integrating state-of-the-art segmentation algorithms into tracking pipelines, Nature Methods 19, 829 (2022)

    work page 2022

  45. [45]

    Fiji: an open-source platform for biological-image analysis, Nature Methods 9, 676 (2012). 38

    work page 2012