####### Estimating the Area of Every burnt class in Hectares or other units when required ########
v<-extract(dNBR_reclass, Empedrado_shp) #Extract raster values to polygons 
v.counts <- lapply(v,table) # Get class counts for each polygon depending on the number of polygons
class.df <- as.data.frame(t(sapply(v.counts,'[',1:length(unique(dNBR_reclass)))))  # Create a data.frame where missing classes are NA
Empedrado@data<- data.frame(Empedrado@data, class.df)# Add back to polygon data
Empedrado_Attributes<-monchique@data
Empedrado_Burn_Severity<-subset(Empedrado_Attributes, select = c(NAME_0, NAME_2, X.1, X1, X2, X3, X4, X5, X6, X7))
##### The size of a landsat pixel is 30m X30m, and the area is 900sqm= 0.9 hectare ####
Empedrado_Burn_Severity$"No Data"<-Empedrado_Burn_Severity$ X.1*0.09 # multiply to get ha
Empedrado_Burn_Severity$"High Enhanced Regrowth hectares"<-Empedrado_Burn_Severity$X1*0.09
Empedrado_Burn_Severity$"Low Enhanced Regrowth hectares"<-Empedrado_Burn_Severity$X2*0.09
Empedrado_Burn_Severity$"Unburned hectares"<-Empedrado_Burn_Severity$X3*0.09
Empedrado_Burn_Severity$"Low Severity hectares"<-Empedrado_Burn_Severity$X4*0.09
Empedrado_Burn_Severity$"Moderate-Low Severity hectares"<-Empedrado_Burn_Severity$X5*0.09
Empedrado_Burn_Severity$"Moderate-High Severity hectares"<-Empedrado_Burn_Severity$X6*0.09
Empedrado_Burn_Severity$"High Severity hectares"<-Empedrado_Burn_Severity$X7*0.09

### Removing unwanted columns and saving information both as csv and as shapefile ###
Empedrado_Burn_Severity<- subset(Empedrado_Burn_Severity, select = -c(X.1, X1, X2, X3, X4, X5, X6, X7)) # To remove unwanted columns 
write.csv(Empedrado_Burn_Severity,"G:/landsat 08/Empedrado_Burn_Severity_21_Aug_2018.csv") # write the csv with values                      
monchiqueBurntHectares<-merge(Empedrado,Empedrado_Burn_Severity,by=("NAME_2")) # merge csv to shapefile
writeOGR(EmpedradoBurntHectares,"EmpedradoBurntHectares_21_Aug_2018.shp", "EmpedradoBurntHectares_21_Aug_2018", driver = "ESRI Shapefile" ) # write shapefile with burnt classes estimated in hectares
plot(Empedrado_Burn_Severity)