Visualizing PRISMA L2D data interactively with prismatools¶

This notebook demonstrates how to visualize PRISMA hyperspectral L2D data interactively with prismatools. For more information about PRISMA, please visit the links below:

• https://www.asi.it/en/earth-science/prisma/
• https://en.wikipedia.org/wiki/PRISMA_(spacecraft)

Note:

• Register and obtain the PRISMA datasets from the ASI PRISMA portal
• Once you have downloaded a dataset, replace the filepath which we will use in this notebook.

In [1]:

# %pip install prismatools

# %pip install prismatools

In [2]:

import prismatools

import prismatools

In [3]:

file = "/path/to/prisma.he5"  # example: r"..\..\data\PRS_L2D_STD_20240429095823_20240429095827_0001.he5"

file = "/path/to/prisma.he5" # example: r"..\..\data\PRS_L2D_STD_20240429095823_20240429095827_0001.he5"

Read the PRISMA data as an xarray.Dataset object.

In [4]:

# ds = prismatools.read_prismaL2D(file_path=file)
# ds

# ds = prismatools.read_prismaL2D(file_path=file) # ds

Create an interactive map.

In [5]:

# m = prismatools.Map()
# m

# m = prismatools.Map() # m

Add the PRISMA data to the map.

In [6]:

# m.add_prisma(ds, wavelengths=[650.0, 550.0, 450.0], vmin=0, vmax=0.2)
# m.add("spectral")
# m

# m.add_prisma(ds, wavelengths=[650.0, 550.0, 450.0], vmin=0, vmax=0.2) # m.add("spectral") # m

In [7]:

# m.add_prisma(ds, wavelengths=[650.0, 550.0, 450.0], vmin=0, vmax=0.2)
# m.add('spectral')

# m.add_prisma(ds, wavelengths=[650.0, 550.0, 450.0], vmin=0, vmax=0.2) # m.add('spectral')