For geospatial analysis, I used a heatmap plots from geopandas point geodataframes. Defined the heatmap function which takes a GeoDataFrame with point geometries and shows a matplotlib plot of heatmap density. This is done using numpy's 2D histogram binning with smoothing from scipy.7