155 lines
5.8 KiB
Bash
155 lines
5.8 KiB
Bash
#!/bin/sh
|
|
##
|
|
# This software was developed and / or modified by Raytheon Company,
|
|
# pursuant to Contract DG133W-05-CQ-1067 with the US Government.
|
|
#
|
|
# U.S. EXPORT CONTROLLED TECHNICAL DATA
|
|
# This software product contains export-restricted data whose
|
|
# export/transfer/disclosure is restricted by U.S. law. Dissemination
|
|
# to non-U.S. persons whether in the United States or abroad requires
|
|
# an export license or other authorization.
|
|
#
|
|
# Contractor Name: Raytheon Company
|
|
# Contractor Address: 6825 Pine Street, Suite 340
|
|
# Mail Stop B8
|
|
# Omaha, NE 68106
|
|
# 402.291.0100
|
|
#
|
|
# See the AWIPS II Master Rights File ("Master Rights File.pdf") for
|
|
# further licensing information.
|
|
#
|
|
#
|
|
# SOFTWARE HISTORY
|
|
# Date Ticket# Engineer Description
|
|
# ------------ ---------- ----------- --------------------------
|
|
# 03/25/2014 #2664 randerso Added support for importing non-WGS84 shape files
|
|
# 10/23/2014 #3685 randerso Fixed bug where .prj was not recognized when shape file
|
|
# was in the current directory (no directory specified)
|
|
# 02/11/2016 #5348 randerso Add code to create a county_names view into the county table
|
|
#
|
|
##
|
|
|
|
function usage()
|
|
{
|
|
echo
|
|
echo usage: `basename $0` shapefile schema table [simplev [dbUser [dbPort [installDir [srid]]]]]
|
|
echo "where: shapefile - pathname of the shape file to be imported"
|
|
echo " schema - database schema where the shape file is to be imported"
|
|
echo " table - database table where the shape file is to be imported"
|
|
echo " simplev - optional list of geometry simplification levels to be created"
|
|
echo " dbUser - optional database user id"
|
|
echo " dbPort - optional database port number"
|
|
echo " installDir- optional directory path to awips installation"
|
|
echo " srid - optional srid of source shape file"
|
|
echo "example: `basename $0` uscounties.shp mapdata County 0.064,0.016,0.004,0.001 awips 5432 /awips2"
|
|
}
|
|
|
|
if [ $# -lt 3 ] ; then
|
|
usage
|
|
exit -1
|
|
fi
|
|
|
|
SHP_PATH=`readlink -f ${1}`
|
|
SHP_DIR="${SHP_PATH%/*}" # shape file dir
|
|
SHP_NAME="${SHP_PATH##*/}" # shape file name with extension
|
|
SHP_BASE="${SHP_NAME%.*}" # shape file name without extension
|
|
SHP_EXT="${SHP_NAME##*.}" # shape file extension
|
|
PRJ_PATH="${SHP_DIR}/${SHP_BASE}.prj"
|
|
|
|
SCHEMA=`echo "${2}" | tr '[:upper:]' '[:lower:]'`
|
|
TABLE=`echo "${3}" | tr '[:upper:]' '[:lower:]'`
|
|
SIMPLEVS=${4}
|
|
|
|
if [ -z $5 ] ; then
|
|
PGUSER=awips
|
|
else
|
|
PGUSER=${5}
|
|
fi
|
|
|
|
if [ -z $6 ] ; then
|
|
PGPORT=5432
|
|
else
|
|
PGPORT=${6}
|
|
fi
|
|
|
|
if [ -z $7 ] ; then
|
|
PGBINDIR=/awips2/postgresql/bin/
|
|
PSQLBINDIR=/awips2/psql/bin/
|
|
else
|
|
PGBINDIR=${7}/postgresql/bin/
|
|
PSQLBINDIR=${7}/psql/bin/
|
|
fi
|
|
|
|
if [ -z $8 ] ; then
|
|
if [ -e $PRJ_PATH ]
|
|
then
|
|
echo
|
|
echo "WARNING, found projection file: ${PRJ_PATH}"
|
|
echo "It is probable that this shape file is not in EPSG:4326 (WGS 84, unprojected lat/lon) format."
|
|
echo "Please determine the correct srid by uploading the .prj file to http://prj2epsg.org/search"
|
|
echo "and re-run `basename $0` supplying the correct srid."
|
|
usage
|
|
exit -1
|
|
fi
|
|
|
|
SRID=4326
|
|
else
|
|
SRID=${8}:4326
|
|
fi
|
|
|
|
if [ ! -r ${SHP_PATH} ]; then
|
|
echo
|
|
echo "ERROR, file not found or is not readable: ${SHP_PATH}"
|
|
exit -1
|
|
fi
|
|
|
|
echo " Importing ${SHP_NAME} into ${SCHEMA}.${TABLE} ..."
|
|
|
|
#
|
|
# If updating county table drop the county names view
|
|
#
|
|
if [ "county" == ${TABLE} ] ; then
|
|
${PSQLBINDIR}psql -d maps -U ${PGUSER} -q -p ${PGPORT} -c "
|
|
DROP VIEW IF EXISTS ${SCHEMA}.county_names;
|
|
"
|
|
fi
|
|
|
|
${PSQLBINDIR}psql -d maps -U ${PGUSER} -q -p ${PGPORT} -c "
|
|
DELETE FROM public.geometry_columns WHERE f_table_schema = '${SCHEMA}' AND f_table_name = '${TABLE}';
|
|
DELETE FROM ${SCHEMA}.map_version WHERE table_name='${TABLE}';
|
|
DROP TABLE IF EXISTS ${SCHEMA}.${TABLE}
|
|
"
|
|
${PGBINDIR}shp2pgsql -W LATIN1 -s ${SRID} -g the_geom -I ${SHP_PATH} ${SCHEMA}.${TABLE} | ${PSQLBINDIR}psql -d maps -U ${PGUSER} -q -p ${PGPORT} -f -
|
|
${PSQLBINDIR}psql -d maps -U ${PGUSER} -q -p ${PGPORT} -c "
|
|
INSERT INTO ${SCHEMA}.map_version (table_name, filename) values ('${TABLE}','${SHP_NAME}');
|
|
SELECT AddGeometryColumn('${SCHEMA}','${TABLE}','the_geom_0','4326',(SELECT type FROM public.geometry_columns WHERE f_table_schema='${SCHEMA}' and f_table_name='${TABLE}' and f_geometry_column='the_geom'),2);
|
|
UPDATE ${SCHEMA}.${TABLE} SET the_geom_0=ST_Segmentize(the_geom,0.1);
|
|
CREATE INDEX ${TABLE}_the_geom_0_gist ON ${SCHEMA}.${TABLE} USING gist(the_geom_0);
|
|
"
|
|
|
|
#
|
|
# If updating county table recreate the county names view
|
|
#
|
|
if [ "county" == ${TABLE} ] ; then
|
|
${PSQLBINDIR}psql -d maps -U ${PGUSER} -q -p ${PGPORT} -c "
|
|
CREATE OR REPLACE VIEW ${SCHEMA}.county_names AS
|
|
SELECT countyname as name, ST_SetSRID(ST_Point(lon,lat), 4326)::geometry(Point, 4326) as the_geom
|
|
FROM ${SCHEMA}.county;
|
|
"
|
|
fi
|
|
|
|
if [ -n "$SIMPLEVS" ] ; then
|
|
echo " Creating simplification levels ${SIMPLEVS}..."
|
|
IFS=", "
|
|
for LEV in $SIMPLEVS ; do
|
|
echo " Creating simplified geometry level $LEV ..."
|
|
IFS="."
|
|
SUFFIX=
|
|
for x in $LEV ; do SUFFIX=${SUFFIX}_${x} ; done
|
|
${PSQLBINDIR}psql -d maps -U ${PGUSER} -q -p ${PGPORT} -c "
|
|
SELECT AddGeometryColumn('${SCHEMA}','${TABLE}','the_geom${SUFFIX}','4326',(SELECT type FROM public.geometry_columns WHERE f_table_schema='${SCHEMA}' and f_table_name='${TABLE}' and f_geometry_column='the_geom'),2);
|
|
UPDATE ${SCHEMA}.${TABLE} SET the_geom${SUFFIX}=ST_Segmentize(ST_Multi(ST_SimplifyPreserveTopology(the_geom,${LEV})),0.1);
|
|
CREATE INDEX ${TABLE}_the_geom${SUFFIX}_gist ON ${SCHEMA}.${TABLE} USING gist(the_geom${SUFFIX});"
|
|
done
|
|
fi
|
|
${PGBINDIR}vacuumdb -d maps -t ${SCHEMA}.${TABLE} -U ${PGUSER} -p ${PGPORT} -vfz
|