Wavelet-based reconstruction for limited-angle X-ray tomography

The aim of X-ray tomography is to reconstruct an unknown physical body from a collection of projection images. When the projection images are only available from a limited angle of view, the reconstruction problem is a severely ill-posed inverse problem. Statistical inversion allows stable solution of the limited-angle tomography problem by complementing the measurement data by a priori information. In this work, the unknown attenuation distribution inside the body is represented as a wavelet expansion, and a Besov space prior distribution together with positivity constraint is used. The wavelet expansion is thresholded before reconstruction to reduce the dimension of the computational problem. Feasibility of the method is demonstrated by numerical examples using in vitro data from mammography and dental radiology.     

Structure-from-motion without correspondence from tomographic projections by Bayesian inversion theory

In conventional tomography, the interior of an object is reconstructed from tomographic projections such as X-ray or transmission electron microscope images. All the current reconstruction methods assume that projection geometry of the imaging device is either known or solved in advance by using e.g., fiducial or nonfiducial feature points in the images. In this paper, we propose a novel approach where the imaging geometry is solved simultaneously with the volume reconstruction problem while no correspondence information is needed. Our approach is a direct application of Bayesian inversion theory and produces the maximum likelihood or maximum a posteriori estimates for the motion parameters under the selected noise and prior distributions. In this paper, the method is implemented for a two-dimensional model problem with one-dimensional affine projection data. The performance of the method is tested with simulated and measured X-ray projection data.     

Reconstruction of 3-D geometry using 2-D profiles and a geometric prior model

A method has been developed to reconstruct three-dimensional (3-D) surfaces from two-dimensional (2-D) projection data. It is used to produce individualized boundary element models, consisting of thorax and lung surfaces, for electro- and magnetocardiographic inverse problems. Two orthogonal projections are utilized. A geometrical prior model, built using segmented magnetic resonance images, is deformed according to profiles segmented from projection images. In the authors' method, virtual X-ray images of the prior model are first constructed by simulating real X-ray imaging. The 2-D profiles of the model are segmented from the projections and elastically matched with the profiles segmented from patient data. The displacement vectors produced by the elastic 2-D matching are back projected onto the 3-D surface of the prior model. Finally, the model is deformed, using the back-projected vectors. Two different deformation methods are proposed. The accuracy of the method is validated by a simulation. The average reconstruction error of a thorax and lungs was 1.22 voxels, corresponding to about 5 mm     

Estimation of the heart respiratory motion with applications for cone beam computed tomography imaging: a simulation study

Computed tomography (CT) reconstruction methods assume imaging of static objects; object movement during projection data acquisition causes tomogram artifacts. The continuously moving heart, therefore, represents a complicated imaging case. The associated problems due to the heart beating can be overcome either by using very short projection acquisition times, during which the heart may be considered static, or by ECG-gated acquisition. In the latter case, however, the acquisition of a large number of projections may not be completed in a single breath hold, thus heart displacement occurs as an additional problem. This problem has been addressed by applying heart motion models in various respiratory motion compensation algorithms. Our paper focuses on cone beam computed tomography (CBCT), performed in conjunction with isocentric, fluoroscopic equipment, and continuous ECG and respiratory monitoring. Such equipment is used primarily for in-theater three-dimensional (3-D) imaging and benefits particularly from the recent developments in flat panel detector technologies. The objectives of this paper are: (i) to develop a model for the motion of the heart due to respiration during the respiratory cycle; (ii) to apply this model to the tomographic reconstruction algorithm, in order to account for heart movement due to respiration in the reconstruction; and (iii) to initially evaluate this method by means of simulation studies. Based on simulation studies, we were able to demonstrate that heart displacement due to respiration can be estimated from the same projection data, required for a CBCT reconstruction. Our paper includes semiautomatic segmentation of the heart on the X-ray projections and reconstruction of a convex 3-D-heart object that performs the same motion as the heart during respiration, and use of this information into the CBCT reconstruction algorithm. The results reveal significant image quality improvements in cardiac image reconstruction.     

Evaluation of iterative sparse object reconstruction from few projections for 3-D rotational coronary angiography

A 3-D reconstruction of the coronary arteries offers great advantages in the diagnosis and treatment of cardiovascular disease, compared to 2-D X-ray angiograms. Besides improved roadmapping, quantitative vessel analysis is possible. Due to the heart's motion, rotational coronary angiography typically provides only 5–10 projections for the reconstruction of each cardiac phase, which leads to a strongly undersampled reconstruction problem. Such an ill-posed problem can be approached with regularized iterative methods. The coronary arteries cover only a small fraction of the reconstruction volume. Therefore, the minimization of the ${mbi L}_1$ norm of the reconstructed image, favoring spatially sparse images, is a suitable regularization. Additional problems are overlaid background structures and projection truncation, which can be alleviated by background reduction using a morphological top-hat filter. This paper quantitatively evaluates image reconstruction based on these ideas on software phantom data, in terms of reconstructed absorption coefficients and vessel radii. Results for different algorithms and different input data sets are compared. First results for electrocardiogram-gated reconstruction from clinical catheter-based rotational X-ray coronary angiography are presented. Excellent 3-D image quality can be achieved.      

A cone-beam reconstruction algorithm for circle-plus-arc data-acquisition geometry

In cone-beam computerized tomography (CT), projections acquired with the focal spot constrained on a planar orbit cannot provide a complete set of data to reconstruct the object function exactly. There are severe distortions in the reconstructed noncentral transverse planes when the cone angle is large. In this work, a new method is proposed which can obtain a complete set of data by acquiring cone-beam projections along a circle-plus-arc orbit. A reconstruction algorithm using this circle-plus-arc orbit is developed, based on the Radon transform and Grangeat's formula. This algorithm first transforms the cone-beam projection data of an object to the first derivative of the three-dimensional (3-D) Radon transform, using Grangeat's formula, and then reconstructs the object using the inverse Radon transform. In order to reduce interpolation errors, new rebinning equations have been derived accurately, which allows one-dimensional (1-D) interpolation to be used in the rebinning process instead of 3-D interpolation. A noise-free Defrise phantom and a Poisson noise-added Shepp-Logan phantom were simulated and reconstructed for algorithm validation. The results from the computer simulation indicate that the new cone-beam data-acquisition scheme can provide a complete set of projection data and the image reconstruction algorithm can achieve exact reconstruction. Potentially, the algorithm can be applied in practice for both a standard CT gantry-based volume tomographic imaging system and a C-arm-based cone-beam tomographic imaging system, with little mechanical modification required.     

Bayesian approach with the maximum entropy principle in image reconstruction from microwave scattered field data

Microwave imaging is of great interest in medical applications owing to its high sensitivity with respect to dielectric properties. It allows detection of very small inhomogeneities. The image reconstruction employing the microwave inverse scattering consists of reconstructing the image of an object from the scattered field measured behind the object. This reconstruction runs up against the nonuniqueness of the solution of the inverse scattering problem. The authors propose to solve the ill-posed inverse problem by a statistical regularization method based on the Bayesian maximum a posteriori estimation where the principle of maximum entropy is used for assigning the a priori laws. The results obtained demonstrate the power and potential of this method in image reconstruction.     

Optimal CT scanning plan for long-bone 3-D reconstruction

Digital computed tomographic (CT) data are widely used in three-dimensional (3-D) construction of bone geometry and density features for 3-D modelling purposes. During in vivo CT data acquisition the number of scans must be limited in order to protect patients from the risks related to X-ray absorption. The aim of this work is to automatically define, given a finite number of CT slices, the scanning plan which returns the optimal 3-D reconstruction of a bone segment from in vivo acquired CT images. An optimization algorithm based on a Discard-Insert-Exchange technique has been developed. In the proposed method the optimal scanning sequence is searched by minimizing the overall reconstruction error of a two-dimensional (2-D) prescanning image: an anterior-posterior (AP) X-ray projection of the bone segment. This approach has been validated in vitro on 3 different femurs. The 3-D reconstruction errors obtained through the optimization of the scanning plan on the 3-D prescanning images and on the corresponding 3-D data sets have been compared. 2-D and 3-D data sets have been reconstructed by linear interpolation along the longitudinal axis. Results show that direct 3-D optimization yields root mean square reconstruction errors which are only 4%-7% lower than the 2-D-optimized plan, thus proving that 2-D-optimization provides a good suboptimal scanning plan for 3-D reconstruction. Further on, 3-D reconstruction errors given by the optimized scanning plan and a standard radiological protocol for long bones have been compared. Results show that the optimized plan yields 20%-50% lower 3-D reconstruction errors     

多数据源电离层层析成像方法

卫星信标电离层层析成像(CIT)是一种高度不适定性问题,有限视角、稀疏布站等原因造成电离层CIT数据采集严重不完整,CIT结果精度不高,尤其是垂直误差可达50km以上。为了解决CIT结果精度不高的问题,提出一种融合地面垂测、斜测数据、顶部探测仪数据和三频卫星信标数据的多数据源电离层联合CIT方法。该方法采用建立在实测数据基础上的迭代初值,与实际偏离较小,能够提高CIT反演结果的精度;同时地基设备和顶部探测设备有较高的垂直分辨率,与卫星信标结合能够有效解决有限视角问题,提高CIT结果的垂直分辨率。     

Iterative Reconstruction?Reprojection: An Algorithm for Limited Data Cardiac-Computed Tomography

Cardiac X-ray computed tomography (CT) has been limited due to scanning times which are considerably longer (1 s) than required to resolve the beating heart (0.1 s). The otherwise attractive convolution-backprojection algorithm is not suited for CT image reconstruction from measurements comprising an incomplete set of projection data. In this paper, an iterative reconstruction-reprojection (IRR) algorithm is proposed for limited projection data CT image reconstruction. At each iteration, the missing views are estimated based on reprojection, which is a software substitute for the scanning process. The standard fan-beam convolution-backprojection algorithm is then used for image reconstruction. The proposed IRR algorithm enables the use of convolution-backprojection in limited angle of view and in limited field of view CT cases. The potential of this method for cardiac CT reconstruction is demonstrated using computer simulated data.     

Polygonal and polyhedral contour reconstruction in computed tomography

This paper is about three-dimensional (3-D) reconstruction of a binary image from its X-ray tomographic data. We study the special case of a compact uniform polyhedron totally included in a uniform background and directly perform the polyhedral surface estimation. We formulate this problem as a nonlinear inverse problem using the Bayesian framework. Vertice estimation is done without using a voxel approximation of the 3-D image. It is based on the construction and optimization of a regularized criterion that accounts for surface smoothness. We investigate original deterministic local algorithms, based on the exact computation of the line projections, their update, and their derivatives with respect to the vertice coordinates. Results are first derived in the two-dimensional (2-D) case, which consists of reconstructing a 2-D object of deformable polygonal contour from its tomographic data. Then, we investigate the 3-D extension that requires technical adaptations. Simulation results illustrate the performance of polygonal and polyhedral reconstruction algorithms in terms of quality and computation time.     

Data Specific Spatially Varying Regularization for Multimodal Fluorescence Molecular Tomography

Fluorescence molecular tomography (FMT) allows in vivo localization and quantification of fluorescence biodistributions in whole animals. The ill-posed nature of the tomographic reconstruction problem, however, limits the attainable resolution. Improvements in resolution and overall imaging performance can be achieved by forming image priors from geometric information obtained by a secondary anatomical or functional high-resolution imaging modality such as X-ray computed tomography or magnetic resonance imaging. A particular challenge in using image priors is to avoid the use of assumptions that may bias the solution and reduced the accuracy of the inverse problem. This is particularly relevant in FMT inversions where there is not an evident link between secondary geometric information and the underlying fluorescence biodistribution. We present here a new, two step approach to incorporating structural priors into the FMT inverse problem. By using the anatomic information to define a low dimensional inverse problem, we obtain a solution which we then use to determine the parameters defining a spatially varying regularization matrix for the full resolution problem. The regularization term is thus customized for each data set and is guided by the data rather than depending only on user defined a priori assumptions. Results are presented for both simulated and experimental data sets, and show significant improvements in image quality as compared to traditional regularization techniques.      

Implementation of uniform and simultaneous ART for 3-D reconstruction in an X-ray imaging system

The authors propose a 3-D volume reconstruction method using X-ray images with a calibration method to implement it in an X-ray imaging system. Previously the authors have proposed an advanced 3-D reconstruction algorithm based on an algebraic reconstruction technique (ART), called a uniform and simultaneous ART (USART). In practice, however, there are two main issues in implementing it in a realised X-ray imaging system. The first one is the huge computation time and memory required in achieving 3-D volume, which is a common limitation in most ART methods. The second issue is the system calibration for determining the geometry of the X-ray imaging conditions needed for the ART method. These two critical problems are addressed. A fast computing model of USART is proposed, where spherical voxel elements are employed in computation to reduce the computation time and memory. Then, a calibration method is proposed to identify the X-ray imaging geometry based on a cone beam projection model. For this purpose, a set of X-ray images of a reference grid pattern is used and the X-ray source positions are determined from the analysis of the image features, the centres of the grid points in the X-ray images. The validity of the proposed 3-D reconstruction method is investigated using a series of experiments.     

基于联合代数重建技术的介观荧光分子层析成像

针对介观荧光分子层析成像重建问题,提出了一种基于联合代数重建技术的介观荧光分子层析成像重建方法.首先应用主成分分析算法对敏感矩阵进行双重降维操作,以消除敏感矩阵中的冗余信息;其次为保持重建结果与目标数据的一致性,对降维后的矩阵进行零填充;最后应用荧光光强测量数据和填充后的敏感矩阵,经带有总变差正则化项的联合代数重建技术...     

基于一种新的非局部二次MRF先验模型的Bayesian图像重建

Bayesian方法或最大后验(maximum a posteriori,MAP)法被认为是解决图像重建中的病态问题的有效方法.根据Bayesian理论,目标图像的先验信息被加诸于图像重建中来抑制噪声.然而,大部分的先验模型提供的先验信息来自于一个较小的局部邻域内灰度值的简单加权差,只能对Bayesian重建提供有限的先验信息.本文提出一个新的非局部的且具有二次的先验能量方程的马尔可夫随机场(Markov Random Fields,MRS)先验,该先验通过选择较大邻域和新的加权方式来充分利用图像的全局信息.文章同时还提出使用该非局部先验的Bayesian重建的迭代算法.最后给出了该先验在PET(正电子发射成像)重建中的应用.实验结果以及同其他先验的比较证明该先验在降低噪声效果和保持边缘方面具有很好的表现.     

Bayesian reconstruction and use of anatomical a priori information for emission tomography

A Bayesian method is presented for simultaneously segmenting and reconstructing emission computed tomography (ECT) images and for incorporating high-resolution, anatomical information into those reconstructions. The anatomical information is often available from other imaging modalities such as computed tomography (CT) or magnetic resonance imaging (MRI). The Bayesian procedure models the ECT radiopharmaceutical distribution as consisting of regions, such that radiopharmaceutical activity is similar throughout each region. It estimates the number of regions, the mean activity of each region, and the region classification and mean activity of each voxel. Anatomical information is incorporated by assigning higher prior probabilities to ECT segmentations in which each ECT region stays within a single anatomical region. This approach is effective because anatomical tissue type often strongly influences radiopharmaceutical uptake. The Bayesian procedure is evaluated using physically acquired single-photon emission computed tomography (SPECT) projection data and MRI for the three-dimensional (3-D) Hoffman brain phantom. A clinically realistic count level is used. A cold lesion within the brain phantom is created during the SPECT scan but not during the MRI to demonstrate that the estimation procedure can detect ECT structure that is not present anatomically.     

基于轮廓先验约束的复杂异形工件CT成像方法研究

在X射线CT成像检测系统中,复杂异形工件由于结构复杂,散射、硬化等现象较严重,而且在部分投影角度上,由于在射线透照方向上的有效厚度差异大,固定能量的射线剂量与厚度不匹配,投影数据质量较差,传统的CT重建算法得到的重建图像质量低,边缘模糊,无法获取工件的完整轮廓信息.为此文章研究了一种基于轮廓先验约束的复杂异形工件CT成像方法.首先利用双目立体视觉技术获取工件的轮廓信息,根据双目坐标系与CT坐标系间的空间位置关系完成先验图像的配准;然后将轮廓先验纳入到CT重建过程中并结合TV正则化进行轮廓约束重建.实验结果表明,该方法能够有效地抑制伪影和噪声,保留重建图像的边缘,改善重建图像质量,有助于提高复杂异形工件缺陷检测的可靠性.     

A Hybrid Reconstruction Algorithm for 3-D Ionospheric Tomography

In this paper, a hybrid reconstruction algorithm (HRA) is presented to solve the ill-posed inverse problem associated with 3-D ionospheric stochastic tomography. In this new method, the ionospheric electron density (IED) can be inverted by using two steps. First, a truncated singular value decomposition (TSVD) method, whose value is independent on any initial estimation, is used to resolve the ill-posed problem of the tomography system. Second, taking into account the "approximation" of its solution, an iterative improvement process of the solution is then implemented by utilizing the conventional algebraic reconstruction algorithm (ART). The HRA, therefore, offers a more reasonable approach to choose an initial approximate for the ART and to improve the quality of the final reconstructed image. A simulated experiment demonstrates that the HRA method is superior to the TSVD or the ART alone for the tomographic inversion of IED. Finally, the HRA is used to perform GPS-based tomographic reconstruction of the IED at mid- and low-latitude regions.     

Improved performance of bayesian solutions for inverse electrocardiography using multiple information sources

The usual goal in inverse electrocardiography (ECG) is to reconstruct cardiac electrical sources from body surface potentials and a mathematical model that relates the sources to the measurements. Due to attenuation and smoothing that occurs in the thorax, the inverse ECG problem is ill-posed and imposition of a priori constraints is needed to combat this ill-posedness. When the problem is posed in terms of reconstructing heart surface potentials, solutions have not yet achieved clinical utility; limitations include the limited availability of good a priori information about the solution and the lack of a "good" error metric. We describe an approach that combines body surface measurements and standard forward models with two additional information sources: statistical prior information about epicardial potential distributions and sparse simultaneous measurements of epicardial potentials made with multielectrode coronary venous catheters. We employ a Bayesian methodology which offers a general way to incorporate these information sources and additionally provides statistical performance analysis tools. In a simulation study, we first compare solutions using one or more of these information sources. Then, we study the effects of varying the number of sparse epicardial potential measurements on reconstruction accuracy. To evaluate accuracy, we used the Bayesian error covariance as well as traditional error metrics such as relative error. Our results show that including even sparsely sampled information from coronary venous catheters can substantially improve the reconstruction of epicardial potential distributions and that a Bayesian framework provides a feasible approach to using this information. Moreover, computing the Bayesian error standard deviations offers a means to indicate confidence in the results even in the absence of validation data.     

Cone-beam reprojection using projection-matrices

This paper addresses reprojection of three-dimensional (3-D) reconstructions obtained from cone-beam scans using a C-arm imaging equipment assisted by a pose-determining system. The emphasis is on reprojecting without decomposing the estimated projection matrix (P-matrix) associated with a pose. Both voxel- and ray-driven methods are considered. The voxel-driven reprojector follows the algorithm for backprojection using a P-matrix. The ray-driven reprojector is derived by extracting from the P-matrix the equation of the line joining a detector-pixel and the X-ray source position. This reprojector can be modified to a ray-driven backprojector. When the geometry is specified explicitly in terms of the physical parameters of the imaging system, the projection matrices can be constructed. The resulting "projection-matrix method" is advantageous, especially when the scanning trajectory is irregular. The algorithms presented are useful in iterative methods of image reconstruction and enhancement procedures, apart from their well-known role in visualization and volume rendering. Reprojections of 3-D patient data compare favorably with the original X-ray projections obtained from a prototype C-arm system. The algorithms for reprojection can be modified to compute perspective maximum intensity projection.