Features such as particles, pores, or cracks are challenging to measure accurately in CT data when they are small relative to the data resolution, characterized as a point-spread function (PSF). These challenges are particularly acute when paired with segmentation, as the PSF distributes some of the signal from a voxel among neighboring ones; effectively dispersing some of the signal from a given object to a region outside of it. Any feature of interest with one or more dimensions on the order of the PSF will be impacted by this effect, and measurements based on global thresholds necessarily fail. Measurements of the same features should be consistent across different instruments and data resolutions. The PVB (partial volume and blurring) method successfully compensates by quantifying features that are small in all three dimensions based on their attenuation anomaly. By calibrating the CT number of the phase of interest (in this case, gold) it is possible to accurately measure particles down to <6 voxels in data acquired on two instruments, 14 years apart, despite severe artifacts. Altogether, the PVB method is accurate, reproducible, resolution-invariant, and objective; it is also notable for its favorable error structure. The principal challenge is the need for representative effective CT numbers, which reflect not only the features of interest themselves, but also the X-ray spectrum, the size, shape and composition of the enclosing sample, and processing details such as beam-hardening correction. Empirical calibration is the most effective approach.
more »
« less
Accurate Measurement of Small Features in X‐Ray CT Data Volumes, Demonstrated Using Gold Grains
Abstract We present a method for measuring small, discrete features near the resolution limit of X‐ray computed tomography (CT) data volumes with the aim of providing consistent answers across instruments and data resolutions. The appearances of small features are impacted by the partial volume effect and blurring due to the data point‐spread function, and we call our approach the partial‐volume and blurring (PVB) method. Features are segmented to encompass their total attenuation signal, which is then converted to a volume, in turn allowing a subset of voxels to be used to measure shape and orientation. We demonstrate the method on a set of gold grains, scanned with two instruments at a range of resolutions and with various surrounding media. We recover volume accurately over a factor of 27 range in grain volume and factor of 5 range in data resolution, successfully characterizing particles as small as 5.4 voxels in true volume. Shape metrics are affected variably by resolution effects and are more reliable when based on image‐based caliper measurements than perimeter length or surface area. Orientations are reproducible when maximum or minimum axis lengths are sufficiently different from the intermediate axis. Calibration requires end‐member CT numbers for the materials of interest, which we obtained empirically; we describe a first‐principles calculation and discuss its challenges. The PVB method is accurate, reproducible, resolution invariant, and objective, all important improvements over any method based on global thresholds.
more »
« less
- Award ID(s):
- 1762458
- PAR ID:
- 10453999
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 124
- Issue:
- 4
- ISSN:
- 2169-9313
- Page Range / eLocation ID:
- p. 3508-3529
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Optimal surface segmentation is a state-of-the-art method used for segmentation of multiple globally optimal surfaces in volumetric datasets. The method is widely used in numerous medical image segmentation applications. However, nodes in the graph based optimal surface segmentation method typically encode uniformly distributed orthogonal voxels of the volume. Thus the segmentation cannot attain an accuracy greater than a single unit voxel, i.e. the distance between two adjoining nodes in graph space. Segmentation accuracy higher than a unit voxel is achievable by exploiting partial volume information in the voxels which shall result in non-equidistant spacing between adjoining graph nodes. This paper reports a generalized graph based multiple surface segmentation method with convex priors which can optimally segment the target surfaces in an irregularly sampled space. The proposed method allows non-equidistant spacing between the adjoining graph nodes to achieve subvoxel segmentation accuracy by utilizing the partial volume information in the voxels. The partial volume information in the voxels is exploited by computing a displacement field from the original volume data to identify the subvoxel-accurate centers within each voxel resulting in non-equidistant spacing between the adjoining graph nodes. The smoothness of each surface modeled as a convex constraint governs the connectivity and regularity of the surface. We employ an edge-based graph representation to incorporate the necessary constraints and the globally optimal solution is obtained by computing a minimum s-t cut. The proposed method was validated on 10 intravascular multi-frame ultrasound image datasets for subvoxel segmentation accuracy. In all cases, the approach yielded highly accurate results. Our approach can be readily extended to higher-dimensional segmentations.more » « less
-
Superresolution of GOES-16 ABI Bands to a Common High Resolution with a Convolutional Neural NetworkAbstract Superresolution is the general task of artificially increasing the spatial resolution of an image. The recent surge in machine learning (ML) research has yielded many promising ML-based approaches for performing single-image superresolution including applications to satellite remote sensing. We develop a convolutional neural network (CNN) to superresolve the 1- and 2-km bands on the GOES-R series Advanced Baseline Imager (ABI) to a common high resolution of 0.5 km. Access to 0.5-km imagery from ABI band 2 enables the CNN to realistically sharpen lower-resolution bands without significant blurring. We first train the CNN on a proxy task, which allows us to only use ABI imagery, namely, degrading the resolution of ABI bands and training the CNN to restore the original imagery. Comparisons at reduced resolution and at full resolution withLandsat-8/Landsat-9observations illustrate that the CNN produces images with realistic high-frequency detail that is not present in a bicubic interpolation baseline. Estimating all ABI bands at 0.5-km resolution allows for more easily combining information across bands without reconciling differences in spatial resolution. However, more analysis is needed to determine impacts on derived products or multispectral imagery that use superresolved bands. This approach is extensible to other remote sensing instruments that have bands with different spatial resolutions and requires only a small amount of data and knowledge of each channel’s modulation transfer function. Significance StatementSatellite remote sensing instruments often have bands with different spatial resolutions. This work shows that we can artificially increase the resolution of some lower-resolution bands by taking advantage of the texture of higher-resolution bands on theGOES-16ABI instrument using a convolutional neural network. This may help reconcile differences in spatial resolution when combining information across bands, but future analysis is needed to precisely determine impacts on derived products that might use superresolved bands.more » « less
-
Abstract Magnetic resonance imaging (MRI) is a technique that scans the anatomical structure of the brain, whereas functional magnetic resonance imaging (fMRI) uses the same basic principles of atomic physics as MRI scans but image metabolic function. A major goal of MRI and fMRI study is to precisely delineate various types of tissues, anatomical structure, pathologies, and detect the brain regions that react to outer stimuli (e.g., viewing an image). As a key feature of these MRI‐based neuroimaging data, voxels (cubic pixels of the brain volume) are highly correlated. However, the associations between voxels are often overlooked in the statistical analysis. We adapt a recently proposed dimension reduction method called the envelope method to analyze neuoimaging data taking into account correlation among voxels. We refer to the modified procedure the envelope chain procedure. Because the envelope chain procedure has not been employed before, we demonstrate in simulations the empirical performance of estimator, and examine its sensitivity when our assumptions are violated. We use the estimator to analyze the MRI data from ADHD‐200 study. Data analyses demonstrate that leveraging the correlations among voxels can significantly increase the efficiency of the regression analysis, thus achieving higher detection power with small sample sizes.more » « less
-
In image-based finite element analysis of bone, partial volume effects (PVEs) arise from image blur at tissue boundaries and as a byproduct of geometric reconstruction and meshing during model creation. In this study, we developed and validated a material assignment approach to mitigate partial volume effects. Our validation data consisted of physical torsion testing of intact tibiae from N = 20 Swiss alpine sheep. We created finite element models from micro-CT scans of these tibiae using three popular element types (10-node tetrahedral, 8-node hexahedral, and 20-node hexahedral). Without partial volume management, the models over-predicted the torsional rigidity compared to physical biomechanical tests. To address this problem, we implemented a dual-zone material model to treat elements that overlap low-density surface voxels as soft tissue rather than bone. After in situ inverse optimization, the dual-zone material model produced strong correlations and high absolute agreement between the virtual and physical tests. This suggests that with appropriate partial volume management, virtual mechanical testing can be a reliable surrogate for physical biomechanical testing. For maximum flexibility in partial volume management regardless of element type, we recommend the use of the following dual-zone material model for ovine tibiae: soft-tissue cutoff density of 665 mgHA/cm3 with a soft tissue modulus of 50 MPa (below cutoff) and a density-modulus conversion slope of 10,225 MPa-cm3/mgHA for bone (above cutoff).more » « less
An official website of the United States government
