#4 imported the codes from kml

Author: cicuser

getDisplacement/getDisplacement.ipynb

@@ -83,16 +83,43 @@
     "    icon_url = \"https://maps.google.com/mapfiles/kml/paddle/grn-blank.png\"\n",
     "    html = '<a target=\"_blank\" href=\"https://sideshow.jpl.nasa.gov/post/links/{site}.html\"><img src=\"https://sideshow.jpl.nasa.gov/post/plots/{site}.jpg\" width=\"300\" height=\"300\"></a>'\n",
     "    for cord, disp, site in zip(zip(df[lat_col], df[lon_col]), df[data_col], df[site_col]):\n",
-    "        folium.Circle(\n",
-    "            location=cord,\n",
-    "            radius=abs(disp)*scale,\n",
+    "        theta = 0\n",
+    "        coords = []\n",
+    "        for k in range(0,31):\n",
+    "            vrad = disp\n",
+    "            lat = cord[0]\n",
+    "            lon = cord[1]\n",
+    "            angle = k/30*2*math.pi\n",
+    "            elon = vrad*math.cos(angle)*math.cos(theta)-vrad*math.sin(angle)*math.sin(theta)\n",
+    "            elat = vrad*math.cos(angle)*math.sin(theta)+vrad*math.sin(angle)*math.cos(theta)\n",
+    "            elon = (elon+0)/scale/math.cos(lat*math.pi/180.)\n",
+    "            elat = (elat+0)/scale\n",
+    "            coords.append([lat + elat,lon + elon])\n",
+    "        if (disp > 0):\n",
+    "            lcolor = 'red'\n",
+    "        else:\n",
+    "            lcolor = 'blue'\n",
+    "        folium.Polygon(\n",
+    "            locations=coords,\n",
+    "            fillColor=lcolor,\n",
     "            fill=True,\n",
-    "            color=colormap(disp),\n",
+    "            fillOpacity=0.9,\n",
+    "            color=lcolor,\n",
+    "            weight=2,\n",
     "            popup=f'Vertical Displacement:{\"{:10.4f}\".format(disp)} mm,   Site:{site}',\n",
-    "            stroke=True,\n",
-    "            weight=1,\n",
-    "            fill_opacity=.8\n",
-    "        ).add_to(vertical)\n",
+    "            ).add_to(vertical)\n",
+    "#         folium.Circle(\n",
+    "#             location=cord,\n",
+    "#             radius=abs(disp)*scale,\n",
+    "#             fill=True,\n",
+    "#             color=colormap(disp),\n",
+    "#             popup=f'Vertical Displacement:{\"{:10.4f}\".format(disp)} mm,   Site:{site}',\n",
+    "#             stroke=True,\n",
+    "#             weight=1,\n",
+    "#             fill_opacity=.8\n",
+    "#         ).add_to(vertical)\n",
+    "        \n",
+    "        \n",
     "        icon = folium.features.CustomIcon(icon_url,icon_size=(16, 16),icon_anchor=(8,16))\n",
     "        #iframe = folium.element.IFrame(html=html, width=500, height=300)\n",
     "        s_html = html.format(site=site)\n",
@@ -342,7 +369,7 @@
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "72fbff5266044cc5819c63cf92614446",
+       "model_id": "8254138517d549418c2db180ac830fa1",
        "version_major": 2,
        "version_minor": 0
       },
@@ -447,7 +474,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.9.2"
+   "version": "3.8.2"
   }
  },
  "nbformat": 4,