# calculates pre-fire NBR
pre_fire_NBR<-(reflectance.jan08NIR-reflectance.jan08SWIR2)/reflectance.jan08NIR+reflectance.jan08SWIR2)
  
# calculates post-fire NBR
post_fire_NBR<-(reflectance.feb25NIR-reflectance.feb25SWIR2)/(reflectance.feb25NIR+reflectance.feb25SWIR2)