Readme
Author: Alex Chansky, some code and methods sourced from Michelle Meyer
Date: 3/26/21

ClippingCTIs.ipynb: Extract CTI, define with WGS 1984 (4326), then project

DroppingNonLiqCells.ipynb: For 1 km, 100x100 m buffer polygon, identifies cells containing 1 or more liq point or poly, then drops all other cells.

DroppingOffshoreNonLiqPoints.ipynb: For all events, some points fall offshore, where DC = 0. These are deleted as we want to keep only onshore points for non-liquefaction sampling. As fishnet polygon is irrelevant for non-liq points, we can ignore polygon altogether.

ExportingData.ipynb: Exporting data from liq 1 km 100x100 m buffer shapefile of liq cells to CSV file.

ExtractLiqPtsData&DropBadData.ipynb: Extracting liquefaction data at centroid of liquefaction cell, or the closest cell where valid data exists for that layer. See full documentation for further details.

LiqEventsFishnets.ipynb: Setting up fishnets for 1km buffers (100x100 m) and space between 1km and 15km buffers (1500x1500 m).

LiqEventsNonLiqCells01.ipynb: Erase 1km buff from 15km buff, then lay fishnet across space between them of 1.5x1.5 km. Results in a poly shapefile and a point shapefile (one point at centroid of each polygon).

LiqEventsNonLiqCellsBufferExtent02.ipynb: Determining X and Y mins and Maxes of buffer extents and place in CSV, which tells us where to place fishnets 

MosaicASTER_03.ipynb: Mosaic together cells representing “attributes” of water bodies, where 0=land, 
    1=ocean, 2=rivers, and 3=lakes.

MosaicClipping02.ipynb: Create new folder of quake name, extract DEM to PGA's extent and project. 
    Then, go into Downloads directly and use GDAL codes to create slope, TPI, TRI
    Then Use the last block of this code to project slope, TPI, and TRI into appropriate folders

MosaicCTIs.ipynb: Creating mosaic of all CTI layers

MosaicDEMs.ipynb: Mosaic together elevetion (aka DEM or GMTED) and associated standard deviation tiles

ReclassEucDist.ipynb: Reclassify and Euclidean distance tools to obtain DC, DR, DWB layers from single attribute layer

SamplingNoLiqPts_NoLiqEvents.ipynb: For non-liquefaction events, this script samples non-liquefaction points using a radius of 40 km and 1.5 km between sampled points. Points over the ocean (DC=0) are later dropped.

WaterDEMEucAllocation.ipynb: Allocate elevation of nearest water body value to every cell in the region


Typical Order of Actions:
- Obtain / Digitize points and polygons representing liquefaction as GIS shapefiles
- Create complementary spreadsheet containing each event listing projection, 1 or 0 for "Liq", and 1 or 0 for "Working"
- Download all appropriate variable data for CTI, DEM, Water Bodies (see full documentation for further details)
- Run all scripts:

Putting rasters together and clipping:
1. MosaicCTIs.ipynb
2. ClippingCTIs.ipynb
3. MosaicDEMs.ipynb
4. MosaicClipping02
5. MosaicASTER_03.ipynb
6. ReclassEucDist.ipynb
7. WaterDEMEucAllocation.ipynb

Establishing buffers and fishnets:
8. LiqEventsNonLiqCellsBufferExtent02.ipynb
9. SamplingNoLiqPts_NoLiqEvents.ipynb
10. LiqEventsNonLiqCells01.ipynb
11. LiqEventsFishnets.ipynb

Eliminating unwanted data:
12. DroppingNonLiqCells.ipynb
13. DroppingOffshoreNonLiqPoints.ipynb

Extracting and exporting values from raster data:
14. ExtractLiqPtsData&DropBadData.ipynb
15. ExportingData.ipynb