SEARCH · Engineering Papers

Results for “3D segmentation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS2017

Ocular Coherence Tomography in the Evaluation of Anterior Eye Injuries in Space Flight

While Ocular Coherence Tomography (OCT) is not a first-line modality to evaluate anterior eye structures terrestrially, it is a resource already available on the International Space Station (ISS) that can be used in medical contingencies that involve the anterior eye. With remote guidance and subject matter expert (SME) support from the ground, a minimally trained crewmember can now use OCT to evaluate anterior eye pathologies on orbit. OCT utilizes low-coherence interferometry to produce detailed cross-sectional and 3D images of the eye in real time. Terrestrially, it has been used to evaluate macular pathologies and glaucoma. Since 2013, OCT has been used onboard the ISS as one part of a suite of hardware to evaluate the Visual Impairment/Intracranial Pressure risk faced by astronauts, specifically assessing changes in the retina and choroid during space flight. The Anterior Segment Module (ASM), an add-on lens, was also flown for research studies, providing an opportunity to evaluate the anterior eye in real time if clinically indicated. Anterior eye pathologies that could be evaluated using OCT were identified. These included corneal abrasions and ulcers, scleritis, and acute angle closure glaucoma. A remote guider script was written to provide ground specialists with step-by-step instructions to guide ISS crewmembers, who do not get trained on the ASM, to evaluate the anterior eye. The instructions were tested on novice subjects and/or operators, whose feedback was incorporated iteratively. The final remote guider script was reviewed by SME optometrists and NASA flight surgeons. The novel application of OCT technology to space flight allows for the acquisition of objective data to diagnose anterior eye pathologies when other modalities are not available. This demonstrates the versatility of OCT and highlights the advantages of using existing hardware and remote guidance skills to expand clinical capabilities in space flight.

Fer, Dan M.↗

DOE OSTI2014

SU-E-J-196: New Visualization Methods for Longitudinal MRI Registrations and Segmentations

Purpose: To develop visualization techniques to facilitate easy assessment of (a) registration and (b) tracking volumetric changes in structures during radiation therapy from MRI. Method: The frequently used method for visualizing registrations between scans is a multi-color overlay technique or deformation vector fields. However, the overlay technique is unintuitive and does not help to appreciate the quality of registration particularly when the registration mismatches are not very large. Similarly, the deformation fields give an indication of extent of deformation but do not help to assess the differences in registration. We present a mirroring and edge-augmented mirroring technique that places the fixed and moving image next to each other and allows the user to quickly assess the small differences in registration. Next, we present a volumetric intersection based 3D model to visualize the changes in diseased lymph node volumes in head and neck cancer. 3D model-based visualization provides more information about the location-specific changes in volume rather than the simplistic one dimensional information obtained from 2D plot of nodal volume changes. Result: We show results comparing our approach with the standard colorbased overlay method for comparing registrations of intra-patient registrations using T2-MRI. Upon comparing the mirroring technique with the color-overlay, one can more easily appreciate the differences in registration. Adding edge-based mirroring seems to further assist in evaluating the registration. Our approach for viewing registrations seems to be more intuitive and easy to use in order to help assess the quality of registration compared to color-based overlays. Similarly, the change volumetric model together with a 2D plot reveals more information including the locations undergoing changes and responding to treatment. Conclusions: Better approaches are necessary for assessing the quality of registrations and changes in diseased structures over time. Our approach presents some novel visualization solutions for registration and segmentation.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

batteryNET v0.0.1

Code to organize, prepare, and segment multi-scan, larger-than-memory, volumetric dataset. Deep learning segmentation uses pytorch-lightning for code organization, logging and metrics, and multi-GPU parallelization of training and prediction. Also uses moani for volumetric data augmentation and model implementation (currently U-Net++). Training can be done slice-by-slice with 2D models trained on 2D patches either aligned with or perpendicular to the electrode plane, or with 3D models on 3D patches. Sparse labels can be used to minimize required hand-labeling -- less than 0.1% of the entire dataset was labeled. Also includes code for data preparation (alignment/cropping/patching).

Ushizima, Daniela↗

NASA NTRS2011

Visual Object Recognition and Tracking of Tools

A method has been created to automatically build an algorithm off-line, using computer-aided design (CAD) models, and to apply this at runtime. The object type is discriminated, and the position and orientation are identified. This system can work with a single image and can provide improved performance using multiple images provided from videos. The spatial processing unit uses three stages: (1) segmentation; (2) initial type, pose, and geometry (ITPG) estimation; and (3) refined type, pose, and geometry (RTPG) calculation. The image segmentation module files all the tools in an image and isolates them from the background. For this, the system uses edge-detection and thresholding to find the pixels that are part of a tool. After the pixels are identified, nearby pixels are grouped into blobs. These blobs represent the potential tools in the image and are the product of the segmentation algorithm. The second module uses matched filtering (or template matching). This approach is used for condensing synthetic images using an image subspace that captures key information. Three degrees of orientation, three degrees of position, and any number of degrees of freedom in geometry change are included. To do this, a template-matching framework is applied. This framework uses an off-line system for calculating template images, measurement images, and the measurements of the template images. These results are used online to match segmented tools against the templates. The final module is the RTPG processor. Its role is to find the exact states of the tools given initial conditions provided by the ITPG module. The requirement that the initial conditions exist allows this module to make use of a local search (whereas the ITPG module had global scope). To perform the local search, 3D model matching is used, where a synthetic image of the object is created and compared to the sensed data. The availability of low-cost PC graphics hardware allows rapid creation of synthetic images. In this approach, a function of orientation, distance, and articulation is defined as a metric on the difference between the captured image and a synthetic image with an object in the given orientation, distance, and articulation. The synthetic image is created using a model that is looked up in an object-model database. A composable software architecture is used for implementation. Video is first preprocessed to remove sensor anomalies (like dead pixels), and then is processed sequentially by a prioritized list of tracker-identifiers.

English, James↗

NASA NTRS2012

NASA Tech Briefs, December 2012

The topics include: Pattern Generator for Bench Test of Digital Boards; 670-GHz Down- and Up-Converting HEMT-Based Mixers; Lidar Electro-Optic Beam Switch with a Liquid Crystal Variable Retarder; Feedback Augmented Sub-Ranging (FASR) Quantizer; Real-Time Distributed Embedded Oscillator Operating Frequency Monitoring; Software Modules for the Proximity-1 Space Link Interleaved Time Synchronization (PITS) Protocol; Description and User Instructions for the Quaternion to Orbit v3 Software; AdapChem; Mars Relay Lander and Orbiter Overflight Profile Estimation; Extended Testability Analysis Tool; Interactive 3D Mars Visualization; Rapid Diagnostics of Onboard Sequences; MER Telemetry Processor; pyam: Python Implementation of YaM; Process for Patterning Indium for Bump Bonding; Archway for Radiation and Micrometeorite Occurrence Resistance; 4D Light Field Imaging System Using Programmable Aperture; Device and Container for Reheating and Sterilization; Radio Frequency Plasma Discharge Lamps for Use as Stable Calibration Light Sources; Membrane Shell Reflector Segment Antenna; High-Speed Transport of Fluid Drops and Solid Particles via Surface Acoustic Waves; Compact Autonomous Hemispheric Vision System; A Distributive, Non-Destructive, Real-Time Approach to Snowpack Monitoring; Wideband Single-Crystal Transducer for Bone Characterization; Numerical Simulation of Rocket Exhaust Interaction With Lunar Soil; Motion Imagery and Robotics Application (MIRA): Standards-Based Robotics; Particle Filtering for Model-Based Anomaly Detection in Sensor Networks; Ka-band Digitally Beamformed Airborne Radar Using SweepSAR Technique; Composite With In Situ Plenums; Multi-Beam Approach for Accelerating Alignment and Calibration of HyspIRI-Like Imaging Spectrometers; JWST Lifting System; Next-Generation Tumbleweed Rover; Pneumatic System for Concentration of Micrometer-Size Lunar Soil.

Source record↗

DOE OSTI2025

Designing Antifouling and Antimicrobial Interfaces: Structural Characterization using CryoEM, Automated Microscopy, and AI Image Segmentation

The design of functionalized surfaces for interactions with biological systems is critical across sectors such as healthcare, energy, and agriculture. Tailoring materials for specific applications, such as antifouling and antimicrobial surfaces, demands a comprehensive understanding of topology and chemistry across multiple length and time scales on both biological and materials systems. This work presents the development and characterization of nanostructured surfaces with controlled topographies and chemistries that enhance bacterial membrane disruption, reduce biofilm formation, and improve antimicrobial and antifouling capabilities. Two specific use cases will be presented - the use of cellulose nanocrystals (CNCs) for bacterial growth inhibition and the development of antifouling surfaces to prevent protein and bacterial adsorption [1-4]. By leveraging large language models (LLMs) for image segmentation and training [5], we enable automated analysis of terabyte-scale cryogenic electron microscopy (cryoEM) datasets. This analysis provides statistical insights into the biotic/abiotic interface and facilitates automated electron microscopy experiments to mitigate time and dose. The integration of cryogenic electron tomography (cryoET) and cryogenic focused ion beam (cryoFIB) milling enables high-resolution, near-native-state imaging and 3D reconstructions of bio/material interfaces [6]. Orthogonal characterization techniques and computational modeling further enhances our understanding, offering a robust platform for the design and optimization of next-generation functional surfaces [7].

Williams, Alexis [ORNL] (ORCID:0000000252835822)↗

DOE OSTI2018

Ultrasonic Evaluation of Segmental Variability in Additively Manufactured Metal Components

One of the primary challenges with additively manufactured metal components is the potential for the material heterogeneity, such as the variations in density, porosity, and elasticity. These heterogeneities are difficult to predict and can cause the deviation in material strength and structural reliability from the specified design. This paper presents a novel approach to the ultrasound characterization of spatially varying material properties and its application to a sample manufactured by the selective laser melting. The method is demonstrated for a 3D-printed aluminum alloy sample. The independent measurements of the elastic modulus, the surface morphology, and the X-ray computed tomography of the sample were consistent with the ultrasound results.

Roy, Manish↗

DOE OSTI2013

Upgrade and benchmarking of a 4D treatment planning system for scanned ion beam therapy

Purpose: Upgrade and benchmarking of a research 4D treatment planning system (4DTPS) suitable for realistic patient treatment planning and treatment simulations taking into account specific requirements for scanned ion beam therapy, i.e., modeling of dose heterogeneities due to interplay effects and range changes caused by patient motion and dynamic beam delivery. Methods: The 4DTPS integrates data interfaces to 4D computed tomography (4DCT), deformable image registration and clinically used motion monitoring devices. The authors implemented a novel data model for 4D image segmentation using Boolean mask volume datasets and developed an algorithm propagating a manually contoured reference contour dataset to all 4DCT phases. They further included detailed treatment simulation and dose reconstruction functionality, based on the irregular patient motion and the temporal structure of the beam delivery. The treatment simulation functionality was validated against experimental data from irradiation of moving radiographic films in air, 3D moving ionization chambers in a water phantom, and moving cells in a biological phantom with a scanned carbon ion beam. The performance of the program was compared to results obtained with predecessor programs. Results: The measured optical density distributions of the radiographic films were reproduced by the simulations to (-2 {+-} 12)%. Compared to earlier versions of the 4DTPS, the mean agreement improved by 2%, standard deviations were reduced by 7%. The simulated dose to the moving ionization chambers in water showed an agreement with the measured dose of (-1 {+-} 4)% for the typical beam configuration. The mean deviation of the simulated from the measured biologically effective dose determined via cell survival was (617 {+-} 538) mGy relative biological effectiveness corresponding to (10 {+-} 9)%. Conclusions: The authors developed a research 4DTPS suitable for realistic treatment planning on patient data and capable of simulating dose delivery to a moving patient geometry for scanned ion beams. The accuracy and reliability of treatment simulations improved considerably with respect to earlier versions of the 4DTPS.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-C-BRA-04: Automated Segmentation of Head-And-Neck CT Images for Radiotherapy Treatment Planning Via Multi-Atlas Machine Learning (MAML)

Purpose: Accurate image segmentation is a crucial step during image guided radiation therapy. This work proposes multi-atlas machine learning (MAML) algorithm for automated segmentation of head-and-neck CT images. Methods: As the first step, the algorithm utilizes normalized mutual information as similarity metric, affine registration combined with multiresolution B-Spline registration, and then fuses together using the label fusion strategy via Plastimatch. As the second step, the following feature selection strategy is proposed to extract five feature components from reference or atlas images: intensity (I), distance map (D), box (B), center of gravity (C) and stable point (S). The box feature B is novel. It describes a relative position from each point to minimum inscribed rectangle of ROI. The center-of-gravity feature C is the 3D Euclidean distance from a sample point to the ROI center of gravity, and then S is the distance of the sample point to the landmarks. Then, we adopt random forest (RF) in Scikit-learn, a Python module integrating a wide range of state-of-the-art machine learning algorithms as classifier. Different feature and atlas strategies are used for different ROIs for improved performance, such as multi-atlas strategy with reference box for brainstem, and single-atlas strategy with reference landmark for optic chiasm. Results: The algorithm was validated on a set of 33 CT images with manual contours using a leave-one-out cross-validation strategy. Dice similarity coefficients between manual contours and automated contours were calculated: the proposed MAML method had an improvement from 0.79 to 0.83 for brainstem and 0.11 to 0.52 for optic chiasm with respect to multi-atlas segmentation method (MA). Conclusion: A MAML method has been proposed for automated segmentation of head-and-neck CT images with improved performance. It provides the comparable result in brainstem and the improved result in optic chiasm compared with MA. Xuhua Ren and Hao Gao were partially supported by the NSFC (#11405105), the 973 Program (#2015CB856000), and the Shanghai Pujiang Talent Program (#14PJ1404500).

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

Rolling hills on the core-mantle boundary

Recent results suggest that an iron-rich oxide may have fractionally crystallized from a primordial magma ocean and settled on the core–mantle boundary (CMB). Based on experimental results, the presence of only a few percent of Fe-rich oxide could slow seismic waves down by several percent. This heavy layer can become highly undulating as predicted from dynamic modeling but can remain as a distinct structure with uniform velocity reductions. Here, we use the large USArray seismic network to search for such structures. Strong constraints on D" are provided by the core-phase SKS where it bifurcates, containing a short segment of P-wave diffractions (P d ) when crossing the CMB, called SKS d . Synthetics from models with moderate velocity drops (less than 10%) involving a layer with variable thickness, perhaps a composite of sharp small structures, with strong variation in thickness can explain both the observed SKS d waveforms and large scatter in differential times between SKKS and SKS. A smooth 3D image is obtained from inverting SKS d waveforms displaying rolling-hills with elongated dome-like structures sitting on the CMB. The most prominent one has an 80-km height, ~8° length, and ~4° width, thus adding still more structural complexity to the lower mantle. We suggest that these results can be explained by a dynamically-stabilized material containing small amounts (~5%) iron-rich (Mg,Fe)O providing a self-consistent physical interpretation.

Sun, Daoyuan↗

NASA NTRS1996

Meshing of a Spiral Bevel Gearset with 3D Finite Element Analysis

Recent advances in spiral bevel gear geometry and finite element technology make it practical to conduct a structural analysis and analytically roll the gearset through mesh. With the advent of user specific programming linked to 3D solid modelers and mesh generators, model generation has become greatly automated. Contact algorithms available in general purpose finite element codes eliminate the need for the use and alignment of gap elements. Once the gearset is placed in mesh, user subroutines attached to the FE code easily roll the gearset through mesh. The method is described in detail. Preliminary results for a gearset segment showing the progression of the contact lineload is given as the gears roll through mesh.

Bibel, George D.↗

DOE OSTI2011

Potential of discrete Gaussian edge feathering method for improving abutment dosimetry in eMLC-delivered segmented-field electron conformal therapy

Purpose: The purpose of this work was to investigate the potential of discrete Gaussian edge feathering of the higher energy electron fields for improving abutment dosimetry in the planning volume when using an electron multileaf collimator (eMLC) to deliver segmented-field electron conformal therapy (ECT). Methods: A discrete (five-step) Gaussian edge spread function was used to match dose penumbras of differing beam energies (6-20 MeV) at a specified depth in a water phantom. Software was developed to define the leaf eMLC positions of an eMLC that most closely fit each electron field shape. The effect of 1D edge feathering of the higher energy field on dose homogeneity was computed and measured for segmented-field ECT treatment plans for three 2D PTVs in a water phantom, i.e., depth from the water surface to the distal PTV surface varied as a function of the x-axis (parallel to leaf motion) and remained constant along the y-axis (perpendicular to leaf motion). Additionally, the effect of 2D edge feathering was computed and measured for one radially symmetric, 3D PTV in a water phantom, i.e., depth from the water surface to the distal PTV surface varied as a function of both axes. For the 3D PTV, the feathering scheme was evaluated for 0.1-1.0-cm leaf widths. Dose calculations were performed using the pencil beam dose algorithm in the Pinnacle{sup 3} treatment planning system. Dose verification measurements were made using a prototype eMLC (1-cm leaf width). Results: 1D discrete Gaussian edge feathering reduced the standard deviation of dose in the 2D PTVs by 34, 34, and 39%. In the 3D PTV, the broad leaf width (1 cm) of the eMLC hindered the 2D application of the feathering solution to the 3D PTV, and the standard deviation of dose increased by 10%. However, 2D discrete Gaussian edge feathering with simulated eMLC leaf widths of 0.1-0.5 cm reduced the standard deviation of dose in the 3D PTV by 33-28%, respectively. Conclusions: A five-step discrete Gaussian edge spread function applied in 2D improves the abutment dosimetry but requires an eMLC leaf resolution better than 1 cm.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2023

10. Path Modifiers

Toolpaths for additive manufacturing are typically broken into three types: closed-loop, open-loop, and space-filling. Each of these toolpaths serve a specific purpose and are often used in combination to construct an object. The standard implementation of each toolpath commands a constant motion speed and extrusion rate throughout the length of the path. The 3D printer will then make any necessary adjustments during interpolation of the motion, such as accelerating and decelerating when traversing a corner. However, sometimes more finite control and additional pathing can be needed to ensure the successful printing of a path. This chapter will outline path modifiers and how they can be used to enhance the printability of a toolpath by modifying parameters or adding segments.

Roschli, Alex↗

NASA NTRS2007

Robotic Vision-Based Localization in an Urban Environment

A system of electronic hardware and software, now undergoing development, automatically estimates the location of a robotic land vehicle in an urban environment using a somewhat imprecise map, which has been generated in advance from aerial imagery. This system does not utilize the Global Positioning System and does not include any odometry, inertial measurement units, or any other sensors except a stereoscopic pair of black-and-white digital video cameras mounted on the vehicle. Of course, the system also includes a computer running software that processes the video image data. The software consists mostly of three components corresponding to the three major image-data-processing functions: Visual Odometry This component automatically tracks point features in the imagery and computes the relative motion of the cameras between sequential image frames. This component incorporates a modified version of a visual-odometry algorithm originally published in 1989. The algorithm selects point features, performs multiresolution area-correlation computations to match the features in stereoscopic images, tracks the features through the sequence of images, and uses the tracking results to estimate the six-degree-of-freedom motion of the camera between consecutive stereoscopic pairs of images (see figure). Urban Feature Detection and Ranging Using the same data as those processed by the visual-odometry component, this component strives to determine the three-dimensional (3D) coordinates of vertical and horizontal lines that are likely to be parts of, or close to, the exterior surfaces of buildings. The basic sequence of processes performed by this component is the following: 1. An edge-detection algorithm is applied, yielding a set of linked lists of edge pixels, a horizontal-gradient image, and a vertical-gradient image. 2. Straight-line segments of edges are extracted from the linked lists generated in step 1. Any straight-line segments longer than an arbitrary threshold (e.g., 30 pixels) are assumed to belong to buildings or other artificial objects. 3. A gradient-filter algorithm is used to test straight-line segments longer than the threshold to determine whether they represent edges of natural or artificial objects. In somewhat oversimplified terms, the test is based on the assumption that the gradient of image intensity varies little along a segment that represents the edge of an artificial object.

Mchenry, Michael↗

DOE OSTI2019

Pore system characterization of organic-rich shales using nanoscale-resolution 3D imaging

The Vaca Muerta Shale in Argentina is the first signifcant commercial shale oil/gas play outside of North America. High-resolution 2D/3D imaging of shale rocks for the purpose of establishing their mineralogy, total, and connected porosities has become more and more sophisticated. In this paper, nanoscale-resolution focused ion beam (FIB)-scanning electron microscopy (SEM) nano-tomography was used to obtain images of pore structures within two organic-rich regions of interest (ROIs), selected based on correlative SEM and automated mineralogy maps. Advanced machine learning classification tools were used to segment the images and assign porosities and other components. Pore size distribution and pore connectivity analyzes revealed that about 95% of all the pores, present within the two ROIs, had a diameter of less than approximately 75 nm, and that most of these pores were poorly connected. In a similar fashion, the flow rate distribution analysis showed that pores with diameters of about 150-330 nm contributed to over 50% of the flow capacity of the connected pore systems. These results suggest that although most of the pores typically found in shales have pore diameter smaller than about 100 nm, most of the hydrocarbon production may be carried by a relatively small number of larger connected pores with pore diameter greater than about 150 nm. This research implies that a large portion of the organic-hosted pores (with diameter typically smaller than about 100 nm) may not provide permeable flow pathways for the oil and/or gas migration, and hence may have very little contribution to the hydrocarbon production.

58 GEOSCIENCES↗

DOE OSTI2021

Xylem network connectivity and embolism spread in grapevine( Vitis vinifera L.)

Abstract Xylem networks are vulnerable to the formation and spread of gas embolisms that reduce water transport. Embolisms spread through interconduit pits, but the three-dimensional (3D) complexity and scale of xylem networks means that the functional implications of intervessel connections are not well understood. Here, xylem networks of grapevine (Vitis vinifera L.) were reconstructed from 3D high-resolution X-ray micro-computed tomography (microCT) images. Xylem network performance was then modeled to simulate loss of hydraulic conductivity under increasingly negative xylem sap pressure simulating drought stress conditions. We also considered the sensitivity of xylem network performance to changes in key network parameters. We found that the mean pit area per intervessel connection was constant across 10 networks from three, 1.5-m stem segments, but short (0.5 cm) segments fail to capture complete network connectivity. Simulations showed that network organization imparted additional resistance to embolism spread beyond the air-seeding threshold of pit membranes. Xylem network vulnerability to embolism spread was most sensitive to variation in the number and location of vessels that were initially embolized and pit membrane vulnerability. Our results show that xylem network organization can increase stem resistance to embolism spread by 40% (0.66 MPa) and challenge the notion that a single embolism can spread rapidly throughout an entire xylem network.

Wason, Jay↗

DOE OSTI2015

Improving anatomical mapping of complexly deformed anatomy for external beam radiotherapy and brachytherapy dose accumulation in cervical cancer

Purpose: In the treatment of cervical cancer, large anatomical deformations, caused by, e.g., tumor shrinkage, bladder and rectum filling changes, organ sliding, and the presence of the brachytherapy (BT) applicator, prohibit the accumulation of external beam radiotherapy (EBRT) and BT dose distributions. This work proposes a structure-wise registration with vector field integration (SW+VF) to map the largely deformed anatomies between EBRT and BT, paving the way for 3D dose accumulation between EBRT and BT. Methods: T2w-MRIs acquired before EBRT and as a part of the MRI-guided BT procedure for 12 cervical cancer patients, along with the manual delineations of the bladder, cervix-uterus, and rectum-sigmoid, were used for this study. A rigid transformation was used to align the bony anatomy in the MRIs. The proposed SW+VF method starts by automatically segmenting features in the area surrounding the delineated organs. Then, each organ and feature pair is registered independently using a feature-based nonrigid registration algorithm developed in-house. Additionally, a background transformation is calculated to account for areas far from all organs and features. In order to obtain one transformation that can be used for dose accumulation, the organ-based, feature-based, and the background transformations are combined into one vector field using a weighted sum, where the contribution of each transformation can be directly controlled by its extent of influence (scope size). The optimal scope sizes for organ-based and feature-based transformations were found by an exhaustive analysis. The anatomical correctness of the mapping was independently validated by measuring the residual distances after transformation for delineated structures inside the cervix-uterus (inner anatomical correctness), and for anatomical landmarks outside the organs in the surrounding region (outer anatomical correctness). The results of the proposed method were compared with the results of the rigid transformation and nonrigid registration of all structures together (AST). Results: The rigid transformation achieved a good global alignment (mean outer anatomical correctness of 4.3 mm) but failed to align the deformed organs (mean inner anatomical correctness of 22.4 mm). Conversely, the AST registration produced a reasonable alignment for the organs (6.3 mm) but not for the surrounding region (16.9 mm). SW+VF registration achieved the best results for both regions (3.5 and 3.4 mm for the inner and outer anatomical correctness, respectively). All differences were significant (p < 0.02, Wilcoxon rank sum test). Additionally, optimization of the scope sizes determined that the method was robust for a large range of scope size values. Conclusions: The novel SW+VF method improved the mapping of large and complex deformations observed between EBRT and BT for cervical cancer patients. Future studies that quantify the mapping error in terms of dose errors are required to test the clinical applicability of dose accumulation by the SW+VF method.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Change-point Detection and Image Segmentation for Time Series of Astrophysical Images

Many astrophysical phenomena are time-varying, in the sense that their intensity, energy spectrum, and/or the spatial distribution of the emission suddenly change. This paper develops a method for modeling a time series of images. Under the assumption that the arrival times of the photons follow a Poisson process, the data are binned into 4D grids of voxels (time, energy band, and x-y coordinates), and viewed as a time series of non-homogeneous Poisson images. The method assumes that at each time point, the corresponding multiband image stack is an unknown 3D piecewise constant function including Poisson noise. It also assumes that all image stacks between any two adjacent change points (in time domain) share the same unknown piecewise constant function. The proposed method is designed to estimate the number and the locations of all of the change points (in time domain), as well as all of the unknown piecewise constant functions between any pairs of the change points. The method applies the minimum description length principle to perform this task. A practical algorithm is also developed to solve the corresponding complicated optimization problem. Simulation experiments and applications to real data sets show that the proposed method enjoys very promising empirical properties. Applications to two real data sets, the XMM observation of a flaring star and an emerging solar coronal loop, illustrate the usage of the proposed method and the scientific insight gained from it.

79 ASTRONOMY AND ASTROPHYSICS↗