Shapely Experiments.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [],
   "source": [
    "import geopandas as gpd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [],
   "source": [
    "from shapely.geometry import Polygon"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [],
   "source": [
    "polygon1 = Polygon([(0, 0), (1, 1), (1, 0)])\n",
    "polygon2 = Polygon([(0.5, 0.5), (1.5, 1.5), (1.5, 0)])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {},
   "outputs": [],
   "source": [
    "gs = gpd.GeoSeries()\n",
    "gs['poly1'] = polygon1\n",
    "gs['poly2'] = polygon2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x112ec1cc0>"
      ]
     },
     "execution_count": 44,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQIAAAD8CAYAAACcoKqNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAELJJREFUeJzt3X2MZXV9x/H3R1DoA1V0MEVwmCVdYpeCAQZLkdSlYli2CdRoG1BasAtbWyVNME1r1tJVbAua1MbUhxJKqCZdRFrJVnbDgyzSlu7KUB4WaIFlQQuY7ghIQxFk8ds/7hn27N25c8/M/Z17nj6vZMKde8/c882Z4b13zt49P0UEZtZtr6l6ADOrnkNgZg6BmTkEZoZDYGY4BGaGQ2BmOARmhkNgZsD+Ve14YmIipqamqtq9WSfcddddP4iIQ4ZtV1kIpqammJmZqWr3Zp0g6btFtvOvBmbmEJiZQ2BmOARmhkNgZjgEZoZDYGY4BGaGQ2DWOA9v2pT8OR0CswbZsn49/3rZZcmf1yEwa4gt69fz7U9+spTndgjMGqDMCIBDYFZ7ZUcAHAKzWhtHBMAhMKutcUUAHAKzWhpnBMAhMKudcUcAHAKzWqkiAlAgBJKukrRL0v1DtjtR0m5J7083nll3VBUBKPaK4Gpg1UIbSNoPuBy4KcFMZp1TZQSgQAgi4nbgmSGbXQT8I7ArxVBmXVJ1BCDBOQJJhwHvBb5UYNu1kmYkzczOzo66a7PGq0MEIM3Jwr8G/jgifjJsw4i4IiKmI2L6kEOGXmrdrNXqEgFIs67BNHCNJIAJYLWk3RFxfYLnNmulOkUAEoQgIpbN3ZZ0NfBNR8BssLpFAAqEQNIGYCUwIekJ4M+A1wJExJdLnc6sZeoYASgQgog4p+iTRcT5I01j1mJ1jQD4nYVmY1HnCIBDYFa6ukcAHAKzUjUhAuAQmJWmKREAh8CsFE2KADgEZsk1LQLgEJgl1cQIgENglkxTIwAOgVkSTY4AOARmI2t6BMAhMBtJGyIADoHZkrUlAuAQmC1JmyIADoHZorUtAuAQmC1KGyMADoFZYW2NADgEZoW0OQKQYKUjSR+UdJ+k7ZLukPT29GOaVaftEYA0Kx09BrwrIo4BLgWuSDCXWS10IQJQ7JqFt0uaWuDxO3KfbgUOH30ss+p1JQKQ/hzBGmBz4uc0G7suRQDSLHACgKRT6YXglAW2WQusBZicnEy1a7OkuhYBSPSKQNKxwJXAWRHx9KDtvOSZ1V0XIwBpFkGdBP4J+O2IeHj0kcyq0dUIQJqVji4B3gR8MVv/cHdETJc1sFkZuhwBSLDSUURcAFyQbCKzMet6BMDvLLSOcwR6HALrLEdgD4fAOskR2JtDYJ3jCOzLIbBOcQTm5xBYZzgCgzkE1gmOwMIcAms9R2A4h8BazREoxiGw1nIEinMIrJUcgcVxCKx1HIHFcwisVRyBpXEIrDUcgaVzCKwVHIHROATWeI7A6BwCazRHII0UKx1J0ucl7chWPDo+/Zhm+3IE0kmx0tEZwPLsYy3wpdHHMluYI5DW0BBExO3AMwtschbwlejZCrxB0qGpBjTr5wikl+IcwWHAf+c+fyK7zyw5R6AcYz1ZKGmtpBlJM7Ozs+PctbXAlks+4QiUJEUIngTemvv88Oy+fXilI1uqLX+6jm9f+udVj9FaKUKwEfid7G8PTgKei4jvJ3heMyCLwKf/ouoxWi3FSkebgNXADuAF4ENlDWvd4wiMR4qVjgL4SLKJzDKOwPj4nYVWS47AeDkEVjuOwPg5BFYrjkA1HAKrDUegOg6B1YIjUC2HwCrnCFTPIbBKOQL14BBYZRyB+nAIrBKOQL04BDZ2jkD9OAQ2Vo5APTkENjaOQH05BDYWjkC9OQRWOkeg/hwCK5Uj0AwOgZXGEWgOh8BK4Qg0i0NgyTkCzVMoBJJWSXooW9bsT+Z5fFLSFkl3Z8uerU4/qjWBI9BMRdY+3A/4Ar2lzVYA50ha0bfZJ4BrI+I44Gzgi6kHtfpzBJqryCuCdwA7ImJnRPwYuIbeMmd5Afxcdvv1wFPpRrQmcASarUgIiixpth44N7vc+SbgovmeyCsdtZMj0HypThaeA1wdEYfTW+Pgq5L2eW6vdNQ+jkA7FAlBkSXN1gDXAkTEvwMHAhMpBrT6cgTao0gI7gSWS1om6XX0TgZu7Nvme8C7AST9Ir0Q+LV/izkC7TI0BBGxG/gocCPwn/T+duABSZ+SdGa22ceACyXdC2wAzs9WQLIWcgTaZ+iSZwARsYneScD8fZfkbj8IvDPtaFZHjkA7+Z2FVpgj0F4OgRXiCLSbQ2BDOQLt5xDYghyBbnAIbCBHoDscApuXI9AtDoHtwxHoHofA9uIIdJNDYK9yBLrLITDAEeg6h8AcAXMIus4RMHAIOs0RsDkOQUc5ApbnEHSQI2D9HIKOcQRsPg5BhzgCNkiSlY6ybX5L0oOSHpD0D2nHtFE5AraQoZcqy6109B56axrcKWljdnmyuW2WAx8H3hkRz0p6c1kD2+I5AjZMqpWOLgS+EBHPAkTErrRj2lI5AlZEqpWOjgKOkvRvkrZKWpVqQFs6R8CKKnQV44LPsxxYSW8BlNslHRMRP8xvJGktsBZgcnIy0a5tPo6ALUaqlY6eADZGxMsR8RjwML0w7MVLno2HI2CLlWqlo+vpvRpA0gS9XxV2JpzTCnIEbClSrXR0I/C0pAeBLcAfRcTTZQ1t83MEbKlSrXQUwMXZh1XAEbBR+J2FLeAI2KgcgoZzBCwFh6DBHAFLxSFoKEfAUnIIGsgRsNQcgoZxBKwMDkGDOAJWFoegIRwBK5ND0ADx4//jsBOO4S0nHFf1KNZSDkHdvfwC+smLHLXqPVz4LzfxgW9scBAsOYegzl5+AV750aufSnIQrBQOQV31RSDPQbDUHII6WiACeQ6CpeIQ1E3BCOQ5CDYqh6BOlhCBPAfBlsohqIsRI5DnINhiOQR1kDACeQ6CFeUQVK2kCOQ5CDZMsiXPsu3eJykkTacbscXGEIE8B8EGGRqC3JJnZwArgHMkrZhnu4OAPwS2pR6ylcYcgTwHwfqlWvIM4FLgcuDFhPO1U4URyHMQbE6SJc8kHQ+8NSJuWOiJJK2VNCNpZnZ2dtHDtkJNIpDnINjIJwslvQb4K+Bjw7bt/EpHNYxAnoPQXSmWPDsI+CXgNkmPAycBG33CsE/NI5DnIHTPyEueRcRzETEREVMRMQVsBc6MiJlSJm6iBkUgz0HojlRLntkgDY1AnoPQfkmWPOu7f+XoY7VECyKQNxeE5aefxiM33sJtn/4MT911d9VjWQJ+Z2FZWhaBPL9CaB+HoAwtjkCeg9AeDkFqHYlAnoPQfA5BSh2MQJ6D0FwOQSodj0Ceg9A8DkEKjsC8HITmcAhG5QgM5SDUn0MwCkdgURyE+nIIlsoRWDIHoX4cgqVwBJJwEOrDIVgsRyA5B6F6DsFiOAKlchCq4xAU5QiMjYMwfg5BEY5AJfYJwvTxVY/UWg7BMI5A5fYE4WY+cMMNvOXEE6seqXUcgoU4ArUiiaNWr+bCbdschMQcgkEcgdpyENJLstKRpIslPSjpPknfknRE+lHHyBFoBAchnVQrHd0NTEfEscB1wGdSDzo2jkDjOAijS7LSUURsiYgXsk+30rvkefM4Ao3mICxdkpWO+qwBNo8yVCUcgdZwEBYv6clCSecC08BnBzxezyXPHIFWchCKS7HSEQCSTgPW0Vvc5KX5nqiWS545Aq3nIAw38kpHAJKOA/6WXgR2pR+zJI5ApzgIg6Va6eizwM8CX5d0j6SNA56uPhyBznIQ9pVkpaOIOC3xXOVyBIw9QVh+xhk8snkzt61fz1N33ln1WJXo3jsLHQHr41cIXQuBI2AL6HIQuhMCR8AK6mIQuhECR8CWoEtBaH8IHAEbUReC0O4QOAKWUJuD0N4QOAJWkjYGoZ0hcARsDNoUhPaFwBGwMWtDENoVAkfAKtTkILQnBI6A1UQTg9COEDgCVkNNCkLzQ+AIWM01IQjNDoEjYA1S5yA0NwSOgDVUHYPQzBA4AtYCdQpC80LgCFjL1CEIqVY6OkDS17LHt0maSj0o4AhYq1UZhFQrHa0Bno2IXwA+B1yeelBHwLqiiiAkWeko+/zvs9vXAe+WpGRTOgLWQeMMQpGLl8630tEvD9omInZLeg54E/CDkSfc/SN45UUgXVesqbr5M9B/kdWdt9ySfB+FrmKciqS1wFqAycnJYl+0/0/1Psw6bi4IR61enfy5U6109Oo2kvYHXg883f9EtVzpyMzSrHSUfX5edvv9wK0REenGNLMyDf3VIPudf26lo/2Aq+ZWOgJmImIj8HfAVyXtAJ6hFwsza4hUKx29CPxm2tHMbFya985CM0vOITAzh8DMHAIzwyEwM0BV/XW/pFnguwU2nSDFW5VH4xk8Q1NnOCIihr57r7IQFCVpJiKmPYNn8AzlzeBfDczMITCzZoTgiqoHwDPM8Qw9rZuh9ucIzKx8TXhFYGYlqzQEo1wUVdLHs/sfknR6iTNcLOlBSfdJ+pakI3KPvSLpnuyj/59mp5zhfEmzuX1dkHvsPEmPZB/n9X9tov1/LrfvhyX9MPdYqmNwlaRdku4f8LgkfT6b8T5Jx+ceG/kYFJzhg9m+t0u6Q9Lbc489nt1/j6SZEmdYKem53DG/JPfYgt/HBUVEJR/0/knzo8CRwOuAe4EVfdv8AfDl7PbZwNey2yuy7Q8AlmXPs19JM5wK/HR2+/fnZsg+f35Mx+F84G/m+do3Ajuz/x6c3T449f77tr+I3j9FT3YMsuf5VeB44P4Bj68GNtO7XtlJwLZUx2ARM5w899z0Lua7LffY48DEGI7DSuCbo34f+z+qfEUwykVRzwKuiYiXIuIxYEf2fMlniIgtEfFC9ulWeldoSqnIcRjkdODmiHgmIp4FbgZWlbz/c4ANi9zHUBFxO71rWQxyFvCV6NkKvEHSoaQ5BoVmiIg7sn1AOT8LRY7DIKP8HFUagvkuinrYoG0iYjcwd1HUIl+baoa8NfT+VJpzoKQZSVsl/cYS9r+YGd6XvSy9TtLcpeNSHIfCz5H9WrQMuDV3d4pjUMSgOVP9LCxW/89CADdJuiu7NmeZfkXSvZI2Szo6u2+k4zDWi5c2maRzgWngXbm7j4iIJyUdCdwqaXtEPFrC7v8Z2BARL0n6PXqvkn6thP0MczZwXUS8krtvXMegNiSdSi8Ep+TuPiU7Dm8Gbpb0X9mf7qn9B71j/ryk1cD1wPJRn7TKVwSjXBS1yNemmgFJpwHrgDMj4qW5+yPiyey/O4HbgOPKmCEins7t90rghMXMP+r+c86m79eCRMegiEFzpvpZKETSsfS+B2dFxKsX6M0dh13AN1jar6pDRcT/RsTz2e1NwGslTTDqcRj15MYIJ0X2p3diZxl7Tm4c3bfNR9j7ZOG12e2j2ftk4U6WdrKwyAzH0TsJs7zv/oOBA7LbE8AjLOLkzCJnODR3+73A1thzouyxbJaDs9tvTL3/bLu30TshptTHIPd8Uww+Sfbr7H2y8DupjsEiZpikdz7q5L77fwY4KHf7DmBVSTP8PHve//MO4HvZMSn0fRy4z6UOm+KD3pngh7P/0dZl932K3p+8AAcCX88O/neAI3Nfuy77uoeAM0qc4Rbgf4B7so+N2f0nA9uzA74dWFPiDH8JPJDtawvwttzX/m52fHYAHypj/9nn64HL+r4u5THYAHwfeJne77drgA8DH84eF72l9x7N9jWd8hgUnOFK4Nncz8JMdv+R2TG4N/s+rStxho/mfha2kovSfN/Hoh9+Z6GZ+Z2FZuYQmBkOgZnhEJgZDoGZ4RCYGQ6BmeEQmBnw/ya4bx0WB/9SAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "gs.plot(cmap='OrRd')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"100.0\" height=\"100.0\" viewBox=\"0.47 0.22 0.56 0.81\" preserveAspectRatio=\"xMinYMin meet\"><g transform=\"matrix(1,0,0,-1,0,1.25)\"><path fill-rule=\"evenodd\" fill=\"#66cc99\" stroke=\"#555555\" stroke-width=\"0.016200000000000003\" opacity=\"0.6\" d=\"M 0.5,0.5 L 1.0,1.0 L 1.0,0.25 L 0.5,0.5 z\" /></g></svg>"
      ],
      "text/plain": [
       "<shapely.geometry.polygon.Polygon at 0x112f57cc0>"
      ]
     },
     "execution_count": 45,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "polygon1.intersection(polygon2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "True"
      ]
     },
     "execution_count": 46,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "polygon1.intersects(polygon2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "False"
      ]
     },
     "execution_count": 47,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "polygon1.within(polygon2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"100.0\" height=\"100.0\" viewBox=\"-0.06 -0.06 1.62 1.62\" preserveAspectRatio=\"xMinYMin meet\"><g transform=\"matrix(1,0,0,-1,0,1.5)\"><path fill-rule=\"evenodd\" fill=\"#66cc99\" stroke=\"#555555\" stroke-width=\"0.032400000000000005\" opacity=\"0.6\" d=\"M 0.0,0.0 L 0.5,0.5 L 1.0,1.0 L 1.5,1.5 L 1.5,0.0 L 1.0,0.25 L 1.0,0.0 L 0.0,0.0 z\" /></g></svg>"
      ],
      "text/plain": [
       "<shapely.geometry.polygon.Polygon at 0x112f57c50>"
      ]
     },
     "execution_count": 48,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "polygon1.union(polygon2)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}