Search NASA⌕ Search

SEARCH · Search NASA

Results for “imaging methods”

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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-D-206-06: Task-Specific Optimization of Scintillator Thickness for CMOS-Detector Based Cone-Beam Breast CT

Purpose: To optimize the cesium iodide (CsI:Tl) scintillator thickness in a complimentary metal-oxide semiconductor (CMOS)-based detector for use in dedicated cone-beam breast CT. Methods: The imaging task considered was the detection of a microcalcification cluster comprising six 220µm diameter calcium carbonate spheres, arranged in the form of a regular pentagon with 2 mm spacing on its sides and a central calcification, similar to that in ACR-recommended mammography accreditation phantom, at a mean glandular dose of 4.5 mGy. Generalized parallel-cascades based linear systems analysis was used to determine Fourier-domain image quality metrics in reconstructed object space, from which the detectability index inclusive of anatomical noise was determined for a non-prewhitening numerical observer. For 300 projections over 2π, magnification-associated focal-spot blur, Monte Carlo derived x-ray scatter, K-fluorescent emission and reabsorption within CsI:Tl, CsI:Tl quantum efficiency and optical blur, fiberoptic plate transmission efficiency and blur, CMOS quantum efficiency, pixel aperture function and additive noise, and filtered back-projection to isotropic 105µm voxel pitch with bilinear interpolation were modeled. Imaging geometry of a clinical prototype breast CT system, a 60 kV Cu/Al filtered x-ray spectrum from 0.3 mm focal spot incident on a 14 cm diameter semi-ellipsoidal breast were used to determine the detectability index for 300–600 µm thick (75µm increments) CsI:Tl. The CsI:Tl thickness that maximized the detectability index was considered optimal. Results: The limiting resolution (10% modulation transfer function, MTF) progressively decreased with increasing CsI:Tl thickness. The zero-frequency detective quantum efficiency, DQE(0), in projection space increased with increasing CsI:Tl thickness. The maximum detectability index was achieved with 525µm thick CsI:Tl scintillator. Reduced MTF at mid-to-high frequencies for 600µm thick CsI:Tl lowered the detectability index than 525µm CsI:Tl. Conclusion: For the x-ray spectrum and imaging conditions considered, a 525µm thick CsI:Tl scintillator integrated with the CMOS detector is optimal for detecting microcalcification cluster. Funding support: Supported in part by NIH R01 CA195512. The contents are solely the responsibility of the authors and do not reflect the official views of the NIH or the NCI. Disclosures: SV, GV and AK - Research collaboration, Koning Corp., West Henrietta, NY.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Facile approach to prepare hollow core–shell NiO microspherers for supercapacitor electrodes

A facile lamellar template method (see image) has been developed for the preparation of uniform hollow core–shell structure NiO (HCS–NiO) with a nanoarchitectured wall structure. The prepared NiO was found to be highly crystalline in uniform microstructures with high specific surface area and pore volume. The results indicated that ethanol interacted with trisodium citrate played an important role for the formation of hollow core–shell spheres. On the basis of the analysis of the composition and the morphology, a possible formation mechanism was investigated. NiO microspheres with hollow core–shell showed excellent capacitive properties. The exceptional cyclic, structural and electrochemical stability with ∼95% coulombic efficiency, and very low ESR value from impedance measurements promised good utility value of hollow core–shell NiO material in fabricating a wide range of high-performance electrochemical supercapacitors. - The hollow core–shell NiO was prepared with a facile lamellar template method. The prepared NiO show higher capacitance, lower ion diffusion resistance and better electroactive surface utilization for Faradaic reactions. - Highlights: • Formation of hollow core–shell NiO via a novel and facile precipitation route. • Exhibited uniform feature sizes and high surface area of hollow core–shell NiO. • Synthesized NiO has high specific capacitance ( 448 F g{sup 1}) and very low ESR value. • Increased 20% of long life cycles capability after 500 charge–discharge cycles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Characterization and Modeling of Cell Wall Imperfections in Aluminum Honeycomb Cores using X-ray CT Imaging

Honeycomb cell imperfections may significantly affect the failure of sandwich composites in core crushing and shear modes. Identification and quantification of the different types of imperfections and evaluation of their effects on the mechanical response of honeycomb cores in flatwise compression is presented. X-Ray computed tomography (CT) scans of honeycomb cores in co-cured composite sandwich panels were performed. Image analysis methods are applied to characterize and quantify core cell geometry and cell wall waviness. An approach is developed for extracting cell wall mid surface planes and waviness for constructing finite element models using shell elements. Analysis results from models with the as-manufactured imperfections are compared to results from models with idealized geometries to evaluate the effect of the as-manufactured imperfections on the initiation of cell wall buckling and subsequent core crushing.

aluminum honeycomb↗
DOE OSTI2016

Coupling sky images with radiative transfer models: a new method to estimate cloud optical depth

A method for retrieving cloud optical depth ( τ c ) using a UCSD developed ground-based sky imager (USI) is presented. The radiance red–blue ratio (RRBR) method is motivated from the analysis of simulated images of various τ c produced by a radiative transfer model (RTM). From these images the basic parameters affecting the radiance and red–blue ratio (RBR) of a pixel are identified as the solar zenith angle ( θ 0 ), τ c , solar pixel angle/scattering angle ( θ s ), and pixel zenith angle/view angle ( θ z ). The effects of these parameters are described and the functions for radiance, I λ τ c , θ 0 , θ s , θ z , and RBR τ c , θ 0 , θ s , θ z are retrieved from the RTM results. RBR, which is commonly used for cloud detection in sky images, provides non-unique solutions for τ c , where RBR increases with τ c up to about τ c = 1 (depending on other parameters) and then decreases. Therefore, the RRBR algorithm uses the measured I λ meas θ s , θ z , in addition to RBR meas θ s , θ z , to obtain a unique solution for τ c . The RRBR method is applied to images of liquid water clouds taken by a USI at the Oklahoma Atmospheric Radiation Measurement (ARM) program site over the course of 220 days and compared against measurements from a microwave radiometer (MWR) and output from the Min et al. (2003) method for overcast skies. τ c values ranged from 0 to 80 with values over 80, being capped and registered as 80. A τ c RMSE of 2.5 between the Min et al. (2003) method and the USI are observed. The MWR and USI have an RMSE of 2.2, which is well within the uncertainty of the MWR. In conclusion, the procedure developed here provides a foundation to test and develop other cloud detection algorithms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

In situ photoacoustic characterization for porous silicon growing: Detection principles

There are a few methodologies for monitoring the in-situ formation of Porous Silicon (PS). One of the methodologies is photoacoustic. Previous works that reported the use of photoacoustic to study the PS formation do not provide the physical explanation of the origin of the signal. In this paper, a physical explanation of the origin of the photoacoustic signal during the PS etching is provided. The incident modulated radiation and changes in the reflectance are taken as thermal sources. In this paper, a useful methodology is proposed to determine the etching rate, porosity, and refractive index of a PS film by the determination of the sample thickness, using scanning electron microscopy images. This method was developed by carrying out two different experiments using the same anodization conditions. The first experiment consisted of growth of the samples with different etching times to prove the periodicity of the photoacoustic signal, while the second one considered the growth samples using three different wavelengths that are correlated with the period of the photoacoustic signal. The last experiment showed that the period of the photoacoustic signal is proportional to the laser wavelength.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

SU-G-206-05: A Comparison of Head Phantoms Used for Dose Determination in Imaging Procedures

Purpose: To determine similarities and differences between various head phantoms that might be used for dose measurements in diagnostic imaging procedures. Methods: We chose four frequently used anthropomorphic head phantoms (SK-150, PBU-50, RS-240T and Alderson Rando), a computational patient phantom (Zubal) and the CTDI head phantom for comparison in our study. We did a CT scan of the head phantoms using the same protocol and compared their dimensions and CT numbers. The scan data was used to calculate dose values for each of the phantoms using EGSnrc Monte Carlo software. An .egsphant file was constructed to describe these phantoms using a Visual C++ program for DOSXYZnrc/EGSnrc simulation. The lens dose was calculated for a simulated CBCT scan using DOSXYZnrc/EGSnrc and the calculated doses were validated with measurements using Gafchromic film and an ionization chamber. Similar calculations and measurements were made for PA radiography to investigate the attenuation and backscatter differences between these phantoms. We used the Zubal phantom as the standard for comparison since it was developed based on a CT scan of a patient. Results: The lens dose for the Alderson Rando phantom is around 9% different than the Zubal phantom, while the lens dose for the PBU-50 phantom was about 50% higher, possibly because its skull thickness and the density of bone and soft tissue are lower than anthropometric values. The lens dose for the CTDI phantom is about 500% higher because of its totally different structure. The entrance dose profiles are similar for the five anthropomorphic phantoms, while that for the CTDI phantom was distinctly different. Conclusion: The CTDI and PBU-50 head phantoms have substantially larger lens dose estimates in CBCT. The other four head phantoms have similar entrance dose with backscatter hence should be preferred for dose measurement in imaging procedures of the head. Partial support from NIH Grant R01-EB002873 and Toshiba Medical Systems Corp.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Core Characterization of Bakken Shale from the Bedwell 33-52-1-1H Well, Sheridan County, Montana

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia were used to characterize core of the Bedwell 33-52-1-1H well (API 25-091-21920), a wildcat field in Sheridan County, Montana. The primary impetus of this work is to provide a unique dataset to researchers at universities, national laboratories, geological surveys, and other locations for additional analyses. The resultant datasets are presented in this report and can be accessed from NETL's Energy Data eXchange (EDX) online system using the following link: https://edx.netl.doe.gov/dataset/bedwell-33-52-1-1h-well. The equipment and techniques used to characterize the full core were non-destructive, enabling future examinations and analyses to be performed on these cores. However, none of the equipment used was suitable for direct visualization of the pore space in fine-grained structures such as those found in shales; fractures, discontinuities, and millimeter scale features were readily detectable with the methods tested. Imaging with the NETL medical CT scanner was performed on the entire core. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), P-wave, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for more detailed analysis. A selection of samples (Table 1) were cored, sliced, and powdered for additional characterization analyses. These analyses did involve destructive subsampling of the full core at discrete sections, and the results are presented here in the context of the larger core description. Selected sections of the core were removed from the full core for detailed analyses using the following destructive techniques: bulk inductively coupled plasma-optical emission spectroscopy (ICP-OES), isolation of carbon and sulfur, Fourier transform infrared (FTIR) spectroscopy, scanning electron microscopy (SEM), and higher resolution CT scan with NETL’s micro-CT scanning systems. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core; the resulting macro and micro descriptions are relevant to many subsurface energy related examinations traditionally performed at NETL.

04 OIL SHALES AND TAR SANDS↗
NASA NTRS1999

Global Geometric Properties of Martian Impact Craters: A Preliminary Assessment Using Mars Orbiter Laser Altimeter (MOLA)

Impact craters on Mars have been used to provide fundamental insights into the properties of the martian crust, the role of volatiles, the relative age of the surface, and on the physics of impact cratering in the Solar System. Before the three-dimensional information provided by the Mars Orbiter Laser Altimeter (MOLA) instrument which is currently operating in Mars orbit aboard the Mars Global Surveyor (MGS), impact features were characterized morphologically using orbital images from Mariner 9 and Viking. Fresh-appearing craters were identified and measurements of their geometric properties were derived from various image-based methods. MOLA measurements can now provide a global sample of topographic cross-sections of martian impact features as small as approx. 2 km in diameter, to basin-scale features. We have previously examined MOLA cross-sections of Northern Hemisphere and North Polar Region impact features, but were unable to consider the global characteristics of these ubiquitous landforms. Here we present our preliminary assessment of the geometric properties of a globally-distributed sample of martian impact craters, most of which were sampled during the initial stages of the MGS mapping mission (i.e., the first 600 orbits). Our aim is to develop a framework for reconsidering theories concerning impact cratering in the martian environment. This first global analysis is focused upon topographically-fresh impact craters, defined here on the basis of MOLA topographic profiles that cross the central cavities of craters that can be observed in Viking-based MDIM global image mosaics. We have considered crater depths, rim heights, ejecta topologies, cross-sectional "shapes", and simple physical models for ejecta emplacement. To date (May, 1999), we have measured the geometric properties of over 1300 impact craters in the 2 to 350 km diameter size interval. A large fraction of these measured craters were sampled with cavity-center cross-sections during the first two months of MGS mapping. Many of these craters are included in Nadine Barlow's Catalogue of Martian Impact Craters, although we have treated simple craters smaller than about 7 km in greater detail than all previous investigations. Additional information is contained in the original extended abstract.

Garvin, J. B.↗
DOE OSTI2017

Image binarization to calculate porosity of porous anodic oxides and derivation of porosity vs current

Highlights: • A facile method based on image binarization is introduced to calculate porosity. • The method is compared with two popular methods, to show higher accuracy. • The method is tested to show robustness under different brightness and contrast. • Mathematical derivation of porosity vs current is provided. - Abstract: Current methods to calculate porosity of tubular or porous structures are either inaccurate or complicated for operation. We introduce a facile method (M7) based on image binarization of top-view specimen images from scanning electron microscopy. M7 is both accurate and convenient, and shows high robustness regardless of external parameters like brightness or contrast. Relation between porosity and applied current during galvanostatic anodization is derived. The porosity-current relation differs in aqueous and organic solution, which can be explained as different conductivity of two electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2015

SU-E-J-27: Shifting Multiple EPID Imager Layers to Improve Image Quality and Resolution in MV CBCT

Purpose: Vertical stacking of four conventional EPID layers can improve DQE for MV-CBCT applications. We hypothesize that shifting each layer laterally by half a pixel relative to the layer above, will improve the contrast-to-noise ratio (CNR) and image resolution. Methods: For CNR assessment, a 20 cm diameter digital phantom with 8 inserts is created. The attenuation coefficient of the phantom is similar to lung at the average energy of a 6 MV photon beam. The inserts have attenuations 1, 2…8 times of lung. One of the inserts is close to soft tissue, resembling the case of a tumor in lung. For resolution assessment, a digital phantom featuring a bar pattern is created. The phantom has an attenuation coefficient similar to soft tissue and the bars have an attenuation coefficient of calcium sulfate. A 2 MeV photon beam is attenuated through these phantoms and hits each of the four stacked detector layers. Each successive layer is shifted by half a pixel in the x only, y only, and x and y (combined) directions, respectively. Blurring and statistical noise are added to the projections. Projections from one, two, three and four layers are used for reconstruction. CNR and image resolution are evaluated and compared. Results: When projections from multiple layers are combined for reconstruction, CNR increases with the number of layers involved. CNR in reconstructions from two, three and four layers are 1.4, 1.7 and 1.99 times that from one layer. The resolution from the shifted four layer detector is also improved from a single layer. In a comparison between one layer versus four layers in this preliminary study, the resolution from four shifted layers is at least 20% better. Conclusion: Layer-shifting in a stacked EPID imager design enhances resolution as well as CNR for half scan MV-CBCT. The project described was supported, in part, by a grant from Varian Medical Systems, Inc., and Award No. R01CA188446-01 from the National Cancer Institute. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Cancer Institute or the National Institutes of Health.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Frequency-modulated atomic force microscopy operation by imaging at the frequency shift minimum: The dip-df mode

In frequency modulated non-contact atomic force microscopy, the change of the cantilever frequency (Δf) is used as the input signal for the topography feedback loop. Around the Δf(z) minimum, however, stable feedback operation is challenging using a standard proportional-integral-derivative (PID) feedback design due to the change of sign in the slope. When operated under liquid conditions, it is furthermore difficult to address the attractive interaction regime due to its often moderate peakedness. Additionally, the Δf signal level changes severely with time in this environment due to drift of the cantilever frequency f{sub 0} and, thus, requires constant adjustment. Here, we present an approach overcoming these obstacles by using the derivative of Δf with respect to z as the input signal for the topography feedback loop. Rather than regulating the absolute value to a preset setpoint, the slope of the Δf with respect to z is regulated to zero. This new measurement mode not only makes the minimum of the Δf(z) curve directly accessible, but it also benefits from greatly increased operation stability due to its immunity against f{sub 0} drift. We present isosurfaces of the Δf minimum acquired on the calcite CaCO{sub 3}(101{sup ¯}4) surface in liquid environment, demonstrating the capability of our method to image in the attractive tip-sample interaction regime.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Large-Scale Image Analysis for Investigating Spatio-Temporal Changes in Nuclear DNA Damage Caused by Nitrogen Atmospheric Pressure Plasma Jets

The effective clinical application of atmospheric pressure plasma jet (APPJ) treatments requires a well-founded methodology that can describe the interactions between the plasma jet and a treated sample and the temporal and spatial changes that result from the treatment. In this study, we developed a large-scale image analysis method to identify the cell-cycle stage and quantify damage to nuclear DNA in single cells. The method was then tested and used to examine spatio-temporal distributions of nuclear DNA damage in two cell lines from the same anatomic location, namely the oral cavity, after treatment with a nitrogen APPJ. One cell line was malignant, and the other, nonmalignant. The results showed that DNA damage in cancer cells was maximized at the plasma jet treatment region, where the APPJ directly contacted the sample, and declined radially outward. As incubation continued, DNA damage in cancer cells decreased slightly over the first 4 h before rapidly decreasing by approximately 60% at 8 h post-treatment. In nonmalignant cells, no damage was observed within 1 h after treatment, but damage was detected 2 h after treatment. Notably, the damage was 5-fold less than that detected in irradiated cancer cells. Moreover, examining damage with respect to the cell cycle showed that S phase cells were more susceptible to DNA damage than either G1 or G2 phase cells. The proposed methodology for large-scale image analysis is not limited to APPJ post-treatment applications and can be utilized to evaluate biological samples affected by any type of radiation, and, more so, the cell-cycle classification can be used on any cell type with any nuclear DNA staining.

cellular imaging↗
DOE OSTI2025

Exploiting correlations in multi-coincidence Coulomb explosion patterns for differentiating molecular structures using machine learning

Coulomb explosion imaging (CEI) is a powerful technique for capturing the real-time motion of individual atoms during ultrafast photochemical reactions. CEI generates high-dimensional data with naturally embedded correlations that allow mapping the coordinated motion of nuclei in molecules. This enables reliable separation of competing reaction pathways and makes this approach uniquely suited for characterizing weak reaction channels. However, rich information contained in experimental CEI patterns remains largely underexploited due to challenges in visualizing correlations between multiple observables in multi-dimensional parameter space. Here we present a new approach to CEI of intermediate-sized polyatomic molecules, detecting up to eight ionic fragments in coincidence and leveraging machine-learning-based analysis to identify patterns and correlations in the resulting high-dimensional momentum-space data, enabling robust molecular structure identification and differentiation. Our approach provides high-dimensional background-free data encoding exceptionally rich structural information and establishes an automated, scalable framework for extracting insightful information from the data. As a demonstration, we apply this method to image and distinguish dichloroethylene isomers, showcasing its potential for broader applications in molecular imaging. Our results pave the way for channel-specific analysis of ultrafast structural dynamics in chemically relevant systems, particularly for disentangling mixed reaction pathways and detecting contributions from weak channels and minority species.

Chemical Physics (physics.chem-ph)↗
DOE OSTI2015

Robust PET-guided intensity-modulated radiation therapy

Purpose: Functional image guided intensity-modulated radiation therapy has the potential to improve cancer treatment quality by basing treatment parameters such as heterogeneous dose distributions information derived from imaging. However, such heterogeneous dose distributions are subject to imaging uncertainty. In this paper, the authors develop a robust optimization model to design plans that are desensitized to imaging uncertainty. Methods: Starting from the pretreatment fluorodeoxyglucose-positron emission tomography scans, the authors use the raw voxel standard uptake values (SUVs) as input into a series of intermediate functions to transform the SUV into a desired dose. The calculated desired doses were used as an input into a robust optimization model to generate beamlet intensities. For each voxel, the authors assume that the true SUV cannot be observed but instead resides in an interval centered on the nominal (i.e., observed) SUV. Then the authors evaluated the nominal and robust solutions through a simulation study. The simulation considered the effect of the true SUV being different from the nominal SUV on the quality of the treatment plan. Treatment plans were compared on the metrics of objective function value and tumor control probability (TCP). Results: Computational results demonstrate the potential for improvements in tumor control probability and deviation from the desired dose distribution compared to a nonrobust model while maintaining acceptable tissue dose. Conclusions: Robust optimization can help design treatment plans that are more stable in the presence of image value uncertainties.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-J-101: Retroactive Calculation of TLD and Film Dose in Anthropomorphic Phantom as Assessment of Updated TPS Performance

Purpose: To demonstrate a quick and comprehensive method verifying the accuracy of the updated dose model by recalculating dose distribution in an anthropomorphic phantom with a new version of the TPS and comparing the results to measured values. Methods: CT images and IMRT plan of an RPC anthropomorphic head phantom, previously calculated by Pinnacle 9.0, was re-computed using Pinnacle 9.2 and 9.6. The dosimeters within the phantom include four TLD capsules representing a primary PTV, two TLD capsules representing a secondary PTV, and two TLD capsules representing an organ at risk. Also included were three sheets of Gafchromic film. Performance of the updated TPS version was assessed by recalculating point doses and dose profiles corresponding to TLD and film position respectively and then comparing the results to reported values by the RPC. Results: Comparing calculated doses to reported measured doses from the RPC yielded an average disagreement of 1.48%, 2.04% and 2.10% for versions 9.0, 9.2, 9.6 respectively. Computed doses points all meet the RPC's passing criteria with the exception of the point representing the superior organ at risk in version 9.6. However, qualitative analysis of the recalculated dose profiles showed improved agreement with those of the RPC, especially in the penumbra region. Conclusion: This work has demonstrated the calculation results of Pinnacle 9.2 and 9.6 vs 9.0 version. Additionally, this study illustrates a method for the user to gain confidence upgrade to a newer version of the treatment planning system.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2003

Numerical Models of Human Circulatory System under Altered Gravity: Brain Circulation

A computational fluid dynamics (CFD) approach is presented to model the blood flow through the human circulatory system under altered gravity conditions. Models required for CFD simulation relevant to major hemodynamic issues are introduced such as non-Newtonian flow models governed by red blood cells, a model for arterial wall motion due to fluid-wall interactions, a vascular bed model for outflow boundary conditions, and a model for auto-regulation mechanism. The three-dimensional unsteady incompressible Navier-Stokes equations coupled with these models are solved iteratively using the pseudocompressibility method and dual time stepping. Moving wall boundary conditions from the first-order fluid-wall interaction model are used to study the influence of arterial wall distensibility on flow patterns and wall shear stresses during the heart pulse. A vascular bed modeling utilizing the analogy with electric circuits is coupled with an auto-regulation algorithm for multiple outflow boundaries. For the treatment of complex geometry, a chimera overset grid technique is adopted to obtain connectivity between arterial branches. For code validation, computed results are compared with experimental data for steady and unsteady non-Newtonian flows. Good agreement is obtained for both cases. In sin-type Gravity Benchmark Problems, gravity source terms are added to the Navier-Stokes equations to study the effect of gravitational variation on the human circulatory system. This computational approach is then applied to localized blood flows through a realistic carotid bifurcation and two Circle of Willis models, one using an idealized geometry and the other model using an anatomical data set. A three- dimensional anatomical Circle of Willis configuration is reconstructed from human-specific magnetic resonance images using an image segmentation method. The blood flow through these Circle of Willis models is simulated to provide means for studying gravitational effects on the brain circulation under auto-regulation.

Kim, Chang Sung↗
DOE OSTI2013

Transconvolution and the virtual positron emission tomograph-A new method for cross calibration in quantitative PET/CT imaging

Purpose: Positron emission tomography (PET)/computed tomography (CT) measurements on small lesions are impaired by the partial volume effect, which is intrinsically tied to the point spread function of the actual imaging system, including the reconstruction algorithms. The variability resulting from different point spread functions hinders the assessment of quantitative measurements in clinical routine and especially degrades comparability within multicenter trials. To improve quantitative comparability there is a need for methods to match different PET/CT systems through elimination of this systemic variability. Consequently, a new method was developed and tested that transforms the image of an object as produced by one tomograph to another image of the same object as it would have been seen by a different tomograph. The proposed new method, termed Transconvolution, compensates for differing imaging properties of different tomographs and particularly aims at quantitative comparability of PET/CT in the context of multicenter trials. Methods: To solve the problem of image normalization, the theory of Transconvolution was mathematically established together with new methods to handle point spread functions of different PET/CT systems. Knowing the point spread functions of two different imaging systems allows determining a Transconvolution function to convert one image into the other. This function is calculated by convolving one point spread function with the inverse of the other point spread function which, when adhering to certain boundary conditions such as the use of linear acquisition and image reconstruction methods, is a numerically accessible operation. For reliable measurement of such point spread functions characterizing different PET/CT systems, a dedicated solid-state phantom incorporating {sup 68}Ge/{sup 68}Ga filled spheres was developed. To iteratively determine and represent such point spread functions, exponential density functions in combination with a Gaussian distribution were introduced. Furthermore, simulation of a virtual PET system provided a standard imaging system with clearly defined properties to which the real PET systems were to be matched. A Hann window served as the modulation transfer function for the virtual PET. The Hann's apodization properties suppressed high spatial frequencies above a certain critical frequency, thereby fulfilling the above-mentioned boundary conditions. The determined point spread functions were subsequently used by the novel Transconvolution algorithm to match different PET/CT systems onto the virtual PET system. Finally, the theoretically elaborated Transconvolution method was validated transforming phantom images acquired on two different PET systems to nearly identical data sets, as they would be imaged by the virtual PET system. Results: The proposed Transconvolution method matched different PET/CT-systems for an improved and reproducible determination of a normalized activity concentration. The highest difference in measured activity concentration between the two different PET systems of 18.2% was found in spheres of 2 ml volume. Transconvolution reduced this difference down to 1.6%. In addition to reestablishing comparability the new method with its parameterization of point spread functions allowed a full characterization of imaging properties of the examined tomographs. Conclusions: By matching different tomographs to a virtual standardized imaging system, Transconvolution opens a new comprehensive method for cross calibration in quantitative PET imaging. The use of a virtual PET system restores comparability between data sets from different PET systems by exerting a common, reproducible, and defined partial volume effect.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Cavitation inception by the backscattering of pressure waves from a bubble interface

The secondary cavitation that occurs by the backscattering of focused ultrasound from a primary cavitation bubble caused by the negative pressure part of the ultrasound (Maxwell, et al., 2011) might be useful for the energy exchange due to bubble oscillations in High Intensity Focused Ultrasound (HIFU). The present study is concerned with the cavitation inception by the backscattering of ultrasound from a bubble. In the present experiment, a laser-induced bubble which is generated by a pulsed focused laser beam with high intensity is utilized as a primary cavitation bubble. After generating the bubble, focused ultrasound is emitted to the bubble. The acoustic field and the bubble motion are observed with a high-speed video camera. It is confirmed that the secondary cavitation bubble clouds are generated by the backscattering from the laser-induced bubble. The growth of cavitation bubble clouds is analyzed with the image processing method. The experimental results show that the height and width of the bubble clouds grow in stepwise during their evolution. The direct numerical simulations are also conducted for the backscattering of incident pressure waves from a bubble in order to evaluate a pressure field near the bubble. It is shown that the ratio of a bubble collapse time t{sub 0} to a characteristic time of wave propagation t{sub S}, η = t{sub 0}/t{sub s}, is an important determinant for generating negative pressure region by backscattering. The minimum pressure location by the backscattering in simulations is in good agreement with the experiment.

42 ENGINEERING↗