基于布雷格曼Split Bregman去噪的荧光漫反射迭代重建光学断层扫描（Matlab代码实现）

 💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、数据、文献

---

💥1 概述

文献来源：

摘要：

荧光漫反射光学断层扫描（fDOT）是一种非侵入性成像技术，可以量化小动物中荧光示踪剂的空间分布。fDOT图像重建通常通过代数重建技术（ART）等迭代方法进行。更高级产生的有用结果l1-用于信号恢复和图像重建的正则化技术，以及最近发表的 Split Bregman （SB） ，使我们提出了一种解决 fDOT 逆问题的新方法，即 ART-SB。该方法将具有成本效益的重建步骤（ART迭代）与基于最小化图像总变化的去噪滤波步骤交替使用SB方法，可以高效快速地求解。我们将这种方法应用于模拟和实验的 fDOT 数据，发现 ART-SB 比传统 ART 具有显着优势。

最近发表的 Split Bregman （SB） 方法19是一种简单高效的求解算法l11基于正则化的问题，可以拆分l11和l22泛 函。在参考文献 19 中，通过将 SB 方法应用于图像去噪和压缩传感，作者表明 SB 在计算上是高效的，因为 SB 公式导致了一个可以使用 Gauss-Seidel 和傅里叶变换方法解决的问题。

SB最近被应用于荧光断层扫描重建。20,21Abascal 等人。20使用 SB 通过施加非负性约束来解决优化问题。使用非线性高斯-牛顿步长更新图像6基于非线性电视泛函的一阶导数和二阶导数的计算。Behrooz 等人。21比较l22基于ROF和SB的正则化方法和带有TV重建方法的ART。作者实现了一种预处理共轭梯度方法6在 SB 的每次迭代中，在某些情况下会导致收敛缓慢。为了验证该方法并比较重建，他们使用了非接触式恒定波透射荧光断层扫描系统，并得出结论，电视正则化有可能提供比传统更高的分辨率和鲁棒性l22正则化算法和 ART。

我们的工作目标是提出一种解决fDOT逆问题的新方法，即ART-SB。基于两步迭代过程，我们将计算效率高的重建方法ART与去噪步骤相结合。去噪步骤基于 SB 公式，并使用 Gauss-Seidel 和收缩运算有效地实现，而无需计算 TV 泛函的一阶导数和二阶导数。

详细文章讲解见第4部分。

📚2 运行结果

部分代码：

for it = 1:numIter
    % ART reconstruction step: Iterative linear solver
    sol = ARTReconstruction_Fast(JacMatrix,data,relaxParam,numIterART,uARTSB(:)); 
    %sol = ARTReconstruction(JacMatrix,data,relaxParam,numIterART,uARTSB(:)); 
    solGrid     = reshape(sol,N);   

    % SB denoising step
    uARTSB      = TV_SB_denoising_3D(solGrid,mu,lambda,alpha,nInner,nOuter);    
    % Uncomment below to do 2D slice-by-slice smoothing instead. It takes
    % similar time but it could be parallelized, which can be faster in some
    % applications for large scale problems (now SB denoising takes less
    % than a second), by changing the for to a parfor loop  
%     for iz = 1:N(3)
%         uARTSB(:,:,iz) = TV_SB_denoising_2D(solGrid(:,:,iz),mu,lambda,nInner,nOuter);
%     end

    % Compute solution error norm
    err(it) = norm(uARTSB(:)-uTarget(:))/norm(uTarget(:));
    
    waitbar(it/numIter);
end % it
toc
close(h);

% Display results
figure; plot(err); ylabel('Solution error'); xlabel('Number of iterations');
title('Convergence of ART-SB');

% Target image
Plot2DMapsGridSolution(reshape(uTarget,N),X,Y,Z,3); 
set(gcf,'name','TARGET','numbertitle','off'); 
colormap gray;

% Reconstructed images0
Plot2DMapsGridSolution(uART,X,Y,Z,3); 
set(gcf,'name','ART reconstruction','numbertitle','off') 
colormap gray;

Plot2DMapsGridSolution(uARTSB,X,Y,Z,3); 
set(gcf,'name','ART-SB reconstruction','numbertitle','off') 
colormap gray;  

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

🌈4 Matlab代码、数据、文献