In this paper, we propose a new wavelet-based reconstruction method suited to three-dimensional (3-D) cone-beam (CB) tomography. It is derived from the Feldkamp algorithm and is valid for the same geometrical conditions. The demonstration is done in the framework of nonseparable wavelets and requires ideally radial wavelets. The proposed inversion formula yields to a filtered backprojection algorithm but the filtering step is implemented using quincunx wavelet filters. The proposed algorithm reconstructs slice by slice both the wavelet and approximation coefficients of the 3-D image directly from the CB projection data. The validity of this multiresolution approach is demonstrated on simulations from both mathematical phantoms and 3-D rotational angiography clinical data. The same quality is achieved compared with the standard Feldkamp algorithm, but in addition, the multiresolution decomposition allows to apply directly image processing techniques in the wavelet domain during the inversion process. As an example, a fast low-resolution reconstruction of the 3-D arterial vessels with the progressive addition of details in a region of interest is demonstrated. Other promising applications are the improvement of image quality by denoising techniques and also the reduction of computing time using the space localization of wavelets.

In this paper, a new multiresolution reconstruction approach for fan-beam tomography is established. The theoretical development assumes radial wavelets. An approximate reconstruction formula based on a near-radial quincunx multiresolution scheme is proposed. This multiresolution algorithm allows to compute both the quincunx approximation and detail coefficients of an image from its fan-beam projections. Simulations on mathematical phantoms show that wavelet decomposition is acceptable for small beam angles but deteriorates at high angles. The main applications of the method are denoising and wavelet-based image analysis.

In this paper, the use of nonseparable wavelets for tomographic reconstruction is investigated. Local tomography is also presented. The algorithm computes both the quincunx approximation and detail coefficients of a function from its projections. Simulation results showed that nonseparable wavelets provide a reconstruction improvement versus separable wavelets.