# !pip install geemap

import ee
import geemap

Map = geemap.Map()
Map

URL = 'https://storage.googleapis.com/pdd-stac/disasters/hurricane-harvey/0831/20170831_172754_101c_3B_AnalyticMS.tif'
# URL = 'gs://pdd-stac/disasters/hurricane-harvey/0831/20170831_172754_101c_3B_AnalyticMS.tif'

image = geemap.load_GeoTIFF(URL)

vis = {"bands": ["B3", "B2", "B1"], "min": 3000, "max": 13500}

Map.addLayer(image, vis, 'Cloud Image')
Map.centerObject(image.geometry(), 12)

Map = geemap.Map()
Map

B3 = 'gs://gcp-public-data-landsat/LC08/01/044/034/LC08_L1TP_044034_20131228_20170307_01_T1/LC08_L1TP_044034_20131228_20170307_01_T1_B3.TIF'
B4 = 'gs://gcp-public-data-landsat/LC08/01/044/034/LC08_L1TP_044034_20131228_20170307_01_T1/LC08_L1TP_044034_20131228_20170307_01_T1_B4.TIF'
B5 = 'gs://gcp-public-data-landsat/LC08/01/044/034/LC08_L1TP_044034_20131228_20170307_01_T1/LC08_L1TP_044034_20131228_20170307_01_T1_B5.TIF'

URLs = [B3, B4, B5]

collection = geemap.load_GeoTIFFs(URLs)

image = collection.toBands().rename(['Green', 'Red', 'NIR']).selfMask()

vis = {'bands': ['NIR', 'Red', 'Green'], 'min': 100, 'max': 12000, 'gamma': 0.8}

Map.addLayer(image, vis, 'Image')
Map.centerObject(image.geometry(), 8)

ndvi = image.normalizedDifference(['NIR', 'Red'])
ndvi_vis = {'min': 0, 'max': 1, 'palette': ['blue', 'green']}
Map.addLayer(ndvi, {}, 'NDVI')