Data Science, GIS, Risk Management, Insuretech.
numpy as np from rasterio.features import rasterize import geopandas as gpd features_buffered = gpd.GeoDataFrame({'geometry': buffered_geometries}) r_total = np.zeros((w,h)) for i in range(features_wgs.shape[0]): print(i/features_wgs.shape[0]) geometry = features_buffered['geometry'].values[i]