This is the current code:
import geopandas as gpd
from shapely.geometry import Polygon
lat_point_list = gdf['latitude']
lon_point_list = gdf['longitude']
polygon_geom = Polygon(zip(lon_point_list, lat_point_list))
crs = {'init': 'epsg:4326'}
polygon = gpd.GeoDataFrame(index=[0], crs=crs, geometry=[polygon_geom])
print(polygon.geometry)
polygon.to_file(filename='polygon.geojson', driver='GeoJSON')
polygon.to_file(filename='polygon.shp', driver="ESRI Shapefile")
m = folium.Map([33.435598, -112.349602], zoom_start=5, tiles='cartodbpositron')
folium.GeoJson(polygon).add_to(m)
folium.LatLngPopup().add_to(m)
m
The map im getting is all over out of bounds:
Could this be because of latitude, longitude mixed places? Thank you. In essence, the goal is to create a Choropleth map of US crime stats of cities(their bounds of coordinates.