X-Git-Url: http://gitweb.michael.orlitzky.com/?p=dead%2Fcensus-tools.git;a=blobdiff_plain;f=makefile;h=7d765c6cdaab32bf2a280020097a7455e5e91183;hp=36543d4e156bb7dd5a6ea8c60e948d613dd29739;hb=f38573f5ab9115c2b4f315cfb7cc2be14e025f54;hpb=d6aa9f176f67ee918f854c1df3ec5b072a502288 diff --git a/makefile b/makefile index 36543d4..7d765c6 100644 --- a/makefile +++ b/makefile @@ -1,89 +1,163 @@ -DB_NAME='census' -DB_USER='postgres' -TIGER_DATA_URL='http://www2.census.gov/geo/tiger/TIGER2009/24_MARYLAND/tl_2009_24_tabblock00.zip' -TIGER_SRID='4269' - -# Starting with PostGIS 1.4.0, these paths are calculated at install -# time using the pg_config utility. Rather than try to guess where -# PostGIS will wind up installed, we can just check the output of -# pg_config ourselves. -PG_BINDIR=`pg_config --bindir` -PG_SHAREDIR=`pg_config --sharedir` +DB_NAME=census +DB_USER=postgres +TIGER_SRID=4269 +SHAPELY_URL=http://pypi.python.org/packages/source/S/Shapely/Shapely-1.0.14.tar.gz + + +# Dark magic. We set these makefile variables to be the result of the +# 'shell' function. The shell function, in turn, executes a Python +# script which determines the locations of these files. +SHP2PGSQL := $(shell bin/find_file_paths --root /usr --single shp2pgsql) +POSTGIS_SQL := $(shell bin/find_file_paths --root /usr lwpostgis.sql postgis.sql) +SPATIAL_REF_SYS_SQL := $(shell bin/find_file_paths --root /usr spatial_ref_sys.sql) # Necessary to run test/data without prerequisites. # -.PHONY : test data +.PHONY : test data lib # The default task, since it comes first in the list. # -all: clean test +all: clean lib test test: ./bin/run_tests +# Download or check out any third-party libraries. +lib: + if [ ! -d lib/Shapely ]; then \ + wget -O shapely.tar.gz $(SHAPELY_URL); \ + tar -xvzf shapely.tar.gz -C lib/ ; \ + rm shapely.tar.gz; \ + mv lib/Shapely* lib/Shapely; \ + fi; + + # Remove byte-compiled python code. # clean: find ./ -name '*.pyc' -print0 | xargs -0 rm -f -# Download the shapefiles from Tiger if they don't already exist. -# For now, we're only dealing with the Census 2000 Maryland Block -# data, so the filenames are hard-coded. Easy enough to change. -# data: - mkdir -p data/census2000/maryland/block - if [ ! -f data/census2000/maryland/block/tl_2009_24_tabblock00.shp ]; then \ - wget -O tmp.zip $(TIGER_DATA_URL); \ - unzip tmp.zip -d ./data/census2000/maryland/block; \ - rm tmp.zip; \ - fi; + bin/download_data + + +# This imports the Tiger data using shp2pgsql. The shapefiles +# should exist, since this task depends on the "data" task, which +# downloads said shapefiles. +# +# After the TIGER import is done, we use the sf1blocks2sql script to +# parse and import the geographic header record information. +# +db: data newdb tiger_blocks_table tiger_lines_table sf1_blocks_table +# All Blocks +# +# The table already exists, so we can append to it, and we don't have +# to create the GiST index. + for state in data/census2000/*; do \ + $(SHP2PGSQL) \ + -a \ + -s $(TIGER_SRID) \ + -D \ + $$state/blocks/*.shp \ + tiger_blocks \ + | psql -U $(DB_USER) -d $(DB_NAME); \ + done; + +# All Lines +# +# Since the table and index already exist, we can utilize -a, +# and leave -I out. + for state in data/census2000/*; do \ + for shapefile in $$state/lines/*.shp; do \ + echo "Importing $$shapefile."; \ + $(SHP2PGSQL) \ + -a \ + -s $(TIGER_SRID) \ + $$shapefile \ + tiger_lines \ + | bin/filter-transactions \ + | psql -U $(DB_USER) -d $(DB_NAME) \ + > /dev/null; \ + done; \ + done; + +# Summary File 1 +# +# Run all of the geo (uf1) files through the import script. This has +# to happen after the blocks import since we impose a foreign key +# restriction. + for state in data/census2000/*; do \ + bin/sf1blocks2sql $$state/sf1/*.uf1 sf1_blocks \ + | psql -U postgres -d $(DB_NAME) \ + > /dev/null; \ + done; + -# This task does a couple of things. First, it drops and re-creates -# the DB_NAME database (or schema, whatever). Then, it adds PL/pgSQL -# support to the database. +# First, we drop and re-create the DB_NAME database (or schema, +# whatever). Then, we add PL/pgSQL support to the database. # # At that point, we import the two PostGIS files, postgis.sql and # spatial_ref_sys.sql. The postgis.sql file contains the geometry # functions, while spatial_ref_sys.sql contains a table of SRIDs, and # their associated properties. PostGIS requires both. # -# Then, we import the Tiger data using shp2pgsql. The shapefiles -# should exist, since this task depends on the "data" task, which -# downloads said shapefiles. -# -# Finally, we create the table for the demographic data (obtained from -# the geographic header records), and populate that table with the output -# of the sf1blocks2sql script. -# -db: data - # Ignore the result of dropdb when it fails. +newdb: +# Ignore the result of dropdb when it fails. dropdb -U $(DB_USER) $(DB_NAME) || true createdb -U $(DB_USER) $(DB_NAME) createlang -U $(DB_USER) plpgsql $(DB_NAME) - psql -d $(DB_NAME) \ - -U $(DB_USER) \ - -f $(PG_SHAREDIR)/contrib/postgis.sql - - psql -d $(DB_NAME) \ - -U $(DB_USER) \ - -f $(PG_SHAREDIR)/contrib/spatial_ref_sys.sql - - $(PG_BINDIR)/shp2pgsql -I \ - -s $(TIGER_SRID) \ - data/census2000/maryland/block/tl_2009_24_tabblock00.shp \ - tiger \ - | psql -U $(DB_USER) -d $(DB_NAME) - - psql -d $(DB_NAME) \ - -U $(DB_USER) \ - -f sql/create-sf1_blocks-table.sql - - bin/sf1blocks2sql src/Tests/Fixtures/SummaryFile1/mdgeo.uf1 sf1_blocks \ - | psql -U postgres -d $(DB_NAME) - + psql -d $(DB_NAME) \ + -U $(DB_USER) \ + -f $(POSTGIS_SQL) \ + > /dev/null + + psql -d $(DB_NAME) \ + -U $(DB_USER) \ + -f $(SPATIAL_REF_SYS_SQL) \ + > /dev/null + + +# This just runs the SQL script to create the sf1_blocks table. +sf1_blocks_table: + psql -d $(DB_NAME) \ + -U $(DB_USER) \ + -f sql/create-sf1_blocks-table.sql \ + > /dev/null + + +# Create the tiger_blocks table, and create its GiST index. Having the +# table already exist makes importing via shp2pgsql much easier. +# Any blocks file will work as an argument. +tiger_blocks_table: + $(SHP2PGSQL) \ + -p \ + -I \ + -s $(TIGER_SRID) \ + data/census2000/maryland/blocks/tl_2009_24_tabblock00.shp \ + tiger_blocks \ + | psql -U postgres -d $(DB_NAME) \ + > /dev/null + + +# Prepare the tiger_lines table, and create the GiST index on its +# geometry column. Any lines shapefile will do here. +tiger_lines_table: + $(SHP2PGSQL) \ + -p \ + -I \ + -s $(TIGER_SRID) \ + data/census2000/maryland/lines/tl_2009_24510_edges.shp \ + tiger_lines \ + | psql -U postgres -d $(DB_NAME) \ + > /dev/null + +# Add a unique index on the "tlid" column. + psql -U postgres \ + -d census \ + -f sql/create_tlid_unique_index.sql