数据插值对正电子发射断层成像设备的图像重建影响的研究

正电子发射断层扫描(positron emission computed tomography, PET)是核医学领域最先进的临床检查影像技术. PET技术是目前临床上用于诊断和指导治疗肿瘤的最佳手段之一. 正电子发射断层成像设备探测器采集到的数据需要进行数据处理, 把原始数据转换成正弦图形式的数据才能用于图像重建. 平行束断层重建和扇形束图像重建是图像重建的两种方法, 分别对应平行束和扇形束形式的数据处理方法. 对原始数据的操作不可避免地破坏了原始数据的完整性. 现今, 正电子发射断层设备在重建过程中普遍采用平行束重建的方法. 平行束的数据分离会对PET数据进行插值操作, 扇形束的数据分离不会对PET数据进行插值操作. 本文通过对比平行束图像重建和扇形束图像重建结果, 研究了数据插值对PET图像重建结果的影响. 关键词： 正电子发射计算机断层扫描 数据插值 图像重建 原始数据     

An image reconstruction algorithm based on Bayesian theorem for electrical resistance tomography

Electrical resistance tomography (ERT) is considered to be one of the most promising process tomography techniques for two-phase/multiphase flow measurement due to the advantages such as high speed, low cost and non-intrusive sensing. The iterative image reconstruction algorithm based on Bayesian theorem has been derived to solve the ERT inverse problem. It is taken into account, which is different to the existing ERT image reconstruction algorithms, including the prior probability of the permittivity distribution and the noise information in the measurement data. Both simulation and experimental data were carried out for typical two-phase flow regimes. The flow regimes are identified according to the reconstructed images of the Bayesian iteration method and conventional methods. Results obtained indicate that the Bayesian iteration method improves the reconstructed image quality with the traditional linear back projection (LBP) and reweighted least squares (RLS). Therefore, the Bayesian iteration method is suitable for identification of online two-phase flow regimes.     

Dictionary learning based sinogram inpainting for CT sparse reconstruction

In CT (computed tomography), reconstruction from undersampling projection data is often ill-posed and suffers from severe artifact in the reconstructed images. To overcome this problem, this paper proposes a sinogram inpainting method based on recently rising sparse representation technology. In this approach, a dictionary learning based inpainting is used to estimate the missing projection data. The final image can be reconstructed by the analytic filtered back projection (FBP) reconstruction. We conduct experiments using both simulated and real phantom data. Compared to the comparative interpolation method, visual and numerical results validate the clinical potential of the proposed method.     

A modified discrete algebraic reconstruction technique for multiple grey image reconstruction for limited angle range tomography

The `missing wedge', which is due to a restricted rotation range, is a major challenge for quantitative analysis of an object using tomography. With prior knowledge of the grey levels, the discrete algebraic reconstruction technique (DART) is able to reconstruct objects accurately with projections in a limited angle range. However, the quality of the reconstructions declines as the number of grey levels increases. In this paper, a modified DART (MDART) was proposed, in which each independent region of homogeneous material was chosen as a research object, instead of the grey values. The grey values of each discrete region were estimated according to the solution of the linear projection equations. The iterative process of boundary pixels updating and correcting the grey values of each region was executed alternately. Simulation experiments of binary phantoms as well as multiple grey phantoms show that MDART is capable of achieving high‐quality reconstructions with projections in a limited angle range. The interesting advancement of MDART is that neither prior knowledge of the grey values nor the number of grey levels is necessary.     

基于最小区域的锥束ART算法快速图像重建

针对锥束ART算法重建速度慢的问题,考虑到工业CT重建目标尺寸差异较大的特点,提出了一种基于最小重建区域的快速三维图像重建方法。由不同视角下的锥束投影构建最小区域包络,对最小区域包络进行膨胀处理,以消除重建目标边界附近所出现的伪影,在此基础上,给出了射线穿过最小重建区域的体素索引及权因子的计算。该方法能够根据重建目标的尺寸自适应地确定最小重建区域,从而最大限度地减少了计算量。实验结果表明,该方法不仅大幅度提高了锥束ART算法的重建速度,而且有效地提高了重建质量。     

An algorithm for computed tomography image reconstruction from limited-view projections

<正>With the development of the compressive sensing theory,the image reconstruction from the projections viewed in limited angles is one of the hot problems in the research of computed tomography technology.This paper develops an iterative algorithm for image reconstruction,which can fit most cases.This method gives an image reconstruction flow with the difference image vector,which is based on the concept that the difference image vector between the reconstructed and the reference image is sparse enough.Then the l_2-norm minimization method is used to reconstruct the difference vector to recover the image for flat subjects in limited angles.The algorithm has been tested with a thin planar phantom and a real object in limited-view projection data.Moreover,all the studies showed the satisfactory results in accuracy at a rather high reconstruction speed.     

PITRE: software for phase‐sensitive X‐ray image processing and tomography reconstruction

Synchrotron‐radiation computed tomography has been applied in many research fields. Here, PITRE (Phase‐sensitive X‐ray Image processing and Tomography REconstruction) and PITRE_BM (PITRE Batch Manager) are presented. PITRE supports phase retrieval for propagation‐based phase‐contrast imaging/tomography (PPCI/PPCT), extracts apparent absorption, refractive and scattering information of diffraction enhanced imaging (DEI), and allows parallel‐beam tomography reconstruction for conventional absorption CT data and for PPCT phase retrieved and DEI‐CT extracted information. PITRE_BM is a batch processing manager for PITRE: it executes a series of tasks, created viaPITRE, without manual intervention. Both PITRE and PITRE_BM are coded in Interactive Data Language (IDL), and have a user‐friendly graphical user interface. They are freeware and can run on Microsoft Windows systems via IDL Virtual Machine, which can be downloaded for free and does not require a license. The data‐processing principle and some examples of application will be presented.     

层析成像图像重建算法综述

介绍了层析成像技术的图像重建算法,并从正向问题数学模型的简化和反向问题数学模型的映射结构的角度比较了各种算法的特点和优劣。研究表明:用本质是线性算法的各种变换方法重建图像存在严重失真,而卷积滤波的引入可以使变换方法的重建效果有所改善;基于导数搜索的迭代算法对初始值依赖性强、收敛速度慢并且容易陷入局部最优解;基于Fourier变换的方法具有本质的局限性;小波变换则可以同时刻画图像时域和频域的细节特征;有限元法通过重建对象像素的智能划分可以简化正问题的复杂性;而具有物理背景的蒙特卡罗法、模拟退火法、遗传算法、粒子滤波法及神经网络法更适合于复杂且非线性的图像重建;智能化、仿生化、并行化以及各种算法的融合是层析成像图像重建算法的发展趋势。     

A multi-phase level set framework for source reconstruction in bioluminescence tomography

We propose a novel multi-phase level set algorithm for solving the inverse problem of bioluminescence tomography. The distribution of unknown interior source is considered as piecewise constant and represented by using multiple level set functions. The localization of interior bioluminescence source is implemented by tracing the evolution of level set function. An alternate search scheme is incorporated to ensure the global optimal of reconstruction. Both numerical and physical experiments are performed to evaluate the developed level set reconstruction method. Reconstruction results show that the proposed method can stably resolve the interior source of bioluminescence tomography.     

Optimization-based image reconstruction in x-ray computed tomography by sparsity exploitation of local continuity and nonlocal spatial self-similarity

The additional sparse prior of images has been the subject of much research in problems of sparse-view computed tomography(CT) reconstruction. A method employing the image gradient sparsity is often used to reduce the sampling rate and is shown to remove the unwanted artifacts while preserve sharp edges, but may cause blocky or patchy artifacts.To eliminate this drawback, we propose a novel sparsity exploitation-based model for CT image reconstruction. In the presented model, the sparse representation and sparsity exploitation of both gradient and nonlocal gradient are investigated.The new model is shown to offer the potential for better results by introducing a similarity prior information of the image structure. Then, an effective alternating direction minimization algorithm is developed to optimize the objective function with a robust convergence result. Qualitative and quantitative evaluations have been carried out both on the simulation and real data in terms of accuracy and resolution properties. The results indicate that the proposed method can be applied for achieving better image-quality potential with the theoretically expected detailed feature preservation.     

深度非线性度量学习在说话人确认中的应用

将非线性度量学习(Nonlinear Metric Learning,NML)应用于说话人确认,提出了一种基于深度独立子空间分析(Independent Subspace Analysis,ISA)网络的说话人确认方法。区别于传统的线性度量学习方法,该方法使用深度独立子空间分析网络来学习一种从说话人原始空间到优化子空间的非线性显式映射,并在此基础上计算两条语音之间的相似性,以获得更好的说话人确认性能.所提方法在NIST SRE 2008数据集上进行了评估。评估结果表明,所提算法的等错误率指标相比传统的基于余弦距离打分的i-vector算法、线性判别分析(Linear Discriminant Analysis,LDA)算法、概率线性判别分析(Probabilistic Linear Discriminant Analysis,PLDA)算法分别下降了11.02%,6.40%和4.579%。     

基于稀疏优化的计算机断层成像图像不完全角度重建综述

计算机断层成像(computed tomography,CT)技术在医学和工业无损检测中都具有非常广泛的应用,CT重建算法是其中的核心,而不完全角度重建问题则是实际应用中重建算法研究领域的一个热点和难点问题.近年来,随着稀疏优化理论与算法的飞速发展,基于稀疏优化的重建算法已经在不完全角度重建问题中得到了较广泛的应用,且表现出了良好的精度与速度性能.本文首先对稀疏优化的基本理论结论与常用算法进行了介绍;而后对稀疏优化理论在CT图像不完全角度重建中的应用进行归纳,分类介绍了其主要研究成果及稀疏优化所发挥的作用;最后对基于稀疏优化的不完全角度重建研究进行了展望.     

Fast parallel algorithm for three-dimensional distance-driven model in iterative computed tomography reconstruction

The projection matrix model is used to describe the physical relationship between reconstructed object and projection.Such a model has a strong influence on projection and backprojection,two vital operations in iterative computed tomographic reconstruction.The distance-driven model(DDM) is a state-of-the-art technology that simulates forward and back projections.This model has a low computational complexity and a relatively high spatial resolution;however,it includes only a few methods in a parallel operation with a matched model scheme.This study introduces a fast and parallelizable algorithm to improve the traditional DDM for computing the parallel projection and backprojection operations.Our proposed model has been implemented on a GPU(graphic processing unit) platform and has achieved satisfactory computational efficiency with no approximation.The runtime for the projection and backprojection operations with our model is approximately 4.5 s and 10.5 s per loop,respectively,with an image size of 256×256×256 and 360 projections with a size of 512×512.We compare several general algorithms that have been proposed for maximizing GPU efficiency by using the unmatched projection/backprojection models in a parallel computation.The imaging resolution is not sacrificed and remains accurate during computed tomographic reconstruction.     

Depth-adaptive regularized reconstruction for reflection diffuse optical tomography

We study reflection diffuse optical tomography using two-dimensional (2D) continuous-wave source-detector arrays on the surface of semi-infinite medium, aiming at imaging the perfusion and the hemoglobin oxygen saturation variation of human cerebral cortex with brain activation. We had previously formulated the inverse problem with Moore-Penrose inversion. When we use simple regularization in this inverse problem, the reconstruction sensitivity decreases markedly with the depth so that the signal in the deep range may be masked by an unwanted signal in the shallow range. In this paper, we propose a depth-adaptive regularized reconstruction, in which we assign a smaller regularization parameter with the depth. We demonstrate improvement of the three-dimensional (3D) reconstruction uniformity using the proposed scheme.     

Deep learning aided beam vector assignment for massive MIMO maritime communication considering location information and handover impact

Location aided beamforming (LAB) has been proposed as a potential solution to fast beam vector (BV) selection for massive multi-input multi-output (M-MIMO) maritime communication systems. However, the existing LAB schemes hardly consider the impact from handover. In this paper, we propose to apply the Gaussian anomaly detection (GAD) to improve the system performance, particularly in the handover zone, by exploiting the information of sea lane and locations of the user equipment (UE). Furthermore, inspired by the rationale of deep learning (DL) aided image processing methods, we first transform the BV assignment problem to a location-based artificial gray value classification problem, and then utilize a convolutional deep belief networks (CDBN) to predict the BV indices from a predefined codebook constituted by location-related feature matrices. In addition, we exploit GAD to rectify the training data of CDBN, such that the prediction on BV assignments can be more accurate to achieve a higher capacity than known methods. Finally, an extensive comparison with selected existing schemes through simulations demonstrates the effectiveness of the proposed scheme.     

Deep learning method for testing the cosmic distance duality relation

The cosmic distance duality relation (DDR) is constrained by a combination of type-Ia supernovae (SNe Ia) and strong gravitational lensing (SGL) systems using the deep learning method. To make use of the full SGL data, we reconstruct the luminosity distance from SNe Ia up to the highest redshift of SGL using deep learning, and then, this luminosity distance is compared with the angular diameter distance obtained from SGL. Considering the influence of the lens mass profile, we constrain the possible violation of the DDR in three lens mass models. The results show that, in the singular isothermal sphere and extended power-law models, the DDR is violated at a high confidence level, with the violation parameter \begin{document}$ \eta_0=-0.193^{+0.021}_{-0.019} $\end{document} and \begin{document}$ \eta_0=-0.247^{+0.014}_{-0.013} $\end{document}, respectively. In the power-law model, however, the DDR is verified within a 1σ confidence level, with the violation parameter \begin{document}$ \eta_0=-0.014^{+0.053}_{-0.045} $\end{document}. Our results demonstrate that the constraints on the DDR strongly depend on the lens mass models. Given a specific lens mass model, the DDR can be constrained at a precision of \begin{document}$O(10^{-2}) $\end{document} using deep learning.     

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

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

Tomographic reconstruction using prior with adaptively binary neighborhood

Bayesian approaches can be used to improve ill-posed image reconstruction by regularizing the inverse solution using spacial or temporal neighborhood information. This paper proposes a prior with adaptively binary neighborhood (ABN prior) for statistical tomographic reconstruction. With binary weight map adapted to local image features, the proposed prior can lead to improved reconstruction by including relevant pixels belonging to similar structures and excluding those not. A two-step algorithm is also put forward for tomographic reconstruction using the proposed prior. Experiments using both simulated and clinical computerized tomography (CT) data are performed to validate the reconstructions with the proposed ABN prior.     

小波分析在层析图像重构中的应用研究

小波分析作为一种非平稳信号分析方法,具有良好的时（ 空）、频局部化特性和多分辨率特性。介绍了小波分析的基本原理和应用,引入小波分析进行图像重构,利用小波分解后得到的多分辨率的稀疏矩阵表示,设计了一迭代重构算法。通过计算机仿真试验,验证了小波分析在计算层析（ＣＴ） 成像光谱技术中能够应用于图像重构,并证实了小波分析的迭代重构算法是稳定、多分辨率和快速收敛的。     

结合透射CT的康普顿背散射图像重建

在康普顿散射成像(CST)技术中可以结合透射成像重建出衰减系数来消除散射重建的非线性,但这样得到的投影矩阵带有误差。而CST重建问题的不适定性对噪声和投影矩阵的误差非常敏感,重建结果会有较大误差。针对此问题,基于压缩感知理论提出了一种新的CST重建算法。新方法将图像重建问题归结为一个图像的全变分(TV)最小化问题,并使用收敛速度较快的基于交替方向法的Split-Bregman方法进行求解。在仿真实验中,通过与代数重建技术(ART)进行比较,在测量数据充足和测量数据不足两种情况下,本文算法都具有更好的重建质量,证明了所提算法在重建精度和抗噪性能方面的优势。