Algorithms and programming tools for image processing on the MPP. Appendix A: VMS command file used to implement pp command
A detailed description of the data mapping and rotation algorithms for the Massively Parallel Processor is given.
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.
A detailed description of the data mapping and rotation algorithms for the Massively Parallel Processor is given.
Purpose: To analyze the performance of the multi-leaf collimators (MLCs) from the log files recorded during the intensity modulated radiotherapy (IMRT) treatment and to construct the relative fluence maps and do the gamma analysis to compare the planned and executed MLCs movement. Methods: We developed a program to extract and analyze the data from dynamic log files (dynalog files) generated from sliding window IMRT delivery treatments. The program extracts the planned and executed (actual or delivered) MLCs movement, calculates and compares the relative planned and executed fluences. The fluence maps were used to perform the gamma analysis (with 3% dose difference and 3 mm distance to agreement) for 3 IMR patients. We compared our gamma analysis results with those obtained from portal dose image prediction (PDIP) algorithm performed using the EPID. Results: For 3 different IMRT patient treatments, the maximum difference between the planned and the executed MCLs positions was 1.2 mm. The gamma analysis results of the planned and delivered fluences were in good agreement with the gamma analysis from portal dosimetry. The maximum difference for number of pixels passing the gamma criteria (3%/3mm) was 0.19% with respect to portal dosimetry results. Conclusion: MLC log files can be used to verify the performance of the MLCs. Patientspecific IMRT QA based on MLC movement log files gives similar results to EPID dosimetry results. This promising method for patient-specific IMRT QA is fast, does not require dose measurements in a phantom, can be done before the treatment and for every fraction, and significantly reduces the IMRT workload. The author would like to thank King Fahd University of petroleum and Minerals for the support.
Our Co-design team is from the University of Texas, working on a Department of Energy-funded project focused on the Beaumont-Port Arthur area. As part of this project, we will be developing climate-resilient design solutions for areas of the region. More on www.caee.utexas.edu.We captured aerial photos in the Port Arthur Coastal Neighborhood Community and the Golf Course on Pleasure Island, Texas, in June 2024.Aerial photos taken were through DroneDeploy autonomous flight, and models were processed through the DroneDeploy engine as well. All aerial photos are in .JPG format and contained in zipped files for each area.The processed data package includes 3D models, geospatial data, mappings, and point clouds. Please be aware that DTM, Elevation toolbox, Point cloud, and Orthomosaic use EPSG: 6588. And 3D Model uses EPSG: 3857.For using these data:- The Adobe Suite gives you great software to open .Tif files.- You can use LASUtility (Windows), ESRI ArcGIS Pro (Windows), or Blaze3D (Windows, Linux) to open a LAS file and view the data it contains.- Open an .OBJ file with a large number of free and commercial applications. Some examples include Microsoft 3D Builder, Apple Preview, Blender, and Autodesk.- You may use ArcGIS, Merkaartor, Blender (with the Google Earth Importer plug-in), Global Mapper, and Marble to open .KML files.- The .tfw world file is a text file used to georeference the GeoTIFF raster images, like the orthomosaic and the DSM. You need suitable software like ArcView to open a .TFW file.This dataset provides researchers with sufficient geometric data and the status quo of the land surface at the locations mentioned above. This dataset could streamline researchers' decision-making processes and enhance the design as well.
Purpose: To present a strategy to integrate the imaging database of a VERO unit with a treatment management system (TMS) to improve clinical workflow and consolidate image data to facilitate clinical quality control and documentation. Methods: A VERO unit is equipped with both kV and MV imaging capabilities for IGRT treatments. It has its own imaging database behind a firewall. It has been a challenge to transfer images on this unit to a TMS in a radiation therapy clinic so that registered images can be reviewed remotely with an approval or rejection record. In this study, a software system, iPump-VERO, was developed to connect VERO and a TMS in our clinic. The patient database folder on the VERO unit was mapped to a read-only folder on a file server outside VERO firewall. The application runs on a regular computer with the read access to the patient database folder. It finds the latest registered images and fuses them in one of six predefined patterns before sends them via DICOM connection to the TMS. The residual image registration errors will be overlaid on the fused image to facilitate image review. Results: The fused images of either registered kV planar images or CBCT images are fully DICOM compatible. A sentinel module is built to sense new registered images with negligible computing resources from the VERO ExacTrac imaging computer. It takes a few seconds to fuse registered images and send them to the TMS. The whole process is automated without any human intervention. Conclusion: Transferring images in DICOM connection is the easiest way to consolidate images of various sources in your TMS. Technically the attending does not have to go to the VERO treatment console to review image registration prior delivery. It is a useful tool for a busy clinic with a VERO unit.
Purpose: A system was developed that automatically calculates the organ and effective dose for individual fluoroscopically-guided procedures using a log of the clinical exposure parameters. Methods: We have previously developed a dose tracking system (DTS) to provide a real-time color-coded 3D- mapping of skin dose. This software produces a log file of all geometry and exposure parameters for every x-ray pulse during a procedure. The data in the log files is input into PCXMC, a Monte Carlo program that calculates organ and effective dose for projections and exposure parameters set by the user. We developed a MATLAB program to read data from the log files produced by the DTS and to automatically generate the definition files in the format used by PCXMC. The processing is done at the end of a procedure after all exposures are completed. Since there are thousands of exposure pulses with various parameters for fluoroscopy, DA and DSA and at various projections, the data for exposures with similar parameters is grouped prior to entry into PCXMC to reduce the number of Monte Carlo calculations that need to be performed. Results: The software developed automatically transfers data from the DTS log file to PCXMC and runs the program for each grouping of exposure pulses. When the dose from all exposure events are calculated, the doses for each organ and all effective doses are summed to obtain procedure totals. For a complicated interventional procedure, the calculations can be completed on a PC without manual intervention in less than 30 minutes depending on the level of data grouping. Conclusion: This system allows organ dose to be calculated for individual procedures for every patient without tedious calculations or data entry so that estimates of stochastic risk can be obtained in addition to the deterministic risk estimate provided by the DTS. Partial support from NIH grant R01EB002873 and Toshiba Medical Systems Corp.
This archive contains data used to support conclusions made in “Colonisation of the alpine tundra by trees: alpine neighbours assist late-seral but not early-seral conifer seedlings”, by Jabis et al., 2020. Data were collected in the alpine field location of the Alpine Treeline Warming Experiment (ATWE), on Niwot Ridge, in the Front Range of the Colorado Rocky Mountains, USA.This package includes survivorship and physiology data for limber pine (Pinus flexilis), Engelmann spruce (Picea engelmannii), and Rocky Mountain snowlover (Chionophila jamesii). Site climate data such as soil moisture and temperature are also included. This data package contains ten comma-separated-values (.csv) files, and two rich-text-format (.rtf) files all compressed within one folder named “Neighbor_data_repository.zip”. Both file types can be opened by text-edit softwares such as TextEdit (Mac) and Notepad (Windows). The files are also compatible with analyses softwares such as R. .csv files can also be opened by Microsoft Excel. Two geospatial datasets are also included in this archive: one keyhole markup language (.kml) file with four points marking the corners of the study site, and a compressed file containing two ESRI shapefiles (.shp). The .kml files can be opened with Google Earth or Google Maps, and the shapefiles can be opened using any geographic information system applications, including the entire ArcGIS suite, and QGIS. -------------------------------------------------------------------------------------------------------------------------------------------------------The elevation mountain treeline is expected to shift upward with climate warming, and seed germination and seedling survival are critical local controls on treeline expansion. Neighboring alpine plants, either through competition for resources or through altering the microclimate, can also affect seedling emergence and survival. We asked whether establishing tree seedlings and an alpine herb are similarly sensitive to alpine plant neighbours under ambient and altered climate. We imposed active heating, watering, and neighbor removal experiments for emerging conifer seedlings and an alpine herb.We compared target plant survival, photosynthetic efficiency, and water use efficiency under ambient and experimental conditions. Picea engelmannii seedlings showed lower survival compared with Pinus flexilis three weeks following neighbour removal, and after 1 year only survived in watered plots. Pinus seedlings responded to neighbour removal by lowering the quantum yield of photosynthesis (ϕPSII). Contrary to expectations from the stress gradient hypothesis, survival was reduced without neighbours near the low-elevation range limit of Chionophila jamesii.
This archive contains data used to draw conclusions in “Warming of alpine tundra enhances belowground production and shifts community towards resource acquisition traits”, by Yang et al. 2020. Data were collected on Niwot Ridge, in an alpine meadow within the Alpine Treeline Warming Experiment (ATWE) field sites in Colorado, USA. Samples were also processed in the U.S. Geological Survey Forest and Rangeland Ecosystem Science Center, in Boise, Idaho. File formats in this archive include comma-separated values (.csv), portable document format (.pdf), Microsoft Excel (.xlsx), and two types of geospatial files: keyhole markup language (.kml), and ESRI shapefiles (.shp). Leaf scans are .jpg images, and root scans are .tiff/.tif images.The .csv files can be opened using R, Microsoft Excel, or any simple text-editing software such as TextEdit and Notepad. Microsoft Excel files can be opened using Microsoft Excel, and .pdf files can be opened with Adobe Acrobat Reader, Preview, or other compatible programs. Scanned images can be opened using any photo and/or picture viewing software.The .kml file can be opened using Google Earth and Google Maps, and the shapefiles can be opened by any programs compatible with shapefiles, such as the ArcGIS Desktop suite, and QGIS.------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------Measures of belowground net primary productivity (BNPP) are required to understand whether aboveground net primary production (ANPP) changes reflect changes in allocation or are indicative of a whole plant NPP response. Plant functional traits provide a key way to scale from the individual plant to the community level, and provide insight into drivers of NPP responses to environmental change. We used infrared heaters to warm an alpine plant community at Niwot Ridge, Colorado, and applied supplemental water to compensate for soil water loss induced by warming. We measured ANPP, BNPP, and leaf and root functional traits across treatments after 5 years of continuous warming. Community-level ANPP and total NPP (ANPP + BNPP) did not respond to heating or watering, but BNPP increased in response to heating. Heating decreased community-level leaf dry matter content and increased total root length, indicating a shift in strategy from resource conservation to acquisition in response to warming.
Plants allocate a substantial portion of their fixed carbon belowground to mycorrhizal fungi in exchange for nutrients and other benefits. However, most current ecosystem models omit mycorrhizal processes, limiting our ability to predict plant–soil carbon dynamics under environmental change. To address this gap, we used a mycorrhiza-explicit soil biogeochemical model, Myco-CORPSE (Mycorrhizal Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment), to simulate tree carbon allocation to arbuscular mycorrhizal (AM) and ectomycorrhizal (ECM) fungi in temperate forests.The dataset includes outputs from two sets of model simulations:1. Perturbation experiments: Simulations across gradients of ECM dominance (0–100%), nitrogen deposition, soil temperature, and net primary productivity (NPP) to test how these factors affect mycorrhizal C allocation and nutrient cycling.2. FIA-based simulations: Model applications to over 1,800 U.S. forest sites using site-specific data from the U.S. Forest Inventory and Analysis (FIA) program, including vegetation composition, mycorrhizal type, climate, litter traits, soil properties, and N deposition.Model outputs include simulated mycorrhizal carbon allocation and related biogeochemical variables, such as soil and microbial carbon and nitrogen stocks. Data are provided in CSV format and organized by experiment type (in separate ZIP files). Python scripts for running simulations, plotting, and spatial mapping are also included and organized similarly. No proprietary software is required. These outputs support a peer-reviewed study and were used to generate figures and tables in the associated publication.
Irregularly shaped satellites, such as Phobos and Amalthea, do not lend themselves to mapping by conventional methods because mathematical projections of their surfaces fail to convey an accurate visual impression of the landforms, and because large and irregular scale changes make their features difficult to measure on maps. A digital mapping technique has therefore been developed by which maps are compiled from digital topographic and spacecraft image files. The digital file is geometrically transformed as desired for human viewing, either on video screens or on hard copy. Digital files of this kind consist of digital images superimposed on another digital file representing the three-dimensional form of a body.
The ezAlign is aimed at clustering coarse grain simulation to find common or uncommon occurrences and convert coarse (CG) grained coordinate and topology files to atomistic formats. We use a PointNet based approach to map individual frames of simulation to points in a latent space. These points are then clustered using a variety of clustering methods. Clusters are analyzed to associate them with states in the simulation. Frames can be chosen from these clusters based on proximity to cluster centers. ezAlign takes CG coordinate and topology files and converts and outputs their corresponding atomistic formats using an alignment and relaxation procedure. ezAlign is designed to convert complex, solvated biological systems including lipid membranes with drug-like molecules using GROMACS. A GROMACS checkpoint (.cpt) file is also outputted to enable continuation simulations that retain the equilibrated atomic velocities. Independent atomistic coordinates and topologies for every molecule must already be included in ezAlign/files. A number of commonly simulated biological molecules are currently provided.
The Southeast Texas Urban Integrated field lab’s Co-design team captured aerial photos in the Port Arthur Coastal Neighborhood Community and the Golf Course on Pleasure Island, Texas, in June 2024. Aerial photos taken were through autonomous flight, and models were processed through the DroneDeploy engine. All aerial photos are in .JPG format and contained in zipped files for each area. The processed data package includes 3D models, geospatial data, mappings, and point clouds. Please be aware that DTM, Elevation toolbox, Point Cloud, and Orthomosaic use EPSG: 6588. And 3D Model uses EPSG: 3857.For using these data:- The Adobe Suite gives you great software to open .Tif files.- You can use LASUtility (Windows), ESRI ArcGIS Pro (Windows), or Blaze3D (Windows, Linux) to open a LAS file and view the data it contains.- Open an .OBJ file with a large number of free and commercial applications. Some examples include Microsoft 3D Builder, Apple Preview, Blender, and Autodesk.- You may use ArcGIS, Merkaartor, Blender (with the Google Earth Importer plug-in), Global Mapper, and Marble to open .KML files.- The .tfw world file is a text file used to georeference the GeoTIFF raster images, like the orthomosaic and the DSM. You need suitable software like ArcView to open a .TFW file.This dataset provides researchers with sufficient geometric data and the status quo of the land surface at the locations mentioned above. This dataset will support researchers' decision-making processes under uncertainties.
To help guide its future data collection efforts, The DOE GTO funded a data gap analysis in FY2012 to identify high potential hydrothermal areas where critical data are needed. This analysis was updated in FY2013 and the resulting datasets are represented by this metadata. The original process was published in FY 2012 and is available here: https://pangea.stanford.edu/ERE/db/GeoConf/papers/SGW/2013/Esposito.pdf Though there are many types of data that can be used for hydrothermal exploration, five types of exploration data were targeted for this analysis. These data types were selected for their regional reconnaissance potential, and include many of the primary exploration techniques currently used by the geothermal industry. The data types include: 1. well data 2. geologic maps 3. fault maps 4. geochemistry data 5. geophysical data To determine data coverage, metadata for exploration data (including data type, data status, and coverage information) were collected and catalogued from nodes on the National Geothermal Data System (NGDS). It is the intention of this analysis that the data be updated from this source in a semi-automated fashion as new datasets are added to the NGDS nodes. In addition to this upload, an online tool was developed to allow all geothermal data providers to access this assessment and to directly add metadata themselves and view the results of the analysis via maps of data coverage in Geothermal Prospector (http://maps.nrel.gov/gt_prospector). A grid of the contiguous U.S. was created with 88,000 10-km by 10-km grid cells, and each cell was populated with the status of data availability corresponding to the five data types. Using these five data coverage maps and the USGS Resource Potential Map, sites were identified for future data collection efforts. These sites signify both that the USGS has indicated high favorability of occurrence of geothermal resources and that data gaps exist. The uploaded data are contained in two data files for each data category. The first file contains the grid and is in the SHP file format (shape file.) Each populated grid cell represents a 10k area within which data is known to exist. The second file is a CSV (comma separated value) file that contains all of the individual layers that intersected with the grid. This CSV can be joined with the map to retrieve a list of datasets that are available at any given site. The attributes in the CSV include: 1. grid_id : The id of the grid cell that the data intersects with 2. title: This represents the name of the WFS service that intersected with this grid cell 3. abstract: This represents the description of the WFS service that intersected with this grid cell 4. gap_type: This represents the category of data availability that these data fall within. As the current processing is pulling data from NGDS, this category universally represents data that are available in the NGDS and are ready for acquisition for analytic purposes. 5. proprietary_type: Whether the data are considered proprietary 6. service_type: The type of service 7. base_url: The service URL
Maps of aspen ploidy (diploid or triploid) developed using National Ecological Observatory Network (NEON) Airborne Observation Platform visible to shortwave infrared imaging spectroscopy data that were acquired over the Upper East River, Colorado in 2018. Maps of aspen distribution, as well as other land cover types that were used to screen for aspen dominated pixels are also provided. Each image file is provided as a GeoTiff with internal tiling and LZW compression. For use of ploidy maps, it is necessary to mask out non-aspen and shaded pixels. Shaded pixels are provided in an associated data package. For full details, please see associated publication and methods description. This dataset is associated with the publication in Ecological Applications.
Atom probe tomography (APT) is a material probing technique that has undergone dramatic improvements in its capability to map individual atoms within a material sample resulting in data files with hundreds of millions of atoms. Understanding the nano-structural features hidden in these massive amounts of atomic data is a crucial analysis task for materials scientists. However, fast analysis capabilities for large APT workloads remains a critical bottleneck. In this paper, we present the design, implementation and detailed performance evaluations of a parallel software capable of efficiently performing extremely time-consuming correlation analyses of massive high density APT data. Starting with shared memory implementations to motivate our design choices, we extend the implementation to hybrid architectures keeping realistic APT workloads in mind. Detailed performance analyses of three different parallel implementations of the software are supported by empirical results on a Cray XC30 and a Cray XC40 architecture. Its usefulness is demonstrated by reducing the turnaround time of an end-to-end APT correlation analysis on 100 millions atoms by three orders of magnitude using 2048 MPI ranks on 1024 nodes (24 cores per node) of a Cray XC30. The software reported here equips material scientists for the first time with a high-speed scalable capability for efficient and timely analyses of massive APT data.
Purpose: Using the graphical processing units (GPU) hardware technology, an extremely fast Monte Carlo (MC) code ARCHER{sub RT} is developed for radiation dose calculations in radiation therapy. This paper describes the detailed software development and testing for three clinical TomoTherapy® cases: the prostate, lung, and head and neck. Methods: To obtain clinically relevant dose distributions, phase space files (PSFs) created from optimized radiation therapy treatment plan fluence maps were used as the input to ARCHER{sub RT}. Patient-specific phantoms were constructed from patient CT images. Batch simulations were employed to facilitate the time-consuming task of loading large PSFs, and to improve the estimation of statistical uncertainty. Furthermore, two different Woodcock tracking algorithms were implemented and their relative performance was compared. The dose curves of an Elekta accelerator PSF incident on a homogeneous water phantom were benchmarked against DOSXYZnrc. For each of the treatment cases, dose volume histograms and isodose maps were produced from ARCHER{sub RT} and the general-purpose code, GEANT4. The gamma index analysis was performed to evaluate the similarity of voxel doses obtained from these two codes. The hardware accelerators used in this study are one NVIDIA K20 GPU, one NVIDIA K40 GPU, and six NVIDIA M2090 GPUs. In addition, to make a fairer comparison of the CPU and GPU performance, a multithreaded CPU code was developed using OpenMP and tested on an Intel E5-2620 CPU. Results: For the water phantom, the depth dose curve and dose profiles from ARCHER{sub RT} agree well with DOSXYZnrc. For clinical cases, results from ARCHER{sub RT} are compared with those from GEANT4 and good agreement is observed. Gamma index test is performed for voxels whose dose is greater than 10% of maximum dose. For 2%/2mm criteria, the passing rates for the prostate, lung case, and head and neck cases are 99.7%, 98.5%, and 97.2%, respectively. Due to specific architecture of GPU, modified Woodcock tracking algorithm performed inferior to the original one. ARCHER{sub RT} achieves a fast speed for PSF-based dose calculations. With a single M2090 card, the simulations cost about 60, 50, 80 s for three cases, respectively, with the 1% statistical error in the PTV. Using the latest K40 card, the simulations are 1.7–1.8 times faster. More impressively, six M2090 cards could finish the simulations in 8.9–13.4 s. For comparison, the same simulations on Intel E5-2620 (12 hyperthreading) cost about 500–800 s. Conclusions: ARCHER{sub RT} was developed successfully to perform fast and accurate MC dose calculation for radiotherapy using PSFs and patient CT phantoms.
This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (Chem. Mater. 2024, 36, 11306-11316). The calculations characterize MnRhBi3, an orthorhombic (Cmmm) van der Waals-layered intermetallic compound that cleaves easily between neighboring Bi layers. The dataset is organized into three calculation types: (i) Bulk: Structural relaxations of the periodic MnRhBi3 crystal in antiferromagnetic (AFM) and ferromagnetic (FM) configurations, using the vdW-DF-optB86b functional. These provide the equilibrium lattice constants, magnetic energy differences (AFM is 0.5 meV/f.u. lower than FM), and magnetic moments (4.4 µB/Mn, 0.17 µB/Rh, 0.18 µB/Bi) reported in Table 1 of the main text. (ii) Slab: Same magnetic configurations computed with an 18 Ang vacuum layer introduced between Bi layers, used to calculate the cleavage energy Ec = 0.56 J/m2 (AFM) and 0.57 J/m2 (FM), establishing MnRhBi3 as a van der Waals-layered material comparable to graphite, MoS2, and CrI3. (iii) ELF: Single-point calculation on the relaxed bulk AFM geometry with LELF=.TRUE., producing the ELFCAR file used to generate electron localization function isosurfaces and contour maps (Fig. 2, main text) showing Bi lone pairs directed into the van der Waals gaps. All folders contain CONTCAR, INCAR, KPOINTS, OUTCAR, and POSCAR. The ELF/ folder additionally contains ELFCAR. Calculations were performed using VASP 6.3.2 with PBE + vdW-DF-optB86b, PAW potentials, and an energy cutoff of 800 eV.
The data measurements of CH4 and CO2 at different patch types (classified by dominant macrophyte vegetation or hydrological state). CH4 and CO2 data include 3 types of measurements: 1) Surface flux measurements from chamber measurements. 2) Leaf flux measurements from plant chambers, providing the rate per leaf area of methane flow through three species of macropyte aerenchemous plants (Water Lily, American Lotus, and Cattail) with the corresponding flux of CO2 at the same leaves. 3) Pore water concentrations of dissolved CO2 and CH4 concentrations. These are measured along vertical profiles at varying depths from the soil-water interface using in-situ porewater dialysis profiles (peepers). Observations of 1 and 3 are available monthly for most months of the growing season (May-September) of 2015, 2016 and 2018. Leaf flux measurements (plant chambers) are available only in June, July and August of 2018. We also included the patch type classification maps during 2015, 2017 and 2018 the relative water elevation map (bathymetry) for the site.In 8/23/2020 we added the data file ChamberFluxes_2018_Separated_Diffusive_Ebullitive_Fluxes.xlsx with observation related to chamber measurements in OWC during September 2018. In these observations, we partitioned the methane flux to its diffusive and ebullitive components, using an analysis approach we describe in section 5 of the methods "Ebullition and diffusion partitioning".
Procedures and software for the digital mosaicking of Landsat data are discussed. It is noted that image processing in support of JPL's planetary program has furnished the software and procedures necessary to achieve digital image mosaicking of vidicon imagery using rigid projective geometry. The incorporation of ground-control points, by either manual or automatic ground-control point file identification, has yielded rms positional accuracies exceeding the National Map Accuracy Standards. Rotation to north vertical is effected at low computational cost. It is possible to cut input data frames in any arbitrary shape to remove cloud cover and accommodate terrain offset effects. In similar fashion, the final digital mosaic can be segmented arbitrarily to conform to user requirements. Information is presented on applications in Pennsylvania and Bolivia.