def _get_modis_learning_bands(domain, b):
    diff        = b['b2'].subtract(b['b1'])
    ratio       = b['b2'].divide(b['b1'])
    modisBands  = b['b1'].addBands(b['b2']).addBands(diff).addBands(ratio).addBands(b['NDVI']).addBands(b['NDWI'])
    return modisBands
Set up features for the classifier to be trained on: [b2, b2/b1, b2/b1, NDVI, NDWI]
dataset/ETHPy150Open nasa/CrisisMappingToolkit/cmt/modis/ee_classifiers.py/_get_modis_learning_bands