Radon Transport from Soil to Air and Monte-Carlo Simulation

The exhalation of geochemical entities from soil to air is significant to understand Lithosphere-Atmospheric relationships. Some of these geochemical entities are capable of modifying the lower atmosphere, and they are employed in various studies. Radon is one of the geochemical gasses widely recognized as a dominant ionization source in near ground regions of the troposphere. The steady state Rn transport equation is considered in many cases for estimating Rn migration from soil to air on the condition that the time evolution is ignored. A method is proposed for estimating radon space-time transport from soil to air. This is achieved by solving the radon transport equation in soil with special boundary conditions. Similar results are obtained with some experimented models, as well as reported radon values in literature for some set of parameter combinations. Strengths and limitations of the method are discussed. The model is useable to study Lithosphere-Atmosphere relationships. It can also be significant in other studies like the Global Electric Circuit or Seismo-Ionospheric studies.

Change detection in sparse repeat CT scans with non-rigid deformations

BACKGROUND: Detecting and interpreting changes in the images of follow-up CT scans by the clinicians is often time-consuming and error-prone due to changes in patient position and non-rigid anatomy deformations. Thus, reconstructed repeat scan images are required, precluding reduced dose sparse-view repeat scanning. OBJECTIVE: To develop a method to automatically detect changes in a region of interest of sparse-view repeat CT scans in the presence of non-rigid deformations of the patient’s anatomy without reconstructing the original images. METHODS: The proposed method uses the sparse sinogram data of two CT scans to distinguish between genuine changes in the repeat scan and differences due to non-rigid anatomic deformations. First, size and contrast level of the changed regions are estimated from the difference between the scans’ sinogram data. The estimated types of changes in the repeat scan help optimize the method’s parameter values. Two scans are then aligned using Radon space non-rigid registration. Rays which crossed changes in the ROI are detected and back-projected onto image space in a two-phase procedure. These rays form a likelihood map from which the binary changed region map is computed. RESULTS: Experimental studies on four pairs of clinical lung and liver CT scans with simulated changed regions yield a mean changed region recall rate > 86%and a mean precision rate > 83%when detecting large changes with low contrast, and high contrast changes, even when small. The new method outperforms image space methods using prior image constrained compressed sensing (PICCS) reconstruction, particularly for small, low contrast changes (recall = 15.8%, precision = 94.7%). CONCLUSION: Our method for automatic change detection in sparse-view repeat CT scans with non-rigid deformations may assist radiologists by highlighting the changed regions and may obviate the need for a high-quality repeat scan image when no changes are detected.

Deformable registration and region-of-interest image reconstruction in sparse repeat CT scanning

BACKGROUND: Repeat CT scanning is ubiquitous in many clinical situations, e.g. to follow disease progression, to evaluate treatment efficacy, and to monitor interventional CT procedures. However, it incurs in cumulative radiation to the patient which can be significantly reduced by using a region of interest (ROI) and the existing baseline scan. OBJECTIVE: To obtain a high-quality reconstruction of a ROI with a significantly reduced X-ray radiation dosage that accounts for deformations. METHODS: We present a new method for deformable registration and image reconstruction inside an ROI in repeat CT scans with a highly reduced X-ray radiation dose based on sparse scanning. Our method uses the existing baseline scan data, a user-defined ROI, and a new sparse repeat scan to compute a high-quality repeat scan ROI image with a significantly reduced radiation dose. Our method first performs rigid registration between the densely scanned baseline and the sparsely scanned repeat CT scans followed by deformable registration with a low-order parametric model, both in 3D Radon space and without reconstructing the repeat scan image. It then reconstructs the repeat scan ROI without computing the entire repeat scan image. RESULTS: Our experimental results on clinical lung and liver CT scans yield a mean × 14 computation speedup and a × 7.6-12.5 radiation dose reduction, with a minor image quality loss of 0.0157 in the NRMSE metric. CONCLUSION: Our method is considerably faster than existing methods, thereby enabling intraoperative online repeat scanning that it is accurate and accounts for position, deformation, and structure changes at a fraction of the radiation dose required by existing methods.

Radon-Augmented Sentinel-2 Satellite Imagery to Derive Wave-Patterns and Regional Bathymetry

Climatological changes occur globally but have local impacts. Increased storminess, sea level rise and more powerful waves are expected to batter the coastal zone more often and more intense. To understand climate change impacts, regional bathymetry information is paramount. A major issue is that the bathymetries are often non-existent or if they do exist, outdated. This sparsity can be overcome by space-borne satellite techniques to derive bathymetry. Sentinel-2 optical imagery is collected continuously and has a revisit-time around a few days depending on the orbital-position around the world. In this work, Sentinel-2 imagery derived wave patterns are extracted using a localized radon transform. A discrete fast-Fourier (DFT) procedure per direction in Radon space (sinogram) is then applied to derive wave spectra. Sentinel-2 time-lag between detector bands is employed to compute the spectral wave-phase shift and depth using the gravity wave linear dispersion. With this novel technique, regional bathymetries are derived at the test-site of Capbreton, France with an root mean squared (RMS)-error of 2.58 m and a correlation coefficient of 0.82 when compared to the survey for depths until 30 m. With the proposed method, the 10 m Sentinel-2 resolution is sufficient to adequately estimate bathymetries for a wave period of 6.5 s or greater. For shorter periods, the pixel resolution does not allow to detect a stable celerity. In addition to the wave-signature enhancement, the capability of the Radon Transform to augment Sentinel-2 20 m resolution imagery to 10 m is demonstrated, increasing the number of suitable bands for the depth inversion.

Multiscale Spectral Methods for Data Indexing in a Radon Space

This paper presents three methods for multiscale indexing of the content of projection data of computed tomography images in the CBIR-systems for medical database search. The feature spaces in the developed algorithms have been created by means of, respectively: Discrete Stationary Wavelet Transform (DSWT), Shearlet Transform (ST) and Repagulum Wavelet Transform (RWT). A comparative analysis and assessment of the proposed algorithms have been carried out based on experimental studies with computed tomography images.