Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D Visualization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52

On the Origin of Sinter‐Resistance and Catalyst Accessibility in Raspberry‐Colloid‐Templated Catalyst Design

Abstract Nanoparticle (NP) sintering is a major cause of the deactivation of supported catalysts. Raspberry‐colloid‐templated (RCT) catalysts are an emerging class of materials that show an unprecedented level of sinter‐resistance and exhibit high catalytic activity. Here a comprehensive study of the origin of NP stability and accessibility in RCT catalysts using theoretical modeling, 3D electron microscopy, and epitaxial overgrowth is reported. The approach is showcased for silica‐based RCT catalysts containing dilute Pd‐in‐Au NPs previously used in hydrogenation and oxidation catalysis. Modeling of the contact line of the silica precursor infiltrating into the assembled raspberry colloids suggests that a large part of the particles must be embedded into silica, which is confirmed by quantitative visualization of >200 individual NPs by dual‐axis electron tomography. The RCT catalysts have a unique structure in which all NPs reside at the pore wall but have >50% of their surface embedded in the matrix, giving rise to the strongly enhanced thermal and mechanical stability. Importantly, epitaxial overgrowth of Ag on the supported NPs reveals that not only the NP surface exposed to the pore but the embedded interface as well remained chemically accessible. This mechanistic understanding provides valuable guidance in the design of stable catalytic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A measure of active interfaces in supported catalysts for high-temperature reactions

Formulating knowledge of structure-function relationships in heterogeneous catalysis is central to the design of efficient catalysts; yet, the elucidation of dominant reaction sites has remained as a challenge. Here, in this paper, we present a methodology that can be used to visualize metal-gas and metal-oxide-gas interfaces in three dimensions and to quantify their catalytic activity levels. As a case study, CH 4 oxidation occurring in a Pt/CeO 2 system is chosen. By employing thermally robust Pt@CeO 2 model catalysts with size-tunable and monodisperse cores, and gas-permeable shells, we reconstruct a series of structures in 3D via electron tomography and match the information to activity data and theoretical calculations. This strategy reveals that two different interfaces catalyze the CH 4 oxidation and that their contribution to the overall rate changes with the Pt size, temperature, and gas atmosphere. Our results provide an analytic platform on which to explore reaction pathways and mechanisms applicable to multiple reactions and materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Real-time visualization of particle evolution during reactive flux-assisted processing of aluminum melts

Here, a multi-modal, multi-scale correlative tomography investigation of Al-TiC metal matrix composites processed via flux-assisted reaction synthesis is reported. Synchrotron X-ray microradiography is utilized to visualize the reaction and particle evolution in real-time. Changes in particle diameter and areal number density suggest that the process is nucleation- rather than growth-dominated. At 950 °C, the bulk of the reaction takes place in a relatively short time span of less than 600 s. The microstructure is imaged at higher resolution in 2D (scanning electron microscopy) and 3D (synchrotron X-ray nanotomography), revealing the formation of carbide particles with a hexagonal platelet morphology. We propose that the morphology arises due to the incorporation of Si impurities during the experiment. It is expected that the correlative tomography workflow and analysis may guide future metal matrix composite (MMC) processing strategies.

36 MATERIALS SCIENCE↗

Three-dimensional X-ray diffraction imaging of dislocations in polycrystalline metals under tensile loading

Abstract The nucleation and propagation of dislocations is an ubiquitous process that accompanies the plastic deformation of materials. Consequently, following the first visualization of dislocations over 50 years ago with the advent of the first transmission electron microscopes, significant effort has been invested in tailoring material response through defect engineering and control. To accomplish this more effectively, the ability to identify and characterize defect structure and strain following external stimulus is vital. Here, using X-ray Bragg coherent diffraction imaging, we describe the first direct 3D X-ray imaging of the strain field surrounding a line defect within a grain of free-standing nanocrystalline material following tensile loading. By integrating the observed 3D structure into an atomistic model, we show that the measured strain field corresponds to a screw dislocation.

36 MATERIALS SCIENCE↗

SU-E-T-141: Automated Dose Point Placement for Cervical Cancer Brachytherapy Using Tandem and Ovoid Applicators

Purpose: To develop a standalone application, which automatically and consistently calculates the coordinates of points A and H based solely on the implanted applicator geometry for cervical cancer HDR brachytherapy. Methods: Manchester point A and ABS point H are both located 2cm lateral from the central tandem plane. While both points are located 2cm above the cervical os, surrogates for the os differ. Point A is defined relative to the anatomical cervical os. Point H is defined relative to the intersection of the tandem with the superior aspects of the ovoids. The application takes an input text file generated by the treatment planning system (TPS, BrachyVision, Varian) that specifies the source geometries. It then outputs the 3D coordinates of points A and H in both the left and right directions. The algorithm was implemented and tested on 34 CT scans of 7 patients treated with HDR brachytherapy delivered using tandem and ovoids. A single experienced user retrospectively and manually placed points A and H on the CT scans, whose coordinates were used as the gold standard for the comparison to the automatically calculated points. Results: The automatically calculated coordinates of points A and H agree within 0.7mm with the gold standard. The averages and standard deviations of the 3D coordinate difference between points placed by the two methods are 0.3±0.1 and 0.4±0.1mm for points A and H, respectively. The maximum difference in 3D magnitude is 0.7mm. Conclusion: The algorithm consistently calculates dose point coordinates independently of the planner for cervical cancer brachytherapy treated with tandem and ovoids. Automated point placement based on the geometry of the implanted applicators agrees in sub-millimeter with careful manual placements by an experienced user. This algorithm expedites the planning process and eliminates dependencies on either user input or TPS visualization tools.

60 APPLIED LIFE SCIENCES↗

WE-AB-207A-12: HLCC Based Quantitative Evaluation Method of Image Artifact in Dental CBCT

Purpose: Image artifacts are usually evaluated qualitatively via visual observation of the reconstructed images, which is susceptible to subjective factors due to the lack of an objective evaluation criterion. In this work, we propose a Helgason-Ludwig consistency condition (HLCC) based evaluation method to quantify the severity level of different image artifacts in dental CBCT. Methods: Our evaluation method consists of four step: 1) Acquire Cone beam CT(CBCT) projection; 2) Convert 3D CBCT projection to fan-beam projection by extracting its central plane projection; 3) Convert fan-beam projection to parallel-beam projection utilizing sinogram-based rebinning algorithm or detail-based rebinning algorithm; 4) Obtain HLCC profile by integrating parallel-beam projection per view and calculate wave percentage and variance of the HLCC profile, which can be used to describe the severity level of image artifacts. Results: Several sets of dental CBCT projections containing only one type of artifact (i.e. geometry, scatter, beam hardening, lag and noise artifact), were simulated using gDRR, a GPU tool developed for efficient, accurate, and realistic simulation of CBCT Projections. These simulated CBCT projections were used to test our proposed method. HLCC profile wave percentage and variance induced by geometry distortion are about 3∼21 times and 16∼393 times as large as that of the artifact-free projection, respectively. The increase factor of wave percentage and variance are 6 and133 times for beam hardening, 19 and 1184 times for scatter, and 4 and16 times for lag artifacts, respectively. In contrast, for noisy projection the wave percentage, variance and inconsistency level are almost the same with those of the noise-free one. Conclusion: We have proposed a quantitative evaluation method of image artifact based on HLCC theory. According to our simulation results, the severity of different artifact types is found to be in a following order: Scatter>Geometry>Beam hardening>Lag>Noise>Artifact-free in dental CBCT.

60 APPLIED LIFE SCIENCES↗

Examination of Multiple Lithologies Within the Primitive Ordinary Chondrite NWA 5717

Northwest Africa 5717 is a primitive (subtype 3.05) ungrouped ordinary chondrite which contains two apparently distinct lithologies. In large cut meteorite slabs, the darker of these, lithology A, looks to host the second, much lighter in color, lithology B (upper left, Fig. 1). The nature of the boundary between the two is uncertain, ranging from abrupt to gradational and not always following particle boundaries. The distinction between the lithologies, beyond the obvious color differences, has been supported by a discrepancy in oxygen isotopes and an incongruity in the magnesium contents of chondrule olivine. Here, quantitative textural analysis and mineralogical methods have been used to investigate the two apparent lithologies within NWA 5717. Olivine grains contained in a thin section from NWA 7402, thought to be paired to 5717, were also measured to re-examine the distinct compositional range among the light and dark areas. Procedure: Particles from a high-resolution mosaic image of a roughly 13x15cm slice of NWA 5717 were traced in Adobe Photoshop. Due to the large size of the sample, visually representative regions of each lithology were chosen to be analyzed. The resulting layers of digitized particles were imported into ImageJ, which was used to measure their area, along with the axes, the angle from horizontal, and the centroid coordinates of ellipses fitted to each particle following the approach. Resulting 2D pixel areas were converted to spherical diameters employing the unfolding algorithm, which outputs a 3D particle size distribution based on digitized 2D size frequency data. Spatstat was used to create kernel density plots of the centroid coordinates for each region. X-ray compositional maps, microprobe analyses, and Mossbauer spectroscopy was conducted on a thin section of NWA 7402, tentatively paired to NWA 5717.

Cato, M. J.↗

GMI-IPS: Python Processing Software for Aircraft Campaigns

NASA's Atmospheric Tomography Mission (ATom) seeks to understand the impact of anthropogenic air pollution on gases in the Earth's atmosphere. Four flight campaigns are being deployed on a seasonal basis to establish a continuous global-scale data set intended to improve the representation of chemically reactive gases in global atmospheric chemistry models. The Global Modeling Initiative (GMI), is creating chemical transport simulations on a global scale for each of the ATom flight campaigns. To meet the computational demands required to translate the GMI simulation data to grids associated with the flights from the ATom campaigns, the GMI ICARTT Processing Software (GMI-IPS) has been developed and is providing key functionality for data processing and analysis in this ongoing effort. The GMI-IPS is written in Python and provides computational kernels for data interpolation and visualization tasks on GMI simulation data. A key feature of the GMI-IPS, is its ability to read ICARTT files, a text-based file format for airborne instrument data, and extract the required flight information that defines regional and temporal grid parameters associated with an ATom flight. Perhaps most importantly, the GMI-IPS creates ICARTT files containing GMI simulated data, which are used in collaboration with ATom instrument teams and other modeling groups. The initial main task of the GMI-IPS is to interpolate GMI model data to the finer temporal resolution (1-10 seconds) of a given flight. The model data includes basic fields such as temperature and pressure, but the main focus of this effort is to provide species concentrations of chemical gases for ATom flights. The software, which uses parallel computation techniques for data intensive tasks, linearly interpolates each of the model fields to the time resolution of the flight. The temporally interpolated data is then saved to disk, and is used to create additional derived quantities. In order to translate the GMI model data to the spatial grid of the flight path as defined by the pressure, latitude, and longitude points at each flight time record, a weighted average is then calculated from the nearest neighbors in two dimensions (latitude, longitude). Using SciPya's Regular Grid Interpolator, interpolation functions are generated for the GMI model grid and the calculated weighted averages. The flight path points are then extracted from the ATom ICARTT instrument file, and are sent to the multi-dimensional interpolating functions to generate GMI field quantities along the spatial path of the flight. The interpolated field quantities are then written to a ICARTT data file, which is stored for further manipulation. The GMI-IPS is aware of a generic ATom ICARTT header format, containing basic information for all flight campaigns. The GMI-IPS includes logic to edit metadata for the derived field quantities, as well as modify the generic header data such as processing dates and associated instrument files. The ICARTT interpolated data is then appended to the modified header data, and the ICARTT processing is complete for the given flight and ready for collaboration. The output ICARTT data adheres to the ICARTT file format standards V1.1. The visualization component of the GMI-IPS uses Matplotlib extensively and has several functions ranging in complexity. First, it creates a model background curtain for the flight (time versus model eta levels) with the interpolated flight data superimposed on the curtain. Secondly, it creates a time-series plot of the interpolated flight data. Lastly, the visualization component creates averaged 2D model slices (longitude versus latitude) with overlaid flight track circles at key pressure levels. The GMI-IPS consists of a handful of classes and supporting functionality that have been generalized to be compatible with any ICARTT file that adheres to the base class definition. The base class represents a generic ICARTT entry, only defining a single time entry and 3D spatial positioning parameters. Other classes inherit from this base class; several classes for input ICARTT instrument files, which contain the necessary flight positioning information as a basis for data processing, as well as other classes for output ICARTT files, which contain the interpolated model data. Utility classes provide functionality for routine procedures such as: comparing field names among ICARTT files, reading ICARTT entries from a data file and storing them in data structures, and returning a reduced spatial grid based on a collection of ICARTT entries. Although the GMI-IPS is compatible with GMI model data, it can be adapted with reasonable effort for any simulation that creates Hierarchical Data Format (HDF) files. The same can be said of its adaptability to ICARTT files outside of the context of the ATom mission. The GMI-IPS contains just under 30,000 lines of code, eight classes, and a dozen drivers and utility programs. It is maintained with GIT source code management and has been used to deliver processed GMI model data for the ATom campaigns that have taken place to date.

Damon, M. R.↗

On the evolution of vortex rings with swirl

A laminar vortex ring with swirl, which has the meridional velocity component inside the vortex core, was experimentally generated by the brief fluid ejection from a rotating outlet. The evolution of the vortex ring was investigated with flow visualizations and particle image velocimetry measurements in order to find the influence of swirling flow in particular upon the transition to turbulence. Immediately after the formation of a vortex ring with swirl, a columnar strong vortex along the symmetric axis is observed in all cases of the present experiment. Then the characteristic fluid discharging from a vortex ring with swirl referred to as “peeling off” appears. The amount of discharging fluid due to the “peeling off” increases with the angular velocity of the rotating outlet. We conjectured that the mechanism generating the “peeling off” is related to the columnar strong vortex by close observations of the spatio-temporal development of the vorticity distribution and the cutting 3D images constructed from the successive cross sections of a vortex ring. While a laminar vortex ring without swirl may develop azimuthal waves around its circumference at some later time and the ring structure subsequently breaks, the swirling flow in a vortex ring core reduces the amplification rate of the azimuthal wavy deformation and preserved its ring structure. Then the traveling distance of a vortex ring can be extended using the swirl flow under certain conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Three-dimensional visualization and characterization of bone structure using reconstructed in-vitro μCT images: A pilot study for bone microarchitecture analysis

Micro Computed Tomography (μCT) has been largely used to perform micrometer scale imaging of specimens, bone biopsies and small animals for the study of porous or cavity-containing objects. One of its favored applications is for assessing structural properties of bone. In this research, we perform a pilot study to visualize and characterize bone structure of a chicken bone thigh, as well as to delineate its cortical and trabecular bone regions. We utilize an In-Vitro μCT scanner Skyscan 1173 to acquire a three dimensional image data of a chicken bone thigh. The thigh was scanned using X-ray voltage of 45 kV and current of 150 μA. The reconstructed images have spatial resolution of 142.50 μm/pixel. Using image processing and analysis e.i segmentation by thresholding the gray values (which represent the pseudo density) and binarizing the images, we were able to visualize each part of the bone, i.e., the cortical and trabecular regions. Total volume of the bone is 4663.63 mm{sup 3}, and the surface area of the bone is 7913.42 mm{sup 2}. The volume of the cortical is approximately 1988.62 mm{sup 3} which is nearly 42.64% of the total bone volume. This pilot study has confirmed that the μCT is capable of quantifying 3D bone structural properties and defining its regions separately. For further development, these results can be improved for understanding the pathophysiology of bone abnormality, testing the efficacy of pharmaceutical intervention, or estimating bone biomechanical properties.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Direct visualization of polaron formation in the thermoelectric SnSe

SnSe is a layered material that currently holds the record for bulk thermoelectric efficiency. The primary determinant of this high efficiency is thought to be the anomalously low thermal conductivity resulting from strong anharmonic coupling within the phonon system. Here we show that the nature of the carrier system in SnSe is also determined by strong coupling to phonons by directly visualizing polaron formation in the material. We employ ultrafast electron diffraction and diffuse scattering to track the response of phonons in both momentum and time to the photodoping of free carriers across the bandgap, observing the bimodal and anisotropic lattice distortions that drive carrier localization. Relatively large (18.7 Å), quasi-one-dimensional (1D) polarons are formed on the 300-fs timescale with smaller (4.2 Å) 3D polarons taking an order of magnitude longer (4 ps) to form. This difference appears to be a consequence of the profoundly anisotropic electron–phonon coupling in SnSe, with strong Fröhlich coupling only to zone-center polar optical phonons. These results demonstrate a high density of polarons in SnSe at optimal doping levels. Strong electron-phonon coupling is critical to the thermoelectric performance of this benchmark material and, potentially, high performance thermoelectrics more generally.

36 MATERIALS SCIENCE↗

WE-D-BRE-07: Variance-Based Sensitivity Analysis to Quantify the Impact of Biological Uncertainties in Particle Therapy

Purpose: In particle therapy, treatment planning and evaluation are frequently based on biological models to estimate the relative biological effectiveness (RBE) or the equivalent dose in 2 Gy fractions (EQD2). In the context of the linear-quadratic model, these quantities depend on biological parameters (α, β) for ions as well as for the reference radiation and on the dose per fraction. The needed biological parameters as well as their dependency on ion species and ion energy typically are subject to large (relative) uncertainties of up to 20–40% or even more. Therefore it is necessary to estimate the resulting uncertainties in e.g. RBE or EQD2 caused by the uncertainties of the relevant input parameters. Methods: We use a variance-based sensitivity analysis (SA) approach, in which uncertainties in input parameters are modeled by random number distributions. The evaluated function is executed 10{sup 4} to 10{sup 6} times, each run with a different set of input parameters, randomly varied according to their assigned distribution. The sensitivity S is a variance-based ranking (from S = 0, no impact, to S = 1, only influential part) of the impact of input uncertainties. The SA approach is implemented for carbon ion treatment plans on 3D patient data, providing information about variations (and their origin) in RBE and EQD2. Results: The quantification enables 3D sensitivity maps, showing dependencies of RBE and EQD2 on different input uncertainties. The high number of runs allows displaying the interplay between different input uncertainties. The SA identifies input parameter combinations which result in extreme deviations of the result and the input parameter for which an uncertainty reduction is the most rewarding. Conclusion: The presented variance-based SA provides advantageous properties in terms of visualization and quantification of (biological) uncertainties and their impact. The method is very flexible, model independent, and enables a broad assessment of uncertainties. Supported by DFG grant WI 3745/1-1 and DFG cluster of excellence: Munich-Centre for Advanced Photonics.

60 APPLIED LIFE SCIENCES↗

Visual feedback and guided balance training in an immersive virtual reality environment for lower extremity rehabilitation

Balance training is essential for physical rehabilitation procedures, as it can improve functional mobility and enhance cognitive coordination. However, conventional balance training methods may have limitations in terms of motivation, real-time objective feedback, and personalization, which a virtual reality (VR) setup may provide a better alternative. In this work, we present an immersive VR training environment for lower extremity balance rehabilitation with real-time guidance and feedback. The VR training environment immerses the user in a 3D ice rink model where a virtual coach (agent) leads them through a series of balance poses, and the user controls a trainee avatar with their own movements. Here we developed two coaching styles: positive-reinforcement and autonomous-supportive, and two viewpoints of the trainee avatar: first-person and third-person. The proposed environment was evaluated in a user study with healthy, non-clinical participants (n = 16, 24.4 ± 5.7 years old, 9 females). Our results show that participants showed stronger performance in the positive-reinforcement style compared to the autonomous-supportive style. Additionally, in the third-person viewpoint, the participants exhibited more stability in the positive-reinforcement style compared to the autonomous-supportive style. For viewpoint, participants exhibited stronger performance in the first-person viewpoint compared to third-person in the autonomous-supportive style, while they were comparable in the positive-reinforcement style. We observed no significant effects on the foot height and number of mistakes. Furthermore, we report the analysis of user performance with balance training poses and subjective measures based on questionnaires to assess the user experience, usability, and task load. The proposed VR balance training could offer an interactive, adaptive, and engaging environment and open new potential research directions for lower extremity rehabilitation.

97 MATHEMATICS AND COMPUTING↗

Machine-Learning-based Algorithms for Automated Image Segmentation Techniques of Transmission X-ray Microscopy (TXM)

Four state-of-the-art Deep Learning-based Convolutional Neural Networks (DCNN) were applied to automate the semantic segmentation of a 3D Transmission x-ray Microscopy (TXM) nanotomography image data. The standard U-Net architecture as baseline along with UNet++, PSPNet, and DeepLab v3+ networks were trained to segment the microstructural features of an AA7075 micropillar. A workflow was established to evaluate and compare the DCNN prediction dataset with the manually segmented features using the Intersection of Union (IoU) scores, time of training, confusion matrix, and visual assessment. Comparing all model segmentation accuracy metrics, it was found that using pre-trained models as a backbone along with appropriate training encoder-decoder architecture of the Unet++ can robustly handle large volumes of x-ray radiographic images in a reasonable amount of time. This opens a new window for handling accurate and efficient image segmentation of in situ time-dependent 4D x-ray microscopy experimental datasets.

36 MATERIALS SCIENCE↗

Measurements of the streamwise vortical structures in a plane mixing layer

The 3D structure of a plane two-stream mixing layer of velocity ratio 0.6 and originating from laminar initial boundary layers was investigated through direct measurements made in a specially constructed mixing-layer wind tunnel. The main objective of the study was to establish quantitatively the presence and the role of the secondary streamwise vortex structure (of the kind that has been shown in past flow visualization investigations to ride among the primary spanwise vortices) in the development of a plane turbulent mixing layer at relatively high Reynolds numbers. Results indicate that the instability leading to the formation of streamwise vortices is initially amplified just downstream of the first spanwise roll-up. The streamwise vortices, which first appear in clusters containing vorticity of both signs, realign further downstream to form counterrotating pairs. Due to the amalgamation of like-sign vortices, the streamwise vortex spacing increases in a stepwise fashion.

Bell, James H.↗

The exploration metaphor

NASA's experience in planetary exploration has demonstrated that the desktop workstation is inadequate for many visualization situations. The primary mission displays for the unmanned Surveyor missions to the moon during the mid-1960's, for example, were environmental images assembled on the inside surfaces of spherical shells. Future exploration missions will greatly benefit from advances in digital computer and display technology, but there remain unmet user interface needs. Alternative user interfaces and metaphors are needed for planetary exploration and other interactions with complex spatial environments. These interfaces and metaphors would enable the user to directly explore environments and naturally manipulate objects in those environments. Personal simulators, virtual workstations, and telepresence user interfaces are systems capable of providing this integration of user space and task space. The Exploration Metaphor is a useful concept for guiding the design of user interfaces for virtual environments and telepresence. To apply the Exploration Metaphor is to assert that computing is like exploration, and to support objects, operations, and contexts comparable to those encountered in the exploration of natural environments. The Exploration Metaphor, under development for user interfaces in support of NASA's planetary exploration missions and goals, will also benefit other applications where complex spatial information must be visualized. Visualization methods and systems for planetary exploration are becoming increasingly integrated and interactive as computing technology improves. These advances will benefit from virtual environment and telepresence interface technology. A key development has been the processing of multiple images and other sensor data to create detailed digital models of the planets and moons. Data from images of the Earth, Mars, and Miranda, for example, have been converted into 3D models, and dynamic virtual fly-overs have been computed as demonstrations. Similar processing of lower altitude photography and the use of computer aided design tools promise to produce very detailed models in the future.

Mcgreevy, Michael W.↗

Experimental study of crystallization fouling by calcium carbonate: Effects of surface structure and material

We report mitigation of crystallization fouling on heat transfer surface is crucial for developing highly efficient desalination techniques, because the tenacious and porous crystal layer can drastically deteriorate system efficiency. Surface modification is a potential approach to decrease scale formation and promote crystal removal, but some of the mechanisms between crystallization and substrate properties are still not fully understood. In this study, we use dual-angle microscopic visualization and image processing to quantify three-dimensional crystallization behavior and analyze the impact of surface structure and material. The transition of crystal size and distance with roughness and wettability are identified from statistical analysis. The determination of crystal porosity and effective thermal conductivity is enabled by visualized geometry and measured deposition mass. This study shed lights on potential strategies of surface modification to mitigate scaling in desalination by enhancing scale removal and localizing crystal growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cardiac motion correction based on partial angle reconstructed images in x-ray CT

Purpose: Cardiac x-ray CT imaging is still challenging due to heart motion, which cannot be ignored even with the current rotation speed of the equipment. In response, many algorithms have been developed to compensate remaining motion artifacts by estimating the motion using projection data or reconstructed images. In these algorithms, accurate motion estimation is critical to the compensated image quality. In addition, since the scan range is directly related to the radiation dose, it is preferable to minimize the scan range in motion estimation. In this paper, the authors propose a novel motion estimation and compensation algorithm using a sinogram with a rotation angle of less than 360°. The algorithm estimates the motion of the whole heart area using two opposite 3D partial angle reconstructed (PAR) images and compensates the motion in the reconstruction process. Methods: A CT system scans the thoracic area including the heart over an angular range of 180° + α + β, where α and β denote the detector fan angle and an additional partial angle, respectively. The obtained cone-beam projection data are converted into cone-parallel geometry via row-wise fan-to-parallel rebinning. Two conjugate 3D PAR images, whose center projection angles are separated by 180°, are then reconstructed with an angular range of β, which is considerably smaller than a short scan range of 180° + α. Although these images include limited view angle artifacts that disturb accurate motion estimation, they have considerably better temporal resolution than a short scan image. Hence, after preprocessing these artifacts, the authors estimate a motion model during a half rotation for a whole field of view via nonrigid registration between the images. Finally, motion-compensated image reconstruction is performed at a target phase by incorporating the estimated motion model. The target phase is selected as that corresponding to a view angle that is orthogonal to the center view angles of two conjugate PAR images. To evaluate the proposed algorithm, digital XCAT and physical dynamic cardiac phantom datasets are used. The XCAT phantom datasets were generated with heart rates of 70 and 100 bpm, respectively, by assuming a system rotation time of 300 ms. A physical dynamic cardiac phantom was scanned using a slowly rotating XCT system so that the effective heart rate will be 70 bpm for a system rotation speed of 300 ms. Results: In the XCAT phantom experiment, motion-compensated 3D images obtained from the proposed algorithm show coronary arteries with fewer motion artifacts for all phases. Moreover, object boundaries contaminated by motion are well restored. Even though object positions and boundary shapes are still somewhat different from the ground truth in some cases, the authors see that visibilities of coronary arteries are improved noticeably and motion artifacts are reduced considerably. The physical phantom study also shows that the visual quality of motion-compensated images is greatly improved. Conclusions: The authors propose a novel PAR image-based cardiac motion estimation and compensation algorithm. The algorithm requires an angular scan range of less than 360°. The excellent performance of the proposed algorithm is illustrated by using digital XCAT and physical dynamic cardiac phantom datasets.

60 APPLIED LIFE SCIENCES↗