A modified equally sloped algorithm based on the total variation algorithm in computed tomography for insufficient data

Computed tomography (CT) has become an important technique for analyzing the inner structures of material, biological and energy fields. However, there are often challenges in the practical application of CT due to insufficient data. For example, the maximum rotation angle of the sample stage is limited by sample space or image reconstruction from the limited number of views required to reduce the X‐ray dose delivered to the sample. Therefore, it is difficult to acquire CT images with complete data. In this work, an iterative reconstruction algorithm based on the minimization of the image total variation (TV) has been utilized to develop equally sloped tomography (EST), and the reconstruction was carried out from limited‐angle, few‐view and noisy data. A synchrotron CT experiment on hydroxyapatite was also carried out to demonstrate the ability of the TV‐EST algorithm. The results indicated that the new TV‐EST algorithm was capable of achieving high‐quality reconstructions from projections with insufficient data.     

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

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

Parallel variable-density spiral imaging using nonlocal total variation reconstruction

The relatively long scan time is still a bottleneck for both clinical applications and research of magnetic resonance imaging. To reduce the data acquisition time, we propose a novel fast magnetic resonance imaging method based on parallel variable-density spiral acquisition, which combines undersampling optimization and nonlocal total variation reconstruction.The undersampling optimization promotes the incoherence of resultant aliasing artifact via the "worst-case" residual error metric, and thus accelerates the data acquisition. Moreover, nonlocal total variation reconstruction is utilized to remove such an incoherent aliasing artifact and so improve image quality. The feasibility of the proposed method is demonstrated by both numerical phantom simulation and in vivo experiment. The experimental results show that the proposed method can achieve high acceleration factor and effectively remove an aliasing artifact from data undersampling with well-preserved image details. The image quality is better than that achieved with the total variation method.     

基于压缩感知与扩展小波树的自适应压缩成像

为了在现有的采样条件下,通过新的压缩采样方式获得计算量小且质量更好的图像,提出了基于压缩感知与扩展小波树的自适应压缩成像方法。首先将图像投影到分区控制的DMD上,获得图像在低分辨率下的测量值,并通过压缩感知重构算法重构出低分辨图像,接着利用扩展小波树预测重要小波位置,通过DMD在小波域采样获取图像的细节信息,最后由小波逆变换恢复高分辨率图像。将该方法与最小化全变分算法（TVAL3）和近来提出的基于扩展小波树的自适应成像算法（EWT-ACS）效果进行对比,实验结果表明,以boat图像为例,在压缩感知采样率为0.75,整体采样率为10%的无噪声条件下,该方法相较于TVAL3、EWT-ACS算法信噪比提高了4.63 dB和2.87 dB,在附加噪声条件下成像效果也较好。该方法能极大地降低压缩感知重建算法的运行时间,同时减少采样次数,具有较好的抗噪性。     

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.     

Fast and accurate X‐ray fluorescence computed tomography imaging with the ordered‐subsets expectation maximization algorithm

The ordered‐subsets expectation maximization algorithm (OSEM) is introduced to X‐ray fluorescence computed tomography (XFCT) and studied; here, simulations and experimental results are presented. The simulation results indicate that OSEM is more accurate than the filtered back‐projection algorithm, and it can efficiently suppress the deterioration of image quality within a large range of angular sampling intervals. Experimental results of both an artificial phantom and cirrhotic liver show that with a satisfying image quality the angular sampling interval could be improved to save on the data‐acquisition time when OSEM is employed. In addition, with an optimum number of subsets, the image reconstruction time of OSEM could be reduced to about half of the time required for one subset. Accordingly, it can be concluded that OSEM is a potential method for fast and accurate XFCT imaging.     

有限角光学计算机层析重建中的投影采样研究

针对光学计算机层析术(Optical Computerized Tomography,OCT)中有限角条件下的严重非完全数据重建问题,提出了正交投影采样,结合基于改进的代数重建术(ART)的先验知识算法,以在尽可能少的投影方向数下较好地重建含遮挡物的三维流场。通过计算机模拟,详细讨论了在有限角条件下含遮挡物的三维流场的非完全数据重建精度及误差分析。结果表明,在有限角及含遮挡物条件下,采用正交投影采样可以极大地减少严重非完全数据重建中的误差,提高重建精度,从而为将光学计算机层析术应用于非完全数据的实测中提供了参考。     

A Compton scattering image reconstruction algorithm based on total variation minimization

Compton scattering imaging is a novel radiation imaging method using scattered photons.Its main characteristics are detectors that do not have to be on the opposite side of the source,so avoiding the rotation process.The reconstruction problem of Compton scattering imaging is the inverse problem to solve electron densities from nonlinear equations,which is ill-posed.This means the solution exhibits instability and sensitivity to noise or erroneous measurements.Using the theory for reconstruction of sparse images,a reconstruction algorithm based on total variation minimization is proposed.The reconstruction problem is described as an optimization problem with nonlinear data-consistency constraint.The simulated results show that the proposed algorithm could reduce reconstruction error and improve image quality,especially when there are not enough measurements.     

基于L1范数的微分相位衬度CT稀疏角度重建算法

微分相位衬度成像及其计算层析(CT)技术是近年出现的无损检测新方法。但是,相位衬度CT往往需要对样品进行多次扫描,这必将导致非常长的辐射时间和巨大的辐射剂量。稀疏角度重建在降低辐射剂量方面有着非常明显的优势,因此,研究针对相位衬度CT的稀疏角度重建算法就显得尤为重要。在分析了相位衬度CT的特点之后,将压缩感知理论引入相位衬度CT重建中,并在该理论框架下将L1约束融入代数迭代重建(ART)算法中,提出了一种微分相位衬度CT重建算法。数值模拟和实际实验表明,该方法可以根据少量投影数据给出较好的重建结果。     

Optimal tube potential and tube current for radiation dose reduction in routine adult head computed tomography scanning

ABSTRACT

We proposed a new, optimised scanning parameter that can reduce the patient’s radiation dose while maintaining image quality in a head computed tomography (CT) scan. We evaluated the clinical CT scan parameters (tube voltage, tube current, slice thickness pitch, scan range, rotation time, and CT dose index) for brain CT examination in a total of 52 multi-detector row spiral CTs (SOMATOM Definition AS+, Siemens Healthcare, Germany). The data were analysed, and the range of valid scan parameters was determined clinically using quartile distribution within the 95% confidence interval. The American Association of Physicists in Medicine performance evaluation phantom was used to acquire images using these scan parameters, and new, optimised CT scan parameters were proposed by analysing CT number accuracy, noise, uniformity, spatial resolution, and contrast resolution. The new CT scan parameters proposed were determined as tube voltage 100?kVp and tube current 300?mAs. Compared with conventional clinical scan parameters, tube voltage was reduced by 16.7% and tube current was decreased by 33.3%. Loss in imaging accuracy and uniformity of CT number was less than 20%, loss in noise was less than 40%, and no change in resolution was observed. Conversely, the CT dose index and effective dose was 20%–50%. A new systematic method for clinically assessing the optimised CT scan parameters were proposed, and the effective dose was decreased, with changing exposure conditions.     

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

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

Magnetic resonance image enhancement using highly sparse input

Lately, the Magnetic Resonance scans have struggled with its own inherent limitations, such as spatial resolution as well as long examination times. In this paper, a novel, rapid compressively-sensed magnetic resonance high resolution image resolution algorithm is presented. This technique addresses these two key issues by employing a highly-sparse sampling scheme and super-resolution reconstruction (SRR) method. Due to highly challenging requirements for the accuracy of diagnostic images registration, the presented technique exploits image priors, deblurring, parallel imaging, and a discrete dense displacement sampling for the deformable human body and motion analysis. The clinical trials as well as phantom based studied have been conducted. It has been proven that the proposed algorithm is able to enhance image spatial resolution, reduce motion artefacts and scan times.     

螺旋锥束计算机断层成像倾斜扇束反投影滤波局部重建算法

螺旋锥束计算机断层成像(CT)作为常用的临床诊断工具,如何尽可能地减少其辐射剂量是热点研究领域之一.局部成像利用准直器减小射线直照区域,能够有效降低CT辐射剂量.然而,局部成像会造成投影数据横向截断,产生局部重建问题.现有螺旋反投影滤波(BPF)算法只能实现局部曲面重建,难以实现局部体区域重建.在圆轨迹扇束BPF算法的基础上,通过加权修正和坐标扩展,提出了螺旋锥束CT倾斜扇束反投影滤波(TFB-BPF)重建算法.该算法把重建区域按层划分,对各层构建倾斜的扇形束几何,并使用经过加权修正的TFB-BPF算法逐层进行重建.该算法最大的优势是滤波线(与原始螺旋锥束BPF算法中PI线等价)在二维平面内选择,算法更加简洁高效,并且能够应用于局部体区域的重建.实验结果表明,算法能够有效重建物体局部体区域,并且重建图像质量较好,没有明显的截断伪影.     

光学投影层析的快速图像重建和影响因素分析

介绍了光学投影层析(Optical Projection Tomography,OPT)图像重建中的滤波反投影(Filtered Back-projection,FBP)算法,针对该算法中运算量最大的射束运算,提出了利用像素点坐标对称性来减少运算量的FBP快速实现方法,可使重建时间缩短约15%.并对影响图像重建质量的投影角数目、旋转台步进精度、CCD中心像素偏离等因素进行了分析,得出的结论对OPT系统设计具有重要参考价值.     

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

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

基于全变分最小化和交替方向法的康普顿散射成像重建算法

康普顿散射成像技术利用射线与物质作用后的散射光子信息对物质的电子密度进行成像.与传统的透射成像方式相比,康普顿散射成像具有系统结构灵活、成像对比度高、辐射剂量低等优势,在无损检测、医疗诊断、安全检查等领域有着广阔的应用前景.但其重建问题是一个非线性的逆问题,通常是不适定的,其解对噪声和测量误差非常敏感.为解决此问题,本文结合全变分最小化正则化方法和交替方向法提出了一种新的康普顿散射成像重建算法.该算法首先将问题对应的TV模型转化为与之等价的带约束的优化问题,然后利用增广拉格朗日乘子法将优化问题分解为两个具有解析解的子问题,并通过交替求解子问题使增广拉格朗日函数达到最小,进而得到重建的图像.在仿真实验中,通过与主流的ASD-POCS方法进行对比,证明了该算法在重建精度和重建效率方面的优势.     

Total variation based gradient descent algorithm for sparse-view photoacoustic image reconstruction

In photoacoustic imaging (PAI), reconstruction from sparse-view sampling data is a remaining challenge in the cases of fast or real-time imaging. In this paper, we present our study on a total variation based gradient descent (TV-GD) algorithm for sparse-view PAI reconstruction. This algorithm involves the total variation (TV) method in compressed sensing (CS) theory. The objective function of the algorithm is modified by adding the TV value of the reconstructed image. With this modification, the reconstructed image could be closer to the real optical energy distribution map. Additionally in the proposed algorithm, the photoacoustic data is processed and the image is updated individually at each detection point. In this way, the calculation with large matrix can be avoided and a more frequent image update can be obtained. Through the numerical simulations, the proposed algorithm is verified and compared with other reconstruction algorithms which have been widely used in PAI. The peak signal-to-noise ratio (PSNR) of the image reconstructed by this algorithm is higher than those by the other algorithms. Additionally, the convergence of the algorithm, the robustness to noise and the tunable parameter are further discussed. The TV-based algorithm is also implemented in the in vitro experiment. The better performance of the proposed method is revealed in the experiments results. From the results, it is seen that the TV-GD algorithm may be a practical and efficient algorithm for sparse-view PAI reconstruction.     

基于同伦l0范数最小化重建的三维动态磁共振成像

高欠采倍数的动态磁共振图像重建具有重要意义,是同时实现高时间分辨率和高空间分辨率动态对比度增强成像的重要环节.本研究提出一种结合黄金角变密度螺旋采样、并行成像和基于同伦l₀范数最小化的压缩感知的图像重建的三维动态磁共振成像方法.黄金角变密度螺旋采样轨迹被用来连续获取k空间数据,具有数据采集效率高、对运动不敏感等优点.在重建算法中,将多线圈稀疏约束应用于时间总变分域,使用基于l₀范数最小化的非线性重建算法代替传统的l₁范数最小化算法,进一步提高了欠采样率.仿真实验和在体实验表明本文所提的方法在保持图像质量的同时,也可以实现较高的空间分辨率和时间分辨率,初步验证了基于同伦l₀范数最小化重建在三维动态磁共振成像上的优势和临床价值.     

压缩感知同步扫描重建及其采样方案的研究

压缩感知(compressed sensing,CS)-磁共振成像(magnetic resonance imaging,MRI)技术使用随机欠采样的k空间数据来重建图像,大大提高了成像速度.但典型的CS重建很费时,这也是CS-MRI临床应用的主要障碍之一.针对这一问题,该文提出了在扫描时同步进行CS图像重建的方案.在同步重建的过程中,可以实时显示重建图像的结果,用户可以根据图像质量来决定何时终止扫描,这样可以在节约扫描和重建时间的同时,更好地控制图像质量.由于预先无法确定最终的采样率,因此传统的变密度随机采样方法并不完全适用.该文设计了适用于同步重建过程的采样模式生成方案,同时提出了分段采样方法,把采样过程分为两个阶段,不同阶段使用不同的概率密度函数(probability density function,PDF)确定待采样的相位编码行.模拟实验的结果表明,与使用单一密度函数的采样方案相比,分段采样方案能够在整个同步扫描重建过程中始终获得更好的图像.     

Compressed sensing MR image reconstruction using a motion-compensated reference

Compressed sensing (CS)-based methods have been proposed for image reconstruction from undersampled magnetic resonance data. Recently, CS-based schemes using reference images have also been proposed to further reduce the sampling requirement. In this study, we propose a new reference-constrained CS reconstruction method that accounts for the misalignment between the reference and the target image to be reconstructed. The proposed method uses a new image model that represents the target image as a linear combination of a motion-dependent reference image and a sparse difference image. We then use an efficient iterative algorithm to jointly estimate the motion parameters and the difference image from sparsely sampled data. Simulation results from a numerical phantom data set and an in vivo data set show that the proposed method can accurately compensate the motion effects between the reference and the target images and improve reconstruction quality. The proposed method should prove useful for several applications such as interventional imaging, longitudinal imaging studies and dynamic contrast-enhanced imaging.