Real-time Tomography Research Articles

Distributed travel-time seismic tomography in large-scale sensor networks

Current geophysical techniques for visualizing seismic activity employ image reconstruction methods that rely on a centralized approach for processing the raw data captured by seismic sensors. The data is either gathered manually, or relayed by expensive broadband stations, and then processed at a base station. This approach is time-consuming (weeks to months) and hazardous as the task involves manual data gathering in extreme conditions. Also, raw seismic samples are typically in the range of 16–24 bit, sampled at 50–200 Hz and transferring this high fidelity sample from large number of sensors to a centralized station results in a bottleneck due to bandwidth limitations. To avoid these issues, a new distributed method is required which processes raw seismic samples inside each node and obtains a high-resolution seismic tomography in real time. In this paper, we present a component-averaged distributed multi-resolution evolving tomography algorithm for processing data and inverting volcano tomography in the network while avoiding centralized computation and costly data collection. The algorithm is first evaluated for the correctness using a synthetic model in a CORE emulator. Later, our proposed algorithm runs using the real data obtained from Mt. St. Helens, WA, USA. The results validate that our distributed algorithm is able to obtain a satisfactory image similar to centralized computation under constraints of network resources, while distributing the computational burden to sensor nodes.

Decentralized multigrid for in-situ big data computing

Modern seismic sensors are capable of recording high precision vibration data continuously for several months. Seismic raw data consists of information regarding earthquake's origin time, location, wave velocity, etc. Currently, these high volume data are gathered manually from each station for analysis. This process restricts us from obtaining high-resolution images in real-time. A new in-network distributed method is required that can obtain a high-resolution seismic tomography in real time. In this paper, we present a distributed multigrid solution to reconstruct seismic image over large dense networks. The algorithm performs in-network computation on large seismic samples and avoids expensive data collection and centralized computation. Our evaluation using synthetic data shows that the proposed method accelerates the convergence and reduces the number of messages exchanged. The distributed scheme balances the computation load and is also tolerant to severe packet loss.

High-definition optical coherence tomography intrinsic skin ageing assessment in women: a pilot study.

Several non-invasive two-dimensional techniques with different lateral resolution and measurable depth range have proved to be useful in assessing and quantifying morphological changes in skin ageing. Among these, only in vivo microscopy techniques permit histometric measurements in vivo. Qualitative and quantitative assessment of chronological (intrinsic) age-related (IAR) morphological changes of epidermis, dermo-epidermal junction (DEJ), papillary dermis (PD), papillary-reticular dermis junction and reticular dermis (RD) have been performed by high-definition optical coherence tomography in real time 3-D. HD-OCT images were taken at the internal site of the right upper arm. Qualitative HD-OCT IAR descriptors were reported at skin surface, at epidermal layer, DEJ, PD and upper RD. Quantitative evaluation of age-related compaction and backscattered intensity or brightness of different skin layers was performed by using the plugin plot z-axis profile of ImageJ® software permitting intensity assessment of HD-OCT (DICOM) images (3-D images). Analysis was in blind from all clinical information. Sixty, fair-skinned (Fitzpatrick types I–III) healthy females were analysed retrospectively in this study. The subjects belonged to three age groups: twenty in group I aged 20–39, twenty in group II aged 40–59 and twenty in group III aged 60–79. Only intrinsic ageing in women has been studied. Significant age-related qualitative and quantitative differences could be noticed. IAR changes in dermal matrix fibers morphology/organisation and in microvasculature were observed. The brightness and compaction of the different skin layers increased significantly with intrinsic skin ageing. The depth of visibility of fibers in RD increased significantly in the older age group. In conclusion, HD-OCT allows 3-D in vivo and real time qualitative and quantitative assessment of chronological (intrinsic) age-related morphological skin changes at high resolution from skin surface to a depth of the superficial reticular dermis.

High speed imaging of dynamic processes with a switched source x-ray CT system

Conventional x-ray computed tomography (CT) scanners are limited in their scanning speed by the mechanical constraints of their rotating gantries and as such do not provide the necessary temporal resolution for imaging of fast-moving dynamic processes, such as moving fluid flows. The Real Time Tomography (RTT) system is a family of fast cone beam CT scanners which instead use multiple fixed discrete sources and complete rings of detectors in an offset geometry. We demonstrate the potential of this system for use in the imaging of such high speed dynamic processes and give results using simulated and real experimental data. The unusual scanning geometry results in some challenges in image reconstruction, which are overcome using algebraic iterative reconstruction techniques and explicit regularisation. Through the use of a simple temporal regularisation term and by optimising the source firing pattern, we show that temporal resolution of the system may be increased at the expense of spatial resolution, which may be advantageous in some situations. Results are given showing temporal resolution of approximately 500 µs with simulated data and 3 ms with real experimental data.

Mapping the structure

The demonstration of real-time and non-destructive Doppler-assisted tomography of the internal structure of photonic-crystal fibres could aid the fabrication of high-quality fibres with enhanced performance.

TU-F-18C-05: Evaluation of a Method to Calculate Patient-Oriented MGD Coefficients Using Estimates of Glandular Tissue Distribution

Purpose: Evaluate mammographic mean glandular dose (MGD) coefficients for particular known tissue distributions using a novel formalism that incorporates the effect of the heterogeneous glandular tissue distribution, by comparing them with MGD coefficients derived from the corresponding anthropomorphic computer breast phantom. Methods: MGD coefficients were obtained using MCNP5 simulations with the currently used homogeneous assumption and the heterogeneously-layered breast (HLB) geometry and compared against those from the computer phantom (ground truth). The tissue distribution for the HLB geometry was estimated using glandularity map image pairs corrected for the presence of non-glandular fibrous tissue. Heterogeneity of tissue distribution was quantified using the glandular tissue distribution index, Idist. The phantom had 5 cm compressed breast thickness (MLO and CC views) and 29% whole breast glandular percentage. Results: Differences as high as 116% were found between the MGD coefficients with the homogeneous breast core assumption and those from the corresponding ground truth. Higher differences were found for cases with more heterogeneous distribution of glandular tissue. The Idist for all cases was in the [−0.8{sup −}+0.3] range. The use of the methods presented in this work results in better agreement with ground truth with an improvement as high as 105 pp. The decrease inmore » difference across all phantom cases was in the [9{sup −}105] pp range, dependent on the distribution of glandular tissue and was larger for the cases with the highest Idist values. Conclusion: Our results suggest that the use of corrected glandularity image pairs, as well as the HLB geometry, improves the estimates of MGD conversion coefficients by accounting for the distribution of glandular tissue within the breast. The accuracy of this approach with respect to ground truth is highly dependent on the particular glandular tissue distribution studied. Predrag Bakic discloses current funding from NIH, NSF, and DoD, former funding from Real Time Tomography, LLC and a current research collaboration with Barco and Hologic.« less

Multi-sheet surface rebinning methods for reconstruction from asymmetrically truncated cone beam projections: I. Approximation and optimality

The mechanical motion of the gantry in conventional cone beam CT scanners restricts the speed of data acquisition in applications with near real time requirements. A possible resolution of this problem is to replace the moving source detector assembly with static parts that are electronically activated. An example of such a system is the Rapiscan Systems RTT80 real time tomography scanner, with a static ring of sources and axially offset static cylinder of detectors. A consequence of such a design is asymmetrical axial truncation of the cone beam projections resulting, in the sense of integral geometry, in severely incomplete data. In particular we collect data only in a fraction of the Tam–Danielsson window, hence the standard cone beam reconstruction techniques do not apply. In this work we propose a family of multi-sheet surface rebinning methods for reconstruction from such truncated projections. The proposed methods combine analytical and numerical ideas utilizing linearity of the ray transform to reconstruct data on multi-sheet surfaces, from which the volumetric image is obtained through deconvolution. In this first paper in the series, we discuss the rebinning to multi-sheet surfaces. In particular we concentrate on the underlying transforms on multi-sheet surfaces and their approximation with data collected by offset multi-source scanning geometries like the RTT. The optimal multi-sheet surface and the corresponding rebinning function are found as a solution of a variational problem. In the case of the quadratic objective, the variational problem for the optimal rebinning pair can be solved by a globally convergent iteration. Examples of optimal rebinning pairs are computed for different trajectories. We formulate the axial deconvolution problem for the recovery of the volumetric image from the reconstructions on multi-sheet surfaces. Efficient and stable solution of the deconvolution problem is the subject of the second paper in this series (Betcke and Lionheart 2013 Inverse Problems 29 115004).

Real-time in vivo computed optical interferometric tomography

High-resolution real-time tomography of scattering tissues is important for many areas of medicine and biology1–6. However, the compromise between transverse resolution and depth-of-field in addition to low sensitivity deep in tissue continue to impede progress towards cellular-level volumetric tomography. Computed imaging has the potential to solve these long-standing limitations. Interferometric synthetic aperture microscopy (ISAM)7–9 is a computed imaging technique enabling high-resolution volumetric tomography with spatially invariant resolution. However, its potential for clinical diagnostics remains largely untapped since full volume reconstructions required lengthy postprocessing, and the phase-stability requirements have been difficult to satisfy in vivo. Here we demonstrate how 3-D Fourier-domain resampling, in combination with high-speed optical coherence tomography (OCT), can achieve high-resolution in vivo tomography. Enhanced depth sensitivity was achieved over a depth-of-field extended in real time by more than an order of magnitude. This work lays the foundation for high-speed volumetric cellular-level tomography.

A Survey of Recent MARTe Based Systems

The Multithreaded Application Real-Time executor (MARTe) is a data driven framework environment for the development and deployment of real-time control algorithms. The main ideas which led to the present version of the framework were to standardize the development of real-time control systems, while providing a set of strictly bounded standard interfaces to the outside world and also accommodating a collection of facilities which promote the speed and ease of development, commissioning and deployment of such systems. At the core of every MARTe based application, is a set of independent inter-communicating software blocks, named Generic Application Modules (GAM), orchestrated by a real-time scheduler. The platform independence of its core library provides MARTe the necessary robustness and flexibility for conveniently testing applications in different environments including non-real-time operating systems. MARTe is already being used in several machines, each with its own peculiarities regarding hardware interfacing, supervisory control configuration, operating system and target control application. This paper presents and compares the most recent results of systems using MARTe: the JET Vertical Stabilization system, which uses the Real Time Application Interface (RTAI) operating system on Intel multi-core processors; the COMPASS plasma control system, driven by Linux RT also on Intel multi-core processors; ISTTOK real-time tomography equilibrium reconstruction which shares the same support configuration of COMPASS; JET error field correction coils based on VME, PowerPC and VxWorks; FTU LH reflected power system running on VME, Intel with RTAI.

Data acquisition and real-time bolometer tomography using LabVIEW RT

Abstract The currently available multi-core PCI Express systems running LabVIEW RT (real-time), equipped with FPGA cards for data acquisition and real-time parallel signal processing, greatly shorten the design and implementation cycles of large-scale, real-time data acquisition and control systems. This paper details a data acquisition and real-time tomography system using LabVIEW RT for the bolometer diagnostic on the ASDEX Upgrade tokamak (Max Planck Institute for Plasma Physics, Garching, Germany). The transformation matrix for tomography is pre-computed based on the geometry of distributed radiation sources and sensors. A parallelized iterative algorithm is adapted to solve a constrained linear system for the reconstruction of the radiated power density. Real-time bolometer tomography is performed with LabVIEW RT. Using multi-core machines to execute the parallelized algorithm, a cycle time well below 1 ms is reached.

ChemInform Abstract: From Solid-Liquid Ratios to Real Time Tomography: The Development of NMR in Food Applications

ChemInform Abstract: From Solid-Liquid Ratios to Real Time Tomography: The Development of NMR in Food Applications

Ring artefact suppression in real-time x-ray tomography

A ring artefact suppression algorithm in x-ray tomography is proposed and allows one to process input data in real time. The algorithm is based on methods of the theory of inverse and ill-posed problems. Its numerical implementation uses minimisation of the smoothing Tikhonov’s functional with the conjugate gradient method.

MIRTO: a prototype for real-time ionospheric imaging over the Mediterranean area

MIRTO (Mediterranean Ionosphere with Real-time TOmography) is a collaborative project between Istituto Nazionale di Geofisica (INGV) of Rome, the University of Bath (U.K.) and the Istituto Fisica Applicata «Nello Carrara»-Consiglio Nazionale delle Ricerche (IFAC-CNR) of Florence. The goal of the project is the development of a prototype for real-time imaging of the ionosphere over the Italian region with extension to the Mediterranean Sea. MIRTO uses an original imaging technique developed at the University of Bath and upgraded for real-time use in cooperation with IFAC. The prototype makes use of the data acquired by the real-time ionospheric and geodetic instrumentation operated by INGV. Such measurements drive the imaging algorithm to produce the image of electron density as well as maps and movies of the Total Electron Content (TEC) over the considered area.

Study on concurrent measurements of velocity and density distributions for high-speed CO2 flow

Velocity and density distributions of a high-speed and initial CO$_{2}$ jet flow have been analyzed concurrently by a developed three-dimensional digital speckle tomography and a particle image velocimetry (PIV). Two high-speed cameras have been used for the tomography and one for the PIV since a shape of a nozzle for the jet flow is bilaterally symmetric and the starting flow is fast and unsteady. The speckle movements between no flow and CO$_{2}$ jet flow have been obtained by a cross-correlation tracking method so that those distances can be transferred to deflection angles of laser rays for density gradients. The three-dimensional density fields for the high-speed CO$_{2}$ jet flow have been reconstructed from the deflection angles by the real-time tomography method, and the two-dimensional velocity fields have been calculated by the PIV method concurrently and instantaneously.

Analysis of Simultaneous Velocity and Density Distributions for High-Speed CO<sub>2</sub> Flow Using Particle Image Velocimetry and Digital Speckle Tomography

Velocity and density distributions of a high-speed and initial CO2 jet flow have been analyzed simultaneously by a developed three-dimensional digital speckle tomography and a particle image velocimetry (PIV). Three high-speed cameras have been used for the tomography and the PIV since a shape of a nozzle for the jet flow is asymmetric and the initial flow is fast and unsteady. The speckle movements between no flow and CO2 jet flow have been obtained by a cross-correlation tracking method so that those distances can be transferred to deflection angles of laser rays for density gradients. The three-dimensional density fields for the high-speed CO2 jet flow have been reconstructed from the deflection angles by the real-time tomography method, and the two-dimensional velocity fields have been calculated by the PIV method simultaneously.

Experimental analysis of high-speed helium jet flow using four-dimensional digital speckle tomography

A high-speed and initial helium jet flow has been analyzed by a developed four-dimensional digital speckle tomography. Multiple high-speed cameras have been used to capture movements of speckles in multiple angles of view simultaneously because a shape of a nozzle for the jet flow is asymmetric and the initial jet flow is fast and unsteady. The speckle movements between no flow and helium jet flow from the asymmetric nozzle controlled by a solenoid valve have been obtained by a cross-correlation tracking method so that those distances can be transferred to deflection angles of laser rays for density gradients. The four-dimensional density fields for the high-speed helium jet flow have been reconstructed from the deflection angles by a developed real-time tomography method.

Recursive travel-time inversion: A tool for real-time seismic tomography

Abstract A new recursive inverse scheme is applied to a currently popular problem named seismic travel-time tomography, in order to enhance the efficiency and reliability in obtaining a new velocity model if a small number of new data are added to a large data set in the past. In comparison with conventional inverse schemes in seismic tomography, either least-squares or iterative types, this scheme does not require large amounts of matrix-type computations but utilizes the amount of modification in model parameters responsible for each new data set. We also introduce the computation of a collocation travel time (i.e., from a given station to every grid point) for the reference velocity model inverted by the data for all the past events, using a ray tracing scheme called the Huygens’ method (Saito, 2001), suitable to computations prior to a new event. Combining the above information already stored with the recursive inverse scheme, we can obtain a new or updated velocity model immediately after a new event takes place, because a temporal interval between two events is usually very long in a given local area. Since the model is revised at each recursive step, we perform ray tracings with the updated reference model to get more accurate ray paths and travel times than the conventional inversion schemes that use all the ray tracings for the same reference model. We first showed the validity and stability of the proposed method with synthetic data. We then applied the new approach to the P-wave travel-time data recorded in the Hidaka, south-central Hokkaido, Japan, region, and compared our result with other previous results. Our result shares the overall feature with the previous ones. In addition, a new low-velocity zone is detected in the east of the Hidaka mountains at the depth of 10 km, corresponding to the collision zone of two arcs, due to the use of the updated reference velocity model at each recursive step. We also confirmed that the order of data does not affect the final result, so that the present approach is shown as an appropriate tool for so-called real-time seismic tomography: a updated velocity model is immediately obtained at each time that a new event takes place, in order to monitor temporal variations of model parameters such as velocity structure on the real-time basis.

Present and perspective roles of soft X-ray tomography in tokamak plasma position measurements

This paper shows the importance and feasibility of real-time tomography in fusion experiments for the example of soft X-ray (SXR) position measurements. The requirement of non-magnetic real-time diagnostics in low frequencies for ITER is discussed. This is illustrated by recent results of rapid tomographic inversion of SXR measurements on tokamak TCV. Comparison with the magnetic reconstruction data has not only shown the valuable resolution capabilities of both techniques, but also revealed a slight dependence of magnetic measurements on toroidal magnetic field and an unnoticed drift of plasma position observer. A feasibility study using current hardware capacities for programmable real-time tomographic system with plasma position feedback output was carried out. A compact solution is found to be tractable opening wide perspectives for development.

An algorithm for real-time tomography of gas concentrations, using prior information about spatial derivatives

We present a new computed tomography method, the low third derivative (LTD) method, that is particularly suited for reconstructing the spatial distribution of gas concentrations from path-integral data for a small number of optical paths. The method finds a spatial distribution of gas concentrations that (1) has path integrals that agree with measured path integrals, and (2) has a low third spatial derivative in each direction, at every point. The trade-off between (1) and (2) is controlled by an adjustable parameter, which can be set based on analysis of the path-integral data. The method produces a set of linear equations, which can be solved with a single matrix multiplication if the constraint that all concentrations must be positive is ignored; the method is therefore extremely rapid. Analysis of experimental data from thousands of concentration distributions shows that the method works nearly as well as smooth basis function minimization (the best method previously available), yet is about 100 times faster.

Computational and interferometric system for real-time limited-view tomography of flow fields

A modified iterative conjugate gradient technique is developed and applied to limited-view interferometric tomography for reconstructing flow fields of simple structure. The approach is also effective in dealing with significant noise present in the data. From a computer simulation of experiments, a multipath interferometric system is built, from which multidirectional projection interferograms can be captured in real time for fast transient phenomena. The system is employed to measure an asymmetric temperature field. The three-dimensional field is then reconstructed with the developed experimental and computational techniques. The results are satisfactory when compared with thermocouple measurements.