skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Numerical Computation of X-ray Computerized Tomography from partial measurements
The present paper proposes a novel numerical scheme to X-ray Computerized Tomography (CT) from partial measurement data. In order to reduce radiation exposure, it is desirable to irradiate X-ray only around region of interest (ROI), while the conventional reconstruction methods such as filtered back projection (FBP) could not work due to its intrinsic limitation of dependency on whole measurement data. The proposed method gives a direct numerical reconstruction employing a Cauchy type boundary integration in $$A$$-analytic theory and a singular integral equation which maps boundary measurement to interior data. Numerical examples using experimental data are also exhibited to show validity of the proposed numerical procedure.  more » « less
Award ID(s):
1907097
PAR ID:
10332388
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Keisan suri kogaku ronbunshu
Volume:
12
Issue:
2021
ISSN:
1348-5245
Page Range / eLocation ID:
37-43
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, Keisan suuri kougaku rombunshuu)
    This paper presents a new algorithm for X-ray Computerized Tomography (CT) based on Bukhgeim’s theory of analytic maps. The reconstruction relies on a Cauchy-type integral formula, where the integration over the boundary replaces the integration in the back- projection operator used in existing algorithms. From the numerical computation stand point, the proposed method recovers the attenuation coefficient at arbitrarily points by utilizing the boundary integration without internal global meshes. This means that it achieves high-parallel efficiency, and it reduces computational resources. Some numerical examples are presented to show feasibility of the proposed algorithm. 
    more » « less
  2. ; (, Advances in mathematical sciences and applications)
    This work concerns the numerical realization of a Cauchy-type integral formula for sequence valued analytic functions in the sense of Bukhgeim, and its applications to the source reconstruction problem in inverse radiative transport through a non-absorbing and non-scattering medium. The inverse source problem is mathematically equivalent to the classical X-ray Computed Tomography (CT), where a function is to be determined from its line integrals. The proposed algorithms have the added advantage to extend to the source determination problems in media with absorbing and scattering properties. Such extensions cannot be achieved in the existing X-ray CT algorithms. The numerical experiments demonstrate the feasibility of our new tomographic algorithms. 
    more » « less
  3. ; ; ; ; (, Optics Express)
    Compressive spectral X-ray imaging (CSXI) introduces a pixelated spectral modulator called K-edge coded aperture (KCA) in front of the X-ray source, which enables both, lower dosage to the subject, as well as the capability of spectral tomography while using low-cost integrating X-ray detectors. CSXI systems generally use hundreds of different spectral modulators, each with a distinct pattern to uniquely modulate the illumination at every view angle. In contrast, this paper introduces the use of a single and static coded aperture placed in a tomosynthesis gantry. The compressive system thus interrogates the subject with a fixed coded illumination pattern on all view angles. The advantages of the system are many including reduced cost and the feasibility of implementation. Given the reduced set of coded measurement and the limited spectral separation ability in the resulting architecture, the nonlinear inverse reconstruction problem results in a highly ill-posed problem. An efficient alternating minimization method with three-dimensional total variation regularization is developed for image reconstruction. Furthermore, rather than simply using a random pattern, the coded aperture is optimized under a uniform sensing criterion that shapes the spatial and spectral pattern of the coded aperture so as to minimize the overall radiation exposure placed on any volumetric area of the patient. This is of particular importance in medical imaging where patients at risk are recommended to have periodical X-ray tomosynthesis screenings. The coded aperture optimization is then posed as a binary programming problem solved by a gradient-based algorithm with equilibrium constraints. Numerical experiments show that spatial and spectral coding used in the proposed system to interrogate the subject not only reduces the radiation dose but it also improves the quality of image reconstruction. Gains close to 5dB in peak signal to noise ratio are observed in simulations. Furthermore, it is shown that the optimization of the KCA can effectively improve the uniformity of X-ray radiation compared to random KCA modulation, thus reducing the radiation dose throughout all volumetric sub-areas of the subject — an objective that is not possible with the use of random KCAs. 
    more » « less
  4. ; (, Medical Physics)
    Abstract PurposeForward and backprojections are the basis of all model‐based iterative reconstruction (MBIR) methods. However, computing these accurately is time‐consuming. In this paper, we present a method for MBIR in parallel X‐ray beam geometry that utilizes a Gram filter to efficiently implement forward and backprojection. MethodsWe propose using voxel‐basis and modeling its footprint in a box spline framework to calculate the Gram filter exactly and improve the performance of backprojection. In the special case of parallel X‐ray beam geometry, the forward and backprojection can be implemented by an estimated Gram filter efficiently if the sinogram signal is bandlimited. In this paper, a specialized sinogram interpolation method is proposed to eliminate the bandlimited prerequisite and thus improve the reconstruction accuracy. We build on this idea by utilizing the continuity of the voxel‐basis' footprint, which provides a more accurate sinogram interpolation and further improves the efficiency and quality of backprojection. In addition, the detector blur effect can be efficiently accounted for in our method to better handle realistic scenarios. ResultsThe proposed method is tested on both phantom and real computed tomography (CT) images under different resolutions, sinogram sampling steps, and noise levels. The proposed method consistently outperforms other state‐of‐the‐art projection models in terms of speed and accuracy for both backprojection and reconstruction. ConclusionsWe proposed a iterative reconstruction methodology for 3D parallel‐beam X‐ray CT reconstruction. Our experimental results demonstrate that the proposed methodology is accurate, fast, and reproducible, and outperforms alternative state‐of‐the‐art projection models on both backprojection and reconstruction results significantly. 
    more » « less
  5. ; (, Optics express)
    Abstract: Coded aperture X-ray computed tomography (CT) has the potential to revolutionize X-ray tomography systems in medical imaging and air and rail transit security - both areas of global importance. It allows either a reduced set of measurements in X-ray CT without degrada- tion in image reconstruction, or measure multiplexed X-rays to simplify the sensing geometry. Measurement reduction is of particular interest in medical imaging to reduce radiation, and airport security often imposes practical constraints leading to limited angle geometries. Coded aperture compressive X-ray CT places a coded aperture pattern in front of the X-ray source in order to obtain patterned projections onto a detector. Compressive sensing (CS) reconstruction algorithms are then used to recover the image. To date, the coded illumination patterns used in conventional CT systems have been random. This paper addresses the code optimization prob- lem for general tomography imaging based on the point spread function (PSF) of the system, which is used as a measure of the sensing matrix quality which connects to the restricted isom- etry property (RIP) and coherence of the sensing matrix. The methods presented are general, simple to use, and can be easily extended to other imaging systems. Simulations are presented where the peak signal to noise ratios (PSNR) of the reconstructed images using optimized coded apertures exhibit significant gain over those attained by random coded apertures. Additionally, results using real X-ray tomography projections are presented. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email