Contained within are copies of the variety of Python- and R-based scripts that were used as a part of the New Mexico statewide landslide susceptibility project. With the exception of Mauro Rossi and Paolo Reichenbach's LAND-SE script, none of these scripts was written with the intention of being used for any other project. They are provided as a part of the deliverable for completeness and transparency. They are provided with no support or guaruntee of operability. They are not intended to be used for any additional assessments of landslide or other geohazard susceptibility. Neither the authors of these scripts nor the agencies that employ them are responsible for the results of any subsequent uses of these scripts. Scripts are presented as-is. Documentation varies, and is mostly contained in comments within each script.


------------------------

The LAND-SE script was written by Mauro Rossi and Paola Reichenbach of the geomorphology group at the Istituto di Ricerca per la Protezione Idrogeologica, Consiglio Nazionale delle Ricerche (Research Institute for Geo-Hydrological Protection of the Italian National Research Council). The script is available at https://github.com/maurorossi/LAND-SE. Documentation is provided in the LAND-SE folder. A description of the script and its functionality was published as:

Rossi, M., and Reichenbach, P., 2016, LAND-SE: a software for statistically based landslide susceptibility zonation, version 1.0: Geoscience Model Development 9, p. 3533-3543.

Modifications to the original LAND-SE script are by Colin Cikoski, New Mexico Bureau of Geology and Mineral Resources. These modifications do not change the core functions of the original script, but remove some aspects of the LAND-SE script that were not used as a part of this particular project.


-----------------------

Original scripts written for this project include:

-ContinuousDataBinningScript.py:
Processes rasters of continuous data for restricted equal-frequency binning

-LANDSE_TrainingData_InputGenerator.py:
Takes as inputs landslide extents, training area extents, geology feature classes, and raster indicator variables and generates shapefiles and text files in the approriate format and structure for use in the LAND-SE script

-LANDSE_ExternalValidation_DataSampler.py:
Takes as inputs landslide extents, external validation area extents, geology feature classes, raster indicator variables, and a set of R-determined model coefficients to generate an external validation shapefile and text table for use in a modified version of the LAND-SE script that only performs the validation.

-SpatialSubsampleScript.py:
Takes the same inputs as above, and creates randomly-located subsamples of the combined training-validation area, samples these data as per the training data sampling method, and exports text data tables for generating subsample models for assessing model coefficient variability within the training-validation area extent.

-SpatialSubsampleStatsCompileScript.R:
Takes as inputs the text tables generated by the above script and processes the data tables for assessing model coefficient variability.

-ModelProbabilityMapGenerator.py:
Takes as inputs rasters of indicator variables and model coefficients to create raster model probability maps.


These scripts were written by Colin Cikoski, New Mexico Bureau of Geology and Mineral Resources, Colin.Cikoski@nmt.edu, for the purposes of this project and this project alone. Additional documentation is available as comments within each script file.