AutoCT: Automated CT registration, segmentation, and quantification
Not Available
SEARCH · Engineering Papers
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.
Not Available
We propose a new method to infer road networks from GPS trace data and accurately segment road regions in high-resolution aerial images. Unlike previous efforts that rely on GPS traces alone, we exploit image features to infer road networks from noisy trace data. The inferred road network is used to guide road segmentation. We show that the number of image segments spanned by the traces and the trace orientation validated with image features are important attributes for identifying GPS traces on road regions. Based on filtered traces , we construct road networks and integrate them with image features to segment road regions. Lastly, our experiments show that the proposed method produces more accurate road networks than the leading method that uses GPS traces alone, and also achieves high accuracy in segmenting road regions even with very noisy GPS data.
Abstract not provided
We propose a new method for joint segmentation of monotonously growing or shrinking shapes in a time sequence of noisy images. The task of segmenting the image time series is expressed as an optimization problem using the spatio-temporal graph of pixels, in which we are able to impose the constraint of shape growth or of shrinkage by introducing monodirectional infinite links connecting pixels at the same spatial locations in successive image frames. The globally optimal solution is computed with a graph cut. The performance of the proposed method is validated on three applications: segmentation of melting sea ice floes and of growing burned areas from time series of 2D satellite images, and segmentation of a growing brain tumor from sequences of 3D medical scans. In the latter application, we impose an additional intersequences inclusion constraint by adding directed infinite links between pixels of dependent image structures.
The Extravehicular Activity Helper/Retriever (EVAHR) is a robotic device currently under development at the NASA Johnson Space Center that is designed to fetch objects or to assist in retrieving an astronaut who may have become inadvertently de-tethered. The EVAHR will be required to exhibit a high degree of intelligent autonomous operation and will base much of its reasoning upon information obtained from one or more three-dimensional sensors that it will carry and control. At the highest level of visual cognition and reasoning, the EVAHR will be required to detect objects, recognize them, and estimate their spatial orientation and location. The recognition phase and estimation of spatial pose will depend on the ability of the vision system to reliably extract geometric features of the objects such as whether the surface topologies observed are planar or curved and the spatial relationships between the component surfaces. In order to achieve these tasks, three-dimensional sensing of the operational environment and objects in the environment will therefore be essential. One of the sensors being considered to provide image data for object recognition and pose estimation is a phase-shift laser scanner. The characteristics of the data provided by this scanner have been studied and algorithms have been developed for segmenting range images into planar surfaces, extracting basic features such as surface area, and recognizing the object based on the characteristics of extracted features. Also, an approach has been developed for estimating the spatial orientation and location of the recognized object based on orientations of extracted planes and their intersection points. This paper presents some of the algorithms that have been developed for the purpose of recognizing and estimating the pose of objects as viewed by the laser scanner, and characterizes the desirability and utility of these algorithms within the context of the scanner itself, considering data quality and noise.
Purpose: To investigate the accuracy and potential limitations of MV image-based dynamic multileaf collimator (DMLC) tracking in a porcine model on a linear accelerator. Methods and Materials: A thermo-expandable NiTi stent designed for kilovoltage (kV) X-ray visualization of lung lesions was inserted into the bronchia of three anaesthetized Goettingen minipigs. A four-dimensional computed tomography scan was used for planning a five-field conformal treatment with circular multileaf collimator (MLC) apertures. A 22.5 Gy single fraction treatment was delivered to the pigs. The peak-to-peak stent motion was 3 to 8 mm, with breathing periods of 1.2 to 4 s. Before treatment, X-ray images were used for image-guided setup based on the stent. During treatment delivery, continuous megavoltage (MV) portal images were acquired at 7.5 Hz. The stent was segmented in the images and used for continuous adaptation of the MLC aperture. Offline, the tracking error in beam's eye view of the treatment beam was calculated for each MV image as the difference between the MLC aperture center and the segmented stent position. The standard deviations of the systematic error {Sigma} and the random error {sigma} were determined and compared with the would-be errors for a nontracking treatment with pretreatment image-guided setup. Results: Reliable stent segmentation was obtained for 11 of 15 fields. Segmentation failures occurred when image contrast was dominated by overlapping anatomical structures (ribs, diaphragm) rather than by the stent, which was designed for kV rather than MV X-ray visibility. For the 11 fields with reliable segmentation, {Sigma} was 0.5 mm/0.4 mm in the two imager directions, whereas {sigma} was 0.5 mm/1.1 mm. Without tracking, {Sigma} and {sigma} would have been 1.7 mm/1.4 mm and 0.8 mm/1.4 mm, respectively. Conclusion: For the first time, in vivo DMLC tracking has been demonstrated on a linear accelerator showing the potential for improved targeting accuracy. The study mimicked the envisioned patient workflow of future patient treatments. Clinical implementation of MV image-based tracking would require markers designed for MV visibility.
Purpose: Investigating a new approach in MRI based treatment planning using the combination of (Ultrashort Echo Time) UTE and T1 weighted spin echo pulse sequences to delineate air, bone and water (soft tissues) in generating pseudo CT images comparable with CT. Methods: A gel phantom containing chicken bones, ping pang balls filled with distilled water and air bubbles, was made. It scanned with MRI using UTE and 2D T1W SE pulse sequences with (in plane resolution= 0.53mm, slice thickness= 2 mm) and CT with (in plane resolution= 0.5 mm and slice thickness= 0.75mm) as a ground truth for geometrical accuracy. The UTE and T1W SE images were registered with CT using mutual information registration algorithm provided by Philips Pinnacle treatment planning system. The phantom boundaries were detected using Canny edge detection algorithm for CT, and MR images. The bone, air bubbles and water in ping pong balls were segmented from CT images using threshold 300HU, - 950HU and 0HU, respectively. These tissue inserts were automatically segmented from combined UTE and T1W SE images using edge detection and relative intensity histograms of the phantom. The obtained segmentations of air, bone and water inserts were evaluated with those obtained from CT. Results: Bone and air can be clearly differentiated in UTE images comparable to CT. Combining UTE and T1W SE images successfully segmented the air, bone and water. The maximum segmentation differences from combine MRI images (UTE and T1W SE) and CT are within 1.3 mm, 1.1mm for bone, air, respectively. The geometric distortion of UTE sequence is small less than 1 pixel (0.53 mm) of MR image resolution. Conclusion: Our approach indicates that MRI can be used solely for treatment planning and its quality is comparable with CT.
At the former nuclear weapon production site in Hanford, WA, caustic radioactive tank waste leaks into subsurface sediments and causes dissolution of quartz and aluminosilicate minerals, and precipitation of sodalite and cancrinite. This work examines changes in pore structure due to these reactions in a previously-conducted column experiment. The column was sectioned and 2D images of the pore space were generated using backscattered electron microscopy and energy dispersive X-ray spectroscopy. A pre-precipitation scenario was created by digitally removing mineral matter identified as secondary precipitates. Porosity, determined by segmenting the images to distinguish pore space from mineral matter, was up to 0.11 less after reaction. Erosion-dilation analysis was used to compute pore and throat size distributions. Images with precipitation had more small and fewer large pores. Precipitation decreased throat sizes and the abundance of large throats. These findings agree with previous findings based on 3D X-ray CMT imaging, observing decreased porosity, clogging of small throats, and little change in large throats. However, 2D imaging found an increase in small pores, mainly in intragranular regions or below the resolution of the 3D images. Also, an increase in large pores observed via 3D imaging was not observed in the 2D analysis. Changes in flow conducting throats that are the key permeability-controlling features were observed in both methods.
Here, the microstructure of plastic bonded explosives (PBXs) is known to influence behavior during mechanical deformation, but characterizing the microstructure can be challenging. For example, the explosive crystals and binder in formulations such as PBX 9501 do not have sufficient X-ray contrast to obtain three-dimensional data by in situ, absorption contrast imaging. To address this difficulty, we have formulated a series of PBXs using octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) crystals and low-density binder systems. The binders were hydroxyl-terminated polybutadiene (HTPB) or glycidyl azide polymer (GAP) cured with a commercial blend of acrylic monomers/oligomers. The binder density is approximately half of the HMX, allowing for excellent contrast using in situ X-ray computed tomography (CT) imaging. The samples were imaged during unaxial compression using micro-scale CT in an interrupted in situ modality. The rigidity of the binder was observed to significantly influence fracture, crystal-binder delamination, and flow. Additionally, 2D slices from the segmented 3D images were meshed for finite element simulation of the mesoscale response. At low stiffness, the binder and crystal do not delaminate and the crystals move with the material flow; at high stiffness, marked delamination is noted between the crystals and the binder, leading to very different mechanical properties. Initial model results exhibit qualitatively similar delamination.
Post processing X-ray computational tomography (CT) inspection data for additively manufactured (AM) components can induce deviations in defect quantification, affecting subsequent fatigue and failure predictions. To assess the influence and potential impact of segmentation-induced measurement deviations, this paper applies several segmentation techniques to X-ray CT data for powder bed fusion Ti-6Al-4V specimens exhibiting porosity conditions. X-ray CT reconstructions were segmented with varying techniques including Otsu’s thresholding, random forest, k-nearest neighbors, and the multilayer perceptron. Metrics such as pore size and global porosity were compared for internal validity. Then, top-down X-ray CT measurements of surface-breaking porosity were compared to optical profilometry for cross-validation.
Explore the source record for details and available documents.
This dataset is associated with a broader study using game camera timeseries photos collected to evaluate stream variable inundation via changes in width (i.e. wet fraction). Four game cameras were deployed along Umtanum Creek (Washington, United States) to track changes in stream inundation over time. Drone imagery was collected at the same location on October 18, 2024 which was used to construct a digital elevation model (DEM) of the streambed topography. The associated paper and data can be found at https://doi.org/10.1016/j.envsoft.2025.106715 (Bao et al., 2025a)) and https://doi.org/10.15485/2589885 (Bao et al., 2025b), respectively. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to this readme, this data package also includes a file-level metadata (FLMD) files that describes each file and a data dictionaries (DD) that describe all column/row headers and variable definitions. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) readme; (4) field metadata; (5) field protocol; and (5) folders containing game camera photos. Game camera photos are organized into folders for each camera (CDL, CUL, CDR, CUR; see readme for information on camera naming) by the month photos were collected. All files are .csv, .jpg, or .pdf.
Experiments on clinical datasets of H&N patients demonstrated the effectiveness of the proposed deep neural network segmentation method for multi-organ segmentation on volumetric CT scans. Additionally, the accuracy and robustness of the segmentation were further increased by incorporating shape priors using SMR. The proposed method showed competitive performance and took shorter time to segment multiple organs in comparison to state of the art methods.
Monte Carlo simulation methods are widely used in medical physics research and are starting to be implemented in clinical applications such as radiation therapy planning systems. Monte Carlo simulations offer the capability to accurately estimate quantities of interest that are challenging to measure experimentally while taking into account the realistic anatomy of an individual patient. Traditionally, practical application of Monte Carlo simulation codes in diagnostic imaging was limited by the need for large computational resources or long execution times. However, recent advancements in high-performance computing hardware, combined with a new generation of Monte Carlo simulation algorithms and novel postprocessing methods, are allowing for the computation of relevant imaging parameters of interest such as patient organ doses and scatter-to-primaryratios in radiographic projections in just a few seconds using affordable computational resources. Programmable Graphics Processing Units (GPUs), for example, provide a convenient, affordable platform for parallelized Monte Carlo executions that yield simulation times on the order of 10{sup 7} xray/ s. Even with GPU acceleration, however, Monte Carlo simulation times can be prohibitive for routine clinical practice. To reduce simulation times further, variance reduction techniques can be used to alter the probabilistic models underlying the x-ray tracking process, resulting in lower variance in the results without biasing the estimates. Other complementary strategies for further reductions in computation time are denoising of the Monte Carlo estimates and estimating (scoring) the quantity of interest at a sparse set of sampling locations (e.g. at a small number of detector pixels in a scatter simulation) followed by interpolation. Beyond reduction of the computational resources required for performing Monte Carlo simulations in medical imaging, the use of accurate representations of patient anatomy is crucial to the virtual generation of medical images and accurate estimation of radiation dose and other imaging parameters. For this, detailed computational phantoms of the patient anatomy must be utilized and implemented within the radiation transport code. Computational phantoms presently come in one of three format types, and in one of four morphometric categories. Format types include stylized (mathematical equation-based), voxel (segmented CT/MR images), and hybrid (NURBS and polygon mesh surfaces). Morphometric categories include reference (small library of phantoms by age at 50th height/weight percentile), patient-dependent (larger library of phantoms at various combinations of height/weight percentiles), patient-sculpted (phantoms altered to match the patient's unique outer body contour), and finally, patient-specific (an exact representation of the patient with respect to both body contour and internal anatomy). The existence and availability of these phantoms represents a very important advance for the simulation of realistic medical imaging applications using Monte Carlo methods. New Monte Carlo simulation codes need to be thoroughly validated before they can be used to perform novel research. Ideally, the validation process would involve comparison of results with those of an experimental measurement, but accurate replication of experimental conditions can be very challenging. It is very common to validate new Monte Carlo simulations by replicating previously published simulation results of similar experiments. This process, however, is commonly problematic due to the lack of sufficient information in the published reports of previous work so as to be able to replicate the simulation in detail. To aid in this process, the AAPM Task Group 195 prepared a report in which six different imaging research experiments commonly performed using Monte Carlo simulations are described and their results provided. The simulation conditions of all six cases are provided in full detail, with all necessary data on material composition, source, geometry, scoring and other parameters provided. The results of these simulations when performed with the four most common publicly available Monte Carlo packages are also provided in tabular form. The Task Group 195 Report will be useful for researchers needing to validate their Monte Carlo work, and for trainees needing to learn Monte Carlo simulation methods. In this symposium we will review the recent advancements in highperformance computing hardware enabling the reduction in computational resources needed for Monte Carlo simulations in medical imaging. We will review variance reduction techniques commonly applied in Monte Carlo simulations of medical imaging systems and present implementation strategies for efficient combination of these techniques with GPU acceleration. Trade-offs involved in Monte Carlo acceleration by means of denoising and “sparse sampling” will be discussed. A method for rapid scatter correction in cone-beam CT (<5 min/scan) will be presented as an illustration of the simulation speeds achievable with optimized Monte Carlo simulations. We will also discuss the development, availability, and capability of the various combinations of computational phantoms for Monte Carlo simulation of medical imaging systems. Finally, we will review some examples of experimental validation of Monte Carlo simulations and will present the AAPM Task Group 195 Report. Learning Objectives: Describe the advances in hardware available for performing Monte Carlo simulations in high performance computing environments. Explain variance reduction, denoising and sparse sampling techniques available for reduction of computational time needed for Monte Carlo simulations of medical imaging. List and compare the computational anthropomorphic phantoms currently available for more accurate assessment of medical imaging parameters in Monte Carlo simulations. Describe experimental methods used for validation of Monte Carlo simulations in medical imaging. Describe the AAPM Task Group 195 Report and its use for validation and teaching of Monte Carlo simulations in medical imaging.
Continuing on the successful progress of NIAC (NASA Innovative Advanced Concepts) Phase I, this report summarizes the technical progress achieved under NIAC Phase II during the performance period September 19, 2014 to June 18, 2017. During this period, the research team has made the following accomplishments: designed and layout a silica photonic integrated circuit (PIC) as a two baselineinterferometric imager; constructed an experiment to utilize the two baselines for complex visibility measurementon a point source and a variable width slit; analyzed and studied the testbed results (in collaboration with Lockheed Martin); designed and layout Si3N4 PICs for the low-resolution and high-resolution SPIDER (Segmented Planar Imaging Detector for Electro-Optical Reconnaissance) telescope; fabricated the multi-layer Si3N4 PIC for low and high resolution SPIDER telescope; characterized the optical throughput and heater response for Si3N4 PIC for low and highresolution SPIDER telescopes; carried out imaging experiments using the Si3N4 PIC low-resolution version (in collaboration with Lockheed Martin); investigated signal-to-noise (SNR) ratio of SPIDER imager compared to the conventional panchromatic imager (in collaboration with Lockheed Martin); fulfilled the SNR simulation upon SPIDER imager (in collaboration with Lockheed Martin).
A computer-assisted method for localizing a rack, including sensing an image of the rack, detecting line segments in the sensed image, recognizing a candidate arrangement of line segments in the sensed image indicative of a predetermined feature of the rack, generating a matrix of correspondence between the candidate arrangement of line segments and an expected position and orientation of the predetermined feature of the rack, and estimating a position and orientation of the rack based on the matrix of correspondence.
The quantification accuracy of dual-energy imaging is influenced by the fundamentals of x-ray physics, system geometry, data acquisition hardware/protocol, system calibration, and image processing technique. This symposium will provide updates on the following advanced application areas: Mammography. Volumetric breast density techniques based on standard mammograms require estimation of breast thickness, which is difficult to accurately measure. By comparison, calculation of breast density using dual energy mammography does not require measurement of breast thickness. Dual energy mammography has been implemented using both energy integrating flat panel detectors in conjunction with beam energy switching and energy resolved photon counting detectors. These techniques have been optimized using simulation studies and validated using physical phantoms and postmortem breasts. Chemical decomposition was used as the gold standard for volumetric breast density measurement in postmortem breasts. Breast density measurements have also been compared with results from four-category BI-RADS density rankings, standard image thresholding and Fuzzy k-mean clustering techniques. These studies indicate that dual energy mammography can be used to accurately measure volumetric breast density. Cardiovascular CT. The predicative accuracy of risk models for recurrent stroke and cardiac arrest depends heavily on accurate differentiation of thrombus or calcium from iodine in left atrial appendage or coronary arteries. The amount of energy separation is constrained by image noise; therefore, optimal kVp, beam filtration, and balanced flux are essential for the quantification accuracy of iodine and calcium. The basis materials are combined linearly to generate monochromatic energy images, where CT# accuracy and CNR are energy dependent. With optimal monochromatic energy, the mean iodine concentration for the thrombus, circulatory stasis, and control groups are significantly different. Risk classification based on calcium scores shows excellent agreement with classification on the basis of conventional coronary artery calcium scoring. These studies demonstrate dual-energy cardiovascular CT can potentially be a noninvasive and sensitive modality in high risk patients. On-board KV/MV Imaging. To enhance soft tissue contrast and reduce metal artifacts, we have developed a dual-energy CBCT technique and a novel on-board kV/MV imaging technique based on hardware available on modern linear accelerators. We have also evaluated the feasibility of these two techniques in various phantom studies. Optimal techniques (energy, beam filtration, # of overlapping projections, etc) have been investigated with unique calibration procedures, which leads to successful decomposition of imaged material into acrylic-aluminum basis material pair. This enables the synthesis of virtual monochromatic (VM) CBCT images that demonstrate much less beam hardening, significantly reduced metal artifacts, and/or higher soft tissue CNR compared to single-energy CBCT. Adaptive Radiation Therapy. DECT could actually contribute to the area of Dose-Guided Radiation Therapy (or Adaptive Therapy). The application of DECT imaging using 80kV and 140 kV combinations could potentially increase the image quality by reducing the bone or high density material artifacts and also increase the soft tissue contrast by a light contrast agent. The result of this higher contrast / quality images is beneficial for deformable image registration / segmentation algorithm to improve its accuracy hence to make adaptive therapy less time consuming in its recontouring process. The real time re-planning prior to per treatment fraction could become more realistic with this improvement especially in hypofractional SBRT cases. Learning Objectives: Learn recent developments of dual-energy imaging in diagnosis and radiation therapy; Understand the unique clinical problem and required quantification accuracy in each application; Understand the different approaches to optimize dual-energy imaging techniques for different applications.
Purpose: With the continuous development of new stent grafts and implantation techniques, it has now become technically feasible to treat abdominal aortic aneurysms (AAA) with challenging anatomy using endovascular repair with standard, fenestrated, or branched stent-grafts. In vitro experimentations are very useful to improve stent-graft design and conformability or imaging guidance for stent-graft delivery or follow-up. Vascular replicas also help to better understand the limitation of endovascular approaches in challenging anatomy and possibly improve surgical planning or training by practicing high risk clinical procedures in the laboratory to improve outcomes in the operating room. Most AAA phantoms available have a very basic anatomy, which is not representative of the clinical reality. This paper presents a method of fabrication of a realistic AAA phantom with a visible thrombus, as well as some mechanical properties characterizing such phantom. Methods: A realistic AAA geometry replica of a real patient anatomy taken from a multidetector computed tomography (CT) scan was manufactured. To demonstrate the multimodality imaging capability of this new phantom with a thrombus visible in magnetic resonance (MR) angiography, CT angiography (CTA), digital subtraction angiography (DSA), and ultrasound, image acquisitions with all these modalities were performed by using standard clinical protocols. Potential use of this phantom for stent deployment was also tested. A rheometer allowed defining hyperelastic and viscoelastic properties of phantom materials. Results: MR imaging measurements of SNR and CNR values on T1 and T2-weighted sequences and MR angiography indicated reasonable agreement with published values of AAA thrombus and abdominal components in vivo. X-ray absorption also lay within normal ranges of AAA patients and was representative of findings observed on CTA, fluoroscopy, and DSA. Ultrasound propagation speeds for developed materials were also in concordance with the literature for vascular and abdominal tissues. Conclusions: The mimicked abdominal tissues, AAA wall, and surrounding thrombus were developed to match imaging features of in vivo MR, CT, and ultrasound examinations. This phantom should be of value for image calibration, segmentation, and testing of endovascular devices for AAA endovascular repair.