In 2013, my family's home burned down in the Black Forest Fire. As the name implies, the Black Forest area is densely wooded, and the fire quickly burned over 14,000 acres and more than 500 homes. Related to the mapping in Project #2, I wanted to create a map of where the fire had burned.
The official fire perimeter map shows the greatest extent of the fire. However, there are many pockets of healthy forest inside the fire perimeter. I needed another source of data that could show the burned areas more directly.
I looked to use remote sensing data as the record covers the time periods of change, and the spatial resolution would be high enough. For processing, I used Google Earth Engine. It's a powerful tool that aggregates remote sensing datasets, and also provides a framework to do large scale data processing in the cloud. Unfortunately, the python API is quite challenging for me to use in my typical mindset of local data processing, and it seems like a second-class citizen to their JavaScript API. So, I used the JavaScript API through the Earth Engine Code Editor for this project.
I searched for an appropriate remote sensing dataset that might capture the change in vegetation. Initially, I anticipated using a normalized difference vegetation index (NDVI) which uses the ratio of infrared to red light to detect the presence and amount of vegetation. Here, you can see a visualization of NDVI in Black Forest throughout 2013. The dark green area in the center is Black Forest, and you can see a lighter scar emerge following the fire.
NDVI could be used for this analysis, but I found the Global Forest Change dataset from Hansen et al. This uses data from Landsat 8, and marks changes in forest density from a reference baseline in the year 2000. The data clearly shows regions of forest loss, and forest gain. It's astonishing to be able to instantly get a snapshot of forest change at a 30m resolution across the globe. After limiting the data to the Black Forest area, the burn scar becomes apparent.The burn scar can also be seen in the Landsat Burned Area Index dataset from August of 2013. White areas represent a burned classification.
Next, I wanted to convert this raster pixel data into a vector format for display on a stylized map. In order to do this, I needed to convert the pixel layer to a vector form in Google Earth Engine. The vector data still had retained the pixel artifact structure, so I exported the GeoJSON data to QGIS, and then filtered the small regions, smoothed the edges using a "snakes" generalization algorithm, and removed any islands. This smoothing and filtering took much longer than the remote sensing part! Below, you can see the vector output before and after filtering.
I was happy with how the map of the fire turned out. The sensed burned areas correspond to what I see in the forest, and the vector format output looks good on the map. I added the burn scar to a map of the Pineries Open Space that I've been working on.
My biggest takeaway from this project is the ease of access for this powerful remote sensing data. After getting spun up on the tools, there is no shortage of data available to analyze, and this project only scratches the surface. It's astonishing to be able to work with data with this level of spatial fidelity, temporal span and resolution, and multimodal sensing.