def calculate_ndbi(swir, nir):

    ndbi = (swir - nir) / (swir + nir + 1e-10)

    return ndbi