Declarative Satellite Data

from datetime import datetime

# Create variables for URL generation
# YOUR CODE GOES HERE

# Construct the data_url string
# YOUR CODE GOES HERE

# Print out your URL and verify it works!
# YOUR CODE GOES HERE

# %load solutions/data_url.py

# Cell content replaced by load magic replacement.
from datetime import datetime

# Create variables for URL generation
image_date = datetime.utcnow().date()
region = 'CONUS'
channel = 8

# We want to match something like:
# https://thredds-test.unidata.ucar.edu/thredds/catalog/satellite/goes16/GOES16/Mesoscale-1/Channel08/20181113/catalog.html

# Construct the data_url string
data_url = ('https://thredds.ucar.edu/thredds/catalog/satellite/goes/east/products/'
            f'CloudAndMoistureImagery/{region}/Channel{channel:02d}/'
            f'{image_date:%Y%m%d}/catalog.xml')

# Print out your URL and verify it works!
print(data_url)

from siphon.catalog import TDSCatalog

image_date = datetime.utcnow().date()
region = 'CONUS'
channel = 8

data_url = ('https://thredds.ucar.edu/thredds/catalog/satellite/goes/east/products/'
            f'CloudAndMoistureImagery/{region}/Channel{channel:02d}/'
            f'{image_date:%Y%m%d}/catalog.xml')

cat = TDSCatalog(data_url)

cat.datasets[:5]

cat.datasets[0]

dataset = cat.datasets[1]

print(dataset)

ds = dataset.remote_access(use_xarray=True)

ds

import xarray as xr

from metpy.plots.declarative import *

img = ImagePlot()
img.data = ds
img.field = 'Sectorized_CMI'

panel = MapPanel()
panel.plots = [img]

pc = PanelContainer()
pc.panels = [panel]

pc.show()

# Create the ImagePlot
# YOUR CODE GOES HERE

# Create the MapPanel
# YOUR CODE GOES HERE

# Create the PanelContainer
# YOUR CODE GOES HERE

# Show it!
# YOUR CODE GOES HERE

# %load solutions/dec_map_modifications.py

# Cell content replaced by load magic replacement.
# Create the ImagePlot
img = ImagePlot()
img.data = ds
img.field = 'Sectorized_CMI'
img.colormap = 'WVCIMSS'

# Create the MapPanel
panel = MapPanel()
panel.plots = [img]
panel.title = "CONUS GOES-East"

# Create the PanelContainer
pc = PanelContainer()
pc.panels = [panel]
pc.size = (10, 8)

# Show it!
pc.show()

#
# Get the data
#

# Get the data for band 7
# YOUR CODE GOES HERE

# Get the data for band 10
# YOUR CODE GOES HERE

#
# Make the plot
#
# YOUR CODE GOES HERE

# %load solutions/dec_map_multipanel.py

# Cell content replaced by load magic replacement.
#
# Get the data
#

# Get the data for band 7
url = ('https://thredds.ucar.edu/thredds/catalog/satellite/goes/east/products/'
            f'CloudAndMoistureImagery/{region}/Channel07/'
            f'{datetime.utcnow():%Y%m%d}/catalog.xml')
cat = TDSCatalog(url)
band_7 = cat.datasets[-2].remote_access(use_xarray=True)

# Get the data for band 10
url = ('https://thredds.ucar.edu/thredds/catalog/satellite/goes/east/products/'
            f'CloudAndMoistureImagery/{region}/Channel10/'
            f'{datetime.utcnow():%Y%m%d}/catalog.xml')
cat = TDSCatalog(url)
band_10 = cat.datasets[-2].remote_access(use_xarray=True)

#
# Make the plot
#
band_7_img = ImagePlot()
band_7_img.data = band_7
band_7_img.field = 'Sectorized_CMI'
band_7_img.colormap = 'ir_rgbv'

band_10_img = ImagePlot()
band_10_img.data = band_10
band_10_img.field = 'Sectorized_CMI'
band_10_img.colormap = 'WVCIMSS'

left_panel = MapPanel()
left_panel.plots = [band_7_img]
left_panel.title = 'IR - GOES EAST'
left_panel.layout = (1, 2, 1)

right_panel = MapPanel()
right_panel.plots = [band_10_img]
right_panel.title = 'Mid-Level Water Vapor - GOES EAST'
right_panel.layout = (1, 2, 2)

pc = PanelContainer()
pc.panels = [left_panel, right_panel]
pc.size = (12, 10)

pc.show()

# %load solutions/dec_map_multipanel_bonus.py

# Cell content replaced by load magic replacement.
def make_channel_image_plot(channel):
    url = ('https://thredds.ucar.edu/thredds/catalog/satellite/goes/east/products/'
                f'CloudAndMoistureImagery/{region}/Channel{channel:02d}/'
                f'{datetime.utcnow():%Y%m%d}/catalog.xml')
    cat = TDSCatalog(url)
    band_data = cat.datasets[-2].remote_access(use_xarray=True)
    
    img = ImagePlot()
    img.data = band_data
    img.field = 'Sectorized_CMI'
    return img

#
# Make the plot
#
band_7_img = make_channel_image_plot(7)
band_10_img = make_channel_image_plot(10)

band_7_img.colormap = 'ir_rgbv'
band_10_img.colormap = 'WVCIMSS'

left_panel = MapPanel()
left_panel.plots = [band_7_img]
left_panel.title = 'IR - GOES EAST'
left_panel.layout = (1, 2, 1)

right_panel = MapPanel()
right_panel.plots = [band_10_img]
right_panel.title = 'Mid-Level Water Vapor - GOES EAST'
right_panel.layout = (1, 2, 2)

pc = PanelContainer()
pc.panels = [left_panel, right_panel]
pc.size = (12, 10)

pc.show()

img = ImagePlot()
img.data = ds
img.field = 'Sectorized_CMI'
img.colormap = 'WVCIMSS'

# Create the MapPanel
panel = MapPanel()
panel.plots = [img]
panel.title = "CONUS GOES-East"
panel.area = 'ma'
panel.layers = ['coastline', 'borders', 'states']

# Create the PanelContainer
pc = PanelContainer()
pc.panels = [panel]
pc.size = (10, 8)

# Show it!
pc.show()

import cartopy.crs as ccrs

# Explore ccrs with tab completion!

from metpy.plots import add_timestamp, add_metpy_logo

print(ds.start_date_time)

start_time = datetime.strptime(ds.start_date_time, '%Y%j%H%M%S')

img = ImagePlot()
img.data = ds
img.field = 'Sectorized_CMI'
img.colormap = 'WVCIMSS'

# Create the MapPanel
panel = MapPanel()
panel.plots = [img]
panel.title = "CONUS GOES-East"
panel.area = 'us'
panel.layers = ['coastline', 'borders', 'states']

# Create the PanelContainer
pc = PanelContainer()
pc.panels = [panel]
pc.size = (10, 8)

# Add Annotations
add_timestamp(panel.ax, start_time)

# Show it!
pc.show()

img = ImagePlot()
img.data = ds
img.field = 'Sectorized_CMI'
img.colormap = 'WVCIMSS'

# Create the MapPanel
panel = MapPanel()
panel.plots = [img]
panel.title = "CONUS GOES-East"
panel.area = 'us'
panel.layers = ['coastline', 'borders', 'states']

# Create the PanelContainer
pc = PanelContainer()
pc.panels = [panel]
pc.size = (10, 8)

# Add Annotations
add_timestamp(panel.ax, start_time, high_contrast=True,
              pretext='Valid: ', y=0.02)

add_metpy_logo(pc.figure, y=10)

# Show it!
pc.show()