| 1 |
11496
|
psarando
|
#!/bin/bash
|
| 2 |
|
|
|
| 3 |
|
|
# Bash script to create a new postgis database and prep it with GADM
|
| 4 |
|
|
# data for geovalidation purposes.
|
| 5 |
|
|
#
|
| 6 |
|
|
# Basic workflow:
|
| 7 |
|
|
# 1. create polygon geometry upon import (loader does this)
|
| 8 |
|
|
# 2. create simplified polygon geometry after import
|
| 9 |
|
|
# 3. create polygon geography after import
|
| 10 |
|
|
# 4. create indexes
|
| 11 |
|
|
# 5. cluster data in geom index order
|
| 12 |
|
|
#
|
| 13 |
|
|
# todo:
|
| 14 |
|
|
# * better handle Antarctica problem so it can be included instead of
|
| 15 |
|
|
# just dropped? probably not a huge deal for plants...
|
| 16 |
|
|
#
|
| 17 |
|
|
# Jim Regetz
|
| 18 |
|
|
# NCEAS
|
| 19 |
|
|
# Created November 2012
|
| 20 |
|
|
#
|
| 21 |
|
|
# refactored and reorganized by
|
| 22 |
|
|
# Paul Sarando
|
| 23 |
|
|
# iPlant Collaborative
|
| 24 |
|
|
# Updated Oct 2013
|
| 25 |
|
|
|
| 26 |
|
|
DB_NAME="geoscrub"
|
| 27 |
|
|
DB_USER="bien"
|
| 28 |
|
|
DB_HOST="localhost"
|
| 29 |
|
|
SCRIPT_DIR="$(dirname $0)"
|
| 30 |
|
|
|
| 31 |
|
|
# GADM data originally available at http://www.gadm.org/data2/gadm_v2_shp.zip
|
| 32 |
|
|
# gadm.org now links this file from biogeo.ucdavis.edu.
|
| 33 |
|
|
GADM_DATA_URL="http://biogeo.ucdavis.edu/data/gadm2/gadm_v2_shp.zip"
|
| 34 |
|
|
GADM_DATA_DIR="${SCRIPT_DIR}/gadm_v2_shp"
|
| 35 |
|
|
|
| 36 |
|
|
#
|
| 37 |
|
|
# assemble input data
|
| 38 |
|
|
#
|
| 39 |
|
|
|
| 40 |
|
|
echo "Updating gadm2 tables from GADM data in ${GADM_DATA_DIR}"
|
| 41 |
|
|
echo -n "Note, to force data to download again from ${GADM_DATA_URL},"
|
| 42 |
|
|
echo " delete the directory ${GADM_DATA_DIR} before running this script."
|
| 43 |
|
|
|
| 44 |
|
|
# Check for GADM data (320MB zip file) unziped in data directory.
|
| 45 |
|
|
if [[ ! -d "$GADM_DATA_DIR" ]]; then
|
| 46 |
|
|
echo "making directory ${GADM_DATA_DIR}"
|
| 47 |
|
|
mkdir -p "$GADM_DATA_DIR"
|
| 48 |
|
|
|
| 49 |
|
|
if [[ $? != 0 ]]; then
|
| 50 |
|
|
echo "Could not create directory ${GADM_DATA_DIR}"
|
| 51 |
|
|
exit 1
|
| 52 |
|
|
fi
|
| 53 |
|
|
|
| 54 |
|
|
pushd "$GADM_DATA_DIR"
|
| 55 |
|
|
wget -O gadm_v2_shp.zip "$GADM_DATA_URL"
|
| 56 |
|
|
unzip gadm_v2_shp.zip
|
| 57 |
|
|
rm -f gadm_v2_shp.zip
|
| 58 |
|
|
popd
|
| 59 |
|
|
fi
|
| 60 |
|
|
|
| 61 |
|
|
#
|
| 62 |
|
|
# create and populate gadm2 table
|
| 63 |
|
|
#
|
| 64 |
|
|
|
| 65 |
|
|
psql -e -U "$DB_USER" -h "$DB_HOST" -d "$DB_NAME" -c "DROP TABLE IF EXISTS gadm2"
|
| 66 |
|
|
if [[ $? != 0 ]]; then
|
| 67 |
|
|
echo "Could not drop GADM2 table in ${DB_NAME} database."
|
| 68 |
|
|
exit 1
|
| 69 |
|
|
fi
|
| 70 |
|
|
|
| 71 |
|
|
echo "Creating gadm2 table with shp2pgsql from ${GADM_DATA_DIR} data."
|
| 72 |
|
|
|
| 73 |
|
|
pushd "$GADM_DATA_DIR"
|
| 74 |
|
|
|
| 75 |
|
|
# load gadm2 data (took 4.7 minutes on willow, 26-Oct-2012)
|
| 76 |
|
|
shp2pgsql -s 4326 -W latin1 gadm2 gadm2 "$DB_NAME" | \
|
| 77 |
|
|
psql --set ON_ERROR_STOP=1 -U "$DB_USER" -h "$DB_HOST" -d "$DB_NAME"
|
| 78 |
|
|
if [[ $? != 0 ]]; then
|
| 79 |
|
|
popd
|
| 80 |
|
|
echo "Could not load GADM2 data into ${DB_NAME} database."
|
| 81 |
|
|
exit 1
|
| 82 |
|
|
fi
|
| 83 |
|
|
|
| 84 |
|
|
popd
|
| 85 |
|
|
|
| 86 |
|
|
# create indexes and additional columns
|
| 87 |
|
|
psql -e -U "$DB_USER" -h "$DB_HOST" -d "$DB_NAME" --set ON_ERROR_STOP=1 < "${SCRIPT_DIR}"/update.gadm2.sql
|
| 88 |
|
|
if [[ $? != 0 ]]; then
|
| 89 |
|
|
echo "Could not update GADM2 data and indexes in ${DB_NAME} database."
|
| 90 |
|
|
exit 1
|
| 91 |
|
|
fi
|
| 92 |
|
|
|
| 93 |
|
|
echo "Update gadm2 tables successfully completed."
|
