ICESat-2 and Sentinel 2 cloud access
ICESat-2 and Sentinel 2 cloud access#
This notebook (
download) builds off of the icepyx IS2_cloud_data_access.ipynb and ICESat-2 Hackweek Data Integration 1 tutorials. It illustrates the use of icepyx for accessing ICESat-2 data currently available through the AWS (Amazon Web Services) us-west2 hub s3 data bucket as well as data integration with Sentinel 2 (cloud-optimized geotiff) and ATM (downloaded csv) datasets.
Identify and locate ICESat-2 and Sentinel 2 data
Acquire data from the cloud
Open data in
xarrayand basic functioning of DataFrames
By the end of this tutorial, you will be able to visualize Sentinel 2 Cloud Optimized Geotiffs with ICESat-2 and ATM data.
ICESat-2 data became publicly available on the cloud on 29 September 2022. Thus, access methods and example workflows are still being developed by NSIDC, and the underlying code in icepyx will need to be updated now that these data (and the associated metadata) are available. We appreciate your patience and contributions (e.g. reporting bugs, sharing your code, etc.) during this transition!
This example and the code it describes are part of ongoing development. Current limitations to using these features are described throughout the example, as appropriate.
You MUST be working within an AWS instance. Otherwise, you will get a permissions error.
Authentication is still more steps than we’d like. We’re working to address this - let us know if you’d like to join the conversation!
For this tutorial, feel free to run the code along with us as we live code by downsizing the zoom window and splitting your screen (or using two screens). Or you can simply watch the zoom walkthrough. Don’t worry if you fall behind on the code. The notebook is standalone and you can easily run the code at your own pace another time to catch anything you missed.
We’ll be using the following open source Python libraries in this notebook:
# Suppress library deprecation warnings import logging logging.captureWarnings(True)
import ipyleaflet from ipyleaflet import Map, basemaps, basemap_to_tiles, Polyline import ipywidgets import datetime import re
%matplotlib widget import pystac_client import geopandas as gpd import h5py import ast import pandas as pd import geoviews as gv import hvplot.pandas from ipywidgets import interact from IPython.display import display, Image import intake # if you've installed intake-STAC, it will automatically import alongside intake import intake_stac import xarray as xr import matplotlib.pyplot as plt import boto3 import rasterio as rio from rasterio.session import AWSSession from rasterio.plot import show import rioxarray as rxr from dask.utils import SerializableLock import os import fiona import hvplot.xarray import numpy as np from pyproj import Proj, transform