Computed tomography (CT) is a powerful tool for reconstruction and analysis of inner structure of objects applied in various fields. Although many classes of objects of interest may have highly absorbent inclusions, leading to a certain type of distortions on reconstructed volume images (metal-like artifacts). The correction of this type of artifacts can’t be considered a solved task, despite all the efforts in this direction. The development and research of methods for suppressing CT artifacts require high-quality synthetic data which allow for numerical assessment of the accuracy of the metal-like artifacts reduction methods and training of neural networks. Although simplified methods considering only beam hardening and Poisson photon distribution are commonly used to simulate the data with type of distortions. In present work we design experiments using the tomographic scanner of the Federal Research Center “Crystallography and Photonics” of the Russian Academy of Sciences to demonstrate that in some cases beam hardening may not be the dominant reason for the arising of metal-like artifacts. These experiments are closely analyzed and modeled within different approaches. The problems in both simplified and state of the art approaches are emphasized and discussed. The provided results show the importance of paying attention to the dark current modeling for synthesized data generation under the conditions of total photon absorption.
The rotation axis position is an important parameter of classical reconstruction algorithms in X-ray computed tomography (CT). The use of incorrect values of the axis position parameters during the reconstruction leads to the appearance of various artifacts distorting the reconstructed image. Therefore, to obtain a reconstruction of better quality, automatic rotation axis position determination and misalignment correction methods are of use. Most of the existing high-precision automatic rotation axis position determination methods are either fast, but suitable only within a parallel-beam geometric scheme, or indifferent to the geometric scheme, but computationally laborious. In this paper, we propose a method for auto-detection of two scalar parameters of rotation axis position — axis shift and tilt in the plane parallel to the detector window plane — using a pixel-wise arithmetically averaged projection image. The described method is highly accurate within both parallel-beam and cone-beam geometric schemes whereas it is characterized by robustness to noise in projection data. The method has performed an increase in reconstruction quality when compared with some well-known and still used in practice methods both on synthetic data and on real data obtained in real laboratory conditions.
In X-ray computed tomography (CT) the real rotation axis position often does not coincide with the assumed one: technical imperfections of the tomographic setup, the fast speed of movement of the gantry and goniometer cause rotation axis displacements and inclinations. At the same time the use of incorrect axis location parameters during reconstruction leads to the appearance of so-called tuning-fork artifacts in the form of stripes and blurs at the object boundary. The existing rotation axis alignment methods for cone-beam CT require a large amount of computing resources, are laborious in implementation, are not able to accurately determine several axis location parameters at once, or are based on the processing of additional equipped with reference markers object post-scans and shots that are not always available. Thus the rotation axis alignment methods development in the cone-beam CT still seems to be relevant. In this paper, the developed model for parameterizing the rotation axis position is described and justified. The novel several-stage automatic method for rotation axis parameters determination is described. The proposed method is based on usage of mean projection image and tested both on synthetic and real data in parallel-beam and cone-beam geometric schemes. The absolute error of that method on the simulated data is no more than 1 pixel and 1 degree, respectively for shift and slope.
In computed tomography (CT), the relative trajectories of a sample, a detector, and a signal source are traditionally considered to be known, since they are caused by the intentional preprogrammed movement of the instrument parts. However, due to the mechanical backlashes, rotation sensor measurement errors, thermal deformations real trajectory differs from desired ones. This negatively affects the resulting quality of tomographic reconstruction. Neither the calibration nor preliminary adjustments of the device completely eliminates the inaccuracy of the trajectory but significantly increase the cost of instrument maintenance. A number of approaches to this problem are based on an automatic refinement of the source and sensor position estimate relative to the sample for each projection (at each time step) during the reconstruction process. A similar problem of position refinement while observing different images of an object from different angles is well known in robotics (particularly, in mobile robots and selfdriving vehicles) and is called Simultaneous Localization And Mapping (SLAM). The scientific novelty of this work is to consider the problem of trajectory refinement in microtomography as a SLAM problem. This is achieved by extracting Speeded Up Robust Features (SURF) features from Xray projections, filtering matches with Random Sample Consensus (RANSAC), calculating angles between projections, and using them in factor graph in combination with stepper motor control signals in order to refine rotation angles.
In Computed tomography (CT) usage of common reconstruction algorithms to the projection data acquired with polychromatic probing radiation leads to the appearance of a cup-like distortion. CT image quality can be improved by adjusting the CT scanner or the reconstruction algorithm, but for this purpose assessment of cupping artifacts evaluation needs to be done. Existing assessment methods either rely on expert opinion or require an object binary mask, which can be unavailable. In this paper, we propose a method for blind assessment of cupping artifacts that do not require any prior information. The main idea of the proposed method is to evaluate the degree of change in intensity near automatically found edges of optically dense objects. We prove the applicability of the method on the collected dataset with cupping artifacts. The results show a monotonic dependency between the severity of cupping artifacts and the calculated with the proposed method value.
The method of Computed Tomography (CT) has progressed throughout the past decade with advances in CT apparatus and program parts that have resulted in an increasing number of CT applications. Today innovative CT Xray detectors have high spatial resolution till a tenth or hundredth of a micron. However, itsfield of view is significantly limited. The object being scanned with a high resolution does not always completely enter in (covered by) the field of view of the detector. The collected projections data may be incomplete. The use of incomplete data in classical reconstruction methods leads to image quality loss. This paper provides a new advanced reconstruction method that demonstrates image quality improvements compared with classical methods when incomplete data collected. The method uses the hypothesis about the consistency of object description in sinogram space and reconstruction space. Input data for the algorithm proposed are incomplete data, and the output data are the reconstructed image and the confidence values for all pixels of the image (reconstruction reliability). A detailed description of the algorithm is presented. Its quality characteristics are based on Shepp-Logan phantom studies.
Despite significant progress in computer vision, pattern recognition, and image analysis, artifacts in imaging still hampers the progress in many scientific fields relying on the results of image analysis. We here present an advanced image-based artifacts suppression algorithm for high-resolution tomography. The algorithm is based on guided filtering of a reconstructed image mapped from the Cartesian to the polar coordinates space. This postprocessing method efficiently reduces both ring- and radial streak artifacts in a reconstructed image. Radial streak artifacts can appear in tomography with an off-center rotation of a large object over 360 degrees used to increase the reconstruction field of view. We successfully applied the developed algorithm for improving x-ray phase-contrast images of human post-mortem pineal gland and olfactory bulbs.
The algorithm for 3D vector image reconstruction from a set of spectral tomographic projections collected with CT set-up completed with an optical element or elements inside the optical path behind the sample is proposed. The purpose of their placement into the optical path is to divide the integral polychromatic projection into a series of monochromatic projections, i.e., to get a multi-channel image. Understanding of the reconstruction results in the monochromatic case is beyond question, the relationship between the reconstructed spatial distribution of the linear attenuation coefficient and the discrete description of the elemental structure of the probed object is linear. In difference with monochromatic case the result of the reconstruction from polychromatic projections is a spatial distribution of the so-called effective or average attenuation coefficient, its connection to a discrete description of the elemental structure is nontrivial. However, if the distribution of the averaged coefficient is supplemented by distributions of linear coefficients for several energies, then it is possible to estimate of the local composition of the object. We present a model for the formation of spectral multi-channel projection based on crystal analyzer usage and describe the steps needed to solve the tomography inverse problem.
Computer vision for biomedical imaging applications is fast developing and at once demanding field of computer science. In particular, computer vision technique provides excellent results for detection and segmentation problems in tomographic imaging. X-ray phase contrast Tomography (XPCT) is a noninvasive 3D imaging technique with high sensitivity for soft tissues. Despite a considerable progress in XPCT data acquisition and data processing methods, the problem in degradation of image quality due to artifacts remains a widespread and often critical issue for computer vision applications. One of the main problems originates from a sample alteration during a long tomographic scan. We proposed and tested Simultaneous Iterative Reconstruction algorithm with Total Variation regularization to reduce the number of projections in high resolution XPCT scans of ex-vivo mouse spinal cord. We have shown that the proposed algorithm allows tenfold reducing the number of projections and, therefore, the exposure time, with conservation of the important morphological information in 3D image with quality acceptable for computer graphics and computer vision applications. Our research paves a way for more effective implementation of advanced computer technologies in phase contrast tomographic research.
Nowadays, microtomography experiments require a lot of time to collect data and process it. In order to observe realtime processes (e.g. fluid flow through porous media), measurements and calculations should be carried out fast enough, therefore optimization task should be solved. Two approaches were developed to solve it. The first one is associated with the search of optimal experimental parameters: number of projection and the quality of the detector. The second one is involved with representative elementary volume determination. Moreover, this determination technique is described in general terms and can be applied not only for porous media studies. Both algorithms are based on comparison methods of pore sizes distribution histograms. On this purposes, apart from common Earth Mover’s Distance (Wasserstein Distance) metric, a new Mean Vector Distance (MVD) metric was designed and described in this paper.
In this work, we propose a method for tomography reconstruction in case of a limited field of view, when the whole image of the investigated sample does not fit on the detector. Proposed technique based on iterative procedure with corrections on each step in sinogram space and reconstruction space. On synthetic and experimental data shown, that proposed technique allows to improve tomography reconstruction quality and extends the field of view.
Usage of common reconstruction algorithms like Filtered Back Projection and Algebraic Reconstruction Technique to the projection data acquired with poly-chromatic probing radiation leads to the appearance of a cup-like distortion of the value profile in reconstructed images. While many methods of the poly-chromatic probing artifacts suppression are suggested, the numerical estimation algorithm of the “Cupping effect” typically is not considered to be important. Described methods imply manual regions selection where the intensity will be compared, or just use experts’ opinion on the effect presence. In this paper, we suggest automatic estimation of the “Cupping effect” method based on utilizing the distance transform built using the objects mask. As a result, we obtain a numeric estimation of the intensity change from the border to the center of the object. As the final image index, a weighted sum of the ratings of all objects is used. While positive value shows the magnitude of the “Cupping effect”, a negative value, on the contrary, shows magnitude of the reverse “Cupping effect”. In the paper, we demonstrate the method used on simulated data and compare it with several different techniques for distortion evaluation due to poly-chromatic probing. Finally, we show method effectiveness on real data acquired with laboratory tomography.
Porous materials are widely used in different applications, in particular they are used to create various filters. Their quality depends on parameters that characterize the internal structure such as porosity, permeability and so on. Сomputed tomography (CT) allows one to see the internal structure of a porous object without destroying it. The result of tomography is a gray image. To evaluate the desired parameters, the image should be segmented. Traditional intensity threshold approaches did not reliably produce correct results due to limitations with CT images quality. Errors in the evaluation of characteristics of porous materials based on segmented images can lead to the incorrect estimation of their quality and consequently to the impossibility of exploitation, financial losses and even to accidents. It is difficult to perform correctly segmentation due to the strong difference in voxel intensities of the reconstructed object and the presence of noise. Image filtering as a preprocessing procedure is used to improve the quality of segmentation. Nevertheless, there is a problem of choosing an optimal filter. In this work, a method for selecting an optimal filter based on attributive indicator of porous objects (should be free from "levitating stones" inside of pores) is proposed. In this paper, we use real data where beam hardening artifacts are removed, which allows us to focus on the noise reduction process.
The present paper is devoted to the solution of a tomographic reconstruction problem of using a regularized algebraic approach for large scale data. The paper explores the issues related to the use of cone beam polychromatic computed tomography. An algorithm for regularized solution of the linear operator equation is described. The minimizing parametric composite function is given and step of the iterative procedure developed is written. The reconstructed volumetric image is about 60 billions voxels. It forces to divide the task of reconstruction of the full volume into subtasks for the efficient implementation of the reconstruction algorithm on the GPU. In each of the subtasks the current solution for the local volume of a given size is calculated. An approach to local volumes selection and solutions crosslinking is described. We compared the image quality of the proposed algorithm with results of Filtered Back Projection (FBP) algorithm.
Artifacts caused by intensely absorbing inclusions are encountered in computed tomography via polychromatic scanning and may obscure or simulate pathologies in medical applications. Тo improve the quality of reconstruction if high-Z inclusions in presence, previously we proposed and tested with synthetic data an iterative technique with soft penalty mimicking linear inequalities on the photon-starved rays. This note reports a test at the tomographic laboratory set-up at the Institute of Crystallography FSRC “Crystallography and Photonics” RAS in which tomographic scans were successfully made of temporary tooth without inclusion and with Pb inclusion.
Digital X-ray imaging became widely used in science, medicine, non-destructive testing. This allows using modern digital images analysis for automatic information extraction and interpretation. We give short review of scientific applications of machine vision in scientific X-ray imaging and microtomography, including image processing, feature detection and extraction, images compression to increase camera throughput, microtomography reconstruction, visualization and setup adjustment.
Motion blur caused by camera vibration is a common source of degradation in photographs. In this paper we study the problem of finding the point spread function (PSF) of a blurred image using the tomography technique. The PSF reconstruction result strongly depends on the particular tomography technique used. We present a tomography algorithm with regularization adapted specifically for this task. We use the algebraic reconstruction technique (ART algorithm) as the starting algorithm and introduce regularization. We use the conjugate gradient method for numerical implementation of the proposed approach. The algorithm is tested using a dataset which contains 9 kernels extracted from real photographs by the Adobe corporation where the point spread function is known. We also investigate influence of noise on the quality of image reconstruction and investigate how the number of projections influence the magnitude change of the reconstruction error.
In this paper we propose a novel method for vanishing points detection based on convolutional neural network (CNN) approach and fast Hough transform algorithm. We show how to determine fast Hough transform neural network layer and how to use it in order to increase usability of the neural network approach to the vanishing point detection task. Our algorithm includes CNN with consequence of convolutional and fast Hough transform layers. We are building estimator for distribution of possible vanishing points in the image. This distribution can be used to find candidates of vanishing point. We provide experimental results from tests of suggested method using images collected from videos of road trips. Our approach shows stable result on test images with different projective distortions and noise. Described approach can be effectively implemented for mobile GPU and CPU.
The presence of errors in tomographic image may lead to misdiagnosis when computed tomography (CT) is used in medicine, or the wrong decision about parameters of technological processes when CT is used in the industrial applications. Two main reasons produce these errors. First, the errors occur on the step corresponding to the measurement, e.g. incorrect calibration and estimation of geometric parameters of the set-up. The second reason is the nature of the tomography reconstruction step. At the stage a mathematical model to calculate the projection data is created. Applied optimization and regularization methods along with their numerical implementations of the method chosen have their own specific errors. Nowadays, a lot of research teams try to analyze these errors and construct the relations between error sources. In this paper, we do not analyze the nature of the final error, but present a new approach for the calculation of its distribution in the reconstructed volume. We hope that the visualization of the error distribution will allow experts to clarify the medical report impression or expert summary given by them after analyzing of CT results. To illustrate the efficiency of the proposed approach we present both the simulation and real data processing results.
Obtaining high quality images from Computed Tomography (CT) is important for correct image interpretation. In this paper, we propose novel procedures that can be used for a quantitative description of the degree of artifact expressiveness in CT images, and show that the use of this type of metric allows to assess the dynamics of image degradation. We perform different image reconstruction tests in order to analyse our approach, and the obtained results confirm the usefulness of the proposed method. We conclude that the use of the proposed estimates allows moving from image quality assessment based on visual scoring to a quantitative approach and consequently to support a CT setup providing high quality reconstructed images obtained by appropriate changes of the reconstruction parameters or algorithms.
KEYWORDS: Video, Detection and tracking algorithms, Mobile devices, Image filtering, Machine vision, Image quality, Sensors, Current controlled current source, Patents, Internet
In this paper we consider a task of finding information fields within document with flexible form for credit card expiration date field as example. We discuss main difficulties and suggest possible solutions. In our case this task is to be solved on mobile devices therefore computational complexity has to be as low as possible. In this paper we provide results of the analysis of suggested algorithm. Error distribution of the recognition system shows that suggested algorithm solves the task with required accuracy.
The artifacts (known as metal-like artifacts) arising from incorrect reconstruction may obscure or simulate pathology in medical applications, hide or mimic cracks and cavities in the scanned objects in industrial tomographic scans. One of the main reasons caused such artifacts is photon starvation on the rays which go through highly absorbing regions. We indroduce a way to suppress such artifacts in the reconstructions using soft penalty mimicing linear inequalities on the photon starved rays. An efficient algorithm to use such information is provided and the effect of those inequalities on the reconstruction quality is studied.
Access to the requested content is limited to institutions that have purchased or subscribe to SPIE eBooks.
You are receiving this notice because your organization may not have SPIE eBooks access.*
*Shibboleth/Open Athens users─please
sign in
to access your institution's subscriptions.
To obtain this item, you may purchase the complete book in print or electronic format on
SPIE.org.
INSTITUTIONAL Select your institution to access the SPIE Digital Library.
PERSONAL Sign in with your SPIE account to access your personal subscriptions or to use specific features such as save to my library, sign up for alerts, save searches, etc.