def calculate_ndwi(green, nir):

    ndwi = (green - nir) / (green + nir + 1e-10)

    return ndwi