38 lines
1.3 KiB
Bash
38 lines
1.3 KiB
Bash
#!/usr/bin/env bash
|
|
|
|
# Assign latitute and longitute to buildings with OpenTOID data
|
|
# - assume postgres connection details are set in the environment using PGUSER, PGHOST etc.
|
|
|
|
: ${1?"Usage: $0 ./path/to/opentoid/dir"}
|
|
|
|
opentoid_dir=$1
|
|
|
|
echo "Create a temporary table for coordinates..."
|
|
psql -c "CREATE TABLE open_toid (
|
|
toid varchar,
|
|
version_number smallint,
|
|
version_date date,
|
|
source_product varchar,
|
|
easting float,
|
|
northing float,
|
|
latitute float,
|
|
longitude float
|
|
);"
|
|
|
|
echo "Loading Open TOID CSV(s) to temporary table..."
|
|
find $opentoid_dir -type f -name '*.csv' \
|
|
-printf "$opentoid_dir/%f\n" | \
|
|
parallel \
|
|
cat {} '|' psql -c "\"COPY open_toid ( toid, version_number, version_date, source_product, easting, northing ) FROM stdin WITH CSV HEADER;\""
|
|
|
|
# Convert the northing/easting coordinates to latitude/longitute with PostGIS
|
|
psql -c "select AddGeometryColumn( 'public', 'open_toid', 'the_geom', 27700, 'POINT', 2);"
|
|
psql -c "update open_toid set the_geom=GeomFromText('POINT('||easting||' '||northing||')',27700);"
|
|
psql -c "update open_toid set longitude=st_x(st_transform(the_geom,4326)), latitude=st_y(st_transform(the_geom,4326));"
|
|
|
|
# Update the buildings table with coordinates
|
|
# psql -c ""
|
|
|
|
|
|
# Delete the temporary table
|
|
# psql -c "DROP TABLE open_toid;" |