demo_jharkhand.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Overview\n",
    "\n",
    "Randomly select 10 villages / wards in Ranchi and Khunti districts of Jharkhand and then select 3 points in each village / ward.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import geopandas as gpd\n",
    "from pathlib import Path\n",
    "import pandas as pd\n",
    "import folium\n",
    "from tqdm import tqdm\n",
    "from shapely.geometry import LineString\n",
    "tqdm.pandas()\n",
    "from shapely import Point\n",
    "from pin_drop_sampling2.utils import get_s2_cell_id, count_neighbors_in_radius, get_nearest_point_on_road, dist_in_meters\n",
    "\n",
    "# set parameter for the number of PSUs to sample\n",
    "NUM_PSUS_TO_SAMPLE = 10"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Set file locations"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "DB_DIR = Path.home() / 'IDinsight Dropbox' / 'Random Walk Testing' \n",
    "RAW_DATA = Path.home() / 'IDinsight Dropbox'/'Random Walk Testing'/'01_Raw data'\n",
    "SHRUG_DIR = RAW_DATA / '02_Admin boundary data' / 'SHRUG' \n",
    "ROOFTOP_DIR = RAW_DATA / '01_Rooftop'/'India'\n",
    "OUTPUT_DIR = DB_DIR / '03_Output' / '08_Jharkhand test pins'"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Import and clean files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# shrug\n",
    "shrug = gpd.read_file(SHRUG_DIR / 'shrug-pc11-village-poly-shp' / 'village.shp')\n",
    "shrug = shrug.rename(columns={'pc11_s_id': 'state_id', 'pc11_d_id': 'district_id', 'pc11_sd_id': 'subdistrict_id', 'pc11_tv_id': 'village_id'})\n",
    "shrug['shrid2'] = '11-' + shrug['state_id'].astype(str) + '-' + shrug['district_id'].astype(str) + '-' + shrug['subdistrict_id'].astype(str) + '-' + shrug['village_id'].astype(str)\n",
    "shrug = shrug.astype({'state_id': int, 'district_id': int, 'subdistrict_id': int, 'village_id': int})\n",
    "\n",
    "# shrug district names and ids\n",
    "shrug_districts = pd.read_csv(SHRUG_DIR / 'shrug-shrid-keys-csv' / 'shrug_district_names_ids.csv')\n",
    "\n",
    "# shrug_loc_names \n",
    "shrug_loc_names = pd.read_csv(SHRUG_DIR / 'shrug-shrid-keys-csv' /  'shrid_loc_names.csv')\n",
    "shrug_loc_names['urban'] = shrug_loc_names['town_name'].notnull().astype(int)\n",
    "shrug_loc_names = shrug_loc_names[['shrid2', 'urban']]\n",
    "\n",
    "# census\n",
    "census = pd.read_csv(SHRUG_DIR/'shrug-pca11-csv'/'pc11_pca_clean_shrid.csv', usecols=['shrid2', 'pc11_pca_no_hh'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Filter SHRUG data for Ranchi and Khunti and merge with census data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "length of ranchi khunti: 2084\n",
      "Num shrug rows with no census data: 24\n",
      "Num shrug rows with no data on urban status: 24\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "urban\n",
       "0.0    0.625107\n",
       "1.0    0.374893\n",
       "Name: pc11_pca_no_hh, dtype: float64"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# filter shrug data for Ranchi and Khunti districts\n",
    "ranchi_khunti_districts =  shrug_districts[shrug_districts['district_name'].isin(['ranchi', 'khunti'])]\n",
    "ranchi_khunti = shrug.merge(ranchi_khunti_districts, on=['state_id', 'district_id'])\n",
    "print(f\"length of ranchi khunti: {len(ranchi_khunti)}\")\n",
    "\n",
    "# left merge census and ranchi_khunti on shrid2\n",
    "ranchi_khunti_census = ranchi_khunti.merge(census, on='shrid2', how='left')\n",
    "\n",
    "# check how many shrug rows couldn't be merged with the census data\n",
    "print(f\"Num shrug rows with no census data: {int(ranchi_khunti_census['pc11_pca_no_hh'].isna().sum())}\")\n",
    "\n",
    "# merge with shrug_loc_names on shrid2 to get urban status of the sampled villages / wards\n",
    "ranchi_khunti_census = ranchi_khunti_census.merge(shrug_loc_names, on='shrid2', how='left')\n",
    "\n",
    "# check how many shrug rows couldn't be merged with the shrug_loc_names\n",
    "print(f\"Num shrug rows with no data on urban status: {int(ranchi_khunti_census['urban'].isna().sum())}\")\n",
    "\n",
    "# check the share of the population urban \n",
    "ranchi_khunti_census.groupby('urban')['pc11_pca_no_hh'].sum() / ranchi_khunti_census['pc11_pca_no_hh'].sum()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/var/folders/_b/2wf3gjnx42s4k_vs5n079mh40000gn/T/ipykernel_57522/676096483.py:6: UserWarning: Geometry is in a geographic CRS. Results from 'centroid' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.\n",
      "\n",
      "  center = gdf.centroid.iloc[0]\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_3c62ec64be52148a1a49e86e22e91445 {\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_3c62ec64be52148a1a49e86e22e91445&quot; &gt;&lt;/div&gt;\n",
       "        \n",
       "&lt;/body&gt;\n",
       "&lt;script&gt;\n",
       "    \n",
       "    \n",
       "            var map_3c62ec64be52148a1a49e86e22e91445 = L.map(\n",
       "                &quot;map_3c62ec64be52148a1a49e86e22e91445&quot;,\n",
       "                {\n",
       "                    center: [23.37148841470954, 85.2588506831631],\n",
       "                    crs: L.CRS.EPSG3857,\n",
       "                    zoom: 12,\n",
       "                    zoomControl: true,\n",
       "                    preferCanvas: false,\n",
       "                }\n",
       "            );\n",
       "\n",
       "            \n",
       "\n",
       "        \n",
       "    \n",
       "            var tile_layer_af219bd0ba39b1f1446ebb5ea1993d0e = 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_af219bd0ba39b1f1446ebb5ea1993d0e.addTo(map_3c62ec64be52148a1a49e86e22e91445);\n",
       "        \n",
       "    \n",
       "\n",
       "        function geo_json_5c1eb2c2f07c487a3f7d86e2aed9b9e3_onEachFeature(feature, layer) {\n",
       "            layer.on({\n",
       "            });\n",
       "        };\n",
       "        var geo_json_5c1eb2c2f07c487a3f7d86e2aed9b9e3 = L.geoJson(null, {\n",
       "                onEachFeature: geo_json_5c1eb2c2f07c487a3f7d86e2aed9b9e3_onEachFeature,\n",
       "            \n",
       "        });\n",
       "\n",
       "        function geo_json_5c1eb2c2f07c487a3f7d86e2aed9b9e3_add (data) {\n",
       "            geo_json_5c1eb2c2f07c487a3f7d86e2aed9b9e3\n",
       "                .addData(data);\n",
       "        }\n",
       "            geo_json_5c1eb2c2f07c487a3f7d86e2aed9b9e3_add({&quot;bbox&quot;: [84.98376454300008, 22.79738033400002, 85.8772120930001, 23.544132385000076], &quot;features&quot;: [{&quot;bbox&quot;: [85.25608253600011, 23.365590747000013, 85.26217985600002, 23.377069695428602], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.26217985600002, 23.376179139], [85.26188487300004, 23.375371636000068], [85.261023546, 23.37326744500005], [85.26012727700001, 23.371104414000047], [85.25982220300008, 23.37020477800004], [85.25976626900007, 23.369798579000093], [85.25974958600003, 23.368940928000008], [85.25976205600008, 23.368583174999983], [85.26006977300005, 23.367612047000012], [85.26067659900008, 23.36597346100006], [85.26029516100004, 23.36598919100003], [85.25846834600006, 23.366066075999996], [85.25781534400004, 23.365999153000054], [85.25608253600011, 23.365590747000013], [85.25609016800011, 23.366260370000084], [85.25618905700013, 23.36669543400002], [85.25638121000013, 23.36728808000003], [85.25706958400006, 23.369142489000076], [85.25723460700003, 23.36982881500002], [85.25735488900004, 23.370658407000093], [85.25733432400003, 23.371413043000093], [85.2572385110001, 23.37375646100001], [85.25720930700007, 23.37573712100002], [85.25726790000003, 23.37602311700006], [85.25749592400007, 23.376359322000077], [85.25776466400002, 23.37662653600007], [85.25830512957151, 23.377069695428602], [85.25909906400003, 23.37655521500003], [85.25939071000009, 23.376501004000076], [85.26029283800003, 23.376380169000115], [85.26217985600002, 23.376179139]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;162&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02684-373529&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2684, &quot;tv_name&quot;: &quot;Bajra&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 373529}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.25806278400007, 23.376179139, 85.26408852400003, 23.37831517500008], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.26408852400003, 23.377991184000095], [85.26312296850011, 23.377357754], [85.26283002400008, 23.377086720000044], [85.26258685300003, 23.376808498000027], [85.26241215800006, 23.376541201000027], [85.26217985600002, 23.376179139], [85.26029283800003, 23.376380169000115], [85.25939071000009, 23.376501004000076], [85.25909906400003, 23.37655521500003], [85.25830512957151, 23.377069695428602], [85.25830856450013, 23.377359556500103], [85.25806278400007, 23.377753842000086], [85.2585241815001, 23.37801068300011], [85.25973198700008, 23.378258560000102], [85.26012429700003, 23.37829932800001], [85.26197287700003, 23.37831517500008], [85.26245957200007, 23.378286699000032], [85.26408852400003, 23.377991184000095]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;163&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02684-373530&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2684, &quot;tv_name&quot;: &quot;Bariatu&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 373530}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.38700112000004, 23.403867955000106, 85.4024487240001, 23.419818337000013], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.4024487240001, 23.412652749000078], [85.40171898000003, 23.41187406100004], [85.40097534950009, 23.411042939000108], [85.40065622900002, 23.410801969000065], [85.40027720050006, 23.410538046500033], [85.39980978100006, 23.410320854000048], [85.399315297, 23.410183922000087], [85.39813490400003, 23.409984505000107], [85.39787352166675, 23.4097750850001], [85.39781844200006, 23.409378575000066], [85.39761902550012, 23.40785583650006], [85.39736983033339, 23.407696917333435], [85.39696893300003, 23.407690608000024], [85.39510891600004, 23.407680931500067], [85.39481769350006, 23.407659742000046], [85.39451572400007, 23.40744905300011], [85.3932513720001, 23.40449837300004], [85.39288509375014, 23.40456518175013], [85.39247960450005, 23.40466543125012], [85.391484394, 23.403867955000106], [85.39152871600002, 23.404183795000037], [85.39151582900011, 23.404540182000037], [85.39159703000009, 23.406257804000045], [85.39151910900003, 23.406662156000017], [85.39116419100006, 23.40779660400009], [85.39048428200003, 23.409454190000105], [85.39029103900002, 23.409864022000086], [85.38975258400006, 23.410536803000063], [85.38880983000004, 23.411684578000084], [85.38823647300002, 23.412172375000033], [85.38700544500004, 23.413240183500054], [85.38700112000004, 23.414463464500017], [85.387256108, 23.41468004000004], [85.38895024300011, 23.415359513000038], [85.38987056100014, 23.41566822300009], [85.39253397400006, 23.415928380000032], [85.392865972, 23.416014498000052], [85.39321949900003, 23.416238289000006], [85.39419917200009, 23.417114093000084], [85.3943785870001, 23.41735223600007], [85.39521758000006, 23.418735496000036], [85.39631106900003, 23.418863490000057], [85.3981238410001, 23.419282617000096], [85.39931149400005, 23.419679748000092], [85.39979215600002, 23.419818337000013], [85.4002059575001, 23.41908090550006], [85.40028919500011, 23.41862057200003], [85.40023233100003, 23.417976909000064], [85.40008761100005, 23.417043514000017], [85.39973495900009, 23.416078441000014], [85.39962244933344, 23.41518006400011], [85.40110720400011, 23.414740727500043], [85.40205923233344, 23.414234273000016], [85.40204324700002, 23.413697910000046], [85.40212206300006, 23.413294072500072], [85.4024487240001, 23.412652749000078]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;194&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02684-373561&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2684, &quot;tv_name&quot;: &quot;Buti&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 373561}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.2972074610001, 23.395634051000062, 85.30842608100009, 23.410600172000066], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.29726477400004, 23.395634051000062], [85.29752824900008, 23.396667346000115], [85.29806020000004, 23.398487885000023], [85.29853363600004, 23.39983242200007], [85.2994562980001, 23.402180296], [85.29988758350011, 23.403142252500146], [85.29996301600009, 23.40385428600004], [85.29988968200007, 23.405460739000088], [85.29931295200004, 23.406033359000038], [85.29886617300002, 23.406302400000094], [85.2972074610001, 23.406584333000062], [85.29797908400008, 23.407598044000093], [85.29816687100008, 23.407901288000048], [85.29875373400004, 23.409033984000065], [85.29902047700008, 23.40942098900011], [85.29942166100007, 23.409836781500008], [85.29966589200002, 23.410010113000094], [85.30058632400007, 23.41022020300007], [85.30154588600006, 23.410420741000053], [85.30244866900006, 23.410531344000056], [85.30323514700002, 23.41056000000003], [85.30611824, 23.410600172000066], [85.30625657300008, 23.40923932300009], [85.30637288700002, 23.408720543000065], [85.30671436100003, 23.408116453999998], [85.30687744100013, 23.407814139500076], [85.3069952380001, 23.40619597700008], [85.3071165410001, 23.40515302800003], [85.30719884100007, 23.404416143000077], [85.3074015185001, 23.403488601500115], [85.30801050999999, 23.402648542000108], [85.30839960800003, 23.4017850003334], [85.30842608100009, 23.401436333000085], [85.30834760800008, 23.40112604600008], [85.30779177650005, 23.399850523000058], [85.30650647800007, 23.398427909000034], [85.30618564000008, 23.39819659000008], [85.3035848440001, 23.396667798000063], [85.30304542300004, 23.39647263299999], [85.30247643600005, 23.396323268000067], [85.30124357100009, 23.396143492000114], [85.29726477400004, 23.395634051000062]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;197&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02684-373564&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2684, &quot;tv_name&quot;: &quot;Misirgonda Alias Pahargonda&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 373564}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.52659886671437, 23.53795941357148, 85.53443560480004, 23.544132385000076], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.53443560480004, 23.541081150000036], [85.53215027400006, 23.539372582000066], [85.53143768000007, 23.53870588500006], [85.53053354228577, 23.53795941357148], [85.52992606900006, 23.538073730000068], [85.52942829300008, 23.538227091000056], [85.52879287800009, 23.53838683700004], [85.52848554700006, 23.538517620000025], [85.52785743100004, 23.53892857000001], [85.52659886671437, 23.539754666571426], [85.52667890600006, 23.540058377000054], [85.52676234450007, 23.541377146500057], [85.52814807900002, 23.54252081400005], [85.52990248000006, 23.54388590000011], [85.53021904000005, 23.544132385000076], [85.53164912900016, 23.542598557000087], [85.53193033500007, 23.542426669000076], [85.53409465700008, 23.541303565000053], [85.53443560480004, 23.541081150000036]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;247&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02685-373613&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2685, &quot;tv_name&quot;: &quot;Jaruadihi&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 373613}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.54133764100003, 23.47466361300005, 85.54930425480005, 23.488588821000082], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.54930425480005, 23.487198324399998], [85.54853401200002, 23.48626870000005], [85.54797081400005, 23.48521980100004], [85.54754394999999, 23.48451923500007], [85.54732221200004, 23.483765250000037], [85.54723399700002, 23.483319278000085], [85.54691987700005, 23.482714207000114], [85.54668787100005, 23.48220353300005], [85.54654449900008, 23.48137494000006], [85.54640130000004, 23.480577948000107], [85.546184022, 23.480051539000073], [85.54609711200001, 23.479676471000037], [85.54593484600008, 23.477551415500024], [85.54541274100005, 23.477214778000103], [85.54492875700006, 23.47683258100001], [85.54473007800014, 23.476216389000097], [85.54508839700004, 23.475925489200115], [85.54375641400006, 23.47466361300005], [85.54314220300003, 23.47545007900002], [85.54293577200009, 23.4762400380001], [85.54236882900008, 23.47761274200008], [85.54164101700009, 23.479187351000068], [85.54142011900011, 23.47999987300006], [85.54133764100003, 23.48097264000008], [85.54136509000008, 23.481818252000057], [85.54152866900003, 23.482315217000053], [85.5416389310001, 23.482647294000085], [85.54236530900006, 23.483739445000033], [85.54262913800004, 23.484200699000006], [85.54271638100006, 23.48471338200006], [85.54309236100009, 23.48596991900011], [85.54357733900007, 23.4867029970001], [85.54391657700008, 23.48695964000011], [85.544653323, 23.48736492400007], [85.54505100300005, 23.487566165000022], [85.5454299545001, 23.487833147000085], [85.54584176200011, 23.488588821000082], [85.54646441600002, 23.488207630999984], [85.54703456800004, 23.487690529000098], [85.5474614870001, 23.487579353000047], [85.54831786400008, 23.4874823010001], [85.54894829700002, 23.4872879980001], [85.54930425480005, 23.487198324399998]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;280&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02685-373646&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2685, &quot;tv_name&quot;: &quot;Sidarmuta&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 373646}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.5481177724286, 23.457409193250093, 85.56314496425006, 23.470600838400017], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.56296235400004, 23.46995722700011], [85.56314496425006, 23.469461895750047], [85.56311216400007, 23.46894724600003], [85.56270890100011, 23.467669238400017], [85.56234945700004, 23.46680296900007], [85.5617157980001, 23.46529092600002], [85.56163077600007, 23.465010746000107], [85.56143863650003, 23.46410177800004], [85.56126422000004, 23.463746537000045], [85.56016229250004, 23.46152921600008], [85.55988876400009, 23.460284804000057], [85.55960623500008, 23.459055517000042], [85.55937291500014, 23.45859789650006], [85.55884695700004, 23.458109810000042], [85.558280741, 23.457624308000018], [85.55775528775003, 23.457409193250093], [85.55641857900012, 23.457602192000053], [85.556052145, 23.457728959333394], [85.5550618817501, 23.459197840250116], [85.55407060500009, 23.459084184000105], [85.552918859, 23.45884305600005], [85.5520482060001, 23.458582864000014], [85.55166263300008, 23.45845470600007], [85.5507988850001, 23.45796331600009], [85.55048770300003, 23.457953578000083], [85.54998203500011, 23.45804047999999], [85.54956965700009, 23.45817787400008], [85.5493196190001, 23.458333629], [85.5489986760001, 23.458555622500114], [85.54856124100006, 23.45991199200003], [85.5481177724286, 23.46099830585721], [85.54852624400009, 23.461789792000058], [85.548753795, 23.463315836000024], [85.54877065700006, 23.46391836800005], [85.54857824200005, 23.46462373500006], [85.54848200900007, 23.465250312000034], [85.54837966400004, 23.465671045000118], [85.54839821250012, 23.466370037500088], [85.54982661100004, 23.46782632900001], [85.55010466000004, 23.468109758000015], [85.5504336482501, 23.468482810500106], [85.55032999800004, 23.46916274200001], [85.55042667650008, 23.46961156100002], [85.55108144800003, 23.470201358000057], [85.55137122240006, 23.470461704200122], [85.55174188700006, 23.470427440000023], [85.55230716200002, 23.470308204000016], [85.55291805400003, 23.470301031000044], [85.5534321260001, 23.470412788000054], [85.55401069460007, 23.470600838400017], [85.55434522600008, 23.470551582999992], [85.55487781300009, 23.47028106500008], [85.55681039300003, 23.47028750400008], [85.55861703800012, 23.47016489100006], [85.56008908600006, 23.46995020500006], [85.56209026000008, 23.469875246000075], [85.56296235400004, 23.46995722700011]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;290&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02685-373656&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2685, &quot;tv_name&quot;: &quot;Hurputoli (Haruptoli)&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 373656}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.65158921600008, 23.380746534833367, 85.65673456900001, 23.384385240000054], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.65615082700002, 23.382369086000036], [85.65556798000003, 23.38168690900008], [85.65505949800006, 23.381142203000042], [85.65468985600002, 23.381031533000034], [85.6535431670001, 23.380929294000012], [85.65240269200005, 23.380790365000063], [85.65175353666675, 23.380746534833367], [85.65158921600008, 23.38105223400006], [85.65166942100005, 23.38320761600005], [85.6518590580001, 23.383975582000087], [85.65407626300006, 23.38341548600002], [85.65477833200006, 23.383244060000045], [85.65506782050008, 23.38322797850009], [85.65532910400007, 23.383340811000057], [85.65630570200005, 23.384082574000058], [85.65673456900001, 23.384385240000054], [85.65642486400004, 23.38323336300005], [85.65615082700002, 23.382369086000036]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;357&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02686-373723&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2686, &quot;tv_name&quot;: &quot;Ikarhatu&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 373723}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.80920679500011, 23.357154621000024, 85.8318292655001, 23.370689329000076], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.82325794700007, 23.36947910600004], [85.82333633500008, 23.36909625699998], [85.82338937800012, 23.36853806800007], [85.82334034700006, 23.368119605000086], [85.82323446, 23.367778485000017], [85.82310613000004, 23.367354723000037], [85.8229577440001, 23.366901158500035], [85.8228110925001, 23.366594972500014], [85.82256508300013, 23.366245463000006], [85.82236751700003, 23.365993835000058], [85.82215115100007, 23.365710527000065], [85.82194951825011, 23.365362977000085], [85.82191684900005, 23.364937911], [85.82201881800013, 23.36376485250006], [85.82277383300008, 23.363294474], [85.82335936400003, 23.363086668000037], [85.82364429900007, 23.362984887000053], [85.82398055500009, 23.362959540000087], [85.82576646300006, 23.362811819000033], [85.82697536799999, 23.362622195000064], [85.82759728600011, 23.362560098000078], [85.8291481310001, 23.362264397000036], [85.83025027400006, 23.361943189000016], [85.83084037850008, 23.361658356000074], [85.83136289500003, 23.361158350000004], [85.8318292655001, 23.360740360000058], [85.83151291100009, 23.360330841000103], [85.83084530400009, 23.359809860000073], [85.83056053200006, 23.35962811250009], [85.83010631500005, 23.359440387000056], [85.8289669220001, 23.359026999000047], [85.82799166100006, 23.358600334000094], [85.8272298860001, 23.35825497200003], [85.82671963100006, 23.358344470000063], [85.82463004600001, 23.358685725000054], [85.82363306100002, 23.35878071500001], [85.8232893900001, 23.358767070000034], [85.821953906, 23.3586218830001], [85.82003169400012, 23.35825499400003], [85.81931804600005, 23.35811866400003], [85.81849648800002, 23.357975231000065], [85.81812952600004, 23.35788022800002], [85.8174489478001, 23.357666375600054], [85.81601462300002, 23.357300388000095], [85.81526824900008, 23.357204014000025], [85.81479218600003, 23.357154621000024], [85.81445290900002, 23.35715518600007], [85.81376824600005, 23.357233824000005], [85.81317369300005, 23.357315992000053], [85.81287060500003, 23.35738657600003], [85.81241313000001, 23.357560909000068], [85.81143760100012, 23.357998054000088], [85.81070941700008, 23.35848492500005], [85.80961042900002, 23.359343656000046], [85.80954439050011, 23.36010876050007], [85.80926502150008, 23.360575874000126], [85.80920679500011, 23.361978592000103], [85.80921689100005, 23.362308481000042], [85.80930692100013, 23.362611033500073], [85.80947695100008, 23.362894376000025], [85.81007839200001, 23.363606288000042], [85.81033840650008, 23.363838749500047], [85.81076315600009, 23.364038916000045], [85.81115105300003, 23.36418544900009], [85.81164534300008, 23.364360142000066], [85.81203368800006, 23.36456340800001], [85.81319289750009, 23.36536095450005], [85.81333199725003, 23.36562001950008], [85.81322657400008, 23.366869063], [85.81315951700005, 23.36723453100004], [85.81296442400006, 23.367811645000096], [85.81219587600006, 23.369846545000065], [85.8119118210001, 23.370689329000076], [85.81724593100002, 23.370607021000104], [85.82120862800006, 23.370565048000095], [85.82152563250006, 23.370520629500106], [85.82232731699999, 23.370035737], [85.82325794700007, 23.36947910600004]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;495&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02688-373861&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2688, &quot;tv_name&quot;: &quot;Kachojara&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 373861}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.865925743, 23.340253373000046, 85.8772120930001, 23.352886864000048], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.86958257500015, 23.35136883350008], [85.87015860300012, 23.350975119250077], [85.87042274900003, 23.350564587000104], [85.87084519700008, 23.349852269000053], [85.87101349100016, 23.349591657500028], [85.87124087800005, 23.349401477000058], [85.87209930800002, 23.34896891800006], [85.87246103550007, 23.348872734500052], [85.87301139100008, 23.348822265000006], [85.87455711000007, 23.348718383000026], [85.87486543350013, 23.348594655000092], [85.875398853, 23.348178310000016], [85.87594360900005, 23.34774602100004], [85.87630547900005, 23.347516948000095], [85.87671879300005, 23.347304876000045], [85.8772120930001, 23.34708381700005], [85.87651686200005, 23.346610139000028], [85.87601575900005, 23.346268727000023], [85.8754945620001, 23.34585935399999], [85.8747333810001, 23.345250507500108], [85.87398819550012, 23.344549610500113], [85.87342174400001, 23.344047206000067], [85.87300009450007, 23.34381654250012], [85.87251522700006, 23.343528578000054], [85.87229357500001, 23.343291945000104], [85.87215056900008, 23.34302020400009], [85.87208880575011, 23.342550137750077], [85.87216551500008, 23.341987486999997], [85.87231735800003, 23.34154533350005], [85.87259378000007, 23.340909867833386], [85.87259146250011, 23.34049157750008], [85.8724258546667, 23.340253373000046], [85.87187532100005, 23.34028922900008], [85.87123258000007, 23.340409747000017], [85.87064737700005, 23.340519476000036], [85.86944694500004, 23.34082259300007], [85.86900105500008, 23.34088447600004], [85.86859655300002, 23.34090142100011], [85.86829034500008, 23.34085450550009], [85.86802187400006, 23.34072979500012], [85.86778197335148, 23.340564743815577], [85.86763694900002, 23.343081440000006], [85.86753893000008, 23.344499613000018], [85.86741851100007, 23.345481442000075], [85.86725855200007, 23.346348156000023], [85.86706195300005, 23.34706314600003], [85.86686229000004, 23.34762209600009], [85.86655916200002, 23.348155194000057], [85.86606269400006, 23.348813990000032], [85.86594959400011, 23.34907864600008], [85.865925743, 23.34937480000007], [85.86607637600012, 23.351047151000063], [85.86623997700008, 23.35190621700009], [85.86633399800003, 23.352239062000006], [85.86654669450009, 23.352564115500034], [85.86696417700004, 23.352886864000048], [85.86769749400005, 23.35243841400007], [85.86910716500007, 23.351626358000026], [85.86958257500015, 23.35136883350008]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;511&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02688-373877&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2688, &quot;tv_name&quot;: &quot;Rangamati&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 373877}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.84462840500001, 23.316905093000077, 85.85833683300004, 23.3263801805001], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.85297898600002, 23.326338656000082], [85.8538880630001, 23.326326134000023], [85.85507108700001, 23.32632859900005], [85.85540223500006, 23.32625231266671], [85.855767578, 23.32590901200006], [85.85619080600011, 23.325308971000062], [85.8569088400001, 23.324369373000035], [85.85761529600002, 23.32354649050003], [85.85789109600002, 23.322965326000038], [85.85833683300004, 23.32207200600002], [85.85812326400004, 23.32185637900003], [85.85780762100003, 23.321626806000044], [85.85744292825008, 23.321446021000096], [85.85707285733334, 23.321223675000084], [85.85700610700007, 23.320741635500035], [85.85711947200012, 23.320166490000098], [85.85728073175012, 23.319343526250062], [85.8571304225, 23.31885063550007], [85.85664229600008, 23.318430348000067], [85.85622681100007, 23.318119379000095], [85.85514423350007, 23.317336734000087], [85.85455405400005, 23.316905093000077], [85.8515867020001, 23.317255179], [85.85093721450006, 23.31734895200012], [85.8498883780001, 23.31795283400008], [85.84941845800003, 23.318222472000098], [85.84854790900003, 23.31868299900001], [85.84783824200011, 23.318943485750136], [85.84620115450006, 23.31855748350006], [85.84584729500004, 23.31855647700011], [85.84510460700005, 23.318751902000027], [85.84462840500001, 23.318886433000056], [85.8449440310001, 23.32015947800005], [85.84511869700007, 23.32122754200005], [85.84517787500008, 23.321853093000115], [85.84524063200003, 23.32263891200003], [85.84523666700007, 23.323334587000037], [85.84512570300002, 23.324210730000054], [85.84495919300011, 23.325333184000023], [85.84496209780006, 23.325761506000024], [85.84632295300005, 23.326336333000068], [85.84663556000008, 23.3263801805001], [85.8478928230001, 23.326373498000034], [85.84902061800007, 23.326368863000006], [85.85101931700005, 23.32635890000007], [85.85297898600002, 23.326338656000082]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;519&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02688-373885&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2688, &quot;tv_name&quot;: &quot;Dumarbera&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 373885}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.31416227900006, 23.276385378000068, 85.33068281300007, 23.288463620200066], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.33068281300007, 23.283274203000108], [85.32966949600008, 23.283184915000092], [85.32899838700006, 23.28295372600006], [85.32839604800004, 23.28254439900003], [85.32786953350005, 23.28197106100011], [85.32771969050009, 23.281363211000087], [85.32756479100001, 23.281108654000075], [85.32735477200005, 23.28089605800007], [85.32708837600003, 23.280683575000076], [85.32687734600006, 23.280484524000038], [85.32666703500007, 23.280203616000065], [85.32641101700003, 23.279417396000028], [85.32615826500007, 23.279097328000102], [85.32573841300007, 23.27883345500003], [85.32493953100004, 23.278429043000052], [85.32418361800003, 23.27798258200006], [85.32379022700007, 23.277610006000053], [85.32302793600006, 23.276660483000093], [85.32177967000008, 23.27667260400003], [85.32087630800002, 23.27667927700005], [85.32046130800006, 23.27667620600006], [85.32002171700002, 23.276602400000083], [85.31949614400006, 23.276470833000033], [85.31917935500005, 23.27639459800003], [85.31886100900009, 23.276389733000073], [85.31847081600004, 23.276385378000068], [85.317494584, 23.27652458800004], [85.31634754400004, 23.276786229000017], [85.31583418500011, 23.276940039000042], [85.315187753, 23.277368405000047], [85.3149162660001, 23.277624740000117], [85.31416227900006, 23.278505505000073], [85.31468154500004, 23.279370964000094], [85.31646939100008, 23.281704530000084], [85.3177898470001, 23.28263647100001], [85.31873813100005, 23.28315997400007], [85.31916372100008, 23.28346617250005], [85.32136459000002, 23.284224307000027], [85.32253232200009, 23.284855972000067], [85.32275646900007, 23.285086274000026], [85.32317146800004, 23.285421549000056], [85.3233354350001, 23.285719907500095], [85.32335931650009, 23.286192002500076], [85.32330113850007, 23.286608923500022], [85.32329476400011, 23.28696655800008], [85.32329599400003, 23.287361113000042], [85.32329712100004, 23.287726826000014], [85.32333541450006, 23.288090838000024], [85.32356675780004, 23.288463620200066], [85.32455514200004, 23.287499174000097], [85.32532747750012, 23.287347456000067], [85.32764010800011, 23.286466946000044], [85.32812576300006, 23.286391595000026], [85.32887808900006, 23.286326690000074], [85.32951755400012, 23.286397404250085], [85.33002409833344, 23.286420799000044], [85.33030792000005, 23.28630052200004], [85.33044311700004, 23.286043047000053], [85.3305737160001, 23.284612753000015], [85.33050960600004, 23.283896099500094], [85.33068281300007, 23.283274203000108]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;602&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02690-373968&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2690, &quot;tv_name&quot;: &quot;Hatia&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 373968}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.34641770800009, 23.21559887900006, 85.35354693200013, 23.227027978000063], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.35253674500012, 23.2183844425001], [85.35286836600002, 23.217289299000036], [85.35354693200013, 23.21559887900006], [85.35266320300002, 23.215644435500053], [85.35222399800008, 23.215757466000035], [85.35137952800002, 23.216061157000027], [85.35097464700002, 23.216211069000085], [85.3505929185001, 23.21641327800003], [85.35025695700006, 23.216480385000075], [85.34883420000001, 23.216831755000015], [85.34757050800009, 23.217074172], [85.3471215505001, 23.217191722000134], [85.34686465700013, 23.217313621000095], [85.3465917740001, 23.217624401500096], [85.34641770800009, 23.218267575000066], [85.34694892400007, 23.218744980000043], [85.34720423000006, 23.219196785000033], [85.34730682950007, 23.21948515100009], [85.34732199400013, 23.21978704200002], [85.34731581800006, 23.220679917000002], [85.34721668000003, 23.22262405400001], [85.3474020410001, 23.22388009900002], [85.34766289300009, 23.22518496500007], [85.348278113, 23.22652495200001], [85.34846114300005, 23.227027978000063], [85.34893355900006, 23.22701817000005], [85.34926660350004, 23.226929523000088], [85.34947970600004, 23.226726272000057], [85.34980966700007, 23.226345062000046], [85.35051451200005, 23.22570741750009], [85.3508247985001, 23.225672921500077], [85.35112228900006, 23.225516689000074], [85.3519311680001, 23.22379660000004], [85.35204947300012, 23.223306257000058], [85.35205346400001, 23.222921012000068], [85.35213044100007, 23.220019298000054], [85.35212570400004, 23.21941240400002], [85.35225324500003, 23.218888765000038], [85.35253674500012, 23.2183844425001]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;639&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02690-374005&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2690, &quot;tv_name&quot;: &quot;Guniyadih&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 374005}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.33725565950013, 23.280753064500097, 85.34965520725011, 23.303365068800012], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.34684198100011, 23.280969350000095], [85.34653023500007, 23.28077845200002], [85.34612645100006, 23.280753064500097], [85.34577675100009, 23.280833865999988], [85.345482786, 23.28097405400004], [85.34497236700005, 23.28129940800005], [85.34430762300012, 23.281792808000016], [85.34356525600006, 23.28254642400002], [85.3429623110001, 23.28317462000004], [85.3426543510001, 23.283632726000107], [85.34249473050009, 23.284141083500113], [85.34253535900002, 23.287558569000037], [85.34259439533336, 23.28881186566675], [85.34193734700003, 23.289136189000008], [85.34164259900005, 23.289174293000052], [85.34134728700005, 23.289178829000072], [85.34100946800008, 23.289175382000053], [85.34047118600003, 23.28914629250009], [85.33895100900011, 23.288787324000097], [85.33842008700006, 23.28865243600004], [85.33809475600005, 23.288632659500045], [85.33779572900009, 23.288662529000074], [85.33749821720001, 23.28884096340004], [85.33725565950013, 23.29083732850003], [85.33736196300008, 23.29113859900013], [85.33759543000004, 23.291339855000086], [85.33858362900008, 23.291624344000088], [85.33922044500011, 23.292101856000073], [85.33949988700004, 23.292426400000103], [85.33991171900001, 23.293203857000037], [85.34015915500007, 23.293873252], [85.34031059100008, 23.294546023000066], [85.34033906800003, 23.295060037000052], [85.34023724600004, 23.29544917500001], [85.33989807800013, 23.29644727100008], [85.33923798000005, 23.296961031000045], [85.34004111000007, 23.297399359000075], [85.34091159000006, 23.29786563100002], [85.34128159700003, 23.29818782900003], [85.34154627000008, 23.299241294000065], [85.34164387200008, 23.299851234000084], [85.34249711900006, 23.301098604000085], [85.34326966500004, 23.302060736000023], [85.34410861600009, 23.302763064000047], [85.3445680532501, 23.302772307250034], [85.34558509300004, 23.302670392000017], [85.34690347000004, 23.302576125000034], [85.34713922100002, 23.30238736833338], [85.3473656740001, 23.301995073666674], [85.34767549650012, 23.301654251250074], [85.34866110340003, 23.303365068800012], [85.34913607900003, 23.302057914000056], [85.34922451100007, 23.30150988400001], [85.3492292500001, 23.301206095000097], [85.34921930300008, 23.300736211000032], [85.34896798400007, 23.299508009000107], [85.34893846200002, 23.29920365600006], [85.34893413550007, 23.29813408950007], [85.34911926400007, 23.29749271200007], [85.34927977500011, 23.29706842200005], [85.34965520725011, 23.29630013475005], [85.34957684600009, 23.29559886200005], [85.34955418200008, 23.29489312500001], [85.34962137300005, 23.2942920300001], [85.3495663510001, 23.2923523100001], [85.34941936400011, 23.289679110000076], [85.34905782200005, 23.287134196000096], [85.34877387600008, 23.285395568000038], [85.34829076800008, 23.283732749000034], [85.34766672900005, 23.282432917000044], [85.34684198100011, 23.280969350000095]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;646&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02690-374012&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2690, &quot;tv_name&quot;: &quot;Hethu&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 374012}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.36575850600012, 23.236368115000065, 85.37202739860008, 23.24562945497765], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.37202739860008, 23.240804681000064], [85.37188877200002, 23.23963588599999], [85.37171839000004, 23.238792440000015], [85.37133052800002, 23.237916618000025], [85.37092300100005, 23.23701212100006], [85.37067100400007, 23.236652365000047], [85.37037606100009, 23.236368115000065], [85.36940538400006, 23.23786956600004], [85.36889133100007, 23.2384038220001], [85.36848136900008, 23.238746090000088], [85.36783970100002, 23.23914457600002], [85.36711206150008, 23.239554979000047], [85.3669063033334, 23.239844389333435], [85.36699075900003, 23.240221999500125], [85.36726872000011, 23.240505315000075], [85.36786962850015, 23.24077334650009], [85.36791678725011, 23.241273068750072], [85.36735262400003, 23.241572211000065], [85.36687628000006, 23.24171656600002], [85.36648444200006, 23.241890527000024], [85.36595586150007, 23.242259028000138], [85.36580156350009, 23.242518293000096], [85.36575850600012, 23.24286779125008], [85.36596959900001, 23.243396855000046], [85.36655059400003, 23.244330388000016], [85.36696392700011, 23.244705619000065], [85.3680098601221, 23.24562945497765], [85.36820227599999, 23.24531159000002], [85.37050614200011, 23.242242758000035], [85.37132816399999, 23.24111365100005], [85.37157051800007, 23.240906087000017], [85.37202739860008, 23.240804681000064]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;655&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02690-374021&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2690, &quot;tv_name&quot;: &quot;Gobindpur&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 374021}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.39273466400006, 23.24528070100004, 85.40443170400012, 23.257203148000098], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.40420428900005, 23.249775436000057], [85.40268878300006, 23.24743401500004], [85.40244907300003, 23.247148093000096], [85.40185111700004, 23.246746340000072], [85.4010017740001, 23.246481689000024], [85.39961285600002, 23.2456399540001], [85.39889668800001, 23.24539613200009], [85.39649673600012, 23.24528070100004], [85.39599450500005, 23.245285873000025], [85.39568073400001, 23.245289128000028], [85.39535554400011, 23.24536172850013], [85.39494339200003, 23.24558483600007], [85.39468141450007, 23.24584485200009], [85.39429556400007, 23.246243938500115], [85.39398437000008, 23.24653509100009], [85.39357154600012, 23.2466645575001], [85.39311969000012, 23.246606544000052], [85.39276537100002, 23.246322444999993], [85.39273466400006, 23.246656906000045], [85.39274944000005, 23.24701611300001], [85.39296507300011, 23.247918547000072], [85.39313924900006, 23.249348185000056], [85.3932149490001, 23.250042615000048], [85.39335847600007, 23.250574786000115], [85.39390047400006, 23.251545309000036], [85.39429991600008, 23.25215601700003], [85.3944011830001, 23.25254360749999], [85.39433306400007, 23.253418458000034], [85.394335556, 23.253731924000036], [85.39451335200006, 23.254262825000087], [85.39490331700009, 23.254935231000033], [85.39538298400007, 23.255522210000017], [85.39571081400001, 23.25584380000007], [85.39647982400005, 23.256428821000068], [85.39686536200007, 23.256704435000046], [85.39779180900007, 23.25715249700008], [85.39868284800008, 23.256826770000043], [85.39921686300004, 23.25678500700005], [85.3996241045001, 23.256801246000062], [85.40014636400011, 23.256930768000032], [85.4009491520001, 23.25715496700002], [85.40142676600009, 23.257203148000098], [85.40194789850008, 23.257142317500083], [85.4025304340001, 23.256667721000028], [85.4029060970001, 23.255974972999997], [85.40331658800005, 23.254823289000097], [85.40344504450006, 23.253732170000056], [85.40366903800007, 23.25337600050003], [85.4039382455001, 23.252937925000065], [85.40420428000004, 23.252036221000065], [85.40436321400004, 23.25130786300002], [85.40443170400012, 23.250625969000087], [85.40420428900005, 23.249775436000057]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;672&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02690-374038&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2690, &quot;tv_name&quot;: &quot;Jordag&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 374038}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.3952803760001, 23.26891153500003, 85.40015923300011, 23.278290334600086], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.39888262699999, 23.26891153500003], [85.3981582480001, 23.269307393000098], [85.3976912950001, 23.269731262999997], [85.39735937000009, 23.27010787300007], [85.39712793400012, 23.270530637000057], [85.39684016200002, 23.271306366000026], [85.39660027000006, 23.27198650700008], [85.39639669400002, 23.272261231000073], [85.39588251350006, 23.272920792500067], [85.39573959350004, 23.273396260000105], [85.39574486500008, 23.27397013600007], [85.39566228600009, 23.27460739900004], [85.39542909600007, 23.275842995000044], [85.3952803760001, 23.277038212999994], [85.39555001400004, 23.27740614000006], [85.39565885400008, 23.27770573500007], [85.39583377400007, 23.278005061000055], [85.39621726850014, 23.278290334600086], [85.39662996000013, 23.278054874500047], [85.39692361000009, 23.277780931000073], [85.39793844400003, 23.27683374800004], [85.39891599500001, 23.27604654600001], [85.39984016450013, 23.275427792500107], [85.40003275400004, 23.275169377000086], [85.40014903100007, 23.27416907900001], [85.40015923300011, 23.2717580940001], [85.40005497200009, 23.270023284750067], [85.39979828100009, 23.26972689500008], [85.39888262699999, 23.26891153500003]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;681&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02690-374047&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2690, &quot;tv_name&quot;: &quot;Arma&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 374047}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.2444536850001, 23.282789181000034, 85.26582984000004, 23.302405091000004], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.26241862600006, 23.292039959000103], [85.26292874800008, 23.291724025000075], [85.26340434600009, 23.29131842800009], [85.26360648100011, 23.290976525000044], [85.2637630290001, 23.290547061000037], [85.2638338420001, 23.290249330000094], [85.26385338600002, 23.28988251300003], [85.263838128, 23.289562584000027], [85.26383782800006, 23.28927695300011], [85.26392349600005, 23.28894037800006], [85.26422831900005, 23.287822018000075], [85.26450277000002, 23.28693316300007], [85.26582984000004, 23.28386795500003], [85.26465662000005, 23.28338658500008], [85.26407262900004, 23.283304078000093], [85.26377254900008, 23.28331060800001], [85.26184530900005, 23.283497787000044], [85.2611878440001, 23.283562119000038], [85.26086287050008, 23.28349128700004], [85.26027358650008, 23.282853248000094], [85.25984665800006, 23.282789181000034], [85.25727805200007, 23.282970127000084], [85.25617978700006, 23.2830314870001], [85.25508454500005, 23.282916819000025], [85.2543445650001, 23.282989596000007], [85.25278210600005, 23.283330371000048], [85.25230570500007, 23.283365735000043], [85.25190738400006, 23.28337438900005], [85.25095431300008, 23.283238359000073], [85.25055020600003, 23.283160592000034], [85.24957682600005, 23.284171204000074], [85.24921301500007, 23.284519482000025], [85.24873339500007, 23.284795796500106], [85.24835254533343, 23.28506067700004], [85.24819339200008, 23.28540418700004], [85.24770334800012, 23.28635432200005], [85.24761200200003, 23.286662057000022], [85.24758817050005, 23.287033709500086], [85.24765748000006, 23.287312993000057], [85.24786800625017, 23.287924805000053], [85.2476784980001, 23.28832750050011], [85.24727483150005, 23.288583838000047], [85.24677815400013, 23.288896535000074], [85.2463388790001, 23.289432989500096], [85.24646586175011, 23.290371126000135], [85.24688367600007, 23.290904844000064], [85.24755906750012, 23.291755443250032], [85.24699192800006, 23.292278487000033], [85.24659438900011, 23.29247363850006], [85.24620122700009, 23.292469611000037], [85.2456935740001, 23.292360468000066], [85.24537140700008, 23.292242384000065], [85.24497569925005, 23.292148562500163], [85.2444536850001, 23.29251603500006], [85.24485815800007, 23.29296594300007], [85.24547370000002, 23.293779599000068], [85.24597104600005, 23.294543891000068], [85.24612493900008, 23.294856863000028], [85.24626826500013, 23.295249940000044], [85.2466317560001, 23.296346922000037], [85.24692237400006, 23.297385766000016], [85.24739359800007, 23.299317936000012], [85.2477074280001, 23.300657254000072], [85.24817050350009, 23.301884797500076], [85.24855224120012, 23.302269740400032], [85.24919186900009, 23.302405091000004], [85.25006865400007, 23.30229102500007], [85.25144219400003, 23.302000886000044], [85.25280532800011, 23.301555189000027], [85.25367257, 23.30114517100003], [85.25413390300008, 23.300873617000033], [85.25436846900003, 23.3006704680001], [85.25476237150016, 23.300268255500082], [85.25509698800009, 23.299799033000028], [85.25568454400003, 23.29877570600009], [85.256229445, 23.29769555600007], [85.25693890900008, 23.296318613000004], [85.25727478900008, 23.295620886000094], [85.25745217900001, 23.295341054000033], [85.25791949300003, 23.29470781000004], [85.25877416900005, 23.29383734500005], [85.2592395070001, 23.2933570550001], [85.25989836000008, 23.29271963500003], [85.260136919, 23.29248744900009], [85.26055705100005, 23.29225721500005], [85.26088386300009, 23.292153030000044], [85.26133363200009, 23.292035535000004], [85.26176501700007, 23.291997373000044], [85.26241862600006, 23.292039959000103]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;776&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02692-374142&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2692, &quot;tv_name&quot;: &quot;Nachiatu&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 374142}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [84.98376454300008, 23.136716548750062, 84.99603274780003, 23.14937236000007], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[84.99499008800002, 23.146897647000035], [84.99507673275015, 23.145820800750087], [84.99555629400007, 23.14538353700008], [84.99603274780003, 23.1448634776], [84.99556504100008, 23.144182692000108], [84.99505272400006, 23.143294047000033], [84.99433909300008, 23.14184722000005], [84.9939840390001, 23.141337401000055], [84.99300748600002, 23.140608175000025], [84.99144073100007, 23.13962147400001], [84.98969801700012, 23.138613601000017], [84.98879477900007, 23.138290079000072], [84.98844241800006, 23.13813739900005], [84.98759955000013, 23.136849810000097], [84.98713142800003, 23.136716548750062], [84.98670056200008, 23.13715340900008], [84.98524512860007, 23.137948787200095], [84.98537064900009, 23.13868089400008], [84.98536514400008, 23.140725795999987], [84.98511360200006, 23.142720157000067], [84.98496635100012, 23.14459443400007], [84.984724408, 23.14536853800008], [84.98449747500005, 23.145694586000058], [84.98439399650006, 23.146090229500064], [84.98445270300009, 23.14711146500005], [84.9844019375001, 23.14747087150008], [84.98386068600004, 23.148301809000042], [84.98376454300008, 23.148767927500106], [84.98401823095242, 23.14924557951014], [84.98638584200008, 23.14937236000007], [84.98701996500004, 23.149238976000046], [84.99003305900004, 23.148250461000057], [84.99057734600004, 23.1480459370001], [84.991175905, 23.147859492000038], [84.99251410000007, 23.147734925000034], [84.99325443550008, 23.14774341850007], [84.99465890700007, 23.147080015000025], [84.99499008800002, 23.146897647000035]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;1080&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02697-374446&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2697, &quot;tv_name&quot;: &quot;Ranidera&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 374446}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.58978568300012, 23.131880887000015, 85.60254006250004, 23.140038223333452], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.60220303300008, 23.138040576000034], [85.60219157000012, 23.136781515500086], [85.60254006250004, 23.135438576500082], [85.60248605900006, 23.135076705000117], [85.60228493100001, 23.134139584000028], [85.60231095000005, 23.13322168900002], [85.60162001300007, 23.133037845999997], [85.59977609800009, 23.132426407000032], [85.59818731500009, 23.13195441300006], [85.5975895950001, 23.131880887000015], [85.59735044300004, 23.132086429000058], [85.59662773800001, 23.133029452000102], [85.59628840400009, 23.133314924000008], [85.5950385335001, 23.13404824000002], [85.59346148800013, 23.134263708000084], [85.59304341000009, 23.134200167000067], [85.5918765630001, 23.13381962900011], [85.59170048600001, 23.134434437], [85.59112670700004, 23.13533068000008], [85.59089666350006, 23.13564034050008], [85.5905467470001, 23.135776439000015], [85.59024028500006, 23.135774649000044], [85.58978568300012, 23.1369025738334], [85.59112864300003, 23.13771598300009], [85.59238070000006, 23.137996086000083], [85.59283125100012, 23.13814501500006], [85.59388554800015, 23.139366739250093], [85.59530627183338, 23.140038223333452], [85.59657665900005, 23.139141466000083], [85.5969765310001, 23.138929403000034], [85.59833404400005, 23.138521983000075], [85.5992974640001, 23.138404724000054], [85.60065224500002, 23.138224271000045], [85.60220303300008, 23.138040576000034]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;1136&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02698-374502&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2698, &quot;tv_name&quot;: &quot;Tisro&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 374502}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.75193350700007, 22.981841440800014, 85.76727733800007, 22.996108112000073], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.76727733800007, 22.98921359633343], [85.76690687900006, 22.98815140800005], [85.76628160700002, 22.987716976000076], [85.76450632600003, 22.986494646000015], [85.76352396900002, 22.98549795300006], [85.76321446800009, 22.985164749000024], [85.7627552860001, 22.984342148000085], [85.76207421400007, 22.983261566000028], [85.76129062900007, 22.981994427000078], [85.76095311140007, 22.981841440800014], [85.7595723260001, 22.98366818800001], [85.75932454400002, 22.983924309000063], [85.75868393000003, 22.98416878800004], [85.7579998660001, 22.984252550500088], [85.75751423350012, 22.984208614500094], [85.75665764500008, 22.98408278100003], [85.75614139300008, 22.984051919000052], [85.75690423500005, 22.98465411100005], [85.75736891200012, 22.985115515000125], [85.75752409700009, 22.98622994500005], [85.75746924000008, 22.98653072300007], [85.75670079600009, 22.98798918100004], [85.75621119466675, 22.989184724833393], [85.75444484000003, 22.98982764200001], [85.75279767299999, 22.99057824200003], [85.75193350700007, 22.99105544600008], [85.75234040900008, 22.99181819700004], [85.75259542100005, 22.992243432000052], [85.75263176200012, 22.99255006700002], [85.75261221400011, 22.993808380200058], [85.75386765900006, 22.994477069000034], [85.75484741700004, 22.995275708000033], [85.75542367700011, 22.995849089000046], [85.75572130700004, 22.996108112000073], [85.7565518220001, 22.995443174000034], [85.75698053500005, 22.995102532000004], [85.75737637800003, 22.994650856000113], [85.75776499200003, 22.994129371000042], [85.75828242900003, 22.993100805000026], [85.75898646650006, 22.991520231000077], [85.75905564300007, 22.99112689400004], [85.75886892125008, 22.98978252950005], [85.75957080400005, 22.988419657500085], [85.759852434, 22.98833187800006], [85.760708084, 22.988279722000097], [85.76254089200005, 22.988228790666724], [85.76301529000004, 22.98862334900002], [85.76361887500006, 22.98926231400004], [85.76400479200002, 22.989403238000065], [85.76481549900004, 22.98953022000009], [85.7660882040001, 22.989499814000062], [85.76727733800007, 22.98921359633343]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;1238&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 364, &quot;district_name&quot;: &quot;ranchi&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-364-02699-374603&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2699, &quot;tv_name&quot;: &quot;Kadamdih&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 374603}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.12446830100006, 22.79738033400002, 85.15060643900003, 22.810195638000057], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.15060643900003, 22.805993426999997], [85.14910477100004, 22.80241674100009], [85.14691989100001, 22.802157372000035], [85.14606261500005, 22.802025459000046], [85.14536666500001, 22.80182466600004], [85.1438744800001, 22.801468552000063], [85.14221335800005, 22.801089142000023], [85.13952879000007, 22.80019612000003], [85.13839278200004, 22.79978951900005], [85.13705071400003, 22.799255379000087], [85.1359739760001, 22.798687845000114], [85.13521204000003, 22.798312404000047], [85.13424485700003, 22.797936126], [85.13326484300003, 22.797634774000112], [85.1322790070001, 22.797422971000085], [85.13199568, 22.797406010000014], [85.13157688100009, 22.79738033400002], [85.13047918000004, 22.797461163000047], [85.12801682400004, 22.798024427000023], [85.12730928300003, 22.798348551000018], [85.12691835800011, 22.798529103000035], [85.12665399100008, 22.79879051100008], [85.12639208000009, 22.799169432999992], [85.12611872400002, 22.799782551000096], [85.12594122100008, 22.800069647000047], [85.125696545, 22.800272652000057], [85.12539566100006, 22.800410498000094], [85.12510763400003, 22.80047546800005], [85.12446830100006, 22.800572513000088], [85.12501065492353, 22.801400396799824], [85.12547641600008, 22.801510941000018], [85.12607687600011, 22.80162524500006], [85.12639796000003, 22.801715087000034], [85.1267122670001, 22.80189441500004], [85.12738034800005, 22.80243782700006], [85.12780556300007, 22.80295155600002], [85.12896943099999, 22.80461171700007], [85.12946975100009, 22.80526410200005], [85.12987287400004, 22.805685891000042], [85.13154927700002, 22.80726193800006], [85.13248290400009, 22.80804354600007], [85.13290749700005, 22.808305040000036], [85.13368344100003, 22.80868843300005], [85.13564607800006, 22.80963496100003], [85.13632744900004, 22.809927527000077], [85.13676143300006, 22.810062598000062], [85.13755123300004, 22.810189773000015], [85.13793710700008, 22.810195638000057], [85.13826793100002, 22.810142560000017], [85.13893412000007, 22.809982454000092], [85.13919841500008, 22.80987291300009], [85.13968303000004, 22.809535084000025], [85.14063525000005, 22.808878864000086], [85.14101446200007, 22.808721768000055], [85.14171520900005, 22.80852944700004], [85.14293023400003, 22.808257218000076], [85.14513484000011, 22.807801719000047], [85.1462834480001, 22.807487160000044], [85.14803059700007, 22.80685511300006], [85.15060643900003, 22.805993426999997]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;1617&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 365, &quot;district_name&quot;: &quot;khunti&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-365-02702-374982&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2702, &quot;tv_name&quot;: &quot;Birhorchuan&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 374982}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.35547660080007, 23.14844669216675, 85.36219715780005, 23.15248306840009], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.36165490620007, 23.150615107600007], [85.36165070200006, 23.1499878270001], [85.36189280200001, 23.149308106999992], [85.36219715780005, 23.148930894200063], [85.36185511400001, 23.148694494000026], [85.36159106950005, 23.148518588000037], [85.36090931683347, 23.14844669216675], [85.36042123600004, 23.148808171000017], [85.36007685500005, 23.149024824], [85.35951324900005, 23.14929067400006], [85.35890719100009, 23.149548250000056], [85.35819958400006, 23.149791906000054], [85.35724382000004, 23.150060393000064], [85.35619625600005, 23.15034132900004], [85.35547660080007, 23.150478340800078], [85.35558126900006, 23.150915217000005], [85.35600187200004, 23.151772454000056], [85.35616196000005, 23.152104126000054], [85.35629862560003, 23.15248306840009], [85.35648778700009, 23.152085262000092], [85.35668826700007, 23.15157165700005], [85.35684971050006, 23.151236641500088], [85.35746192900005, 23.15109406400006], [85.35837193400002, 23.150996822000007], [85.3593839380001, 23.150894764], [85.36040650400005, 23.150825758000064], [85.3611691770001, 23.150765768000074], [85.36165490620007, 23.150615107600007]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;1908&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 365, &quot;district_name&quot;: &quot;khunti&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-365-02704-375273&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2704, &quot;tv_name&quot;: &quot;Chhota Banri&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 375273}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [85.57004344840003, 22.985635219200066, 85.58483503350007, 22.998818353000054], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[85.58412205560006, 22.985635219200066], [85.58355379750012, 22.985901161000026], [85.58300448900009, 22.986059669000085], [85.58257090700002, 22.986105847999998], [85.58208952700006, 22.98612819199998], [85.58152744200001, 22.98614763700008], [85.58055928100002, 22.98611374200003], [85.57914514800007, 22.986030071000034], [85.57825153200012, 22.986000378000064], [85.5770253330001, 22.985975226000026], [85.57605191700011, 22.98596574900004], [85.57525222500003, 22.98597264900008], [85.57442485700004, 22.98589081800006], [85.57288133200005, 22.985789184000087], [85.57246701200006, 22.98578350300006], [85.57209027300014, 22.985766974000114], [85.5717031470001, 22.985676598000097], [85.57117618050006, 22.985705897000035], [85.57081264800004, 22.985844532000048], [85.57037117100002, 22.986126248000087], [85.57004344840003, 22.986337222399992], [85.57045354200008, 22.98656890800003], [85.57090726400006, 22.986921198000047], [85.57173637400008, 22.987605338999987], [85.57231586500008, 22.987958942000024], [85.57304026700005, 22.98834225700003], [85.57350919400004, 22.988563985999992], [85.57392106000002, 22.988687843000026], [85.57433383600006, 22.988798063000086], [85.57466356100008, 22.988922683000094], [85.57505079500005, 22.989026953000106], [85.57539367200009, 22.98915830100005], [85.57591445800001, 22.989525274000037], [85.57642647600001, 22.989889134000034], [85.57678823500001, 22.99015647500004], [85.57722586600008, 22.99047912700009], [85.57758072100003, 22.99070456200009], [85.57828654200009, 22.99108504100004], [85.5788362910001, 22.991398415000035], [85.57913970300007, 22.99155879700007], [85.57940218300011, 22.99167008300004], [85.57980504300014, 22.991834872000087], [85.5800654160001, 22.99198806200006], [85.58023764600014, 22.99228666950012], [85.58045020800004, 22.992682022000015], [85.58062537000006, 22.993060616000076], [85.5806912180001, 22.993546405000036], [85.58086201500002, 22.99423694800006], [85.58110891700002, 22.995261803000037], [85.58122867000007, 22.996302680000007], [85.58119655400007, 22.997132068000063], [85.58116189900005, 22.997885063000105], [85.58115561100006, 22.998459250000096], [85.58161149200004, 22.998539963000024], [85.58204233500007, 22.998652915000044], [85.58234024400011, 22.998725673000024], [85.58290756400007, 22.99876215400011], [85.58386955499999, 22.998797968000055], [85.58439909200011, 22.998818353000054], [85.58407446950002, 22.998082921000044], [85.58396535900002, 22.997562074000026], [85.58400210800004, 22.996838649000033], [85.58396988800004, 22.995231627000027], [85.58404634900006, 22.994449544000066], [85.5840714180001, 22.99407685500005], [85.58412274100007, 22.993694330000018], [85.58427214500009, 22.992678133000048], [85.58426818800002, 22.991813989000047], [85.58428484700008, 22.991132551000064], [85.5842904850001, 22.990429449000064], [85.58441308000002, 22.98995849300009], [85.58457655000005, 22.989443048000073], [85.58478201700008, 22.98890496900005], [85.58483233800007, 22.988599673000103], [85.58483503350007, 22.988215672000038], [85.58476499300009, 22.987909737000038], [85.58460763800005, 22.987581843999997], [85.58447028700004, 22.98723001700002], [85.58439908500002, 22.986867005000093], [85.58430239500011, 22.9862358980001], [85.58412205560006, 22.985635219200066]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;2033&quot;, &quot;properties&quot;: {&quot;district_id&quot;: 365, &quot;district_name&quot;: &quot;khunti&quot;, &quot;pc11_pca_no_hh&quot;: null, &quot;shrid2&quot;: &quot;11-20-365-02705-375397&quot;, &quot;state_id&quot;: 20, &quot;state_name&quot;: &quot;jharkhand&quot;, &quot;subdistrict_id&quot;: 2705, &quot;tv_name&quot;: &quot;Banribera&quot;, &quot;urban&quot;: null, &quot;village_id&quot;: 375397}, &quot;type&quot;: &quot;Feature&quot;}], &quot;type&quot;: &quot;FeatureCollection&quot;});\n",
       "        geo_json_5c1eb2c2f07c487a3f7d86e2aed9b9e3.setStyle(function(feature) {return feature.properties.style;});\n",
       "\n",
       "        \n",
       "    \n",
       "    geo_json_5c1eb2c2f07c487a3f7d86e2aed9b9e3.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;tv_name&quot;];\n",
       "    let aliases = [&quot;tv_name&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_5c1eb2c2f07c487a3f7d86e2aed9b9e3.addTo(map_3c62ec64be52148a1a49e86e22e91445);\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 0x12fb7b950>"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# create a folium map for the rows with no census data\n",
    "ranchi_khunti_no_census = ranchi_khunti_census[ranchi_khunti_census['pc11_pca_no_hh'].isna()]\n",
    "\n",
    "def create_folium_map(gdf):\n",
    "    # get the center of the delhi map\n",
    "    center = gdf.centroid.iloc[0]\n",
    "\n",
    "    # create a folium map showing the villages in the buffer\n",
    "    m = folium.Map(location=[center.y, center.x], zoom_start=12)\n",
    "    folium.GeoJson(gdf, tooltip=folium.features.GeoJsonTooltip(fields=['tv_name'])).add_to(m)\n",
    "    return m\n",
    "\n",
    "m = create_folium_map(ranchi_khunti_no_census)\n",
    "m"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Sample SHRUG using PPS\n",
    "I sample WITH replacement since Ranchi is not broken up by ward so if I sampled without replacement Ranchi would be very underrepresented. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "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>state_id</th>\n",
       "      <th>district_id</th>\n",
       "      <th>subdistrict_id</th>\n",
       "      <th>village_id</th>\n",
       "      <th>tv_name</th>\n",
       "      <th>geometry</th>\n",
       "      <th>shrid2</th>\n",
       "      <th>state_name</th>\n",
       "      <th>district_name</th>\n",
       "      <th>pc11_pca_no_hh</th>\n",
       "      <th>urban</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>200</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2684</td>\n",
       "      <td>801794</td>\n",
       "      <td>Ranchi</td>\n",
       "      <td>POLYGON ((85.38199 23.43523, 85.38327 23.43388...</td>\n",
       "      <td>11-20-364-02684-801794</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>207636.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1871</th>\n",
       "      <td>20</td>\n",
       "      <td>365</td>\n",
       "      <td>2704</td>\n",
       "      <td>375236</td>\n",
       "      <td>Chinchal</td>\n",
       "      <td>POLYGON ((85.461 23.0822, 85.46087 23.08191, 8...</td>\n",
       "      <td>11-20-365-02704-375236</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>khunti</td>\n",
       "      <td>42.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>980</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2695</td>\n",
       "      <td>374346</td>\n",
       "      <td>Tero</td>\n",
       "      <td>POLYGON ((85.00861 23.33943, 85.00877 23.33827...</td>\n",
       "      <td>11-20-364-02695-374346</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>465.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>692</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2690</td>\n",
       "      <td>374058</td>\n",
       "      <td>Tati</td>\n",
       "      <td>POLYGON ((85.4482 23.3896, 85.44715 23.38811, ...</td>\n",
       "      <td>11-20-364-02690-374058</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>2544.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>200</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2684</td>\n",
       "      <td>801794</td>\n",
       "      <td>Ranchi</td>\n",
       "      <td>POLYGON ((85.38199 23.43523, 85.38327 23.43388...</td>\n",
       "      <td>11-20-364-02684-801794</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>207636.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>200</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2684</td>\n",
       "      <td>801794</td>\n",
       "      <td>Ranchi</td>\n",
       "      <td>POLYGON ((85.38199 23.43523, 85.38327 23.43388...</td>\n",
       "      <td>11-20-364-02684-801794</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>207636.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>115</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2684</td>\n",
       "      <td>373482</td>\n",
       "      <td>Uparkonki</td>\n",
       "      <td>POLYGON ((85.2784 23.5295, 85.27796 23.52918, ...</td>\n",
       "      <td>11-20-364-02684-373482</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>297.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1441</th>\n",
       "      <td>20</td>\n",
       "      <td>365</td>\n",
       "      <td>2700</td>\n",
       "      <td>374806</td>\n",
       "      <td>Mehan</td>\n",
       "      <td>POLYGON ((85.12808 23.17615, 85.12972 23.17556...</td>\n",
       "      <td>11-20-365-02700-374806</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>khunti</td>\n",
       "      <td>212.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>692</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2690</td>\n",
       "      <td>374058</td>\n",
       "      <td>Tati</td>\n",
       "      <td>POLYGON ((85.4482 23.3896, 85.44715 23.38811, ...</td>\n",
       "      <td>11-20-364-02690-374058</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>2544.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>903</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2694</td>\n",
       "      <td>374269</td>\n",
       "      <td>Kanjgi</td>\n",
       "      <td>POLYGON ((84.91115 23.42465, 84.90996 23.42472...</td>\n",
       "      <td>11-20-364-02694-374269</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>132.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      state_id  district_id  subdistrict_id  village_id    tv_name  \\\n",
       "200         20          364            2684      801794     Ranchi   \n",
       "1871        20          365            2704      375236   Chinchal   \n",
       "980         20          364            2695      374346       Tero   \n",
       "692         20          364            2690      374058       Tati   \n",
       "200         20          364            2684      801794     Ranchi   \n",
       "200         20          364            2684      801794     Ranchi   \n",
       "115         20          364            2684      373482  Uparkonki   \n",
       "1441        20          365            2700      374806      Mehan   \n",
       "692         20          364            2690      374058       Tati   \n",
       "903         20          364            2694      374269     Kanjgi   \n",
       "\n",
       "                                               geometry  \\\n",
       "200   POLYGON ((85.38199 23.43523, 85.38327 23.43388...   \n",
       "1871  POLYGON ((85.461 23.0822, 85.46087 23.08191, 8...   \n",
       "980   POLYGON ((85.00861 23.33943, 85.00877 23.33827...   \n",
       "692   POLYGON ((85.4482 23.3896, 85.44715 23.38811, ...   \n",
       "200   POLYGON ((85.38199 23.43523, 85.38327 23.43388...   \n",
       "200   POLYGON ((85.38199 23.43523, 85.38327 23.43388...   \n",
       "115   POLYGON ((85.2784 23.5295, 85.27796 23.52918, ...   \n",
       "1441  POLYGON ((85.12808 23.17615, 85.12972 23.17556...   \n",
       "692   POLYGON ((85.4482 23.3896, 85.44715 23.38811, ...   \n",
       "903   POLYGON ((84.91115 23.42465, 84.90996 23.42472...   \n",
       "\n",
       "                      shrid2 state_name district_name  pc11_pca_no_hh  urban  \n",
       "200   11-20-364-02684-801794  jharkhand        ranchi        207636.0    1.0  \n",
       "1871  11-20-365-02704-375236  jharkhand        khunti            42.0    0.0  \n",
       "980   11-20-364-02695-374346  jharkhand        ranchi           465.0    0.0  \n",
       "692   11-20-364-02690-374058  jharkhand        ranchi          2544.0    1.0  \n",
       "200   11-20-364-02684-801794  jharkhand        ranchi        207636.0    1.0  \n",
       "200   11-20-364-02684-801794  jharkhand        ranchi        207636.0    1.0  \n",
       "115   11-20-364-02684-373482  jharkhand        ranchi           297.0    0.0  \n",
       "1441  11-20-365-02700-374806  jharkhand        khunti           212.0    0.0  \n",
       "692   11-20-364-02690-374058  jharkhand        ranchi          2544.0    1.0  \n",
       "903   11-20-364-02694-374269  jharkhand        ranchi           132.0    0.0  "
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "probabilities = ranchi_khunti_census['pc11_pca_no_hh'] / ranchi_khunti_census['pc11_pca_no_hh'].sum()\n",
    "sampled_group = ranchi_khunti_census.sample(n=NUM_PSUS_TO_SAMPLE, weights=probabilities, replace=True, random_state=42)\n",
    "sampled_group"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "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>state_id</th>\n",
       "      <th>district_id</th>\n",
       "      <th>subdistrict_id</th>\n",
       "      <th>village_id</th>\n",
       "      <th>tv_name</th>\n",
       "      <th>geometry</th>\n",
       "      <th>shrid2</th>\n",
       "      <th>state_name</th>\n",
       "      <th>district_name</th>\n",
       "      <th>pc11_pca_no_hh</th>\n",
       "      <th>urban</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>200</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2684</td>\n",
       "      <td>801794</td>\n",
       "      <td>Ranchi</td>\n",
       "      <td>POLYGON ((85.38199 23.43523, 85.38327 23.43388...</td>\n",
       "      <td>11-20-364-02684-801794</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>207636.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1871</th>\n",
       "      <td>20</td>\n",
       "      <td>365</td>\n",
       "      <td>2704</td>\n",
       "      <td>375236</td>\n",
       "      <td>Chinchal</td>\n",
       "      <td>POLYGON ((85.461 23.0822, 85.46087 23.08191, 8...</td>\n",
       "      <td>11-20-365-02704-375236</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>khunti</td>\n",
       "      <td>42.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>980</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2695</td>\n",
       "      <td>374346</td>\n",
       "      <td>Tero</td>\n",
       "      <td>POLYGON ((85.00861 23.33943, 85.00877 23.33827...</td>\n",
       "      <td>11-20-364-02695-374346</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>465.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>692</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2690</td>\n",
       "      <td>374058</td>\n",
       "      <td>Tati</td>\n",
       "      <td>POLYGON ((85.4482 23.3896, 85.44715 23.38811, ...</td>\n",
       "      <td>11-20-364-02690-374058</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>2544.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>115</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2684</td>\n",
       "      <td>373482</td>\n",
       "      <td>Uparkonki</td>\n",
       "      <td>POLYGON ((85.2784 23.5295, 85.27796 23.52918, ...</td>\n",
       "      <td>11-20-364-02684-373482</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>297.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1441</th>\n",
       "      <td>20</td>\n",
       "      <td>365</td>\n",
       "      <td>2700</td>\n",
       "      <td>374806</td>\n",
       "      <td>Mehan</td>\n",
       "      <td>POLYGON ((85.12808 23.17615, 85.12972 23.17556...</td>\n",
       "      <td>11-20-365-02700-374806</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>khunti</td>\n",
       "      <td>212.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>903</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2694</td>\n",
       "      <td>374269</td>\n",
       "      <td>Kanjgi</td>\n",
       "      <td>POLYGON ((84.91115 23.42465, 84.90996 23.42472...</td>\n",
       "      <td>11-20-364-02694-374269</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>132.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      state_id  district_id  subdistrict_id  village_id    tv_name  \\\n",
       "200         20          364            2684      801794     Ranchi   \n",
       "1871        20          365            2704      375236   Chinchal   \n",
       "980         20          364            2695      374346       Tero   \n",
       "692         20          364            2690      374058       Tati   \n",
       "115         20          364            2684      373482  Uparkonki   \n",
       "1441        20          365            2700      374806      Mehan   \n",
       "903         20          364            2694      374269     Kanjgi   \n",
       "\n",
       "                                               geometry  \\\n",
       "200   POLYGON ((85.38199 23.43523, 85.38327 23.43388...   \n",
       "1871  POLYGON ((85.461 23.0822, 85.46087 23.08191, 8...   \n",
       "980   POLYGON ((85.00861 23.33943, 85.00877 23.33827...   \n",
       "692   POLYGON ((85.4482 23.3896, 85.44715 23.38811, ...   \n",
       "115   POLYGON ((85.2784 23.5295, 85.27796 23.52918, ...   \n",
       "1441  POLYGON ((85.12808 23.17615, 85.12972 23.17556...   \n",
       "903   POLYGON ((84.91115 23.42465, 84.90996 23.42472...   \n",
       "\n",
       "                      shrid2 state_name district_name  pc11_pca_no_hh  urban  \n",
       "200   11-20-364-02684-801794  jharkhand        ranchi        207636.0    1.0  \n",
       "1871  11-20-365-02704-375236  jharkhand        khunti            42.0    0.0  \n",
       "980   11-20-364-02695-374346  jharkhand        ranchi           465.0    0.0  \n",
       "692   11-20-364-02690-374058  jharkhand        ranchi          2544.0    1.0  \n",
       "115   11-20-364-02684-373482  jharkhand        ranchi           297.0    0.0  \n",
       "1441  11-20-365-02700-374806  jharkhand        khunti           212.0    0.0  \n",
       "903   11-20-364-02694-374269  jharkhand        ranchi           132.0    0.0  "
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# drop duplicates\n",
    "sampled_group_no_dupes = sampled_group.drop_duplicates(subset='shrid2')\n",
    "sampled_group_no_dupes"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Get rooftop data for these PSUs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\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>state_id</th>\n",
       "      <th>district_id</th>\n",
       "      <th>subdistrict_id</th>\n",
       "      <th>village_id</th>\n",
       "      <th>tv_name</th>\n",
       "      <th>geometry</th>\n",
       "      <th>shrid2</th>\n",
       "      <th>state_name</th>\n",
       "      <th>district_name</th>\n",
       "      <th>pc11_pca_no_hh</th>\n",
       "      <th>urban</th>\n",
       "      <th>s2_cell_id</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>200</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2684</td>\n",
       "      <td>801794</td>\n",
       "      <td>Ranchi</td>\n",
       "      <td>POLYGON ((85.38199 23.43523, 85.38327 23.43388...</td>\n",
       "      <td>11-20-364-02684-801794</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>207636.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>4176244229456003072</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1871</th>\n",
       "      <td>20</td>\n",
       "      <td>365</td>\n",
       "      <td>2704</td>\n",
       "      <td>375236</td>\n",
       "      <td>Chinchal</td>\n",
       "      <td>POLYGON ((85.461 23.0822, 85.46087 23.08191, 8...</td>\n",
       "      <td>11-20-365-02704-375236</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>khunti</td>\n",
       "      <td>42.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>4176244229456003072</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>980</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2695</td>\n",
       "      <td>374346</td>\n",
       "      <td>Tero</td>\n",
       "      <td>POLYGON ((85.00861 23.33943, 85.00877 23.33827...</td>\n",
       "      <td>11-20-364-02695-374346</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>465.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>4176244229456003072</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>692</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2690</td>\n",
       "      <td>374058</td>\n",
       "      <td>Tati</td>\n",
       "      <td>POLYGON ((85.4482 23.3896, 85.44715 23.38811, ...</td>\n",
       "      <td>11-20-364-02690-374058</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>2544.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>4176244229456003072</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>115</th>\n",
       "      <td>20</td>\n",
       "      <td>364</td>\n",
       "      <td>2684</td>\n",
       "      <td>373482</td>\n",
       "      <td>Uparkonki</td>\n",
       "      <td>POLYGON ((85.2784 23.5295, 85.27796 23.52918, ...</td>\n",
       "      <td>11-20-364-02684-373482</td>\n",
       "      <td>jharkhand</td>\n",
       "      <td>ranchi</td>\n",
       "      <td>297.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>4176244229456003072</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      state_id  district_id  subdistrict_id  village_id    tv_name  \\\n",
       "200         20          364            2684      801794     Ranchi   \n",
       "1871        20          365            2704      375236   Chinchal   \n",
       "980         20          364            2695      374346       Tero   \n",
       "692         20          364            2690      374058       Tati   \n",
       "115         20          364            2684      373482  Uparkonki   \n",
       "\n",
       "                                               geometry  \\\n",
       "200   POLYGON ((85.38199 23.43523, 85.38327 23.43388...   \n",
       "1871  POLYGON ((85.461 23.0822, 85.46087 23.08191, 8...   \n",
       "980   POLYGON ((85.00861 23.33943, 85.00877 23.33827...   \n",
       "692   POLYGON ((85.4482 23.3896, 85.44715 23.38811, ...   \n",
       "115   POLYGON ((85.2784 23.5295, 85.27796 23.52918, ...   \n",
       "\n",
       "                      shrid2 state_name district_name  pc11_pca_no_hh  urban  \\\n",
       "200   11-20-364-02684-801794  jharkhand        ranchi        207636.0    1.0   \n",
       "1871  11-20-365-02704-375236  jharkhand        khunti            42.0    0.0   \n",
       "980   11-20-364-02695-374346  jharkhand        ranchi           465.0    0.0   \n",
       "692   11-20-364-02690-374058  jharkhand        ranchi          2544.0    1.0   \n",
       "115   11-20-364-02684-373482  jharkhand        ranchi           297.0    0.0   \n",
       "\n",
       "               s2_cell_id  \n",
       "200   4176244229456003072  \n",
       "1871  4176244229456003072  \n",
       "980   4176244229456003072  \n",
       "692   4176244229456003072  \n",
       "115   4176244229456003072  "
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sampled_group_no_dupes['s2_cell_id'] = sampled_group_no_dupes.apply(lambda x: get_s2_cell_id(x.geometry.centroid), axis=1)\n",
    "sampled_group_no_dupes.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "# create empty gdf to store rooftops\n",
    "all_rooftops = gpd.GeoDataFrame()\n",
    "\n",
    "# since all the PSUs have the same s2_cell_id we can just use the first one\n",
    "s2_cell_id = sampled_group_no_dupes['s2_cell_id'].iloc[0]\n",
    "\n",
    "rooftops_gdf = gpd.read_parquet(ROOFTOP_DIR /f'{s2_cell_id}.parquet')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "."
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n",
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "."
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n",
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "."
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n",
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "."
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n",
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "."
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n",
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "."
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n",
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "."
     ]
    }
   ],
   "source": [
    "temp_rooftops = gpd.GeoDataFrame()\n",
    "all_rooftops = gpd.GeoDataFrame()\n",
    "\n",
    "for item, row in sampled_group_no_dupes.iterrows():\n",
    "    # print a single dot without the newline character\n",
    "    print('.', end='')\n",
    "    # filter rooftops to only include those within the psu\n",
    "    temp_rooftops = rooftops_gdf[rooftops_gdf.geometry.within(row.geometry)]\n",
    "    # set the psid for the rooftops\n",
    "    temp_rooftops['shrid2'] = row['shrid2']\n",
    "    temp_rooftops['urban'] = row['urban']\n",
    "    temp_rooftops.to_crs(epsg=4326, inplace=True)\n",
    "    # try to append temp_rooftops to all_rooftops and catch a value error. if there is an error, print the crs of the two dataframes\n",
    "    try:\n",
    "        all_rooftops = gpd.GeoDataFrame(pd.concat([all_rooftops, temp_rooftops], ignore_index=True))\n",
    "    except ValueError:\n",
    "        print(f\"Error with s2 cell {s2_cell_id}\")\n",
    "        print(f\"temp_rooftops crs: {temp_rooftops.crs}\")\n",
    "        print(f\"all_rooftops crs: {all_rooftops.crs}\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "# save all rooftops to a parquet file in case I close this notebook or the kernel gets messed up\n",
    "all_rooftops.to_parquet(OUTPUT_DIR / f'rooftops_sampled_psus_ranchi_khunti.parquet')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Sample 4 rooftops from each PSU (except Ranchi where I sample 12)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "all_rooftops['neighbors'] = count_neighbors_in_radius(all_rooftops)\n",
    "all_rooftops['isolated'] = (all_rooftops['neighbors'] < 10)\n",
    "\n",
    "# create dataframe of all non-isolated rooftops\n",
    "non_isolated_rooftops = all_rooftops[~all_rooftops['isolated']]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "ranchi_rooftops = non_isolated_rooftops[non_isolated_rooftops['shrid2'].str.contains('11-20-364-02684-801794')]\n",
    "other_rooftops =   non_isolated_rooftops[~non_isolated_rooftops['shrid2'].str.contains('11-20-364-02684-801794')]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/var/folders/_b/2wf3gjnx42s4k_vs5n079mh40000gn/T/ipykernel_57522/96762827.py:2: DeprecationWarning: DataFrameGroupBy.apply operated on the grouping columns. This behavior is deprecated, and in a future version of pandas the grouping columns will be excluded from the operation. Either pass `include_groups=False` to exclude the groupings or explicitly select the grouping columns after groupby to silence this warning.\n",
      "  sampled_points_other = other_rooftops.groupby('shrid2', group_keys=False).apply(lambda x: x.sample(n=min(4, x.shape[0]),random_state=42))\n"
     ]
    }
   ],
   "source": [
    "# sample 4 rooftops in each PSU for other rooftops\n",
    "sampled_points_other = other_rooftops.groupby('shrid2', group_keys=False).apply(lambda x: x.sample(n=min(4, x.shape[0]),random_state=42))\n",
    "\n",
    "# sample 12 rooftops in ranchi_rooftops\n",
    "sampled_points_ranchi = ranchi_rooftops.sample(n=12, random_state=42)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/var/folders/_b/2wf3gjnx42s4k_vs5n079mh40000gn/T/ipykernel_57522/3276395114.py:5: UserWarning: Geometry is in a geographic CRS. Results from 'centroid' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.\n",
      "\n",
      "  sampled_points['geometry'] = sampled_points.geometry.centroid\n"
     ]
    }
   ],
   "source": [
    "# concat the two dataframes\n",
    "sampled_points = pd.concat([sampled_points_other, sampled_points_ranchi])\n",
    "\n",
    "# get the centroid of each rooftop\n",
    "sampled_points['geometry'] = sampled_points.geometry.centroid"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Get nearest points on road "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "100%|██████████| 32/32 [00:04<00:00,  7.55it/s]\n"
     ]
    }
   ],
   "source": [
    "# for each point in the sampled points, get the nearest point on the road\n",
    "sampled_points['nearest_point_on_road'] = sampled_points.progress_apply(lambda x: get_nearest_point_on_road(x.geometry), axis=1)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Manual revision of points\n",
    "\n",
    "Some of the rural points that Google says lie on a road don't really seem to be on a road. For these, I set the \"nearest_point_on_road\" to null. I use the geohash for this because checking for equality of lat and long is really annoying because Google maps truncates the precision of the point."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "points_not_really_on_road = [\n",
    "    'tuhky4j4',\n",
    "    'tuhkwntc',\n",
    "    'tuhec1j9',\n",
    "    'tuhg23zu',\n",
    "    'tuhe9p4w',\n",
    "    'tuhg298p',\n",
    "    'tuhg2992' \n",
    "]\n",
    "\n",
    "# for points that are not really on the road, set nearest_point_on_road to None\n",
    "sampled_points.loc[sampled_points['geohash'].isin(points_not_really_on_road), 'nearest_point_on_road'] = None\n",
    "\n",
    "# remove one point that is in the airport\n",
    "sampled_points = sampled_points[sampled_points['geohash'] != 'tuhstt20']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "sampled_points['distance_to_road_m'] = sampled_points.apply(lambda x: dist_in_meters(x.geometry, x.nearest_point_on_road), axis=1)\n",
    "\n",
    "# create binary variable called \"on_road\" that is 1 if nearest_point_on_road is not null\n",
    "sampled_points['on_road'] = sampled_points['nearest_point_on_road'].notnull().astype(int)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Display map of points"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "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_949c60b3b945d722c556e79b10c74465 {\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",
       "&lt;/head&gt;\n",
       "&lt;body&gt;\n",
       "    \n",
       "    \n",
       "            &lt;div class=&quot;folium-map&quot; id=&quot;map_949c60b3b945d722c556e79b10c74465&quot; &gt;&lt;/div&gt;\n",
       "        \n",
       "&lt;/body&gt;\n",
       "&lt;script&gt;\n",
       "    \n",
       "    \n",
       "            var map_949c60b3b945d722c556e79b10c74465 = L.map(\n",
       "                &quot;map_949c60b3b945d722c556e79b10c74465&quot;,\n",
       "                {\n",
       "                    center: [23.516696608670312, 85.27140160918752],\n",
       "                    crs: L.CRS.EPSG3857,\n",
       "                    zoom: 12,\n",
       "                    zoomControl: true,\n",
       "                    preferCanvas: false,\n",
       "                }\n",
       "            );\n",
       "\n",
       "            \n",
       "\n",
       "        \n",
       "    \n",
       "            var tile_layer_28cfdd6c391cb6f578a63a92a89ffaa8 = L.tileLayer(\n",
       "                &quot;https://server.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer/tile/{z}/{y}/{x}&quot;,\n",
       "                {&quot;attribution&quot;: &quot;Esri&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_28cfdd6c391cb6f578a63a92a89ffaa8.addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var marker_ca6287c06488ebd58706ed2124021248 = L.marker(\n",
       "                [23.516696608670312, 85.27140160918752],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_ccb590f65aad5a154a0cc30e4766eb2d = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;green&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_ca6287c06488ebd58706ed2124021248.setIcon(icon_ccb590f65aad5a154a0cc30e4766eb2d);\n",
       "        \n",
       "    \n",
       "        var popup_9a399cc276813230830309d502312bff = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_2cb07efc00ef55b3fa22ef425a1d6c67 = $(`&lt;div id=&quot;html_2cb07efc00ef55b3fa22ef425a1d6c67&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_9a399cc276813230830309d502312bff.setContent(html_2cb07efc00ef55b3fa22ef425a1d6c67);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_ca6287c06488ebd58706ed2124021248.bindPopup(popup_9a399cc276813230830309d502312bff)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_875ab6e80bcbb23685fe3391089b805a = L.marker(\n",
       "                [23.513345590769276, 85.26910756164905],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_3fb0be6e733a98760b813aceb187bfe3 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;green&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_875ab6e80bcbb23685fe3391089b805a.setIcon(icon_3fb0be6e733a98760b813aceb187bfe3);\n",
       "        \n",
       "    \n",
       "        var popup_87803c9ff79bdef1ce40b96c659c2e89 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_162299b4d5464920fa587c25e6a4ae8e = $(`&lt;div id=&quot;html_162299b4d5464920fa587c25e6a4ae8e&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_87803c9ff79bdef1ce40b96c659c2e89.setContent(html_162299b4d5464920fa587c25e6a4ae8e);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_875ab6e80bcbb23685fe3391089b805a.bindPopup(popup_87803c9ff79bdef1ce40b96c659c2e89)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_197310c591a2b46e55a1849d92de4698 = L.marker(\n",
       "                [23.516188233204613, 85.27203417020021],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_bc477000080fa48cc80c296aeb725b09 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;green&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_197310c591a2b46e55a1849d92de4698.setIcon(icon_bc477000080fa48cc80c296aeb725b09);\n",
       "        \n",
       "    \n",
       "        var popup_f6474fd716943d5e101ef3e4b6733c05 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_3594aab4ba1b8d7d6354619c0bdc5e15 = $(`&lt;div id=&quot;html_3594aab4ba1b8d7d6354619c0bdc5e15&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_f6474fd716943d5e101ef3e4b6733c05.setContent(html_3594aab4ba1b8d7d6354619c0bdc5e15);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_197310c591a2b46e55a1849d92de4698.bindPopup(popup_f6474fd716943d5e101ef3e4b6733c05)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_2577bd84543680566d84a497d42a12ed = L.marker(\n",
       "                [23.516171774590465, 85.2705781306977],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_760b76d1ad872a49839a8a2425ead8d6 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;green&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_2577bd84543680566d84a497d42a12ed.setIcon(icon_760b76d1ad872a49839a8a2425ead8d6);\n",
       "        \n",
       "    \n",
       "        var popup_f74fb4bbc1c080c627ec7b7d18b3d376 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_2e5a28d624bf82a829986c1046dc3491 = $(`&lt;div id=&quot;html_2e5a28d624bf82a829986c1046dc3491&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_f74fb4bbc1c080c627ec7b7d18b3d376.setContent(html_2e5a28d624bf82a829986c1046dc3491);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_2577bd84543680566d84a497d42a12ed.bindPopup(popup_f74fb4bbc1c080c627ec7b7d18b3d376)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_0956bf3b7148a4fbf0828a4368f7f2b5 = L.marker(\n",
       "                [23.372241298416583, 85.43897965725373],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_e11dc490d2a05e4e21a85499b90c3458 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_0956bf3b7148a4fbf0828a4368f7f2b5.setIcon(icon_e11dc490d2a05e4e21a85499b90c3458);\n",
       "        \n",
       "    \n",
       "        var popup_619a360df085e79fc3514411848c24fa = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_287a83828475ad5144c61e2108738f85 = $(`&lt;div id=&quot;html_287a83828475ad5144c61e2108738f85&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_619a360df085e79fc3514411848c24fa.setContent(html_287a83828475ad5144c61e2108738f85);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_0956bf3b7148a4fbf0828a4368f7f2b5.bindPopup(popup_619a360df085e79fc3514411848c24fa)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_32fd52bebb78e9395ff10468b8fe2a28 = L.marker(\n",
       "                [23.373972109571895, 85.43999555353629],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_c6e11cc57524a838480ee5841aaa3e59 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_32fd52bebb78e9395ff10468b8fe2a28.setIcon(icon_c6e11cc57524a838480ee5841aaa3e59);\n",
       "        \n",
       "    \n",
       "        var popup_72bec4515c0494e2e30c8bcdf0bf0b78 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_21d9b0635dd812f6389348ed5c4d3421 = $(`&lt;div id=&quot;html_21d9b0635dd812f6389348ed5c4d3421&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_72bec4515c0494e2e30c8bcdf0bf0b78.setContent(html_21d9b0635dd812f6389348ed5c4d3421);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_32fd52bebb78e9395ff10468b8fe2a28.bindPopup(popup_72bec4515c0494e2e30c8bcdf0bf0b78)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_b851a903e3ce0824c68f78c1069291d4 = L.marker(\n",
       "                [23.374838677462858, 85.43339249620735],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_5bc3a4c366a59e168ace45af168b925f = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_b851a903e3ce0824c68f78c1069291d4.setIcon(icon_5bc3a4c366a59e168ace45af168b925f);\n",
       "        \n",
       "    \n",
       "        var popup_ea0e5c89a2ea66aedb1069d86fe287bb = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_837900f773e469721d47c8744bdf7f7b = $(`&lt;div id=&quot;html_837900f773e469721d47c8744bdf7f7b&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_ea0e5c89a2ea66aedb1069d86fe287bb.setContent(html_837900f773e469721d47c8744bdf7f7b);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_b851a903e3ce0824c68f78c1069291d4.bindPopup(popup_ea0e5c89a2ea66aedb1069d86fe287bb)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_29d11ef08b6a475a6bf02c6d0162cb10 = L.marker(\n",
       "                [23.377310439024242, 85.42472082350446],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_d9b27ca32912f1d274dbdbffb3dbdf5d = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_29d11ef08b6a475a6bf02c6d0162cb10.setIcon(icon_d9b27ca32912f1d274dbdbffb3dbdf5d);\n",
       "        \n",
       "    \n",
       "        var popup_c254e6c8834da7714c9b5d95725edafa = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_acfaa860cebdeb3feda51e6f342754e1 = $(`&lt;div id=&quot;html_acfaa860cebdeb3feda51e6f342754e1&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_c254e6c8834da7714c9b5d95725edafa.setContent(html_acfaa860cebdeb3feda51e6f342754e1);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_29d11ef08b6a475a6bf02c6d0162cb10.bindPopup(popup_c254e6c8834da7714c9b5d95725edafa)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_adbe0adc9e4eaf29948c32d8c675f6a1 = L.marker(\n",
       "                [23.336287073080115, 85.00040518638738],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_96e140bdd0876fa32a94863b487a5762 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;green&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_adbe0adc9e4eaf29948c32d8c675f6a1.setIcon(icon_96e140bdd0876fa32a94863b487a5762);\n",
       "        \n",
       "    \n",
       "        var popup_238d09424f3ce965eef1c46cb7a25499 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_7cd6daf36335176e6c229328eb9fe90d = $(`&lt;div id=&quot;html_7cd6daf36335176e6c229328eb9fe90d&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_238d09424f3ce965eef1c46cb7a25499.setContent(html_7cd6daf36335176e6c229328eb9fe90d);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_adbe0adc9e4eaf29948c32d8c675f6a1.bindPopup(popup_238d09424f3ce965eef1c46cb7a25499)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_fd2af0e74e2cdb7e1b9ab3542ca9dfdc = L.marker(\n",
       "                [23.346431044894366, 84.9974142806709],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_c3e97f66a21fe9ce16381c45ee133bc9 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;gray&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_fd2af0e74e2cdb7e1b9ab3542ca9dfdc.setIcon(icon_c3e97f66a21fe9ce16381c45ee133bc9);\n",
       "        \n",
       "    \n",
       "        var popup_20b5127e120f3b00650d77c17ff360d3 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_e2c3c3064b2a9eadbba3b1c90079fe90 = $(`&lt;div id=&quot;html_e2c3c3064b2a9eadbba3b1c90079fe90&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 0&lt;/div&gt;`)[0];\n",
       "                popup_20b5127e120f3b00650d77c17ff360d3.setContent(html_e2c3c3064b2a9eadbba3b1c90079fe90);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_fd2af0e74e2cdb7e1b9ab3542ca9dfdc.bindPopup(popup_20b5127e120f3b00650d77c17ff360d3)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_383e5c5d289e2efcfa21fe3f29299f9e = L.marker(\n",
       "                [23.326990651506627, 84.99831989013455],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_9d71ad4c0eb3576bac0a72d503e3616e = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;gray&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_383e5c5d289e2efcfa21fe3f29299f9e.setIcon(icon_9d71ad4c0eb3576bac0a72d503e3616e);\n",
       "        \n",
       "    \n",
       "        var popup_ec3f56897d0d6d80ed31cb083ed6c022 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_074d2bb961fd48f74ed99a10bb569845 = $(`&lt;div id=&quot;html_074d2bb961fd48f74ed99a10bb569845&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 0&lt;/div&gt;`)[0];\n",
       "                popup_ec3f56897d0d6d80ed31cb083ed6c022.setContent(html_074d2bb961fd48f74ed99a10bb569845);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_383e5c5d289e2efcfa21fe3f29299f9e.bindPopup(popup_ec3f56897d0d6d80ed31cb083ed6c022)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_7db0cd8a615e4b50614ff3448ed5ad3e = L.marker(\n",
       "                [23.33341009062788, 84.9953755216113],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_bb34a01ae8cce81d5b1e7a5cec5aa52c = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;green&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_7db0cd8a615e4b50614ff3448ed5ad3e.setIcon(icon_bb34a01ae8cce81d5b1e7a5cec5aa52c);\n",
       "        \n",
       "    \n",
       "        var popup_5e0920d3fb3ed8241e5a37d5d93b0236 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_fcafb59d7309da28c43dcb9068e564c6 = $(`&lt;div id=&quot;html_fcafb59d7309da28c43dcb9068e564c6&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_5e0920d3fb3ed8241e5a37d5d93b0236.setContent(html_fcafb59d7309da28c43dcb9068e564c6);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_7db0cd8a615e4b50614ff3448ed5ad3e.bindPopup(popup_5e0920d3fb3ed8241e5a37d5d93b0236)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_e36616eefc1389fd9b93fd74491adc59 = L.marker(\n",
       "                [23.166867948215128, 85.13220109702206],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_929f22f28cecefddc596fda300c8d299 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;gray&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_e36616eefc1389fd9b93fd74491adc59.setIcon(icon_929f22f28cecefddc596fda300c8d299);\n",
       "        \n",
       "    \n",
       "        var popup_bf4ccf80348246beb938511301e0fe9c = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_7729bad04e875b8e05c9fc90e859a3fe = $(`&lt;div id=&quot;html_7729bad04e875b8e05c9fc90e859a3fe&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 0&lt;/div&gt;`)[0];\n",
       "                popup_bf4ccf80348246beb938511301e0fe9c.setContent(html_7729bad04e875b8e05c9fc90e859a3fe);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_e36616eefc1389fd9b93fd74491adc59.bindPopup(popup_bf4ccf80348246beb938511301e0fe9c)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_8f086fbe3c421192db7f79dc6b3d2d7e = L.marker(\n",
       "                [23.164953659628566, 85.12970232016603],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_be9f64c4d47ec919d477c4a4e4c97cd1 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;gray&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_8f086fbe3c421192db7f79dc6b3d2d7e.setIcon(icon_be9f64c4d47ec919d477c4a4e4c97cd1);\n",
       "        \n",
       "    \n",
       "        var popup_a505316aed8aeb1daac97fa751c8a745 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_d6b5614bac2a1fed8c4c9315fd6759bf = $(`&lt;div id=&quot;html_d6b5614bac2a1fed8c4c9315fd6759bf&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 0&lt;/div&gt;`)[0];\n",
       "                popup_a505316aed8aeb1daac97fa751c8a745.setContent(html_d6b5614bac2a1fed8c4c9315fd6759bf);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_8f086fbe3c421192db7f79dc6b3d2d7e.bindPopup(popup_a505316aed8aeb1daac97fa751c8a745)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_a3d4cfd70f41608ad5edd78c0620d5ef = L.marker(\n",
       "                [23.154841289925145, 85.12571620642944],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_1a51881fa6e9c044bb4878b014926536 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;gray&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_a3d4cfd70f41608ad5edd78c0620d5ef.setIcon(icon_1a51881fa6e9c044bb4878b014926536);\n",
       "        \n",
       "    \n",
       "        var popup_64480e554d5990cd98c915fcf5ab1c97 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_a3f8161f557007b2ef9b198f5da91766 = $(`&lt;div id=&quot;html_a3f8161f557007b2ef9b198f5da91766&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 0&lt;/div&gt;`)[0];\n",
       "                popup_64480e554d5990cd98c915fcf5ab1c97.setContent(html_a3f8161f557007b2ef9b198f5da91766);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_a3d4cfd70f41608ad5edd78c0620d5ef.bindPopup(popup_64480e554d5990cd98c915fcf5ab1c97)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_4066375e552e60f64c4302b4a2b9ddd7 = L.marker(\n",
       "                [23.16662239064889, 85.13615745933186],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_81e4e56e67e77425731885d4230d4fe2 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;gray&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_4066375e552e60f64c4302b4a2b9ddd7.setIcon(icon_81e4e56e67e77425731885d4230d4fe2);\n",
       "        \n",
       "    \n",
       "        var popup_482f13f349f29f1a54db0367bc5fdd44 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_6d0ce4e9d206d26fc842c4a7e995ea4a = $(`&lt;div id=&quot;html_6d0ce4e9d206d26fc842c4a7e995ea4a&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 0&lt;/div&gt;`)[0];\n",
       "                popup_482f13f349f29f1a54db0367bc5fdd44.setContent(html_6d0ce4e9d206d26fc842c4a7e995ea4a);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_4066375e552e60f64c4302b4a2b9ddd7.bindPopup(popup_482f13f349f29f1a54db0367bc5fdd44)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_a433fcbe1ab8e09ae264036b8d6349ec = L.marker(\n",
       "                [23.08163893878751, 85.45154714394273],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_8951b60450daa437b759ae04260cf333 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;gray&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_a433fcbe1ab8e09ae264036b8d6349ec.setIcon(icon_8951b60450daa437b759ae04260cf333);\n",
       "        \n",
       "    \n",
       "        var popup_c6842acdb281299b5cb59faff2b45525 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_7e7a6df0d3f9b4b312899715441bee2b = $(`&lt;div id=&quot;html_7e7a6df0d3f9b4b312899715441bee2b&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 0&lt;/div&gt;`)[0];\n",
       "                popup_c6842acdb281299b5cb59faff2b45525.setContent(html_7e7a6df0d3f9b4b312899715441bee2b);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_a433fcbe1ab8e09ae264036b8d6349ec.bindPopup(popup_c6842acdb281299b5cb59faff2b45525)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_a69963506976efb16c303cdc3e0538ad = L.marker(\n",
       "                [23.08080603817383, 85.45175065128267],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_0f1fc5990845d8c98a81507570e5ee12 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;gray&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_a69963506976efb16c303cdc3e0538ad.setIcon(icon_0f1fc5990845d8c98a81507570e5ee12);\n",
       "        \n",
       "    \n",
       "        var popup_b2ecaab9513fa124463f7f6a91c700b2 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_631267ed06ca1647b376ed3bfd43852b = $(`&lt;div id=&quot;html_631267ed06ca1647b376ed3bfd43852b&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 0&lt;/div&gt;`)[0];\n",
       "                popup_b2ecaab9513fa124463f7f6a91c700b2.setContent(html_631267ed06ca1647b376ed3bfd43852b);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_a69963506976efb16c303cdc3e0538ad.bindPopup(popup_b2ecaab9513fa124463f7f6a91c700b2)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_f0f4c74b42ae1a6370d0f2ffa1a7cac6 = L.marker(\n",
       "                [23.079334620991418, 85.44561403673399],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_38faba45497a34219f84206a7b244520 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;gray&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_f0f4c74b42ae1a6370d0f2ffa1a7cac6.setIcon(icon_38faba45497a34219f84206a7b244520);\n",
       "        \n",
       "    \n",
       "        var popup_23e3924b6470b7a69c41c7f7d58012fc = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_892656377560b6fcf6480f709360d8f6 = $(`&lt;div id=&quot;html_892656377560b6fcf6480f709360d8f6&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 0&lt;/div&gt;`)[0];\n",
       "                popup_23e3924b6470b7a69c41c7f7d58012fc.setContent(html_892656377560b6fcf6480f709360d8f6);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_f0f4c74b42ae1a6370d0f2ffa1a7cac6.bindPopup(popup_23e3924b6470b7a69c41c7f7d58012fc)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_80a54b70947c439c6fd2ad7453e9b515 = L.marker(\n",
       "                [23.07969236302457, 85.45347161228239],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_d8647b83550b94ead6b08b88991003e9 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;gray&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_80a54b70947c439c6fd2ad7453e9b515.setIcon(icon_d8647b83550b94ead6b08b88991003e9);\n",
       "        \n",
       "    \n",
       "        var popup_78a10dd09b9ecc8f3aedf23f75584946 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_d6f5fce8c229107b587eca55cacff2ea = $(`&lt;div id=&quot;html_d6f5fce8c229107b587eca55cacff2ea&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 0.0, On Road: 0&lt;/div&gt;`)[0];\n",
       "                popup_78a10dd09b9ecc8f3aedf23f75584946.setContent(html_d6f5fce8c229107b587eca55cacff2ea);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_80a54b70947c439c6fd2ad7453e9b515.bindPopup(popup_78a10dd09b9ecc8f3aedf23f75584946)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_58d1db4dd581976f0cdf55ce249a95d0 = L.marker(\n",
       "                [23.299515870569923, 85.31872650736022],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_9bfe2e9b49550c52209606137032f42d = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_58d1db4dd581976f0cdf55ce249a95d0.setIcon(icon_9bfe2e9b49550c52209606137032f42d);\n",
       "        \n",
       "    \n",
       "        var popup_dbff41274d9e55e9d09a8886752ecac8 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_bd1be09c9a6d387dd00043684b8e2e93 = $(`&lt;div id=&quot;html_bd1be09c9a6d387dd00043684b8e2e93&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_dbff41274d9e55e9d09a8886752ecac8.setContent(html_bd1be09c9a6d387dd00043684b8e2e93);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_58d1db4dd581976f0cdf55ce249a95d0.bindPopup(popup_dbff41274d9e55e9d09a8886752ecac8)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_3b7790bd50c6e3cf5939e27633c987d1 = L.marker(\n",
       "                [23.403447205226765, 85.3534967753574],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_3b44a66a7fc6918a51e7a212a661b1b0 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_3b7790bd50c6e3cf5939e27633c987d1.setIcon(icon_3b44a66a7fc6918a51e7a212a661b1b0);\n",
       "        \n",
       "    \n",
       "        var popup_922a3a90a529fe72c370169a3a7a0515 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_c08427c0f793b7fa2fd5fe6d2e0a2f3a = $(`&lt;div id=&quot;html_c08427c0f793b7fa2fd5fe6d2e0a2f3a&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_922a3a90a529fe72c370169a3a7a0515.setContent(html_c08427c0f793b7fa2fd5fe6d2e0a2f3a);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_3b7790bd50c6e3cf5939e27633c987d1.bindPopup(popup_922a3a90a529fe72c370169a3a7a0515)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_4b2e336bf3b272427685633d5d583567 = L.marker(\n",
       "                [23.3834366695394, 85.3021144665907],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_ce40f7e22f2d0570521fd3d053f8b915 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_4b2e336bf3b272427685633d5d583567.setIcon(icon_ce40f7e22f2d0570521fd3d053f8b915);\n",
       "        \n",
       "    \n",
       "        var popup_9e0e0b88b52900733a6b3e37072842e5 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_4f8085dc479142812d3f7a6709849f37 = $(`&lt;div id=&quot;html_4f8085dc479142812d3f7a6709849f37&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_9e0e0b88b52900733a6b3e37072842e5.setContent(html_4f8085dc479142812d3f7a6709849f37);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_4b2e336bf3b272427685633d5d583567.bindPopup(popup_9e0e0b88b52900733a6b3e37072842e5)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_b794298ef41bb271b3e0bbbf65202400 = L.marker(\n",
       "                [23.40888574880101, 85.31579577594854],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_8e54d8288ca14fa5b4cabeaaf5f0bbaf = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_b794298ef41bb271b3e0bbbf65202400.setIcon(icon_8e54d8288ca14fa5b4cabeaaf5f0bbaf);\n",
       "        \n",
       "    \n",
       "        var popup_421a9fcaf136c02ab8a6fe45b501d923 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_85657a3443862ba06f1c51a78c570d2f = $(`&lt;div id=&quot;html_85657a3443862ba06f1c51a78c570d2f&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_421a9fcaf136c02ab8a6fe45b501d923.setContent(html_85657a3443862ba06f1c51a78c570d2f);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_b794298ef41bb271b3e0bbbf65202400.bindPopup(popup_421a9fcaf136c02ab8a6fe45b501d923)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_c1e63f41d1e63533dcd9f5cf2d6b691b = L.marker(\n",
       "                [23.36787142768376, 85.32240039125793],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_a36b198af6a85e6aba4b3aefc082d31a = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_c1e63f41d1e63533dcd9f5cf2d6b691b.setIcon(icon_a36b198af6a85e6aba4b3aefc082d31a);\n",
       "        \n",
       "    \n",
       "        var popup_fdf3c6e099b38ee16bbfd5dfb7f5ea5d = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_8467ad7271b2dd5cd50bfc282ffad7d9 = $(`&lt;div id=&quot;html_8467ad7271b2dd5cd50bfc282ffad7d9&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_fdf3c6e099b38ee16bbfd5dfb7f5ea5d.setContent(html_8467ad7271b2dd5cd50bfc282ffad7d9);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_c1e63f41d1e63533dcd9f5cf2d6b691b.bindPopup(popup_fdf3c6e099b38ee16bbfd5dfb7f5ea5d)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_13c7edd6bb00641e5a243e78af28c43f = L.marker(\n",
       "                [23.371722793370505, 85.31780960668829],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_1aaddf931c661800de2bb3240a236400 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_13c7edd6bb00641e5a243e78af28c43f.setIcon(icon_1aaddf931c661800de2bb3240a236400);\n",
       "        \n",
       "    \n",
       "        var popup_18645423b2481efe9d0e3bfdf7375903 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_238740bc8303e70eac3221936930b40a = $(`&lt;div id=&quot;html_238740bc8303e70eac3221936930b40a&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_18645423b2481efe9d0e3bfdf7375903.setContent(html_238740bc8303e70eac3221936930b40a);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_13c7edd6bb00641e5a243e78af28c43f.bindPopup(popup_18645423b2481efe9d0e3bfdf7375903)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_7612904fdff3201e75d8721e1423a059 = L.marker(\n",
       "                [23.29700761291445, 85.27580875027589],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_bebb3c1cba0f2f6fe5c6c16549249fb0 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_7612904fdff3201e75d8721e1423a059.setIcon(icon_bebb3c1cba0f2f6fe5c6c16549249fb0);\n",
       "        \n",
       "    \n",
       "        var popup_ee86a19dc2fd6dcf17626d30ee66e68d = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_a8e2fedb50954608fc8242f31c3d0a26 = $(`&lt;div id=&quot;html_a8e2fedb50954608fc8242f31c3d0a26&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_ee86a19dc2fd6dcf17626d30ee66e68d.setContent(html_a8e2fedb50954608fc8242f31c3d0a26);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_7612904fdff3201e75d8721e1423a059.bindPopup(popup_ee86a19dc2fd6dcf17626d30ee66e68d)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_65aee89fe20cefb268634e08c3e1b7c4 = L.marker(\n",
       "                [23.381000156327076, 85.30650260308568],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_73f293264b0e71f6498033d3616901b4 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_65aee89fe20cefb268634e08c3e1b7c4.setIcon(icon_73f293264b0e71f6498033d3616901b4);\n",
       "        \n",
       "    \n",
       "        var popup_dc72514afff785304262daf9287aa262 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_888a9f1ff12b6b304be717fc54d533b5 = $(`&lt;div id=&quot;html_888a9f1ff12b6b304be717fc54d533b5&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_dc72514afff785304262daf9287aa262.setContent(html_888a9f1ff12b6b304be717fc54d533b5);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_65aee89fe20cefb268634e08c3e1b7c4.bindPopup(popup_dc72514afff785304262daf9287aa262)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_f91522c30a107a53acad3b8b3b3f5d73 = L.marker(\n",
       "                [23.33849146390648, 85.32708675945723],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_0d0b79cee6443bca3e65eeca2c3d7c4f = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_f91522c30a107a53acad3b8b3b3f5d73.setIcon(icon_0d0b79cee6443bca3e65eeca2c3d7c4f);\n",
       "        \n",
       "    \n",
       "        var popup_ecbeb8dd347b90d01c469e7baf18b809 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_1ccc4e5c9211396783e472e104013f0b = $(`&lt;div id=&quot;html_1ccc4e5c9211396783e472e104013f0b&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_ecbeb8dd347b90d01c469e7baf18b809.setContent(html_1ccc4e5c9211396783e472e104013f0b);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_f91522c30a107a53acad3b8b3b3f5d73.bindPopup(popup_ecbeb8dd347b90d01c469e7baf18b809)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_1590b37aa54b80cd4c28b391e06c81f5 = L.marker(\n",
       "                [23.336064373449627, 85.30800989821344],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_05ee732ae53b58c5824ad138ba477364 = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_1590b37aa54b80cd4c28b391e06c81f5.setIcon(icon_05ee732ae53b58c5824ad138ba477364);\n",
       "        \n",
       "    \n",
       "        var popup_df6e583c2b5a1e3081f38e7472413f1c = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_9880b2913441da3da7e04e2d76fc07c7 = $(`&lt;div id=&quot;html_9880b2913441da3da7e04e2d76fc07c7&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_df6e583c2b5a1e3081f38e7472413f1c.setContent(html_9880b2913441da3da7e04e2d76fc07c7);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_1590b37aa54b80cd4c28b391e06c81f5.bindPopup(popup_df6e583c2b5a1e3081f38e7472413f1c)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            var marker_e1b4e9315cbf478748baaa2bf9e8c6cf = L.marker(\n",
       "                [23.318979900644862, 85.29468369629942],\n",
       "                {}\n",
       "            ).addTo(map_949c60b3b945d722c556e79b10c74465);\n",
       "        \n",
       "    \n",
       "            var icon_6d1cdd1deb95cdd57a6522e59b9f222f = L.AwesomeMarkers.icon(\n",
       "                {&quot;extraClasses&quot;: &quot;fa-rotate-0&quot;, &quot;icon&quot;: &quot;info-sign&quot;, &quot;iconColor&quot;: &quot;white&quot;, &quot;markerColor&quot;: &quot;red&quot;, &quot;prefix&quot;: &quot;glyphicon&quot;}\n",
       "            );\n",
       "            marker_e1b4e9315cbf478748baaa2bf9e8c6cf.setIcon(icon_6d1cdd1deb95cdd57a6522e59b9f222f);\n",
       "        \n",
       "    \n",
       "        var popup_c2cd64bfb868052c7c8d169fc4de6257 = L.popup({&quot;maxWidth&quot;: &quot;100%&quot;});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_3a27b01ef31b679a045b8367c8463cb3 = $(`&lt;div id=&quot;html_3a27b01ef31b679a045b8367c8463cb3&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;Urban: 1.0, On Road: 1&lt;/div&gt;`)[0];\n",
       "                popup_c2cd64bfb868052c7c8d169fc4de6257.setContent(html_3a27b01ef31b679a045b8367c8463cb3);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_e1b4e9315cbf478748baaa2bf9e8c6cf.bindPopup(popup_c2cd64bfb868052c7c8d169fc4de6257)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            tile_layer_28cfdd6c391cb6f578a63a92a89ffaa8.addTo(map_949c60b3b945d722c556e79b10c74465);\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 0x168977090>"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# create folium map of points\n",
    "center = sampled_points.geometry.iloc[0]\n",
    "\n",
    "# create a folium map showing the villages in the buffer\n",
    "m = folium.Map(location=[center.y, center.x], zoom_start=12, tiles='https://server.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer/tile/{z}/{y}/{x}', attr='Esri')\n",
    "\n",
    "# Function to determine marker color based on 'urban' and 'on_road' values\n",
    "def get_marker_color(urban, on_road):\n",
    "    if urban and on_road:\n",
    "        return 'red'\n",
    "    elif urban and not on_road:\n",
    "        return 'blue'\n",
    "    elif not urban and on_road:\n",
    "        return 'green'\n",
    "    else:\n",
    "        return 'gray'\n",
    "\n",
    "# Iterate through the GeoDataFrame and add points to the map\n",
    "for idx, row in sampled_points.iterrows():\n",
    "    folium.Marker(\n",
    "        location=[row.geometry.y, row.geometry.x],\n",
    "        popup=f\"Urban: {row['urban']}, On Road: {row['on_road']}\",\n",
    "        icon=folium.Icon(color=get_marker_color(row['urban'], row['on_road']))\n",
    "    ).add_to(m)\n",
    "\n",
    "# save the map as sampled points\n",
    "m.save(OUTPUT_DIR / 'sampled_points.html')\n",
    "m"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Save sample outputs as CSV"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n",
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n",
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n",
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n",
      "/var/folders/_b/2wf3gjnx42s4k_vs5n079mh40000gn/T/ipykernel_57522/524848326.py:14: FutureWarning: Series.__getitem__ treating keys as positions is deprecated. In a future version, integer keys will always be treated as labels (consistent with DataFrame behavior). To access a value by position, use `ser.iloc[pos]`\n",
      "  ].apply(lambda x: f\"https://www.google.com/maps?q={x[0]},{x[1]}\", axis=1)\n"
     ]
    }
   ],
   "source": [
    "temp_samp_points = sampled_points[['shrid2', 'geometry', 'nearest_point_on_road', 'distance_to_road_m', 'urban']]\n",
    "temp_samp_points['orig_lat'] = temp_samp_points['geometry'].y\n",
    "temp_samp_points['orig_lon'] = temp_samp_points['geometry'].x\n",
    "temp_samp_points['road_lat'] = temp_samp_points['nearest_point_on_road'].y\n",
    "temp_samp_points['road_lon'] = temp_samp_points['nearest_point_on_road'].x\n",
    "temp_samp_points = temp_samp_points.drop(columns=['geometry', 'nearest_point_on_road'])\n",
    "temp_samp_points['google_maps_directions_link'] = temp_samp_points.apply(\n",
    "    lambda row: f\"https://www.google.com/maps/dir/{row['orig_lat']},{row['orig_lon']}/{row['road_lat']},{row['road_lon']}\",\n",
    "    axis=1\n",
    ")\n",
    "# replace google_maps_directions_link with a link to just the original point if nearest_point_on_road is null\n",
    "temp_samp_points.loc[temp_samp_points['road_lat'].isna(), 'google_maps_directions_link'] = temp_samp_points.loc[\n",
    "    temp_samp_points['road_lat'].isna(), ['orig_lat', 'orig_lon']\n",
    "].apply(lambda x: f\"https://www.google.com/maps?q={x[0]},{x[1]}\", axis=1)\n",
    "\n",
    "\n",
    "# save the sampled points as csv\n",
    "temp_samp_points.to_csv(OUTPUT_DIR / f'samp_points.csv', index=False)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Save points as kmls\n",
    "I create several different kml files:\n",
    "\n",
    "1. a kml files of lines between the original point and the nearest point on road\n",
    "2. a kml file of points on road that are urban\n",
    "3. a kml file of points on road that are rural\n",
    "4. a kml file of points NOT on road (all are rural)\n",
    "\n",
    "I then upload each of these kml files to a single Google My Map separately. While possible to combine all of these in python it is a bit of a pain.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/geopandas/geodataframe.py:1819: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  super().__setitem__(key, value)\n",
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/pyogrio/geopandas.py:662: UserWarning: 'crs' was not provided.  The output dataset will not have projection information defined and may not be usable in other systems.\n",
      "  write(\n",
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/pyogrio/geopandas.py:662: UserWarning: 'crs' was not provided.  The output dataset will not have projection information defined and may not be usable in other systems.\n",
      "  write(\n",
      "/Users/douglasjohnson/Documents/code/pin_drop_sampling2/.venv/lib/python3.11/site-packages/pyogrio/geopandas.py:662: UserWarning: 'crs' was not provided.  The output dataset will not have projection information defined and may not be usable in other systems.\n",
      "  write(\n"
     ]
    }
   ],
   "source": [
    "sampled_points_on_road = sampled_points[~(sampled_points['nearest_point_on_road'].isnull())]\n",
    "sampled_points_on_road_urban = sampled_points_on_road[sampled_points_on_road['urban'] == 1]\n",
    "sampled_points_on_road_rural = sampled_points_on_road[sampled_points_on_road['urban'] == 0]\n",
    "sampled_points_off_road = sampled_points[sampled_points['nearest_point_on_road'].isnull()]\n",
    "\n",
    "# save lines between original points and nearest points on road\n",
    "sampled_points_on_road['line'] = sampled_points_on_road.apply(lambda row: LineString([row['geometry'], row['nearest_point_on_road']]), axis=1)\n",
    "lines_gdf = sampled_points_on_road[['line', 'shrid2']]\n",
    "lines_gdf.to_file(OUTPUT_DIR / 'lines.kml', driver='KML')\n",
    "\n",
    "# save sampled points on road and urban as kml\n",
    "sampled_points_on_road_urban[['nearest_point_on_road', 'shrid2']].to_file(OUTPUT_DIR / 'points_on_road_urban.kml', driver='KML')\n",
    "\n",
    "# save sampled points on road and rural as kml\n",
    "sampled_points_on_road_rural[['nearest_point_on_road', 'shrid2']].to_file(OUTPUT_DIR / 'points_on_road_rural.kml', driver='KML')\n",
    "\n",
    "# save sampled points off road as kml\n",
    "sampled_points_off_road[['geometry', 'shrid2']].to_file(OUTPUT_DIR / 'points_off_road.kml', driver='KML')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "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>boundary_id</th>\n",
       "      <th>bf_source</th>\n",
       "      <th>confidence</th>\n",
       "      <th>area_in_meters</th>\n",
       "      <th>s2_id</th>\n",
       "      <th>country_iso</th>\n",
       "      <th>geometry</th>\n",
       "      <th>geohash</th>\n",
       "      <th>shrid2</th>\n",
       "      <th>urban</th>\n",
       "      <th>neighbors</th>\n",
       "      <th>isolated</th>\n",
       "      <th>nearest_point_on_road</th>\n",
       "      <th>distance_to_road_m</th>\n",
       "      <th>on_road</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>347252</th>\n",
       "      <td>56</td>\n",
       "      <td>google</td>\n",
       "      <td>0.8294</td>\n",
       "      <td>32.1143</td>\n",
       "      <td>4176244229456003072</td>\n",
       "      <td>IND</td>\n",
       "      <td>POINT (85.2714 23.5167)</td>\n",
       "      <td>tuhtu3hd</td>\n",
       "      <td>11-20-364-02684-373482</td>\n",
       "      <td>0.0</td>\n",
       "      <td>54</td>\n",
       "      <td>False</td>\n",
       "      <td>POINT (85.27137 23.51662)</td>\n",
       "      <td>8.897919</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>346684</th>\n",
       "      <td>56</td>\n",
       "      <td>google</td>\n",
       "      <td>0.7639</td>\n",
       "      <td>18.1215</td>\n",
       "      <td>4176244229456003072</td>\n",
       "      <td>IND</td>\n",
       "      <td>POINT (85.26911 23.51335)</td>\n",
       "      <td>tuhtu27p</td>\n",
       "      <td>11-20-364-02684-373482</td>\n",
       "      <td>0.0</td>\n",
       "      <td>10</td>\n",
       "      <td>False</td>\n",
       "      <td>POINT (85.26904 23.5135)</td>\n",
       "      <td>18.728227</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>346934</th>\n",
       "      <td>56</td>\n",
       "      <td>google</td>\n",
       "      <td>0.8076</td>\n",
       "      <td>44.6793</td>\n",
       "      <td>4176244229456003072</td>\n",
       "      <td>IND</td>\n",
       "      <td>POINT (85.27203 23.51619)</td>\n",
       "      <td>tuhtu2vp</td>\n",
       "      <td>11-20-364-02684-373482</td>\n",
       "      <td>0.0</td>\n",
       "      <td>53</td>\n",
       "      <td>False</td>\n",
       "      <td>POINT (85.2721 23.51617)</td>\n",
       "      <td>6.582941</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>346857</th>\n",
       "      <td>56</td>\n",
       "      <td>google</td>\n",
       "      <td>0.8746</td>\n",
       "      <td>61.9966</td>\n",
       "      <td>4176244229456003072</td>\n",
       "      <td>IND</td>\n",
       "      <td>POINT (85.27058 23.51617)</td>\n",
       "      <td>tuhtu2up</td>\n",
       "      <td>11-20-364-02684-373482</td>\n",
       "      <td>0.0</td>\n",
       "      <td>59</td>\n",
       "      <td>False</td>\n",
       "      <td>POINT (85.27071 23.51613)</td>\n",
       "      <td>14.341226</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>342876</th>\n",
       "      <td>56</td>\n",
       "      <td>google</td>\n",
       "      <td>0.6514</td>\n",
       "      <td>63.0674</td>\n",
       "      <td>4176244229456003072</td>\n",
       "      <td>IND</td>\n",
       "      <td>POINT (85.00041 23.33629)</td>\n",
       "      <td>tuhky0pr</td>\n",
       "      <td>11-20-364-02695-374346</td>\n",
       "      <td>0.0</td>\n",
       "      <td>20</td>\n",
       "      <td>False</td>\n",
       "      <td>POINT (85.00047 23.33629)</td>\n",
       "      <td>6.287872</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>342179</th>\n",
       "      <td>56</td>\n",
       "      <td>microsoft</td>\n",
       "      <td>NaN</td>\n",
       "      <td>24.6467</td>\n",
       "      <td>4176244229456003072</td>\n",
       "      <td>IND</td>\n",
       "      <td>POINT (84.99538 23.33341)</td>\n",
       "      <td>tuhkwpew</td>\n",
       "      <td>11-20-364-02695-374346</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16</td>\n",
       "      <td>False</td>\n",
       "      <td>POINT (84.99538 23.33335)</td>\n",
       "      <td>6.574817</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        boundary_id  bf_source  confidence  area_in_meters  \\\n",
       "347252           56     google      0.8294         32.1143   \n",
       "346684           56     google      0.7639         18.1215   \n",
       "346934           56     google      0.8076         44.6793   \n",
       "346857           56     google      0.8746         61.9966   \n",
       "342876           56     google      0.6514         63.0674   \n",
       "342179           56  microsoft         NaN         24.6467   \n",
       "\n",
       "                      s2_id country_iso                   geometry   geohash  \\\n",
       "347252  4176244229456003072         IND    POINT (85.2714 23.5167)  tuhtu3hd   \n",
       "346684  4176244229456003072         IND  POINT (85.26911 23.51335)  tuhtu27p   \n",
       "346934  4176244229456003072         IND  POINT (85.27203 23.51619)  tuhtu2vp   \n",
       "346857  4176244229456003072         IND  POINT (85.27058 23.51617)  tuhtu2up   \n",
       "342876  4176244229456003072         IND  POINT (85.00041 23.33629)  tuhky0pr   \n",
       "342179  4176244229456003072         IND  POINT (84.99538 23.33341)  tuhkwpew   \n",
       "\n",
       "                        shrid2  urban  neighbors  isolated  \\\n",
       "347252  11-20-364-02684-373482    0.0         54     False   \n",
       "346684  11-20-364-02684-373482    0.0         10     False   \n",
       "346934  11-20-364-02684-373482    0.0         53     False   \n",
       "346857  11-20-364-02684-373482    0.0         59     False   \n",
       "342876  11-20-364-02695-374346    0.0         20     False   \n",
       "342179  11-20-364-02695-374346    0.0         16     False   \n",
       "\n",
       "            nearest_point_on_road  distance_to_road_m  on_road  \n",
       "347252  POINT (85.27137 23.51662)            8.897919        1  \n",
       "346684   POINT (85.26904 23.5135)           18.728227        1  \n",
       "346934   POINT (85.2721 23.51617)            6.582941        1  \n",
       "346857  POINT (85.27071 23.51613)           14.341226        1  \n",
       "342876  POINT (85.00047 23.33629)            6.287872        1  \n",
       "342179  POINT (84.99538 23.33335)            6.574817        1  "
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sampled_points_on_road_rural.head(10)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": ".venv",
   "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.11.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}