Quasi-SVD把矩阵分解速度提了3-20倍,还能在GPU上并行跑,医学实时成像终于有戏了。
传统SVD算法串行,制约GPU吞吐量,难以满足临床实时需求。Quasi-SVD是一种可微分、完全并行的矩阵分解框架,通过单个李参数化因子保证正交性,软约束恢复其余部分。在SSIM=0.89-0.94下,相比cuSOLVER和随机SVD加速3-20倍,吞吐量超过25 FPS。在超声定位显微镜和Mueller矩阵偏振测量两个医学成像任务上验证了鲁棒性,使结构化矩阵分解适用于实时图像引导工作流。
Quasi-SVD: Learning a Lie-constrained matrix factorisation for real-time imaging
Singular Value Decomposition (SVD) underlies matrix factorisation tasks across computational imaging, with medical applications increasingly demanding real-time processing. Yet SVD algorithms are inherently sequential, constraining real-time GPU throughput and limit online deployment in clinical pipelines. This study introduces Quasi-SVD, a differentiable, fully parallelized matrix factorization framework for GPUs. Rather than enforcing orthogonality on both factors, it guarantees exact orthogonality for a single Lie-parameterized factor while recovering the remaining components through soft constraints, enabling efficient parallel decomposition without iterative singular-vector optimization. This asymmetric design, provably sufficient for valid factorisation, achieves reconstruction fidelity of SSIM = 0.89-0.94 and accelerates computation by 3-20x relative to cuSOLVER and randomised SVD, enabling throughput above 25 FPS. Performance is evaluated on two medical imaging tasks spanning complementary computational regimes: (1) spatio-temporal background subtraction for ultrasound localisation microscopy, requiring high-dimensional matrix separation, and (2) Mueller matrix polarimetry for neurosurgical tissue characterisation, requiring massive batch processing of small matrices. Across both regimes and multiple imaging instruments, the proposed framework demonstrates robust domain transfer and throughput exceeding 25 FPS at clinical matrix scales, a rate sufficient for live image-guided workflows that classical solvers cannot currently support in these settings. By prioritising downstream reconstruction fidelity over exact spectral recovery, Quasi-SVD makes structured matrix factorisation practical for real-time imaging.