SEARCH · Search NASA
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.
Earth's atmosphere reduces a sensors ability in currently discriminating targets. Using radiometric correction to reduce the atmospheric effects may improve considerably the performance of an automatic image interpreter. Several methods for radiometric correction from the open literature are compared leading to the development of an atmospheric correction system.
Purpose: To estimate and remove the scatter contamination in the acquired projection of cone-beam CT (CBCT), to suppress the shading artifacts and improve the image quality without prior information. Methods: The uncorrected CBCT images containing shading artifacts are reconstructed by applying the standard FDK algorithm on CBCT raw projections. The uncorrected image is then segmented to generate an initial template image. To estimate scatter signal, the differences are calculated by subtracting the simulated projections of the template image from the raw projections. Since scatter signals are dominantly continuous and low-frequency in the projection domain, they are estimated by low-pass filtering the difference signals and subtracted from the raw CBCT projections to achieve the scatter correction. Finally, the corrected CBCT image is reconstructed from the corrected projection data. Since an accurate template image is not readily segmented from the uncorrected CBCT image, the proposed scheme is iterated until the produced template is not altered. Results: The proposed scheme is evaluated on the Catphan©600 phantom data and CBCT images acquired from a pelvis patient. The result shows that shading artifacts have been effectively suppressed by the proposed method. Using multi-detector CT (MDCT) images as reference, quantitative analysis is operated to measure the quality of corrected images. Compared to images without correction, the method proposed reduces the overall CT number error from over 200 HU to be less than 50 HU and can increase the spatial uniformity. Conclusion: An iterative strategy without relying on the prior information is proposed in this work to remove the shading artifacts due to scatter contamination in the projection domain. The method is evaluated in phantom and patient studies and the result shows that the image quality is remarkably improved. The proposed method is efficient and practical to address the poor image quality issue of CBCT images. This work is supported by the Zhejiang Provincial Natural Science Foundation of China (Grant No. LR16F010001), National High-tech R&D Program for Young Scientists by the Ministry of Science and Technology of China (Grant No. 2015AA020917).
A new image processing method for automated horizon detection is presented in this paper. This method was developed for the Mars Exploration Rovers mission’s surface operations to aid in determining orbiter visibility periods for downlink data, and is now being considered for all Mars surface missions and for a variety of applications. Starting with images received from Mars rovers or landers, the horizon pixel location is found by (1) obtaining the gradient images using the Sobel operator, (2) calculating the preliminary horizon locations by optimizing the maximum difference between sky and ground regions in the gradient domain, and (3) applying a multi-variable thresholding method. The results are analyzed using a variety of images from the Mars rovers, as an example one dataset’s results were categorized as 94.2% Good, 4.11% Okay, and 1.65% Poor. The pixel location corresponding to the found horizon in an image is converted to azimuth and elevation values using the camera model. The elevation of the terrain surrounding the rover or lander directly affects the duration orbiters are visible for downlink data volume, and when the sun is visible for image quality or solar energy consumption. Having an automated process to detect the elevation of the terrain allows the operation team to better predict the outcome of the plan, reducing heath and safety risks and allowing new operational limits to be defined.
Purpose: Accounting for tumor motion during radiation therapy is important to ensure that the tumor receives the prescribed dose. Increasing the field size to account for this motion exposes the surrounding healthy tissues to unnecessary radiation. In contrast to using motion-encompassing techniques to treat moving tumors, conformal radiation therapy (RT) uses a smaller field to track the tumor and adapts the beam aperture according to the motion detected. This work investigates and compares the performance of three markerless, EPID based, optical flow methods to track tumor motion with conformal RT. Methods: Three techniques were used to track the motions of a 3D printed lung tumor programmed to move according to the tumor of seven lung cancer patients. These techniques utilized a multi-resolution optical flow algorithm as the core computation for image registration. The first method (DIR) registers the incoming images with an initial reference frame, while the second method (RFSF) uses an adaptive reference frame and the third method (CU) uses preceding image frames for registration. The patient traces and errors were evaluated for the seven patients. Results: The average position errors for all patient traces were 0.12 ± 0.33 mm, −0.05 ± 0.04 mm and −0.28 ± 0.44 mm for CU, DIR and RFSF method respectively. The position errors distributed within 1 standard deviation are 0.74 mm, 0.37 mm and 0.96 mm respectively. The CU and RFSF algorithms are sensitive to the characteristics of the patient trace and produce a wider distribution of errors amongst patients. Although the mean error for the DIR method is negatively biased (−0.05 mm) for all patients, it has the narrowest distribution of position error, which can be corrected using an offset calibration. Conclusion: Three techniques of image registration and position update were studied. Using direct comparison with an initial frame yields the best performance. The authors would like to thank Dr.YeLin Suh for making the Cyberknife dataset available to us. Scholarship funding from the Natural Sciences and Engineering Research Council of Canada (NSERC) and CancerCare Manitoba Foundation is acknowledged.
Purpose: To assess regional changes in human lung ventilation and mechanics using four-dimensional computed tomography (4DCT) and deformable image registration. This work extends our prior analysis of the entire lung to a lobe-based analysis. Methods: 4DCT images acquired from 20 patients prior to radiation therapy (RT) were used for this analysis. Jacobian ventilation and motion maps were computed from the displacement field after deformable image registration between the end of expiration breathing phase and the end of inspiration breathing phase. The lobes were manually segmented on the reference phase by a medical physicist expert. The voxel-by-voxel ventilation and motion magnitude for all subjects were grouped by lobes and plotted into cumulative voxel frequency curves respectively. In addition, to eliminate the effect of different breathing efforts across subjects, we applied the inter-subject equivalent lung volume (ELV) method on a subset of the cohort and reevaluated the lobar ventilation. Results: 95% of voxels in the lung are expanding during inspiration. However, some local regions of lung tissue show far more expansion than others. The greatest expansion with respiration occurs within the lower lobes; between exhale and inhale the median expansion in lower lobes is approximately 15%, while the median expansion in upper lobes is 10%. This appears to be driven by a subset of lung tissues within the lobe that have greater expansion; twice the number of voxels in the lower lobes (20%) expand by > 30% when compared to the upper lobes (10%). Conclusion: Lung ventilation and motion show significant difference on the lobar level. There are different lobar fractions of driving voxels that contribute to the major expansion of the lung. This work was supported by NIH grant CA166703.
Nanoparticle self-assembly has been well studied theoretically, but it remains challenging to directly observe and quantify individual nanoparticle interactions. With our custom image analysis method, we track the trajectories of nanoparticle movement with high precision from a stack of relatively noisy images obtained using liquid cell transmission electron microscopy. In a time frame of minutes, Pt-Fe nanoparticles self-assembled into a loosely packed hcp lattice. The energetics and stability of the dynamic assembly were studied quantitatively. From velocity and diffusion measurements, we experimentally determined the magnitude of forces between single particles and the related physical properties. The results illustrate that long-range anisotropic forces drive the formation of chains, which then clump and fold to maximize close range van der Waals interactions.
Abstract Improving the spatial and spectral resolution of 2D X-ray near-edge absorption structure (XANES) has been a decade-long pursuit to probe local chemical reactions at the nanoscale. However, the poor signal-to-noise ratio in the measured images poses significant challenges in quantitative analysis, especially when the element of interest is at a low concentration. In this work, we developed a post-imaging processing method using deep neural network to reliably improve the signal-to-noise ratio in the XANES images. The proposed neural network model could be trained to adapt to new datasets by incorporating the physical features inherent in the latent space of the XANES images and self-supervised to detect new features in the images and achieve self-consistency. Two examples are presented in this work to illustrate the model’s robustness in determining the valence states of Ni and Co in the LiNi x Mn y Co 1-x-y O 2 systems with high confidence.
Purpose: Noise characteristics of clinical multidetector CT (MDCT) systems can be quantified by the noise power spectrum (NPS). Although the NPS of CT has been extensively studied in the past few decades, the joint impact of the bowtie filter and object position on the NPS has not been systematically investigated. This work studies the interplay of these two factors on the two dimensional (2D) local NPS of a clinical CT system that uses the filtered backprojection algorithm for image reconstruction. Methods: A generalized NPS model was developed to account for the impact of the bowtie filter and image object location in the scan field-of-view (SFOV). For a given bowtie filter, image object, and its location in the SFOV, the shape and rotational symmetries of the 2D local NPS were directly computed from the NPS model without going through the image reconstruction process. The obtained NPS was then compared with the measured NPSs from the reconstructed noise-only CT images in both numerical phantom simulation studies and experimental phantom studies using a clinical MDCT scanner. The shape and the associated symmetry of the 2D NPS were classified by borrowing the well-known atomic spectral symbols s, p, and d, which correspond to circular, dumbbell, and cloverleaf symmetries, respectively, of the wave function of electrons in an atom. Finally, simulated bar patterns were embedded into experimentally acquired noise backgrounds to demonstrate the impact of different NPS symmetries on the visual perception of the object. Results: (1) For a central region in a centered cylindrical object, an s-wave symmetry was always present in the NPS, no matter whether the bowtie filter was present or not. In contrast, for a peripheral region in a centered object, the symmetry of its NPS was highly dependent on the bowtie filter, and both p-wave symmetry and d-wave symmetry were observed in the NPS. (2) For a centered region-ofinterest (ROI) in an off-centered object, the symmetry of its NPS was found to be different from that of a peripheral ROI in the centered object, even when the physical positions of the two ROIs relative to the isocenter were the same. (3) The potential clinical impact of the highly anisotropic NPS, caused by the interplay of the bowtie filter and position of the image object, was highlighted in images of specific bar patterns oriented at different angles. The visual perception of the bar patterns was found to be strongly dependent on their orientation. Conclusions: The NPS of CT depends strongly on the bowtie filter and object position. Even if the location of the ROI with respect to the isocenter is fixed, there can be different symmetries in the NPS, which depend on the object position and the size of the bowtie filter. For an isolated off-centered object, the NPS of its CT images cannot be represented by the NPS measured from a centered object.
Purpose: This study evaluated the efficiency of 4D lung radiation treatment planning using Monte Carlo simulation on the cloud. The EGSnrc Monte Carlo code was used in dose calculation on the 4D-CT image set. Methods: 4D lung radiation treatment plan was created by the DOSCTP linked to the cloud, based on the Amazon elastic compute cloud platform. Dose calculation was carried out by Monte Carlo simulation on the 4D-CT image set on the cloud, and results were sent to the FFD4D image deformation program for dose reconstruction. The dependence of computing time for treatment plan on the number of compute node was optimized with variations of the number of CT image set in the breathing cycle and dose reconstruction time of the FFD4D. Results: It is found that the dependence of computing time on the number of compute node was affected by the diminishing return of the number of node used in Monte Carlo simulation. Moreover, the performance of the 4D treatment planning could be optimized by using smaller than 10 compute nodes on the cloud. The effects of the number of image set and dose reconstruction time on the dependence of computing time on the number of node were not significant, as more than 15 compute nodes were used in Monte Carlo simulations. Conclusion: The issue of long computing time in 4D treatment plan, requiring Monte Carlo dose calculations in all CT image sets in the breathing cycle, can be solved using the cloud computing technology. It is concluded that the optimized number of compute node selected in simulation should be between 5 and 15, as the dependence of computing time on the number of node is significant.
‘Radiomics’ refers to studies that extract a large amount of quantitative information from medical imaging studies as a basis for characterizing a specific aspect of patient health. Radiomics models can be built to address a wide range of outcome predictions, clinical decisions, basic cancer biology, etc. For example, radiomics models can be built to predict the aggressiveness of an imaged cancer, cancer gene expression characteristics (radiogenomics), radiation therapy treatment response, etc. Technically, radiomics brings together quantitative imaging, computer vision/image processing, and machine learning. In this symposium, speakers will discuss approaches to radiomics investigations, including: longitudinal radiomics, radiomics combined with other biomarkers (‘pan-omics’), radiomics for various imaging modalities (CT, MRI, and PET), and the use of registered multi-modality imaging datasets as a basis for radiomics. There are many challenges to the eventual use of radiomics-derived methods in clinical practice, including: standardization and robustness of selected metrics, accruing the data required, building and validating the resulting models, registering longitudinal data that often involve significant patient changes, reliable automated cancer segmentation tools, etc. Despite the hurdles, results achieved so far indicate the tremendous potential of this general approach to quantifying and using data from medical images. Specific applications of radiomics to be presented in this symposium will include: the longitudinal analysis of patients with low-grade gliomas; automatic detection and assessment of patients with metastatic bone lesions; image-based monitoring of patients with growing lymph nodes; predicting radiotherapy outcomes using multi-modality radiomics; and studies relating radiomics with genomics in lung cancer and glioblastoma. Learning Objectives: Understanding the basic image features that are often used in radiomic models. Understanding requirements for reliable radiomic models, including robustness of metrics, adequate predictive accuracy, and generalizability. Understanding the methodology behind radiomic-genomic (’radiogenomics’) correlations. Research supported by NIH (US), CIHR (Canada), and NSERC (Canada)
A critical factor in the implementation of direct energy deposition is the ability to maintain the standoff distance between the nozzle and the build surface, as this influences powder capture efficiency and overall part quality. Due to process-related variations, layer height may vary, causing unintended variation in standoff distance and poor build quality. While prior work has utilized contact probing to qualify standoff distance during processing, in situ methods for qualification of standoff distance are of major interest. The present work seeks to understand efficacy of image-based methods for classifying standoff distance variation in real-time using support vector machines (SVMs). It was hypothesized that the size of the melt pool and the amount of spatter will have significant correlations with deviations in the standoff distance; thus, SVMs were used on a dataset that is comprised of morphological features of melt pool size and image entropy. The SVM model was used to classify melt pool images into categories according to standoff distance variation from nominal. K-folds cross validation was used to find the optimal hyperparameters for the SVM model. To understand the impact of the selected features on the classification performance and inference speed, multiple models were trained with differing numbers of included features. Results for classification score, inference time, and image preprocessing/feature extraction from these data are reported. The present results show that the SVM model was able to predict the standoff distance classification with an accuracy of 97 percent and a speed of 0.122 s per image, making it a viable solution for real-time control of standoff distance.
Purpose: To test a method to reconstruct a four-dimensional (4D) dose distribution using the correlation of pre-calculated 4D electronic portal imaging device (EPID) images and measured cine-EPID images. Methods: 1. A phantom designed to simulate a tumor in lung (a polystyrene block with 3.0 cm diameter embedded in cork) was placed on a sinusoidally moving platform with 2 cm amplitude and 4 sec/cycle. Ten-phase 4D CT images were acquired for treatment planning and dose reconstruction. A 6MV photon beam was irradiated on the phantom with static (field size=5×8.5 cm{sup 2}) and dynamic fields (sliding windows, 10×10 cm{sup 2}, X1 MLC closing in parallel with the tumor movement). 2. 4D and 3D doses were calculated forwardly on PTV (1 cm margin). 3. Dose images on EPID under the fields were calculated for 10 phases. 4. Cine EPID images were acquired during irradiation. 5. Their acquisition times were correlated to the phases of the phantom at which irradiation occurred by inter-comparing calculated “reference” EPID images with measured images (2D gamma comparison). For the dynamic beam, the tumor was hidden under MLCs during a portion of irradiation time; the correlation performed when the tumor was visible was extrapolated. 6. Dose for each phase was reconstructed on the 4D CT images and summed over all phases. The summation was compared with forwardly calculated 4D and 3D dose distributions. Monte Carlo methods were used for all calculations. Results: For the open and dynamic beams, the 4D reconstructed doses showed the pass rates of 92.7 % and 100 %, respectively, at the isocenter plane given 3% / 3 mm criteria. The better agreement of the dynamic beam was from its dose gradient which blurred the otherwise sharp difference between forward and reconstructed doses. This also contributed slightly better agreement in DVH of PTV. Conclusion: The feasibility of 4D reconstruction was demonstrated.
Abstract Background The findings of the 2023 AAPM Grand Challenge on Deep Generative Modeling for Learning Medical Image Statistics are reported in this Special Report. Purpose The goal of this challenge was to promote the development of deep generative models for medical imaging and to emphasize the need for their domain‐relevant assessments via the analysis of relevant image statistics. Methods As part of this Grand Challenge, a common training dataset and an evaluation procedure was developed for benchmarking deep generative models for medical image synthesis. To create the training dataset, an established 3D virtual breast phantom was adapted. The resulting dataset comprised about 108 000 images of size 512 512. For the evaluation of submissions to the Challenge, an ensemble of 10 000 DGM‐generated images from each submission was employed. The evaluation procedure consisted of two stages. In the first stage, a preliminary check for memorization and image quality (via the Fréchet Inception Distance [FID]) was performed. Submissions that passed the first stage were then evaluated for the reproducibility of image statistics corresponding to several feature families including texture, morphology, image moments, fractal statistics, and skeleton statistics. A summary measure in this feature space was employed to rank the submissions. Additional analyses of submissions was performed to assess DGM performance specific to individual feature families, the four classes in the training data, and also to identify various artifacts. Results Fifty‐eight submissions from 12 unique users were received for this Challenge. Out of these 12 submissions, 9 submissions passed the first stage of evaluation and were eligible for ranking. The top‐ranked submission employed a conditional latent diffusion model, whereas the joint runners‐up employed a generative adversarial network, followed by another network for image superresolution. In general, we observed that the overall ranking of the top 9 submissions according to our evaluation method (i) did not match the FID‐based ranking, and (ii) differed with respect to individual feature families. Another important finding from our additional analyses was that different DGMs demonstrated similar kinds of artifacts. Conclusions This Grand Challenge highlighted the need for domain‐specific evaluation to further DGM design as well as deployment. It also demonstrated that the specification of a DGM may differ depending on its intended use.
The research effort described in this report represents a first attempt to investigate the radial redistribution of Zr in ternary fuel alloys U-xPu-10Zr (x = 0, 8, 19) irradiated in the in-reactor fuel experiments during the operation of the U.S. Department of Energy’s (DOE’s) Experimental Breeder Reactor II (EBR-II) at Idaho National Laboratory (INL) using methods of image processing on post-irradiation neutron radiographs. Approximately 130,000 metal fuel pins were irradiated in EBR-II during its 30 years of operation to develop and characterize existing and prospective fuels. For many of the metal fuel irradiation experiments, neutron radiography imaging was performed historically now allowing application of modern image analysis techniques to characterize fuel behavior, such as fuel swelling, fluff formation, and now, fuel alloy constituent redistribution. The redistribution of fuel components depends on the temperature field, radially, within the fuel. Specifically, Zr is expected to redistribute radially towards the center of the pin, as well as towards the outer zones. Currently, direct imaging of a pin cross-section through optical methods or scanning electron microscopy (SEM) is used to study the constituents’ redistribution, which is very time-consuming and can only be applied to a limited number of pins. An automated image processing technique allowing for the investigation of fuel radial redistribution zones would significantly accelerate data analysis. In general, if the fuel temperature is hot enough, three redistribution zones are expected, corresponding to the main fuel components (e.g., U, Pu, Zr). While more assessment will be performed in fiscal year (FY)-2023, it seems possible to differentiate the pins according to their fuel composition using image analysis techniques on neutron radiographs of metallic fuel pins based on the analysis to date.
We have developed a method for image reconstruction and analysis of dual-energy X-ray scans that has been demonstrated to yield a physically intuitive characterization of materials (by electron density and effective atomic number) in a system-independent way. Our methodology stands in contrast to common practice of reconstructing images to yield linear attenuation coefficients, which depend strongly on system spectral response. At the core of our method, a technique called dual-energy decomposition is used to generate spectral energy-independent characterizations of the attenuation properties of a material. Inputs to this decomposition include system spectral response estimates and dual-energy X-ray projection data. X-ray cross-section tables are used to convert decomposition results to physical feature estimates (electron density, effective atomic number). Our method currently works best for materials with atomic numbers below 22 (titanium). In prior work for DHS, we demonstrated in a controlled, no-clutter environment that the decomposition technique could be applied to scans of a well-characterized materials on multiple systems with different X-ray spectral responses. In cases where LAC values showed a standard deviation of as much as 20% of mean value, application of our method yielded results with standard deviations as low as 1–3% of mean value.
NASA's planetary missions have collected, and continue to collect, massive volumes of orbital imagery. The volume is such that it is difficult to manually review all of the data and determine its significance. As a result, images are indexed and searchable by location and date but generally not by their content. A new automated method analyzes images and identifies "landmarks," or visually salient features such as gullies, craters, dust devil tracks, and the like. This technique uses a statistical measure of salience derived from information theory, so it is not associated with any specific landmark type. It identifies regions that are unusual or that stand out from their surroundings, so the resulting landmarks are context-sensitive areas that can be used to recognize the same area when it is encountered again. A machine learning classifier is used to identify the type of each discovered landmark. Using a specified window size, an intensity histogram is computed for each such window within the larger image (sliding the window across the image). Next, a salience map is computed that specifies, for each pixel, the salience of the window centered at that pixel. The salience map is thresholded to identify landmark contours (polygons) using the upper quartile of salience values. Descriptive attributes are extracted for each landmark polygon: size, perimeter, mean intensity, standard deviation of intensity, and shape features derived from an ellipse fit.
I describe a systematic method for characterizing the effects of differences in u-v coverage in terms of inferred spectral gradients in interferometric images. This method is directly applicable to optical and radio interferometry; it is useful in situations when observations using scale arrays are not feasible.
Electrostatic modes of atomic force microscopy have shown to be non-destructive and relatively simple methods for imaging conductors embedded in insulating polymers. Here we use electrostatic force microscopy to image the dispersion of carbon nanotubes in a latex-based conductive composite, which brings forth features not observed in previously studied systems employing linear polymer films. A fixed-potential model of the probe-nanotube electrostatics is presented which in principle gives access to the conductive nanoparticle's depth and radius, and the polymer film dielectric constant. Comparing this model to the data results in nanotube depths that appear to be slightly above the film–air interface. Furthermore, this result suggests that water-mediated charge build-up at the film–air interface may be the source of electrostatic phase contrast in ambient conditions.