Burn severity in R

In this project, I wanted to calculate the burn severity after the fire in Ljusdal, Sweden in 2018 (the exact duration was: 14.07.2018. - 09.08.2018.) by using RStudio.  First, I needed two Landsat-8 images (from USGS), one before (06.26.2018.) and one after the fire (15.08.2018.) to assess the changes caused by the fire.  I calculated the dNBR (difference normalized burn ratio) for this, then assigned the standard ranges for burn severity and created the first map. 

I wanted to show the land cover as well, to see which types were mainly affected by the fire. I used the WorldCover dataset (second map) and calculated the distribution of land cover types in only the fire-affected area (it was defined in another research). The third map shows the result of this. The maps were created in RStudio.