No Arabic abstract
The purpose of the present work is the study of reconstruction properties of a new Molecular Breast Imaging (MBI) device for the early diagnosis of breast cancer, in Limited Angle Tomography (LAT), by using two asymmetric detector heads with different collimators. The detectors face each other in anti-parallel viewing direction and, mild-compressing the breast phantom, they are able to reconstruct the inner tumour of the phantoms with only a limited number of projections using a dedicated maximum-likelihood expectation maximization (ML-EM) algorithm. Phantoms, MBI system, as well as Monte Carlo simulator using Geant 4 Application for Tomographic Emission (GATE) software, are briefly described. MBI systems model has been implemented in IDL (Interactive Data Visualization), in order to evaluate the best LAT configuration of the system and its reconstruction ability by varying tumours size, depth and uptake. LAT setup in real and simulated configurations, as well as the ML-EM method and the preliminary reconstruction results, are discussed.
A new set of signals for studying detectability of an x-ray imaging system is presented. The results obtained with these signals are intended to complement the NEQ results. The signals are generated from line spread profiles by progressively removing their lower frequency components and the resulting high frequency residues (HFRs) form the set of signals to be used in detectability studies. Detectability indexes for these HFRs are obtained using a non-prewhitening (NPW) observer and a series of edge images are used to obtain the HFRs, the covariance matrices required by the NPW model and the MTF and NPS used in NEQ calculations. The template used in the model is obtained by simulating the processes of blurring and sampling of the edge images. Comparison between detectability indexes for the HFRs and NEQ are carried out for different acquisition techniques using different beam qualities and doses. The relative sensitivity shown by detectability indexes using HFRs is higher than that of NEQ, especially at lower doses. Also, the different observers produce different results at high doses: while the ideal Bayesian observer used by NEQ distinguishes between beam qualities, the NPW used with the HFRs produces no differences between them. Delta functions used in HFR are the opposite of complex exponential functions in terms of their support in the spatial and frequency domains. Since NEQ can be interpreted as detectability of these complex exponential functions, detectability of HFRs is presented as a natural complement to NEQ in the performance assessment of an imaging system.
This paper considers the reconstruction problem in Acousto-Electrical Tomography, i.e., the problem of estimating a spatially varying conductivity in a bounded domain from measurements of the internal power densities resulting from different prescribed boundary conditions. Particular emphasis is placed on the limited angle scenario, in which the boundary conditions are supported only on a part of the boundary. The reconstruction problem is formulated as an optimization problem in a Hilbert space setting and solved using Landweber iteration. The resulting algorithm is implemented numerically in two spatial dimensions and tested on simulated data. The results quantify the intuition that features close to the measurement boundary are stably reconstructed and features further away are less well reconstructed. Finally, the ill-posedness of the limited angle problem is quantified numerically using the singular value decomposition of the corresponding linearized problem.
Parallel magnetic resonance imaging (MRI) is a technique of image acceleration which takes advantage of the localization of the field of view (FOV) of coils in an array. In this letter we show that metamaterial lenses based on capacitively-loaded rings can provide higher localization of the FOV. Several lens designs are systematically analyzed in order to find the structure providing higher signal-to-noise-ratio. The magnetoinductive (MI) lens is find to be the optimum structure and an experiment is developed to show it. The ability of the fabricated MI lenses to accelerate the image is quantified by means of the parameter known in the MRI community as g-factor.
We investigate the reconstruction problem of limited angle tomography. Such problems arise naturally in applications like digital breast tomosynthesis, dental tomography, electron microscopy etc. Since the acquired tomographic data is highly incomplete, the reconstruction problem is severely ill-posed and the traditional reconstruction methods, such as filtered backprojection (FBP), do not perform well in such situations. To stabilize the reconstruction procedure additional prior knowledge about the unknown object has to be integrated into the reconstruction process. In this work, we propose the use of the sparse regularization technique in combination with curvelets. We argue that this technique gives rise to an edge-preserving reconstruction. Moreover, we show that the dimension of the problem can be significantly reduced in the curvelet domain. To this end, we give a characterization of the kernel of limited angle Radon transform in terms of curvelets and derive a characterization of solutions obtained through curvelet sparse regularization. In numerical experiments, we will present the practical relevance of these results.
In this work we introduce a new method that combines Parallel MRI and Compressed Sensing (CS) for accelerated image reconstruction from subsampled k-space data. The method first computes a convolved image, which gives the convolution between a user-defined kernel and the unknown MR image, and then reconstructs the image by CS-based image deblurring, in which CS is applied for removing the inherent blur stemming from the convolution process. This method is hence termed CORE-Deblur. Retrospective subsampling experiments with data from a numerical brain phantom and in-vivo 7T brain scans showed that CORE-Deblur produced high-quality reconstructions, comparable to those of a conventional CS method, while reducing the number of iterations by a factor of 10 or more. The average Normalized Root Mean Square Error (NRMSE) obtained by CORE-Deblur for the in-vivo datasets was 0.016. CORE-Deblur also exhibited robustness regarding the chosen kernel and compatibility with various k-space subsampling schemes, ranging from regular to random. In summary, CORE-Deblur enables high quality reconstructions and reduction of the CS iterations number by 10-fold.