Layered Rasters with custom colors

[markobud]

Hi!
I'm trying to overlay different rasters, each with a different color scale, but I'm unable to make it work (below my current approach). Do you have any advice on how to move forward?

Thanks,
Marko

import os
import leafmap.leafmap as leafmap

landsat = 'landsat.tif'

url1 = "https://open.gishub.org/data/raster/landsat7.tif"
satellite = leafmap.download_file(url1, "landsat.tif")

m = leafmap.Map()

for b,name,cmap_name in zip([2,1,3],['red','green','blue'],['Reds','Greens','Blues']):
    m.add_raster(landsat, bands=b, layer_name=f'Landsat_{name}',cmap = cmap_name)

m

[giswqs]

Thank you for reporting. This bug has been fixed in #674. You can update the package using leafmap.update_package() after the PR is merged.

[markobud]

Thanks a lot! Your solution works perfectly now. However, it does not seem to process correctly custom colormaps:

import os
import leafmap.leafmap as leafmap
from matplotlib.colors import ListedColormap

landsat = 'landsat.tif'

url1 = "https://open.gishub.org/data/raster/landsat7.tif"
satellite = leafmap.download_file(url1, "landsat.tif")

cmap1 = ListedColormap(['red','blue'],name='test_1')
cmap2 = ListedColormap(['blue','green'],name='test_2')
cmap3 = ListedColormap(['green','red'],name='test_3')

m = leafmap.Map()

for b,name,cmap in zip([2,1,3],['test_1','test_2','test_3'],[cmap1,cmap2,cmap3]):
    m.add_raster(landsat, bands=b, layer_name=f'Landsat_{name}',cmap = cmap)

m

[giswqs]

@markobud This seems to be an issue with localtileserver or rio-tiler. The following example does not use leafmap. It can't render the image, either. rio-tiler appears to support custom colormaps. See https://cogeotiff.github.io/rio-tiler/colormap/#default-rio-tilers-colormaps. I have not figured out how to use custom colormaps with localtileserver. Please report back if you figure it out.

cc @banesullivan

from localtileserver import TileClient, get_leaflet_tile_layer
from matplotlib.colors import ListedColormap
from ipyleaflet import Map

client = TileClient('https://github.com/giswqs/data/raw/main/raster/landsat7.tif')
colormap = ListedColormap(['red','blue'],name='test_1')
t = get_leaflet_tile_layer(client, indexes=1, colormap=colormap)
m = Map(center=client.center(), zoom=client.default_zoom)
m.add(t)
m

[banesullivan]

Ah yes, this is the one thing that I was unable to handle in localtileserver's transition from large-image to rio-tiler. Its possible to do this, I just need to account for in localtileserver.

Sorry for the loss of functionality! I didn't think many people were using custom colormaps like this

[markobud]

Thanks both for your time! @giswqs and @banesullivan!