A subspace-based resolution-enhancing image reconstruction method for few-view differential phase-contrast tomography

It is well known that properly designed image reconstruction methods can facilitate reductions in imaging doses and data-acquisition times in tomographic imaging. The ability to do so is particularly important for emerging modalities, such as differential x-ray phase-contrast tomography (D-XPCT), which are currently limited by these factors. An important application of D-XPCT is high-resolution imaging of biomedical samples. However, reconstructing high-resolution images from few-view tomographic measurements remains a challenging task due to the high-frequency information loss caused by data incompleteness. In this work, a subspace-based reconstruction strategy is proposed and investigated for use in few-view D-XPCT image reconstruction. By adopting a two-step approach, the proposed method can simultaneously recover high-frequency details within a certain region of interest while suppressing noise and/or artifacts globally. The proposed method is investigated by the use of few-view experimental data acquired by an edge-illumination D-XPCT scanner

Error analysis for filtered back projection reconstructions in Besov spaces

Filtered back projection (FBP) methods are the most widely used reconstruction algorithms in computerized tomography (CT). The ill-posedness of this inverse problem allows only an approximate reconstruction for given noisy data. Studying the resulting reconstruction error has been a most active field of research in the 1990s and has recently been revived in terms of optimal filter design and estimating the FBP approximation errors in general Sobolev spaces. However, the choice of Sobolev spaces is suboptimal for characterizing typical CT reconstructions. A widely used model are sums of characteristic functions, which are better modelled in terms of Besov spaces $\mathrm{B}^{\alpha,p}_q(\mathbb{R}^2)$. In particular $\mathrm{B}^{\alpha,1}_1(\mathbb{R}^2)$ with $\alpha \approx 1$ is a preferred model in image analysis for describing natural images. In case of noisy Radon data the total FBP reconstruction error $$\|f-f_L^\delta\| \le \|f-f_L\|+ \|f_L - f_L^\delta\|$$ splits into an approximation error and a data error, where $L$ serves as regularization parameter. In this paper, we study the approximation error of FBP reconstructions for target functions $f \in \mathrm{L}^1(\mathbb{R}^2) \cap \mathrm{B}^{\alpha,p}_q(\mathbb{R}^2)$ with positive $\alpha \not\in \mathbb{N}$ and $1 \leq p,q \leq \infty$. We prove that the $\mathrm{L}^p$-norm of the inherent FBP approximation error $f-f_L$ can be bounded above by \begin{equation*} \|f - f_L\|_{\mathrm{L}^p(\mathbb{R}^2)} \leq c_{\alpha,q,W} \, L^{-\alpha} \, |f|_{\mathrm{B}^{\alpha,p}_q(\mathbb{R}^2)} \end{equation*} under suitable assumptions on the utilized low-pass filter's window function $W$. This then extends by classical methods to estimates for the total reconstruction error.Comment: 32 pages, 8 figure

System Characterizations and Optimized Reconstruction Methods for Novel X-ray Imaging

In the past decade there have been many new emerging X-ray based imaging technologies developed for different diagnostic purposes or imaging tasks. However, there exist one or more specific problems that prevent them from being effectively or efficiently employed. In this dissertation, four different novel X-ray based imaging technologies are discussed, including propagation-based phase-contrast (PB-XPC) tomosynthesis, differential X-ray phase-contrast tomography (D-XPCT), projection-based dual-energy computed radiography (DECR), and tetrahedron beam computed tomography (TBCT). System characteristics are analyzed or optimized reconstruction methods are proposed for these imaging modalities. In the first part, we investigated the unique properties of propagation-based phase-contrast imaging technique when combined with the X-ray tomosynthesis. Fourier slice theorem implies that the high frequency components collected in the tomosynthesis data can be more reliably reconstructed. It is observed that the fringes or boundary enhancement introduced by the phase-contrast effects can serve as an accurate indicator of the true depth position in the tomosynthesis in-plane image. In the second part, we derived a sub-space framework to reconstruct images from few-view D-XPCT data set. By introducing a proper mask, the high frequency contents of the image can be theoretically preserved in a certain region of interest. A two-step reconstruction strategy is developed to mitigate the risk of subtle structures being oversmoothed when the commonly used total-variation regularization is employed in the conventional iterative framework. In the thirt part, we proposed a practical method to improve the quantitative accuracy of the projection-based dual-energy material decomposition. It is demonstrated that applying a total-projection-length constraint along with the dual-energy measurements can achieve a stabilized numerical solution of the decomposition problem, thus overcoming the disadvantages of the conventional approach that was extremely sensitive to noise corruption. In the final part, we described the modified filtered backprojection and iterative image reconstruction algorithms specifically developed for TBCT. Special parallelization strategies are designed to facilitate the use of GPU computing, showing demonstrated capability of producing high quality reconstructed volumetric images with a super fast computational speed. For all the investigations mentioned above, both simulation and experimental studies have been conducted to demonstrate the feasibility and effectiveness of the proposed methodologies

Third International Conference on Inverse Design Concepts and Optimization in Engineering Sciences (ICIDES-3)

Papers from the Third International Conference on Inverse Design Concepts and Optimization in Engineering Sciences (ICIDES) are presented. The papers discuss current research in the general field of inverse, semi-inverse, and direct design and optimization in engineering sciences. The rapid growth of this relatively new field is due to the availability of faster and larger computing machines