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.

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.

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.

Virtual unrolling or unfolding, digital unwrapping, flattening or unfurling - all these terms are used to describe the process of surface straightening of a tomographically reconstructed digital object. For many objects of historical heritage, tomography is the only way to obtain a hidden image of the original object without its destruction. Digital flattening is no longer considered a unique met hodology. It being applied by many research group, but AI-based methods are used insignificantly in such projects, despite the amazing success of AI in computer vision, in particular optical text recognition. It can be explained by the fact that the success of AI depends on large, broad and high quality datasets, but there are very few published CT-based datasets relevant to the task of digital flattening. Accumulation of a sufficient amount of data necessary for training models is a key point for the next technological breakthrough. In this paper, we present open and cumulative dataset CT-OCR-2022. Dataset includes 6 packages data for different model objects that help to enrich tomographic solutions and to train machine learning models. Each package contains optically scanned image of model objects, 400 measured X-ray projections, 2687 CT- reconstructed cross-sections of 3D reconstructed image, segmentation markups. We believe that CT-OCR-2022 dataset will serve as a benchmark for reconstructed object digital flattening and recognition systems, and that it will prove invaluable for advancement of the field of CT-reconstruction, symbols analysis and recognition. The data presented are openly available in Zenodo at doi:10.5281/zenodo.7123495 and linked repositories.

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 self­driving 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 X­ray 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.

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 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.

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.

In this paper we introduce a novel neural network architecture based on Fast Hough Transform layer. The layer of this type allows our neural network to accumulate features from linear areas across the entire image instead of local areas. We demonstrate its potential by solving the problem of vanishing points detection in the images of documents. Such problem occurs when dealing with camera shots of the documents in uncontrolled conditions. In this case, the document image can suffer several specific distortions including projective transform. To train our model, we use MIDV-500 dataset and provide testing results. Strong generalization ability of the suggested method is proven with its applying to a completely different ICDAR 2011 dewarping contest. In previously published papers considering this dataset authors measured quality of vanishing point detection by counting correctly recognized words with open OCR engine Tesseract. To compare with them, we reproduce this experiment and show that our method outperforms the state-of-the-art result.

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.

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.

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.

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.

Artifacts caused by intensely absorbing areas are encountered in computed tomography and may obscure or simulate pathology in medical applications, hide or mimic the cracks and cavities in the devices at industrial applications. We simulated sinograms with different levels of absorption to demonstrate the artifacts dynamics. If the analysis of the measured data shows the presence of strongly absorbing areas in the object under study we propose to use quadratic programming technique for solving the inverse problem. Although the technique is time-consuming it allows us to avoid the typical artifacts. We compare the images reconstructed with different techniques including the proposed one.

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.

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.