loadshapedata

# Takes in the table name, directory containing a bunch of ZIP’ed shape files and a temp working directory
# Flattens all of the zip files in the shape directory into a table named “public.DATASET”
#/bin/bash
-set e 

DATASET=$1
SHAPEDIR=$2
TEMP=$3

if [ “$#” -eq 0 ];
  then
    echo “Please enter the dataset, shape directory and staging directory” 
    exit 1  
fi

if [ -z “$TEMP” ]
  then 
    echo attempting to set $TEMP directory to system temp directory
    $TEMP=$TMPDIR
fi

if [ -z “$TEMP” ]
  then 
     echo Could not set temp directory.  Bailing
     exit 1
fi

#assumes just default install of PostGIS
DB=postgres
USER_NAME=postgres
TABLE=public.$DATASET

#Clean up the staging directory, otherwise we will potentially end up with the wrong schema for each shape file
cd $TEMP 
rm $TEMP/*.shp
rm $TEMP/*.zip
rm $TEMP/*.dbf
rm $TEMP/*.prj
rm $TEMP/*.shx


cd $SHAPEDIR

#unzip files into temp working directory
for z in *.zip; do unzip -o -d $TEMP $z; done

cd $TMPDIR

#Get the first file.  Used to set the schema in PostGIS
FIRSTFILE=$(ls *.shp | head -1)

echo $FIRST FILE

#prepare the tables don't load data
shp2pgsql -p -s 4326 -W LATIN1 $FIRSTFILE $TABLE  | psql -U $USER_NAME -d $DB


#loop thru shapes and put them into SQL as 4326
 for z in *${t}.shp;
  do 
	shp2pgsql -s 4326 -W LATIN1 -a $z $TABLE | psql -d $DB -U $USER_NAME;  
  done

psql -d $DB -U $USER_NAME “CREATE INDEX geom_gix ON $TABLE USING gist\(geom\)” 

cd $TMPDIR 
rm $TMPDIR/*