Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- ogr2ogr -f "PostGreSQL" PG:"host=127.0.0.1 user=youruser dbname=yourdb password=yourpass" "E:path_toaddresses.shp" -nln mcdon_addresses -nlt geometry
- ogr2ogr -f "PostGreSQL" PG:"host=127.0.0.1 user=youruser dbname=yourdb password=yourpass" "E:path_tost_tract10_spMoWest.shp" -nln mo_tracts_2010 -nlt geometry
- ogr2ogr -f "PostGreSQL" PG:"host=127.0.0.1 user=youruser dbname=yourdb password=yourpass" "E:path_tost_block10_spMoWest.shp" -nln mo_blocks_2010 -nlt geometry
- vacuum mcdon_addresses;
- vacuum mo_tracts_2010;
- vacuum mo_blocks_2010;
- SELECT COUNT(*) FROM mcdon_addresses;
- -- 11979
- CREATE TABLE mcdon_addresses_wtract AS
- SELECT
- a.wkb_geometry,
- a.route AS addr_route,
- a.box AS addr_box,
- a.new_add AS addr_new_add,
- a.prefix AS addr_prefix,
- a.rdname AS addr_rdname,
- a.road_name AS addr_road_name,
- a.city AS addr_city,
- a.state AS addr_state,
- a.zip AS addr_zip,
- t.statefp10 AS tr_statefp10,
- t.countyfp10 AS tr_countyfp10,
- t.tractce10 AS tr_tractce10,
- t.name10 AS tr_name10,
- t.pop90 AS tr_pop90,
- t.white90 AS tr_white90,
- t.black90 AS tr_black90,
- t.asian90 AS tr_asian90,
- t.amind90 AS tr_amind90,
- t.other90 AS tr_other90,
- t.hisp90 AS tr_hisp90
- FROM
- mcdon_addresses AS a,
- mo_tracts_2010 AS t
- WHERE
- ST_Contains(t.wkb_geometry, a.wkb_geometry);
- vacuum mcdon_addresses_wtract;
- select count(*) from mcdon_addresses_wtract;
- -- returns 11848
- CREATE TABLE mcdon_addr_trct_and_blk AS
- SELECT
- a.*,
- b.pop90 AS blk_pop90,
- b.white90 AS blk_white90,
- b.black90 AS blk_black90,
- b.asian90 AS blk_asian90,
- b.amind90 AS blk_amind90,
- b.other90 AS blk_other90,
- b.hisp90 AS blk_hisp90
- FROM
- mcdon_addresses_wtract AS a,
- mo_blocks_2010 AS b
- WHERE
- ST_Contains(b.wkb_geometry, a.wkb_geometry);
- vacuum mcdon_addr_trct_and_blk;
- t.pop90 AS tr_pop90 -- would have been simply pop90
- b.pop90 AS blk_pop90 -- also would have been pop90 !
- select count(*) from mcdon_addr_trct_and_blk;
- -- 11848 (thumbs up!)
- ogr2ogr -f "ESRI Shapefile" "E:path_toaddr_trct_blk.shp" PG:"host=127.0.0.1 user=youruser dbname=yourdb password=yourpass" "mcdon_addr_trct_and_blk"
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement