Mapping the Border Markers

import requests

# URL of the ArcGIS Feature Server query
url = "https://gisserver.ibwc.gov/agsserver/rest/services/USIBWCPublic/US_MexicoIBLandMonuments/FeatureServer/0/query?f=geojson&resultOffset=0&resultRecordCount=300&where=1%3D1&orderByFields=&outFields=*&returnGeometry=true&spatialRel=esriSpatialRelIntersects"

# Send a GET request to the URL
response = requests.get(url)

# Check if the request was successful
if response.status_code == 200:
    # Save the data to a file
    with open("monuments_ibwc.geojson", "w") as file:
        file.write(response.text)
    print("Data downloaded successfully.")
else:
    print(f"Failed to download data. Status code: {response.status_code}")

import geopandas as gpd

# Load the GeoJSON file
gdf = gpd.read_file('monuments_ibwc.geojson')

# Display the first few rows of the GeoDataFrame
print(gdf.head())

import geopandas as gpd
import folium
from folium import GeoJson, GeoJsonPopup, GeoJsonTooltip

# Load the GeoJSON file
gdf = gpd.read_file('monuments_ibwc.geojson')

# Create a base map centered around the GeoJSON data
m = folium.Map(location=[gdf.geometry.centroid.y.mean(), gdf.geometry.centroid.x.mean()], zoom_start=7)

# Define the popup and tooltip
popup = GeoJsonPopup(
    fields=["Monument"],
    aliases=["Monument:"],
    localize=True,
    labels=True,
    style="background-color: yellow;",
)

# Define the style function
def style_function(feature):
    return {
        "fillColor": "#000000",
        "color": "black",
        "weight": 2,
        "fillOpacity": 0.5,
    }

# Add the GeoJSON data to the map with the defined popup and tooltip
GeoJson(
    gdf,
    style_function=style_function,
    popup=popup,
).add_to(m)

# Save the map to an HTML file
m.save('map.html')

# Optionally, display the map in a Jupyter notebook
m

## CONVERT TO KML
from pykml.factory import KML_ElementMaker as KML
from pykml.factory import GX_ElementMaker as GX
from lxml import etree

# Read the GeoJSON file
geojson_file = 'monuments_ibwc.geojson'
gdf = gpd.read_file(geojson_file)

# Create a KML document
kml_doc = KML.kml(
    KML.Document(
        *[KML.Placemark(
            KML.name(gdf.loc[i, 'Monument']),
            KML.Point(
                KML.coordinates(f"{gdf.geometry[i].x},{gdf.geometry[i].y}")
            )
        ) for i in range(len(gdf))]
    )
)

# Write the KML document to a file
kml_file = 'monuments_ibwc.kml'
with open(kml_file, 'wb') as f:
    f.write(etree.tostring(kml_doc, pretty_print=True, xml_declaration=True, encoding='UTF-8'))

print(f'GeoJSON file has been converted to KML file: {kml_file}')