calculate_indices function

Calculate indices from the bands of a raster