The pipe wall thickness was estimated based on three-dimensional images of the pipe recovered from several X-ray projections, which were made in a limited angle of view. Since the effects of scattered radiation and beam hardening are up to 50 % of the main radiation, ignoring them leads to blur of the image and inaccuracy in determining dimensions. To restore pipe images from projections, a volume and/or shell representation of the pipe is used, as well as iterative Bayesian methods. Using these methods, the error in estimating the pipe wall thickness from the projection data can be equal to or less than 300 μm. It has been shown that standard X-ray projections on the film or imaging plates used to obtain data can be used to restore pipe wall thickness profiles in factory conditions.
The results of a study of a biaxial stress-strain state using neutron diffraction, standardized in the field of stress measurement, and the Barkhausen effect method are presented. The possibility of using a device with a biaxial load as a standard of stress is evaluated. The article gives an assessment of an error in measuring voltage using strain gauges in relation to measurements made by the neutron diffraction method. The Barkhausen effect method is used as one of the non-destructive testing methods that are sensitive to the stress-strain state. The studies describe the possibility of adapting this method to the proposed stress standard.
OBJECTIVE:To compare image quality and breast density of two reconstruction methods, the widely-used filtered-back projection (FBP) reconstruction and the iterative heuristic Bayesian inference reconstruction (Bayesian inference reconstruction plus the method of total variation applied, HBI). METHODS:Thirty-two clinical DBT data sets with malignant and benign findings, n = 27 and 17, respectively, were reconstructed using FBP and HBI. Three experienced radiologists evaluated the images independently using a 5-point visual grading scale and classified breast density according to the American College of Radiology Breast Imaging-Reporting And Data System Atlas, fifth edition. Image quality metrics included lesion conspicuity, clarity of lesion borders and spicules, noise level, artifacts surrounding the lesion, visibility of parenchyma and breast density. RESULTS:For masses, the image quality of HBI reconstructions was superior to that of FBP in terms of conspicuity,clarity of lesion borders and spicules (p < 0.01). HBI and FBP were not significantly different in calcification conspicuity. Overall, HBI reduced noise and supressed artifacts surrounding the lesions better (p < 0.01). The visibility of fibroglandular parenchyma increased using the HBI method (p < 0.01). On average, five cases per radiologist were downgraded from BI-RADS breast density category C/D to A/B. CONCLUSION:HBI significantly improves lesion visibility compared to FBP. HBI-visibility of breast parenchyma increased, leading to a lower breast density rating. Applying the HBIR algorithm should improve the diagnostic performance of DBT and decrease the need for additional imaging in patients with dense breasts. ADVANCES IN KNOWLEDGE:Iterative heuristic Bayesian inference (HBI) image reconstruction substantially improves the image quality of breast tomosynthesis leading to a better visibility of breast carcinomas and reduction of the perceived breast density compared to the widely-used filtered-back projection (FPB) reconstruction. Applying HBI should improve the accuracy of breast tomosynthesis and reduce the number of unnecessary breast biopsies. It may also reduce the radiation dose for the patients, which is especially important in the screening context.
Experimental investigation of bi-axial steel deformation is presented using new equipment for bi-axial loading. Main conclusion obtained from displayed results is the condition of the invariance of Barkhausen noise intensity relative to any changes of spherical (isotropic) strain/stress tensor, the last being influenced only by material microstructure. This conclusion is supported by independent measurement of Barkhausen noise (BN) intensity on cross-shaped specimens of optimized shape using novel biaxial loading equipment, and tubular specimens like pipe and a balloon loaded by means of oil pressure. Finite Element Modeling (FEM) simulation was also investigated. Thus the BN intensity depends only on the deviatoric (shear) stress tensor value. The presence of this symmetric effect yields too many uncertainties in strain/stress evaluation via BN. Further investigation should identify whether this condition is also present in other magnetic parameters, such as coercive force, remanence, permeability, associated with BN.
Одни из наиболее важных вычислительных задач разработка и адаптация итерационных методов для решения сверхбольших разреженных систем алгебраических уравнений. К таким вычислительным задачам приводит задача итерационной параллельной реконструкции трехмерных изображений промышленных изделий. Важно, что итерационные методы решения вычислительных задач большой размерности реализуются на параллельных структурах намного эффективнее, чем прямые методы их решения. В этой работе описан синхронный параллельный алгоритм, основанный на использовании системы MPI для решения задачи реконструкции трехмерных изображений промышленных изделий. Purpose. Currently, one of the most important tasks is the development and adaptation of iterative methods for solving ultra-large sparse systems of algebraic equations. Such computational problems are caused by the iterative parallel reconstruction of three-dimensional images of industrial products. It is important that iterative methods for solving computational problems of large size are implemented on parallel structures much more efficiently than direct methods for solving them. This paper describes a synchronous parallel algorithm based on the MPI system for solving the problem of reconstruction of three-dimensional images of industrial products. Methodology. It is important that iterative methods for solving computational problems of larger dimensionality on parallel structures are implemented more efficiently than the direct ones. Most algorithms based on direct solving methods have a significant hereditary sequential structure and require a large number of processors interactions which cannot be executed in parallel mode. Iterative methods, for the most part, require a significantly smaller number of interactions of this type and are relatively easily mapped onto parallel computational structures. Equally important, in most cases parallel implementations of classical iterative methods are more effective in terms of computational speed. The parallel execution of the algorithm is based on the distribution of the process in some way between various groups of processors. Depending on the interaction method between local processors, two different types of parallel iterative algorithms execution are distinguished: synchronous and asynchronous. In the former case, it is assumed that the processors complete the calculations and exchange all the necessary results before the start of a new iteration. The main disadvantage of synchronous parallel algorithms is that they require synchronization of iterations. This is a very difficult task, especially with large number of processors. In addition, the overall calculation speed is limited by the speed of the slowest processor. At the same time, faster processors spend most of their time in the waiting mode. But, nevertheless, the implementation of these parallel algorithms can be effectively achieved using the MPI standard. Findings. Synchronous parallel computing algorithms and program codes for threedimensional tomographic reconstruction in a conical beam were developed. Program code debugging and numerical calculations were performed on a hybrid cluster based on the OpenPOWER architecture using the MPI system. For designing a parallel threedimensional tomographic reconstruction, a voxel form of parallelism was used. Originality. Implemented parallel iterative technology of three-dimensional images reconstruction has undeniable advantages over traditional sequential iterative tomographic reconstruction. It allows reducing the time of tomographic reconstruction as many as tens of times, provides the ability to reconstruct products with sizes of 5123 to 10243 voxels with simultaneous storage of the submatrix node of the projection matrix in RAM, which eliminates the need for recalculation of matrix coefficients for each new iteration.
Now one of the most important tasks is development and adaptation of iterative methods for the solution of the big rarefied systems of the algebraic equations. The problem of iterative parallel reconstruction of three-dimensional images of industrial facilities leads to such computing tasks. The fact that iterative methods of the solution of computing problems of big dimension are implemented on parallel structures much more effectively is important than direct methods of their decision. In this work the synchronous parallel algorithm is described, based on use of the MPI system for the solution of a problem of reconstruction of three-dimensional images of industrial facilities.
The technology of three-dimensional Bayesian tomographic reconstruction of homogeneous objects with high-density inclusions is developed. The approach is based on preliminary correction of projections by extracting the data corresponding to X-rays passing through a high-density region, and replacing it with synthesized data obtained by two-dimensional interpolation. An original method for selecting interpolation points is proposed and a mathematical algorithm is described that ensures the implementation of two-dimensional interpolation correction of projections.
The key problem in increasing efficiency of radiotherapy of malignant tumors in brain and other dangerous neoplasms is the problem of increasing the quality of 3D positioning of a patient before radiotherapy. We consider the principles of development of fast parallel iterative algorithms based on graphic accelerators and the OpenGL library. The proposed approaches provide simultaneous residual minimization for the sought solution and total variation of the reconstructed 3D image. In this case, the number of required initial data, i.e., conic X-ray projections, can be reduced several times, and therefore, the radiation load on the patient can also be accordingly reduced with preservation of the necessary contrast and spatial resolution of the 3D image of the patient. The new heuristic iterative algorithm can be used as an alternative to the known 3D Feldkamp algorithm.
The problem of accident prevention for complex structures is considered as multifactor forecasting problem, covering the tasks of Structural Health Monitoring, Non-Destructive Testing, situational modelling for joint behaviour of the object and multi sensor system in extreme cases, risk and residual life time evaluation. This task is considered as an inverse problem being solved in the Bayesian formulation. We introduce three principal characteristics of object behaviour, namely: impact, risk of direct consequences (vulnerability), and risk of indirect consequences. Interaction of the object under observation and monitoring system respectively is displayed in the form of graph, described with a matrix of transition probabilities, given opportunity to evaluate condition probability. The example is given for the calculation of residual life time for the steel cable roof of large sport building, basing on the estimate of operating reliability as a probability integral of critical parameters of detected defects.
The further development of the new iterative reconstruction algorithms to improve three-dimensional breast images quality restored from incomplete and noisy mammograms, is provided. The algebraic reconstruction method with simultaneous iterations - Simultaneous Algebraic Reconstruction Technique (SART) and the iterative method of statistical reconstruction Bayesian Iterative Reconstruction (BIR) are referred here as the preferable iterative methods suitable to improve the image quality. For better processing we use the Graphics Processing Unit (GPU). Method of minimizing the Total Variation (TV) is used as a priori support for regularization of iteration process and to reduce the level of noise in the reconstructed image. Preliminary results with physical phantoms show that all examined methods are capable to reconstruct structures layer-by-layer and to separate layers which images are overlapped in the Z-direction. It was found that the method of traditional Shift-And-Add tomosynthesis (SAA) is worse than iterative methods SART and BIR in terms of suppression of the anatomical noise and image blurring in between the adjacent layers. Despite of the fact that the measured contrast/noise ratio in the presence of low contrast internal structures is higher for the method of tomosynthesis SAA than for SART and BIR methods, its effectiveness in the presence of structured background is rather poor. In our opinion the optimal results can be achieved using Bayesian iterative reconstruction BIR.
The magnetization of ferromagnetic material is considered as periodically inhomogeneous Markov process. The theory assumes both statistically independent and correlated Barkhausen discontinuities. The model, based on the chain evolution-type process theory, assumes that the domain structure of a ferromagnet passes successively the steps of: linear growing, exponential acceleration and domains annihilation to zero density at magnetic saturation. The solution of stochastic differential Kolmogorov equation enables the hysteresis loop calculus.
New parallel iteration algorithms that provide real-time reconstruction of the 3D breast images restored from an incomplete set of noisy mammograms are studied. The simultaneous algebraic reconstruction technique (SART) and Bayesian inference reconstruction (BIR) are considered as advantageous iteration methods that are most suitable for improving the quality of the reconstructed 3D images. The graphics processing unit (GPU) is used to accelerate the reconstruction. The minimization of total variation (TV) is used as a priori support for the regularization of the iteration process and decrease of the noise level in the reconstructed images. Preliminary results for medical physical phantoms show that all the methods are sufficient for the layer-by-layer reconstruction of medical model objects and separation of layers whose images are overlapped on a mammogram that corresponds to vertical transmission (direction along the OZ axis). The traditional shift-and-add (SAA) tomosynthesis is established to be less efficient than SART and BIR in terms of the anatomical-noise reduction and blurring of reconstructed 3D images between conjugate layers. Despite the fact that the estimated contrast-noise ratio, given internal structures with low contrast, is higher for SAA as compared to SART and BIR, its efficiency is very low given the highly structured background. In our opinion, optimal results can be achieved using BIR.
The technique for biaxial stress state quantitative non destructive testing using magnetic, namely Barkhausen Noise, measurements is developed and checked experimentally. The main elaboration concerns the application of uni-axial calibration data for bi-axial stress measurement in the material which treatment pre-history is not definitely known. The article is aimed to get over difficulties, accompanying factual nondestructive stress evaluation, implied from its tensor nature. The developed technique of stress calibration and measurement assumes the bi-axial stress components recovery from uni-axial magnetic and Barkhausen noise measurement results. The complete technology, based on new calibration procedure with grid diagrams is considered in the article.
Conventional X-ray tomography is an effective tool of medical diagnostics, nondestructive testing of engineering structures, and technical diagnostics of manufactured products. Tomographic reconstruction of objects is a strongly incorrect inverse problem in the case of limited vision angle, small number of projections, and/or insufficient X-ray source energy. Hull and hull-voxel methods of tomographic reconstruction of piecewise-uniform objects are considered in this paper. These approaches are quite effective for image reconstruction of industrial and biological objects under conditions of sparse input data and for analyzing binary and piecewise-uniform images. It will be a priori assumed that the object being reconstructed consists of areas with constant attenuation coefficients.
The effect of repeated elasto-plastic deformation by tension of steel specimens on the change of residual stress rearrangement in a crack domain is investigated experimentally and analytically. Reloading can lead to appearance of residual compressive stress in a crack tip, reduce stress concentration and the probability of crack opening when tensile stress is applied. The stress intensity factor calculation should take into consideration the wide stress range in a crack tip, resulting from its deformation pre-history.
A Bayesian iterative method can be the basis for a wide range of technologies in the field of pattern recognition and image reconstruction. It involves finding the most probable solutions for images or patterns, if functionals describing the likelihood function and a priori information, respectively, are already known. The article describes the basic principles and recent advances in the development of BIM and its applications in various fields, mainly in tomography and restoration of functions from incomplete and noisy data.