mirror of
https://github.com/Unidata/python-awips.git
synced 2025-02-23 22:57:56 -05:00
deploy: 8db6e4951c
This commit is contained in:
tiffanycmeyer13
parent
f48336965d
commit
ee756bc4d3
28 changed files with 1312 additions and 441 deletions
|
|
@ -1,4 +1,4 @@
|
|||
# Sphinx build info version 1
|
||||
# This file hashes the configuration used when building these files. When it is not found, a full rebuild will be done.
|
||||
config: f0ebb3dbf3d9374c81a3df5a307e45ee
|
||||
config: f9ff9fe8bb30e963804480f00bc57e27
|
||||
tags: 645f666f9bcd5a90fca523b33c5a78b7
|
||||
|
|
|
|||
BIN
_images/Watch_Warning_and_Advisory_Plotting_34_0.png
Normal file
BIN
_images/Watch_Warning_and_Advisory_Plotting_34_0.png
Normal file
Binary file not shown.
|
After Width: | Height: | Size: 182 KiB |
Binary file not shown.
|
Before Width: | Height: | Size: 109 KiB |
|
|
@ -0,0 +1,576 @@
|
|||
===================================
|
||||
Watch Warning and Advisory Plotting
|
||||
===================================
|
||||
`Notebook <http://nbviewer.ipython.org/github/Unidata/python-awips/blob/master/examples/notebooks/Watch_Warning_and_Advisory_Plotting.ipynb>`_
|
||||
Python-AWIPS Tutorial Notebook
|
||||
|
||||
--------------
|
||||
|
||||
Objectives
|
||||
==========
|
||||
|
||||
- Create a colorized plot with `Warnings, Watches, Advisories and
|
||||
Statements (WWAs) <https://weather.cod.edu/notes/criteria/>`__
|
||||
- Use python-awips to connect to an EDEX server
|
||||
- Create and filter the data request specifically for a warning data
|
||||
type
|
||||
- Create and use accurate time filter for data requests
|
||||
- Define and use functions
|
||||
- Define and use dictionaries
|
||||
- Colorize shapes based on a dictionary
|
||||
- Overlay warnings, watches, and advisories with state and political
|
||||
maps
|
||||
|
||||
--------------
|
||||
|
||||
Table of Contents
|
||||
-----------------
|
||||
|
||||
| `1
|
||||
Imports <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#imports>`__\
|
||||
| `2 Functions:
|
||||
make_map() <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#functions-make-map>`__\
|
||||
| `3 Functions:
|
||||
get_color() <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#functions-get-color>`__\
|
||||
| `4 Functions:
|
||||
get_title() <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#functions-get-title>`__\
|
||||
| `5 Initial
|
||||
Setup <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#initial-setup>`__\
|
||||
| `5.1 EDEX
|
||||
Connection <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#edex-connection>`__\
|
||||
| `5.2 Significance (Sig)
|
||||
Constraints <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#significance-sig-constraints>`__\
|
||||
| `6 Filter by
|
||||
Time <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#filter-by-time>`__\
|
||||
| `7 Use the
|
||||
Data! <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#use-the-data>`__\
|
||||
| `7.1 Get the
|
||||
Data <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#get-the-data>`__\
|
||||
| `7.2 Extract Phensigs, Geometries, and
|
||||
Times <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#extract-phensigs-geometries-and-times>`__\
|
||||
| `8 Plot the
|
||||
Data! <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#plot-the-data>`__\
|
||||
| `8.1 Create State and Political
|
||||
Boundaries <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#create-state-and-political-boundaries>`__\
|
||||
| `8.2 Draw the Plot and Legend for
|
||||
WWAs <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#draw-the-plot-and-legend-for-wwas>`__\
|
||||
| `9 See
|
||||
Also <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#see-also>`__\
|
||||
| `9.1 Related
|
||||
Notebooks <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#related-notebooks>`__\
|
||||
| `9.2 Additional
|
||||
Documentation <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#additional-documentation>`__\
|
||||
|
||||
1 Imports
|
||||
---------
|
||||
|
||||
The imports below are used throughout the notebook. The python-awips
|
||||
imports allow us to connect to an EDEX server, use the warning lookup
|
||||
dictionary, and define a TimeRange. The additional imports are for data
|
||||
manipulation and visualization.
|
||||
|
||||
.. code:: ipython3
|
||||
|
||||
from datetime import datetime, timedelta
|
||||
|
||||
import numpy as np
|
||||
import matplotlib.patches as mpatches
|
||||
import matplotlib.pyplot as plt
|
||||
import cartopy.crs as ccrs
|
||||
import cartopy.feature as cfeature
|
||||
from cartopy.feature import ShapelyFeature, NaturalEarthFeature
|
||||
from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
|
||||
from shapely.geometry import MultiPolygon, Polygon
|
||||
|
||||
from awips.dataaccess import DataAccessLayer
|
||||
from awips.tables import vtec
|
||||
from dynamicserialize.dstypes.com.raytheon.uf.common.time import TimeRange
|
||||
|
||||
`Top <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html>`__
|
||||
|
||||
--------------
|
||||
|
||||
2 Function: make_map()
|
||||
----------------------
|
||||
|
||||
In order to plot more than one image, it’s easiest to define common
|
||||
logic in a function. However, for this notebook we only use it in one
|
||||
place. It is a function you will find in most of our example notebooks.
|
||||
|
||||
Here, a new function called **make_map** is defined. This function uses
|
||||
the `matplotlib.pyplot package
|
||||
(plt) <https://matplotlib.org/3.3.3/api/_as_gen/matplotlib.pyplot.html>`__
|
||||
to create a figure and axis. The lat/lon grids are added.
|
||||
|
||||
.. code:: ipython3
|
||||
|
||||
def make_map(bbox, projection=ccrs.PlateCarree()):
|
||||
fig, ax = plt.subplots(figsize=(20,12),
|
||||
subplot_kw=dict(projection=projection))
|
||||
ax.set_extent(bbox)
|
||||
gl = ax.gridlines(draw_labels=True)
|
||||
gl.top_labels = gl.right_labels = False
|
||||
gl.xformatter = LONGITUDE_FORMATTER
|
||||
gl.yformatter = LATITUDE_FORMATTER
|
||||
return fig, ax
|
||||
|
||||
`Top <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html>`__
|
||||
|
||||
--------------
|
||||
|
||||
3 Function: get_color()
|
||||
-----------------------
|
||||
|
||||
Since we’ll be needing to access the color using the vtec lookup table
|
||||
in several places, creating an easily recognizable function is useful.
|
||||
|
||||
.. code:: ipython3
|
||||
|
||||
def get_color(phensig):
|
||||
return vtec[phensig]['color']
|
||||
|
||||
`Top <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html>`__
|
||||
|
||||
--------------
|
||||
|
||||
4 Function get_title()
|
||||
----------------------
|
||||
|
||||
Similar to the color function just defined, accessing the full name for
|
||||
the phensig will also be necessary, so this function will be helpful.
|
||||
|
||||
.. code:: ipython3
|
||||
|
||||
def get_title(phensig):
|
||||
return vtec[phensig]['hdln']
|
||||
|
||||
`Top <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html>`__
|
||||
|
||||
--------------
|
||||
|
||||
5 Inital Setup
|
||||
--------------
|
||||
|
||||
5.1 EDEX Connection
|
||||
~~~~~~~~~~~~~~~~~~~
|
||||
|
||||
First we establish a connection to Unidata’s public EDEX server. With
|
||||
that connection made, we can create a `new data request
|
||||
object <http://unidata.github.io/python-awips/api/IDataRequest.html>`__
|
||||
and set the data type to **warning**, and set the Parameters to
|
||||
**phensig** and **sig**.
|
||||
|
||||
.. container:: alert-info
|
||||
|
||||
Note: Remember, to see all available parameters use the
|
||||
DataAccess.getAvailableParameters() method as shown in the Grid
|
||||
Levels and Parameters Notebook.
|
||||
|
||||
.. code:: ipython3
|
||||
|
||||
DataAccessLayer.changeEDEXHost("edex-cloud.unidata.ucar.edu")
|
||||
request = DataAccessLayer.newDataRequest()
|
||||
request.setDatatype("warning")
|
||||
params = ["phensig", "sig"]
|
||||
request.setParameters(*(params))
|
||||
|
||||
5.2 Significance (Sig) Constants
|
||||
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
|
||||
|
||||
The two parameters we’re requesting for our warning objects are
|
||||
**phensig** and **sig** where phensig is styled “XX.Y” and sig is “Y”.
|
||||
Phen stands for “Phenomena” and sig stands for “Significance”. `A more
|
||||
detailed description of phensigs and how they’re used is provided with
|
||||
this NWS
|
||||
pamphlet <https://www.weather.gov/media/vtec/VTEC_explanation4-20.pdf>`__.
|
||||
|
||||
The constants in this section correlate the **sig** to what type of
|
||||
message it is (what significance it is).
|
||||
|
||||
.. code:: ipython3
|
||||
|
||||
WATCH_SIG = 'A'
|
||||
WARN_SIG = 'W'
|
||||
ADVIS_SIG = 'Y'
|
||||
STATEM_SIG = 'S'
|
||||
|
||||
`Top <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html>`__
|
||||
|
||||
--------------
|
||||
|
||||
6 Filter by Time
|
||||
----------------
|
||||
|
||||
Here we decide how much data we want to pull from EDEX. By default we’ll
|
||||
request 12 hours, but that value can easily be modified by `adjusting
|
||||
the
|
||||
``timedelta(hours = 12)`` <https://docs.python.org/3/library/datetime.html#timedelta-objects>`__
|
||||
in line ``2``. The more data we request, the longer the next section
|
||||
will take to run.
|
||||
|
||||
.. code:: ipython3
|
||||
|
||||
# Get records from the last 12 hours
|
||||
lastHourDateTime = datetime.utcnow() - timedelta(hours = 12)
|
||||
start = lastHourDateTime.strftime('%Y-%m-%d %H:%M:%S')
|
||||
end = datetime.utcnow().strftime('%Y-%m-%d %H:%M:%S')
|
||||
|
||||
beginRange = datetime.strptime( start , "%Y-%m-%d %H:%M:%S")
|
||||
endRange = datetime.strptime( end , "%Y-%m-%d %H:%M:%S")
|
||||
timerange = TimeRange(beginRange, endRange)
|
||||
|
||||
`Top <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html>`__
|
||||
|
||||
--------------
|
||||
|
||||
7 Use the Data!
|
||||
---------------
|
||||
|
||||
7.1 Get the Data
|
||||
~~~~~~~~~~~~~~~~
|
||||
|
||||
Now that we have our ``request`` and TimeRange ``timerange`` objects
|
||||
ready, it’s time to request the data array from EDEX.
|
||||
|
||||
.. container:: alert-info
|
||||
|
||||
Note: Above we set timerange to be 12 hours worth of data. This can
|
||||
return on the order of ~2000 records and can take a little while to
|
||||
run.
|
||||
|
||||
.. code:: ipython3
|
||||
|
||||
# Get response
|
||||
response = DataAccessLayer.getGeometryData(request, timerange)
|
||||
print("Using " + str(len(response)) + " records")
|
||||
|
||||
|
||||
.. parsed-literal::
|
||||
|
||||
Using 1502 records
|
||||
|
||||
|
||||
7.2 Extract Phensigs, Geometries, and Times
|
||||
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
|
||||
|
||||
In this section we start gathering all the information we’ll need to
|
||||
properly display our data. First we create an array to keep track of
|
||||
unique phensigs. This is useful summary information and will be used to
|
||||
help create the legend which we’ll display along with our plot.
|
||||
|
||||
Next, we create arrays for each of the 4 types of significance a
|
||||
statement can have. We will group our records this way, so we can easily
|
||||
toggle which records to display or not.
|
||||
|
||||
Then, we create two time variables to keep track of the earliest time
|
||||
from our records and the latest time, and will display that information
|
||||
in the title of our plot.
|
||||
|
||||
This section has optional print statements at lines ``65`` and ``85``.
|
||||
The first prints out the title, phensig, ref time, and shape for each
|
||||
unique phensig, and the second prints out a sum of how many unique
|
||||
phensigs there are.
|
||||
|
||||
We cycle through all the data produced from our ``response`` object,
|
||||
access its geometries, and create a new
|
||||
`ShapelyFeature <https://scitools.org.uk/cartopy/docs/latest/reference/generated/cartopy.feature.ShapelyFeature.html>`__
|
||||
with the corresponding color. Then we place this new feature in the
|
||||
appropriate ``shapes`` array. During this process we also populate the
|
||||
phensigs array with all unique phensig entries.
|
||||
|
||||
Finally, after we’re done looping through all the ``response`` data, we
|
||||
create a mapping of phensigs to their corresponding titles. This will be
|
||||
used later to sort the legend alphabetically by titles (which differs
|
||||
from simply sorting by phensig). Ex. *Blizzard Warning (BZ.W)* would
|
||||
come before *Areal Flood Advisory (FA.Y)* if we simply sorted by
|
||||
phensig.
|
||||
|
||||
.. code:: ipython3
|
||||
|
||||
# Keep track of unique phensigs, to use in legend
|
||||
phensigs = []
|
||||
|
||||
# Sort the geometries based on their sig
|
||||
watch_shapes = []
|
||||
warning_shapes = []
|
||||
advisory_shapes = []
|
||||
statement_shapes = []
|
||||
|
||||
# Keep track of the earliest and latest reftime for title
|
||||
# start with the first time from the first object in the response
|
||||
time_str = str(response[0].getDataTime().getRefTime())
|
||||
# truncate the decimal seconds for datetime parsing
|
||||
time_str = time_str[:-4]
|
||||
# convert to datetime object for easy comparison
|
||||
first_time = datetime.strptime(time_str, '%Y-%m-%d %H:%M:%S')
|
||||
last_time = datetime.strptime(time_str, '%Y-%m-%d %H:%M:%S')
|
||||
|
||||
for ob in response:
|
||||
|
||||
# get the geometry for the object
|
||||
poly = ob.getGeometry()
|
||||
# get the reftime for the object
|
||||
ref = ob.getDataTime().getRefTime()
|
||||
|
||||
# do not plot if phensig is blank (SPS)
|
||||
if ob.getString('phensig'):
|
||||
|
||||
# look at the reftime
|
||||
# convert reftime to a string and parse the decimal seconds
|
||||
ref_str = str(ref)
|
||||
ref_str = ref_str[:-4]
|
||||
# convert reftime to a datetime object for comparison
|
||||
ref_time = datetime.strptime(ref_str, '%Y-%m-%d %H:%M:%S')
|
||||
# compare current time with first and last times and set if appropriate
|
||||
if ref_time is not None:
|
||||
if ref_time < first_time:
|
||||
first_time = ref_time
|
||||
elif ref_time > last_time:
|
||||
last_time = ref_time
|
||||
|
||||
# get the phensig and sig values from object
|
||||
phensigString = ob.getString('phensig')
|
||||
sig = ob.getString('sig')
|
||||
|
||||
# set the geometries based on whether it's a MultiPolygon or Polygon
|
||||
if poly.geom_type == 'MultiPolygon':
|
||||
geometries = np.array([])
|
||||
geometries = np.append(geometries,MultiPolygon(poly))
|
||||
else:
|
||||
geometries = np.array([])
|
||||
geometries = np.append(geometries,Polygon(poly))
|
||||
|
||||
for geom in geometries:
|
||||
bounds = Polygon(geom)
|
||||
intersection = bounds.intersection
|
||||
geoms = (intersection(geom) for geom in geometries if bounds.intersects(geom))
|
||||
|
||||
# Store the unique phensigs
|
||||
if not phensigString in phensigs:
|
||||
phensigs.append(phensigString)
|
||||
# Optional printout of unique Phensigs
|
||||
# print(get_title(phensigString) + " (" + phensigString + ")
|
||||
|
||||
# get the corresponding color using the dictionary
|
||||
color = get_color(phensigString)
|
||||
# create a new shape feature for the object
|
||||
shape_feature = ShapelyFeature(geoms,ccrs.PlateCarree(),
|
||||
facecolor=color, edgecolor=color)
|
||||
# store the shape feature in the correct array
|
||||
if sig is WARN_SIG:
|
||||
warning_shapes.append(shape_feature)
|
||||
elif sig is WATCH_SIG:
|
||||
watch_shapes.append(shape_feature)
|
||||
elif sig is ADVIS_SIG:
|
||||
advisory_shapes.append(shape_feature)
|
||||
elif sig is STATEM_SIG:
|
||||
statement_shapes.append(shape_feature)
|
||||
|
||||
# Optional printout for the number of unique phensigs
|
||||
print(len(phensigs), " Unique Phensigs")
|
||||
|
||||
# Map phensigs to their titles (used for displaying alphabetically by
|
||||
# title in legend)
|
||||
phensig_titles = {}
|
||||
for phensig in phensigs:
|
||||
key = get_title(phensig)
|
||||
phensig_titles[key] = phensig
|
||||
|
||||
|
||||
.. parsed-literal::
|
||||
|
||||
14 Unique Phensigs
|
||||
|
||||
|
||||
`Top <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html>`__
|
||||
|
||||
--------------
|
||||
|
||||
8 Plot the Data!
|
||||
----------------
|
||||
|
||||
8.1 Create State and Political Boundaries
|
||||
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
|
||||
|
||||
Define the state and political boundaries that we’ll use in our plot to
|
||||
give more of a frame of reference. These objects are standard method
|
||||
calls in the `Cartopy Feature package, using the NaturalEarthFeature
|
||||
function <https://scitools.org.uk/cartopy/docs/v0.14/matplotlib/feature_interface.html#cartopy.feature.NaturalEarthFeature>`__.
|
||||
|
||||
.. code:: ipython3
|
||||
|
||||
# Define the state and political boundaries for the plot
|
||||
states_provinces = cfeature.NaturalEarthFeature(
|
||||
category='cultural',
|
||||
name='admin_1_states_provinces_lines',
|
||||
scale='50m',
|
||||
facecolor='none')
|
||||
political_boundaries = cfeature.NaturalEarthFeature(category='cultural',
|
||||
name='admin_0_boundary_lines_land',
|
||||
scale='50m', facecolor='none')
|
||||
|
||||
8.2 Draw the Plot and Legend for WWAs
|
||||
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
|
||||
|
||||
Here is where we finally get ot draw something! The very first few lines
|
||||
of this section are constants that we can manually “switch on and off”
|
||||
for what records we want displayed. By default we have all significance
|
||||
types drawn. If we want to “turn off” any of the significance records,
|
||||
simply set it’s corresponding constant to false, and re-run this cell to
|
||||
see how that plot compares.
|
||||
|
||||
The next step involves creating the objects that are used to define the
|
||||
legend. We use the ``phensig_titles`` dictionary to loop through all the
|
||||
phensigs in alphabetical (by title) order. Then, we compare if the
|
||||
phensig will be displayed or not based on the display constants from the
|
||||
previous lines. If the significance will be drawn then we create a new
|
||||
`Patch
|
||||
object <https://matplotlib.org/stable/api/_as_gen/matplotlib.patches.Patch.html#matplotlib.patches.Patch>`__
|
||||
of the corresponding color with the corresponding label and add it to
|
||||
our ``handles`` array.
|
||||
|
||||
After that we define our bounding box and create our new plot with its
|
||||
figure and axes.
|
||||
|
||||
Our next step is to create our Title for our plot. We create a title
|
||||
based on the draw variables to accurately describe what is being drawn
|
||||
in our plot. Here is where we use the first and last times defined in a
|
||||
previous cell.
|
||||
|
||||
Finally, we create and show our plot. We add the title to the plot, add
|
||||
all the features to the axes, and add the legend as well.
|
||||
|
||||
.. code:: ipython3
|
||||
|
||||
# Set these variables for which records to draw
|
||||
DRAW_ADVISORY = True
|
||||
DRAW_WATCH = True
|
||||
DRAW_WARNING = True
|
||||
DRAW_STATEMENT = True
|
||||
|
||||
# Create handles for legend and add items alphabetically by title and
|
||||
# only display based on the display values above
|
||||
handles = []
|
||||
for title in sorted(phensig_titles):
|
||||
phensig = phensig_titles[title]
|
||||
# check draw booleans
|
||||
if ( "."+ADVIS_SIG in phensig and DRAW_ADVISORY or
|
||||
"."+WATCH_SIG in phensig and DRAW_WATCH or
|
||||
"."+WARN_SIG in phensig and DRAW_WARNING or
|
||||
"."+STATEM_SIG in phensig and DRAW_STATEMENT ):
|
||||
entry = mpatches.Patch(color=get_color(phensig), label=title)
|
||||
handles.append(entry)
|
||||
|
||||
# Create the plot
|
||||
bbox=[-127,-64,24,49]
|
||||
fig, ax = make_map(bbox=bbox)
|
||||
|
||||
# Add the title
|
||||
# Construct the title based on which record types are being displayed
|
||||
title_string = ""
|
||||
if DRAW_WATCH:
|
||||
title_string += "Watches, "
|
||||
if DRAW_WARNING:
|
||||
title_string += "Warnings, "
|
||||
if DRAW_ADVISORY:
|
||||
title_string += "Advisories, "
|
||||
if DRAW_STATEMENT:
|
||||
title_string += "Statements, "
|
||||
# remove the last comma and space
|
||||
title_string = title_string[:-2]
|
||||
# add the time range
|
||||
title_string += " from " + str(first_time)[:-3] + " to " + str(last_time)[:-3] + " UTC"
|
||||
# set the title on the plot, give it a bigger font size, and increase
|
||||
# the vertical padding between the title and the figure
|
||||
plt.title(title_string, fontsize=20, pad=10)
|
||||
|
||||
# Draw all features on the plot
|
||||
ax.add_feature(cfeature.LAND)
|
||||
ax.add_feature(cfeature.COASTLINE)
|
||||
ax.add_feature(states_provinces, edgecolor='black')
|
||||
ax.add_feature(political_boundaries, edgecolor='black')
|
||||
|
||||
# Draw WWAs in order: Advisory -> Watch > Warning > Statement
|
||||
if DRAW_ADVISORY:
|
||||
for shape in advisory_shapes:
|
||||
ax.add_feature(shape)
|
||||
if DRAW_WATCH:
|
||||
for shape in watch_shapes:
|
||||
ax.add_feature(shape)
|
||||
if DRAW_WARNING:
|
||||
for shape in warning_shapes:
|
||||
ax.add_feature(shape)
|
||||
if DRAW_STATEMENT:
|
||||
for shape in statement_shapes:
|
||||
ax.add_feature(shape)
|
||||
|
||||
# Draw the legend
|
||||
# use the handles defined earlier for the color associations to
|
||||
# phensig titles, set the location to the lower center, give it
|
||||
# 5 columns so it uses all the horizonatal space, place it under
|
||||
# the actual figure, and give it a larger fontsize
|
||||
bottom = 0.12 + (len(handles)//5 *.04)
|
||||
ax.legend(handles=handles, loc='lower center', ncol=5, bbox_to_anchor=(0.5, -bottom), fontsize=16)
|
||||
|
||||
# Show the plot
|
||||
plt.show()
|
||||
|
||||
|
||||
|
||||
.. image:: Watch_Warning_and_Advisory_Plotting_files/Watch_Warning_and_Advisory_Plotting_34_0.png
|
||||
|
||||
|
||||
`Top <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html>`__
|
||||
|
||||
--------------
|
||||
|
||||
9 See Also
|
||||
----------
|
||||
|
||||
- `National Weather Service WWA Definitions (Baltimore
|
||||
Office) <https://www.weather.gov/lwx/warningsdefined#:~:text=A%20Winter%20Storm%20Warning%20is%20issued%20when%20a,combination%20of%20snow%20and%2For%20ice%20accumulation%20with%20wind.>`__
|
||||
- `College of Dupage WWA
|
||||
Definitions <https://weather.cod.edu/notes/criteria/>`__
|
||||
- `Phensig
|
||||
Explanation <https://www.weather.gov/media/vtec/VTEC_explanation4-20.pdf>`__
|
||||
|
||||
9.1 Related Notebooks
|
||||
~~~~~~~~~~~~~~~~~~~~~
|
||||
|
||||
- `Grid Levels and
|
||||
Parameters <http://unidata.github.io/python-awips/examples/generated/Grid_Levels_and_Parameters.html>`__
|
||||
|
||||
9.2 Additional Documentation
|
||||
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
|
||||
|
||||
**python-awips**
|
||||
|
||||
- `DataAccessLayer.changeEDEXHost() <http://unidata.github.io/python-awips/api/DataAccessLayer.html#awips.dataaccess.DataAccessLayer.changeEDEXHost>`__
|
||||
- `DataAccessLayer.newDataRequest() <http://unidata.github.io/python-awips/api/DataAccessLayer.html#awips.dataaccess.DataAccessLayer.newDataRequest>`__
|
||||
- `DataAccessLayer.getAvailableParameters() <http://unidata.github.io/python-awips/api/DataAccessLayer.html#awips.dataaccess.DataAccessLayer.getAvailableParameters>`__
|
||||
- `IDataRequest <http://unidata.github.io/python-awips/api/IDataRequest.html>`__
|
||||
- `GeometryData <http://unidata.github.io/python-awips/api/PyGeometryData.html>`__
|
||||
|
||||
**datetime**
|
||||
|
||||
- `datetime.datetime <https://docs.python.org/3/library/datetime.html#datetime-objects>`__
|
||||
- `datetime.timedelta <https://docs.python.org/3/library/datetime.html#timedelta-objects>`__
|
||||
|
||||
**cartopy**
|
||||
|
||||
- `cartopy feature
|
||||
interface <https://scitools.org.uk/cartopy/docs/v0.14/matplotlib/feature_interface.html>`__
|
||||
- `cartopy.feature.ShaeplyFeature <https://scitools.org.uk/cartopy/docs/latest/reference/generated/cartopy.feature.ShapelyFeature.html>`__
|
||||
|
||||
**matplotlib**
|
||||
|
||||
- `matplotlib.pyplot() <https://matplotlib.org/stable/api/_as_gen/matplotlib.pyplot.html>`__
|
||||
- `matplotlib.pyplot.legend() <https://matplotlib.org/3.1.1/api/_as_gen/matplotlib.pyplot.legend.html>`__
|
||||
- `matplotlib.pyplot.axes() <https://matplotlib.org/stable/api/_as_gen/matplotlib.pyplot.axes.html>`__
|
||||
- `matplotlib.pyplot.figure() <https://matplotlib.org/stable/api/_as_gen/matplotlib.pyplot.figure.html>`__
|
||||
- `matplotlib.pyplot.title() <https://matplotlib.org/3.1.1/api/_as_gen/matplotlib.pyplot.title.html>`__
|
||||
- `matplotlib.pathes.Patch <https://matplotlib.org/stable/api/_as_gen/matplotlib.patches.Patch.html#matplotlib.patches.Patch>`__
|
||||
|
||||
`Top <https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html>`__
|
||||
|
||||
--------------
|
||||
|
|
@ -1,139 +0,0 @@
|
|||
==========================
|
||||
Watch and Warning Polygons
|
||||
==========================
|
||||
`Notebook <http://nbviewer.ipython.org/github/Unidata/python-awips/blob/master/examples/notebooks/Watch_and_Warning_Polygons.ipynb>`_
|
||||
This example uses matplotlib, cartopy, shapely, and python-awips to plot
|
||||
watch and warning polygons requested from a real-time AWIPS EDEX server.
|
||||
|
||||
First, set up our imports and define functions to be used later:
|
||||
|
||||
.. code:: ipython3
|
||||
|
||||
from awips.dataaccess import DataAccessLayer
|
||||
from awips.tables import vtec
|
||||
from datetime import datetime
|
||||
import numpy as np
|
||||
import matplotlib.pyplot as plt
|
||||
import cartopy.crs as ccrs
|
||||
import cartopy.feature as cfeature
|
||||
from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
|
||||
from cartopy.feature import ShapelyFeature,NaturalEarthFeature
|
||||
from shapely.geometry import MultiPolygon,Polygon
|
||||
|
||||
def warning_color(phensig):
|
||||
return vtec[phensig]['color']
|
||||
|
||||
def make_map(bbox, projection=ccrs.PlateCarree()):
|
||||
fig, ax = plt.subplots(figsize=(20,12),
|
||||
subplot_kw=dict(projection=projection))
|
||||
ax.set_extent(bbox)
|
||||
gl = ax.gridlines(draw_labels=True)
|
||||
gl.top_labels = gl.right_labels = False
|
||||
gl.xformatter = LONGITUDE_FORMATTER
|
||||
gl.yformatter = LATITUDE_FORMATTER
|
||||
return fig, ax
|
||||
|
||||
Next, we create a request for the “warning” data type:
|
||||
|
||||
.. code:: ipython3
|
||||
|
||||
DataAccessLayer.changeEDEXHost("edex-cloud.unidata.ucar.edu")
|
||||
request = DataAccessLayer.newDataRequest()
|
||||
request.setDatatype("warning")
|
||||
request.setParameters('phensig')
|
||||
times = DataAccessLayer.getAvailableTimes(request)
|
||||
|
||||
# Get records for last 50 available times
|
||||
response = DataAccessLayer.getGeometryData(request, times[-50:-1])
|
||||
print("Using " + str(len(response)) + " records")
|
||||
|
||||
# Each record will have a numpy array the length of the number of "parameters"
|
||||
# Default is 1 (request.setParameters('phensig'))
|
||||
parameters = {}
|
||||
for x in request.getParameters():
|
||||
parameters[x] = np.array([])
|
||||
print(parameters)
|
||||
|
||||
|
||||
.. parsed-literal::
|
||||
|
||||
Using 109 records
|
||||
{'phensig': array([], dtype=float64)}
|
||||
|
||||
|
||||
Now loop through each record and plot it as either Polygon or
|
||||
MultiPolygon, with appropriate colors
|
||||
|
||||
.. code:: ipython3
|
||||
|
||||
%matplotlib inline
|
||||
bbox=[-127,-64,24,49]
|
||||
fig, ax = make_map(bbox=bbox)
|
||||
|
||||
siteids=np.array([])
|
||||
periods=np.array([])
|
||||
reftimes=np.array([])
|
||||
|
||||
for ob in response:
|
||||
|
||||
poly = ob.getGeometry()
|
||||
site = ob.getLocationName()
|
||||
pd = ob.getDataTime().getValidPeriod()
|
||||
ref = ob.getDataTime().getRefTime()
|
||||
|
||||
# do not plot if phensig is blank (SPS)
|
||||
if ob.getString('phensig'):
|
||||
|
||||
phensigString = ob.getString('phensig')
|
||||
|
||||
siteids = np.append(siteids,site)
|
||||
periods = np.append(periods,pd)
|
||||
reftimes = np.append(reftimes,ref)
|
||||
|
||||
for parm in parameters:
|
||||
parameters[parm] = np.append(parameters[parm],ob.getString(parm))
|
||||
|
||||
if poly.geom_type == 'MultiPolygon':
|
||||
geometries = np.array([])
|
||||
geometries = np.append(geometries,MultiPolygon(poly))
|
||||
geom_count = ", " + str(len(geometries)) +" geometries"
|
||||
else:
|
||||
geometries = np.array([])
|
||||
geometries = np.append(geometries,Polygon(poly))
|
||||
geom_count=""
|
||||
|
||||
for geom in geometries:
|
||||
bounds = Polygon(geom)
|
||||
intersection = bounds.intersection
|
||||
geoms = (intersection(geom)
|
||||
for geom in geometries
|
||||
if bounds.intersects(geom))
|
||||
|
||||
#print(vtec[phensigString]['hdln']
|
||||
# + " (" + phensigString + ") issued at " + str(ref)
|
||||
# + " ("+str(poly.geom_type) + geom_count + ")")
|
||||
|
||||
color = warning_color(phensigString)
|
||||
shape_feature = ShapelyFeature(geoms,ccrs.PlateCarree(),
|
||||
facecolor=color, edgecolor=color)
|
||||
ax.add_feature(shape_feature)
|
||||
|
||||
states_provinces = cfeature.NaturalEarthFeature(
|
||||
category='cultural',
|
||||
name='admin_1_states_provinces_lines',
|
||||
scale='50m',
|
||||
facecolor='none')
|
||||
political_boundaries = cfeature.NaturalEarthFeature(category='cultural',
|
||||
name='admin_0_boundary_lines_land',
|
||||
scale='50m', facecolor='none')
|
||||
ax.add_feature(cfeature.LAND)
|
||||
ax.add_feature(cfeature.COASTLINE)
|
||||
ax.add_feature(states_provinces, edgecolor='black')
|
||||
ax.add_feature(political_boundaries, edgecolor='black')
|
||||
|
||||
plt.show()
|
||||
|
||||
|
||||
|
||||
.. image:: Watch_and_Warning_Polygons_files/Watch_and_Warning_Polygons_5_0.png
|
||||
|
||||
|
|
@ -154,9 +154,7 @@ var Documentation = {
|
|||
this.fixFirefoxAnchorBug();
|
||||
this.highlightSearchWords();
|
||||
this.initIndexTable();
|
||||
if (DOCUMENTATION_OPTIONS.NAVIGATION_WITH_KEYS) {
|
||||
this.initOnKeyListeners();
|
||||
}
|
||||
this.initOnKeyListeners();
|
||||
},
|
||||
|
||||
/**
|
||||
|
|
@ -269,6 +267,13 @@ var Documentation = {
|
|||
window.history.replaceState({}, '', url);
|
||||
},
|
||||
|
||||
/**
|
||||
* helper function to focus on search bar
|
||||
*/
|
||||
focusSearchBar : function() {
|
||||
$('input[name=q]').first().focus();
|
||||
},
|
||||
|
||||
/**
|
||||
* make the url absolute
|
||||
*/
|
||||
|
|
@ -291,27 +296,54 @@ var Documentation = {
|
|||
},
|
||||
|
||||
initOnKeyListeners: function() {
|
||||
// only install a listener if it is really needed
|
||||
if (!DOCUMENTATION_OPTIONS.NAVIGATION_WITH_KEYS &&
|
||||
!DOCUMENTATION_OPTIONS.ENABLE_SEARCH_SHORTCUTS)
|
||||
return;
|
||||
|
||||
$(document).keydown(function(event) {
|
||||
var activeElementType = document.activeElement.tagName;
|
||||
// don't navigate when in search box, textarea, dropdown or button
|
||||
if (activeElementType !== 'TEXTAREA' && activeElementType !== 'INPUT' && activeElementType !== 'SELECT'
|
||||
&& activeElementType !== 'BUTTON' && !event.altKey && !event.ctrlKey && !event.metaKey
|
||||
&& !event.shiftKey) {
|
||||
switch (event.keyCode) {
|
||||
case 37: // left
|
||||
var prevHref = $('link[rel="prev"]').prop('href');
|
||||
if (prevHref) {
|
||||
window.location.href = prevHref;
|
||||
return false;
|
||||
}
|
||||
break;
|
||||
case 39: // right
|
||||
var nextHref = $('link[rel="next"]').prop('href');
|
||||
if (nextHref) {
|
||||
window.location.href = nextHref;
|
||||
return false;
|
||||
}
|
||||
break;
|
||||
&& activeElementType !== 'BUTTON') {
|
||||
if (event.altKey || event.ctrlKey || event.metaKey)
|
||||
return;
|
||||
|
||||
if (!event.shiftKey) {
|
||||
switch (event.key) {
|
||||
case 'ArrowLeft':
|
||||
if (!DOCUMENTATION_OPTIONS.NAVIGATION_WITH_KEYS)
|
||||
break;
|
||||
var prevHref = $('link[rel="prev"]').prop('href');
|
||||
if (prevHref) {
|
||||
window.location.href = prevHref;
|
||||
return false;
|
||||
}
|
||||
break;
|
||||
case 'ArrowRight':
|
||||
if (!DOCUMENTATION_OPTIONS.NAVIGATION_WITH_KEYS)
|
||||
break;
|
||||
var nextHref = $('link[rel="next"]').prop('href');
|
||||
if (nextHref) {
|
||||
window.location.href = nextHref;
|
||||
return false;
|
||||
}
|
||||
break;
|
||||
case 'Escape':
|
||||
if (!DOCUMENTATION_OPTIONS.ENABLE_SEARCH_SHORTCUTS)
|
||||
break;
|
||||
Documentation.hideSearchWords();
|
||||
return false;
|
||||
}
|
||||
}
|
||||
|
||||
// some keyboard layouts may need Shift to get /
|
||||
switch (event.key) {
|
||||
case '/':
|
||||
if (!DOCUMENTATION_OPTIONS.ENABLE_SEARCH_SHORTCUTS)
|
||||
break;
|
||||
Documentation.focusSearchBar();
|
||||
return false;
|
||||
}
|
||||
}
|
||||
});
|
||||
|
|
|
|||
|
|
@ -8,5 +8,7 @@ var DOCUMENTATION_OPTIONS = {
|
|||
LINK_SUFFIX: '.html',
|
||||
HAS_SOURCE: true,
|
||||
SOURCELINK_SUFFIX: '.txt',
|
||||
NAVIGATION_WITH_KEYS: false
|
||||
NAVIGATION_WITH_KEYS: false,
|
||||
SHOW_SEARCH_SUMMARY: true,
|
||||
ENABLE_SEARCH_SHORTCUTS: true,
|
||||
};
|
||||
|
|
@ -172,10 +172,6 @@ var Search = {
|
|||
}
|
||||
// stem the word
|
||||
var word = stemmer.stemWord(tmp[i].toLowerCase());
|
||||
// prevent stemmer from cutting word smaller than two chars
|
||||
if(word.length < 3 && tmp[i].length >= 3) {
|
||||
word = tmp[i];
|
||||
}
|
||||
var toAppend;
|
||||
// select the correct list
|
||||
if (word[0] == '-') {
|
||||
|
|
@ -276,7 +272,7 @@ var Search = {
|
|||
setTimeout(function() {
|
||||
displayNextItem();
|
||||
}, 5);
|
||||
} else if (DOCUMENTATION_OPTIONS.HAS_SOURCE) {
|
||||
} else if (DOCUMENTATION_OPTIONS.SHOW_SEARCH_SUMMARY) {
|
||||
$.ajax({url: requestUrl,
|
||||
dataType: "text",
|
||||
complete: function(jqxhr, textstatus) {
|
||||
|
|
@ -293,7 +289,7 @@ var Search = {
|
|||
}, 5);
|
||||
}});
|
||||
} else {
|
||||
// no source available, just display title
|
||||
// just display title
|
||||
Search.output.append(listItem);
|
||||
setTimeout(function() {
|
||||
displayNextItem();
|
||||
|
|
|
|||
4
dev.html
4
dev.html
|
|
@ -20,7 +20,7 @@
|
|||
<link rel="index" title="Index" href="genindex.html" />
|
||||
<link rel="search" title="Search" href="search.html" />
|
||||
<link rel="next" title="About Unidata AWIPS" href="about.html" />
|
||||
<link rel="prev" title="Watch and Warning Polygons" href="examples/generated/Watch_and_Warning_Polygons.html" />
|
||||
<link rel="prev" title="Watch Warning and Advisory Plotting" href="examples/generated/Watch_Warning_and_Advisory_Plotting.html" />
|
||||
</head>
|
||||
|
||||
<body class="wy-body-for-nav">
|
||||
|
|
@ -640,7 +640,7 @@ all factories may support this.</p></li>
|
|||
</div>
|
||||
</div>
|
||||
<footer><div class="rst-footer-buttons" role="navigation" aria-label="Footer">
|
||||
<a href="examples/generated/Watch_and_Warning_Polygons.html" class="btn btn-neutral float-left" title="Watch and Warning Polygons" accesskey="p" rel="prev"><span class="fa fa-arrow-circle-left" aria-hidden="true"></span> Previous</a>
|
||||
<a href="examples/generated/Watch_Warning_and_Advisory_Plotting.html" class="btn btn-neutral float-left" title="Watch Warning and Advisory Plotting" accesskey="p" rel="prev"><span class="fa fa-arrow-circle-left" aria-hidden="true"></span> Previous</a>
|
||||
<a href="about.html" class="btn btn-neutral float-right" title="About Unidata AWIPS" accesskey="n" rel="next">Next <span class="fa fa-arrow-circle-right" aria-hidden="true"></span></a>
|
||||
</div>
|
||||
|
||||
|
|
|
|||
|
|
@ -72,7 +72,7 @@
|
|||
<li class="toctree-l2"><a class="reference internal" href="Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
|
|
|
|||
|
|
@ -71,7 +71,7 @@
|
|||
<li class="toctree-l2"><a class="reference internal" href="Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
|
|
|
|||
|
|
@ -63,7 +63,7 @@
|
|||
<li class="toctree-l2"><a class="reference internal" href="Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
|
|
|
|||
|
|
@ -71,7 +71,7 @@
|
|||
<li class="toctree-l2"><a class="reference internal" href="Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
|
|
|
|||
|
|
@ -62,7 +62,7 @@
|
|||
<li class="toctree-l2"><a class="reference internal" href="Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
|
|
|
|||
|
|
@ -75,7 +75,7 @@
|
|||
<li class="toctree-l2"><a class="reference internal" href="Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
|
|
|
|||
|
|
@ -62,7 +62,7 @@
|
|||
<li class="toctree-l2"><a class="reference internal" href="Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
|
|
|
|||
|
|
@ -76,7 +76,7 @@
|
|||
<li class="toctree-l2"><a class="reference internal" href="Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
|
|
|
|||
|
|
@ -65,7 +65,7 @@
|
|||
<li class="toctree-l2"><a class="reference internal" href="Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
|
|
|
|||
|
|
@ -59,7 +59,7 @@
|
|||
<li class="toctree-l2"><a class="reference internal" href="Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
|
|
|
|||
|
|
@ -59,7 +59,7 @@
|
|||
<li class="toctree-l2"><a class="reference internal" href="Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
|
|
|
|||
|
|
@ -64,7 +64,7 @@
|
|||
</li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
|
|
|
|||
|
|
@ -65,7 +65,7 @@
|
|||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
|
|
|
|||
|
|
@ -19,7 +19,7 @@
|
|||
<link rel="author" title="About these documents" href="../../about.html" />
|
||||
<link rel="index" title="Index" href="../../genindex.html" />
|
||||
<link rel="search" title="Search" href="../../search.html" />
|
||||
<link rel="next" title="Watch and Warning Polygons" href="Watch_and_Warning_Polygons.html" />
|
||||
<link rel="next" title="Watch Warning and Advisory Plotting" href="Watch_Warning_and_Advisory_Plotting.html" />
|
||||
<link rel="prev" title="Satellite Imagery" href="Satellite_Imagery.html" />
|
||||
</head>
|
||||
|
||||
|
|
@ -59,7 +59,7 @@
|
|||
<li class="toctree-l2"><a class="reference internal" href="Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2 current"><a class="current reference internal" href="#">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
|
|
@ -240,7 +240,7 @@ used to plot the wind profile.</p>
|
|||
</div>
|
||||
<footer><div class="rst-footer-buttons" role="navigation" aria-label="Footer">
|
||||
<a href="Satellite_Imagery.html" class="btn btn-neutral float-left" title="Satellite Imagery" accesskey="p" rel="prev"><span class="fa fa-arrow-circle-left" aria-hidden="true"></span> Previous</a>
|
||||
<a href="Watch_and_Warning_Polygons.html" class="btn btn-neutral float-right" title="Watch and Warning Polygons" accesskey="n" rel="next">Next <span class="fa fa-arrow-circle-right" aria-hidden="true"></span></a>
|
||||
<a href="Watch_Warning_and_Advisory_Plotting.html" class="btn btn-neutral float-right" title="Watch Warning and Advisory Plotting" accesskey="n" rel="next">Next <span class="fa fa-arrow-circle-right" aria-hidden="true"></span></a>
|
||||
</div>
|
||||
|
||||
<hr/>
|
||||
|
|
|
|||
657
examples/generated/Watch_Warning_and_Advisory_Plotting.html
Normal file
657
examples/generated/Watch_Warning_and_Advisory_Plotting.html
Normal file
|
|
@ -0,0 +1,657 @@
|
|||
<!DOCTYPE html>
|
||||
<html class="writer-html5" lang="en" >
|
||||
<head>
|
||||
<meta charset="utf-8" /><meta name="generator" content="Docutils 0.17.1: http://docutils.sourceforge.net/" />
|
||||
|
||||
<meta name="viewport" content="width=device-width, initial-scale=1.0" />
|
||||
<title>Watch Warning and Advisory Plotting — python-awips documentation</title>
|
||||
<link rel="stylesheet" href="../../_static/pygments.css" type="text/css" />
|
||||
<link rel="stylesheet" href="../../_static/css/theme.css" type="text/css" />
|
||||
<!--[if lt IE 9]>
|
||||
<script src="../../_static/js/html5shiv.min.js"></script>
|
||||
<![endif]-->
|
||||
|
||||
<script data-url_root="../../" id="documentation_options" src="../../_static/documentation_options.js"></script>
|
||||
<script src="../../_static/jquery.js"></script>
|
||||
<script src="../../_static/underscore.js"></script>
|
||||
<script src="../../_static/doctools.js"></script>
|
||||
<script src="../../_static/js/theme.js"></script>
|
||||
<link rel="author" title="About these documents" href="../../about.html" />
|
||||
<link rel="index" title="Index" href="../../genindex.html" />
|
||||
<link rel="search" title="Search" href="../../search.html" />
|
||||
<link rel="next" title="Development Guide" href="../../dev.html" />
|
||||
<link rel="prev" title="Upper Air BUFR Soundings" href="Upper_Air_BUFR_Soundings.html" />
|
||||
</head>
|
||||
|
||||
<body class="wy-body-for-nav">
|
||||
<div class="wy-grid-for-nav">
|
||||
<nav data-toggle="wy-nav-shift" class="wy-nav-side">
|
||||
<div class="wy-side-scroll">
|
||||
<div class="wy-side-nav-search" >
|
||||
<a href="../../index.html" class="icon icon-home"> python-awips
|
||||
</a>
|
||||
<div class="version">
|
||||
18.1.8
|
||||
</div>
|
||||
<div role="search">
|
||||
<form id="rtd-search-form" class="wy-form" action="../../search.html" method="get">
|
||||
<input type="text" name="q" placeholder="Search docs" />
|
||||
<input type="hidden" name="check_keywords" value="yes" />
|
||||
<input type="hidden" name="area" value="default" />
|
||||
</form>
|
||||
</div>
|
||||
</div><div class="wy-menu wy-menu-vertical" data-spy="affix" role="navigation" aria-label="Navigation menu">
|
||||
<ul class="current">
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../api/index.html">API Documentation</a></li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../datatypes.html">Available Data Types</a></li>
|
||||
<li class="toctree-l1 current"><a class="reference internal" href="../index.html">Data Plotting Examples</a><ul class="current">
|
||||
<li class="toctree-l2"><a class="reference internal" href="Colored_Surface_Temperature_Plot.html">Colored Surface Temperature Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Colorized_Grid_Data.html">Colorized Grid Data</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Forecast_Model_Vertical_Sounding.html">Forecast Model Vertical Sounding</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="GOES_CIRA_Product_Writer.html">GOES CIRA Product Writer</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="GOES_Geostationary_Lightning_Mapper.html">GOES Geostationary Lightning Mapper</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Grid_Levels_and_Parameters.html">Grid Levels and Parameters</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="METAR_Station_Plot_with_MetPy.html">METAR Station Plot with MetPy</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Map_Resources_and_Topography.html">Map Resources and Topography</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Model_Sounding_Data.html">Model Sounding Data</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="NEXRAD_Level3_Radar.html">NEXRAD Level3 Radar</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Precip_Accumulation-Region_Of_Interest.html">Precip Accumulation-Region Of Interest</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2 current"><a class="current reference internal" href="#">Watch Warning and Advisory Plotting</a><ul>
|
||||
<li class="toctree-l3"><a class="reference internal" href="#objectives">Objectives</a><ul>
|
||||
<li class="toctree-l4"><a class="reference internal" href="#table-of-contents">Table of Contents</a></li>
|
||||
<li class="toctree-l4"><a class="reference internal" href="#imports">1 Imports</a></li>
|
||||
<li class="toctree-l4"><a class="reference internal" href="#function-make-map">2 Function: make_map()</a></li>
|
||||
<li class="toctree-l4"><a class="reference internal" href="#function-get-color">3 Function: get_color()</a></li>
|
||||
<li class="toctree-l4"><a class="reference internal" href="#function-get-title">4 Function get_title()</a></li>
|
||||
<li class="toctree-l4"><a class="reference internal" href="#inital-setup">5 Inital Setup</a></li>
|
||||
<li class="toctree-l4"><a class="reference internal" href="#filter-by-time">6 Filter by Time</a></li>
|
||||
<li class="toctree-l4"><a class="reference internal" href="#use-the-data">7 Use the Data!</a></li>
|
||||
<li class="toctree-l4"><a class="reference internal" href="#plot-the-data">8 Plot the Data!</a></li>
|
||||
<li class="toctree-l4"><a class="reference internal" href="#see-also">9 See Also</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
</ul>
|
||||
</li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
<li class="toctree-l1"><a class="reference external" href="http://unidata.github.io/awips2/appendix/appendix-grid-parameters/">AWIPS Grid Parameters</a></li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../about.html">About Unidata AWIPS</a></li>
|
||||
</ul>
|
||||
|
||||
</div>
|
||||
</div>
|
||||
</nav>
|
||||
|
||||
<section data-toggle="wy-nav-shift" class="wy-nav-content-wrap"><nav class="wy-nav-top" aria-label="Mobile navigation menu" >
|
||||
<i data-toggle="wy-nav-top" class="fa fa-bars"></i>
|
||||
<a href="../../index.html">python-awips</a>
|
||||
</nav>
|
||||
|
||||
<div class="wy-nav-content">
|
||||
<div class="rst-content">
|
||||
<div role="navigation" aria-label="Page navigation">
|
||||
<ul class="wy-breadcrumbs">
|
||||
<li><a href="../../index.html" class="icon icon-home"></a> »</li>
|
||||
<li><a href="../index.html">Data Plotting Examples</a> »</li>
|
||||
<li>Watch Warning and Advisory Plotting</li>
|
||||
<li class="wy-breadcrumbs-aside">
|
||||
<a href="../../_sources/examples/generated/Watch_Warning_and_Advisory_Plotting.rst.txt" rel="nofollow"> View page source</a>
|
||||
</li>
|
||||
</ul>
|
||||
<hr/>
|
||||
</div>
|
||||
<div role="main" class="document" itemscope="itemscope" itemtype="http://schema.org/Article">
|
||||
<div itemprop="articleBody">
|
||||
|
||||
<section id="watch-warning-and-advisory-plotting">
|
||||
<h1>Watch Warning and Advisory Plotting<a class="headerlink" href="#watch-warning-and-advisory-plotting" title="Permalink to this headline"></a></h1>
|
||||
<p><a class="reference external" href="http://nbviewer.ipython.org/github/Unidata/python-awips/blob/master/examples/notebooks/Watch_Warning_and_Advisory_Plotting.ipynb">Notebook</a>
|
||||
Python-AWIPS Tutorial Notebook</p>
|
||||
<hr class="docutils" />
|
||||
<section id="objectives">
|
||||
<h2>Objectives<a class="headerlink" href="#objectives" title="Permalink to this headline"></a></h2>
|
||||
<ul class="simple">
|
||||
<li><p>Create a colorized plot with <a class="reference external" href="https://weather.cod.edu/notes/criteria/">Warnings, Watches, Advisories and
|
||||
Statements (WWAs)</a></p></li>
|
||||
<li><p>Use python-awips to connect to an EDEX server</p></li>
|
||||
<li><p>Create and filter the data request specifically for a warning data
|
||||
type</p></li>
|
||||
<li><p>Create and use accurate time filter for data requests</p></li>
|
||||
<li><p>Define and use functions</p></li>
|
||||
<li><p>Define and use dictionaries</p></li>
|
||||
<li><p>Colorize shapes based on a dictionary</p></li>
|
||||
<li><p>Overlay warnings, watches, and advisories with state and political
|
||||
maps</p></li>
|
||||
</ul>
|
||||
<hr class="docutils" />
|
||||
<section id="table-of-contents">
|
||||
<h3>Table of Contents<a class="headerlink" href="#table-of-contents" title="Permalink to this headline"></a></h3>
|
||||
<div class="line-block">
|
||||
<div class="line"><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#imports">1
|
||||
Imports</a></div>
|
||||
<div class="line"><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#functions-make-map">2 Functions:
|
||||
make_map()</a></div>
|
||||
<div class="line"><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#functions-get-color">3 Functions:
|
||||
get_color()</a></div>
|
||||
<div class="line"><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#functions-get-title">4 Functions:
|
||||
get_title()</a></div>
|
||||
<div class="line"><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#initial-setup">5 Initial
|
||||
Setup</a></div>
|
||||
<div class="line"> <a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#edex-connection">5.1 EDEX
|
||||
Connection</a></div>
|
||||
<div class="line"> <a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#significance-sig-constraints">5.2 Significance (Sig)
|
||||
Constraints</a></div>
|
||||
<div class="line"><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#filter-by-time">6 Filter by
|
||||
Time</a></div>
|
||||
<div class="line"><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#use-the-data">7 Use the
|
||||
Data!</a></div>
|
||||
<div class="line"> <a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#get-the-data">7.1 Get the
|
||||
Data</a></div>
|
||||
<div class="line"> <a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#extract-phensigs-geometries-and-times">7.2 Extract Phensigs, Geometries, and
|
||||
Times</a></div>
|
||||
<div class="line"><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#plot-the-data">8 Plot the
|
||||
Data!</a></div>
|
||||
<div class="line"> <a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#create-state-and-political-boundaries">8.1 Create State and Political
|
||||
Boundaries</a></div>
|
||||
<div class="line"> <a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#draw-the-plot-and-legend-for-wwas">8.2 Draw the Plot and Legend for
|
||||
WWAs</a></div>
|
||||
<div class="line"><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#see-also">9 See
|
||||
Also</a></div>
|
||||
<div class="line"> <a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#related-notebooks">9.1 Related
|
||||
Notebooks</a></div>
|
||||
<div class="line"> <a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html#additional-documentation">9.2 Additional
|
||||
Documentation</a></div>
|
||||
</div>
|
||||
</section>
|
||||
<section id="imports">
|
||||
<h3>1 Imports<a class="headerlink" href="#imports" title="Permalink to this headline"></a></h3>
|
||||
<p>The imports below are used throughout the notebook. The python-awips
|
||||
imports allow us to connect to an EDEX server, use the warning lookup
|
||||
dictionary, and define a TimeRange. The additional imports are for data
|
||||
manipulation and visualization.</p>
|
||||
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="kn">from</span> <span class="nn">datetime</span> <span class="kn">import</span> <span class="n">datetime</span><span class="p">,</span> <span class="n">timedelta</span>
|
||||
|
||||
<span class="kn">import</span> <span class="nn">numpy</span> <span class="k">as</span> <span class="nn">np</span>
|
||||
<span class="kn">import</span> <span class="nn">matplotlib.patches</span> <span class="k">as</span> <span class="nn">mpatches</span>
|
||||
<span class="kn">import</span> <span class="nn">matplotlib.pyplot</span> <span class="k">as</span> <span class="nn">plt</span>
|
||||
<span class="kn">import</span> <span class="nn">cartopy.crs</span> <span class="k">as</span> <span class="nn">ccrs</span>
|
||||
<span class="kn">import</span> <span class="nn">cartopy.feature</span> <span class="k">as</span> <span class="nn">cfeature</span>
|
||||
<span class="kn">from</span> <span class="nn">cartopy.feature</span> <span class="kn">import</span> <span class="n">ShapelyFeature</span><span class="p">,</span> <span class="n">NaturalEarthFeature</span>
|
||||
<span class="kn">from</span> <span class="nn">cartopy.mpl.gridliner</span> <span class="kn">import</span> <span class="n">LONGITUDE_FORMATTER</span><span class="p">,</span> <span class="n">LATITUDE_FORMATTER</span>
|
||||
<span class="kn">from</span> <span class="nn">shapely.geometry</span> <span class="kn">import</span> <span class="n">MultiPolygon</span><span class="p">,</span> <span class="n">Polygon</span>
|
||||
|
||||
<span class="kn">from</span> <span class="nn">awips.dataaccess</span> <span class="kn">import</span> <span class="n">DataAccessLayer</span>
|
||||
<span class="kn">from</span> <span class="nn">awips.tables</span> <span class="kn">import</span> <span class="n">vtec</span>
|
||||
<span class="kn">from</span> <span class="nn">dynamicserialize.dstypes.com.raytheon.uf.common.time</span> <span class="kn">import</span> <span class="n">TimeRange</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
<p><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html">Top</a></p>
|
||||
</section>
|
||||
<hr class="docutils" />
|
||||
<section id="function-make-map">
|
||||
<h3>2 Function: make_map()<a class="headerlink" href="#function-make-map" title="Permalink to this headline"></a></h3>
|
||||
<p>In order to plot more than one image, it’s easiest to define common
|
||||
logic in a function. However, for this notebook we only use it in one
|
||||
place. It is a function you will find in most of our example notebooks.</p>
|
||||
<p>Here, a new function called <strong>make_map</strong> is defined. This function uses
|
||||
the <a class="reference external" href="https://matplotlib.org/3.3.3/api/_as_gen/matplotlib.pyplot.html">matplotlib.pyplot package
|
||||
(plt)</a>
|
||||
to create a figure and axis. The lat/lon grids are added.</p>
|
||||
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="k">def</span> <span class="nf">make_map</span><span class="p">(</span><span class="n">bbox</span><span class="p">,</span> <span class="n">projection</span><span class="o">=</span><span class="n">ccrs</span><span class="o">.</span><span class="n">PlateCarree</span><span class="p">()):</span>
|
||||
<span class="n">fig</span><span class="p">,</span> <span class="n">ax</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">subplots</span><span class="p">(</span><span class="n">figsize</span><span class="o">=</span><span class="p">(</span><span class="mi">20</span><span class="p">,</span><span class="mi">12</span><span class="p">),</span>
|
||||
<span class="n">subplot_kw</span><span class="o">=</span><span class="nb">dict</span><span class="p">(</span><span class="n">projection</span><span class="o">=</span><span class="n">projection</span><span class="p">))</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">set_extent</span><span class="p">(</span><span class="n">bbox</span><span class="p">)</span>
|
||||
<span class="n">gl</span> <span class="o">=</span> <span class="n">ax</span><span class="o">.</span><span class="n">gridlines</span><span class="p">(</span><span class="n">draw_labels</span><span class="o">=</span><span class="kc">True</span><span class="p">)</span>
|
||||
<span class="n">gl</span><span class="o">.</span><span class="n">top_labels</span> <span class="o">=</span> <span class="n">gl</span><span class="o">.</span><span class="n">right_labels</span> <span class="o">=</span> <span class="kc">False</span>
|
||||
<span class="n">gl</span><span class="o">.</span><span class="n">xformatter</span> <span class="o">=</span> <span class="n">LONGITUDE_FORMATTER</span>
|
||||
<span class="n">gl</span><span class="o">.</span><span class="n">yformatter</span> <span class="o">=</span> <span class="n">LATITUDE_FORMATTER</span>
|
||||
<span class="k">return</span> <span class="n">fig</span><span class="p">,</span> <span class="n">ax</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
<p><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html">Top</a></p>
|
||||
</section>
|
||||
<hr class="docutils" />
|
||||
<section id="function-get-color">
|
||||
<h3>3 Function: get_color()<a class="headerlink" href="#function-get-color" title="Permalink to this headline"></a></h3>
|
||||
<p>Since we’ll be needing to access the color using the vtec lookup table
|
||||
in several places, creating an easily recognizable function is useful.</p>
|
||||
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="k">def</span> <span class="nf">get_color</span><span class="p">(</span><span class="n">phensig</span><span class="p">):</span>
|
||||
<span class="k">return</span> <span class="n">vtec</span><span class="p">[</span><span class="n">phensig</span><span class="p">][</span><span class="s1">'color'</span><span class="p">]</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
<p><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html">Top</a></p>
|
||||
</section>
|
||||
<hr class="docutils" />
|
||||
<section id="function-get-title">
|
||||
<h3>4 Function get_title()<a class="headerlink" href="#function-get-title" title="Permalink to this headline"></a></h3>
|
||||
<p>Similar to the color function just defined, accessing the full name for
|
||||
the phensig will also be necessary, so this function will be helpful.</p>
|
||||
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="k">def</span> <span class="nf">get_title</span><span class="p">(</span><span class="n">phensig</span><span class="p">):</span>
|
||||
<span class="k">return</span> <span class="n">vtec</span><span class="p">[</span><span class="n">phensig</span><span class="p">][</span><span class="s1">'hdln'</span><span class="p">]</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
<p><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html">Top</a></p>
|
||||
</section>
|
||||
<hr class="docutils" />
|
||||
<section id="inital-setup">
|
||||
<h3>5 Inital Setup<a class="headerlink" href="#inital-setup" title="Permalink to this headline"></a></h3>
|
||||
<section id="edex-connection">
|
||||
<h4>5.1 EDEX Connection<a class="headerlink" href="#edex-connection" title="Permalink to this headline"></a></h4>
|
||||
<p>First we establish a connection to Unidata’s public EDEX server. With
|
||||
that connection made, we can create a <a class="reference external" href="http://unidata.github.io/python-awips/api/IDataRequest.html">new data request
|
||||
object</a>
|
||||
and set the data type to <strong>warning</strong>, and set the Parameters to
|
||||
<strong>phensig</strong> and <strong>sig</strong>.</p>
|
||||
<div class="alert-info docutils container">
|
||||
<p>Note: Remember, to see all available parameters use the
|
||||
DataAccess.getAvailableParameters() method as shown in the Grid
|
||||
Levels and Parameters Notebook.</p>
|
||||
</div>
|
||||
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="n">DataAccessLayer</span><span class="o">.</span><span class="n">changeEDEXHost</span><span class="p">(</span><span class="s2">"edex-cloud.unidata.ucar.edu"</span><span class="p">)</span>
|
||||
<span class="n">request</span> <span class="o">=</span> <span class="n">DataAccessLayer</span><span class="o">.</span><span class="n">newDataRequest</span><span class="p">()</span>
|
||||
<span class="n">request</span><span class="o">.</span><span class="n">setDatatype</span><span class="p">(</span><span class="s2">"warning"</span><span class="p">)</span>
|
||||
<span class="n">params</span> <span class="o">=</span> <span class="p">[</span><span class="s2">"phensig"</span><span class="p">,</span> <span class="s2">"sig"</span><span class="p">]</span>
|
||||
<span class="n">request</span><span class="o">.</span><span class="n">setParameters</span><span class="p">(</span><span class="o">*</span><span class="p">(</span><span class="n">params</span><span class="p">))</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
</section>
|
||||
<section id="significance-sig-constants">
|
||||
<h4>5.2 Significance (Sig) Constants<a class="headerlink" href="#significance-sig-constants" title="Permalink to this headline"></a></h4>
|
||||
<p>The two parameters we’re requesting for our warning objects are
|
||||
<strong>phensig</strong> and <strong>sig</strong> where phensig is styled “XX.Y” and sig is “Y”.
|
||||
Phen stands for “Phenomena” and sig stands for “Significance”. <a class="reference external" href="https://www.weather.gov/media/vtec/VTEC_explanation4-20.pdf">A more
|
||||
detailed description of phensigs and how they’re used is provided with
|
||||
this NWS
|
||||
pamphlet</a>.</p>
|
||||
<p>The constants in this section correlate the <strong>sig</strong> to what type of
|
||||
message it is (what significance it is).</p>
|
||||
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="n">WATCH_SIG</span> <span class="o">=</span> <span class="s1">'A'</span>
|
||||
<span class="n">WARN_SIG</span> <span class="o">=</span> <span class="s1">'W'</span>
|
||||
<span class="n">ADVIS_SIG</span> <span class="o">=</span> <span class="s1">'Y'</span>
|
||||
<span class="n">STATEM_SIG</span> <span class="o">=</span> <span class="s1">'S'</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
<p><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html">Top</a></p>
|
||||
</section>
|
||||
</section>
|
||||
<hr class="docutils" />
|
||||
<section id="filter-by-time">
|
||||
<h3>6 Filter by Time<a class="headerlink" href="#filter-by-time" title="Permalink to this headline"></a></h3>
|
||||
<p>Here we decide how much data we want to pull from EDEX. By default we’ll
|
||||
request 12 hours, but that value can easily be modified by <cite>adjusting
|
||||
the
|
||||
``timedelta(hours = 12)`</cite> <<a class="reference external" href="https://docs.python.org/3/library/datetime.html#timedelta-objects">https://docs.python.org/3/library/datetime.html#timedelta-objects</a>>`__
|
||||
in line <code class="docutils literal notranslate"><span class="pre">2</span></code>. The more data we request, the longer the next section
|
||||
will take to run.</p>
|
||||
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="c1"># Get records from the last 12 hours</span>
|
||||
<span class="n">lastHourDateTime</span> <span class="o">=</span> <span class="n">datetime</span><span class="o">.</span><span class="n">utcnow</span><span class="p">()</span> <span class="o">-</span> <span class="n">timedelta</span><span class="p">(</span><span class="n">hours</span> <span class="o">=</span> <span class="mi">12</span><span class="p">)</span>
|
||||
<span class="n">start</span> <span class="o">=</span> <span class="n">lastHourDateTime</span><span class="o">.</span><span class="n">strftime</span><span class="p">(</span><span class="s1">'%Y-%m-</span><span class="si">%d</span><span class="s1"> %H:%M:%S'</span><span class="p">)</span>
|
||||
<span class="n">end</span> <span class="o">=</span> <span class="n">datetime</span><span class="o">.</span><span class="n">utcnow</span><span class="p">()</span><span class="o">.</span><span class="n">strftime</span><span class="p">(</span><span class="s1">'%Y-%m-</span><span class="si">%d</span><span class="s1"> %H:%M:%S'</span><span class="p">)</span>
|
||||
|
||||
<span class="n">beginRange</span> <span class="o">=</span> <span class="n">datetime</span><span class="o">.</span><span class="n">strptime</span><span class="p">(</span> <span class="n">start</span> <span class="p">,</span> <span class="s2">"%Y-%m-</span><span class="si">%d</span><span class="s2"> %H:%M:%S"</span><span class="p">)</span>
|
||||
<span class="n">endRange</span> <span class="o">=</span> <span class="n">datetime</span><span class="o">.</span><span class="n">strptime</span><span class="p">(</span> <span class="n">end</span> <span class="p">,</span> <span class="s2">"%Y-%m-</span><span class="si">%d</span><span class="s2"> %H:%M:%S"</span><span class="p">)</span>
|
||||
<span class="n">timerange</span> <span class="o">=</span> <span class="n">TimeRange</span><span class="p">(</span><span class="n">beginRange</span><span class="p">,</span> <span class="n">endRange</span><span class="p">)</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
<p><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html">Top</a></p>
|
||||
</section>
|
||||
<hr class="docutils" />
|
||||
<section id="use-the-data">
|
||||
<h3>7 Use the Data!<a class="headerlink" href="#use-the-data" title="Permalink to this headline"></a></h3>
|
||||
<section id="get-the-data">
|
||||
<h4>7.1 Get the Data<a class="headerlink" href="#get-the-data" title="Permalink to this headline"></a></h4>
|
||||
<p>Now that we have our <code class="docutils literal notranslate"><span class="pre">request</span></code> and TimeRange <code class="docutils literal notranslate"><span class="pre">timerange</span></code> objects
|
||||
ready, it’s time to request the data array from EDEX.</p>
|
||||
<div class="alert-info docutils container">
|
||||
<p>Note: Above we set timerange to be 12 hours worth of data. This can
|
||||
return on the order of ~2000 records and can take a little while to
|
||||
run.</p>
|
||||
</div>
|
||||
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="c1"># Get response</span>
|
||||
<span class="n">response</span> <span class="o">=</span> <span class="n">DataAccessLayer</span><span class="o">.</span><span class="n">getGeometryData</span><span class="p">(</span><span class="n">request</span><span class="p">,</span> <span class="n">timerange</span><span class="p">)</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="s2">"Using "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">response</span><span class="p">))</span> <span class="o">+</span> <span class="s2">" records"</span><span class="p">)</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-default notranslate"><div class="highlight"><pre><span></span><span class="n">Using</span> <span class="mi">1502</span> <span class="n">records</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
</section>
|
||||
<section id="extract-phensigs-geometries-and-times">
|
||||
<h4>7.2 Extract Phensigs, Geometries, and Times<a class="headerlink" href="#extract-phensigs-geometries-and-times" title="Permalink to this headline"></a></h4>
|
||||
<p>In this section we start gathering all the information we’ll need to
|
||||
properly display our data. First we create an array to keep track of
|
||||
unique phensigs. This is useful summary information and will be used to
|
||||
help create the legend which we’ll display along with our plot.</p>
|
||||
<p>Next, we create arrays for each of the 4 types of significance a
|
||||
statement can have. We will group our records this way, so we can easily
|
||||
toggle which records to display or not.</p>
|
||||
<p>Then, we create two time variables to keep track of the earliest time
|
||||
from our records and the latest time, and will display that information
|
||||
in the title of our plot.</p>
|
||||
<p>This section has optional print statements at lines <code class="docutils literal notranslate"><span class="pre">65</span></code> and <code class="docutils literal notranslate"><span class="pre">85</span></code>.
|
||||
The first prints out the title, phensig, ref time, and shape for each
|
||||
unique phensig, and the second prints out a sum of how many unique
|
||||
phensigs there are.</p>
|
||||
<p>We cycle through all the data produced from our <code class="docutils literal notranslate"><span class="pre">response</span></code> object,
|
||||
access its geometries, and create a new
|
||||
<a class="reference external" href="https://scitools.org.uk/cartopy/docs/latest/reference/generated/cartopy.feature.ShapelyFeature.html">ShapelyFeature</a>
|
||||
with the corresponding color. Then we place this new feature in the
|
||||
appropriate <code class="docutils literal notranslate"><span class="pre">shapes</span></code> array. During this process we also populate the
|
||||
phensigs array with all unique phensig entries.</p>
|
||||
<p>Finally, after we’re done looping through all the <code class="docutils literal notranslate"><span class="pre">response</span></code> data, we
|
||||
create a mapping of phensigs to their corresponding titles. This will be
|
||||
used later to sort the legend alphabetically by titles (which differs
|
||||
from simply sorting by phensig). Ex. <em>Blizzard Warning (BZ.W)</em> would
|
||||
come before <em>Areal Flood Advisory (FA.Y)</em> if we simply sorted by
|
||||
phensig.</p>
|
||||
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="c1"># Keep track of unique phensigs, to use in legend</span>
|
||||
<span class="n">phensigs</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
|
||||
<span class="c1"># Sort the geometries based on their sig</span>
|
||||
<span class="n">watch_shapes</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
<span class="n">warning_shapes</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
<span class="n">advisory_shapes</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
<span class="n">statement_shapes</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
|
||||
<span class="c1"># Keep track of the earliest and latest reftime for title</span>
|
||||
<span class="c1"># start with the first time from the first object in the response</span>
|
||||
<span class="n">time_str</span> <span class="o">=</span> <span class="nb">str</span><span class="p">(</span><span class="n">response</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span><span class="o">.</span><span class="n">getDataTime</span><span class="p">()</span><span class="o">.</span><span class="n">getRefTime</span><span class="p">())</span>
|
||||
<span class="c1"># truncate the decimal seconds for datetime parsing</span>
|
||||
<span class="n">time_str</span> <span class="o">=</span> <span class="n">time_str</span><span class="p">[:</span><span class="o">-</span><span class="mi">4</span><span class="p">]</span>
|
||||
<span class="c1"># convert to datetime object for easy comparison</span>
|
||||
<span class="n">first_time</span> <span class="o">=</span> <span class="n">datetime</span><span class="o">.</span><span class="n">strptime</span><span class="p">(</span><span class="n">time_str</span><span class="p">,</span> <span class="s1">'%Y-%m-</span><span class="si">%d</span><span class="s1"> %H:%M:%S'</span><span class="p">)</span>
|
||||
<span class="n">last_time</span> <span class="o">=</span> <span class="n">datetime</span><span class="o">.</span><span class="n">strptime</span><span class="p">(</span><span class="n">time_str</span><span class="p">,</span> <span class="s1">'%Y-%m-</span><span class="si">%d</span><span class="s1"> %H:%M:%S'</span><span class="p">)</span>
|
||||
|
||||
<span class="k">for</span> <span class="n">ob</span> <span class="ow">in</span> <span class="n">response</span><span class="p">:</span>
|
||||
|
||||
<span class="c1"># get the geometry for the object</span>
|
||||
<span class="n">poly</span> <span class="o">=</span> <span class="n">ob</span><span class="o">.</span><span class="n">getGeometry</span><span class="p">()</span>
|
||||
<span class="c1"># get the reftime for the object</span>
|
||||
<span class="n">ref</span> <span class="o">=</span> <span class="n">ob</span><span class="o">.</span><span class="n">getDataTime</span><span class="p">()</span><span class="o">.</span><span class="n">getRefTime</span><span class="p">()</span>
|
||||
|
||||
<span class="c1"># do not plot if phensig is blank (SPS)</span>
|
||||
<span class="k">if</span> <span class="n">ob</span><span class="o">.</span><span class="n">getString</span><span class="p">(</span><span class="s1">'phensig'</span><span class="p">):</span>
|
||||
|
||||
<span class="c1"># look at the reftime</span>
|
||||
<span class="c1"># convert reftime to a string and parse the decimal seconds</span>
|
||||
<span class="n">ref_str</span> <span class="o">=</span> <span class="nb">str</span><span class="p">(</span><span class="n">ref</span><span class="p">)</span>
|
||||
<span class="n">ref_str</span> <span class="o">=</span> <span class="n">ref_str</span><span class="p">[:</span><span class="o">-</span><span class="mi">4</span><span class="p">]</span>
|
||||
<span class="c1"># convert reftime to a datetime object for comparison</span>
|
||||
<span class="n">ref_time</span> <span class="o">=</span> <span class="n">datetime</span><span class="o">.</span><span class="n">strptime</span><span class="p">(</span><span class="n">ref_str</span><span class="p">,</span> <span class="s1">'%Y-%m-</span><span class="si">%d</span><span class="s1"> %H:%M:%S'</span><span class="p">)</span>
|
||||
<span class="c1"># compare current time with first and last times and set if appropriate</span>
|
||||
<span class="k">if</span> <span class="n">ref_time</span> <span class="ow">is</span> <span class="ow">not</span> <span class="kc">None</span><span class="p">:</span>
|
||||
<span class="k">if</span> <span class="n">ref_time</span> <span class="o"><</span> <span class="n">first_time</span><span class="p">:</span>
|
||||
<span class="n">first_time</span> <span class="o">=</span> <span class="n">ref_time</span>
|
||||
<span class="k">elif</span> <span class="n">ref_time</span> <span class="o">></span> <span class="n">last_time</span><span class="p">:</span>
|
||||
<span class="n">last_time</span> <span class="o">=</span> <span class="n">ref_time</span>
|
||||
|
||||
<span class="c1"># get the phensig and sig values from object</span>
|
||||
<span class="n">phensigString</span> <span class="o">=</span> <span class="n">ob</span><span class="o">.</span><span class="n">getString</span><span class="p">(</span><span class="s1">'phensig'</span><span class="p">)</span>
|
||||
<span class="n">sig</span> <span class="o">=</span> <span class="n">ob</span><span class="o">.</span><span class="n">getString</span><span class="p">(</span><span class="s1">'sig'</span><span class="p">)</span>
|
||||
|
||||
<span class="c1"># set the geometries based on whether it's a MultiPolygon or Polygon</span>
|
||||
<span class="k">if</span> <span class="n">poly</span><span class="o">.</span><span class="n">geom_type</span> <span class="o">==</span> <span class="s1">'MultiPolygon'</span><span class="p">:</span>
|
||||
<span class="n">geometries</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([])</span>
|
||||
<span class="n">geometries</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">geometries</span><span class="p">,</span><span class="n">MultiPolygon</span><span class="p">(</span><span class="n">poly</span><span class="p">))</span>
|
||||
<span class="k">else</span><span class="p">:</span>
|
||||
<span class="n">geometries</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([])</span>
|
||||
<span class="n">geometries</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">geometries</span><span class="p">,</span><span class="n">Polygon</span><span class="p">(</span><span class="n">poly</span><span class="p">))</span>
|
||||
|
||||
<span class="k">for</span> <span class="n">geom</span> <span class="ow">in</span> <span class="n">geometries</span><span class="p">:</span>
|
||||
<span class="n">bounds</span> <span class="o">=</span> <span class="n">Polygon</span><span class="p">(</span><span class="n">geom</span><span class="p">)</span>
|
||||
<span class="n">intersection</span> <span class="o">=</span> <span class="n">bounds</span><span class="o">.</span><span class="n">intersection</span>
|
||||
<span class="n">geoms</span> <span class="o">=</span> <span class="p">(</span><span class="n">intersection</span><span class="p">(</span><span class="n">geom</span><span class="p">)</span> <span class="k">for</span> <span class="n">geom</span> <span class="ow">in</span> <span class="n">geometries</span> <span class="k">if</span> <span class="n">bounds</span><span class="o">.</span><span class="n">intersects</span><span class="p">(</span><span class="n">geom</span><span class="p">))</span>
|
||||
|
||||
<span class="c1"># Store the unique phensigs</span>
|
||||
<span class="k">if</span> <span class="ow">not</span> <span class="n">phensigString</span> <span class="ow">in</span> <span class="n">phensigs</span><span class="p">:</span>
|
||||
<span class="n">phensigs</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">phensigString</span><span class="p">)</span>
|
||||
<span class="c1"># Optional printout of unique Phensigs</span>
|
||||
<span class="c1"># print(get_title(phensigString) + " (" + phensigString + ")</span>
|
||||
|
||||
<span class="c1"># get the corresponding color using the dictionary</span>
|
||||
<span class="n">color</span> <span class="o">=</span> <span class="n">get_color</span><span class="p">(</span><span class="n">phensigString</span><span class="p">)</span>
|
||||
<span class="c1"># create a new shape feature for the object</span>
|
||||
<span class="n">shape_feature</span> <span class="o">=</span> <span class="n">ShapelyFeature</span><span class="p">(</span><span class="n">geoms</span><span class="p">,</span><span class="n">ccrs</span><span class="o">.</span><span class="n">PlateCarree</span><span class="p">(),</span>
|
||||
<span class="n">facecolor</span><span class="o">=</span><span class="n">color</span><span class="p">,</span> <span class="n">edgecolor</span><span class="o">=</span><span class="n">color</span><span class="p">)</span>
|
||||
<span class="c1"># store the shape feature in the correct array</span>
|
||||
<span class="k">if</span> <span class="n">sig</span> <span class="ow">is</span> <span class="n">WARN_SIG</span><span class="p">:</span>
|
||||
<span class="n">warning_shapes</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">shape_feature</span><span class="p">)</span>
|
||||
<span class="k">elif</span> <span class="n">sig</span> <span class="ow">is</span> <span class="n">WATCH_SIG</span><span class="p">:</span>
|
||||
<span class="n">watch_shapes</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">shape_feature</span><span class="p">)</span>
|
||||
<span class="k">elif</span> <span class="n">sig</span> <span class="ow">is</span> <span class="n">ADVIS_SIG</span><span class="p">:</span>
|
||||
<span class="n">advisory_shapes</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">shape_feature</span><span class="p">)</span>
|
||||
<span class="k">elif</span> <span class="n">sig</span> <span class="ow">is</span> <span class="n">STATEM_SIG</span><span class="p">:</span>
|
||||
<span class="n">statement_shapes</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">shape_feature</span><span class="p">)</span>
|
||||
|
||||
<span class="c1"># Optional printout for the number of unique phensigs</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">phensigs</span><span class="p">),</span> <span class="s2">" Unique Phensigs"</span><span class="p">)</span>
|
||||
|
||||
<span class="c1"># Map phensigs to their titles (used for displaying alphabetically by</span>
|
||||
<span class="c1"># title in legend)</span>
|
||||
<span class="n">phensig_titles</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
<span class="k">for</span> <span class="n">phensig</span> <span class="ow">in</span> <span class="n">phensigs</span><span class="p">:</span>
|
||||
<span class="n">key</span> <span class="o">=</span> <span class="n">get_title</span><span class="p">(</span><span class="n">phensig</span><span class="p">)</span>
|
||||
<span class="n">phensig_titles</span><span class="p">[</span><span class="n">key</span><span class="p">]</span> <span class="o">=</span> <span class="n">phensig</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-default notranslate"><div class="highlight"><pre><span></span><span class="mi">14</span> <span class="n">Unique</span> <span class="n">Phensigs</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
<p><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html">Top</a></p>
|
||||
</section>
|
||||
</section>
|
||||
<hr class="docutils" />
|
||||
<section id="plot-the-data">
|
||||
<h3>8 Plot the Data!<a class="headerlink" href="#plot-the-data" title="Permalink to this headline"></a></h3>
|
||||
<section id="create-state-and-political-boundaries">
|
||||
<h4>8.1 Create State and Political Boundaries<a class="headerlink" href="#create-state-and-political-boundaries" title="Permalink to this headline"></a></h4>
|
||||
<p>Define the state and political boundaries that we’ll use in our plot to
|
||||
give more of a frame of reference. These objects are standard method
|
||||
calls in the <a class="reference external" href="https://scitools.org.uk/cartopy/docs/v0.14/matplotlib/feature_interface.html#cartopy.feature.NaturalEarthFeature">Cartopy Feature package, using the NaturalEarthFeature
|
||||
function</a>.</p>
|
||||
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="c1"># Define the state and political boundaries for the plot</span>
|
||||
<span class="n">states_provinces</span> <span class="o">=</span> <span class="n">cfeature</span><span class="o">.</span><span class="n">NaturalEarthFeature</span><span class="p">(</span>
|
||||
<span class="n">category</span><span class="o">=</span><span class="s1">'cultural'</span><span class="p">,</span>
|
||||
<span class="n">name</span><span class="o">=</span><span class="s1">'admin_1_states_provinces_lines'</span><span class="p">,</span>
|
||||
<span class="n">scale</span><span class="o">=</span><span class="s1">'50m'</span><span class="p">,</span>
|
||||
<span class="n">facecolor</span><span class="o">=</span><span class="s1">'none'</span><span class="p">)</span>
|
||||
<span class="n">political_boundaries</span> <span class="o">=</span> <span class="n">cfeature</span><span class="o">.</span><span class="n">NaturalEarthFeature</span><span class="p">(</span><span class="n">category</span><span class="o">=</span><span class="s1">'cultural'</span><span class="p">,</span>
|
||||
<span class="n">name</span><span class="o">=</span><span class="s1">'admin_0_boundary_lines_land'</span><span class="p">,</span>
|
||||
<span class="n">scale</span><span class="o">=</span><span class="s1">'50m'</span><span class="p">,</span> <span class="n">facecolor</span><span class="o">=</span><span class="s1">'none'</span><span class="p">)</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
</section>
|
||||
<section id="draw-the-plot-and-legend-for-wwas">
|
||||
<h4>8.2 Draw the Plot and Legend for WWAs<a class="headerlink" href="#draw-the-plot-and-legend-for-wwas" title="Permalink to this headline"></a></h4>
|
||||
<p>Here is where we finally get ot draw something! The very first few lines
|
||||
of this section are constants that we can manually “switch on and off”
|
||||
for what records we want displayed. By default we have all significance
|
||||
types drawn. If we want to “turn off” any of the significance records,
|
||||
simply set it’s corresponding constant to false, and re-run this cell to
|
||||
see how that plot compares.</p>
|
||||
<p>The next step involves creating the objects that are used to define the
|
||||
legend. We use the <code class="docutils literal notranslate"><span class="pre">phensig_titles</span></code> dictionary to loop through all the
|
||||
phensigs in alphabetical (by title) order. Then, we compare if the
|
||||
phensig will be displayed or not based on the display constants from the
|
||||
previous lines. If the significance will be drawn then we create a new
|
||||
<a class="reference external" href="https://matplotlib.org/stable/api/_as_gen/matplotlib.patches.Patch.html#matplotlib.patches.Patch">Patch
|
||||
object</a>
|
||||
of the corresponding color with the corresponding label and add it to
|
||||
our <code class="docutils literal notranslate"><span class="pre">handles</span></code> array.</p>
|
||||
<p>After that we define our bounding box and create our new plot with its
|
||||
figure and axes.</p>
|
||||
<p>Our next step is to create our Title for our plot. We create a title
|
||||
based on the draw variables to accurately describe what is being drawn
|
||||
in our plot. Here is where we use the first and last times defined in a
|
||||
previous cell.</p>
|
||||
<p>Finally, we create and show our plot. We add the title to the plot, add
|
||||
all the features to the axes, and add the legend as well.</p>
|
||||
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="c1"># Set these variables for which records to draw</span>
|
||||
<span class="n">DRAW_ADVISORY</span> <span class="o">=</span> <span class="kc">True</span>
|
||||
<span class="n">DRAW_WATCH</span> <span class="o">=</span> <span class="kc">True</span>
|
||||
<span class="n">DRAW_WARNING</span> <span class="o">=</span> <span class="kc">True</span>
|
||||
<span class="n">DRAW_STATEMENT</span> <span class="o">=</span> <span class="kc">True</span>
|
||||
|
||||
<span class="c1"># Create handles for legend and add items alphabetically by title and</span>
|
||||
<span class="c1"># only display based on the display values above</span>
|
||||
<span class="n">handles</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
<span class="k">for</span> <span class="n">title</span> <span class="ow">in</span> <span class="nb">sorted</span><span class="p">(</span><span class="n">phensig_titles</span><span class="p">):</span>
|
||||
<span class="n">phensig</span> <span class="o">=</span> <span class="n">phensig_titles</span><span class="p">[</span><span class="n">title</span><span class="p">]</span>
|
||||
<span class="c1"># check draw booleans</span>
|
||||
<span class="k">if</span> <span class="p">(</span> <span class="s2">"."</span><span class="o">+</span><span class="n">ADVIS_SIG</span> <span class="ow">in</span> <span class="n">phensig</span> <span class="ow">and</span> <span class="n">DRAW_ADVISORY</span> <span class="ow">or</span>
|
||||
<span class="s2">"."</span><span class="o">+</span><span class="n">WATCH_SIG</span> <span class="ow">in</span> <span class="n">phensig</span> <span class="ow">and</span> <span class="n">DRAW_WATCH</span> <span class="ow">or</span>
|
||||
<span class="s2">"."</span><span class="o">+</span><span class="n">WARN_SIG</span> <span class="ow">in</span> <span class="n">phensig</span> <span class="ow">and</span> <span class="n">DRAW_WARNING</span> <span class="ow">or</span>
|
||||
<span class="s2">"."</span><span class="o">+</span><span class="n">STATEM_SIG</span> <span class="ow">in</span> <span class="n">phensig</span> <span class="ow">and</span> <span class="n">DRAW_STATEMENT</span> <span class="p">):</span>
|
||||
<span class="n">entry</span> <span class="o">=</span> <span class="n">mpatches</span><span class="o">.</span><span class="n">Patch</span><span class="p">(</span><span class="n">color</span><span class="o">=</span><span class="n">get_color</span><span class="p">(</span><span class="n">phensig</span><span class="p">),</span> <span class="n">label</span><span class="o">=</span><span class="n">title</span><span class="p">)</span>
|
||||
<span class="n">handles</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">entry</span><span class="p">)</span>
|
||||
|
||||
<span class="c1"># Create the plot</span>
|
||||
<span class="n">bbox</span><span class="o">=</span><span class="p">[</span><span class="o">-</span><span class="mi">127</span><span class="p">,</span><span class="o">-</span><span class="mi">64</span><span class="p">,</span><span class="mi">24</span><span class="p">,</span><span class="mi">49</span><span class="p">]</span>
|
||||
<span class="n">fig</span><span class="p">,</span> <span class="n">ax</span> <span class="o">=</span> <span class="n">make_map</span><span class="p">(</span><span class="n">bbox</span><span class="o">=</span><span class="n">bbox</span><span class="p">)</span>
|
||||
|
||||
<span class="c1"># Add the title</span>
|
||||
<span class="c1"># Construct the title based on which record types are being displayed</span>
|
||||
<span class="n">title_string</span> <span class="o">=</span> <span class="s2">""</span>
|
||||
<span class="k">if</span> <span class="n">DRAW_WATCH</span><span class="p">:</span>
|
||||
<span class="n">title_string</span> <span class="o">+=</span> <span class="s2">"Watches, "</span>
|
||||
<span class="k">if</span> <span class="n">DRAW_WARNING</span><span class="p">:</span>
|
||||
<span class="n">title_string</span> <span class="o">+=</span> <span class="s2">"Warnings, "</span>
|
||||
<span class="k">if</span> <span class="n">DRAW_ADVISORY</span><span class="p">:</span>
|
||||
<span class="n">title_string</span> <span class="o">+=</span> <span class="s2">"Advisories, "</span>
|
||||
<span class="k">if</span> <span class="n">DRAW_STATEMENT</span><span class="p">:</span>
|
||||
<span class="n">title_string</span> <span class="o">+=</span> <span class="s2">"Statements, "</span>
|
||||
<span class="c1"># remove the last comma and space</span>
|
||||
<span class="n">title_string</span> <span class="o">=</span> <span class="n">title_string</span><span class="p">[:</span><span class="o">-</span><span class="mi">2</span><span class="p">]</span>
|
||||
<span class="c1"># add the time range</span>
|
||||
<span class="n">title_string</span> <span class="o">+=</span> <span class="s2">" from "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">first_time</span><span class="p">)[:</span><span class="o">-</span><span class="mi">3</span><span class="p">]</span> <span class="o">+</span> <span class="s2">" to "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">last_time</span><span class="p">)[:</span><span class="o">-</span><span class="mi">3</span><span class="p">]</span> <span class="o">+</span> <span class="s2">" UTC"</span>
|
||||
<span class="c1"># set the title on the plot, give it a bigger font size, and increase</span>
|
||||
<span class="c1"># the vertical padding between the title and the figure</span>
|
||||
<span class="n">plt</span><span class="o">.</span><span class="n">title</span><span class="p">(</span><span class="n">title_string</span><span class="p">,</span> <span class="n">fontsize</span><span class="o">=</span><span class="mi">20</span><span class="p">,</span> <span class="n">pad</span><span class="o">=</span><span class="mi">10</span><span class="p">)</span>
|
||||
|
||||
<span class="c1"># Draw all features on the plot</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">add_feature</span><span class="p">(</span><span class="n">cfeature</span><span class="o">.</span><span class="n">LAND</span><span class="p">)</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">add_feature</span><span class="p">(</span><span class="n">cfeature</span><span class="o">.</span><span class="n">COASTLINE</span><span class="p">)</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">add_feature</span><span class="p">(</span><span class="n">states_provinces</span><span class="p">,</span> <span class="n">edgecolor</span><span class="o">=</span><span class="s1">'black'</span><span class="p">)</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">add_feature</span><span class="p">(</span><span class="n">political_boundaries</span><span class="p">,</span> <span class="n">edgecolor</span><span class="o">=</span><span class="s1">'black'</span><span class="p">)</span>
|
||||
|
||||
<span class="c1"># Draw WWAs in order: Advisory -> Watch > Warning > Statement</span>
|
||||
<span class="k">if</span> <span class="n">DRAW_ADVISORY</span><span class="p">:</span>
|
||||
<span class="k">for</span> <span class="n">shape</span> <span class="ow">in</span> <span class="n">advisory_shapes</span><span class="p">:</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">add_feature</span><span class="p">(</span><span class="n">shape</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="n">DRAW_WATCH</span><span class="p">:</span>
|
||||
<span class="k">for</span> <span class="n">shape</span> <span class="ow">in</span> <span class="n">watch_shapes</span><span class="p">:</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">add_feature</span><span class="p">(</span><span class="n">shape</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="n">DRAW_WARNING</span><span class="p">:</span>
|
||||
<span class="k">for</span> <span class="n">shape</span> <span class="ow">in</span> <span class="n">warning_shapes</span><span class="p">:</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">add_feature</span><span class="p">(</span><span class="n">shape</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="n">DRAW_STATEMENT</span><span class="p">:</span>
|
||||
<span class="k">for</span> <span class="n">shape</span> <span class="ow">in</span> <span class="n">statement_shapes</span><span class="p">:</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">add_feature</span><span class="p">(</span><span class="n">shape</span><span class="p">)</span>
|
||||
|
||||
<span class="c1"># Draw the legend</span>
|
||||
<span class="c1"># use the handles defined earlier for the color associations to</span>
|
||||
<span class="c1"># phensig titles, set the location to the lower center, give it</span>
|
||||
<span class="c1"># 5 columns so it uses all the horizonatal space, place it under</span>
|
||||
<span class="c1"># the actual figure, and give it a larger fontsize</span>
|
||||
<span class="n">bottom</span> <span class="o">=</span> <span class="mf">0.12</span> <span class="o">+</span> <span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">handles</span><span class="p">)</span><span class="o">//</span><span class="mi">5</span> <span class="o">*</span><span class="mf">.04</span><span class="p">)</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">legend</span><span class="p">(</span><span class="n">handles</span><span class="o">=</span><span class="n">handles</span><span class="p">,</span> <span class="n">loc</span><span class="o">=</span><span class="s1">'lower center'</span><span class="p">,</span> <span class="n">ncol</span><span class="o">=</span><span class="mi">5</span><span class="p">,</span> <span class="n">bbox_to_anchor</span><span class="o">=</span><span class="p">(</span><span class="mf">0.5</span><span class="p">,</span> <span class="o">-</span><span class="n">bottom</span><span class="p">),</span> <span class="n">fontsize</span><span class="o">=</span><span class="mi">16</span><span class="p">)</span>
|
||||
|
||||
<span class="c1"># Show the plot</span>
|
||||
<span class="n">plt</span><span class="o">.</span><span class="n">show</span><span class="p">()</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
<img alt="../../_images/Watch_Warning_and_Advisory_Plotting_34_0.png" src="../../_images/Watch_Warning_and_Advisory_Plotting_34_0.png" />
|
||||
<p><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html">Top</a></p>
|
||||
</section>
|
||||
</section>
|
||||
<hr class="docutils" />
|
||||
<section id="see-also">
|
||||
<h3>9 See Also<a class="headerlink" href="#see-also" title="Permalink to this headline"></a></h3>
|
||||
<ul class="simple">
|
||||
<li><p><a class="reference external" href="https://www.weather.gov/lwx/warningsdefined#:~:text=A%20Winter%20Storm%20Warning%20is%20issued%20when%20a,combination%20of%20snow%20and%2For%20ice%20accumulation%20with%20wind.">National Weather Service WWA Definitions (Baltimore
|
||||
Office)</a></p></li>
|
||||
<li><p><a class="reference external" href="https://weather.cod.edu/notes/criteria/">College of Dupage WWA
|
||||
Definitions</a></p></li>
|
||||
<li><p><a class="reference external" href="https://www.weather.gov/media/vtec/VTEC_explanation4-20.pdf">Phensig
|
||||
Explanation</a></p></li>
|
||||
</ul>
|
||||
<section id="related-notebooks">
|
||||
<h4>9.1 Related Notebooks<a class="headerlink" href="#related-notebooks" title="Permalink to this headline"></a></h4>
|
||||
<ul class="simple">
|
||||
<li><p><a class="reference external" href="http://unidata.github.io/python-awips/examples/generated/Grid_Levels_and_Parameters.html">Grid Levels and
|
||||
Parameters</a></p></li>
|
||||
</ul>
|
||||
</section>
|
||||
<section id="additional-documentation">
|
||||
<h4>9.2 Additional Documentation<a class="headerlink" href="#additional-documentation" title="Permalink to this headline"></a></h4>
|
||||
<p><strong>python-awips</strong></p>
|
||||
<ul class="simple">
|
||||
<li><p><a class="reference external" href="http://unidata.github.io/python-awips/api/DataAccessLayer.html#awips.dataaccess.DataAccessLayer.changeEDEXHost">DataAccessLayer.changeEDEXHost()</a></p></li>
|
||||
<li><p><a class="reference external" href="http://unidata.github.io/python-awips/api/DataAccessLayer.html#awips.dataaccess.DataAccessLayer.newDataRequest">DataAccessLayer.newDataRequest()</a></p></li>
|
||||
<li><p><a class="reference external" href="http://unidata.github.io/python-awips/api/DataAccessLayer.html#awips.dataaccess.DataAccessLayer.getAvailableParameters">DataAccessLayer.getAvailableParameters()</a></p></li>
|
||||
<li><p><a class="reference external" href="http://unidata.github.io/python-awips/api/IDataRequest.html">IDataRequest</a></p></li>
|
||||
<li><p><a class="reference external" href="http://unidata.github.io/python-awips/api/PyGeometryData.html">GeometryData</a></p></li>
|
||||
</ul>
|
||||
<p><strong>datetime</strong></p>
|
||||
<ul class="simple">
|
||||
<li><p><a class="reference external" href="https://docs.python.org/3/library/datetime.html#datetime-objects">datetime.datetime</a></p></li>
|
||||
<li><p><a class="reference external" href="https://docs.python.org/3/library/datetime.html#timedelta-objects">datetime.timedelta</a></p></li>
|
||||
</ul>
|
||||
<p><strong>cartopy</strong></p>
|
||||
<ul class="simple">
|
||||
<li><p><a class="reference external" href="https://scitools.org.uk/cartopy/docs/v0.14/matplotlib/feature_interface.html">cartopy feature
|
||||
interface</a></p></li>
|
||||
<li><p><a class="reference external" href="https://scitools.org.uk/cartopy/docs/latest/reference/generated/cartopy.feature.ShapelyFeature.html">cartopy.feature.ShaeplyFeature</a></p></li>
|
||||
</ul>
|
||||
<p><strong>matplotlib</strong></p>
|
||||
<ul class="simple">
|
||||
<li><p><a class="reference external" href="https://matplotlib.org/stable/api/_as_gen/matplotlib.pyplot.html">matplotlib.pyplot()</a></p></li>
|
||||
<li><p><a class="reference external" href="https://matplotlib.org/3.1.1/api/_as_gen/matplotlib.pyplot.legend.html">matplotlib.pyplot.legend()</a></p></li>
|
||||
<li><p><a class="reference external" href="https://matplotlib.org/stable/api/_as_gen/matplotlib.pyplot.axes.html">matplotlib.pyplot.axes()</a></p></li>
|
||||
<li><p><a class="reference external" href="https://matplotlib.org/stable/api/_as_gen/matplotlib.pyplot.figure.html">matplotlib.pyplot.figure()</a></p></li>
|
||||
<li><p><a class="reference external" href="https://matplotlib.org/3.1.1/api/_as_gen/matplotlib.pyplot.title.html">matplotlib.pyplot.title()</a></p></li>
|
||||
<li><p><a class="reference external" href="https://matplotlib.org/stable/api/_as_gen/matplotlib.patches.Patch.html#matplotlib.patches.Patch">matplotlib.pathes.Patch</a></p></li>
|
||||
</ul>
|
||||
<p><a class="reference external" href="https://unidata.github.io/python-awips/examples/generated/Watch_Warning_and_Advisory_Plotting.html">Top</a></p>
|
||||
<hr class="docutils" />
|
||||
</section>
|
||||
</section>
|
||||
</section>
|
||||
</section>
|
||||
|
||||
|
||||
</div>
|
||||
</div>
|
||||
<footer><div class="rst-footer-buttons" role="navigation" aria-label="Footer">
|
||||
<a href="Upper_Air_BUFR_Soundings.html" class="btn btn-neutral float-left" title="Upper Air BUFR Soundings" accesskey="p" rel="prev"><span class="fa fa-arrow-circle-left" aria-hidden="true"></span> Previous</a>
|
||||
<a href="../../dev.html" class="btn btn-neutral float-right" title="Development Guide" accesskey="n" rel="next">Next <span class="fa fa-arrow-circle-right" aria-hidden="true"></span></a>
|
||||
</div>
|
||||
|
||||
<hr/>
|
||||
|
||||
<div role="contentinfo">
|
||||
<p>© Copyright 2018, Unidata.</p>
|
||||
</div>
|
||||
|
||||
Built with <a href="https://www.sphinx-doc.org/">Sphinx</a> using a
|
||||
<a href="https://github.com/readthedocs/sphinx_rtd_theme">theme</a>
|
||||
provided by <a href="https://readthedocs.org">Read the Docs</a>.
|
||||
|
||||
|
||||
</footer>
|
||||
</div>
|
||||
</div>
|
||||
</section>
|
||||
</div>
|
||||
<script>
|
||||
jQuery(function () {
|
||||
SphinxRtdTheme.Navigation.enable(true);
|
||||
});
|
||||
</script>
|
||||
|
||||
</body>
|
||||
</html>
|
||||
|
|
@ -1,253 +0,0 @@
|
|||
<!DOCTYPE html>
|
||||
<html class="writer-html5" lang="en" >
|
||||
<head>
|
||||
<meta charset="utf-8" /><meta name="generator" content="Docutils 0.17.1: http://docutils.sourceforge.net/" />
|
||||
|
||||
<meta name="viewport" content="width=device-width, initial-scale=1.0" />
|
||||
<title>Watch and Warning Polygons — python-awips documentation</title>
|
||||
<link rel="stylesheet" href="../../_static/pygments.css" type="text/css" />
|
||||
<link rel="stylesheet" href="../../_static/css/theme.css" type="text/css" />
|
||||
<!--[if lt IE 9]>
|
||||
<script src="../../_static/js/html5shiv.min.js"></script>
|
||||
<![endif]-->
|
||||
|
||||
<script data-url_root="../../" id="documentation_options" src="../../_static/documentation_options.js"></script>
|
||||
<script src="../../_static/jquery.js"></script>
|
||||
<script src="../../_static/underscore.js"></script>
|
||||
<script src="../../_static/doctools.js"></script>
|
||||
<script src="../../_static/js/theme.js"></script>
|
||||
<link rel="author" title="About these documents" href="../../about.html" />
|
||||
<link rel="index" title="Index" href="../../genindex.html" />
|
||||
<link rel="search" title="Search" href="../../search.html" />
|
||||
<link rel="next" title="Development Guide" href="../../dev.html" />
|
||||
<link rel="prev" title="Upper Air BUFR Soundings" href="Upper_Air_BUFR_Soundings.html" />
|
||||
</head>
|
||||
|
||||
<body class="wy-body-for-nav">
|
||||
<div class="wy-grid-for-nav">
|
||||
<nav data-toggle="wy-nav-shift" class="wy-nav-side">
|
||||
<div class="wy-side-scroll">
|
||||
<div class="wy-side-nav-search" >
|
||||
<a href="../../index.html" class="icon icon-home"> python-awips
|
||||
</a>
|
||||
<div class="version">
|
||||
18.1.8
|
||||
</div>
|
||||
<div role="search">
|
||||
<form id="rtd-search-form" class="wy-form" action="../../search.html" method="get">
|
||||
<input type="text" name="q" placeholder="Search docs" />
|
||||
<input type="hidden" name="check_keywords" value="yes" />
|
||||
<input type="hidden" name="area" value="default" />
|
||||
</form>
|
||||
</div>
|
||||
</div><div class="wy-menu wy-menu-vertical" data-spy="affix" role="navigation" aria-label="Navigation menu">
|
||||
<ul class="current">
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../api/index.html">API Documentation</a></li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../datatypes.html">Available Data Types</a></li>
|
||||
<li class="toctree-l1 current"><a class="reference internal" href="../index.html">Data Plotting Examples</a><ul class="current">
|
||||
<li class="toctree-l2"><a class="reference internal" href="Colored_Surface_Temperature_Plot.html">Colored Surface Temperature Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Colorized_Grid_Data.html">Colorized Grid Data</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Forecast_Model_Vertical_Sounding.html">Forecast Model Vertical Sounding</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="GOES_CIRA_Product_Writer.html">GOES CIRA Product Writer</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="GOES_Geostationary_Lightning_Mapper.html">GOES Geostationary Lightning Mapper</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Grid_Levels_and_Parameters.html">Grid Levels and Parameters</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="METAR_Station_Plot_with_MetPy.html">METAR Station Plot with MetPy</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Map_Resources_and_Topography.html">Map Resources and Topography</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Model_Sounding_Data.html">Model Sounding Data</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="NEXRAD_Level3_Radar.html">NEXRAD Level3 Radar</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Precip_Accumulation-Region_Of_Interest.html">Precip Accumulation-Region Of Interest</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2 current"><a class="current reference internal" href="#">Watch and Warning Polygons</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../dev.html">Development Guide</a></li>
|
||||
<li class="toctree-l1"><a class="reference external" href="http://unidata.github.io/awips2/appendix/appendix-grid-parameters/">AWIPS Grid Parameters</a></li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../../about.html">About Unidata AWIPS</a></li>
|
||||
</ul>
|
||||
|
||||
</div>
|
||||
</div>
|
||||
</nav>
|
||||
|
||||
<section data-toggle="wy-nav-shift" class="wy-nav-content-wrap"><nav class="wy-nav-top" aria-label="Mobile navigation menu" >
|
||||
<i data-toggle="wy-nav-top" class="fa fa-bars"></i>
|
||||
<a href="../../index.html">python-awips</a>
|
||||
</nav>
|
||||
|
||||
<div class="wy-nav-content">
|
||||
<div class="rst-content">
|
||||
<div role="navigation" aria-label="Page navigation">
|
||||
<ul class="wy-breadcrumbs">
|
||||
<li><a href="../../index.html" class="icon icon-home"></a> »</li>
|
||||
<li><a href="../index.html">Data Plotting Examples</a> »</li>
|
||||
<li>Watch and Warning Polygons</li>
|
||||
<li class="wy-breadcrumbs-aside">
|
||||
<a href="../../_sources/examples/generated/Watch_and_Warning_Polygons.rst.txt" rel="nofollow"> View page source</a>
|
||||
</li>
|
||||
</ul>
|
||||
<hr/>
|
||||
</div>
|
||||
<div role="main" class="document" itemscope="itemscope" itemtype="http://schema.org/Article">
|
||||
<div itemprop="articleBody">
|
||||
|
||||
<section id="watch-and-warning-polygons">
|
||||
<h1>Watch and Warning Polygons<a class="headerlink" href="#watch-and-warning-polygons" title="Permalink to this headline"></a></h1>
|
||||
<p><a class="reference external" href="http://nbviewer.ipython.org/github/Unidata/python-awips/blob/master/examples/notebooks/Watch_and_Warning_Polygons.ipynb">Notebook</a>
|
||||
This example uses matplotlib, cartopy, shapely, and python-awips to plot
|
||||
watch and warning polygons requested from a real-time AWIPS EDEX server.</p>
|
||||
<p>First, set up our imports and define functions to be used later:</p>
|
||||
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="kn">from</span> <span class="nn">awips.dataaccess</span> <span class="kn">import</span> <span class="n">DataAccessLayer</span>
|
||||
<span class="kn">from</span> <span class="nn">awips.tables</span> <span class="kn">import</span> <span class="n">vtec</span>
|
||||
<span class="kn">from</span> <span class="nn">datetime</span> <span class="kn">import</span> <span class="n">datetime</span>
|
||||
<span class="kn">import</span> <span class="nn">numpy</span> <span class="k">as</span> <span class="nn">np</span>
|
||||
<span class="kn">import</span> <span class="nn">matplotlib.pyplot</span> <span class="k">as</span> <span class="nn">plt</span>
|
||||
<span class="kn">import</span> <span class="nn">cartopy.crs</span> <span class="k">as</span> <span class="nn">ccrs</span>
|
||||
<span class="kn">import</span> <span class="nn">cartopy.feature</span> <span class="k">as</span> <span class="nn">cfeature</span>
|
||||
<span class="kn">from</span> <span class="nn">cartopy.mpl.gridliner</span> <span class="kn">import</span> <span class="n">LONGITUDE_FORMATTER</span><span class="p">,</span> <span class="n">LATITUDE_FORMATTER</span>
|
||||
<span class="kn">from</span> <span class="nn">cartopy.feature</span> <span class="kn">import</span> <span class="n">ShapelyFeature</span><span class="p">,</span><span class="n">NaturalEarthFeature</span>
|
||||
<span class="kn">from</span> <span class="nn">shapely.geometry</span> <span class="kn">import</span> <span class="n">MultiPolygon</span><span class="p">,</span><span class="n">Polygon</span>
|
||||
|
||||
<span class="k">def</span> <span class="nf">warning_color</span><span class="p">(</span><span class="n">phensig</span><span class="p">):</span>
|
||||
<span class="k">return</span> <span class="n">vtec</span><span class="p">[</span><span class="n">phensig</span><span class="p">][</span><span class="s1">'color'</span><span class="p">]</span>
|
||||
|
||||
<span class="k">def</span> <span class="nf">make_map</span><span class="p">(</span><span class="n">bbox</span><span class="p">,</span> <span class="n">projection</span><span class="o">=</span><span class="n">ccrs</span><span class="o">.</span><span class="n">PlateCarree</span><span class="p">()):</span>
|
||||
<span class="n">fig</span><span class="p">,</span> <span class="n">ax</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">subplots</span><span class="p">(</span><span class="n">figsize</span><span class="o">=</span><span class="p">(</span><span class="mi">20</span><span class="p">,</span><span class="mi">12</span><span class="p">),</span>
|
||||
<span class="n">subplot_kw</span><span class="o">=</span><span class="nb">dict</span><span class="p">(</span><span class="n">projection</span><span class="o">=</span><span class="n">projection</span><span class="p">))</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">set_extent</span><span class="p">(</span><span class="n">bbox</span><span class="p">)</span>
|
||||
<span class="n">gl</span> <span class="o">=</span> <span class="n">ax</span><span class="o">.</span><span class="n">gridlines</span><span class="p">(</span><span class="n">draw_labels</span><span class="o">=</span><span class="kc">True</span><span class="p">)</span>
|
||||
<span class="n">gl</span><span class="o">.</span><span class="n">top_labels</span> <span class="o">=</span> <span class="n">gl</span><span class="o">.</span><span class="n">right_labels</span> <span class="o">=</span> <span class="kc">False</span>
|
||||
<span class="n">gl</span><span class="o">.</span><span class="n">xformatter</span> <span class="o">=</span> <span class="n">LONGITUDE_FORMATTER</span>
|
||||
<span class="n">gl</span><span class="o">.</span><span class="n">yformatter</span> <span class="o">=</span> <span class="n">LATITUDE_FORMATTER</span>
|
||||
<span class="k">return</span> <span class="n">fig</span><span class="p">,</span> <span class="n">ax</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
<p>Next, we create a request for the “warning” data type:</p>
|
||||
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="n">DataAccessLayer</span><span class="o">.</span><span class="n">changeEDEXHost</span><span class="p">(</span><span class="s2">"edex-cloud.unidata.ucar.edu"</span><span class="p">)</span>
|
||||
<span class="n">request</span> <span class="o">=</span> <span class="n">DataAccessLayer</span><span class="o">.</span><span class="n">newDataRequest</span><span class="p">()</span>
|
||||
<span class="n">request</span><span class="o">.</span><span class="n">setDatatype</span><span class="p">(</span><span class="s2">"warning"</span><span class="p">)</span>
|
||||
<span class="n">request</span><span class="o">.</span><span class="n">setParameters</span><span class="p">(</span><span class="s1">'phensig'</span><span class="p">)</span>
|
||||
<span class="n">times</span> <span class="o">=</span> <span class="n">DataAccessLayer</span><span class="o">.</span><span class="n">getAvailableTimes</span><span class="p">(</span><span class="n">request</span><span class="p">)</span>
|
||||
|
||||
<span class="c1"># Get records for last 50 available times</span>
|
||||
<span class="n">response</span> <span class="o">=</span> <span class="n">DataAccessLayer</span><span class="o">.</span><span class="n">getGeometryData</span><span class="p">(</span><span class="n">request</span><span class="p">,</span> <span class="n">times</span><span class="p">[</span><span class="o">-</span><span class="mi">50</span><span class="p">:</span><span class="o">-</span><span class="mi">1</span><span class="p">])</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="s2">"Using "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">response</span><span class="p">))</span> <span class="o">+</span> <span class="s2">" records"</span><span class="p">)</span>
|
||||
|
||||
<span class="c1"># Each record will have a numpy array the length of the number of "parameters"</span>
|
||||
<span class="c1"># Default is 1 (request.setParameters('phensig'))</span>
|
||||
<span class="n">parameters</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
<span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">request</span><span class="o">.</span><span class="n">getParameters</span><span class="p">():</span>
|
||||
<span class="n">parameters</span><span class="p">[</span><span class="n">x</span><span class="p">]</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([])</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="n">parameters</span><span class="p">)</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-default notranslate"><div class="highlight"><pre><span></span><span class="n">Using</span> <span class="mi">109</span> <span class="n">records</span>
|
||||
<span class="p">{</span><span class="s1">'phensig'</span><span class="p">:</span> <span class="n">array</span><span class="p">([],</span> <span class="n">dtype</span><span class="o">=</span><span class="n">float64</span><span class="p">)}</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
<p>Now loop through each record and plot it as either Polygon or
|
||||
MultiPolygon, with appropriate colors</p>
|
||||
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="o">%</span><span class="k">matplotlib</span> inline
|
||||
<span class="n">bbox</span><span class="o">=</span><span class="p">[</span><span class="o">-</span><span class="mi">127</span><span class="p">,</span><span class="o">-</span><span class="mi">64</span><span class="p">,</span><span class="mi">24</span><span class="p">,</span><span class="mi">49</span><span class="p">]</span>
|
||||
<span class="n">fig</span><span class="p">,</span> <span class="n">ax</span> <span class="o">=</span> <span class="n">make_map</span><span class="p">(</span><span class="n">bbox</span><span class="o">=</span><span class="n">bbox</span><span class="p">)</span>
|
||||
|
||||
<span class="n">siteids</span><span class="o">=</span><span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([])</span>
|
||||
<span class="n">periods</span><span class="o">=</span><span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([])</span>
|
||||
<span class="n">reftimes</span><span class="o">=</span><span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([])</span>
|
||||
|
||||
<span class="k">for</span> <span class="n">ob</span> <span class="ow">in</span> <span class="n">response</span><span class="p">:</span>
|
||||
|
||||
<span class="n">poly</span> <span class="o">=</span> <span class="n">ob</span><span class="o">.</span><span class="n">getGeometry</span><span class="p">()</span>
|
||||
<span class="n">site</span> <span class="o">=</span> <span class="n">ob</span><span class="o">.</span><span class="n">getLocationName</span><span class="p">()</span>
|
||||
<span class="n">pd</span> <span class="o">=</span> <span class="n">ob</span><span class="o">.</span><span class="n">getDataTime</span><span class="p">()</span><span class="o">.</span><span class="n">getValidPeriod</span><span class="p">()</span>
|
||||
<span class="n">ref</span> <span class="o">=</span> <span class="n">ob</span><span class="o">.</span><span class="n">getDataTime</span><span class="p">()</span><span class="o">.</span><span class="n">getRefTime</span><span class="p">()</span>
|
||||
|
||||
<span class="c1"># do not plot if phensig is blank (SPS)</span>
|
||||
<span class="k">if</span> <span class="n">ob</span><span class="o">.</span><span class="n">getString</span><span class="p">(</span><span class="s1">'phensig'</span><span class="p">):</span>
|
||||
|
||||
<span class="n">phensigString</span> <span class="o">=</span> <span class="n">ob</span><span class="o">.</span><span class="n">getString</span><span class="p">(</span><span class="s1">'phensig'</span><span class="p">)</span>
|
||||
|
||||
<span class="n">siteids</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">siteids</span><span class="p">,</span><span class="n">site</span><span class="p">)</span>
|
||||
<span class="n">periods</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">periods</span><span class="p">,</span><span class="n">pd</span><span class="p">)</span>
|
||||
<span class="n">reftimes</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">reftimes</span><span class="p">,</span><span class="n">ref</span><span class="p">)</span>
|
||||
|
||||
<span class="k">for</span> <span class="n">parm</span> <span class="ow">in</span> <span class="n">parameters</span><span class="p">:</span>
|
||||
<span class="n">parameters</span><span class="p">[</span><span class="n">parm</span><span class="p">]</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">parameters</span><span class="p">[</span><span class="n">parm</span><span class="p">],</span><span class="n">ob</span><span class="o">.</span><span class="n">getString</span><span class="p">(</span><span class="n">parm</span><span class="p">))</span>
|
||||
|
||||
<span class="k">if</span> <span class="n">poly</span><span class="o">.</span><span class="n">geom_type</span> <span class="o">==</span> <span class="s1">'MultiPolygon'</span><span class="p">:</span>
|
||||
<span class="n">geometries</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([])</span>
|
||||
<span class="n">geometries</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">geometries</span><span class="p">,</span><span class="n">MultiPolygon</span><span class="p">(</span><span class="n">poly</span><span class="p">))</span>
|
||||
<span class="n">geom_count</span> <span class="o">=</span> <span class="s2">", "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">geometries</span><span class="p">))</span> <span class="o">+</span><span class="s2">" geometries"</span>
|
||||
<span class="k">else</span><span class="p">:</span>
|
||||
<span class="n">geometries</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([])</span>
|
||||
<span class="n">geometries</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">geometries</span><span class="p">,</span><span class="n">Polygon</span><span class="p">(</span><span class="n">poly</span><span class="p">))</span>
|
||||
<span class="n">geom_count</span><span class="o">=</span><span class="s2">""</span>
|
||||
|
||||
<span class="k">for</span> <span class="n">geom</span> <span class="ow">in</span> <span class="n">geometries</span><span class="p">:</span>
|
||||
<span class="n">bounds</span> <span class="o">=</span> <span class="n">Polygon</span><span class="p">(</span><span class="n">geom</span><span class="p">)</span>
|
||||
<span class="n">intersection</span> <span class="o">=</span> <span class="n">bounds</span><span class="o">.</span><span class="n">intersection</span>
|
||||
<span class="n">geoms</span> <span class="o">=</span> <span class="p">(</span><span class="n">intersection</span><span class="p">(</span><span class="n">geom</span><span class="p">)</span>
|
||||
<span class="k">for</span> <span class="n">geom</span> <span class="ow">in</span> <span class="n">geometries</span>
|
||||
<span class="k">if</span> <span class="n">bounds</span><span class="o">.</span><span class="n">intersects</span><span class="p">(</span><span class="n">geom</span><span class="p">))</span>
|
||||
|
||||
<span class="c1">#print(vtec[phensigString]['hdln']</span>
|
||||
<span class="c1"># + " (" + phensigString + ") issued at " + str(ref)</span>
|
||||
<span class="c1"># + " ("+str(poly.geom_type) + geom_count + ")")</span>
|
||||
|
||||
<span class="n">color</span> <span class="o">=</span> <span class="n">warning_color</span><span class="p">(</span><span class="n">phensigString</span><span class="p">)</span>
|
||||
<span class="n">shape_feature</span> <span class="o">=</span> <span class="n">ShapelyFeature</span><span class="p">(</span><span class="n">geoms</span><span class="p">,</span><span class="n">ccrs</span><span class="o">.</span><span class="n">PlateCarree</span><span class="p">(),</span>
|
||||
<span class="n">facecolor</span><span class="o">=</span><span class="n">color</span><span class="p">,</span> <span class="n">edgecolor</span><span class="o">=</span><span class="n">color</span><span class="p">)</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">add_feature</span><span class="p">(</span><span class="n">shape_feature</span><span class="p">)</span>
|
||||
|
||||
<span class="n">states_provinces</span> <span class="o">=</span> <span class="n">cfeature</span><span class="o">.</span><span class="n">NaturalEarthFeature</span><span class="p">(</span>
|
||||
<span class="n">category</span><span class="o">=</span><span class="s1">'cultural'</span><span class="p">,</span>
|
||||
<span class="n">name</span><span class="o">=</span><span class="s1">'admin_1_states_provinces_lines'</span><span class="p">,</span>
|
||||
<span class="n">scale</span><span class="o">=</span><span class="s1">'50m'</span><span class="p">,</span>
|
||||
<span class="n">facecolor</span><span class="o">=</span><span class="s1">'none'</span><span class="p">)</span>
|
||||
<span class="n">political_boundaries</span> <span class="o">=</span> <span class="n">cfeature</span><span class="o">.</span><span class="n">NaturalEarthFeature</span><span class="p">(</span><span class="n">category</span><span class="o">=</span><span class="s1">'cultural'</span><span class="p">,</span>
|
||||
<span class="n">name</span><span class="o">=</span><span class="s1">'admin_0_boundary_lines_land'</span><span class="p">,</span>
|
||||
<span class="n">scale</span><span class="o">=</span><span class="s1">'50m'</span><span class="p">,</span> <span class="n">facecolor</span><span class="o">=</span><span class="s1">'none'</span><span class="p">)</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">add_feature</span><span class="p">(</span><span class="n">cfeature</span><span class="o">.</span><span class="n">LAND</span><span class="p">)</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">add_feature</span><span class="p">(</span><span class="n">cfeature</span><span class="o">.</span><span class="n">COASTLINE</span><span class="p">)</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">add_feature</span><span class="p">(</span><span class="n">states_provinces</span><span class="p">,</span> <span class="n">edgecolor</span><span class="o">=</span><span class="s1">'black'</span><span class="p">)</span>
|
||||
<span class="n">ax</span><span class="o">.</span><span class="n">add_feature</span><span class="p">(</span><span class="n">political_boundaries</span><span class="p">,</span> <span class="n">edgecolor</span><span class="o">=</span><span class="s1">'black'</span><span class="p">)</span>
|
||||
|
||||
<span class="n">plt</span><span class="o">.</span><span class="n">show</span><span class="p">()</span>
|
||||
</pre></div>
|
||||
</div>
|
||||
<img alt="../../_images/Watch_and_Warning_Polygons_5_0.png" src="../../_images/Watch_and_Warning_Polygons_5_0.png" />
|
||||
</section>
|
||||
|
||||
|
||||
</div>
|
||||
</div>
|
||||
<footer><div class="rst-footer-buttons" role="navigation" aria-label="Footer">
|
||||
<a href="Upper_Air_BUFR_Soundings.html" class="btn btn-neutral float-left" title="Upper Air BUFR Soundings" accesskey="p" rel="prev"><span class="fa fa-arrow-circle-left" aria-hidden="true"></span> Previous</a>
|
||||
<a href="../../dev.html" class="btn btn-neutral float-right" title="Development Guide" accesskey="n" rel="next">Next <span class="fa fa-arrow-circle-right" aria-hidden="true"></span></a>
|
||||
</div>
|
||||
|
||||
<hr/>
|
||||
|
||||
<div role="contentinfo">
|
||||
<p>© Copyright 2018, Unidata.</p>
|
||||
</div>
|
||||
|
||||
Built with <a href="https://www.sphinx-doc.org/">Sphinx</a> using a
|
||||
<a href="https://github.com/readthedocs/sphinx_rtd_theme">theme</a>
|
||||
provided by <a href="https://readthedocs.org">Read the Docs</a>.
|
||||
|
||||
|
||||
</footer>
|
||||
</div>
|
||||
</div>
|
||||
</section>
|
||||
</div>
|
||||
<script>
|
||||
jQuery(function () {
|
||||
SphinxRtdTheme.Navigation.enable(true);
|
||||
});
|
||||
</script>
|
||||
|
||||
</body>
|
||||
</html>
|
||||
|
|
@ -59,7 +59,7 @@
|
|||
<li class="toctree-l2"><a class="reference internal" href="generated/Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="generated/Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="generated/Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="generated/Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l2"><a class="reference internal" href="generated/Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="../dev.html">Development Guide</a></li>
|
||||
|
|
@ -109,7 +109,7 @@
|
|||
<li class="toctree-l1"><a class="reference internal" href="generated/Regional_Surface_Obs_Plot.html">Regional Surface Obs Plot</a></li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="generated/Satellite_Imagery.html">Satellite Imagery</a></li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="generated/Upper_Air_BUFR_Soundings.html">Upper Air BUFR Soundings</a></li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="generated/Watch_and_Warning_Polygons.html">Watch and Warning Polygons</a></li>
|
||||
<li class="toctree-l1"><a class="reference internal" href="generated/Watch_Warning_and_Advisory_Plotting.html">Watch Warning and Advisory Plotting</a></li>
|
||||
</ul>
|
||||
</div>
|
||||
</section>
|
||||
|
|
|
|||
BIN
objects.inv
BIN
objects.inv
Binary file not shown.
File diff suppressed because one or more lines are too long
Loading…
Add table
Reference in a new issue