notebooks/unit06/unit06_04.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This notebook was prepared by Cayetano Benavent, 2016."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Vector analysis with Shapely and NumPy"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We review in this examples how to use Shapely with vector data and NumPy:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from shapely.wkt import loads\n",
    "from shapely.geometry import LineString"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def runTest01():\n",
    "    wkt_ln = \"LINESTRING(178044.505034367 4080523.32981373,185051.202770193 4115968.97717771,218436.056686834 4125036.46836944,231625.134776712 4158833.48097265,229564.341325156 4174495.51120244,216787.42192546 4189333.22405021,212253.676331916 4213238.42808471,216787.421925395 4231373.41045699)\"\n",
    "\n",
    "    ln = loads(wkt_ln)\n",
    "    ln_np = np.array(ln)\n",
    "\n",
    "    ln_buf = ln.buffer(20000)\n",
    "\n",
    "    return(ln_np, ln_buf.wkt)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def runTest02():\n",
    "    lat = np.array([88.2598888, 87.8656022505306, 87.47131570106119, 87.07729258677767, 86.68423509786567, 86.29117760895366, 85.89812012004164, 85.50506263112963, 85.11200514221763, 84.71894765330562, 84.35408870391326, 84.00746812090851, 83.66121815097686, 83.31589924664645, 82.96712093784359, 82.6290776997575, 82.25749423519197, 81.87712508822979, 81.49675594126761, 81.09931880902852, 80.6998987769223, 80.38219308946239, 80.06394742324582, 79.70348526693819, 79.58169647523991, 79.7679112, 79.7679112])\n",
    "    lon = np.array([49.7934288, 49.86079434688794, 49.92815989377588, 49.99699057427323, 50.07119172258097, 50.14539287088871, 50.21959401919645, 50.2937951675042, 50.367996315811936, 50.44219746411967, 50.59280312395832, 50.792436728275504, 50.99271393197662, 51.19459015098207, 51.390376419073085, 51.60421432299298, 51.746632193203595, 51.87040340164587, 51.99417461008816, 52.02591753929958, 52.04744979306631, 51.826426464592416, 51.58531405704919, 51.43357869745786, 51.39167315475912, 51.3637112, 51.3637112])\n",
    "\n",
    "    coor = np.vstack((lon, lat))\n",
    "    coor_t = coor.T\n",
    "\n",
    "    ln = LineString(coor_t)\n",
    "\n",
    "    ln_wkt = ln.wkt\n",
    "\n",
    "    return(ln_wkt)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "res01 = runTest01()\n",
    "res02 = runTest02()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Test01 array: \n",
      "[[  178044.50503437  4080523.32981373]\n",
      " [  185051.20277019  4115968.97717771]\n",
      " [  218436.05668683  4125036.46836944]\n",
      " [  231625.13477671  4158833.48097265]\n",
      " [  229564.34132516  4174495.51120244]\n",
      " [  216787.42192546  4189333.22405021]\n",
      " [  212253.67633192  4213238.42808471]\n",
      " [  216787.42192539  4231373.41045699]]\n"
     ]
    }
   ],
   "source": [
    "print(\"Test01 array: \\n{}\".format(res01[0]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Test01 buff wkt: \n",
      "POLYGON ((165430.8635830167 4119847.416319962, 165921.6067607095 4121804.94809213, 166607.1255483623 4123703.058599259, 167480.4400630784 4125522.421472323, 168532.6583048801 4127244.512142734, 169753.066694152 4128851.796457624, 171129.2391562242 4130327.909210815, 172647.1636424063 4131657.820771689, 174291.3847992491 4132827.990115369, 176045.1613333884 4133826.502696061, 177890.6364696993 4134643.191759754, 179809.0197671729 4135269.74186108, 203468.2671049733 4141695.710277001, 211125.8141482724 4161318.174574235, 210513.957774473 4165968.283014506, 201632.0458799098 4176282.76134125, 200351.1938988329 4177938.031719375, 199250.3392348049 4179718.093749422, 198341.53762748 4181603.453518976, 197634.7415946002 4183573.463973271, 197137.6914394404 4185606.551026347, 192603.9458458964 4209511.755060846, 192316.5713361304 4211653.548156443, 192261.9556375852 4213813.844177521, 192440.736369668 4215967.42241456, 192850.8263291123 4218089.140585848, 197384.5719225913 4236224.122958128, 197953.4549166156 4238102.577323119, 198703.7193086281 4239916.226132652, 199628.1396423975 4241647.602957251, 200717.81324477 4243280.033691626, 201962.2459632786 4244797.79713543, 203349.453230568 4246186.276397005, 204866.0754823294 4247432.09966204, 206497.5068172085 4248523.268971459, 208228.0356596191 4249449.275768317, 210040.9960708057 4250201.202100954, 211918.9282509435 4250771.806507742, 213843.7466865555 4251155.593756319, 215796.9143238942 4251348.867765682, 217759.6210909019 4251349.767201466, 219712.9650484876 4251158.283401614, 221638.1344265334 4250776.260459794, 223516.5887915241 4250207.37746577, 225330.2376010571 4249457.113073757, 227061.6144256555 4248532.692739988, 228694.0451600314 4247443.019137615, 230211.8086038354 4246198.586419106, 231600.2878654097 4244811.379151817, 232846.1111304456 4243294.756900055, 233937.2804398643 4241663.325565176, 234863.2872367223 4239932.796722766, 235615.2135693592 4238119.83631158, 236185.8179761468 4236241.904131441, 236569.6052247238 4234317.085695829, 236762.8792340868 4232363.91805849, 236763.7786698714 4230401.211291483, 236572.2948700193 4228447.867333897, 236190.2719281987 4226522.697955851, 232721.9215191716 4212649.296320994, 235438.830068518 4198323.778517858, 244719.7173707063 4187545.9739114, 245908.7016030077 4186022.084288478, 246945.0326530277 4184390.538689276, 247819.0313880259 4182666.575435976, 248522.5348278711 4180866.296013875, 249048.9723857055 4179006.514686618, 249393.427235963 4177104.601454201, 251454.220687519 4161442.571224411, 251615.9858384124 4159438.356024965, 251575.6911367771 4157428.026980566, 251333.7438659973 4155431.90372588, 250892.5895362567 4153470.162308903, 250256.68716626 4151562.631259203, 237067.609076382 4117765.618655994, 236308.7205333188 4116060.506747945, 235394.1915413952 4114433.560171473, 234331.9860898644 4112898.946859104, 233131.3541783169 4111470.03067796, 231802.7512649446 4110159.255053247, 230357.7472173158 4108978.034607226, 228808.9255585422 4107936.655757299, 227169.7738862334 4107044.187138842, 225454.5664185037 4106308.400632846, 223678.2396898541 4105735.70368607, 202265.7300857639 4099919.960333928, 197664.8442215433 4076644.890671478, 197190.2134119071 4074740.43687661, 196531.1989644863 4072891.675501978, 195694.1475486635 4071116.411130254, 194687.1204271521 4069431.740527797, 193519.8158216292 4067853.887993443, 192203.4755135595 4066398.049109639, 190750.7765796946 4065078.244400681, 189175.7093048937 4063907.184307392, 187493.4424480323 4062896.146778638, 185720.1771585628 4062054.868658511, 183872.990950588 4061391.451915196, 181969.6732370655 4060912.285614592, 180028.5540080358 4060621.984390114, 178068.3273027983 4060523.344001247, 176107.8711760931 4060617.314408858, 174166.0658921149 4060902.990626554, 172261.612097247 4061377.62143619, 170412.8507226145 4062036.635883611, 168637.5863508907 4062873.687299433, 166952.9157484337 4063880.714420945, 165375.0632140797 4065048.019026468, 163919.2243302764 4066364.359334537, 162599.4196213174 4067817.058268402, 161428.3595280286 4069392.125543203, 160417.321999275 4071074.392400065, 159576.0438791477 4072847.657689534, 158912.627135833 4074694.843897509, 158433.4608352295 4076598.161611031, 158143.1596107509 4078539.280840061, 158044.5192218838 4080499.507545299, 158138.4896294954 4082459.963672004, 158424.1658471907 4084401.768955982, 165430.8635830167 4119847.416319962))\n"
     ]
    }
   ],
   "source": [
    "print(\"Test01 buff wkt: \\n{}\".format(res01[1]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Test02 wkt: \n",
      "LINESTRING (49.7934288 88.2598888, 49.86079434688794 87.8656022505306, 49.92815989377588 87.47131570106119, 49.99699057427323 87.07729258677767, 50.07119172258097 86.68423509786567, 50.14539287088871 86.29117760895366, 50.21959401919645 85.89812012004164, 50.2937951675042 85.50506263112963, 50.36799631581194 85.11200514221763, 50.44219746411967 84.71894765330562, 50.59280312395832 84.35408870391326, 50.7924367282755 84.00746812090851, 50.99271393197662 83.66121815097686, 51.19459015098207 83.31589924664645, 51.39037641907309 82.96712093784359, 51.60421432299298 82.6290776997575, 51.7466321932036 82.25749423519197, 51.87040340164587 81.87712508822979, 51.99417461008816 81.49675594126761, 52.02591753929958 81.09931880902852, 52.04744979306631 80.6998987769223, 51.82642646459242 80.38219308946239, 51.58531405704919 80.06394742324582, 51.43357869745786 79.70348526693819, 51.39167315475912 79.58169647523991, 51.3637112 79.7679112, 51.3637112 79.7679112)\n"
     ]
    }
   ],
   "source": [
    "print(\"Test02 wkt: \\n{}\".format(res02))"
   ]
  }
 ],
 "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.4.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}