New surface temperature colormap

[anayeaye]

What

Add a new custom linear colormap for surface temperature datasets. The example below uses specific values for transition stops in the mapping--ultimately what is needed is a continuous linear transition between the color stops.

Matplotlib

#Create colormap:
def interpolate_color(color1, color2, fraction):
    r1, g1, b1 = mcolors.hex2color(color1)
    r2, g2, b2 = mcolors.hex2color(color2)
    r = r1 + fraction * (r2 - r1)
    g = g1 + fraction * (g2 - g1)
    b = b1 + fraction * (b2 - b1)
    return (r, g, b)
def temp_colormap():
    temp_colormap = []
    stops = [
        "#3CCBCE",   # -30
        "#C5F8FF",   # -20
        "#FEC5FF",   # -10
        "#E079FB",   # 0
        "#094FC9",   # 10
        "#009FFF",   # 20
        "#44E2FF",   # 30
        "#147F4F",   # 40
        "#79B32C",   # 50
        "#FDFE00",   # 60
        "#FF8700",   # 70
        "#FF0F00",   # 80
        "#9D0F2B",   # 90
        "#4D0000",   # 100
        "#BA2E6D",   # 110
    ]
    for i in range(len(stops) - 1):
        color1 = stops[i]
        color2 = stops[i + 1]
        for j in range(10):
            fraction = j / 10.0
            interpolated_color = interpolate_color(color1, color2, fraction)
            temp_colormap.append(interpolated_color)
    return mcolors.ListedColormap(temp_colormap)

Sample visualization

AC

• named custom surface temperature colormap can be used in tile xyz requests

[smohiudd]

@acblackford what collection is this colormap associated with? Do you have a sample file I can test out?

[acblackford]

@smohiudd This colormap will be associated with MERRA2 and ERA5 2 meter temperature reanalysis data for five select blizzard case dates. I do not have them ingested into STAC yet (will be done soon) but they are in the s3 bucket veda-data-store-staging/Blizzards/ with the general name structures ERA5_2m_Temp_cog_YYYY-MM-DDTHH:MM:SSZ.tif and MERRA2_2m_Temp_cog_YYYY-MM-DDTHH:MM:SSZ.tif .

[acblackford]

@smohiudd Here are the STAC collections that this colormap will be used on:

https://staging.openveda.cloud/stac/collections/blizzard-merra2-2m-temp

https://staging.openveda.cloud/stac/collections/blizzard-era5-2m-temp