Meeting radiation dosimetry capacity requirements of population-scale exposures by geostatistical sampling

This is the data repository for the PLOS ONE Manuscript: "Meeting radiation dosimetry capacity requirements of population-scale exposures by geostatistical sampling". This repository contains the following data:

1. "State-and-Subdivision-Boundary-Files.Edited-for-ArcMap-10.4.KML-Format.zip":

This file contains modified U.S. state and sub-division boundary files [in KML format], which can be imported into ArcMap using its KMLtoLayer function. These files have been modified to prevent sub-division naming issues that we encountered when importing boundary data into ArcMap: A) State sub-divisions with identical names are considered a single sub-division by ArcMap (corrected by adding a letter after each sub-division of the same name, i.e. CenterA, CenterB, etc), and; B) ArcMap would only identify the sub-division by its first word if sub-division name contained spaces (corrected by converting all spaces into dashes). 

2. "HPAC-Plumes.Processed.zip" and "HPAC-Plumes.Unprocessed.zip":

These files contain HPAC plume coordinate (WGS1984) and dose (in cGy) values for all scenarios discussed in the manuscript. We provide "processed" and "unprocessed" HPAC plume data files. The "unprocessed" HPAC plume data is provided in its original XML format, which cannot be imported into ArcMap directly. The "processed" HPAC plumes are provided in tab-delimited X,Y,Z format (Latitude, Longitude, and Dose). We have also added a "0 cGy" contour in the "processed" plumes (surrounding the HPAC plume), as the presence of unirradiated data points adjacent to the plume was found to be crucial for accurate kriging, since these points served as boundaries for kriging.

3. "Final-Derived-Plumes.Data-Points.zip":

This file contains geostatistically-derived plume coordinate (WGS1984) and dose (in cGy) values for all scenarios discussed in the manuscript. Data is in comma-delimited format (Latitude, Longitude, and Dose). Data points consist of a set of initial coordinates generated at random locations within each Census sub-division using the ArcMap tool, ‘CreateRandomPoints_management’, and subsequent points generated by densification (the geostatistical procedure that targets and localizes an additional small cohort of irradiated individuals to mitigate uncertainty in environmental measurements). These data points were assigned radiation level values corresponding to the adjacent outer HPAC contour by a script comparing each sample with its location within the HPAC plume of the same scenario.

4. “Intermediate-Derived-Plumes.Data-Points.zip”

This archive contains coordinate data (WGS1984) and dose values (in cGy) for all intermediary steps of plume development (using our geostatistical method) for all scenarios. Like (3), the data is comma-delimited (Latitude, Longitude, and Dose), and were assigned radiation level values by a script comparing sampling locations with the location of the HPAC plume of the same scenario. Scenario replicate folders contains text files for each iteration step of the plume derivation process, including a file containing just the initial random sampling (“Iteration-1”), a file containing initial sampling and sampling locations selected by the first densification step (“Iteration-2”), a file containing initial sampling and sampling locations selected by the first and second densification steps (“Iteration-3”), and so on.

This archive also contains a Table (“Progression-of-New-Densification-Selected-Sampling-Locations-For-All-Scenarios.xslx”) which provides a categorical breakdown of how many unique densification-selected sampling locations occur within the irradiated region (i.e. overlap the HPAC plume) for each iteration of all scenario replicates. The fraction of irradiated to unirradiated sampling locations varies among each scenario and individual replicates for the same scenario. Our analysis shows that these results depend on the population densities and exact topography of the HPAC plume which is different among each scenario.

5. "Geostatistical-Sampling-Project.All-Scripts.zip"

This archive contains all programs required for this project. This includes Python scripts meant to be run within the ArcMap software environment (for random point generation and data extraction), and Perl scripts used to process HPAC and U.S. State and Sub-division boundary files, and to assign radiation values to sample locations based on a modified HPAC plume. A java program, “CompareReplicates.jar”, compares the overlapping areas between a pair of polygons that overlap one other using the ArcMap software environment, and requires access to the ArcGIS Runtime SDK (https://developers.arcgis.com/arcgis-runtime/).