# Import the following libraries
# For fetching from the Raster API
import requests
# For making maps
import folium
import folium.plugins
from folium import Map, TileLayer
# For talking to the STAC API
from pystac_client import Client
# For working with data
import pandas as pd
# For making time series
import matplotlib.pyplot as plt
# For formatting date/time data
import datetime
# Custom functions for working with GHGC data via the API
import ghgc_utils
Geostationary Satellite Observations of Extreme and Transient Methane Emissions from Oil and Gas Infrastructure
Access this Notebook
You can launch this notebook in the US GHG Center JupyterHub by clicking the link below. If you are a new user, you should first sign up for the hub by filling out this request form and providing the required information.
Access the Geostationary Satellite Observations of Extreme and Transient Methane Emissions from Oil and Gas Infrastructure notebook in the US GHG Center JupyterHub.
Table of Contents
Data Summary and Application
- Spatial coverage: North and South America
- Spatial resolution: 1 km x 1 km (nadir)
- Temporal extent: April 7, 2019 - July 26, 2023
- Temporal resolution: 5 minutes
- Unit: Total column methane enhancement above background, moles of methane per square meter
- Utility: Climate Research
For more, visit the Geostationary Satellite Observations of Extreme and Transient Methane Emissions from Oil and Gas Infrastructure data overview page.
Approach
- Identify available dates and temporal frequency of observations for the given collection using the GHGC API
/stac
endpoint. The collection processed in this notebook is the gridded methane emissions data product. - Pass the STAC item into the raster API
/collections/{collection_id}/items/{item_id}/tilejson.json
endpoint. - Using
folium.plugins.DualMap
, we will visualize two tiles (side-by-side), allowing us to compare time points.
About the Data
Geostationary Satellite Observations of Extreme and Transient Methane Emissions from Oil and Gas Infrastructure
A sample of methane plumes from point sources observed since 2019 by the U.S. Geostationary Operational Environmental Satellites (GOES) over North and South America.
For more information regarding this dataset, please visit the Geostationary Satellite Observations of Extreme and Transient Methane Emissions from Oil and Gas Infrastructure data overview page.
Terminology
Navigating data via the GHGC API, you will encounter terminology that is different from browsing in a typical filesystem. We’ll define some terms here which are used throughout this notebook. - catalog
: All datasets available at the /stac
endpoint - collection
: A specific dataset, e.g. Geostationary Satellite Observations of Extreme and Transient Methane Emissions from Oil and Gas Infrastructure - item
: One granule in the dataset, e.g. one file of GOES observations - asset
: A variable available within the granule, e.g. CO₂ fluxes - STAC API
: SpatioTemporal Asset Catalogs - Endpoint for fetching metadata about available datasets - Raster API
: Endpoint for fetching data itself, for imagery and statistics
Install the Required Libraries
Required libraries are pre-installed on the GHG Center Hub. If you need to run this notebook elsewhere, please install them with this line in a code cell:
%pip install requests folium rasterstats pystac_client pandas matplotlib –quiet
Query the STAC API
STAC API Collection Names
Now, you must fetch the dataset from the STAC API by defining its associated STAC API collection ID as a variable. The collection ID, also known as the collection name, for the Geostationary Satellite Observations of Extreme and Transient Methane Emissions from Oil and Gas Infrastructure dataset is goes-ch4plume-v1.*
**You can find the collection name of any dataset on the GHGC data portal by navigating to the dataset landing page within the data catalog. The collection name is the last portion of the dataset landing page’s URL, and is also listed in the pop-up box after clicking “ACCESS DATA.”*
# Provide STAC and RASTER API endpoints
= "https://earth.gov/ghgcenter/api/stac"
STAC_API_URL = "https://earth.gov/ghgcenter/api/raster"
RASTER_API_URL
# Please use the collection name similar to the one used in STAC collection.
# Name of the collection for gosat budget methane.
= "goes-ch4plume-v1" collection_name
# Fetch the collection from the STAC API using the appropriate endpoint
# The 'pystac_client' library enables an HTTP request
= Client.open(STAC_API_URL)
catalog = catalog.get_collection(collection_name)
collection
# Print the properties of the collection to the console
collection
- type "Collection"
- id "goes-ch4plume-v1"
- stac_version "1.0.0"
- description "A sample of methane plumes from point sources observed since 2019 by the U.S. Geostationary Operational Environmental Satellites (GOES) over North and South America"
links[] 5 items
0
- rel "items"
- href "https://earth.gov/ghgcenter/api/stac/collections/goes-ch4plume-v1/items"
- type "application/geo+json"
1
- rel "parent"
- href "https://earth.gov/ghgcenter/api/stac/"
- type "application/json"
2
- rel "root"
- href "https://earth.gov/ghgcenter/api/stac"
- type "application/json"
- title "US GHG Center STAC API"
3
- rel "self"
- href "https://earth.gov/ghgcenter/api/stac/collections/goes-ch4plume-v1"
- type "application/json"
4
- rel "http://www.opengis.net/def/rel/ogc/1.0/queryables"
- href "https://earth.gov/ghgcenter/api/stac/collections/goes-ch4plume-v1/queryables"
- type "application/schema+json"
- title "Queryables"
stac_extensions[] 2 items
- 0 "https://stac-extensions.github.io/render/v1.0.0/schema.json"
- 1 "https://stac-extensions.github.io/item-assets/v1.0.0/schema.json"
renders
rad
assets[] 1 items
- 0 "rad"
rescale[] 1 items
0[] 2 items
- 0 0.0
- 1 0.4
- colormap_name "plasma"
dashboard
assets[] 1 items
- 0 "rad"
rescale[] 1 items
0[] 2 items
- 0 0.0
- 1 0.4
- colormap_name "plasma"
item_assets
rad
- type "image/tiff; application=geotiff; profile=cloud-optimized"
roles[] 2 items
- 0 "data"
- 1 "layer"
- title "Methane Column Enhancement"
- description "Sample of methane plumes from point sources."
- dashboard:is_periodic True
- dashboard:time_density "daily"
- title "Geostationary Satellite Observations of Extreme and Transient Methane Emissions from Oil and Gas Infrastructure"
extent
spatial
bbox[] 1 items
0[] 4 items
- 0 -104.72692922440127
- 1 25.251622016105113
- 2 -86.82596069124111
- 3 41.12229053684302
temporal
interval[] 1 items
0[] 2 items
- 0 "2019-04-07T16:16:00Z"
- 1 "2023-07-26T18:01:00Z"
- license "CC-BY-4.0"
summaries
datetime[] 2 items
- 0 "2019-04-07T16:16:00Z"
- 1 "2023-07-26T18:01:00Z"
Examining the contents of our collection
under the temporal
variable, we see that the data is available starting April 7, 2019 . By looking at the dashboard:time density
, we observe that the data is available for multiple days.
# The search function lets you search for items within a specific date/time range
= catalog.search(
search =collection_name,
collections=['2019-05-11T00:00:00Z','2019-05-13T00:00:00Z']
datetime
)= search.item_collection()
items # Print how many items we found in our search
print(f"# items found: {len(items)}")
# items found: 52
for item in items:
print(item.properties["datetime"])
2019-05-12T19:41:00+00:00
2019-05-12T19:36:00+00:00
2019-05-12T19:31:00+00:00
2019-05-12T19:26:00+00:00
2019-05-12T19:21:00+00:00
2019-05-12T19:16:00+00:00
2019-05-12T19:11:00+00:00
2019-05-12T19:06:00+00:00
2019-05-12T19:01:00+00:00
2019-05-12T18:56:00+00:00
2019-05-12T18:51:00+00:00
2019-05-12T18:46:00+00:00
2019-05-12T18:41:00+00:00
2019-05-12T18:36:00+00:00
2019-05-12T18:31:00+00:00
2019-05-12T18:26:00+00:00
2019-05-12T18:21:00+00:00
2019-05-12T18:16:00+00:00
2019-05-12T18:11:00+00:00
2019-05-12T18:06:00+00:00
2019-05-12T18:01:00+00:00
2019-05-12T17:56:00+00:00
2019-05-12T17:51:00+00:00
2019-05-12T17:46:00+00:00
2019-05-12T17:41:00+00:00
2019-05-12T17:36:00+00:00
2019-05-12T17:31:00+00:00
2019-05-12T17:26:00+00:00
2019-05-12T17:21:00+00:00
2019-05-12T17:16:00+00:00
2019-05-12T17:11:00+00:00
2019-05-12T17:06:00+00:00
2019-05-12T17:01:00+00:00
2019-05-12T16:56:00+00:00
2019-05-12T16:51:00+00:00
2019-05-12T16:46:00+00:00
2019-05-12T16:41:00+00:00
2019-05-11T16:11:00+00:00
2019-05-11T16:06:00+00:00
2019-05-11T16:01:00+00:00
2019-05-11T15:56:00+00:00
2019-05-11T15:51:00+00:00
2019-05-11T15:46:00+00:00
2019-05-11T15:41:00+00:00
2019-05-11T15:36:00+00:00
2019-05-11T15:31:00+00:00
2019-05-11T15:26:00+00:00
2019-05-11T15:21:00+00:00
2019-05-11T15:16:00+00:00
2019-05-11T15:11:00+00:00
2019-05-11T15:06:00+00:00
2019-05-11T15:01:00+00:00
# Examining the first item in the collection
0] items[
- type "Feature"
- stac_version "1.0.0"
stac_extensions[] 2 items
- 0 "https://stac-extensions.github.io/raster/v1.1.0/schema.json"
- 1 "https://stac-extensions.github.io/projection/v1.1.0/schema.json"
- id "goes-ch4plume-v1-GOES-CH4_Mexico_Durango_BV2_BV2-5_2019-05-12T19:41:00Z"
geometry
- type "Polygon"
coordinates[] 1 items
0[] 5 items
0[] 2 items
- 0 -104.40938696034158
- 1 26.06828954903695
1[] 2 items
- 0 -104.11184817858141
- 1 26.06828954903695
2[] 2 items
- 0 -104.11184817858141
- 1 26.266648736877066
3[] 2 items
- 0 -104.40938696034158
- 1 26.266648736877066
4[] 2 items
- 0 -104.40938696034158
- 1 26.06828954903695
bbox[] 4 items
- 0 -104.40938696034158
- 1 26.06828954903695
- 2 -104.11184817858141
- 3 26.266648736877066
properties
- datetime "2019-05-12T19:41:00Z"
links[] 5 items
0
- rel "collection"
- href "https://earth.gov/ghgcenter/api/stac/collections/goes-ch4plume-v1"
- type "application/json"
1
- rel "parent"
- href "https://earth.gov/ghgcenter/api/stac/collections/goes-ch4plume-v1"
- type "application/json"
2
- rel "root"
- href "https://earth.gov/ghgcenter/api/stac"
- type "application/json"
- title "US GHG Center STAC API"
3
- rel "self"
- href "https://earth.gov/ghgcenter/api/stac/collections/goes-ch4plume-v1/items/goes-ch4plume-v1-GOES-CH4_Mexico_Durango_BV2_BV2-5_2019-05-12T19:41:00Z"
- type "application/geo+json"
4
- rel "preview"
- href "https://earth.gov/ghgcenter/api/raster/collections/goes-ch4plume-v1/items/goes-ch4plume-v1-GOES-CH4_Mexico_Durango_BV2_BV2-5_2019-05-12T19:41:00Z/map?assets=rad&rescale=0.0%2C0.4&colormap_name=plasma"
- type "text/html"
- title "Map of Item"
assets
rad
- href "s3://ghgc-data-store/goes-ch4plume-v1/GOES-CH4_Mexico_Durango_BV2_BV2-5_2019-05-12T19:41:00Z.tif"
- type "image/tiff; application=geotiff"
- title "Methane Column Enhancement"
- description "Sample of methane plumes from point sources"
proj:bbox[] 4 items
- 0 -104.40938696034158
- 1 26.06828954903695
- 2 -104.11184817858141
- 3 26.266648736877066
- proj:wkt2 "GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4326"]]"
proj:shape[] 2 items
- 0 14
- 1 21
raster:bands[] 1 items
0
- scale 1.0
- nodata 0.0
- offset 0.0
- sampling "area"
- data_type "float64"
histogram
- max 0.3214902157450974
- min -0.545543008550919
- count 11
buckets[] 10 items
- 0 2
- 1 1
- 2 0
- 3 1
- 4 0
- 5 6
- 6 60
- 7 40
- 8 7
- 9 1
statistics
- mean 0.03760829813286496
- stddev 0.10854316578818861
- maximum 0.3214902157450974
- minimum -0.545543008550919
- valid_percent 40.136054421768705
proj:geometry
- type "Polygon"
coordinates[] 1 items
0[] 5 items
0[] 2 items
- 0 -104.40938696034158
- 1 26.06828954903695
1[] 2 items
- 0 -104.11184817858141
- 1 26.06828954903695
2[] 2 items
- 0 -104.11184817858141
- 1 26.266648736877066
3[] 2 items
- 0 -104.40938696034158
- 1 26.266648736877066
4[] 2 items
- 0 -104.40938696034158
- 1 26.06828954903695
proj:transform[] 9 items
- 0 0.01416851341715102
- 1 0.0
- 2 -104.40938696034158
- 3 0.0
- 4 -0.014168513417151129
- 5 26.266648736877066
- 6 0.0
- 7 0.0
- 8 1.0
roles[] 2 items
- 0 "data"
- 1 "layer"
rendered_preview
- href "https://earth.gov/ghgcenter/api/raster/collections/goes-ch4plume-v1/items/goes-ch4plume-v1-GOES-CH4_Mexico_Durango_BV2_BV2-5_2019-05-12T19:41:00Z/preview.png?assets=rad&rescale=0.0%2C0.4&colormap_name=plasma"
- type "image/png"
- title "Rendered preview"
- rel "preview"
roles[] 1 items
- 0 "overview"
- collection "goes-ch4plume-v1"
# Restructure our items into a dictionary where keys are the datetime items
# Then we can query more easily by date/time, e.g. "2020"
= {item.properties["datetime"]: item for item in items} items_dict
# Before we go further, let's pick which asset to focus on for the remainder of the notebook.
# This dataset only has one variable, so our asset is:
= "rad" asset_name
Creating Maps Using Folium
You will now visualize the evolution of a methane plume across for three different time steps using folium
.
Fetch Imagery Using Raster API
Here we get information from the Raster API
which we will add to our map in the next section.
# Take a look at the datetime options we have to choose from based on our search
list(items_dict.keys())
['2019-05-12T19:41:00Z',
'2019-05-12T19:36:00+00:00',
'2019-05-12T19:31:00+00:00',
'2019-05-12T19:26:00+00:00',
'2019-05-12T19:21:00+00:00',
'2019-05-12T19:16:00+00:00',
'2019-05-12T19:11:00+00:00',
'2019-05-12T19:06:00+00:00',
'2019-05-12T19:01:00+00:00',
'2019-05-12T18:56:00+00:00',
'2019-05-12T18:51:00+00:00',
'2019-05-12T18:46:00+00:00',
'2019-05-12T18:41:00+00:00',
'2019-05-12T18:36:00+00:00',
'2019-05-12T18:31:00+00:00',
'2019-05-12T18:26:00+00:00',
'2019-05-12T18:21:00+00:00',
'2019-05-12T18:16:00+00:00',
'2019-05-12T18:11:00+00:00',
'2019-05-12T18:06:00+00:00',
'2019-05-12T18:01:00+00:00',
'2019-05-12T17:56:00+00:00',
'2019-05-12T17:51:00+00:00',
'2019-05-12T17:46:00+00:00',
'2019-05-12T17:41:00+00:00',
'2019-05-12T17:36:00+00:00',
'2019-05-12T17:31:00+00:00',
'2019-05-12T17:26:00+00:00',
'2019-05-12T17:21:00+00:00',
'2019-05-12T17:16:00+00:00',
'2019-05-12T17:11:00+00:00',
'2019-05-12T17:06:00+00:00',
'2019-05-12T17:01:00+00:00',
'2019-05-12T16:56:00+00:00',
'2019-05-12T16:51:00+00:00',
'2019-05-12T16:46:00+00:00',
'2019-05-12T16:41:00+00:00',
'2019-05-11T16:11:00+00:00',
'2019-05-11T16:06:00+00:00',
'2019-05-11T16:01:00+00:00',
'2019-05-11T15:56:00+00:00',
'2019-05-11T15:51:00+00:00',
'2019-05-11T15:46:00+00:00',
'2019-05-11T15:41:00+00:00',
'2019-05-11T15:36:00+00:00',
'2019-05-11T15:31:00+00:00',
'2019-05-11T15:26:00+00:00',
'2019-05-11T15:21:00+00:00',
'2019-05-11T15:16:00+00:00',
'2019-05-11T15:11:00+00:00',
'2019-05-11T15:06:00+00:00',
'2019-05-11T15:01:00+00:00']
# Specify dates/times to visualize
# We'll select three observations, each 15 minutes apart
= ["2019-05-12T17:21:00+00:00","2019-05-12T17:36:00+00:00","2019-05-12T17:51:00+00:00"] datetimes
# Extract collection name and item ID for the first date
= items_dict[datetimes[0]]
observation_date_1 = observation_date_1.collection_id
collection_id = observation_date_1.id
item_id # Select relevant asset and get raster bands
object = observation_date_1.assets[asset_name]
= object.extra_fields.get("raster:bands", [{}])
raster_bands # Print the raster bands' information
raster_bands
[{'scale': 1.0,
'nodata': 0.0,
'offset': 0.0,
'sampling': 'area',
'data_type': 'float64',
'histogram': {'max': 0.27649252674814284,
'min': -0.10735534384080142,
'count': 11,
'buckets': [3, 4, 11, 19, 28, 19, 18, 6, 6, 4]},
'statistics': {'mean': 0.08132970903405966,
'stddev': 0.07654679032735916,
'maximum': 0.27649252674814284,
'minimum': -0.10735534384080142,
'valid_percent': 51.75438596491229}}]
# Use min and max to define desired colorbar range
= {
rescale_values "max": raster_bands[0]['statistics']['maximum'],
"min": raster_bands[0]['statistics']['minimum'],
}
print(rescale_values)
{'max': 0.27649252674814284, 'min': -0.10735534384080142}
Now, you will pass the item id
, collection name
, asset name
, and the rescale values
to the Raster API endpoint, along with a colormap. This step is required for every date/time you want to visualize, and tells the Raster API which collection, item, and asset you want to view, specifying the colormap and colorbar ranges to use for visualization. The API returns a JSON with information about the requested image. Each image will be referred to as a tile.
# Choose a colormap for displaying the data
# Make sure capitilization follows Matplotlib standard practice
# For more information on Colormaps in Matplotlib, please visit https://matplotlib.org/stable/users/explain/colors/colormaps.html
= "plasma" # please select the color ramp from matplotlib library. color_map
items.keys()
# Make a GET request to retrieve information for your first date/time
= requests.get(
observation_1_tile f"{RASTER_API_URL}/collections/{collection_id}/items/{item_id}/tilejson.json?"
f"&assets={asset_name}"
f"&color_formula=gamma+r+1.05&colormap_name={color_map}"
f"&rescale={rescale_values['min']},{rescale_values['max']}",
).json()# Print the resulting JSON
observation_1_tile
{'tilejson': '2.2.0',
'version': '1.0.0',
'scheme': 'xyz',
'tiles': ['https://earth.gov/ghgcenter/api/raster/collections/goes-ch4plume-v1/items/goes-ch4plume-v1-GOES-CH4_Mexico_Durango_BV2_BV2-5_2019-05-12T17:21:00Z/tiles/WebMercatorQuad/{z}/{x}/{y}@1x?assets=rad&color_formula=gamma+r+1.05&colormap_name=plasma&rescale=-0.10735534384080142%2C0.27649252674814284'],
'minzoom': 0,
'maxzoom': 24,
'bounds': [-104.38926051070443,
26.073879500563976,
-104.11975109944395,
26.24409597083376],
'center': [-104.25450580507419, 26.15898773569887, 0]}
# Repeat the previous step, fetching the tiles for the second and third datetimes of interest
= items_dict[datetimes[1]]
observation_date_2 = observation_date_2.collection_id
collection_id = observation_date_2.id
item_id = requests.get(
observation_2_tile f"{RASTER_API_URL}/collections/{collection_id}/items/{item_id}/tilejson.json?"
f"&assets={asset_name}"
f"&color_formula=gamma+r+1.05&colormap_name={color_map}"
f"&rescale={rescale_values['min']},{rescale_values['max']}",
).json()
= items_dict[datetimes[2]]
observation_date_3 = observation_date_3.collection_id
collection_id = observation_date_3.id
item_id = requests.get(
observation_3_tile f"{RASTER_API_URL}/collections/{collection_id}/items/{item_id}/tilejson.json?"
f"&assets={asset_name}"
f"&color_formula=gamma+r+1.05&colormap_name={color_map}"
f"&rescale={rescale_values['min']},{rescale_values['max']}",
).json()
Generate Map
# Set initial zoom and center of map
# We'll use the "center" variable from our first methane plume tile to set the center of the map
# Note that we specify "tiles=None" because in the next step we're going to set a custom tile to serve as our underlying world map.
= folium.Map(
map_ =(observation_1_tile["center"][1], observation_1_tile["center"][0]),
location=9,
zoom_start=None
tiles
)
# Specify a custom imagery source for the underlying map
folium.TileLayer(='https://server.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer/tile/{z}/{y}/{x}.png',
tiles='ESRI World Imagery',
name='Tiles © Esri — Source: Esri, i-cubed, USDA, USGS, AEX, GeoEye, Getmapping, Aerogrid, IGN, IGP, UPR-EGP, and the GIS User Community',
attr='True').add_to(map_)
overlay
# Add place labels on top
='https://server.arcgisonline.com/arcgis/rest/services/Reference/World_Boundaries_and_Places/MapServer/tile/{z}/{y}/{x}.png',
folium.TileLayer(tiles='ESRI World Boundaries and Places',
name='Tiles © Esri — Source: Esri, i-cubed, USDA, USGS, AEX, GeoEye, Getmapping, Aerogrid, IGN, IGP, UPR-EGP, and the GIS User Community',
attr='True').add_to(map_)
overlay
# Use the 'TileLayer' library to create the map layer for your first date/time tile
= TileLayer(
map_layer =observation_1_tile["tiles"][0], # Path to retrieve the tile
tiles=f'{items_dict[datetimes[0]].assets[asset_name].title}, {items_dict[datetimes[0]].properties["datetime"]}', # Give this layer a title
name='True', # The layer can be overlaid on the map
overlay="GHG", # Set the attribution
attr=1, # Adjust the transparency of the layer
opacity
)# Add layer to the map
map_layer.add_to(map_)
# Repeat for second date/time tile
= TileLayer(
map_layer =observation_2_tile["tiles"][0], # Path to retrieve the tile
tiles=f'{items_dict[datetimes[1]].assets[asset_name].title}, {items_dict[datetimes[1]].properties["datetime"]}', # Give this layer a title
name='True', # The layer can be overlaid on the map
overlay="GHG", # Set the attribution
attr=1, # Adjust the transparency of the layer
opacity
)
map_layer.add_to(map_)
# And finally for the third date/time tile
= TileLayer(
map_layer =observation_3_tile["tiles"][0], # Path to retrieve the tile
tiles=f'{items_dict[datetimes[2]].assets[asset_name].title}, {items_dict[datetimes[2]].properties["datetime"]}', # Give this layer a title
name='True', # The layer can be overlaid on the map
overlay="GHG", # Set the attribution
attr=1, # Adjust the transparency of the layer
opacity
)
map_layer.add_to(map_)
# Add controls to switch elements of the map on/off
=False, position='topright').add_to(map_)
folium.LayerControl(collapsed
# Add colorbar
# We can use 'generate_html_colorbar' from the 'ghgc_utils' module
# to create an HTML colorbar representation.
= ghgc_utils.generate_html_colorbar(
legend_html
color_map,
rescale_values,='Methane Column Enhancement (mol/m2)'
label
)
# Add colorbar to the map
map_.get_root().html.add_child(folium.Element(legend_html))
# Visualizing the map
map_
Click through the three methane layers to see the evolution of this plume over time!
# Set initial zoom and center of map for CH₄ Layer
# Centre of map [latitude,longitude]
= folium.Map(location=(26.29, -104.53), zoom_start=10)
map_
# April 2019
= TileLayer(
map_layer_2019 =april_2019_tile["tiles"][0],
tiles="GHG",
attr=0.7,
opacity
)
map_layer_2019.add_to(map_) map_
Summary
In this notebook we have successfully completed the following steps for the STAC collection for the GOSAT-based Top-down Total and Natural Methane Emissions dataset.
- Install and import the necessary libraries
- Fetch the collection from STAC collections using the appropriate endpoints
- Count the number of existing granules within the collection
- Map the methane emission levels
If you have any questions regarding this user notebook, please contact us using the feedback form.