Day04_geometric_analysis.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 211,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>name</th>\n",
       "      <th>address</th>\n",
       "      <th>tel</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>新北市立圖書館總館</td>\n",
       "      <td>板橋區貴興路139號</td>\n",
       "      <td>02-2953-7868</td>\n",
       "      <td>POINT (295924.3464126067 2766044.031022599)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>八里分館</td>\n",
       "      <td>八里區舊城里舊城路19號5~6樓</td>\n",
       "      <td>02-2610-3385</td>\n",
       "      <td>POINT (290314.0604052807 2782072.96750804)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>三芝分館</td>\n",
       "      <td>三芝區淡金路1段37號</td>\n",
       "      <td>02-2636-2331</td>\n",
       "      <td>POINT (300455.5626360322 2794611.54889705)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>三重五常分館</td>\n",
       "      <td>三重區五華街7巷30號2樓</td>\n",
       "      <td>02-2989-0559</td>\n",
       "      <td>POINT (299953.9171605469 2774884.029527623)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>三重田中分館</td>\n",
       "      <td>三重區忠孝路3段40巷51號2~4樓</td>\n",
       "      <td>02-2988-5482</td>\n",
       "      <td>POINT (298742.620464512 2773123.990503977)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>三重崇德分館</td>\n",
       "      <td>三重區仁孝街62號2樓</td>\n",
       "      <td>02-2973-6987</td>\n",
       "      <td>POINT (300169.5892529833 2774117.205006743)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>三重分館</td>\n",
       "      <td>三重區自強路1段158號</td>\n",
       "      <td>02-2981-4887</td>\n",
       "      <td>POINT (299220.9035412617 2774224.515142027)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>三重東區分館</td>\n",
       "      <td>三重區溪尾街303號</td>\n",
       "      <td>02-2857-5270</td>\n",
       "      <td>POINT (299191.3936854166 2775133.805940764)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>三重南區分館</td>\n",
       "      <td>三重區重安街70號3樓</td>\n",
       "      <td>02-2976-8729</td>\n",
       "      <td>POINT (299582.3565128865 2772199.898055936)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>三重培德分館</td>\n",
       "      <td>三重區忠孝路2段35號2樓</td>\n",
       "      <td>02-2984-2179</td>\n",
       "      <td>POINT (298522.4940371005 2773176.373897539)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>三峽分館</td>\n",
       "      <td>三峽區永安街9巷5號2樓</td>\n",
       "      <td>02-2673-5234</td>\n",
       "      <td>POINT (287895.7292293567 2757921.508451493)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>三峽北大分館</td>\n",
       "      <td>三峽區學成路396號3-4樓</td>\n",
       "      <td>02-8671-5735</td>\n",
       "      <td>POINT (287870.5766551858 2759350.284659256)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>土城柑林埤圖書閱覽室</td>\n",
       "      <td>土城區學享街55號3樓</td>\n",
       "      <td>02-2263-7037</td>\n",
       "      <td>POINT (295626.9146780901 2764425.864276677)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>土城祖田圖書閱覽室</td>\n",
       "      <td>土城區中央路4段283-1號4樓</td>\n",
       "      <td>02-2267-4776</td>\n",
       "      <td>POINT (294322.6471886858 2762537.471766521)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>土城清水圖書閱覽室</td>\n",
       "      <td>土城區清水路243號2樓</td>\n",
       "      <td>02-2265-3293</td>\n",
       "      <td>POINT (296592.4177212388 2764015.961402074)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>土城分館</td>\n",
       "      <td>土城區中正路18號7樓</td>\n",
       "      <td>02-2260-2952</td>\n",
       "      <td>POINT (294704.4472692583 2763112.473279943)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>土城親子分館</td>\n",
       "      <td>土城區中州路28巷30號4樓</td>\n",
       "      <td>02-2267-8901</td>\n",
       "      <td>POINT (291922.6859968934 2761207.394835655)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>中和大同圖書閱覽室</td>\n",
       "      <td>中和區圓通路121巷5號1樓</td>\n",
       "      <td>02-2248-3766</td>\n",
       "      <td>POINT (300048.7626940823 2765237.719257739)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>中和分館</td>\n",
       "      <td>中和區南山路236號6樓</td>\n",
       "      <td>02-2949-2767</td>\n",
       "      <td>POINT (300947.6053299593 2764587.515893403)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>中和員山分館</td>\n",
       "      <td>中和區民享街37號2樓</td>\n",
       "      <td>02-8221-1152</td>\n",
       "      <td>POINT (297765.1688677093 2765857.600809997)</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "          name             address           tel  \\\n",
       "0    新北市立圖書館總館          板橋區貴興路139號  02-2953-7868   \n",
       "1         八里分館    八里區舊城里舊城路19號5~6樓  02-2610-3385   \n",
       "2         三芝分館         三芝區淡金路1段37號  02-2636-2331   \n",
       "3       三重五常分館       三重區五華街7巷30號2樓  02-2989-0559   \n",
       "4       三重田中分館  三重區忠孝路3段40巷51號2~4樓  02-2988-5482   \n",
       "5       三重崇德分館         三重區仁孝街62號2樓  02-2973-6987   \n",
       "6         三重分館        三重區自強路1段158號  02-2981-4887   \n",
       "7       三重東區分館          三重區溪尾街303號  02-2857-5270   \n",
       "8       三重南區分館         三重區重安街70號3樓  02-2976-8729   \n",
       "9       三重培德分館       三重區忠孝路2段35號2樓  02-2984-2179   \n",
       "10        三峽分館        三峽區永安街9巷5號2樓  02-2673-5234   \n",
       "11      三峽北大分館      三峽區學成路396號3-4樓  02-8671-5735   \n",
       "12  土城柑林埤圖書閱覽室         土城區學享街55號3樓  02-2263-7037   \n",
       "13   土城祖田圖書閱覽室    土城區中央路4段283-1號4樓  02-2267-4776   \n",
       "14   土城清水圖書閱覽室        土城區清水路243號2樓  02-2265-3293   \n",
       "15        土城分館         土城區中正路18號7樓  02-2260-2952   \n",
       "16      土城親子分館      土城區中州路28巷30號4樓  02-2267-8901   \n",
       "17   中和大同圖書閱覽室      中和區圓通路121巷5號1樓  02-2248-3766   \n",
       "18        中和分館        中和區南山路236號6樓  02-2949-2767   \n",
       "19      中和員山分館         中和區民享街37號2樓  02-8221-1152   \n",
       "\n",
       "                                       geometry  \n",
       "0   POINT (295924.3464126067 2766044.031022599)  \n",
       "1    POINT (290314.0604052807 2782072.96750804)  \n",
       "2    POINT (300455.5626360322 2794611.54889705)  \n",
       "3   POINT (299953.9171605469 2774884.029527623)  \n",
       "4    POINT (298742.620464512 2773123.990503977)  \n",
       "5   POINT (300169.5892529833 2774117.205006743)  \n",
       "6   POINT (299220.9035412617 2774224.515142027)  \n",
       "7   POINT (299191.3936854166 2775133.805940764)  \n",
       "8   POINT (299582.3565128865 2772199.898055936)  \n",
       "9   POINT (298522.4940371005 2773176.373897539)  \n",
       "10  POINT (287895.7292293567 2757921.508451493)  \n",
       "11  POINT (287870.5766551858 2759350.284659256)  \n",
       "12  POINT (295626.9146780901 2764425.864276677)  \n",
       "13  POINT (294322.6471886858 2762537.471766521)  \n",
       "14  POINT (296592.4177212388 2764015.961402074)  \n",
       "15  POINT (294704.4472692583 2763112.473279943)  \n",
       "16  POINT (291922.6859968934 2761207.394835655)  \n",
       "17  POINT (300048.7626940823 2765237.719257739)  \n",
       "18  POINT (300947.6053299593 2764587.515893403)  \n",
       "19  POINT (297765.1688677093 2765857.600809997)  "
      ]
     },
     "execution_count": 211,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import geopandas as gpd\n",
    "gdf_Lib=gpd.read_file('output/Library.shp',encoding='utf-8')\n",
    "gdf_Lib.crs = {'init' :'epsg:4326'} # 避免資料沒設，這邊再重新給一次\n",
    "gdf_Lib=gdf_Lib.to_crs(epsg=3826)\n",
    "gdf_Lib"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "base=gdf_Lib.head(1).buffer(100).plot()\n",
    "\n",
    "gdf_Lib.head(1).plot(ax=base, marker='o', color='red', markersize=10);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "buffer=gdf_Lib.head(1).buffer(100)\n",
    "area=buffer.area\n",
    "print(area)\n",
    "\n",
    "envelope=buffer.envelope\n",
    "print(envelope)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "convexhull=buffer.convex_hull\n",
    "convexhull\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "base=envelope.plot()\n",
    "convexhull.plot(ax=base,color='brown')\n",
    "gdf_Lib.head(1).plot(ax=base,  color='red');\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "base=buffer.scale(10,5).plot(color='yellow')\n",
    "buffer.plot(ax=base,color='brown')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "p1=gpd.read_file('output/Library.shp',encoding='utf-8').head(1)\n",
    "p1.crs = {'init' :'epsg:4326'} # 避免資料沒設，這邊再重新給一次\n",
    "p1=p1.to_crs(epsg=3826)\n",
    "p1['geometry']=p1.buffer(30).translate(xoff=20.0, yoff=0.0)\n",
    "\n",
    "p2=gpd.read_file('output/Library.shp',encoding='utf-8').head(1)\n",
    "p2.crs = {'init' :'epsg:4326'} # 避免資料沒設，這邊再重新給一次\n",
    "p2=p2.to_crs(epsg=3826)\n",
    "p2['geometry']=p2.buffer(30)\n",
    "intersection = gpd.overlay(p1,p2,  how='intersection')\n",
    "base=p1.plot(color='blue')\n",
    "p2.plot(ax=base,color='brown')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "intersection = gpd.overlay(p1,p2,  how='intersection')\n",
    "base=p1.plot(color='blue')\n",
    "p2.plot(ax=base,color='brown')\n",
    "intersection.plot(ax=base,color='yellow')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "union = gpd.overlay(p1,p2,  how='union')\n",
    "base=p1.plot(color='blue')\n",
    "p2.plot(ax=base,color='brown')\n",
    "union.plot(ax=base,color='yellow')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "difference = gpd.overlay(p1,p2,  how='difference')\n",
    "base=p1.plot(color='blue')\n",
    "p2.plot(ax=base,color='brown')\n",
    "difference.plot(ax=base,color='yellow')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "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
}