wards and pubs.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'id': 1142096055,\n",
       " 'name': \"Finch's\",\n",
       " 'o': (-0.0872945, 51.5217532),\n",
       " 'p': (-64.49651771904955, 83.83026593644367),\n",
       " 'point': <shapely.geometry.point.Point at 0x1118cdd90>}"
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from pyproj import Proj\n",
    "from shapely.geometry import Point, Polygon, LineString, MultiPolygon\n",
    "\n",
    "\n",
    "#pr = Proj('+init=EPSG:27700')\n",
    "\n",
    "# This is a stereographic projection which uses Finsbury Square as origin\n",
    "pr = Proj('+proj=stere +lat_0=51.5209993 +lon_0=-0.0863623 +k=1 +x_0=0 +y_0=0 +a=6371000 +b=6371000 +units=m +no_defs')\n",
    "\n",
    "file_name=\"data/pubs.tsv\"\n",
    "pubs = dict([(int(x[0]), {\n",
    "                'name': x[4], \n",
    "                'o': (float(x[1]), float(x[2])),\n",
    "                'p': pr(float(x[1]), float(x[2])),\n",
    "                'point': Point(float(x[1]), float(x[2])),\n",
    "                'id': int(x[0])}) \n",
    "             for x in [line.strip().split('\\t') \n",
    "                       for line in file(file_name).readlines()] if len(x)>4])\n",
    "len(pubs)\n",
    "pubs[1142096055]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{u'geometry': {u'coordinates': [[[-0.330679062942453, 51.329011010602905],\n",
       "    [-0.330594448736231, 51.329088039878236],\n",
       "    [-0.330506117851447, 51.329148829506465],\n",
       "    [-0.330409211490547, 51.32920859757066],\n",
       "    [-0.330292121849854, 51.32927077559537],\n",
       "    [-0.329913461974029, 51.32946051543075],\n",
       "    [-0.329633594619526, 51.32965705855314],\n",
       "    [-0.329428784834638, 51.32984208201042],\n",
       "    [-0.329320915687227, 51.32996643925324],\n",
       "    [-0.32921448062474, 51.330090816854394],\n",
       "    [-0.32911632355624, 51.33026387227744],\n",
       "    [-0.329112955750983, 51.3303960148678],\n",
       "    [-0.329171167711951, 51.33072957000635],\n",
       "    [-0.329191302316349, 51.33096546228584],\n",
       "    [-0.329302910546211, 51.331528189994685],\n",
       "    [-0.329362771534731, 51.33161896942137],\n",
       "    [-0.329554552333886, 51.33187169784059],\n",
       "    [-0.329578988615119, 51.33218678585108],\n",
       "    [-0.329594210640675, 51.33228142478945],\n",
       "    [-0.329690000087309, 51.332489619653416],\n",
       "    [-0.32971930401878, 51.332670787926034],\n",
       "    [-0.329718711471269, 51.33276610057187],\n",
       "    [-0.329713651327465, 51.332865845784674],\n",
       "    [-0.329643374963557, 51.33349612100018],\n",
       "    [-0.329628972694324, 51.33385292004089],\n",
       "    [-0.329550496490717, 51.33415574920923],\n",
       "    [-0.329484384839356, 51.33435534029025],\n",
       "    [-0.329546433368672, 51.334465035124275],\n",
       "    [-0.329431949069611, 51.33485008214592],\n",
       "    [-0.329410734967829, 51.33511775766491],\n",
       "    [-0.329359230972181, 51.33531036295659],\n",
       "    [-0.329208223769788, 51.33563373943428],\n",
       "    [-0.329197343121298, 51.33593303622248],\n",
       "    [-0.329212414905641, 51.33630824062638],\n",
       "    [-0.329192134820813, 51.33670811987971],\n",
       "    [-0.329041537123449, 51.337572853475976],\n",
       "    [-0.328900444301525, 51.33817603921421],\n",
       "    [-0.328869226896171, 51.33867917658232],\n",
       "    [-0.328882308957264, 51.339030072608615],\n",
       "    [-0.328865515201647, 51.339412915553105],\n",
       "    [-0.328870784697311, 51.33958384913406],\n",
       "    [-0.328876250331204, 51.339749389979964],\n",
       "    [-0.328959320052221, 51.3401102769569],\n",
       "    [-0.328979269794097, 51.34023286018948],\n",
       "    [-0.32890917222982, 51.34077860696534],\n",
       "    [-0.328780550871525, 51.34107802221065],\n",
       "    [-0.328542180279971, 51.34174996116788],\n",
       "    [-0.328525414020483, 51.3417766995562],\n",
       "    [-0.328346139166793, 51.34204751460602],\n",
       "    [-0.328201270790769, 51.3423197196975],\n",
       "    [-0.328010180987364, 51.34295456362381],\n",
       "    [-0.327944192320452, 51.34334749508946],\n",
       "    [-0.327906806306604, 51.34374353255681],\n",
       "    [-0.327796155987095, 51.34390381877997],\n",
       "    [-0.327480812238024, 51.34400722717033],\n",
       "    [-0.328177391612962, 51.34404145085986],\n",
       "    [-0.328483545945542, 51.344545805252174],\n",
       "    [-0.328666818347213, 51.344915315970134],\n",
       "    [-0.328898463686148, 51.345297206761195],\n",
       "    [-0.329205768953539, 51.34588880312093],\n",
       "    [-0.329345842260152, 51.34614349124877],\n",
       "    [-0.329970793224286, 51.34728006667857],\n",
       "    [-0.330327280093989, 51.34790203628056],\n",
       "    [-0.330533862924573, 51.34842115203906],\n",
       "    [-0.330367267889462, 51.34857974401673],\n",
       "    [-0.33021561010705, 51.34876192926616],\n",
       "    [-0.330102635557406, 51.34898603141627],\n",
       "    [-0.329936977146583, 51.349237259224054],\n",
       "    [-0.329752002020655, 51.349545763584665],\n",
       "    [-0.329563268511802, 51.34991806103497],\n",
       "    [-0.329472964225532, 51.350071444835216],\n",
       "    [-0.329311211981746, 51.350333518494104],\n",
       "    [-0.32903599281449, 51.35071447339452],\n",
       "    [-0.328765267701086, 51.351129663380995],\n",
       "    [-0.328353106101917, 51.35172178590493],\n",
       "    [-0.328046517372798, 51.352135563199745],\n",
       "    [-0.327947015118701, 51.35218629969403],\n",
       "    [-0.327759754660003, 51.35228074605213],\n",
       "    [-0.327275614659091, 51.352522027863074],\n",
       "    [-0.326841800917161, 51.352800895868654],\n",
       "    [-0.326479986169664, 51.353034928954614],\n",
       "    [-0.325766123154619, 51.35335295611828],\n",
       "    [-0.324837712030583, 51.353725463393076],\n",
       "    [-0.324453936901546, 51.35389173226],\n",
       "    [-0.324158786694635, 51.35410872665401],\n",
       "    [-0.323752729174304, 51.354413159467434],\n",
       "    [-0.323535077753691, 51.35467082833569],\n",
       "    [-0.323184853046418, 51.355176591634965],\n",
       "    [-0.322861962024969, 51.35552357741356],\n",
       "    [-0.322576101914298, 51.35583961869267],\n",
       "    [-0.322259755120477, 51.35620378233469],\n",
       "    [-0.321953015303272, 51.35661934015992],\n",
       "    [-0.321529264210917, 51.35721037359525],\n",
       "    [-0.32121760273532, 51.35768161266449],\n",
       "    [-0.320948200338684, 51.35809680365024],\n",
       "    [-0.320744893153151, 51.35839334081395],\n",
       "    [-0.320259445798069, 51.359099488898266],\n",
       "    [-0.320079350800145, 51.35942873082388],\n",
       "    [-0.320178339947984, 51.35954975087501],\n",
       "    [-0.320165178404721, 51.35959542404538],\n",
       "    [-0.319957006726642, 51.35963200602732],\n",
       "    [-0.31979241444019, 51.359929995977105],\n",
       "    [-0.319546632385726, 51.36048400740659],\n",
       "    [-0.319408617907466, 51.360919064444396],\n",
       "    [-0.319187512741341, 51.36162270500862],\n",
       "    [-0.319060245805666, 51.362313303504756],\n",
       "    [-0.318285409078045, 51.36350628265066],\n",
       "    [-0.318304733628236, 51.36368461201952],\n",
       "    [-0.318275986396557, 51.36368509862089],\n",
       "    [-0.316898872391401, 51.36367251664156],\n",
       "    [-0.316629141652739, 51.363702812882245],\n",
       "    [-0.316560227340392, 51.363701822656395],\n",
       "    [-0.316498590473247, 51.36369823920765],\n",
       "    [-0.316562845042351, 51.36359125220441],\n",
       "    [-0.315497844471691, 51.36356785084888],\n",
       "    [-0.31497437125262, 51.36354503571201],\n",
       "    [-0.313400209206441, 51.36346123506726],\n",
       "    [-0.312318547134839, 51.363422277604045],\n",
       "    [-0.311025097733082, 51.36336137277011],\n",
       "    [-0.309737891203448, 51.363287054899544],\n",
       "    [-0.308904082696092, 51.363266029571946],\n",
       "    [-0.30842535520976, 51.363237536228844],\n",
       "    [-0.307869697178102, 51.36319174323868],\n",
       "    [-0.30752859878378, 51.36317062951251],\n",
       "    [-0.307264299029795, 51.36317040803606],\n",
       "    [-0.306728222798701, 51.36317794851318],\n",
       "    [-0.306303898268055, 51.36319339694987],\n",
       "    [-0.30589399672871, 51.36320725390605],\n",
       "    [-0.305299165359023, 51.36324990835698],\n",
       "    [-0.303890867602543, 51.36330416625408],\n",
       "    [-0.303040911715541, 51.36333142510303],\n",
       "    [-0.302168894334417, 51.363411414059925],\n",
       "    [-0.30154404068005, 51.363450017355646],\n",
       "    [-0.300459502192416, 51.363527811636686],\n",
       "    [-0.300006596738068, 51.36357789357511],\n",
       "    [-0.299541738708729, 51.363640389840285],\n",
       "    [-0.29930672081789, 51.36366395603959],\n",
       "    [-0.298553351485872, 51.36372046131762],\n",
       "    [-0.298527475541803, 51.36372098477151],\n",
       "    [-0.298036801778111, 51.3637435329516],\n",
       "    [-0.297627895987718, 51.363769065781824],\n",
       "    [-0.297096652708767, 51.363762244361276],\n",
       "    [-0.296894618783457, 51.36374851647035],\n",
       "    [-0.296697025464506, 51.36373125577347],\n",
       "    [-0.296301070620289, 51.363678737112814],\n",
       "    [-0.296265214629983, 51.36375465232255],\n",
       "    [-0.296165655370688, 51.36388269709672],\n",
       "    [-0.296012907593596, 51.36401176669595],\n",
       "    [-0.295845135812233, 51.36419727059796],\n",
       "    [-0.294916710600628, 51.36556321848823],\n",
       "    [-0.294162672239175, 51.36540476410622],\n",
       "    [-0.294139532980609, 51.36544759140917],\n",
       "    [-0.293923130240196, 51.36585899663533],\n",
       "    [-0.293794244827275, 51.365886795090645],\n",
       "    [-0.29395060588565, 51.36624068066671],\n",
       "    [-0.293888999604844, 51.36635219034912],\n",
       "    [-0.292691668215232, 51.36610004252014],\n",
       "    [-0.292694106727627, 51.36607310041628],\n",
       "    [-0.292229806843134, 51.36561760521599],\n",
       "    [-0.292136824804918, 51.365452585421025],\n",
       "    [-0.291784225373345, 51.36496814196446],\n",
       "    [-0.291766532855606, 51.36490313763907],\n",
       "    [-0.291755183006384, 51.364860707165974],\n",
       "    [-0.290845677059296, 51.36467298665931],\n",
       "    [-0.290228593389514, 51.36454078458316],\n",
       "    [-0.290030367119905, 51.36450192128183],\n",
       "    [-0.289874041028482, 51.364495143202916],\n",
       "    [-0.28976325232328, 51.36453848871602],\n",
       "    [-0.289618232355378, 51.36465237553821],\n",
       "    [-0.289399110385555, 51.364865896813704],\n",
       "    [-0.289203330982519, 51.36503119877188],\n",
       "    [-0.288606273977691, 51.36547930020749],\n",
       "    [-0.288098218373541, 51.36523537339522],\n",
       "    [-0.287606732287988, 51.36497100365596],\n",
       "    [-0.287423970183523, 51.36490268654534],\n",
       "    [-0.287270258489396, 51.36482580129217],\n",
       "    [-0.286959931669728, 51.36467288697118],\n",
       "    [-0.286735475223243, 51.36456708963517],\n",
       "    [-0.286602966151974, 51.364499506108466],\n",
       "    [-0.286548381980651, 51.36446093906808],\n",
       "    [-0.286620248021387, 51.364382855637636],\n",
       "    [-0.285726159243823, 51.36405233901505],\n",
       "    [-0.285769002706404, 51.3640205927138],\n",
       "    [-0.285955209819343, 51.36384256680158],\n",
       "    [-0.286026706790759, 51.36377437014037],\n",
       "    [-0.28632801543817, 51.363514396747945],\n",
       "    [-0.286534492170347, 51.363293502142234],\n",
       "    [-0.28667000872075, 51.36312642402789],\n",
       "    [-0.286780715453229, 51.36296977401728],\n",
       "    [-0.286839724771548, 51.36288970392315],\n",
       "    [-0.286889583342469, 51.362823888081806],\n",
       "    [-0.28704479405179, 51.362667888586465],\n",
       "    [-0.287223336967318, 51.36254100618897],\n",
       "    [-0.287372556498594, 51.362468549454675],\n",
       "    [-0.287497631853592, 51.36242721355338],\n",
       "    [-0.28770664517004, 51.36236911926102],\n",
       "    [-0.287918328182864, 51.362316459147074],\n",
       "    [-0.28815434972806, 51.36226595280225],\n",
       "    [-0.288178459397226, 51.362235730398744],\n",
       "    [-0.288301979065533, 51.36208196395354],\n",
       "    [-0.287966523957395, 51.36194846863353],\n",
       "    [-0.287960814914869, 51.361947485947475],\n",
       "    [-0.28762255756466, 51.36181215016591],\n",
       "    [-0.287595682747195, 51.36180096627366],\n",
       "    [-0.287237423338822, 51.361701307229175],\n",
       "    [-0.287027028931216, 51.36164247722088],\n",
       "    [-0.287021353478179, 51.36164059572393],\n",
       "    [-0.286524015278682, 51.36149573652897],\n",
       "    [-0.286165525813508, 51.36140236560232],\n",
       "    [-0.285847127288408, 51.36135094595909],\n",
       "    [-0.285425834961252, 51.36128543000663],\n",
       "    [-0.285101065668824, 51.361212332938514],\n",
       "    [-0.285360539271765, 51.3610722499134],\n",
       "    [-0.286054514398571, 51.360297354333106],\n",
       "    [-0.286413732473064, 51.35998607105688],\n",
       "    [-0.286611840908848, 51.359796527839464],\n",
       "    [-0.286702931819206, 51.359703438157446],\n",
       "    [-0.286857468005672, 51.359526746129355],\n",
       "    [-0.286996708635952, 51.35937500941728],\n",
       "    [-0.288097795281945, 51.359043991775295],\n",
       "    [-0.288996388351274, 51.358711805058384],\n",
       "    [-0.289236119546597, 51.358445529329344],\n",
       "    [-0.289368845131762, 51.35823704155323],\n",
       "    [-0.289486674644875, 51.35808139223847],\n",
       "    [-0.289610312356681, 51.357924029093816],\n",
       "    [-0.289794813382644, 51.3577135989836],\n",
       "    [-0.289881821437393, 51.35761415245672],\n",
       "    [-0.290223167540334, 51.35728011531095],\n",
       "    [-0.290623199278082, 51.35691276177518],\n",
       "    [-0.290898628300749, 51.35669016784578],\n",
       "    [-0.290939791630382, 51.35666468995205],\n",
       "    [-0.291108817798991, 51.356561043510425],\n",
       "    [-0.291130622667594, 51.3564768315063],\n",
       "    [-0.291203689906098, 51.35609571345367],\n",
       "    [-0.291209032976251, 51.35606791444844],\n",
       "    [-0.291243095200363, 51.35588586228224],\n",
       "    [-0.291251555687757, 51.3556584738552],\n",
       "    [-0.291224371747488, 51.355308266633536],\n",
       "    [-0.291224011430227, 51.35520214912961],\n",
       "    [-0.291346773062895, 51.35495214796771],\n",
       "    [-0.291499877648405, 51.354697193672216],\n",
       "    [-0.291725018288241, 51.35443609560903],\n",
       "    [-0.292003299521232, 51.35417487266716],\n",
       "    [-0.292081180469964, 51.35412744830833],\n",
       "    [-0.292307511714818, 51.35398866536462],\n",
       "    [-0.29242017706709, 51.35389408725863],\n",
       "    [-0.292524229628563, 51.35379938350807],\n",
       "    [-0.292684698582744, 51.35361647552649],\n",
       "    [-0.292782176732627, 51.353428153061806],\n",
       "    [-0.292772712159673, 51.35291004339951],\n",
       "    [-0.292783734219189, 51.3526907854388],\n",
       "    [-0.292799759208744, 51.35260738812107],\n",
       "    [-0.29281995655589, 51.35252764863183],\n",
       "    [-0.292930449011733, 51.35237548595083],\n",
       "    [-0.293176897917912, 51.35219742721934],\n",
       "    [-0.293319465282057, 51.35210957872811],\n",
       "    [-0.293470677489115, 51.352020956745264],\n",
       "    [-0.294172454767205, 51.35164719547864],\n",
       "    [-0.294349966797501, 51.35154636575084],\n",
       "    [-0.294519266614886, 51.351434625117086],\n",
       "    [-0.294840606711624, 51.35121179114641],\n",
       "    [-0.294991178411417, 51.35110157573786],\n",
       "    [-0.295264581002036, 51.35085466220706],\n",
       "    [-0.295638114849462, 51.35054086144621],\n",
       "    [-0.296039348929474, 51.350254440134904],\n",
       "    [-0.296476626418785, 51.349964045341835],\n",
       "    [-0.297268233435071, 51.34956729207264],\n",
       "    [-0.297558910099395, 51.349396161873656],\n",
       "    [-0.297795048302183, 51.34926290649267],\n",
       "    [-0.29782144896842, 51.34924800276254],\n",
       "    [-0.298158715314547, 51.34905956294252],\n",
       "    [-0.29841540582426, 51.34887534728589],\n",
       "    [-0.298674228528117, 51.34867227765281],\n",
       "    [-0.298993199461607, 51.34851234593999],\n",
       "    [-0.29919952204469, 51.348369662229715],\n",
       "    [-0.299452692805216, 51.348125143030934],\n",
       "    [-0.299949117814642, 51.3476710300914],\n",
       "    [-0.300100658211388, 51.34749517650769],\n",
       "    [-0.300184183318634, 51.347372291169826],\n",
       "    [-0.300270744619537, 51.347244953545705],\n",
       "    [-0.30054770319883, 51.34678405633227],\n",
       "    [-0.300672895356338, 51.34654397252912],\n",
       "    [-0.300788671674013, 51.346286666128485],\n",
       "    [-0.301045714577185, 51.345742747212135],\n",
       "    [-0.301176010266212, 51.34551982275818],\n",
       "    [-0.301337977184917, 51.345294659669236],\n",
       "    [-0.301550942339489, 51.345088220956185],\n",
       "    [-0.301790840551982, 51.34489116489336],\n",
       "    [-0.301944916023182, 51.34480167432813],\n",
       "    [-0.302021103110804, 51.34476051366637],\n",
       "    [-0.302311133706791, 51.34460554904585],\n",
       "    [-0.302536516309658, 51.344451445413924],\n",
       "    [-0.303287638580883, 51.34378608668764],\n",
       "    [-0.303431454738399, 51.34366317308032],\n",
       "    [-0.303699366969436, 51.34344583598246],\n",
       "    [-0.303806342011612, 51.343348466846045],\n",
       "    [-0.303903204776606, 51.34325274970374],\n",
       "    [-0.304048568303663, 51.3430489246438],\n",
       "    [-0.304146133878932, 51.34281742970631],\n",
       "    [-0.304215536949321, 51.34260980690898],\n",
       "    [-0.304329184654699, 51.34237045132315],\n",
       "    [-0.304396795146305, 51.342211362331014],\n",
       "    [-0.30446170092791, 51.34204773786317],\n",
       "    [-0.304518327734426, 51.34187500095417],\n",
       "    [-0.304569214026287, 51.34170218092605],\n",
       "    [-0.30461846292663, 51.341379161268684],\n",
       "    [-0.304663930759492, 51.3411586021463],\n",
       "    [-0.304753897793453, 51.340860451676285],\n",
       "    [-0.304810962200981, 51.34071469865024],\n",
       "    [-0.30493950599659, 51.34042160259911],\n",
       "    [-0.305010918982431, 51.34027605711182],\n",
       "    [-0.305087673569484, 51.34014137995753],\n",
       "    [-0.305341258510224, 51.339727793771026],\n",
       "    [-0.305437971784711, 51.339480098747764],\n",
       "    [-0.305439437739343, 51.33920674578781],\n",
       "    [-0.305381001051333, 51.3389226343555],\n",
       "    [-0.305379739714676, 51.33868431294534],\n",
       "    [-0.305321910283473, 51.33842269167236],\n",
       "    [-0.305296398005862, 51.33821909040641],\n",
       "    [-0.305245152864049, 51.338129322485926],\n",
       "    [-0.305080313565159, 51.33788863413644],\n",
       "    [-0.304980761122958, 51.33770734206322],\n",
       "    [-0.304883283574739, 51.33754766210264],\n",
       "    [-0.304956216984938, 51.33736077276589],\n",
       "    [-0.304961019291707, 51.33707487838324],\n",
       "    [-0.304949995413292, 51.33675098614315],\n",
       "    [-0.305052831836229, 51.3364539207886],\n",
       "    [-0.305196979115416, 51.33616554640919],\n",
       "    [-0.305279365948208, 51.336033648349265],\n",
       "    [-0.305710121335248, 51.33540860122101],\n",
       "    [-0.305888097848156, 51.33525380545954],\n",
       "    [-0.306031956488495, 51.33516775884644],\n",
       "    [-0.306215856862412, 51.33508588813518],\n",
       "    [-0.30662399586401, 51.334918232759314],\n",
       "    [-0.307014814945081, 51.33475302340159],\n",
       "    [-0.307491652950809, 51.33459085427453],\n",
       "    [-0.30763597191626, 51.33460912623915],\n",
       "    [-0.30772952941814, 51.334485481194314],\n",
       "    [-0.308351632141362, 51.33408170848809],\n",
       "    [-0.308757264199398, 51.33374764687651],\n",
       "    [-0.309137156732768, 51.33348874993342],\n",
       "    [-0.309644387193975, 51.33308511057124],\n",
       "    [-0.309881033907439, 51.33293475253331],\n",
       "    [-0.310419197417255, 51.33266554557943],\n",
       "    [-0.310953248609539, 51.332390881324805],\n",
       "    [-0.311603686293558, 51.332152064234116],\n",
       "    [-0.312009022683188, 51.332098555544526],\n",
       "    [-0.312422142948081, 51.3319894036455],\n",
       "    [-0.312603825725588, 51.33192817395268],\n",
       "    [-0.312776932709508, 51.331865921186534],\n",
       "    [-0.313094451004898, 51.331741001605465],\n",
       "    [-0.313716437388347, 51.3314163336281],\n",
       "    [-0.313955964617329, 51.33126510964038],\n",
       "    [-0.314268361300155, 51.331044791889816],\n",
       "    [-0.315047048628603, 51.33051734144213],\n",
       "    [-0.315210035016298, 51.33037850263434],\n",
       "    [-0.31597891155355, 51.329687240071486],\n",
       "    [-0.316142536657492, 51.329569991338744],\n",
       "    [-0.316765190071445, 51.32914729784507],\n",
       "    [-0.316850274716521, 51.32909726293057],\n",
       "    [-0.317148217576959, 51.328918095695236],\n",
       "    [-0.317831324514175, 51.32852953770007],\n",
       "    [-0.318394041495096, 51.32825345298935],\n",
       "    [-0.318978993789515, 51.32795790107236],\n",
       "    [-0.319307299160142, 51.32781153743493],\n",
       "    [-0.319608586178474, 51.32773672583376],\n",
       "    [-0.320082318825753, 51.32765629370738],\n",
       "    [-0.321236031641869, 51.3275073727147],\n",
       "    [-0.321599582250038, 51.32745682913346],\n",
       "    [-0.321774381195321, 51.32742606141491],\n",
       "    [-0.32210214037995, 51.327333637594705],\n",
       "    [-0.322648093964269, 51.32720117357336],\n",
       "    [-0.323094304718351, 51.327087961896],\n",
       "    [-0.323095739294378, 51.32708798243569],\n",
       "    [-0.32355043925606, 51.32689935298098],\n",
       "    [-0.323836596008216, 51.32676676153044],\n",
       "    [-0.323970443820737, 51.32671741933837],\n",
       "    [-0.32402864555579, 51.326695770773995],\n",
       "    [-0.32426001048832, 51.326648722906576],\n",
       "    [-0.324448789957882, 51.32662804307077],\n",
       "    [-0.324638938320394, 51.3266091810266],\n",
       "    [-0.325113522723552, 51.326583595487605],\n",
       "    [-0.325604895834516, 51.32656993838298],\n",
       "    [-0.326003933138142, 51.326569347625934],\n",
       "    [-0.326706410732271, 51.32659197548505],\n",
       "    [-0.326981711509264, 51.32663907265173],\n",
       "    [-0.327231681204934, 51.32667231853357],\n",
       "    [-0.327323035627756, 51.32668621284312],\n",
       "    [-0.327597597238846, 51.32631963991207],\n",
       "    [-0.328314309328093, 51.32646386104324],\n",
       "    [-0.328926092286686, 51.32680981267641],\n",
       "    [-0.329251043715115, 51.32699160209771],\n",
       "    [-0.329588647712956, 51.32718076510295],\n",
       "    [-0.329590082296634, 51.32718078556312],\n",
       "    [-0.330139023428105, 51.327518640399624],\n",
       "    [-0.330157246912207, 51.32757015782663],\n",
       "    [-0.330339655624747, 51.32788300173333],\n",
       "    [-0.33039873068244, 51.327995351641576],\n",
       "    [-0.330480430272456, 51.328196151124025],\n",
       "    [-0.33059043133774, 51.328408145009035],\n",
       "    [-0.330674185780984, 51.32863145502641],\n",
       "    [-0.330756473741289, 51.32885564332769],\n",
       "    [-0.330719366232742, 51.328928853465115],\n",
       "    [-0.330679062942453, 51.329011010602905]]],\n",
       "  u'type': u'Polygon'},\n",
       " u'properties': {u'BOROUGH': u'Kingston upon Thames',\n",
       "  u'GSS_CODE': u'E05000405',\n",
       "  u'HECTARES': 755.173,\n",
       "  u'LB_GSS_CD': u'E09000021',\n",
       "  u'NAME': u'Chessington South',\n",
       "  u'NONLD_AREA': 0.0,\n",
       "  u'POLY_ID': 50840.0},\n",
       " u'type': u'Feature'}"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from shapely.geometry import asShape, mapping\n",
    "import json\n",
    "\n",
    "wards = json.load(file('data/wards.json'))['features']\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Chessington South 3\n"
     ]
    }
   ],
   "source": [
    "for ward in wards[0:1]:\n",
    "    ward_shape = asShape(ward['geometry'])\n",
    "    count = 0\n",
    "    for pub in pubs.values():\n",
    "        if ward_shape.contains(pub['point']):\n",
    "            count += 1\n",
    "    print ward['properties']['NAME'], count"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<rtree.index.Index at 0x115ada150>"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from rtree import index\n",
    "tree = index.Index()\n",
    "for pub in pubs.values():\n",
    "    tree.insert(pub['id'], (pub['o'][0], pub['o'][1],pub['o'][0], pub['o'][1]))\n",
    "tree"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(-0.330756473741289, 51.32631963991207, -0.285101065668824, 51.36635219034912)"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "asShape(wards[0]['geometry']).bounds"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Chessington South 3\n",
      "Tolworth and Hook Rise 1\n",
      "Berrylands 5\n",
      "Alexandra 2\n",
      "Beverley 2\n",
      "Coombe Hill 1\n",
      "Chessington North and Hook 2\n",
      "Surbiton Hill 1\n",
      "Old Malden 2\n",
      "St. Mark's 12\n",
      "Grove 16\n",
      "Canbury 7\n",
      "Norbiton 6\n",
      "Coombe Vale 2\n",
      "St. James 2\n",
      "Tudor 3\n",
      "Coulsdon East 2\n",
      "Selsdon and Ballards 1\n",
      "Coulsdon West 2\n",
      "Waddon 3\n",
      "Kenley 2\n",
      "Purley 3\n",
      "Sanderstead 0\n",
      "Heathfield 6\n",
      "Fairfield 24\n",
      "Broad Green 7\n",
      "West Thornton 5\n",
      "Bensham Manor 2\n",
      "Norbury 0\n",
      "New Addington 0\n",
      "Croham 2\n",
      "Fieldway 0\n",
      "Shirley 1\n",
      "Selhurst 9\n",
      "Ashburton 0\n",
      "Woodside 2\n",
      "Thornton Heath 1\n",
      "Upper Norwood 5\n",
      "South Norwood 7\n",
      "Addiscombe 7\n",
      "Darwin 7\n",
      "Hayes and Coney Hall 2\n",
      "Bromley Common and Keston 5\n",
      "Chelsfield and Pratts Bottom 7\n",
      "Biggin Hill 1\n",
      "West Wickham 4\n",
      "Clock House 2\n",
      "Kelsey and Eden Park 3\n",
      "Farnborough and Crofton 4\n",
      "Shortlands 0\n",
      "Bromley Town 14\n",
      "Bickley 3\n",
      "Petts Wood and Knoll 4\n",
      "Crystal Palace 5\n",
      "Penge and Cator 10\n",
      "Copers Cope 10\n",
      "Plaistow and Sundridge 8\n",
      "Chislehurst 10\n",
      "Mottingham and Chislehurst North 1\n",
      "Orpington 1\n",
      "Cray Valley West 0\n",
      "Cray Valley East 3\n",
      "Bedfont 3\n",
      "Hanworth 2\n",
      "Cranford 2\n",
      "Syon 15\n",
      "Heston West 2\n",
      "Heston East 2\n",
      "Osterley and Spring Grove 2\n",
      "Brentford 6\n",
      "Chiswick Homefields 5\n",
      "Feltham West 2\n",
      "Hanworth Park 4\n",
      "Feltham North 2\n",
      "Hounslow Heath 4\n",
      "Hounslow West 4\n",
      "Heston Central 2\n",
      "Hounslow South 1\n",
      "Isleworth 8\n",
      "Hounslow Central 10\n",
      "Chiswick Riverside 4\n",
      "Turnham Green 11\n",
      "Norwood Green 3\n",
      "Southall Green 4\n",
      "Northolt West End 1\n",
      "Dormers Wells 1\n",
      "Greenford Broadway 2\n",
      "North Greenford 2\n",
      "Northfield 6\n",
      "East Acton 4\n",
      "Southall Broadway 3\n",
      "Elthorne 8\n",
      "Walpole 6\n",
      "Lady Margaret 0\n",
      "Northolt Mandeville 2\n",
      "Hobbayne 1\n",
      "Cleveland 2\n",
      "Greenford Green 3\n",
      "Perivale 1\n",
      "Ealing Common 7\n",
      "Ealing Broadway 11\n",
      "South Acton 8\n",
      "Southfield 4\n",
      "Acton Central 7\n",
      "Hanger Hill 2\n",
      "Upminster 7\n",
      "Rainham and Wennington 5\n",
      "South Hornchurch 1\n",
      "Elm Park 1\n",
      "Brooklands 0\n",
      "Romford Town 4\n",
      "Mawneys 1\n",
      "Pettits 2\n",
      "Hacton 2\n",
      "St. Andrew's 4\n",
      "Emerson Park 1\n",
      "Squirrel's Heath 2\n",
      "Harold Wood 2\n",
      "Cranham 1\n",
      "Havering Park 3\n",
      "Heaton 1\n",
      "Gooshays 4\n",
      "Hylands 1\n",
      "Heathrow Villages 15\n",
      "Harefield 1\n",
      "West Drayton 0\n",
      "Yiewsley 4\n",
      "Uxbridge South 10\n",
      "Brunel 7\n",
      "Uxbridge North 6\n",
      "Hillingdon East 0\n",
      "Ickenham 3\n",
      "West Ruislip 4\n",
      "Northwood 2\n",
      "South Ruislip 2\n",
      "Manor 2\n",
      "Eastcote and East Ruislip 2\n",
      "Northwood Hills 0\n",
      "Pinkwell 1\n",
      "Botwell 1\n",
      "Charville 0\n",
      "Townfield 2\n",
      "Barnhill 0\n",
      "Yeading 1\n",
      "Cavendish 1\n",
      "Roxeth 3\n",
      "Harrow on the Hill 7\n",
      "Pinner 2\n",
      "Pinner South 1\n",
      "Greenhill 7\n",
      "Headstone North 0\n",
      "Marlborough 5\n",
      "Harrow Weald 3\n",
      "Stanmore Park 2\n",
      "Canons 0\n",
      "Rayners Lane 0\n",
      "Roxbourne 1\n",
      "West Harrow 0\n",
      "Hatch End 2\n",
      "Headstone South 0\n",
      "Kenton West 2\n",
      "Wealdstone 1\n",
      "Belmont 1\n",
      "Kenton East 0\n",
      "Queensbury 0\n",
      "Edgware 2\n",
      "Northwick Park 1\n",
      "Wembley Central 6\n",
      "Preston 2\n",
      "Stonebridge 1\n",
      "Welsh Harp 0\n",
      "Barnhill 2\n",
      "Fryent 0\n",
      "Sudbury 2\n",
      "Alperton 4\n",
      "Kensal Green 5\n",
      "Harlesden 2\n",
      "Willesden Green 3\n",
      "Queens Park 8\n",
      "Brondesbury Park 2\n",
      "Kilburn 10\n",
      "Tokyngton 5\n",
      "Kenton 0\n",
      "Dudden Hill 0\n",
      "Dollis Hill 2\n",
      "Mapesbury 2\n",
      "Queensbury 0\n",
      "Underhill 3\n",
      "High Barnet 8\n",
      "West Hendon 0\n",
      "Golders Green 0\n",
      "Colindale 3\n",
      "Childs Hill 8\n",
      "Garden Suburb 2\n",
      "Finchley Church End 2\n",
      "East Finchley 3\n",
      "Edgware 1\n",
      "Mill Hill 2\n",
      "Hale 2\n",
      "Totteridge 3\n",
      "Oakleigh 3\n",
      "Woodhouse 4\n",
      "Coppetts 2\n",
      "Brunswick Park 2\n",
      "East Barnet 7\n",
      "Burnt Oak 1\n",
      "Hendon 3\n",
      "West Finchley 2\n",
      "Streatham South 3\n",
      "St. Leonard's 7\n",
      "Streatham Wells 1\n",
      "Knight's Hill 3\n",
      "Thornton 1\n",
      "Streatham Hill 1\n",
      "Clapham Common 4\n",
      "Brixton Hill 7\n",
      "Tulse Hill 1\n",
      "Coldharbour 10\n",
      "Ferndale 6\n",
      "Larkhall 9\n",
      "Oval 14\n",
      "Vassall 4\n",
      "Gipsy Hill 5\n",
      "Thurlow Park 4\n",
      "Herne Hill 7\n",
      "Bishop's 24\n",
      "Clapham Town 16\n",
      "Stockwell 5\n",
      "Prince's 16\n",
      "College 3\n",
      "Riverside 20\n",
      "Village 2\n",
      "South Camberwell 4\n",
      "East Dulwich 10\n",
      "Peckham Rye 4\n",
      "Camberwell Green 12\n",
      "Brunswick Park 6\n",
      "The Lane 11\n",
      "Nunhead 5\n",
      "Peckham 5\n",
      "Newington 6\n",
      "Faraday 7\n",
      "East Walworth 8\n",
      "Livesey 6\n",
      "South Bermondsey 4\n",
      "Cathedrals 43\n",
      "Grange 16\n",
      "Rotherhithe 13\n",
      "Surrey Docks 10\n",
      "Chaucer 21\n",
      "Bellingham 3\n",
      "Telegraph Hill 2\n",
      "Downham 1\n",
      "Whitefoot 0\n",
      "Blackheath 12\n",
      "Sydenham 5\n",
      "Forest Hill 10\n",
      "Perry Vale 5\n",
      "Rushey Green 3\n",
      "Catford South 0\n",
      "Crofton Park 6\n",
      "Ladywell 2\n",
      "Lewisham Central 3\n",
      "Brockley 15\n",
      "New Cross 7\n",
      "Grove Park 1\n",
      "Lee Green 6\n",
      "Evelyn 6\n",
      "Middle Park and Sutcliffe 1\n",
      "Coldharbour and New Eltham 1\n",
      "Eltham South 9\n",
      "Shooters Hill 2\n",
      "Peninsula 15\n",
      "Woolwich Riverside 8\n",
      "Greenwich West 28\n",
      "Eltham West 0\n",
      "Blackheath Westcombe 1\n",
      "Kidbrooke with Hornfair 3\n",
      "Eltham North 0\n",
      "Charlton 0\n",
      "Woolwich Common 4\n",
      "Plumstead 3\n",
      "Glyndon 2\n",
      "Thamesmead Moorings 2\n",
      "Abbey Wood 4\n",
      "Longlands 3\n",
      "Blackfen and Lamorbey 2\n",
      "Cray Meadows 1\n",
      "Sidcup 3\n",
      "St. Mary's 6\n",
      "Crayford 5\n",
      "North End 4\n",
      "Erith 5\n",
      "Belvedere 7\n",
      "Thamesmead East 1\n",
      "Falconwood and Welling 2\n",
      "East Wickham 2\n",
      "Blendon and Penhill 2\n",
      "Danson Park 1\n",
      "Christchurch 5\n",
      "St. Michael's 1\n",
      "Brampton 1\n",
      "Northumberland Heath 6\n",
      "Barnehurst 5\n",
      "Colyers 1\n",
      "Lesnes Abbey 1\n",
      "Chase 8\n",
      "Winchmore Hill 5\n",
      "Cockfosters 1\n",
      "Highlands 2\n",
      "Upper Edmonton 2\n",
      "Palmers Green 1\n",
      "Edmonton Green 3\n",
      "Lower Edmonton 1\n",
      "Grange 2\n",
      "Jubilee 1\n",
      "Ponders End 4\n",
      "Enfield Highway 3\n",
      "Bowes 1\n",
      "Southgate Green 2\n",
      "Southgate 5\n",
      "Haselbury 1\n",
      "Bush Hill Park 3\n",
      "Town 11\n",
      "Southbury 2\n",
      "Turkey Street 2\n",
      "Enfield Lock 3\n",
      "Leyton 5\n",
      "High Street 5\n",
      "Higham Hill 3\n",
      "Valley 5\n",
      "Chingford Green 7\n",
      "Cathall 4\n",
      "Lea Bridge 2\n",
      "Markhouse 1\n",
      "Grove Green 2\n",
      "Forest 1\n",
      "William Morris 1\n",
      "Hoe Street 6\n",
      "Wood Street 5\n",
      "Chapel End 3\n",
      "Hale End and Highams Park 2\n",
      "Cann Hall 4\n",
      "Larkswood 3\n",
      "Endlebury 1\n",
      "Hatch Lane 1\n",
      "Leytonstone 6\n",
      "Wanstead 1\n",
      "Cranbrook 0\n",
      "Newbury 2\n",
      "Roding 4\n",
      "Fairlop 2\n",
      "Goodmayes 0\n",
      "Aldborough 0\n",
      "Bridge 1\n",
      "Hainault 0\n",
      "Loxford 0\n",
      "Clementswood 3\n",
      "Mayfield 0\n",
      "Valentines 2\n",
      "Seven Kings 2\n",
      "Snaresbrook 5\n",
      "Clayhall 1\n",
      "Church End 4\n",
      "Barkingside 0\n",
      "Fullwell 0\n",
      "Chadwell 0\n",
      "Monkhams 3\n",
      "Carshalton South and Clockhouse 1\n",
      "Cheam 5\n",
      "Beddington South 0\n",
      "St. Helier 0\n",
      "Wandle Valley 2\n",
      "Belmont 1\n",
      "Nonsuch 2\n",
      "Worcester Park 2\n",
      "Sutton South 0\n",
      "Sutton West 1\n",
      "Sutton Central 16\n",
      "Carshalton Central 9\n",
      "Sutton North 1\n",
      "Stonecot 2\n",
      "The Wrythe 1\n",
      "Wallington South 3\n",
      "Wallington North 4\n",
      "Beddington North 2\n",
      "Hampton 5\n",
      "Teddington 10\n",
      "Hampton Wick 4\n",
      "Heathfield 2\n",
      "Twickenham Riverside 16\n",
      "Ham, Petersham and Richmond Riverside 8\n",
      "North Richmond 3\n",
      "Kew 4\n",
      "East Sheen 4\n",
      "Mortlake and Barnes Common 10\n",
      "Hampton North 2\n",
      "Fulwell and Hampton Hill 3\n",
      "West Twickenham 8\n",
      "Whitton 3\n",
      "South Twickenham 1\n",
      "St. Margarets and North Twickenham 3\n",
      "South Richmond 19\n",
      "Barnes 5\n",
      "Cannon Hill 0\n",
      "Cricket Green 6\n",
      "Village 7\n",
      "Wimbledon Park 1\n",
      "Lower Morden 1\n",
      "West Barnes 1\n",
      "Raynes Park 2\n",
      "Dundonald 4\n",
      "St. Helier 2\n",
      "Ravensbury 1\n",
      "Merton Park 1\n",
      "Abbey 8\n",
      "Figge's Marsh 3\n",
      "Pollards Hill 1\n",
      "Graveney 3\n",
      "Longthornton 0\n",
      "Hillside 3\n",
      "Trinity 3\n",
      "Lavender Fields 0\n",
      "Colliers Wood 5\n",
      "Roehampton and Putney Heath 4\n",
      "Thamesfield 15\n",
      "Wandsworth Common 1\n",
      "West Hill 0\n",
      "West Putney 2\n",
      "East Putney 0\n",
      "Graveney 2\n",
      "Tooting 3\n",
      "Southfields 4\n",
      "Earlsfield 9\n",
      "Nightingale 8\n",
      "Furzedown 1\n",
      "Bedford 1\n",
      "Balham 6\n",
      "Fairfield 18\n",
      "Northcote 11\n",
      "Latchmere 7\n",
      "St. Mary's Park 10\n",
      "Queenstown 9\n",
      "Shaftesbury 3\n",
      "Palace Riverside 3\n",
      "Munster 2\n",
      "Sands End 1\n",
      "Parsons Green and Walham 11\n",
      "Fulham Broadway 13\n",
      "Ravenscourt Park 14\n",
      "Hammersmith Broadway 15\n",
      "Avonmore and Brook Green 10\n",
      "Askew 5\n",
      "Wormholt and White City 2\n",
      "Shepherd's Bush Green 6\n",
      "College Park and Old Oak 1\n",
      "Town 4\n",
      "Fulham Reach 3\n",
      "North End 9\n",
      "Addison 7\n",
      "Chelsea Riverside 4\n",
      "Holland 3\n",
      "Norland 5\n",
      "Redcliffe 8\n",
      "Stanley 8\n",
      "Royal Hospital 8\n",
      "Notting Dale 5\n",
      "St. Helen's 0\n",
      "Colville 7\n",
      "Dalgarno 0\n",
      "Golborne 5\n",
      "Earl's Court 6\n",
      "Courtfield 5\n",
      "Abingdon 7\n",
      "Queen's Gate 5\n",
      "Brompton & Hans Town 10\n",
      "Campden 9\n",
      "Pembridge 4\n",
      "St. James's 111\n",
      "Lancaster Gate 10\n",
      "Knightsbridge and Belgravia 18\n",
      "Hyde Park 19\n",
      "Churchill 7\n",
      "Vincent Square 7\n",
      "West End 105\n",
      "Harrow Road 4\n",
      "Bryanston and Dorset Square 23\n",
      "Regent's Park 2\n",
      "Abbey Road 5\n",
      "Marylebone High Street 33\n",
      "Queen's Park 1\n",
      "Bayswater 9\n",
      "Warwick 14\n",
      "Tachbrook 4\n",
      "Westbourne 7\n",
      "Little Venice 6\n",
      "Maida Vale 1\n",
      "Church Street 9\n",
      "Swiss Cottage 1\n",
      "Highgate 10\n",
      "Kilburn 6\n",
      "Bloomsbury 28\n",
      "Holborn and Covent Garden 56\n",
      "St. Pancras and Somers Town 11\n",
      "Camden Town with Primrose Hill 32\n",
      "Cantelowes 6\n",
      "King's Cross 21\n",
      "Fortune Green 3\n",
      "Frognal and Fitzjohns 1\n",
      "Gospel Oak 5\n",
      "Hampstead Town 12\n",
      "Kentish Town 14\n",
      "West Hampstead 5\n",
      "Belsize 3\n",
      "Haverstock 5\n",
      "Regent's Park 24\n",
      "St. Katharine's & Wapping 9\n",
      "Shadwell 4\n",
      "Island Gardens 4\n",
      "Canary Wharf 13\n",
      "Blackwall & Cubitt Town 2\n",
      "Poplar 4\n",
      "Spitalfields & Banglatown 19\n",
      "St. Peter's 12\n",
      "Bethnal Green 10\n",
      "St. Dunstan's 6\n",
      "Mile End 2\n",
      "Lansbury 5\n",
      "Bromley South 3\n",
      "Bow West 9\n",
      "Bow East 4\n",
      "Bromley North 2\n",
      "Whitechapel 11\n",
      "Stepney Green 3\n",
      "Limehouse 3\n",
      "Weavers 9\n",
      "Clerkenwell 31\n",
      "Caledonian 16\n",
      "St. Mary's 28\n",
      "Holloway 8\n",
      "Highbury East 6\n",
      "Highbury West 9\n",
      "Tollington 13\n",
      "Bunhill 26\n",
      "St. Peter's 19\n",
      "Canonbury 9\n",
      "Barnsbury 19\n",
      "St. George's 2\n",
      "Junction 12\n",
      "Finsbury Park 7\n",
      "Hillrise 1\n",
      "Mildmay 8\n",
      "Hoxton East & Shoreditch 27\n",
      "Haggerston 4\n",
      "De Beauvoir 5\n",
      "London Fields 9\n",
      "Hackney Wick 5\n",
      "Clissold 6\n",
      "Hackney Downs 5\n",
      "Lea Bridge 4\n",
      "Cazenove 1\n",
      "Springfield 3\n",
      "King's Park 2\n",
      "Brownswood 1\n",
      "Hoxton West 14\n",
      "Dalston 5\n",
      "Victoria 6\n",
      "Hackney Central 7\n",
      "Homerton 5\n",
      "Shacklewell 1\n",
      "Stoke Newington 10\n",
      "Woodberry Down 2\n",
      "Stamford Hill West 0\n",
      "Highgate 7\n",
      "Harringay 4\n",
      "Fortis Green 0\n",
      "Bounds Green 3\n",
      "Crouch End 7\n",
      "Stroud Green 2\n",
      "Muswell Hill 5\n",
      "Hornsey 5\n",
      "Alexandra 2\n",
      "Noel Park 5\n",
      "Woodside 4\n",
      "St. Ann's 2\n",
      "Seven Sisters 1\n",
      "Tottenham Green 2\n",
      "West Green 3\n",
      "Tottenham Hale 6\n",
      "White Hart Lane 0\n",
      "Bruce Grove 4\n",
      "Northumberland Park 8\n",
      "Royal Docks 2\n",
      "Canning Town North 5\n",
      "Beckton 1\n",
      "East Ham South 1\n",
      "Stratford and New Town 17\n",
      "Canning Town South 6\n",
      "Custom House 1\n",
      "Plaistow South 0\n",
      "West Ham 5\n",
      "Plaistow North 5\n",
      "Forest Gate South 3\n",
      "Boleyn 3\n",
      "East Ham Central 3\n",
      "Green Street West 0\n",
      "Green Street East 2\n",
      "East Ham North 2\n",
      "Wall End 2\n",
      "Forest Gate North 1\n",
      "Manor Park 3\n",
      "Little Ilford 0\n",
      "Gascoigne 1\n",
      "Thames 3\n",
      "River 2\n",
      "Village 2\n",
      "Abbey 5\n",
      "Longbridge 0\n",
      "Eastbury 1\n",
      "Goresbrook 1\n",
      "Mayesbrook 1\n",
      "Becontree 0\n",
      "Alibon 1\n",
      "Valence 0\n",
      "Heath 2\n",
      "Whalebone 1\n",
      "Eastbrook 1\n",
      "Chadwell Heath 3\n",
      "Parsloes 1\n",
      "City of London 159\n"
     ]
    }
   ],
   "source": [
    "for ward in wards:\n",
    "    ward_shape = asShape(ward['geometry'])\n",
    "    count = 0\n",
    "    for pub_id in tree.intersection(ward_shape.bounds):\n",
    "        pub = pubs[pub_id]\n",
    "        if ward_shape.contains(pub['point']):\n",
    "            count += 1\n",
    "    ward['properties']['pubcount'] = count\n",
    "    ward['properties']['pubcount_per_hectare'] = count / ward['properties']['HECTARES']\n",
    "    print ward['properties']['NAME'], count"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "feature_collection = {\"type\": \"FeatureCollection\",\n",
    "                      \"features\": wards\n",
    "                      }\n",
    "with open('data/wards_with_pubs.json', 'w') as f:\n",
    "    json.dump(feature_collection, f)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}