shapely_demo.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from shapely.geometry import Point, Polygon, LineString\n",
    "import descartes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "circle = Point(5.0, 0.0).buffer(10.0)\n",
    "clip_poly = Polygon([[-9.5, -2], [2, 2], [3, 4], [-1, 3]])\n",
    "clipped_shape = circle.difference(clip_poly)\n",
    "\n",
    "line = LineString([[-10, -5], [15, 5]])\n",
    "line2 = LineString([[-10, -5], [-5, 0], [2, 3]])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Blue line intersects clipped shape: True\n",
      "Green line intersects clipped shape: False\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEACAYAAACj0I2EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VOXd//H3TPaFfUnIRoIKosiiuG8BEde6tdW6VK3W\n9le11q5W7aN59LGrPlfdtVrrjrVUpdYFWUTEigoIsskiBEIICUtYQhJCkvn9cZLHIcxkPWfuM+d8\nXtd1LrIMc77DzHz4zn3u+xwQEREREREREREREREREREREREREYkbzwCVwNKwn/UHZgCrgfeAvgbq\nEhHxlaAN9/E34Ow2P/s1VqAPB2a1fC8iInGgkAM79C+BrJavs1u+FxERB9nRoUeShTUMQ8ufWe3c\nVkREbOBUoIcLtWwiIuKgRIfutxJrqGULMASoinCbxcAYh/YvIuJVS4CxkX7hVKD/C7gG+EPLn29E\nuM0YIODQ/iMpadm8qgQ9PicEgP6QUAS5R0L/wyAxCwYFoSAERUHIa4K8Osivg4J6SOzGJ9JLi+HV\nOV37OzsToTQNytJgcyqUhmBDCMqCsLcG9pXBxqVQsxYoA/Z1vS7blODd12cJsX1sUV9fdgT6FOB0\nYCDWi+Yu4PfAq8D1QClwqQ37EYmF1gAvhLwjod8YyO4H44FxDTBiDxRVQHqz4TqBvo0wdo+1hWsG\ntibDmjxYdjgsAJaHYM9a2Pg57PkK2Ag0xL5mcZIdgX55lJ9PsuG+RWKhpQOPFOBH74SCjbE53GSX\nIJDVYG2nVFs/qw3CkmxYcqkC3rucGnJxozmmC3DYHNMFOGyOzffXG5JHw2ETIbcIjguZDfBJpc7e\nf3oznLjT2uDggF+6Hyo+grL5wDqgyeYC5th8f24yx3QBrWI5ht1WyPD+xX9SIfEIOOQ0yB5tfYg8\nYzccvzO+OnAnbEuCmYNgehIs2Q1l78O2BUA5mqXmNlGzU4EuXpcIHApDT4TsE+HUJJi0F07fAaku\nGAd3o3VpMGMgvBuErzbDutmwdzGw3XRlAijQxYeGQNbJkHc6jMuAyfvhjK3Qv9F0YfGjGfiiF8zs\nB+8FoHwVrJoNTYvReLtJCnTxhQAwFIadByOOhUsa4aytkG9yup5HNAbgo37wdm+YuQfWvgG7PwZq\nTVfmQwp08bQAMByGXwCjjoTL98H5lRpSccraNJgyGP5dD1+9BdvnArtNV+UjCnTxpCAkjoLhF8O4\nIrhiL0ze1r2FPdJ15SnwSha83gRfzYAts9A4eywo0MVTkiBlHAy/BE7Mhit2w6k7NFPFlG1JMDUb\nXg3AV3Nh43SgwnRVHqZAF08IQsp4GHk5TOgHV+yA8fqo7xo1CfBaFkxJglWfwPqpwDbTVXmQAl3i\nXj4MvwqKR8KPq2BUjemCJJr6IPwjG54OwtJ/QvUsNCvGTgp0iVsZkH0+jDkHbqqF86o0tBIvKpLh\nsRyYVgUrn4PGFWiRkh0U6BJ3gpByLBx1FXwnA24oh952L0eXmJjfBx7qDx9/CqV/R8MwPaVAl7jS\nMrwy4XD4SRWM3Gu6IOmpxgBMGdIyDPMaVM9EwzDdpUCXuJABOefDURpe8azKlmGYN7bCimc1DNMt\nCnRxu6TDYNSNcGUfDa/4wfw+8OAAmDcbNr0C1JuuKI4o0MW1EqDfWXDCpXBXNZywy3RBEiv1QXgw\nD16ogOWPYZ3ZUTqmQBdX6gvDrodvjIbfbIKB+00XJCZMHwi/S4VFz8Kej9AQTEcU6OI2SSNh9I/g\n5lS4erPGyv2uNBXuzYEZH0LZS+ikX+1RoItrJMLA8+Gki+Gu7XCMVnpKi4YAPJYHf90Oyx4FNpiu\nyKUU6OIKA+DQG+DikXBHmXWRY5G23u8P92XCZy/A7g+wTswuX1Ogi2lJI2HczXBrAly2RUMs0r7y\nFLgnF6Z/AhueAepMV+QiCnQxKf0EOO2H8LsdMHaP6WokXjQG4M/58NQ6WP0QOud6KwW6GBGAvmfC\nGVfB/ZuhUHONpRv+lgsPbIflDwBbTVfjAgp0ibkgZF0CZ18AfyyDwZqSKD3wxmC4twEW3Q9sMl2N\nYQp0iakkyP0uXFwM923Qqk+xx/v94c5kWPAA7F9juhqDFOgSM6lQ+EO4chzctQGStUhEbLSwN9zW\nFz56COoXm67GEAW6xEQvGP4TuGEY/GyjZrKIM1alwy+yYM5foWau6WoMUKCL4wbCkT+Hnw6A72nl\npzisPAV+mQvvvQrb38JfpwtQoIuj+sJRd8DdveCblaaLEb/YkQi/HArTXmkJdb+Imp1qo6SnMmDE\nT+FXfRTmElv9G+FPG2DSZZB5qulq3ECBLj2RDIfcCDflwVUVposRP+rfCA+UQ/H3IXWs6WpMU6BL\ndyVA/vVw9ZFwU5npYsTPcvfB/ZVw0i3WhVL8S4Eu3RGA7Mvg2yfCHRv0MhLzRtTCH3bC+J8Deaar\nMUXvROmGAefBN86Ge0sh0U+zC8TVxu+G/9kPR/8CGGi6GhMU6NJFmafC5Evh9xshXac1FZeZuB1+\nkwKjfg70Nl1NrCnQpQtSx8KE6+FP5dbBKBE3urgSfjYQht8CpJmuJpYU6NJZeXDCzfBAlXUQSsTN\nrim3ViwPvQYfrXdRoEtnpMGom+DOOjhM13qUOBAEbi2DySdBn9NMVxMrCnTpSADyL4ers2DSdtPF\niHReYgjuLofxVwP5pquJBQW6dCD9eJhYDD/2+zmoJS7l7oNf18Kom/HBeLoCXdozBMZfD3dVQKpm\ntEicmrQdrh0EeVfi8fF0BbpEkwIjb4TbGmCYLtArce7HZTDpNEg/0XQlTlKgSyQBGHIpXJUL5+oa\njuIBySG4ezMc+z0gx3Q1TlGgSwQpR0PxmXCrxs3FQwrr4dcNcMRNQKrpapygQJe2BsP4H8DdW7QS\nVLzn7K1wVQ7kXIYHx9MV6BIuAIdeA7eErJMdiXjRT8vglImQdLjpSuymQJcwKcfA6UfBt7aYrkTE\nOanNcGs1jLgWSDZdjZ2cDvRS4Avgc+BTh/clPZMOR34XbqnS//PifSfuhAuyoM8E05XYyel3bggo\nBsYBxzm8L+mRQefApb1gdI3pSkRi4+bNMPpbwADTldglFq2Y5w48eFAOjDkPbthsuhCR2BnSANcB\nhZfikZyKRYc+E1gA3ODwvqR7AnDolfCjBp0SV/znis1w3AldP0AayobQ7RD6CkJNELrJmfq6JtHh\n+z8ZqAAGATOAL4EPw35fEvb1nJZNYirlaDhtFFy03nQlIrGXHIJbqmHFtbDsN8D+6LcNBYFJwA+A\nCzkwP78NPOpQkcUtW4di+THjbqAGeKDl+1CM9y8HS4Nxv4dn92vsXPztjkJ4/GXY+d7BvwtlA9/D\nGmUoivCXq4FLITDTyQrDCyJKdjo55JIO9Gr5OgOYDCx1cH/SZYPOhW/3VpiL3FQBo8IOkIaCEJoM\noalAGfBbDg7zecB3gdwYhnm7nBxyyQJeD9vPS0CE//3EkMFw1Hnww3LThYiYl7sPrgtAv2vhX8m0\n340/BzwFgRWxrLAzTA55aMjFqNzvwj2nw3U6X4v4XQAYBs3HQOhwSIiUS/OAJ4F/QsD02UejZqfT\nB0XFnfrBIcVwSYXpQkQMygTGAscAfSOMQLu6G49Ege5LgybCt4C+mqYoftPSjXMMMIKIxxFLKyHp\nLsh9wQXdeJco0P0nEw6dDJfpfC3iJ2268YPUA4uBRfBmKty7HbbGVZiDAt2H+p0GFybD4Hbm24p4\nQie6cTYCC4EVQMsn1suSYMpk2Doda6p13FCg+0saHHI+XF5puhARB3WhGyfCFbkG77eani9Pg+q3\nHazTdgp0X8k4Ac5Jg4Iq05WI2Kyb3Xg0l1fC1PNhwftA3Ay9KND9IwmGXwRX6hqhndYQgC96wZI+\n1umIklPhwTWmq5ID9LAbj6agHs7NhpUnwN73bagzJhTovpFyDEzsAyNKTVfiXm0DfCmwewOUvw/7\nh8I9R5uuUADbu/ForqiCNy+Cz+fR7jle3EOB7g8BGH4hXL7DdCHu0l6A71gDbMD6uN0HTr4frtC8\nfbMc6sajGVELEwbDF6OhaWHP7895CnR/yIfh2TBug+lCzOp0gLcx6Az4VoJmBhkRo248mnP2wDvF\nsFKBLm6RdSyc3eS/S8t1N8AP0AeGn6PuPOZi3I1Hc9oOyD4SVvYFdjq3H3so0L0vAXKLYfI204U4\nz5YAb0PdeQwZ7sYjSQ7BpADMHwN1Hzi/v55RoHvfoTAuEwq2my7Efk4E+AHUnceGS7rxaCbthCkT\nYdlcrBNjuZYC3fOGngRnNZiuwh6OB3gb6s4d5MJuPJrxu6BwKCzLxroCm2sp0L0tFbJOhDPidCFR\nrAP8AOrOneHybjySIHB2COYdAzv/bbqa9ijQPS1xJJySFD8XfzYa4G2oO7dRHHXj0UzaBk9NgJ1v\nA82mq4lGge5pw06DM2tNVxGdwQAv4XDgh0AOcCclrA37rbpze8RhNx7NiFo4sgCWFALrTFcTjQLd\nu3pD9mgoduEl5kpT4YksmNcc0w68hCBwDnAL1jVuWwWAS7/+Vt15D3igG4/mrCZ4/zioUKBLrCUf\nBWcEIdVFHw9rg/BsDrywH5Y+D3vnE4shlBL6YF21/WbgkAi3mBv2tbrz7vFQNx7NmVvh4dOgYiou\n/Y9Ige5ZhUfDiXtNV2FpBt4dBI9lwOLZUD4N2OX4bq1hlZuBa4GMNr8NAf8CHqKE2V//WN15F3i4\nG49kSAMcngIL8oBS09VEokD3pgTodQQc44K55ysz4OHBMHsNrHoRp98I0YdVWu0EngYeo4T1bX6n\n7rxzfNCNR3NsAKYOg/pS05VEokD3phw4NMXs7JYdifB0DrxaA8seg32f4eTsgI6HVZYDDwEvUUKU\nTy7qztvhs248mjE1UDQWVs7u+Laxp0D3pOQiGG9o340BeCMLnkiGZW9B5TuAczNtOh5WmQY8DLxP\nSbur/NSdR+bjbjyScbsg83Cs7HTdf1oKdE8qGgtjDExXXNgbHhoIHy2Gr6bg1Kq6ng2rRKHuPIy6\n8ah6N8GIBPjMlePoCnTvaR0/j+HJuCqT4fEh8Pp2WHE/NC7FiXNe2DKsEpG6c4u68U45NgCvHwJ7\nS01X0pYC3Xty4JDk2IyfNwTg5Rx4Blj2D6ieBeyzfTf2DatE4evuXN14l42pgYKxsHKW6UraUqB7\nTqzGz+f2g4f7wmcfw4apgL2fCBwZVonIr925uvFuG7cLeo0AknDZpekU6J5TNBbGOjj/vDQVHsmG\nd8th9WOwfzV2Dq84N6wSha+6c3XjtmgdR/80D3rUTNhOge4tCZB5BBztwPh5+CrPFc/C7nnY+YZ3\nfFjlIKmQeAQcfq4PunN147YbH4TXhsFeBbo4JguKbD67ooOrPGM3rNIqCRgBRSdC1nFwcgJctMej\n3bm6cUcdtQcKjnLbOLoC3VsGWu9huzi0yjO2wyoJ1j5yj4chJ8FxqTBpH0yohL5eDDF14zFRVAup\neaaraEuB7impg6DQhitB/98qz722rvKM3bBKECiArPGQezqMzoTJjdaFPtSNo27cBnn7ICULlx0Y\nVaB7ypB8yK3v/t9vXeX5ZDIstWeVZ+yGVQJADgwYB3kTYGR/OLPJujCBF6+nCqgbNygxBEOagf5A\npelqWinQPaVXPuR3M4BtXuUZu2GVQdBrLBRNgKIhMLkZztgOIzb04D7dTN24awwFGIACXRwQgMQc\nGNrFbrR1lecbO2B5z1d5xmZYpS+kjYZhxZBfBGcCE3fA6A2R880T1I27TmEQEgZAk+lC/o8C3Tsy\noW8XzrBo4yrP2AyrZELCkTDidMg6HCYBE3bD8V4OcXXjrpbbALkF1lPgDgp07xgA+Z08cGnTKs/O\nDas8DLzYzWGVVOtC10WnQtYYOD0IZ+yBk8sg2f5zxbiHuvG4kFcLffJNVxFOge4dA6Aw0P5NbFrl\n2dkrAXVvWKXNXPGTEuGMWuvaqOkuupye7dSNx52iOkjKNV1FOAW6Z/TOgqFRAt2GVZ7ODqu0zBXP\nOw6yT7bmip+xDyZ6da54OHXjcWtQA2RkAuk4ec7/LlCge0Z2AeS0eVHZsMrTuWGVKHPFJ26FrIYu\n1Rh/1I17QhDIa526qEAXO6VlQVbYQc0ervJ0ZljFj3PFw6kb95wcsJ7LTYYLARToHhJIh8zGHq3y\ndG5YxW9zxcOpG/e0TIBk01W0UqB7RjAVPsiEVwKw9G2ofJvOfgx0ZljFj3PFw6kb94X0AAp0sd/O\nBvj9Clj7Mp1d5Wn/sEomJBwBI4p9NFc8nLpx38kMQkKKWxYXKdA9Y93vgGo6moZo/7BKy1zxYafA\n4LE+miseTt24b6U1Q1oa1JguBFCge8mOdn9r77BKEjAcik7y2VzxcOrGBUhtgox0BbrEhn3DKn6e\nKx5O3biESWuE9EzTVbRyMtDPBv6MFQRPA39wcF8Szr5hlQAw1KdzxcOpG5co0pogKc10Fa2cCvQE\n4BGso2LlwGdYneBKh/YnYNewit/niodTNy4dSGuCRM8H+nHAWr5ezPIKcCEKdGeUcChwKz0bVvHz\nXPFw6salCzKaICHddBWtnAr0XKAs7PtNwPEO7cvfSjgXeA1IafObzgyr+H2ueDh149IN6b4I9M5O\nVysJ+3pOyyadZYX56xy4sKGjYRW/zxUPp25ceigxhPNvnOKWrUNOBXo5EH6e4Hwin+ugxKH9e9/B\nYV4KfB+YHWFYRXPFD9SZbnwJVpCrG5d27E2AkNNzFudwYLN7d7QbOhXoC4DDgEJgM3AZcLlD+/Kf\nyGFeTAnh492aK34gdePigLoEaHbFmRbBuUBvxJppMR1rxstf0QFRe7Qf5porfjB14+KgvQmw3/OB\nDvBOyyZ2iR7mG2HIhTBkko/niodTNy4xUpcA++tMV9FKK0XjRfudeToccgm8VOHDueLh1I1LjNUn\nwD5fdOhil47HzBtgfzMU1JsozzB142JQbSLUdOcC6I5QoLtd5w6ANkJjEzQEfDRrRd24uMC+BNir\nDl06oXNh3iJUD7sTYeD+mNUXe+rGxWX2hmD/vo5vFxsKdLfqUpiDFei1CYAXA13duLhUTQhwzeQD\nBbobdTnMwQr0vW2X/8czdeMSB2pDgDp0iaJbYQ6wfydsz4dOX/PTrdSNSxzZHqKz1+6NAQW6m3Q7\nzAGqN8Km4ViXoYs36sYlTpUFANdMFVagu0WPwhygchOUJzlSm3PUjUscqw3C1mZgl+lKWinQ3aDH\nYQ7QtB3Wx8M5WtSNi0esT4PGLYBr3ncKdNNsCXMAtls56FrqxsVjytKhdq3pKsIp0E2yL8wBqqES\nly0uUjcuHlaeBhWuuqKXAt0Ue8McoAn2b4UNaXCY6aPu6sbFBzaEoKbKdBXhFOgm2B/mLeo2wYYR\nhgJd3bj4TGkIF81wAQV67DkW5gBVG6BibM/vp0vUjYtPbQqiQPcxR8McYMcW2BiLC4OqGxef25EI\nu+px2UI+BXqsOB7mAGyH9U4eEFU3LgLA+nRo3AwHXb/XKAV6LMQmzAEqYVUAGgMtVyO3g7pxkYN8\n2Quq5puuoi0FutNiF+YAe2F3GSzNhHF7enhf6sZFolqUAFtcd51kBbqTYhvmLTZ/DkvO72agqxsX\n6VAzsCgATr6Pu0eB7hQjYQ6wbTUsDMK1XflL6sZFOm1lBlRvBnabrqQtBboTjIU5ABtgMZ0YR1c3\nLtIti/tC1Xumq4hEgW43s2EOHY+jqxsX6ZFFQaj40nQVkSjQ7WQ+zFscNI6ublzEFs1Yge6+8XNQ\noNvHNWEO1jj6ggBcq25cxFbuHT8HBbo9XBXmoSCUHwI7zoBQPgTUjYvY5vN+sOVd01VEo0DvKdeE\neSgb+B5wA+QWQW7bG6gbF+mxRUGoXGW6imgU6D1hPMxDQWAS8EPgAiI/n+rGRWzREIBFYL3PXUmB\n3l1Gwzy8G6cowg2qoep1eKQflHwR+RioiHTNvP5QuQLo6Spsx+id3h1GwjwUhNBkCP0TKAN+y8Fh\nPg+4GsiFrO/DtLWwItO5mkT8ZEYmrPnAdBXtUYfeVTEP88504zwP/AUCKw781YZZMPMyGFXjTG0i\nfrE7AT5ogqblpitpjzr0rohZmHe1Gw/cenCYA+xaBNOxVo2KSPfNGgSVnwKmL+/YLnXonRWTMO9J\nNx7RVihbD58MgJN32leniN+8lwLrPjJdRUcU6J3haJh3aqbKPOAvwFQI1HXt/tfMhpnfV6CLdFdF\nMnxWB6w2XUlHFOgdcSzMbe/Go2j4AmY3wy+DkN7c/fsR8asZg6BiJrDfdCUdUaC3x/Ywd7obj2gP\nVCyBOSPgXC0oEumydxNh82emq+gMBXo0toZ5rLrxaNbPhZlj4Vz771rE01alw4oduHgxUTgFeiS2\nhLmRbjyKxpXwUSNsS4KBrv/YKOIeMwbCxn9inWbR9TRtsa0eh3koG0K3A2ux5gxewoFhXg08CBwJ\ngVMh8IKzYQ7APtg4C6ZmO7wfEQ+pSYDXgOpPTFfSWerQw3U7zN3UjUezZRa8ehZclQCZTbHfv0i8\neS0b1n4IbDNdSWcp0Ft1K8xNj413yTbrxfnaSXB1ueFaRFyuPgivJEHZdNOVdIUCHboY5vHQjUdT\nNh1eOR0uDUJqXIwJipjx7yxY9Rmw2XQlXaFA73SYx1U3Hs1m60X679HwrQrTxYi4U2MApqTAurdN\nV9JV/g70DsM8nrvxaNa9DVOOhYsCkBgyXY2I+7w3EJYtJ06mKobzb6C3G+ae6MajKbVerO8VaqGR\nSFvNwMsZsPpNIO4aHn9OW4wU5vt6TaAkNKLnZzh0vZD1Yn05I06m1orE0If9YdF6wLWXmWuP/zr0\ntmHenFDGM3NfZdNJs/FeNx7NKvh8PXw4AE7fYboYEfd4uTeseZo47M7BuQ69BNgEfN6yne3Qfrqm\nbZjXZNXy4NohbDrpV3izG48mBKtfhxf7qEsXafVpH/h4CzQuM11JdznVoYeA/23Z3KGEcwnxOoGW\nMK8uhGfnpLNraPitvNiNR9G4DOauhndy4bwq09WImNUYgEcGwMr7ieMux8kxdPdcJefKc35IU+Kb\nbcKcsDD3cjceTTOsfgEeTbcuryXiZ1OzYd5CaFxqupKecDLQfwwsAf4K9HVwP+279Jup5Cx4nIRG\n67F+HeYmzqniNhth8XR4Jtd0ISLmVCXBXxJh/SvE6dh5q54MucwAIp3s6U7gceCelu/vBR4Aro9w\n25Kwr+e0bPaq6xcg1PJhoboQXnznM3YNfZi4mTfutIo34eVT4Jx0GOHq6yWKOOOJHFjyGlBpupIo\nilu2DsViWKQQeBM4qs3PQzHaP0y84yx6lV/OykueYPWF82Oyz7iSfiJc8yN4ZJ1fZ7KKXy3sDdeH\nYMmdQL3pajopanY6FahDgNal5T8FjgWu6GxREnNBGH47/K8OkIqPNAbguiKYcj80fmG6mi6Imp1O\ntWN/AL7AGkM/HSvUxb2aYfWLOkAq/jI1G/6zKN4PhIYz2SGrQ3edIZfDr86EWzearkTEWVVJcFk2\nzLkD946dRxPzDl3iUsWb8HKddR1FES97MgeWTyP+wrxdCnQJVwPLX4DfZVsn+BfxotkD4NUq2Pqe\n6UrspjettFH7CcyaCw/nma5ExH7lKXBfBix7lPiZ1dJpCnRpKwSbXoLnqmDmANPFiNinMQD/nQsL\nngc8eZxIgS6R1MHyR+G36VCWYroYEXs8kQcz5sPuD0xX4hQFukRTBp89D/fkQYNmI0mcm9cPnt4J\npc8T58v726NAl3bUzIUZH8GTGk+XOFaZDP/TG5Y8CtSYrsZJCnRpTwg2vABP74AP+psuRqTrGgNw\nbx588jKwznQ1TlOgS0f2whePwn29oCK545uLuMkzOfDWQtg5y3QlsaBAl84ohc9esjodjadLvPi4\nLzy5F0r/RhxftKIrFOjSSTtnw5vz4LdDffLekLi2MgPu7AXLHgJ2m64mVhTo0lnNsOlZeGEpPJxv\nuhiR6MpS4BeDYf6D0OD5cfNwCnTpigZY9zg8vhGeyzFdjMjBtiXBz3Phw6egLp5OiWsLBbp0VS2s\n+jPcvxOmDTZdjMjXahLgtgKY+QrsmWe6GhMU6NIdu2DZ/XBvI8ztZ7oYEetg/X8Nhbfegup3TFdj\nigJduqsKFt4Pd6TB571MFyN+1gzcNxSmzoPKf+DhlaAdUaBLT2yETx+A2/rDGp1DXQx5qMA6WL/p\nWaDJdDUmKdClh/avgrmPwC+yrVOTisTSsznwxAZY/zjQYLoa0xToYoN9C+GDp+DGXHXqEhvNwON5\n8MdK6yA9taYrcgNdU1RslHo0nHIz/LEaxu0xXY14VWMA/lgAf1sFax/B4yfciiBqdirQxWZJI+C4\nn8F9++D0HaarEa+pD8JdQ+HVz2DDU8A+0xUZoECXmCqAY34B/5UIF1aZLka8Ymci3F4A02ZCxctA\no+mKDFGgS8wNhlE/h1/0g2s2my5G4l1FMvwqD959DbZNw8dTE1GgiyF9YMRP4P8NhVvKdAxeumdt\nGvwyG2Y/C7vn4O8wBwW6GJQORT+Cq4+COzZAst/fjNIli3vBr/pbU2P3LTRdjUso0MWoZMi/Bs47\nFX5TDrl+PJAlXfZ6FvwpAAv+DPu/NF2NiyjQxbgg9C6G466CO2pggmbASBQ1CfCnfPj7elj1BFBp\nuiKXUaCLaxTCUTfCdQPhxjINwciBlmfAf2fB3Leh8jW0+jMSBbq4SgYUXAWTToa7N0NBvemCxLRm\n4B9D4M/NsPhJqF9suiIXU6CL6wQg8xQYfw38uh7O2ma6IDFlZyL8IR/+uRrWPAnotdA+Bbq4Vh4c\neRNcnQW3bIJUXbDUVz7vBfcMhI+mwdY3gf2mK4oDCnRxtTTI/Q5MnAC3V8LIvaYLEqc1BuCFHHhs\nHyx+HBqXm64ojijQxfUCkDIejroavpMBN5RDb1+f29q75veBh/rDx59C6ctAtemK4owCXeJGBmSf\nD2POhhvr4fxKrTD1iopkeCwH/lUFK56DxhVo1Wd3KNAl7uTD8KugeCTcUglHahgmbjUE4JUh8HQA\nlr0G1bMTrSdqAAAGKElEQVTQdMSeUKBLXApawzCjroLvZML3N0Nfv55hL07N7wMPDoD5n0Dp39EM\nFjso0CWuaRgm7rQOr0yrhJXPQeNKNLxiFwW6eEIeHHYVnDISrtgNE7cr2N2mKgn+ng3/aNbwimMU\n6OIZQUgcCYdeBGOGw5V74ZytkKjuz6iNqTAlC6Y1wNrpsPV9NHvFKQp08ZwAcCgcdj6MHANX7IML\nK7UwKdZWpcNLg+CdOvjq31A9F/9d4zPWFOjiaQVQdA4cdgJ8pxG+uUVz2J32eS94cQDM3glrpsHe\njwGdkyc2FOjiC1mQOxkOLYZvAd+uhCyN39qmGZjXD6b0hY+2wOrXWy46oeX6saVAF1/pD1kToeAM\nOD4FJtfDhO2Qqa69W1alw3sDYHoQ1pfC6mnQuBTQv6cZCnTxpWRgBBxyCmQfA8UJcOYeOLlaB1E7\nUpEM7w2GdxNg5XbYOBuqF2FdbEL/dmYp0MX3MiBhFAwvhtwRMDkAk3bAmD2a+thqdwLMGgTTk2Fh\nHWyeC5s/BTZgjbeIOyjQRcIMgIyxMGwiHJIDk0MwbheM3eO/WTLlKbCwL/wnFT5ohKpPYd1HwBo0\nNu5WjgT6t4ES4HDgWGBR2O9uB67DGmO7BXivK0WJxEgAyIF+YyF3LPQaBqOA8cDo3d4M+NYAX5IK\nC0JQVgO7v4DSz6FpBVBnukLpkCOBfjjWx7AngZ/zdaAfAbyMFfK5wExgOAd/ZIt1oBcDc2K4v1gr\nRo+vp1KBAuh32MEBP3Y3jHYw4P9SCD8otf9+owb4UmgqBaqIzZh4Md59fRYT28cWNTsTe3CnX0b5\n+YXAFKyPa6XAWuA4YH4P9mWHYrz7ggI9PjvUA6uhejVUvwWkwscFMLUl4HsfAsNCUBiAoQHIqYWC\nOiis6/kMmpmF3Q/0ZqAiBUrTYFM6lCdBaTOsD0B5DexeaCDA2yrGu6/PYlzy2HoS6NHkcGB4b8Lq\n1EXiTYSA/89gYAD0GgzZQyEjDxKGwMAEyG+GwiAUNENuHfRrsII+vQkymqyvuzK7phmoD0JNAtQm\nwN4E2JsIm1OtrbTZOl5ZHoS6XdBQDjs3QsUmaNoObMc6u6FmpfhER4E+A8iO8PM7gDe7sB+9oMQL\n6oGN1rYHawOsaTK9gQHWNnAIDCyA5D4QTINgOtZwTgqkAGlAWgjSQ5COtS3uDVfmQS1QG4D6ANQF\noGk/sA9C9dBcC811sKcCyjfC/m1Yob0DrdIUm7wPHB32/a9btlbvAsdH+HuLsYJemzZt2rR1fluM\ng94Hjgn7/oiWHSYDRcBXaDaLiIirXQyUYU1z2gK8E/a7O7AOhn4JnBX70kRERERExLW+DSzHWuR0\ndJvf3Y61Iu5LYHKM63JCCdasos9btrONVmOPs7GenzXAbYZrcUIp8AXW8/Wp2VJs8QzW+V6Whv2s\nP9YEi9VYiwz7GqjLLpEeXwnee9+51uFYC5ve58BAbx3rTwIKsYaI4v2kHncDPzNdhI0SsJ6XQqzn\naTEw0mRBDliPFXhecSowjgMD74/Ar1q+vg34fayLslGkx+ea9128B1hnfInVGbQVbQFUvPPSAejj\nsJ6XUqzn6RWs581rvPScfcjBl567AHiu5evngItiWpG9Ij0+cMlz6IdAjyYH62NSK68sgPoxsAT4\nK/H90Ras56Ms7HuvPEfhQlinx1gA3GC4FqdkYQ1T0PJnlsFanOKK951XAn0G1kegtts3ung/IZvr\nckK0x3oB8DjWVNGxQAXwgKEa7RIPz0dPnYz1Ef4c4Casj/Re1jqX2ktc875zYum/CWd24++UA/lh\n3+e1/MztOvtYn6Zrq3ndqO1zlM+Bn6q8oKLlz63A61jDTB+aK8cRlVgrzrcAQ7DOJ+Ml4Y/H6PvO\nKx16Z4WPc/0L+A5fL4A6jPifZTAk7OuLOfDATTxagPW8FGI9T5dhPW9ekQ70avk6A2umVbw/Z5H8\nC7im5etrgDcM1uIEr73vXM1PC6Cex5oCtwTrTeOFscpzgFVYz9PthmuxWxHWzJ3FwDK88fimAJuB\nBqz33fewZvHMxBvTFts+vuvw5vtOREREREREREREREREREREREREREREREREROLF/wduwsRoyODF\nrwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1109655d0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "print 'Blue line intersects clipped shape:', line.intersects(clipped_shape)\n",
    "print 'Green line intersects clipped shape:', line2.intersects(clipped_shape)\n",
    "\n",
    "fig = plt.figure()\n",
    "ax = fig.add_subplot(111)\n",
    "\n",
    "ax.plot(*np.array(line).T, color='blue', linewidth=3, solid_capstyle='round')\n",
    "ax.plot(*np.array(line2).T, color='green', linewidth=3, solid_capstyle='round')\n",
    "ax.add_patch(descartes.PolygonPatch(clipped_shape, fc='blue', alpha=0.5))\n",
    "ax.axis('equal')\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAEACAYAAACTXJylAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAADt9JREFUeJzt3W2MXNV9gPFn7bXrDXZoUkLcYldbGiCp0vJSB1GC6qXI\nyEmjIKT0CwpRUgVVrZTQlLqJoVE3ipqmQim0SI1UlbSJW5qUNIE44UVuwuZNhChhbQhvsQ24phhc\n200xuDQ23n64Y3s8O7Ozs/fcOfec+/ykK2ZmZ+89ls3DnzOzHpAkSZIkSZIkSZIkSZIkSQJgMTAN\nbI69EEnK3aJA57kWeBSYCXQ+SVIPIcK9Cng78PfASIDzSZLmECLcNwEbgKMBziVJ6qNsuN8B7KXY\n33balqQhKBvbTwBXA0eAZcCrgX8D3tP2nK3AuSWvI0lNsw04r+qLrKX7u0pSecFyMvYCKjYZewEV\nm4y9gApNxl5AxSZLfv8S4Jdg2WXwK38MF/4dvOdW+Jtb4Zt/DQc/DjOT8OVPwxs+xPB3ByYX+H09\n2zm6wBMOfCFJCmQJsAqWnQlnng/Lz4Y3LoY1wLkvwgXPw/JXTv6WIyOwaQx23EEGnQoZ7m+2DkkK\naQGh7vTV0+GhbcDT1S+3eqEn7pRNxV5AxaZiL6BiU7EXUKGp2Auo2FTH/QChbhd92p4KfcJh7PXM\nDOk6ktLUL9QvDBbqTne8Hjb8GHbcTFrbJD3b6cQtadgCT9RziT5tV8JwS6raEEPdKa+97WMMt6TQ\nIoa6XZ7TNhhuSeXVJNSd8py2wXBLGlxNQ90u32kbDLek/hIIdad8p20w3JJmSzDU7fKetsFwS0o+\n1J3ynrbBcEtNlFmo2+U/bYPhlpqgS6jPWQxrRopQ/3rCoe6U/7QNhlvKUYNC3a4Z0zYYbikHDQ11\np2ZM22C4pRQZ6lmaM22D4ZZSYKj7as60DYZbqiNDPZBmTdtguKU6MNSlNGvaBsMtxWCog2netA3l\nw72M4nMmfwZYCtwJbCy7KCkzhroyzZu2oXy4XwYuBQ61zvUd4JLWP6WmMtRD0cxpG8JslRxq/XMp\nsBg4EOCcUkoMdRTNnLYhTLgXAQ8Cvwx8Gng0wDmlOjPU0TV32oYw4T4KnAecCtwLTFDBx9FLNTEG\nb5yEt5xmqGNq7rQNYd9V8j/A1yj+irGpjq9Ntt2e6vJ1KRErLoGrToeP7oq9kubKdtqeaB19jZS8\n0GnAEeAnwBjFxP0x4Ottz5kJcB2pDsZgzafg9hdg/OXYi2mur7wervsx7LiZvMLdqWc7y07cPw98\nlmKfexGwiZOjLWVkxSXwzjEY3xt7Jc11ZAQ+m+O0PZCy4X4YuCDEQqSaG4NzroSrn4+9kGa7q9F7\n28csir0AKQ3Hp223SKJx2j7GcEv9OW3XgtP2MYZb6stpOz6n7XaGW5qb03YtOG23M9zSnJy243Pa\n7mS4pd6ctmvBabuT4ZZ6ctqOz2m7G8Mtdee0XQtO290Ybqkrp+34nLZ7MdzSbE7bteC03YvhlmZx\n2o7PaXsuhls6mdN2LThtz8VwSydx2o7Pabsfwy2d4LRdC07b/Rhu6Tin7fictufDcEsFp+1acNqe\nD8MtAU7bdeC0PV+GW3Largmn7fky3JLTdg04bQ8iRLhXA/cBjwA/Aj4Y4JzSsDht14LT9iDKflgw\nwGHgQ8BWYDnwQ2AL8FiAc0sV85Pb43PaHlSIifs5imgDvEgR7F8IcF6pak7bteC0PajQe9zjwPnA\nA4HPK1XAve34nLYXIsRWyTHLgS8C11JM3u0m225PtQ4pJqftKI4Cj50CW38WHlwE0yPw5Hdx2gaY\naB19jQS64BLgq8DdwM0dX5sJeB0pkBXrYMNV8NFdsVeSt26hPrAHnpuG558AdgEvRF5kXfVsZ4iJ\newS4FXiU2dGW6shpuzI9Q32voQ4nRLjfCrwbeAiYbj22EbgnwLmlCvhOknAMdQzD2MJwq0R1MgZr\nPgW3v+CLkgvh1scQVbpVIiXEaXswTtR1ZLjVJO5t92WoU2C41SBO27MZ6hQZbjWF0zZgqPNguNUQ\nTZ22DXWODLeaoEHTtqFuAsOtBsh52jbUTWS4lbvMpm1DLcOt7KU+bRtqzWa4lbMEp21Drf4MtzKW\nwrRtqDU4w61c1XTaNtQqz3ArU3WZtg21wjPcylHEadtQq3qGWxka5rRtqDV8hlu5qXjaNtSKz3Ar\nM6GnbUOt+jHcykmAadtQq/4MtzKykGnbUCs9IcL9GeC3gb3ArwY4n7QQ85y2DbXSFyLc/wDcAnwu\nwLmkBeo1bRtq5SdEuL8NjAc4j7RQbdO2oVb+3ONWBl71djjjDXDLbkOtJhhWuCfbbk+1DimQU38C\n3/4C3Pl9DLXSNdE6hmYceLjH12aGuA5JykXPdi4a5iokSfXwL8CzwP8Bu4H3dXzdiVuSBhe1nYZb\nkgbnVokk5cJwS1JiDLckJcZwS1JiDLckJcZwS1JiDLckJcZwS1JiDLckJcZwS1JiDLckJcZwS1Ji\nDLckJcZwS1JiDLckJcZwS1JiDLckJcZwS1JiDLckJSZEuNcDjwPbgQ8HOJ8kqUKLgR3AOLAE2Aq8\nqeM5fliwJA2usg8LvpAi3E8Dh4HPA1eUPKckaQ5lw30GsLvt/jOtxyRJFRkt+f3z3QaZbLs91Tok\nSSdMtI7KXQTc03Z/I7NfoHSPW5IGV1k7R4GdFC9OLsUXJyUplErb+TbgCYoXKTcO++KSlKmo7TTc\nkjS4yt4OKEkaMsMtSYkx3JKUGMMtSYkx3JKUGMMtSYkx3JKUGMMtSYkx3JKUGMMtSYkx3JKUGMMt\nSYkx3JKUGMMtSYkx3JKUGMMtSYkx3JKUGMMtSYkx3JKUmDLh/h3gEeAV4IIwy5Ek9VMm3A8DVwLf\nCrQWSdI8jJb43seDrUKSNG/ucUtSYvpN3FuAlV0evx7YPMB1JttuT7UOSdIJE61jKO5j7hcnZ4a1\nEEnKSM92htoqGQl0HklSha4EdgP/CzwH3N3jeU7ckjS4qO003JI0uMq3SiRJQ2K4JSkxhluSEmO4\nJSkxhluSEmO4JSkxhluSEmO4JSkxhluSEmO4JSkxhluSEmO4JSkxhluSEmO4JSkxhluSEmO4JSkx\nhluSEmO4JSkxhluSElM23DcCjwHbgC8Bp5ZekSSpUus4Ef9Pto5OfliwJA2usg8L3gIcbd1+AFhV\n8nySpCHaDFzV5fFUJu7zgNXAktgLkSTmaOfoPL55C7Cyy+PXU8Qa4Abgp8BtPc4x2XZ7qnXUyVs5\nmz9lKc/y38ARnuEQj3OQJ4FngeeAw3GXKClzE62jr5EAF3svcA1wGfByl6/PBLpOlS7mcq7hYnZx\nhEXs4xT2sYL9jPJfzBhzSRH0bOd8Ju65rAc2AGvpHu30jHKUlRxkJQePP1bEfDn7uJT9rDPmkmIq\nG+5bgKUU2ykA9wN/UPKc9WPMJdVI2XCfFWQVKTLmkiIpG261M+aShsBwV82YSwrMcMdgzCWVYLjr\nwphLmifDXWcLi/lOYA/GXMqW4U6NMZcaz3DnwJhLjWK4c2XMpWwZ7iYx5lIWDHfTGXMpOYZbsxlz\nqdYMt+bHmEu1Ybi1cMZcisJwKyxjLlXOcKt6xlwKynArDmMuLZjhVn0Yc2leDLfqzZhLs5QJ98eB\nd1J8EvF+ik973x1gTdLcjLkarutHv8/TCjj+L84HgHOB93d5Xs+PmK+Ri7mca7iYXbEXooCKmJ/C\nPlawn1FjrsT0bGeZiftg2+3lwL4S55LCczJXpsrucf85cDVwCLio/HKkihlzZaDfFsYWYGWXx68H\nNrfd/whwDvC+Ls91q0TpmWub5UWmeYk7Kf5sS1VZ8FbJunle4Dbgrjm+Ptl2e6p1SPXVOZkfYIwH\nWclOXsth9mO0Fd5E6+irzFbJWcD21u0rgOk5njtZ4jpSPAcYYyuvYzsvsY9/5DD3Ay/HXpayNMXJ\nQ+2f9XpimXD/BcX2yCvATuD3S5xLqpeTg73JYKtOyoT7XcFWIdWFwVYC/MlJCQy2kmK41WwGWwky\n3Gomg62EGW41i8FWBgy3msFgKyOGW3kz2MqQ4VaeDLYyZriVF4OtBjDcyoPBVoMYbqXNYKuBDLfS\nZLDVYIZbaTHYkuFWIgy2dJzhVr0ZbGkWw616MthST4Zb9WKwpb4Mt+rBYEvzZrgVl8GWBma4FYfB\nlhYsRLivA24ETgMOBDifcmawpdIWlfz+1cA6YFeAtcS1h5Wxl1CpHzAe9foHGOMb/CK38yruZxN7\n2MBh7iNctCcCnaeOJmIvoGITsRdQsYnQJywb7r8C/iTEQqLbn3m4n4wU7uqDfcxE4PPVyUTsBVRs\nIvYCKjYR+oRltkquAJ4BHgq0FuXELRGpMv3CvQW6TqI3ABuBy9seGwm1qAhmOMRi7mJV7IVUZi+v\nHtqv7xVGeJZDBluqxkJj+2bg68Ch1v1VwH8CFwJ7O567FTh3gdeRpKbaBpxX5QWeAl5b5QUkSYWy\nL04eMxPoPJIkSZJy9wHgMeBHwF9GXksVrgOOkt/W1o0Uv2/bgC8Bp8ZdTjDrgceB7cCHI68ltNXA\nfcAjFP++fTDuciqxGJgGNsdeSM4upXgnzZLW/ddFXEsVVgP3kOdrEus4sfX3ydaRusXADmCc4s/k\nVuBNMRcU2EpOvPi2HHiCvH59AH8E/DPwldgLydm/Ar8VexEVuh34NfIMd7srgX+KvYgAfoPiP7TH\nfKR15OoO4LLYiwhoFfDvFANh0Ik71IuTuTgL+E3ge8AUsCbqasJq0g9M/S5wV+xFBHAGsLvt/jOt\nx3I0DpwPPBB5HSHdBGyg2JoMqol/O+BcP1Q0CrwGuAh4C8UEfubwllZa7j8w1evXdz0nJpobgJ8C\ntw1rURVqyru1lgNfBK4FXoy8llDeQfEzLdPk/yP90d0NrG27vwP4uUhrCenNwPMUWyRPAYeBp4HT\nI66pCu8Fvgssi7yOUC7i5K2SjeT3AuUS4F7gD2MvJLBPUPzf0lPAHuAl4HNRV5Sx3wM+1rp9NvAf\nEddSpRz3uNdTvDvhtNgLCWgU2EmxjbCU/F6cHKGI2U2xF1KxtfiukkotATYBDwM/JN//xXmS/MK9\nneKvF55uHX8bdznBvI3i3RY7KCbunFxCsf+7lRO/b+ujrqgaa/FdJZIkSZIkSZIkSZIkSZIkSZIk\nSdLw/D+Dqw+3qgzMcwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x114ab75d0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "start = Polygon([[-2, 1], [2, 2], [3, 4], [-1, 3]])\n",
    "end = Polygon([[-6, -2], [-6, -4], [-3, -4], [-1, -3]])\n",
    "line = LineString([[-8, -5], [5, 5]])\n",
    "\n",
    "fig = plt.figure()\n",
    "ax = fig.add_subplot(111)\n",
    "\n",
    "#ax.plot(*np.array(line).T, color='red', linewidth=1, solid_capstyle='round')\n",
    "ax.add_patch(descartes.PolygonPatch(start, fc='blue', alpha=0.5))\n",
    "ax.add_patch(descartes.PolygonPatch(end, fc='green', alpha=0.5))\n",
    "ax.axis('equal')\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": []
  }
 ],
 "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
}