Advertisement
Guest User

Untitled

a guest
Feb 20th, 2016
97
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.12 KB | None | 0 0
  1. ogr2ogr -f "PostGreSQL" PG:"host=127.0.0.1 user=youruser dbname=yourdb password=yourpass" "E:path_toaddresses.shp" -nln mcdon_addresses -nlt geometry
  2.  
  3. 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
  4.  
  5. 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
  6.  
  7. vacuum mcdon_addresses;
  8. vacuum mo_tracts_2010;
  9. vacuum mo_blocks_2010;
  10.  
  11. SELECT COUNT(*) FROM mcdon_addresses;
  12. -- 11979
  13.  
  14. CREATE TABLE mcdon_addresses_wtract AS
  15. SELECT
  16. a.wkb_geometry,
  17. a.route AS addr_route,
  18. a.box AS addr_box,
  19. a.new_add AS addr_new_add,
  20. a.prefix AS addr_prefix,
  21. a.rdname AS addr_rdname,
  22. a.road_name AS addr_road_name,
  23. a.city AS addr_city,
  24. a.state AS addr_state,
  25. a.zip AS addr_zip,
  26. t.statefp10 AS tr_statefp10,
  27. t.countyfp10 AS tr_countyfp10,
  28. t.tractce10 AS tr_tractce10,
  29. t.name10 AS tr_name10,
  30. t.pop90 AS tr_pop90,
  31. t.white90 AS tr_white90,
  32. t.black90 AS tr_black90,
  33. t.asian90 AS tr_asian90,
  34. t.amind90 AS tr_amind90,
  35. t.other90 AS tr_other90,
  36. t.hisp90 AS tr_hisp90
  37. FROM
  38. mcdon_addresses AS a,
  39. mo_tracts_2010 AS t
  40. WHERE
  41. ST_Contains(t.wkb_geometry, a.wkb_geometry);
  42.  
  43. vacuum mcdon_addresses_wtract;
  44.  
  45. select count(*) from mcdon_addresses_wtract;
  46. -- returns 11848
  47.  
  48. CREATE TABLE mcdon_addr_trct_and_blk AS
  49. SELECT
  50. a.*,
  51. b.pop90 AS blk_pop90,
  52. b.white90 AS blk_white90,
  53. b.black90 AS blk_black90,
  54. b.asian90 AS blk_asian90,
  55. b.amind90 AS blk_amind90,
  56. b.other90 AS blk_other90,
  57. b.hisp90 AS blk_hisp90
  58. FROM
  59. mcdon_addresses_wtract AS a,
  60. mo_blocks_2010 AS b
  61. WHERE
  62. ST_Contains(b.wkb_geometry, a.wkb_geometry);
  63.  
  64. vacuum mcdon_addr_trct_and_blk;
  65.  
  66. t.pop90 AS tr_pop90 -- would have been simply pop90
  67. b.pop90 AS blk_pop90 -- also would have been pop90 !
  68.  
  69. select count(*) from mcdon_addr_trct_and_blk;
  70. -- 11848 (thumbs up!)
  71.  
  72. 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