High total variation-based method for sparse-view photoacoustic reconstruction

We propose a novel method by combining the total variation（TV） with the high-degree TV（HDTV） to improve the reconstruction quality of sparse-view sampling photoacoustic imaging（PAI）. A weighing function is adaptively updated in an iterative way to combine the solutions of the TV and HDTV minimizations. The fast iterative shrinkage/thresholding algorithm is implemented to solve both the TV and the HDTV minimizations with better convergence rate. Numerical results demonstrate the superiority and efficiency of the proposed method on sparse-view PAI. In vitro experiments also illustrate that the method can be used in practical sparse-view PAI.     

非精确交替方向总变分最小化重建算法

CT (computed tomography)系统实际应用当中, 经常会出现扫描数据不满足数据完备性条件的情况. 针对不完全角度重建问题的研究, 是目前迭代型算法研究中的一个热点. 一系列基于带有约束的总变分最小化的重建算法近年来在不完全角度重建中取得了较好的效果, 这其中基于交替方向法 (alternating direction method, ADM) 的重建算法表现出更好的性能. 然而, ADM方法在求解过程中对矩阵求逆的处理效率不高, 导致极大的计算开销. 本文针对该问题, 使用非精确ADM方法, 利用线性近似的方式替换掉计算开销较大的项, 使得矩阵求逆问题可以通过快速傅里叶变换加速实现. 实验结果表明, 本文提出的非精确交替方向总变分最小化重建算法与精确ADM重建算法相比, 没有明显的精度损失, 计算时间缩减30%左右. 关键词： 不完全角度重建 总变分最小化 非精确交替方向法     

Novel Fourier-based iterative reconstruction for sparse fan projection using alternating direction total variation minimization

Sparse-view x-ray computed tomography(CT) imaging is an interesting topic in CT field and can efficiently decrease radiation dose. Compared with spatial reconstruction, a Fourier-based algorithm has advantages in reconstruction speed and memory usage. A novel Fourier-based iterative reconstruction technique that utilizes non-uniform fast Fourier transform(NUFFT) is presented in this work along with advanced total variation(TV) regularization for a fan sparse-view CT. The proposition of a selective matrix contributes to improve reconstruction quality. The new method employs the NUFFT and its adjoin to iterate back and forth between the Fourier and image space. The performance of the proposed algorithm is demonstrated through a series of digital simulations and experimental phantom studies. Results of the proposed algorithm are compared with those of existing TV-regularized techniques based on compressed sensing method, as well as basic algebraic reconstruction technique. Compared with the existing TV-regularized techniques, the proposed Fourier-based technique significantly improves convergence rate and reduces memory allocation, respectively.     

总变差约束的数据分离最小图像重建模型及其Chambolle-Pock求解算法

基于优化的迭代法,可以结合压缩感知和低秩矩阵等稀疏优化技术高精度地重建图像.其中,总变差最小（total variation minimization,TV）模型是一种简单有效的优化模型.传统的约束TV模型,使用数据保真项为约束项,TV正则项为目标函数.本文研究TV约束的、数据分离最小（TV constrained,data divergence minimization,TVcDM）新型TV模型及其求解算法.详细推导了TVcDM模型的Chambolle-Pock（CP）算法,验证了模型及算法的正确性;分析了算法的收敛行为;评估了模型的稀疏重建能力;分析了模型参数的选择对重建的影响及算法参数对收敛速率的影响.研究表明,TVcDM模型有高精度稀疏重建能力;TVcDM-CP算法确保收敛,但迭代过程中有振荡现象;TV限对重建有重要影响,参数值过大会引入噪声而过小会模糊图像细节;算法参数的不同选取会导致不同的收敛速率.     

Fast multi-contrast MRI reconstruction

Multi-contrast magnetic resonance imaging (MRI) is a useful technique to aid clinical diagnosis. This paper proposes an efficient algorithm to jointly reconstruct multiple T1/T2-weighted images of the same anatomical cross section from partially sampled k-space data. The joint reconstruction problem is formulated as minimizing a linear combination of three terms, corresponding to a least squares data fitting, joint total variation (TV) and group wavelet-sparsity regularization. It is rooted in two observations: 1) the variance of image gradients should be similar for the same spatial position across multiple contrasts; 2) the wavelet coefficients of all images from the same anatomical cross section should have similar sparse modes. To efficiently solve this problem, we decompose it into joint TV regularization and group sparsity subproblems, respectively. Finally, the reconstructed image is obtained from the weighted average of solutions from the two subproblems, in an iterative framework. Experiments demonstrate the efficiency and effectiveness of the proposed method compared to existing multi-contrast MRI methods.     

基于成像差分吸收光谱技术和压缩感知理论的烟囱烟羽断层重建研究

烟羽断层重建质量受两方面条件限制：其中一个限制条件是遥感设备的时间分辨率。以往的研究多使用多轴差分吸收光谱仪（MAX-DOAS）进行CT重建,受采集数据速度的限制,重建图像的时间分辨率较低。另一个限制条件是,采集到的数据量有限,是典型的不完全角度重建。过去多使用代数迭代重建算法或统计迭代重建算法,重建图像受测量误差的影响比较大,分辨率较低且伪影较多。构造了基于成像差分吸收光谱技术（IDOAS）的光谱数据采集系统,与多轴差分吸收光谱仪构造的系统相比,数据采集的时间分辨率提高了160多倍,基本解决了时间分辨率的问题。提出了一种基于压缩感知理论和低三阶导数模型的烟羽断层重建算法--投影凸函数集低三阶导数法,简称为POCS-LTD。在投影的过程中,使用代数重建算法使重建图像符合投影方程;在全变分迭代的过程中使用了优化算法,将低三阶导数模型的全变分归一化值作为优化算法的迭代方向,前次迭代运算结果与本次投影运算的差值的模作为迭代步长。对重建算法进行了数值模拟,并以重建图像的接近度和一致性相关因子为指标,对重建结果进行了分析。数值模拟表明,算法具有良好的抗误差能力,与传统的低三阶导数法相比,本文提出的算法将重建接近度减小了80%以上。使用烟羽数据采集系统进行了外场实验,用POCS-LTD算法对外场实验的数据进行了烟羽重建,重建图像显示烟羽图像清晰,伪影得到了较好的抑制。介绍的烟羽断层数据采集系统和烟羽断层重建算法,提高了烟羽断层重建图像的时间分辨率,减少了重建图像的伪影,扩大了光谱测量技术的应用范围。     

基于lp范数的压缩感知图像重建算法研究

图像重建是光学成像、光声成像、声纳成像、核磁共振成像、 天体成像等物理成像领域中的关键技术之一. 近年来提出的压缩感知理论指出: 对稀疏或者可压缩信号进行少量非自适应线性投影,投影信号含有足够的信息, 从而能对信号进行高概率重建. 压缩感知已被应用于多种物理成像系统. 将罚函数法和修正Hesse阵序列二次规划方法相结合, 并采用了分块压缩感知思想, 提出一种基于l_p范数的压缩感知图像重建算法. 以cameraman, barbara和mandrill图像为例, 采用该算法进行图像重建. 首先, 在不同采样率下对图像重建. 即便采样率低至0.3时, 也能获得高达32.23dB的信噪比, 重建图像清晰可辨. 验证了该算法的正确性. 其次, 将该算法与正交匹配追踪算法进行对比, 在采样率达到0.5以上时, 能够获得高信噪比的重建图像, 成像时间也大为减少, 特别是采样率为0.7时, 成像时间减少88%. 最后, 与现有基于l_p 范数的压缩感知图像重建算法进行对比, 计算结果表明在成像质量有所提高的基础上, 成像时间大为缩短. 关键词： 图像重建 压缩感知 罚函数 修正Hesse阵序列二次规划     

An Efficient Proximity Point Algorithm for Total-Variation-Based Image Restoration

In this paper, we propose a fast proximity point algorithm and apply it to total variation (TV) based image restoration. The novel method is derived from the idea of establishing a general proximity point operator framework based on which new first-order schemes for total variation (TV) based image restoration have been proposed. Many current algorithms for TV-based image restoration, such as Chambolle's projection algorithm, the split Bregman algorithm, the Bermúdez-Moreno algorithm, the Jia-Zhao denoising algorithm, and the fixed point algorithm, can be viewed as special cases of the new first-order schemes. Moreover, the convergence of the new algorithm has been analyzed at length. Finally, we make comparisons with the split Bregman algorithm which is one of the best algorithms for solving TV-based image restoration at present. Numerical experiments illustrate the efficiency of the proposed algorithms.     

用于光声层析成像的虚拟阵列探测及逆拉东变换

提出用固定孔径的换能器虚拟阵列探测器进行数据采集,该换能器只需对样品进行180扫描。实验中,使用调Q Nd: YAG激光器作为光源,激光器输出激光波长为532 nm,重复频率为10 Hz,脉宽为7 ns;使用非聚焦换能器进行光声信号探测,其接收面直径为6 mm,中心频率为2.25 MHz,接收到光声信号用逆拉东变换进行图像重建。以小白鼠的离体爪子作为成像对象,换能器由旋转电机控制并以2作为旋转步长进行180旋转扫描获取实验数据,运用逆拉东变换进行了图像重建。所重建的图像比较清晰地反映了爪子内部的骨骼分布。由此表明,逆拉东变换可以用于光声成像,并且实现了有限角度的光声成像。     

Reconstruction method of computational ghost imaging based on non-local generalized total variation北大核心CSCD

鬼成像是一种能够透过大雾等恶劣环境的成像技术。针对传统鬼成像重建图像存在噪声较多、图像对比度较低等问题,将非局部广义全变分方法用于鬼成像的图像重建之中,提出基于非局部广义全变分的计算鬼成像重建方法。所提方法构造了一种非局部相关性权重设计梯度算子,将其代入全变分重建算法中,使得重建的图像能有效去除噪声的同时实现细节较好的还原。首先在不同条件下进行仿真模拟,得到所提方法的峰值信噪比相对其他方法提升1 dB左右,且具有更好的主观视觉效果,进而设计并搭建实验平台对算法的有效性进行验证,实验结果证明了所提方法在去除噪声和细节重建等方面的优越性。     

Anisotropic total variation minimization approach in in-line phase-contrast tomography and its application to correction of ring artifacts

In-line phase-contrast computed tomography(IL-PC-CT) imaging is a new physical and biochemical imaging method.IL-PC-CT has advantages compared to absorption CT when imaging soft tissues. In practical applications, ring artifacts which will reduce the image quality are commonly encountered in IL-PC-CT, and numerous correction methods exist to either pre-process the sinogram or post-process the reconstructed image. In this study, we develop an IL-PC-CT reconstruction method based on anisotropic total variation(TV) minimization. Using this method, the ring artifacts are corrected during the reconstruction process. This method is compared with two methods: a sinogram preprocessing correction technique based on wavelet-FFT filter and a reconstruction method based on isotropic TV. The correction results show that the proposed method can reduce visible ring artifacts while preserving the liver section details for real liver section synchrotron data.     

Photoacoustic image reconstruction based on Bayesian compressive sensing algorithm

The photoacoustic tomography (PAT) method, based on compressive sensing (CS) theory, requires that, for the CS reconstruction, the desired image should have a sparse representation in a known transform domain. However, the sparsity of photoacoustic signals is destroyed because noises always exist. Therefore, the original sparse signal cannot be effectively recovered using the general reconstruction algorithm. In this study, Bayesian compressive sensing (BCS) is employed to obtain highly sparse representations of photoacoustic images based on a set of noisy CS measurements. Results of simulation demonstrate that the BCS-reconstructed image can achieve superior performance than other state-of-the-art CS-reconstruction algorithms.     

双树小波变换与小波树稀疏联合的低场CS-MRI算法

压缩感知理论常用在磁共振快速成像上,仅采样少量的K空间数据即可重建出高质量的磁共振图像.压缩感知磁共振成像技术的原理是将磁共振图像重建问题建模成一个包含数据保真项、稀疏先验项和全变分项的线性组合最小化问题,显著减少磁共振扫描时间.稀疏表示是压缩感知理论的一个关键假设,重建结果很大程度上依赖于稀疏变换.本文将双树复小波变换和小波树稀疏联合作为压缩感知磁共振成像中的稀疏变换,提出了基于双树小波变换和小波树稀疏的压缩感知低场磁共振图像重建算法.实验表明,本文所提算法可以在某些磁共振图像客观评价指标中表现出一定的优势.     

Remote photoacoustic imaging on solid material using a two-wave mixing interferometer

We report on remote and contactless photoacoustic imaging (PAI) for the inspection of solid materials using a two-wave mixing interferometer. In this Letter, a semitransparent sample was excited with picosecond laser pulses. The local absorption of the electromagnetic radiation led to generation of broadband ultrasonic waves inside the sample. Ultrasonic waves arriving at the sample surface were detected utilizing a two-wave mixing interferometer. After data acquisition, the initial pressure distribution was reconstructed using a Fourier space synthetic aperture technique algorithm. We show the potential of PAI for the inspection of semitransparent solid materials.     

Singular Value Decomposition Using Jacobi Algorithm in pMRI and CS

Parallel magnetic resonance imaging (pMRI) and compressed sensing (CS) have been recently used to accelerate data acquisition process in MRI. Matrix inversion (for rectangular matrices) is required to reconstruct images from the acquired under-sampled data in various pMRI algorithms (e.g., SENSE, GRAPPA) and CS. Singular value decomposition (SVD) provides a mechanism to accurately estimate pseudo-inverse of a rectangular matrix. This work proposes the use of Jacobi SVD algorithm to reconstruct MR images from the acquired under-sampled data both in pMRI and in CS. The use of Jacobi SVD algorithm is proposed in advance MRI reconstruction algorithms, including SENSE, GRAPPA, and low-rank matrix estimation in L + S model for matrix inversion and estimation of singular values. Experiments are performed on 1.5T human head MRI data and 3T cardiac perfusion MRI data for different acceleration factors. The reconstructed images are analyzed using artifact power and central line profiles. The results show that the Jacobi SVD algorithm successfully reconstructs the images in SENSE, GRAPPA, and L + S algorithms. The benefit of using Jacobi SVD algorithm for MRI image reconstruction is its suitability for parallel computation on GPUs, which may be a great help in reducing the image reconstruction time.     

无透镜显微成像的重构算法研究

针对单幅图像进行了无透镜显微成像的重构算法研究，介绍了无透镜显微成像系统实验装置和ASM(angle spectrum method)、改编后的L-R(Lucy-Richardson)两种重构算法。对比两种算法重构后的USAF分辨率板图像的分辨率，利用瑞利判据得出ASM获得的振幅图分辨率最高(即3.10 μm)，且计算用时最少(即0.9 s)，证明了ASM为最佳的单幅无透镜显微重构算法。其次，利用无透镜显微成像系统结合ASM重构的方法，进行细胞成像实验。该无透镜成像视场为5×显微镜的4.4倍，且分辨率介于5×及10×光学显微镜之间，统计学优势明显，在生物医学领域具有广阔的应用前景。     

等离子体图像诊断中带吸收校正的轴对称重建

在轴对称等离子体的图像诊断中,根据等离子体在某个方向的X射线强度图像可以重建出辐射系数的空间分布.由于X射线在等离子体的传播过程中存在一定的吸收,需要在重建过程中考虑吸收校正.本文在滤波反投影法(FBP)的基础上,结合X射线在等离子体中的吸收规律,提出了带吸收校正的FBP迭代算法.通过对钠等离子体辐射系数的轴对称重建模拟,发现该算法重建精度高,收敛性好,取得了很好的重建结果.     

基于广义变分正则化的闪光照相图像重建算法

 针对闪光照相图像信噪比低的特点,提出了一种基于广义变分正则化的图像重建算法,该方法采用p-范数取代目前广泛采用的全变分范数作为正则项,构造了用于图像重建的展平泛函,将图像重建问题转化为目标泛函最优化问题,采用固定点迭代法求解图像重建的最优解。数值计算结果表明,该算法在重建过程中能够有效抑制图像噪声,并加大对图像边缘的保持能力,从而提高了图像重建质量,是一种有效且性能优良的闪光照相图像重建算法。     

面向近红外脑功能成像的光学自导引漫射光断层成像方法

为降低近红外脑功能漫射光断层成像(DOT)固有的逆问题病态性,并避免多模态方法的图像配准等问题,提出了基于光学自导引提供先验功能信息的脑功能DOT方法(OT-DOT),并发展了图像重构方法.模拟验证表明:上皮厚度(TLT)已知时,OT-DOT获得的重构量化度(QR)约为传统DOT的4.2倍;当TLT的估计误差小于±10%时,OT-DOT重构的QR值可达92%以上,远远优于传统DOT;噪声鲁棒性测试表明,OT-DOT与传统DOT的噪声鲁棒性相近.利用连续光DOT测量系统的仿体实验重构结果表明,所发展的OT-DOT算法获得的重构结果优于传统DOT算法.     

不同旋转扫描次数对光声层析成像的影响

为了实现单元探测器高质量的快速光声图像重建,提出了不同旋转扫描次数对光声层析成像的影响方案。实验采用的光源为YAG激光器,波长为1 064 nm,重复频率为20 Hz,脉宽为7 ns,探测器为针状的PVDF膜水听器,接收直径为1 mm,得到了26个字母、12根头发丝、树叶骨架和模拟血管的光声重建图像。由仿真和实验结果表明,在不牺牲光声重建图像质量的前提下,单元探测器环形扫描一圈,均匀采集100个位置的信号,图像重建时间为5.903 s。该研究结果对于单元探测器的快速旋转扫描成像和环形阵列探测器阵元数的设计具有一定的指导意义。