This project was a fun way to practice some of my geospatial Python skills. I chose to use data from my home city of Cambridge, Massachusetts available through the fantastic Cambridge GIS website.
I used the city’s base zoning districts and grouped them into the broad categories provided here. I also obtained point data representing street trees. I then used the OSMnx and GeoPandas Python libraries to calculate the total road length in each district using OpenStreetMap data. Finally, I calculated the number of street trees per meter of road for each zone category and extracted the most common tree species.
This was a fun way to practice working with vector layers and creating geospatial graphics using only Python. I took inspiration for this project from the great Python Foundation series by Spatial Thoughts and Luigi Carboni’s project Retrieving and Visualizing Data from OpenStreetMap.
