new_sites.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<pre>✔️ 803 ms (2024-12-02T16:22:16/2024-12-02T16:22:17)</pre>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%load_ext autotime\n",
    "import geopandas as gpd\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "from scipy.spatial import KDTree\n",
    "from tqdm.auto import tqdm\n",
    "from glob import glob\n",
    "from shapely.geometry import LineString, Point\n",
    "import folium\n",
    "pd.set_option('display.max_columns', None)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<pre>✔️ 17.3 ms (2024-12-02T16:22:17/2024-12-02T16:22:17)</pre>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "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>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>LINESTRING (176.9377 -39.33628, 176.93654 -39....</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry\n",
       "0  LINESTRING (176.9377 -39.33628, 176.93654 -39...."
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "new_shoreline = gpd.GeoDataFrame.from_features({\"type\":\"FeatureCollection\",\"features\":[{\"type\":\"Feature\",\"properties\":{},\"geometry\":{\"type\":\"LineString\",\"coordinates\":[[176.937701,-39.336281],[176.936542,-39.336397],[176.935598,-39.33731]]}}]}, crs=\"EPSG:4326\")\n",
    "new_shoreline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<pre>✔️ 2.94 s (2024-12-02T16:22:17/2024-12-02T16:22:20)</pre>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc=\"&lt;!DOCTYPE html&gt;\n",
       "&lt;html&gt;\n",
       "&lt;head&gt;\n",
       "    \n",
       "    &lt;meta http-equiv=&quot;content-type&quot; content=&quot;text/html; charset=UTF-8&quot; /&gt;\n",
       "    \n",
       "        &lt;script&gt;\n",
       "            L_NO_TOUCH = false;\n",
       "            L_DISABLE_3D = false;\n",
       "        &lt;/script&gt;\n",
       "    \n",
       "    &lt;style&gt;html, body {width: 100%;height: 100%;margin: 0;padding: 0;}&lt;/style&gt;\n",
       "    &lt;style&gt;#map {position:absolute;top:0;bottom:0;right:0;left:0;}&lt;/style&gt;\n",
       "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://code.jquery.com/jquery-3.7.1.min.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/js/bootstrap.bundle.min.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/css/bootstrap.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap-glyphicons.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6.2.0/css/all.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css&quot;/&gt;\n",
       "    \n",
       "            &lt;meta name=&quot;viewport&quot; content=&quot;width=device-width,\n",
       "                initial-scale=1.0, maximum-scale=1.0, user-scalable=no&quot; /&gt;\n",
       "            &lt;style&gt;\n",
       "                #map_2bd106d30d5ce78aadf7430853bbc651 {\n",
       "                    position: relative;\n",
       "                    width: 100.0%;\n",
       "                    height: 100.0%;\n",
       "                    left: 0.0%;\n",
       "                    top: 0.0%;\n",
       "                }\n",
       "                .leaflet-container { font-size: 1rem; }\n",
       "            &lt;/style&gt;\n",
       "        \n",
       "    \n",
       "                    &lt;style&gt;\n",
       "                        .foliumtooltip {\n",
       "                            \n",
       "                        }\n",
       "                       .foliumtooltip table{\n",
       "                            margin: auto;\n",
       "                        }\n",
       "                        .foliumtooltip tr{\n",
       "                            text-align: left;\n",
       "                        }\n",
       "                        .foliumtooltip th{\n",
       "                            padding: 2px; padding-right: 8px;\n",
       "                        }\n",
       "                    &lt;/style&gt;\n",
       "            \n",
       "&lt;/head&gt;\n",
       "&lt;body&gt;\n",
       "    \n",
       "    \n",
       "            &lt;div class=&quot;folium-map&quot; id=&quot;map_2bd106d30d5ce78aadf7430853bbc651&quot; &gt;&lt;/div&gt;\n",
       "        \n",
       "&lt;/body&gt;\n",
       "&lt;script&gt;\n",
       "    \n",
       "    \n",
       "            var map_2bd106d30d5ce78aadf7430853bbc651 = L.map(\n",
       "                &quot;map_2bd106d30d5ce78aadf7430853bbc651&quot;,\n",
       "                {\n",
       "                    center: [-39.3367955, 176.9366495],\n",
       "                    crs: L.CRS.EPSG3857,\n",
       "                    zoom: 10,\n",
       "                    zoomControl: true,\n",
       "                    preferCanvas: false,\n",
       "                }\n",
       "            );\n",
       "            L.control.scale().addTo(map_2bd106d30d5ce78aadf7430853bbc651);\n",
       "\n",
       "            \n",
       "\n",
       "        \n",
       "    \n",
       "            var tile_layer_8b0e0be3f9836e045606aa40214c073f = L.tileLayer(\n",
       "                &quot;https://server.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer/tile/{z}/{y}/{x}&quot;,\n",
       "                {&quot;attribution&quot;: &quot;Tiles \\u0026copy; Esri \\u0026mdash; Source: Esri, i-cubed, USDA, USGS, AEX, GeoEye, Getmapping, Aerogrid, IGN, IGP, UPR-EGP, and the GIS User Community&quot;, &quot;detectRetina&quot;: false, &quot;maxZoom&quot;: 18, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
       "            );\n",
       "        \n",
       "    \n",
       "            tile_layer_8b0e0be3f9836e045606aa40214c073f.addTo(map_2bd106d30d5ce78aadf7430853bbc651);\n",
       "        \n",
       "    \n",
       "            map_2bd106d30d5ce78aadf7430853bbc651.fitBounds(\n",
       "                [[-39.33731, 176.935598], [-39.336281, 176.937701]],\n",
       "                {}\n",
       "            );\n",
       "        \n",
       "    \n",
       "        function geo_json_707c6f579b80ed126e4d2ec55ba977cd_styler(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.5, &quot;weight&quot;: 2};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_707c6f579b80ed126e4d2ec55ba977cd_highlighter(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.75};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_707c6f579b80ed126e4d2ec55ba977cd_pointToLayer(feature, latlng) {\n",
       "            var opts = {&quot;bubblingMouseEvents&quot;: true, &quot;color&quot;: &quot;#3388ff&quot;, &quot;dashArray&quot;: null, &quot;dashOffset&quot;: null, &quot;fill&quot;: true, &quot;fillColor&quot;: &quot;#3388ff&quot;, &quot;fillOpacity&quot;: 0.2, &quot;fillRule&quot;: &quot;evenodd&quot;, &quot;lineCap&quot;: &quot;round&quot;, &quot;lineJoin&quot;: &quot;round&quot;, &quot;opacity&quot;: 1.0, &quot;radius&quot;: 2, &quot;stroke&quot;: true, &quot;weight&quot;: 3};\n",
       "            \n",
       "            let style = geo_json_707c6f579b80ed126e4d2ec55ba977cd_styler(feature)\n",
       "            Object.assign(opts, style)\n",
       "            \n",
       "            return new L.CircleMarker(latlng, opts)\n",
       "        }\n",
       "\n",
       "        function geo_json_707c6f579b80ed126e4d2ec55ba977cd_onEachFeature(feature, layer) {\n",
       "            layer.on({\n",
       "                mouseout: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                            geo_json_707c6f579b80ed126e4d2ec55ba977cd.resetStyle(e.target);\n",
       "                    }\n",
       "                },\n",
       "                mouseover: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                        const highlightStyle = geo_json_707c6f579b80ed126e4d2ec55ba977cd_highlighter(e.target.feature)\n",
       "                        e.target.setStyle(highlightStyle);\n",
       "                    }\n",
       "                },\n",
       "            });\n",
       "        };\n",
       "        var geo_json_707c6f579b80ed126e4d2ec55ba977cd = L.geoJson(null, {\n",
       "                onEachFeature: geo_json_707c6f579b80ed126e4d2ec55ba977cd_onEachFeature,\n",
       "            \n",
       "                style: geo_json_707c6f579b80ed126e4d2ec55ba977cd_styler,\n",
       "                pointToLayer: geo_json_707c6f579b80ed126e4d2ec55ba977cd_pointToLayer,\n",
       "        });\n",
       "\n",
       "        function geo_json_707c6f579b80ed126e4d2ec55ba977cd_add (data) {\n",
       "            geo_json_707c6f579b80ed126e4d2ec55ba977cd\n",
       "                .addData(data);\n",
       "        }\n",
       "            geo_json_707c6f579b80ed126e4d2ec55ba977cd_add({&quot;bbox&quot;: [176.935598, -39.33731, 176.937701, -39.336281], &quot;features&quot;: [{&quot;bbox&quot;: [176.935598, -39.33731, 176.937701, -39.336281], &quot;geometry&quot;: {&quot;coordinates&quot;: [[176.937701, -39.336281], [176.936542, -39.336397], [176.935598, -39.33731]], &quot;type&quot;: &quot;LineString&quot;}, &quot;id&quot;: &quot;0&quot;, &quot;properties&quot;: {}, &quot;type&quot;: &quot;Feature&quot;}], &quot;type&quot;: &quot;FeatureCollection&quot;});\n",
       "\n",
       "        \n",
       "    \n",
       "    geo_json_707c6f579b80ed126e4d2ec55ba977cd.bindTooltip(\n",
       "    function(layer){\n",
       "    let div = L.DomUtil.create(&#x27;div&#x27;);\n",
       "    \n",
       "    return div\n",
       "    }\n",
       "    ,{&quot;className&quot;: &quot;foliumtooltip&quot;, &quot;sticky&quot;: true});\n",
       "                     \n",
       "    \n",
       "            geo_json_707c6f579b80ed126e4d2ec55ba977cd.addTo(map_2bd106d30d5ce78aadf7430853bbc651);\n",
       "        \n",
       "&lt;/script&gt;\n",
       "&lt;/html&gt;\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
      ],
      "text/plain": [
       "<folium.folium.Map at 0x7fd52ca38910>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "new_shoreline.explore(tiles=\"ESRI.WorldImagery\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<pre>✔️ 4.56 s (2024-12-02T16:22:20/2024-12-02T16:22:24)</pre>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "shorelines = gpd.read_file(\"shorelines.geojson\")\n",
    "transects = gpd.read_file(\"transects_extended.geojson\")\n",
    "poly = gpd.read_file(\"polygons.geojson\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<pre>✔️ 14.2 ms (2024-12-02T16:22:25/2024-12-02T16:22:25)</pre>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Latest siteid is nzd0562, so this new site will be nzd0563\n"
     ]
    },
    {
     "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>geometry</th>\n",
       "      <th>id</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>LINESTRING (176.9377 -39.33628, 176.93654 -39....</td>\n",
       "      <td>nzd0563</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry       id\n",
       "0  LINESTRING (176.9377 -39.33628, 176.93654 -39....  nzd0563"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "nz_shorelines = shorelines[shorelines.id.str.startswith(\"nzd\")]\n",
    "latest_siteid = nz_shorelines.id.max()\n",
    "new_siteid = f\"nzd0{int(latest_siteid[4:])+1}\"\n",
    "print(f\"Latest siteid is {latest_siteid}, so this new site will be {new_siteid}\")\n",
    "new_shoreline[\"id\"] = new_siteid\n",
    "new_shoreline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<pre>✔️ 1.27 s (2024-12-02T16:22:25/2024-12-02T16:22:26)</pre>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "pd.concat((shorelines, new_shoreline)).to_file(\"shorelines.geojson\", driver=\"GeoJSON\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<pre>✔️ 56.1 ms (2024-12-02T16:22:26/2024-12-02T16:22:26)</pre>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc=\"&lt;!DOCTYPE html&gt;\n",
       "&lt;html&gt;\n",
       "&lt;head&gt;\n",
       "    \n",
       "    &lt;meta http-equiv=&quot;content-type&quot; content=&quot;text/html; charset=UTF-8&quot; /&gt;\n",
       "    \n",
       "        &lt;script&gt;\n",
       "            L_NO_TOUCH = false;\n",
       "            L_DISABLE_3D = false;\n",
       "        &lt;/script&gt;\n",
       "    \n",
       "    &lt;style&gt;html, body {width: 100%;height: 100%;margin: 0;padding: 0;}&lt;/style&gt;\n",
       "    &lt;style&gt;#map {position:absolute;top:0;bottom:0;right:0;left:0;}&lt;/style&gt;\n",
       "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://code.jquery.com/jquery-3.7.1.min.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/js/bootstrap.bundle.min.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/css/bootstrap.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap-glyphicons.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6.2.0/css/all.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css&quot;/&gt;\n",
       "    \n",
       "            &lt;meta name=&quot;viewport&quot; content=&quot;width=device-width,\n",
       "                initial-scale=1.0, maximum-scale=1.0, user-scalable=no&quot; /&gt;\n",
       "            &lt;style&gt;\n",
       "                #map_dbd833a77680d77f8a419b432d63bff6 {\n",
       "                    position: relative;\n",
       "                    width: 100.0%;\n",
       "                    height: 100.0%;\n",
       "                    left: 0.0%;\n",
       "                    top: 0.0%;\n",
       "                }\n",
       "                .leaflet-container { font-size: 1rem; }\n",
       "            &lt;/style&gt;\n",
       "        \n",
       "    \n",
       "                    &lt;style&gt;\n",
       "                        .foliumtooltip {\n",
       "                            \n",
       "                        }\n",
       "                       .foliumtooltip table{\n",
       "                            margin: auto;\n",
       "                        }\n",
       "                        .foliumtooltip tr{\n",
       "                            text-align: left;\n",
       "                        }\n",
       "                        .foliumtooltip th{\n",
       "                            padding: 2px; padding-right: 8px;\n",
       "                        }\n",
       "                    &lt;/style&gt;\n",
       "            \n",
       "&lt;/head&gt;\n",
       "&lt;body&gt;\n",
       "    \n",
       "    \n",
       "            &lt;div class=&quot;folium-map&quot; id=&quot;map_dbd833a77680d77f8a419b432d63bff6&quot; &gt;&lt;/div&gt;\n",
       "        \n",
       "&lt;/body&gt;\n",
       "&lt;script&gt;\n",
       "    \n",
       "    \n",
       "            var map_dbd833a77680d77f8a419b432d63bff6 = L.map(\n",
       "                &quot;map_dbd833a77680d77f8a419b432d63bff6&quot;,\n",
       "                {\n",
       "                    center: [-39.33663417765773, 176.93651847639467],\n",
       "                    crs: L.CRS.EPSG3857,\n",
       "                    zoom: 10,\n",
       "                    zoomControl: true,\n",
       "                    preferCanvas: false,\n",
       "                }\n",
       "            );\n",
       "            L.control.scale().addTo(map_dbd833a77680d77f8a419b432d63bff6);\n",
       "\n",
       "            \n",
       "\n",
       "        \n",
       "    \n",
       "            var tile_layer_08c2f8e8e3d853af0be1019a36ccbfa5 = L.tileLayer(\n",
       "                &quot;https://tile.openstreetmap.org/{z}/{x}/{y}.png&quot;,\n",
       "                {&quot;attribution&quot;: &quot;\\u0026copy; \\u003ca href=\\&quot;https://www.openstreetmap.org/copyright\\&quot;\\u003eOpenStreetMap\\u003c/a\\u003e contributors&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 19, &quot;maxZoom&quot;: 19, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
       "            );\n",
       "        \n",
       "    \n",
       "            tile_layer_08c2f8e8e3d853af0be1019a36ccbfa5.addTo(map_dbd833a77680d77f8a419b432d63bff6);\n",
       "        \n",
       "    \n",
       "            map_dbd833a77680d77f8a419b432d63bff6.fitBounds(\n",
       "                [[-39.33855011832538, 176.9337392505123], [-39.33471823699008, 176.93929770227703]],\n",
       "                {}\n",
       "            );\n",
       "        \n",
       "    \n",
       "        function geo_json_c79781e85236c459198fca1d4a2a3ec0_styler(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.5, &quot;weight&quot;: 2};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_c79781e85236c459198fca1d4a2a3ec0_highlighter(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.75};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_c79781e85236c459198fca1d4a2a3ec0_pointToLayer(feature, latlng) {\n",
       "            var opts = {&quot;bubblingMouseEvents&quot;: true, &quot;color&quot;: &quot;#3388ff&quot;, &quot;dashArray&quot;: null, &quot;dashOffset&quot;: null, &quot;fill&quot;: true, &quot;fillColor&quot;: &quot;#3388ff&quot;, &quot;fillOpacity&quot;: 0.2, &quot;fillRule&quot;: &quot;evenodd&quot;, &quot;lineCap&quot;: &quot;round&quot;, &quot;lineJoin&quot;: &quot;round&quot;, &quot;opacity&quot;: 1.0, &quot;radius&quot;: 2, &quot;stroke&quot;: true, &quot;weight&quot;: 3};\n",
       "            \n",
       "            let style = geo_json_c79781e85236c459198fca1d4a2a3ec0_styler(feature)\n",
       "            Object.assign(opts, style)\n",
       "            \n",
       "            return new L.CircleMarker(latlng, opts)\n",
       "        }\n",
       "\n",
       "        function geo_json_c79781e85236c459198fca1d4a2a3ec0_onEachFeature(feature, layer) {\n",
       "            layer.on({\n",
       "                mouseout: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                            geo_json_c79781e85236c459198fca1d4a2a3ec0.resetStyle(e.target);\n",
       "                    }\n",
       "                },\n",
       "                mouseover: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                        const highlightStyle = geo_json_c79781e85236c459198fca1d4a2a3ec0_highlighter(e.target.feature)\n",
       "                        e.target.setStyle(highlightStyle);\n",
       "                    }\n",
       "                },\n",
       "            });\n",
       "        };\n",
       "        var geo_json_c79781e85236c459198fca1d4a2a3ec0 = L.geoJson(null, {\n",
       "                onEachFeature: geo_json_c79781e85236c459198fca1d4a2a3ec0_onEachFeature,\n",
       "            \n",
       "                style: geo_json_c79781e85236c459198fca1d4a2a3ec0_styler,\n",
       "                pointToLayer: geo_json_c79781e85236c459198fca1d4a2a3ec0_pointToLayer,\n",
       "        });\n",
       "\n",
       "        function geo_json_c79781e85236c459198fca1d4a2a3ec0_add (data) {\n",
       "            geo_json_c79781e85236c459198fca1d4a2a3ec0\n",
       "                .addData(data);\n",
       "        }\n",
       "            geo_json_c79781e85236c459198fca1d4a2a3ec0_add({&quot;bbox&quot;: [176.9337392505123, -39.33855011832538, 176.93929770227703, -39.33471823699008], &quot;features&quot;: [{&quot;bbox&quot;: [176.9337392505123, -39.33855011832538, 176.93929770227703, -39.33471823699008], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[176.93780118519558, -39.33471823699008], [176.9337392505123, -39.33670605529121], [176.93523571183061, -39.33855011832538], [176.93929770227703, -39.33656224488286], [176.93780118519558, -39.33471823699008]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;0&quot;, &quot;properties&quot;: {&quot;id&quot;: &quot;nzd0563&quot;}, &quot;type&quot;: &quot;Feature&quot;}], &quot;type&quot;: &quot;FeatureCollection&quot;});\n",
       "\n",
       "        \n",
       "    \n",
       "    geo_json_c79781e85236c459198fca1d4a2a3ec0.bindTooltip(\n",
       "    function(layer){\n",
       "    let div = L.DomUtil.create(&#x27;div&#x27;);\n",
       "    \n",
       "    let handleObject = feature=&gt;typeof(feature)==&#x27;object&#x27; ? JSON.stringify(feature) : feature;\n",
       "    let fields = [&quot;id&quot;];\n",
       "    let aliases = [&quot;id&quot;];\n",
       "    let table = &#x27;&lt;table&gt;&#x27; +\n",
       "        String(\n",
       "        fields.map(\n",
       "        (v,i)=&gt;\n",
       "        `&lt;tr&gt;\n",
       "            &lt;th&gt;${aliases[i]}&lt;/th&gt;\n",
       "            \n",
       "            &lt;td&gt;${handleObject(layer.feature.properties[v])}&lt;/td&gt;\n",
       "        &lt;/tr&gt;`).join(&#x27;&#x27;))\n",
       "    +&#x27;&lt;/table&gt;&#x27;;\n",
       "    div.innerHTML=table;\n",
       "    \n",
       "    return div\n",
       "    }\n",
       "    ,{&quot;className&quot;: &quot;foliumtooltip&quot;, &quot;sticky&quot;: true});\n",
       "                     \n",
       "    \n",
       "            geo_json_c79781e85236c459198fca1d4a2a3ec0.addTo(map_dbd833a77680d77f8a419b432d63bff6);\n",
       "        \n",
       "&lt;/script&gt;\n",
       "&lt;/html&gt;\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
      ],
      "text/plain": [
       "<folium.folium.Map at 0x7fd52ba2a380>"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "new_poly = gpd.GeoDataFrame([{\"id\": new_siteid}], geometry=new_shoreline.to_crs(2193).buffer(100).minimum_rotated_rectangle().to_crs(4326), crs=4326)\n",
    "new_poly.explore()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<pre>✔️ 124 ms (2024-12-02T16:22:26/2024-12-02T16:22:26)</pre>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "pd.concat((poly, new_poly)).to_file(\"polygons.geojson\", driver=\"GeoJSON\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<pre>✔️ 79.6 ms (2024-12-02T16:22:27/2024-12-02T16:22:27)</pre>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "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>id</th>\n",
       "      <th>site_id</th>\n",
       "      <th>orientation</th>\n",
       "      <th>along_dist</th>\n",
       "      <th>along_dist_norm</th>\n",
       "      <th>beach_slope</th>\n",
       "      <th>cil</th>\n",
       "      <th>ciu</th>\n",
       "      <th>trend</th>\n",
       "      <th>n_points</th>\n",
       "      <th>n_points_nonan</th>\n",
       "      <th>r2_score</th>\n",
       "      <th>mae</th>\n",
       "      <th>mse</th>\n",
       "      <th>rmse</th>\n",
       "      <th>intercept</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>95050</th>\n",
       "      <td>nzd0001-0000</td>\n",
       "      <td>nzd0001</td>\n",
       "      <td>359.037136</td>\n",
       "      <td>3197.737936</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>0.050</td>\n",
       "      <td>0.0397</td>\n",
       "      <td>0.0679</td>\n",
       "      <td>-0.136498</td>\n",
       "      <td>184.0</td>\n",
       "      <td>143.0</td>\n",
       "      <td>0.001809</td>\n",
       "      <td>19.590361</td>\n",
       "      <td>572.322498</td>\n",
       "      <td>23.923263</td>\n",
       "      <td>325.377492</td>\n",
       "      <td>LINESTRING (1596662.375 6190263.089, 1596646.6...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>95051</th>\n",
       "      <td>nzd0001-0001</td>\n",
       "      <td>nzd0001</td>\n",
       "      <td>359.037136</td>\n",
       "      <td>3097.737936</td>\n",
       "      <td>0.968728</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-0.316164</td>\n",
       "      <td>184.0</td>\n",
       "      <td>157.0</td>\n",
       "      <td>0.011511</td>\n",
       "      <td>17.845204</td>\n",
       "      <td>490.058427</td>\n",
       "      <td>22.137263</td>\n",
       "      <td>330.472572</td>\n",
       "      <td>LINESTRING (1596744.903 6190264.499, 1596729.1...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>95052</th>\n",
       "      <td>nzd0001-0002</td>\n",
       "      <td>nzd0001</td>\n",
       "      <td>359.037136</td>\n",
       "      <td>2997.737936</td>\n",
       "      <td>0.937456</td>\n",
       "      <td>0.060</td>\n",
       "      <td>0.0478</td>\n",
       "      <td>0.0776</td>\n",
       "      <td>-0.412752</td>\n",
       "      <td>184.0</td>\n",
       "      <td>160.0</td>\n",
       "      <td>0.024944</td>\n",
       "      <td>16.023083</td>\n",
       "      <td>374.900676</td>\n",
       "      <td>19.362352</td>\n",
       "      <td>336.564138</td>\n",
       "      <td>LINESTRING (1596827.431 6190265.908, 1596811.6...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>95053</th>\n",
       "      <td>nzd0001-0003</td>\n",
       "      <td>nzd0001</td>\n",
       "      <td>359.037136</td>\n",
       "      <td>2897.737936</td>\n",
       "      <td>0.906184</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-0.469545</td>\n",
       "      <td>184.0</td>\n",
       "      <td>161.0</td>\n",
       "      <td>0.040228</td>\n",
       "      <td>14.359386</td>\n",
       "      <td>299.451117</td>\n",
       "      <td>17.304656</td>\n",
       "      <td>345.186311</td>\n",
       "      <td>LINESTRING (1596909.959 6190267.317, 1596894.2...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>95054</th>\n",
       "      <td>nzd0001-0004</td>\n",
       "      <td>nzd0001</td>\n",
       "      <td>359.037136</td>\n",
       "      <td>2797.737936</td>\n",
       "      <td>0.874912</td>\n",
       "      <td>0.055</td>\n",
       "      <td>0.0450</td>\n",
       "      <td>0.0664</td>\n",
       "      <td>-0.531422</td>\n",
       "      <td>184.0</td>\n",
       "      <td>162.0</td>\n",
       "      <td>0.055976</td>\n",
       "      <td>13.451432</td>\n",
       "      <td>276.975375</td>\n",
       "      <td>16.642577</td>\n",
       "      <td>357.606234</td>\n",
       "      <td>LINESTRING (1596992.486 6190268.724, 1596976.7...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>127421</th>\n",
       "      <td>nzd0562-0014</td>\n",
       "      <td>nzd0562</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.100</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.262019</td>\n",
       "      <td>583.0</td>\n",
       "      <td>447.0</td>\n",
       "      <td>0.043431</td>\n",
       "      <td>6.676504</td>\n",
       "      <td>80.625540</td>\n",
       "      <td>8.979173</td>\n",
       "      <td>193.845690</td>\n",
       "      <td>LINESTRING (1770293.383 5917676.929, 1770648.7...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>127422</th>\n",
       "      <td>nzd0562-0015</td>\n",
       "      <td>nzd0562</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.100</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.158424</td>\n",
       "      <td>583.0</td>\n",
       "      <td>448.0</td>\n",
       "      <td>0.015499</td>\n",
       "      <td>6.952012</td>\n",
       "      <td>84.615902</td>\n",
       "      <td>9.198690</td>\n",
       "      <td>196.070192</td>\n",
       "      <td>LINESTRING (1770254.942 5917751.29, 1770610.27...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>127423</th>\n",
       "      <td>nzd0562-0016</td>\n",
       "      <td>nzd0562</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.100</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.130158</td>\n",
       "      <td>583.0</td>\n",
       "      <td>439.0</td>\n",
       "      <td>0.011180</td>\n",
       "      <td>6.628118</td>\n",
       "      <td>79.216921</td>\n",
       "      <td>8.900389</td>\n",
       "      <td>205.098465</td>\n",
       "      <td>LINESTRING (1770214.898 5917828.821, 1770573.4...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>127424</th>\n",
       "      <td>nzd0562-0017</td>\n",
       "      <td>nzd0562</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.100</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.005730</td>\n",
       "      <td>583.0</td>\n",
       "      <td>405.0</td>\n",
       "      <td>0.000014</td>\n",
       "      <td>8.118279</td>\n",
       "      <td>126.880731</td>\n",
       "      <td>11.264135</td>\n",
       "      <td>209.805326</td>\n",
       "      <td>LINESTRING (1770180.582 5917901.338, 1770542.2...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>127425</th>\n",
       "      <td>nzd0562-0018</td>\n",
       "      <td>nzd0562</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.100</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.424099</td>\n",
       "      <td>583.0</td>\n",
       "      <td>365.0</td>\n",
       "      <td>0.045716</td>\n",
       "      <td>9.717575</td>\n",
       "      <td>202.294073</td>\n",
       "      <td>14.223012</td>\n",
       "      <td>208.659953</td>\n",
       "      <td>LINESTRING (1770147.813 5917970.658, 1770509.4...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>32376 rows × 17 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                  id  site_id  orientation   along_dist  along_dist_norm  \\\n",
       "95050   nzd0001-0000  nzd0001   359.037136  3197.737936         1.000000   \n",
       "95051   nzd0001-0001  nzd0001   359.037136  3097.737936         0.968728   \n",
       "95052   nzd0001-0002  nzd0001   359.037136  2997.737936         0.937456   \n",
       "95053   nzd0001-0003  nzd0001   359.037136  2897.737936         0.906184   \n",
       "95054   nzd0001-0004  nzd0001   359.037136  2797.737936         0.874912   \n",
       "...              ...      ...          ...          ...              ...   \n",
       "127421  nzd0562-0014  nzd0562          NaN          NaN              NaN   \n",
       "127422  nzd0562-0015  nzd0562          NaN          NaN              NaN   \n",
       "127423  nzd0562-0016  nzd0562          NaN          NaN              NaN   \n",
       "127424  nzd0562-0017  nzd0562          NaN          NaN              NaN   \n",
       "127425  nzd0562-0018  nzd0562          NaN          NaN              NaN   \n",
       "\n",
       "        beach_slope     cil     ciu     trend  n_points  n_points_nonan  \\\n",
       "95050         0.050  0.0397  0.0679 -0.136498     184.0           143.0   \n",
       "95051           NaN     NaN     NaN -0.316164     184.0           157.0   \n",
       "95052         0.060  0.0478  0.0776 -0.412752     184.0           160.0   \n",
       "95053           NaN     NaN     NaN -0.469545     184.0           161.0   \n",
       "95054         0.055  0.0450  0.0664 -0.531422     184.0           162.0   \n",
       "...             ...     ...     ...       ...       ...             ...   \n",
       "127421        0.100     NaN     NaN  0.262019     583.0           447.0   \n",
       "127422        0.100     NaN     NaN  0.158424     583.0           448.0   \n",
       "127423        0.100     NaN     NaN  0.130158     583.0           439.0   \n",
       "127424        0.100     NaN     NaN  0.005730     583.0           405.0   \n",
       "127425        0.100     NaN     NaN  0.424099     583.0           365.0   \n",
       "\n",
       "        r2_score        mae         mse       rmse   intercept  \\\n",
       "95050   0.001809  19.590361  572.322498  23.923263  325.377492   \n",
       "95051   0.011511  17.845204  490.058427  22.137263  330.472572   \n",
       "95052   0.024944  16.023083  374.900676  19.362352  336.564138   \n",
       "95053   0.040228  14.359386  299.451117  17.304656  345.186311   \n",
       "95054   0.055976  13.451432  276.975375  16.642577  357.606234   \n",
       "...          ...        ...         ...        ...         ...   \n",
       "127421  0.043431   6.676504   80.625540   8.979173  193.845690   \n",
       "127422  0.015499   6.952012   84.615902   9.198690  196.070192   \n",
       "127423  0.011180   6.628118   79.216921   8.900389  205.098465   \n",
       "127424  0.000014   8.118279  126.880731  11.264135  209.805326   \n",
       "127425  0.045716   9.717575  202.294073  14.223012  208.659953   \n",
       "\n",
       "                                                 geometry  \n",
       "95050   LINESTRING (1596662.375 6190263.089, 1596646.6...  \n",
       "95051   LINESTRING (1596744.903 6190264.499, 1596729.1...  \n",
       "95052   LINESTRING (1596827.431 6190265.908, 1596811.6...  \n",
       "95053   LINESTRING (1596909.959 6190267.317, 1596894.2...  \n",
       "95054   LINESTRING (1596992.486 6190268.724, 1596976.7...  \n",
       "...                                                   ...  \n",
       "127421  LINESTRING (1770293.383 5917676.929, 1770648.7...  \n",
       "127422  LINESTRING (1770254.942 5917751.29, 1770610.27...  \n",
       "127423  LINESTRING (1770214.898 5917828.821, 1770573.4...  \n",
       "127424  LINESTRING (1770180.582 5917901.338, 1770542.2...  \n",
       "127425  LINESTRING (1770147.813 5917970.658, 1770509.4...  \n",
       "\n",
       "[32376 rows x 17 columns]"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "nz_transects = transects[transects.id.str.startswith(\"nzd\")].to_crs(2193)\n",
    "nz_transects"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<pre>✔️ 74.1 ms (2024-12-02T16:22:27/2024-12-02T16:22:27)</pre>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "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>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>LINESTRING (176.9374 -39.3345, 176.938 -39.33807)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>LINESTRING (176.93648 -39.33459, 176.93708 -39...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>LINESTRING (176.93431 -39.33569, 176.93792 -39...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry\n",
       "0  LINESTRING (176.9374 -39.3345, 176.938 -39.33807)\n",
       "1  LINESTRING (176.93648 -39.33459, 176.93708 -39...\n",
       "2  LINESTRING (176.93431 -39.33569, 176.93792 -39..."
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Make sure the origin is inland\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc=\"&lt;!DOCTYPE html&gt;\n",
       "&lt;html&gt;\n",
       "&lt;head&gt;\n",
       "    \n",
       "    &lt;meta http-equiv=&quot;content-type&quot; content=&quot;text/html; charset=UTF-8&quot; /&gt;\n",
       "    \n",
       "        &lt;script&gt;\n",
       "            L_NO_TOUCH = false;\n",
       "            L_DISABLE_3D = false;\n",
       "        &lt;/script&gt;\n",
       "    \n",
       "    &lt;style&gt;html, body {width: 100%;height: 100%;margin: 0;padding: 0;}&lt;/style&gt;\n",
       "    &lt;style&gt;#map {position:absolute;top:0;bottom:0;right:0;left:0;}&lt;/style&gt;\n",
       "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://code.jquery.com/jquery-3.7.1.min.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/js/bootstrap.bundle.min.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/css/bootstrap.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap-glyphicons.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6.2.0/css/all.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css&quot;/&gt;\n",
       "    \n",
       "            &lt;meta name=&quot;viewport&quot; content=&quot;width=device-width,\n",
       "                initial-scale=1.0, maximum-scale=1.0, user-scalable=no&quot; /&gt;\n",
       "            &lt;style&gt;\n",
       "                #map_753b1d50ae52750956566793d919183a {\n",
       "                    position: relative;\n",
       "                    width: 100.0%;\n",
       "                    height: 100.0%;\n",
       "                    left: 0.0%;\n",
       "                    top: 0.0%;\n",
       "                }\n",
       "                .leaflet-container { font-size: 1rem; }\n",
       "            &lt;/style&gt;\n",
       "        \n",
       "    \n",
       "                    &lt;style&gt;\n",
       "                        .foliumtooltip {\n",
       "                            \n",
       "                        }\n",
       "                       .foliumtooltip table{\n",
       "                            margin: auto;\n",
       "                        }\n",
       "                        .foliumtooltip tr{\n",
       "                            text-align: left;\n",
       "                        }\n",
       "                        .foliumtooltip th{\n",
       "                            padding: 2px; padding-right: 8px;\n",
       "                        }\n",
       "                    &lt;/style&gt;\n",
       "            \n",
       "    \n",
       "                    &lt;style&gt;\n",
       "                        .foliumtooltip {\n",
       "                            \n",
       "                        }\n",
       "                       .foliumtooltip table{\n",
       "                            margin: auto;\n",
       "                        }\n",
       "                        .foliumtooltip tr{\n",
       "                            text-align: left;\n",
       "                        }\n",
       "                        .foliumtooltip th{\n",
       "                            padding: 2px; padding-right: 8px;\n",
       "                        }\n",
       "                    &lt;/style&gt;\n",
       "            \n",
       "    \n",
       "                    &lt;style&gt;\n",
       "                        .foliumtooltip {\n",
       "                            \n",
       "                        }\n",
       "                       .foliumtooltip table{\n",
       "                            margin: auto;\n",
       "                        }\n",
       "                        .foliumtooltip tr{\n",
       "                            text-align: left;\n",
       "                        }\n",
       "                        .foliumtooltip th{\n",
       "                            padding: 2px; padding-right: 8px;\n",
       "                        }\n",
       "                    &lt;/style&gt;\n",
       "            \n",
       "    \n",
       "                    &lt;style&gt;\n",
       "                        .foliumtooltip {\n",
       "                            \n",
       "                        }\n",
       "                       .foliumtooltip table{\n",
       "                            margin: auto;\n",
       "                        }\n",
       "                        .foliumtooltip tr{\n",
       "                            text-align: left;\n",
       "                        }\n",
       "                        .foliumtooltip th{\n",
       "                            padding: 2px; padding-right: 8px;\n",
       "                        }\n",
       "                    &lt;/style&gt;\n",
       "            \n",
       "&lt;/head&gt;\n",
       "&lt;body&gt;\n",
       "    \n",
       "    \n",
       "            &lt;div class=&quot;folium-map&quot; id=&quot;map_753b1d50ae52750956566793d919183a&quot; &gt;&lt;/div&gt;\n",
       "        \n",
       "&lt;/body&gt;\n",
       "&lt;script&gt;\n",
       "    \n",
       "    \n",
       "            var map_753b1d50ae52750956566793d919183a = L.map(\n",
       "                &quot;map_753b1d50ae52750956566793d919183a&quot;,\n",
       "                {\n",
       "                    center: [-39.336327006965575, 176.93615190470987],\n",
       "                    crs: L.CRS.EPSG3857,\n",
       "                    zoom: 10,\n",
       "                    zoomControl: true,\n",
       "                    preferCanvas: false,\n",
       "                }\n",
       "            );\n",
       "            L.control.scale().addTo(map_753b1d50ae52750956566793d919183a);\n",
       "\n",
       "            \n",
       "\n",
       "        \n",
       "    \n",
       "            var tile_layer_c8e0f0fd51e222d6377ef3d680b26354 = L.tileLayer(\n",
       "                &quot;https://tile.openstreetmap.org/{z}/{x}/{y}.png&quot;,\n",
       "                {&quot;attribution&quot;: &quot;\\u0026copy; \\u003ca href=\\&quot;https://www.openstreetmap.org/copyright\\&quot;\\u003eOpenStreetMap\\u003c/a\\u003e contributors&quot;, &quot;detectRetina&quot;: false, &quot;maxNativeZoom&quot;: 19, &quot;maxZoom&quot;: 19, &quot;minZoom&quot;: 0, &quot;noWrap&quot;: false, &quot;opacity&quot;: 1, &quot;subdomains&quot;: &quot;abc&quot;, &quot;tms&quot;: false}\n",
       "            );\n",
       "        \n",
       "    \n",
       "            tile_layer_c8e0f0fd51e222d6377ef3d680b26354.addTo(map_753b1d50ae52750956566793d919183a);\n",
       "        \n",
       "    \n",
       "            map_753b1d50ae52750956566793d919183a.fitBounds(\n",
       "                [[-39.33815780276676, 176.93430595052536], [-39.33449621116439, 176.9379978588944]],\n",
       "                {}\n",
       "            );\n",
       "        \n",
       "    \n",
       "        function geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb_styler(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.5, &quot;weight&quot;: 2};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb_highlighter(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.75};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb_pointToLayer(feature, latlng) {\n",
       "            var opts = {&quot;bubblingMouseEvents&quot;: true, &quot;color&quot;: &quot;#3388ff&quot;, &quot;dashArray&quot;: null, &quot;dashOffset&quot;: null, &quot;fill&quot;: true, &quot;fillColor&quot;: &quot;#3388ff&quot;, &quot;fillOpacity&quot;: 0.2, &quot;fillRule&quot;: &quot;evenodd&quot;, &quot;lineCap&quot;: &quot;round&quot;, &quot;lineJoin&quot;: &quot;round&quot;, &quot;opacity&quot;: 1.0, &quot;radius&quot;: 2, &quot;stroke&quot;: true, &quot;weight&quot;: 3};\n",
       "            \n",
       "            let style = geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb_styler(feature)\n",
       "            Object.assign(opts, style)\n",
       "            \n",
       "            return new L.CircleMarker(latlng, opts)\n",
       "        }\n",
       "\n",
       "        function geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb_onEachFeature(feature, layer) {\n",
       "            layer.on({\n",
       "                mouseout: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                            geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb.resetStyle(e.target);\n",
       "                    }\n",
       "                },\n",
       "                mouseover: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                        const highlightStyle = geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb_highlighter(e.target.feature)\n",
       "                        e.target.setStyle(highlightStyle);\n",
       "                    }\n",
       "                },\n",
       "            });\n",
       "        };\n",
       "        var geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb = L.geoJson(null, {\n",
       "                onEachFeature: geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb_onEachFeature,\n",
       "            \n",
       "                style: geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb_styler,\n",
       "                pointToLayer: geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb_pointToLayer,\n",
       "        });\n",
       "\n",
       "        function geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb_add (data) {\n",
       "            geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb\n",
       "                .addData(data);\n",
       "        }\n",
       "            geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb_add({&quot;bbox&quot;: [176.93430595052536, -39.33815780276676, 176.9379978588944, -39.33449621116439], &quot;features&quot;: [{&quot;bbox&quot;: [176.9374041599855, -39.33806578689432, 176.9379978588944, -39.33449621116439], &quot;geometry&quot;: {&quot;coordinates&quot;: [[176.9374041599855, -39.33449621116439], [176.9379978588944, -39.33806578689432]], &quot;type&quot;: &quot;LineString&quot;}, &quot;id&quot;: &quot;0&quot;, &quot;properties&quot;: {}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [176.93648477096534, -39.33815780276676, 176.93707851525744, -39.33458822915117], &quot;geometry&quot;: {&quot;coordinates&quot;: [[176.93648477096534, -39.33458822915117], [176.93707851525744, -39.33815780276676]], &quot;type&quot;: &quot;LineString&quot;}, &quot;id&quot;: &quot;1&quot;, &quot;properties&quot;: {}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [176.93430595052536, -39.33793874657172, 176.93791992888504, -39.335685225885015], &quot;geometry&quot;: {&quot;coordinates&quot;: [[176.93430595052536, -39.335685225885015], [176.93791992888504, -39.33793874657172]], &quot;type&quot;: &quot;LineString&quot;}, &quot;id&quot;: &quot;2&quot;, &quot;properties&quot;: {}, &quot;type&quot;: &quot;Feature&quot;}], &quot;type&quot;: &quot;FeatureCollection&quot;});\n",
       "\n",
       "        \n",
       "    \n",
       "    geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb.bindTooltip(\n",
       "    function(layer){\n",
       "    let div = L.DomUtil.create(&#x27;div&#x27;);\n",
       "    \n",
       "    return div\n",
       "    }\n",
       "    ,{&quot;className&quot;: &quot;foliumtooltip&quot;, &quot;sticky&quot;: true});\n",
       "                     \n",
       "    \n",
       "            geo_json_f76eb587e5c4a7ab3fdde76bd87b98fb.addTo(map_753b1d50ae52750956566793d919183a);\n",
       "        \n",
       "    \n",
       "        function geo_json_2ebc09d5989928c238eb2ebfdd350424_styler(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.5, &quot;weight&quot;: 2};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_2ebc09d5989928c238eb2ebfdd350424_highlighter(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.75};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_2ebc09d5989928c238eb2ebfdd350424_pointToLayer(feature, latlng) {\n",
       "            var opts = {&quot;bubblingMouseEvents&quot;: true, &quot;color&quot;: &quot;#3388ff&quot;, &quot;dashArray&quot;: null, &quot;dashOffset&quot;: null, &quot;fill&quot;: true, &quot;fillColor&quot;: &quot;#3388ff&quot;, &quot;fillOpacity&quot;: 0.2, &quot;fillRule&quot;: &quot;evenodd&quot;, &quot;lineCap&quot;: &quot;round&quot;, &quot;lineJoin&quot;: &quot;round&quot;, &quot;opacity&quot;: 1.0, &quot;radius&quot;: 2, &quot;stroke&quot;: true, &quot;weight&quot;: 3};\n",
       "            \n",
       "            let style = geo_json_2ebc09d5989928c238eb2ebfdd350424_styler(feature)\n",
       "            Object.assign(opts, style)\n",
       "            \n",
       "            return new L.CircleMarker(latlng, opts)\n",
       "        }\n",
       "\n",
       "        function geo_json_2ebc09d5989928c238eb2ebfdd350424_onEachFeature(feature, layer) {\n",
       "            layer.on({\n",
       "                mouseout: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                            geo_json_2ebc09d5989928c238eb2ebfdd350424.resetStyle(e.target);\n",
       "                    }\n",
       "                },\n",
       "                mouseover: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                        const highlightStyle = geo_json_2ebc09d5989928c238eb2ebfdd350424_highlighter(e.target.feature)\n",
       "                        e.target.setStyle(highlightStyle);\n",
       "                    }\n",
       "                },\n",
       "            });\n",
       "        };\n",
       "        var geo_json_2ebc09d5989928c238eb2ebfdd350424 = L.geoJson(null, {\n",
       "                onEachFeature: geo_json_2ebc09d5989928c238eb2ebfdd350424_onEachFeature,\n",
       "            \n",
       "                style: geo_json_2ebc09d5989928c238eb2ebfdd350424_styler,\n",
       "                pointToLayer: geo_json_2ebc09d5989928c238eb2ebfdd350424_pointToLayer,\n",
       "        });\n",
       "\n",
       "        function geo_json_2ebc09d5989928c238eb2ebfdd350424_add (data) {\n",
       "            geo_json_2ebc09d5989928c238eb2ebfdd350424\n",
       "                .addData(data);\n",
       "        }\n",
       "            geo_json_2ebc09d5989928c238eb2ebfdd350424_add({&quot;bbox&quot;: [176.935598, -39.33731, 176.937701, -39.336281], &quot;features&quot;: [{&quot;bbox&quot;: [176.935598, -39.33731, 176.937701, -39.336281], &quot;geometry&quot;: {&quot;coordinates&quot;: [[176.937701, -39.336281], [176.936542, -39.336397], [176.935598, -39.33731]], &quot;type&quot;: &quot;LineString&quot;}, &quot;id&quot;: &quot;0&quot;, &quot;properties&quot;: {&quot;id&quot;: &quot;nzd0563&quot;}, &quot;type&quot;: &quot;Feature&quot;}], &quot;type&quot;: &quot;FeatureCollection&quot;});\n",
       "\n",
       "        \n",
       "    \n",
       "    geo_json_2ebc09d5989928c238eb2ebfdd350424.bindTooltip(\n",
       "    function(layer){\n",
       "    let div = L.DomUtil.create(&#x27;div&#x27;);\n",
       "    \n",
       "    let handleObject = feature=&gt;typeof(feature)==&#x27;object&#x27; ? JSON.stringify(feature) : feature;\n",
       "    let fields = [&quot;id&quot;];\n",
       "    let aliases = [&quot;id&quot;];\n",
       "    let table = &#x27;&lt;table&gt;&#x27; +\n",
       "        String(\n",
       "        fields.map(\n",
       "        (v,i)=&gt;\n",
       "        `&lt;tr&gt;\n",
       "            &lt;th&gt;${aliases[i]}&lt;/th&gt;\n",
       "            \n",
       "            &lt;td&gt;${handleObject(layer.feature.properties[v])}&lt;/td&gt;\n",
       "        &lt;/tr&gt;`).join(&#x27;&#x27;))\n",
       "    +&#x27;&lt;/table&gt;&#x27;;\n",
       "    div.innerHTML=table;\n",
       "    \n",
       "    return div\n",
       "    }\n",
       "    ,{&quot;className&quot;: &quot;foliumtooltip&quot;, &quot;sticky&quot;: true});\n",
       "                     \n",
       "    \n",
       "            geo_json_2ebc09d5989928c238eb2ebfdd350424.addTo(map_753b1d50ae52750956566793d919183a);\n",
       "        \n",
       "    \n",
       "        function geo_json_65ff5a27f20413935a299a3c4bc35b0d_styler(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.5, &quot;weight&quot;: 2};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_65ff5a27f20413935a299a3c4bc35b0d_highlighter(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.75};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_65ff5a27f20413935a299a3c4bc35b0d_pointToLayer(feature, latlng) {\n",
       "            var opts = {&quot;bubblingMouseEvents&quot;: true, &quot;color&quot;: &quot;#3388ff&quot;, &quot;dashArray&quot;: null, &quot;dashOffset&quot;: null, &quot;fill&quot;: true, &quot;fillColor&quot;: &quot;#3388ff&quot;, &quot;fillOpacity&quot;: 0.2, &quot;fillRule&quot;: &quot;evenodd&quot;, &quot;lineCap&quot;: &quot;round&quot;, &quot;lineJoin&quot;: &quot;round&quot;, &quot;opacity&quot;: 1.0, &quot;radius&quot;: 2, &quot;stroke&quot;: true, &quot;weight&quot;: 3};\n",
       "            \n",
       "            let style = geo_json_65ff5a27f20413935a299a3c4bc35b0d_styler(feature)\n",
       "            Object.assign(opts, style)\n",
       "            \n",
       "            return new L.CircleMarker(latlng, opts)\n",
       "        }\n",
       "\n",
       "        function geo_json_65ff5a27f20413935a299a3c4bc35b0d_onEachFeature(feature, layer) {\n",
       "            layer.on({\n",
       "                mouseout: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                            geo_json_65ff5a27f20413935a299a3c4bc35b0d.resetStyle(e.target);\n",
       "                    }\n",
       "                },\n",
       "                mouseover: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                        const highlightStyle = geo_json_65ff5a27f20413935a299a3c4bc35b0d_highlighter(e.target.feature)\n",
       "                        e.target.setStyle(highlightStyle);\n",
       "                    }\n",
       "                },\n",
       "            });\n",
       "        };\n",
       "        var geo_json_65ff5a27f20413935a299a3c4bc35b0d = L.geoJson(null, {\n",
       "                onEachFeature: geo_json_65ff5a27f20413935a299a3c4bc35b0d_onEachFeature,\n",
       "            \n",
       "                style: geo_json_65ff5a27f20413935a299a3c4bc35b0d_styler,\n",
       "                pointToLayer: geo_json_65ff5a27f20413935a299a3c4bc35b0d_pointToLayer,\n",
       "        });\n",
       "\n",
       "        function geo_json_65ff5a27f20413935a299a3c4bc35b0d_add (data) {\n",
       "            geo_json_65ff5a27f20413935a299a3c4bc35b0d\n",
       "                .addData(data);\n",
       "        }\n",
       "            geo_json_65ff5a27f20413935a299a3c4bc35b0d_add({&quot;bbox&quot;: [176.9337392505123, -39.33855011832538, 176.93929770227703, -39.33471823699008], &quot;features&quot;: [{&quot;bbox&quot;: [176.9337392505123, -39.33855011832538, 176.93929770227703, -39.33471823699008], &quot;geometry&quot;: {&quot;coordinates&quot;: [[176.93780118519558, -39.33471823699008], [176.9337392505123, -39.33670605529121], [176.93523571183061, -39.33855011832538], [176.93929770227703, -39.33656224488286], [176.93780118519558, -39.33471823699008]], &quot;type&quot;: &quot;LineString&quot;}, &quot;id&quot;: &quot;0&quot;, &quot;properties&quot;: {}, &quot;type&quot;: &quot;Feature&quot;}], &quot;type&quot;: &quot;FeatureCollection&quot;});\n",
       "\n",
       "        \n",
       "    \n",
       "    geo_json_65ff5a27f20413935a299a3c4bc35b0d.bindTooltip(\n",
       "    function(layer){\n",
       "    let div = L.DomUtil.create(&#x27;div&#x27;);\n",
       "    \n",
       "    return div\n",
       "    }\n",
       "    ,{&quot;className&quot;: &quot;foliumtooltip&quot;, &quot;sticky&quot;: true});\n",
       "                     \n",
       "    \n",
       "            geo_json_65ff5a27f20413935a299a3c4bc35b0d.addTo(map_753b1d50ae52750956566793d919183a);\n",
       "        \n",
       "    \n",
       "        function geo_json_8307a5ea85f51539e0d7dc8100084c82_styler(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;color&quot;: &quot;red&quot;, &quot;fillColor&quot;: &quot;red&quot;, &quot;fillOpacity&quot;: 0.5, &quot;weight&quot;: 2};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_8307a5ea85f51539e0d7dc8100084c82_highlighter(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.75};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_8307a5ea85f51539e0d7dc8100084c82_pointToLayer(feature, latlng) {\n",
       "            var opts = {&quot;bubblingMouseEvents&quot;: true, &quot;color&quot;: &quot;#3388ff&quot;, &quot;dashArray&quot;: null, &quot;dashOffset&quot;: null, &quot;fill&quot;: true, &quot;fillColor&quot;: &quot;#3388ff&quot;, &quot;fillOpacity&quot;: 0.2, &quot;fillRule&quot;: &quot;evenodd&quot;, &quot;lineCap&quot;: &quot;round&quot;, &quot;lineJoin&quot;: &quot;round&quot;, &quot;opacity&quot;: 1.0, &quot;radius&quot;: 2, &quot;stroke&quot;: true, &quot;weight&quot;: 3};\n",
       "            \n",
       "            let style = geo_json_8307a5ea85f51539e0d7dc8100084c82_styler(feature)\n",
       "            Object.assign(opts, style)\n",
       "            \n",
       "            return new L.CircleMarker(latlng, opts)\n",
       "        }\n",
       "\n",
       "        function geo_json_8307a5ea85f51539e0d7dc8100084c82_onEachFeature(feature, layer) {\n",
       "            layer.on({\n",
       "                mouseout: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                            geo_json_8307a5ea85f51539e0d7dc8100084c82.resetStyle(e.target);\n",
       "                    }\n",
       "                },\n",
       "                mouseover: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                        const highlightStyle = geo_json_8307a5ea85f51539e0d7dc8100084c82_highlighter(e.target.feature)\n",
       "                        e.target.setStyle(highlightStyle);\n",
       "                    }\n",
       "                },\n",
       "            });\n",
       "        };\n",
       "        var geo_json_8307a5ea85f51539e0d7dc8100084c82 = L.geoJson(null, {\n",
       "                onEachFeature: geo_json_8307a5ea85f51539e0d7dc8100084c82_onEachFeature,\n",
       "            \n",
       "                style: geo_json_8307a5ea85f51539e0d7dc8100084c82_styler,\n",
       "                pointToLayer: geo_json_8307a5ea85f51539e0d7dc8100084c82_pointToLayer,\n",
       "        });\n",
       "\n",
       "        function geo_json_8307a5ea85f51539e0d7dc8100084c82_add (data) {\n",
       "            geo_json_8307a5ea85f51539e0d7dc8100084c82\n",
       "                .addData(data);\n",
       "        }\n",
       "            geo_json_8307a5ea85f51539e0d7dc8100084c82_add({&quot;bbox&quot;: [176.93430595052536, -39.335685225885015, 176.9374041599855, -39.33449621116439], &quot;features&quot;: [{&quot;bbox&quot;: [176.9374041599855, -39.33449621116439, 176.9374041599855, -39.33449621116439], &quot;geometry&quot;: {&quot;coordinates&quot;: [176.9374041599855, -39.33449621116439], &quot;type&quot;: &quot;Point&quot;}, &quot;id&quot;: &quot;0&quot;, &quot;properties&quot;: {&quot;__folium_color&quot;: &quot;red&quot;}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [176.93648477096534, -39.33458822915117, 176.93648477096534, -39.33458822915117], &quot;geometry&quot;: {&quot;coordinates&quot;: [176.93648477096534, -39.33458822915117], &quot;type&quot;: &quot;Point&quot;}, &quot;id&quot;: &quot;1&quot;, &quot;properties&quot;: {&quot;__folium_color&quot;: &quot;red&quot;}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [176.93430595052536, -39.335685225885015, 176.93430595052536, -39.335685225885015], &quot;geometry&quot;: {&quot;coordinates&quot;: [176.93430595052536, -39.335685225885015], &quot;type&quot;: &quot;Point&quot;}, &quot;id&quot;: &quot;2&quot;, &quot;properties&quot;: {&quot;__folium_color&quot;: &quot;red&quot;}, &quot;type&quot;: &quot;Feature&quot;}], &quot;type&quot;: &quot;FeatureCollection&quot;});\n",
       "\n",
       "        \n",
       "    \n",
       "    geo_json_8307a5ea85f51539e0d7dc8100084c82.bindTooltip(\n",
       "    function(layer){\n",
       "    let div = L.DomUtil.create(&#x27;div&#x27;);\n",
       "    \n",
       "    return div\n",
       "    }\n",
       "    ,{&quot;className&quot;: &quot;foliumtooltip&quot;, &quot;sticky&quot;: true});\n",
       "                     \n",
       "    \n",
       "            geo_json_8307a5ea85f51539e0d7dc8100084c82.addTo(map_753b1d50ae52750956566793d919183a);\n",
       "        \n",
       "&lt;/script&gt;\n",
       "&lt;/html&gt;\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
      ],
      "text/plain": [
       "<folium.folium.Map at 0x7fd52ba295a0>"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def create_transects(line, spacing=80, transect_length=400):\n",
    "    transects = []\n",
    "    distances = np.arange(0, line.length, spacing)\n",
    "\n",
    "    for distance in distances:\n",
    "        # Find point along the average line\n",
    "        point = line.interpolate(distance)\n",
    "        \n",
    "        # Find the direction of the line at this point (tangent direction)\n",
    "        nearest_point_ahead = line.interpolate(min(distance + 1e-6, line.length))\n",
    "        direction = np.arctan2(nearest_point_ahead.y - point.y, nearest_point_ahead.x - point.x)\n",
    "        \n",
    "        # Rotate 90 degrees (perpendicular) and extend to create a transect\n",
    "        transect = LineString([\n",
    "            Point(\n",
    "                point.x - transect_length / 2 * np.cos(direction + np.pi / 2),\n",
    "                point.y - transect_length / 2 * np.sin(direction + np.pi / 2)\n",
    "            ),\n",
    "            Point(\n",
    "                point.x + transect_length / 2 * np.cos(direction + np.pi / 2),\n",
    "                point.y + transect_length / 2 * np.sin(direction + np.pi / 2)\n",
    "            )\n",
    "        ])\n",
    "        \n",
    "        transects.append(transect)\n",
    "    \n",
    "    transects = gpd.GeoDataFrame(geometry=transects, crs=2193).to_crs(4326)\n",
    "    return transects\n",
    "\n",
    "new_transects = create_transects(new_shoreline.to_crs(2193).geometry.iloc[0])\n",
    "display(new_transects)\n",
    "m = new_transects.explore()\n",
    "new_shoreline.explore(m=m)\n",
    "new_poly.boundary.explore(m=m)\n",
    "gpd.GeoSeries(new_transects.geometry.apply(lambda line: Point(line.coords[0])), crs=new_transects.crs).explore(m=m, color=\"red\", name=\"transect start\")\n",
    "print(\"Make sure the origin is inland\")\n",
    "m"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<pre>✔️ 11.5 ms (2024-12-02T16:22:27/2024-12-02T16:22:27)</pre>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "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>geometry</th>\n",
       "      <th>id</th>\n",
       "      <th>site_id</th>\n",
       "      <th>beach_slope</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>LINESTRING (176.9374 -39.3345, 176.938 -39.33807)</td>\n",
       "      <td>nzd0563-0000</td>\n",
       "      <td>nzd0563</td>\n",
       "      <td>0.1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>LINESTRING (176.93648 -39.33459, 176.93708 -39...</td>\n",
       "      <td>nzd0563-0001</td>\n",
       "      <td>nzd0563</td>\n",
       "      <td>0.1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>LINESTRING (176.93431 -39.33569, 176.93792 -39...</td>\n",
       "      <td>nzd0563-0002</td>\n",
       "      <td>nzd0563</td>\n",
       "      <td>0.1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry            id  site_id  \\\n",
       "0  LINESTRING (176.9374 -39.3345, 176.938 -39.33807)  nzd0563-0000  nzd0563   \n",
       "1  LINESTRING (176.93648 -39.33459, 176.93708 -39...  nzd0563-0001  nzd0563   \n",
       "2  LINESTRING (176.93431 -39.33569, 176.93792 -39...  nzd0563-0002  nzd0563   \n",
       "\n",
       "   beach_slope  \n",
       "0          0.1  \n",
       "1          0.1  \n",
       "2          0.1  "
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "new_transects[\"id\"] = new_siteid + \"-\" + new_transects.index.astype(str).str.pad(4, fillchar=\"0\")\n",
    "new_transects[\"site_id\"] = new_siteid\n",
    "#new_transects[\"beach_slope\"] = .1\n",
    "new_transects"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<pre>✔️ 4.33 s (2024-12-02T16:22:27/2024-12-02T16:22:31)</pre>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "pd.concat((transects, new_transects)).to_file(\"transects_extended.geojson\", driver=\"GeoJSON\")"
   ]
  }
 ],
 "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.10.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}