Working with NASA EMIT data in HyperCoast¶

This notebook demonstrates how to work with NASA Earth Surface Mineral Dust Source Investigation (EMIT) data in HyperCoast.

Environment setup¶

Uncomment and run the following cell to install the required packages.

In [1]:

# %pip install "hypercoast[extra]"

# %pip install "hypercoast[extra]"

Import library.

In [2]:

import hypercoast

import hypercoast

To download and access the data, you will need to create an Earthdata login. You can register for an account at urs.earthdata.nasa.gov. Once you have an account, run the following cell and enter your NASA Earthdata login credentials.

In [3]:

hypercoast.nasa_earth_login()

hypercoast.nasa_earth_login()

Search for EMIT data¶

Search for EMIT data programmatically. Specify the bounding box and time range of interest. Set count=-1 to return all results or set count=10 to return the first 10 results.

In [4]:

results, gdf = hypercoast.search_emit(
    bbox=(-83, 25, -81, 28),
    temporal=("2024-04-01", "2024-05-16"),
    count=10,  # use -1 to return all datasets
    return_gdf=True,
)

results, gdf = hypercoast.search_emit( bbox=(-83, 25, -81, 28), temporal=("2024-04-01", "2024-05-16"), count=10, # use -1 to return all datasets return_gdf=True, )

Plot the footprints of the returned datasets on a map.

In [5]:

gdf.explore()

gdf.explore()

Out[5]:

Make this Notebook Trusted to load map: File -> Trust Notebook

Uncomment the following cell to download the first dataset from the search results. Note that the download may take some time.

In [6]:

# hypercoast.download_emit(results[:1], out_dir="data")

# hypercoast.download_emit(results[:1], out_dir="data")

Search for EMIT data interactively. Specify pan and zoom to the area of interest. Specify the time range of interest from the search dialog, then click on the Search button.

In [7]:

m = hypercoast.Map(center=[30.0262, -90.1345], zoom=8)
m.search_emit()
m

m = hypercoast.Map(center=[30.0262, -90.1345], zoom=8) m.search_emit() m

Out[7]:

Uncomment the following cell to display the GeoDataFrame of the search results.

In [8]:

# m._NASA_DATA_GDF.head()

# m._NASA_DATA_GDF.head()

Similarly, you can download the first dataset from the search results by uncommenting the following cell.

In [9]:

# hypercoast.download_emit(results[:1], out_dir="data")

# hypercoast.download_emit(results[:1], out_dir="data")

Download a sample EMIT dataset¶

Let's download a sample EMIT dataset for the demonstration.

In [10]:

url = "https://github.com/opengeos/datasets/releases/download/hypercoast/EMIT_L2A_RFL_001_20240404T161230_2409511_009.nc"
filepath = "../examples/data/EMIT_L2A_RFL_001_20240404T161230_2409511_009.nc"
hypercoast.download_file(url, filepath, quiet=True)

url = "https://github.com/opengeos/datasets/releases/download/hypercoast/EMIT_L2A_RFL_001_20240404T161230_2409511_009.nc" filepath = "../examples/data/EMIT_L2A_RFL_001_20240404T161230_2409511_009.nc" hypercoast.download_file(url, filepath, quiet=True)

../examples/data/EMIT_L2A_RFL_001_20240404T161230_2409511_009.nc already exists. Skip downloading. Set overwrite=True to overwrite.

Out[10]:

'/home/runner/work/HyperCoast/HyperCoast/docs/examples/data/EMIT_L2A_RFL_001_20240404T161230_2409511_009.nc'

Read EMIT data¶

Read the downloaded EMIT data and process it as an xarray.Dataset. Note that the dataset has 285 bands.

In [11]:

dataset = hypercoast.read_emit(filepath)
# dataset

dataset = hypercoast.read_emit(filepath) # dataset

Visualize EMIT data¶

Visualize the EMIT data on an interactive map. You can change the band combination and extract spectral profiles interactively. You can also export the spectral profiles as a CSV file.

In [12]:

m = hypercoast.Map()
m.add_basemap("SATELLITE")
m.add_emit(dataset, wavelengths=[1000, 600, 500], vmin=0, vmax=0.3, layer_name="EMIT")
m.add("spectral")
m

m = hypercoast.Map() m.add_basemap("SATELLITE") m.add_emit(dataset, wavelengths=[1000, 600, 500], vmin=0, vmax=0.3, layer_name="EMIT") m.add("spectral") m

Out[12]:

Create an image cube¶

First, select a subset of the data to avoid nodata areas.

In [13]:

ds = dataset.sel(longitude=slice(-90.1482, -89.7321), latitude=slice(30.0225, 29.7451))

ds = dataset.sel(longitude=slice(-90.1482, -89.7321), latitude=slice(30.0225, 29.7451))

Visualize the EMIT data in 3D with an RGB image overlaid on top of the 3D plot.

In [14]:

p = hypercoast.image_cube(
    ds,
    variable="reflectance",
    cmap="jet",
    clim=(0, 0.4),
    rgb_wavelengths=[1000, 700, 500],
    rgb_gamma=2,
    title="EMIT Reflectance",
)

p = hypercoast.image_cube( ds, variable="reflectance", cmap="jet", clim=(0, 0.4), rgb_wavelengths=[1000, 700, 500], rgb_gamma=2, title="EMIT Reflectance", )

Uncomment the following cell to create an image cube. Note that this function does not work in the Google Colab environment.

In [15]:

# p.show()

# p.show()

Interactive slicing¶

First, select a subset of the data for demonstration purposes.

In [16]:

ds = dataset.sel(longitude=slice(-90.05, -89.99), latitude=slice(30.00, 29.93))

ds = dataset.sel(longitude=slice(-90.05, -89.99), latitude=slice(30.00, 29.93))

Drag the plane up and down to slice the data in 3D.

In [17]:

p = hypercoast.image_cube(
    ds,
    variable="reflectance",
    cmap="jet",
    clim=(0, 0.5),
    rgb_wavelengths=[1000, 700, 500],
    rgb_gamma=2,
    title="EMIT Reflectance",
    widget="plane",
)
p.add_text("Band slicing", position="upper_right", font_size=14)

p = hypercoast.image_cube( ds, variable="reflectance", cmap="jet", clim=(0, 0.5), rgb_wavelengths=[1000, 700, 500], rgb_gamma=2, title="EMIT Reflectance", widget="plane", ) p.add_text("Band slicing", position="upper_right", font_size=14)

Out[17]:

<CornerAnnotation(0x5627f3116910) at 0x7fb9f4704a60>

Uncomment the following cell to display the interactive slicing widget. Note that this function does not work in the Google Colab environment.

In [18]:

# p.show()

# p.show()

Interactive thresholding¶

Drag the threshold slider to threshold the data in 3D.

In [19]:

p = hypercoast.image_cube(
    ds,
    variable="reflectance",
    cmap="jet",
    clim=(0, 0.5),
    rgb_wavelengths=[1000, 700, 500],
    rgb_gamma=2,
    title="EMIT Reflectance",
    widget="threshold",
)
p.add_text("Thresholding", position="upper_right", font_size=14)

p = hypercoast.image_cube( ds, variable="reflectance", cmap="jet", clim=(0, 0.5), rgb_wavelengths=[1000, 700, 500], rgb_gamma=2, title="EMIT Reflectance", widget="threshold", ) p.add_text("Thresholding", position="upper_right", font_size=14)

Out[19]:

<CornerAnnotation(0x5627f2870bb0) at 0x7fb9f4467220>

Uncomment the following cell to display the thresholded data. Note that this function does not work in the Google Colab environment.

In [20]:

# p.show()

# p.show()