{ "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.9" }, "name": "", "signature": "sha256:c5056cab2e638465802653766bf6a7c98bbb86edd67a12e0fa417523c382e697" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "import os, sys\n", "from glob import glob\n", "import networkx as nx\n", "import numpy as np\n", "import pickle\n", "import matplotlib.pyplot as plt\n", "sys.path.append('../');\n", "import Holes as ho" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "import geojson as gj" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "%matplotlib inline" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "max_dimensions=1;" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "import scipy.io\n", "ld = os.listdir('../stores/phom/country_high_act_low_entropy_matrices/');\n", "countries_to_cells={}\n", "for country in ld:\n", " if country == '.gitignore':\n", " continue\n", " countries_to_cells[country[:-4]]=pickle.load(open('../stores/phom/country_high_act_low_entropy_matrices/'+country))\n", " " ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 5 }, { "cell_type": "code", "collapsed": false, "input": [ "grid=gj.load(open('../data/poi/milano-grid.geojson'))" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [ "import pandas as pd\n", "\n", "cell_position = {}\n", "for i,cell in enumerate(grid['features']):\n", " cell_position[cell[\"properties\"][\"cellId\"]] = cell[\"geometry\"][\"coordinates\"][0][0]" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 7 }, { "cell_type": "code", "collapsed": false, "input": [ "import networkx as nx\n", "country_graphs={}\n", "for country in countries_to_cells:\n", " coords = np.nonzero(countries_to_cells[country]);\n", " country_graphs[country]=nx.Graph();\n", " coords = zip(coords[0],coords[1]);\n", " for i, cell in enumerate(coords):\n", " for j, other_cell in enumerate(coords):\n", " if i