add fake cellsites cp

README.md

@@ -1,2 +1,4 @@
 # sds4gdsp
 intro to spatial data science for the globe data science program
+
+1. Setup

scripts/make_fake_cellsites.py

@@ -0,0 +1,24 @@
+# Download Instructions
+# 1. Goto https://gadm.org/download_country.html#google_vignette
+# 2. Select Philippines then download GeoJSON and select level2
+# 3. Put it in your root folder `data/shapefiles/<file>`
+
+# REMOVE CODE BELOW TEMP ONLY
+# import os
+# os.chdir("../")
+
+import geopandas as gpd
+
+filepath = "data/shapefile/gadm41_PHL_2.json.zip"
+gadm = gpd.read_file(filepath)
+
+# filter the town of interest
+# let's use Taguig City for now
+# change when you want to explore
+town_keyword = "Taguig"
+town_filter = gadm.NAME_2==town_keyword
+gadm = gadm.loc[town_filter]
+
+# save file to local disk
+filepath = f"data/shapefile/gadm_{town_keyword.lower()}.gpkg"
+gadm.reset_index(drop=True).to_file(filepath, driver="GPKG")