code: def ndsi(self):
          
          green, swir1 = self.reflectance(3), self.reflectance(6)
          ndsi = self._divide_zero((green - swir1), (green + swir1), nan)
  
          return ndsi
docs: Normalized difference snow index.
          :return: NDSI