.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "examples/recon.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_examples_recon.py>`
        to download the full example code

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_examples_recon.py:


============================
Aircraft Recon Data Analysis
============================
This sample script shows how to use Tropycal to retrieve, analyze and plot aircraft reconnaissance missions.

.. GENERATED FROM PYTHON SOURCE LINES 7-12

.. code-block:: default


    from tropycal import tracks, recon
    import numpy as np
    import datetime as dt








.. GENERATED FROM PYTHON SOURCE LINES 13-16

Reading In Recon Dataset
------------------------
There are several ways of retrieving a ReconDataset object. The first step is to retrieve a TrackDataset object, which for this example we'll use the HURDATv2 database for the North Atlantic basin.

.. GENERATED FROM PYTHON SOURCE LINES 16-19

.. code-block:: default


    basin = tracks.TrackDataset('north_atlantic')





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    --> Starting to read in HURDAT2 data
    --> Completed reading in HURDAT2 data (1.99 seconds)




.. GENERATED FROM PYTHON SOURCE LINES 20-21

We'll now retrieve a Storm object for Hurricane Irma (2017), which will be used for this example code.

.. GENERATED FROM PYTHON SOURCE LINES 21-24

.. code-block:: default


    storm = basin.get_storm(('Irma',2017))








.. GENERATED FROM PYTHON SOURCE LINES 25-54

Now that we have a Storm object, there are several ways of retrieving recon data.

1. The first and easiest way is through the Storm object, which already contains an empty instance of ReconDataset stored as an attribute ``storm.recon``. This method will be highlighted in this script. Its methods can then be accessed as ``storm.recon.plot_summary()``, for example.

2. The second method is through retrieving an instance of ReconDataset, providing the storm object as an argument (e.g., ``recon_obj = recon.ReconDataset(storm)``). Its methods can then be accessed as ``recon_obj.plot_summary()``, for example.

3. The individual sub-classes (to be discussed later) can also be created individually of a ReconDataset object, providing the storm object as an argument (e.g., ``dropsondes = recon.dropsondes(storm)``).

For the rest of this script, we'll be using the first method.

-----------------
Recon Sub-Classes
-----------------
The three primary sub-classes of the recon module are:

.. list-table:: 
   :widths: 25 75
   :header-rows: 1

   * - Class
     - Description
   * - hdobs
     - Class containing all High Density Observations (HDOBs) for this Storm.
   * - dropsondes
     - Class containing all dropsondes for this Storm.
   * - vdms
     - Class containing all Vortex Data Messages (VDMs) for this Storm.

The following functionality is used to retrieve data for each sub-class. Each class has a ``to_pickle()`` method, which can be used to save the data once it's been read in as a local pickle file, which can be re-read in later.

.. GENERATED FROM PYTHON SOURCE LINES 54-67

.. code-block:: default


    storm.recon.get_vdms()
    #Save pickle file of VDM data (list of dictionaries)
    storm.recon.vdms.to_pickle(f'{storm.name}{storm.year}_vdms.pickle')

    storm.recon.get_dropsondes()
    #Save pickle file of Dropsonde data (list of dictionaries)
    storm.recon.dropsondes.to_pickle(f'{storm.name}{storm.year}_dropsondes.pickle')

    storm.recon.get_hdobs()
    #Save pickle file of HDOB data (Pandas dataframe)
    storm.recon.hdobs.to_pickle(f'{storm.name}{storm.year}_hdobs.pickle')





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    Searching through recon VDM files between 20170830 and 20170913 ...
    --> Completed reading in recon VDM files (17.0 seconds)
    Read 82 files
    Searching through recon dropsonde files between 201708290000 and 201709141200 ...
    --> Completed reading in recon dropsonde files (85.4 seconds)
    Read 460 files
    Searching through recon HDOB files between 20170829 and 20170914 ...
    --> Completed reading in recon HDOB files (325.7 seconds)
    Read 1491 files
    Unable to decode 0 files




.. GENERATED FROM PYTHON SOURCE LINES 68-74

These can be initialized again any time later with the saved pickle files:

>>> storm.recon.get_vdms(f'{storm.name}{storm.year}_vdms.pickle')
>>> storm.recon.get_dropsondes(f'{storm.name}{storm.year}_dropsondes.pickle')
>>> storm.recon.get_hdobs(f'{storm.name}{storm.year}_hdobs.pickle')


.. GENERATED FROM PYTHON SOURCE LINES 76-79

Visualizing ReconDataset
------------------------
Print the recon object to see a summary of the data in recon from the three objects:

.. GENERATED FROM PYTHON SOURCE LINES 79-82

.. code-block:: default


    storm.recon





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <tropycal.recon.hdobs>
    Dataset Summary:
        Storm:                        IRMA 2017
        Missions:                     31
        Time range:                   Sep-03 19:11 — Sep-10 22:37
        Max 30sec flight level wind:  162.0 knots
        Max 10sec flight level wind:  167.0 knots
        Max SFMR wind:                172.0 knots
        Min surface pressure:         909.1 hPa
        Source:                       National Hurricane Center (NHC)
    <tropycal.recon.dropsondes>
    Dataset Summary:
        Storm:                   IRMA 2017
        Missions:                31
        Dropsondes:              458
        Max 500m-avg wind:       170.0 knots
        Max 150m-avg wind:       163.0 knots
        Min sea level pressure:  915.0 hPa
        Source:                  National Hurricane Center (NHC)
    <tropycal.recon.vdms>
    Dataset Summary:
        Storm:                   IRMA 2017
        Missions:                23
        VDMs:                    82
        Min sea level pressure:  915.0 hPa
        Source:                  National Hurricane Center (NHC)



.. GENERATED FROM PYTHON SOURCE LINES 83-84

A summary of recon data for this storm can also be plotted:

.. GENERATED FROM PYTHON SOURCE LINES 84-87

.. code-block:: default


    storm.recon.plot_summary()




.. image-sg:: /examples/images/sphx_glr_recon_001.png
   :alt: Hurricane IRMA Recon summary, Start ... 19:11 UTC 03 Sep 2017 End ... 22:37 UTC 10 Sep 2017
   :srcset: /examples/images/sphx_glr_recon_001.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <GeoAxes: title={'left': 'Hurricane IRMA\nRecon summary', 'right': 'Start ... 19:11 UTC 03 Sep 2017\nEnd ... 22:37 UTC 10 Sep 2017'}>



.. GENERATED FROM PYTHON SOURCE LINES 88-89

Find the mission numbers that were active during a given time, within a distance (in km) from the storm:

.. GENERATED FROM PYTHON SOURCE LINES 89-92

.. code-block:: default


    storm.recon.find_mission(dt.datetime(2017,9,7,12), distance=200)





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    ['17']



.. GENERATED FROM PYTHON SOURCE LINES 93-94

Then plot a summary from just that mission:

.. GENERATED FROM PYTHON SOURCE LINES 94-97

.. code-block:: default


    storm.recon.plot_summary(mission=17)




.. image-sg:: /examples/images/sphx_glr_recon_002.png
   :alt: Hurricane IRMA Recon summary for mission 17, Start ... 08:54 UTC 07 Sep 2017 End ... 19:03 UTC 07 Sep 2017
   :srcset: /examples/images/sphx_glr_recon_002.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <GeoAxes: title={'left': 'Hurricane IRMA\nRecon summary for mission 17', 'right': 'Start ... 08:54 UTC 07 Sep 2017\nEnd ... 19:03 UTC 07 Sep 2017'}>



.. GENERATED FROM PYTHON SOURCE LINES 98-101

High Density Observations (HDOBs)
---------------------------------
The first class we'll be reviewing is the HDOBs class, which is the largest containing the most data given the high frequency of observations. Let's start by viewing the HDOB summary:

.. GENERATED FROM PYTHON SOURCE LINES 101-104

.. code-block:: default


    storm.recon.hdobs





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <tropycal.recon.hdobs>
    Dataset Summary:
        Storm:                        IRMA 2017
        Missions:                     31
        Time range:                   Sep-03 19:11 — Sep-10 22:37
        Max 30sec flight level wind:  162.0 knots
        Max 10sec flight level wind:  167.0 knots
        Max SFMR wind:                172.0 knots
        Min surface pressure:         909.1 hPa
        Source:                       National Hurricane Center (NHC)



.. GENERATED FROM PYTHON SOURCE LINES 105-106

And view HDOB Pandas DataFrame data:

.. GENERATED FROM PYTHON SOURCE LINES 106-109

.. code-block:: default


    storm.recon.hdobs.data






.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">
    <div>
    <style scoped>
        .dataframe tbody tr th:only-of-type {
            vertical-align: middle;
        }

        .dataframe tbody tr th {
            vertical-align: top;
        }

        .dataframe thead th {
            text-align: right;
        }
    </style>
    <table border="1" class="dataframe">
      <thead>
        <tr style="text-align: right;">
          <th></th>
          <th>index</th>
          <th>time</th>
          <th>lat</th>
          <th>lon</th>
          <th>plane_p</th>
          <th>plane_z</th>
          <th>p_sfc</th>
          <th>temp</th>
          <th>dwpt</th>
          <th>wdir</th>
          <th>wspd</th>
          <th>pkwnd</th>
          <th>sfmr</th>
          <th>rain</th>
          <th>flag</th>
          <th>mission</th>
          <th>mission_id</th>
          <th>xdist</th>
          <th>ydist</th>
          <th>distance</th>
          <th>iscenter</th>
        </tr>
      </thead>
      <tbody>
        <tr>
          <th>0</th>
          <td>0</td>
          <td>2017-09-03 19:11:30</td>
          <td>13.07</td>
          <td>-59.48</td>
          <td>1000.5</td>
          <td>112.0</td>
          <td>1013.4</td>
          <td>28.4</td>
          <td>21.4</td>
          <td>86.0</td>
          <td>5.0</td>
          <td>7.0</td>
          <td>NaN</td>
          <td>NaN</td>
          <td>[sfmr, rain]</td>
          <td>01</td>
          <td>NOAA2-0111A-IRMA</td>
          <td>-1064.102798</td>
          <td>-494.372330</td>
          <td>1173.336595</td>
          <td>0</td>
        </tr>
        <tr>
          <th>1</th>
          <td>1</td>
          <td>2017-09-03 19:12:00</td>
          <td>13.07</td>
          <td>-59.45</td>
          <td>980.3</td>
          <td>294.0</td>
          <td>1013.3</td>
          <td>26.0</td>
          <td>21.8</td>
          <td>66.0</td>
          <td>5.0</td>
          <td>6.0</td>
          <td>NaN</td>
          <td>NaN</td>
          <td>[sfmr, rain]</td>
          <td>01</td>
          <td>NOAA2-0111A-IRMA</td>
          <td>-1060.748502</td>
          <td>-494.302829</td>
          <td>1170.266069</td>
          <td>0</td>
        </tr>
        <tr>
          <th>2</th>
          <td>2</td>
          <td>2017-09-03 19:12:30</td>
          <td>13.08</td>
          <td>-59.43</td>
          <td>964.9</td>
          <td>435.0</td>
          <td>1013.8</td>
          <td>24.6</td>
          <td>21.6</td>
          <td>57.0</td>
          <td>4.0</td>
          <td>4.0</td>
          <td>NaN</td>
          <td>NaN</td>
          <td>[sfmr, rain]</td>
          <td>01</td>
          <td>NOAA2-0111A-IRMA</td>
          <td>-1058.454254</td>
          <td>-493.121318</td>
          <td>1167.687476</td>
          <td>0</td>
        </tr>
        <tr>
          <th>3</th>
          <td>3</td>
          <td>2017-09-03 19:13:00</td>
          <td>13.08</td>
          <td>-59.40</td>
          <td>942.8</td>
          <td>639.0</td>
          <td>1013.6</td>
          <td>23.4</td>
          <td>19.8</td>
          <td>39.0</td>
          <td>0.0</td>
          <td>3.0</td>
          <td>NaN</td>
          <td>NaN</td>
          <td>[sfmr, rain]</td>
          <td>01</td>
          <td>NOAA2-0111A-IRMA</td>
          <td>-1055.099988</td>
          <td>-493.051701</td>
          <td>1164.618377</td>
          <td>0</td>
        </tr>
        <tr>
          <th>4</th>
          <td>4</td>
          <td>2017-09-03 19:13:30</td>
          <td>13.10</td>
          <td>-59.37</td>
          <td>916.3</td>
          <td>887.0</td>
          <td>1013.4</td>
          <td>22.1</td>
          <td>18.1</td>
          <td>245.0</td>
          <td>6.0</td>
          <td>7.0</td>
          <td>NaN</td>
          <td>NaN</td>
          <td>[sfmr, rain]</td>
          <td>01</td>
          <td>NOAA2-0111A-IRMA</td>
          <td>-1051.745735</td>
          <td>-490.758124</td>
          <td>1160.608732</td>
          <td>0</td>
        </tr>
        <tr>
          <th>...</th>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
        </tr>
        <tr>
          <th>15</th>
          <td>15</td>
          <td>2017-09-10 22:35:30</td>
          <td>29.05</td>
          <td>-85.88</td>
          <td>410.3</td>
          <td>7353.0</td>
          <td>NaN</td>
          <td>-14.5</td>
          <td>-45.9</td>
          <td>351.0</td>
          <td>20.0</td>
          <td>21.0</td>
          <td>35.0</td>
          <td>0.0</td>
          <td>[]</td>
          <td>31</td>
          <td>AF307-3111A-IRMA</td>
          <td>-414.301548</td>
          <td>283.314691</td>
          <td>501.909341</td>
          <td>0</td>
        </tr>
        <tr>
          <th>16</th>
          <td>16</td>
          <td>2017-09-10 22:36:00</td>
          <td>29.07</td>
          <td>-85.93</td>
          <td>410.5</td>
          <td>7352.0</td>
          <td>NaN</td>
          <td>-14.5</td>
          <td>-45.5</td>
          <td>347.0</td>
          <td>21.0</td>
          <td>21.0</td>
          <td>36.0</td>
          <td>0.0</td>
          <td>[]</td>
          <td>31</td>
          <td>AF307-3111A-IRMA</td>
          <td>-419.299612</td>
          <td>285.347534</td>
          <td>507.183773</td>
          <td>0</td>
        </tr>
        <tr>
          <th>17</th>
          <td>17</td>
          <td>2017-09-10 22:36:30</td>
          <td>29.08</td>
          <td>-85.97</td>
          <td>410.4</td>
          <td>7356.0</td>
          <td>NaN</td>
          <td>-14.5</td>
          <td>-45.2</td>
          <td>346.0</td>
          <td>21.0</td>
          <td>21.0</td>
          <td>35.0</td>
          <td>0.0</td>
          <td>[sfmr, rain]</td>
          <td>31</td>
          <td>AF307-3111A-IRMA</td>
          <td>-423.302335</td>
          <td>286.268366</td>
          <td>511.013154</td>
          <td>0</td>
        </tr>
        <tr>
          <th>18</th>
          <td>18</td>
          <td>2017-09-10 22:37:00</td>
          <td>29.12</td>
          <td>-86.02</td>
          <td>410.5</td>
          <td>7350.0</td>
          <td>NaN</td>
          <td>-14.5</td>
          <td>-44.9</td>
          <td>350.0</td>
          <td>21.0</td>
          <td>22.0</td>
          <td>34.0</td>
          <td>0.0</td>
          <td>[]</td>
          <td>31</td>
          <td>AF307-3111A-IRMA</td>
          <td>-428.299609</td>
          <td>290.524992</td>
          <td>517.537753</td>
          <td>0</td>
        </tr>
        <tr>
          <th>19</th>
          <td>19</td>
          <td>2017-09-10 22:37:30</td>
          <td>29.13</td>
          <td>-86.07</td>
          <td>410.5</td>
          <td>7351.0</td>
          <td>NaN</td>
          <td>-14.5</td>
          <td>-44.7</td>
          <td>348.0</td>
          <td>22.0</td>
          <td>22.0</td>
          <td>34.0</td>
          <td>0.0</td>
          <td>[]</td>
          <td>31</td>
          <td>AF307-3111A-IRMA</td>
          <td>-433.296478</td>
          <td>291.445707</td>
          <td>522.193870</td>
          <td>0</td>
        </tr>
      </tbody>
    </table>
    <p>29694 rows × 21 columns</p>
    </div>
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 110-111

Plot a summary of the recon data for this storm, using peak wind gusts with a custom colormap:

.. GENERATED FROM PYTHON SOURCE LINES 111-114

.. code-block:: default


    storm.recon.hdobs.plot_points('pkwnd',prop={'cmap':{1:'dodgerblue',2:'gold',3:'firebrick'},'levels':np.arange(20,161,10)})




.. image-sg:: /examples/images/sphx_glr_recon_003.png
   :alt: Hurricane IRMA Recon: 10s FL wind speed (kt), Start ... 19:11 UTC 03 Sep 2017 End ... 22:37 UTC 10 Sep 2017
   :srcset: /examples/images/sphx_glr_recon_003.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <GeoAxes: title={'left': 'Hurricane IRMA\nRecon: 10s FL wind speed (kt)', 'right': 'Start ... 19:11 UTC 03 Sep 2017\nEnd ... 22:37 UTC 10 Sep 2017'}>



.. GENERATED FROM PYTHON SOURCE LINES 115-116

Plot a hovmoller from recon data interpolated to time and radius - note the eyewall replacement cycles:

.. GENERATED FROM PYTHON SOURCE LINES 116-119

.. code-block:: default


    storm.recon.hdobs.plot_hovmoller(varname='pkwnd',prop={'cmap':{1:'dodgerblue',2:'gold',3:'firebrick'},'levels':np.arange(20,161,10)})




.. image-sg:: /examples/images/sphx_glr_recon_004.png
   :alt: Hurricane IRMA Recon: 10s FL wind speed (kt), Start ... 06:00 UTC 04 Sep 2017 End ... 19:29 UTC 10 Sep 2017
   :srcset: /examples/images/sphx_glr_recon_004.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    --> Starting interpolation
    --> Completed interpolation (34.38 seconds)

    <Axes: title={'left': 'Hurricane IRMA\nRecon: 10s FL wind speed (kt)', 'right': 'Start ... 06:00 UTC 04 Sep 2017\nEnd ... 19:29 UTC 10 Sep 2017'}, xlabel='Radius (km)', ylabel='UTC Time (MM-DD HH)'>



.. GENERATED FROM PYTHON SOURCE LINES 120-121

Plot a map valid at 1200 UTC 6 September 2017 interpolated to time and space:

.. GENERATED FROM PYTHON SOURCE LINES 121-125

.. code-block:: default


    time = dt.datetime(2017,9,6,12)
    storm.recon.hdobs.plot_maps(time=time,varname='pkwnd',prop={'cmap':{1:'dodgerblue',2:'gold',3:'firebrick'},'levels':np.arange(20,161,10)})




.. image-sg:: /examples/images/sphx_glr_recon_005.png
   :alt: Hurricane IRMA Recon: 10s FL wind speed (kt), 12:00 UTC 06 Sep 2017
   :srcset: /examples/images/sphx_glr_recon_005.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    --> Starting interpolation
    2017-09-06 12:00:00
    --> Completed interpolation (11.83 seconds)

    <GeoAxes: title={'left': 'Hurricane IRMA\nRecon: 10s FL wind speed (kt)', 'right': '12:00 UTC 06 Sep 2017'}>



.. GENERATED FROM PYTHON SOURCE LINES 126-129

Dropsonde Data
--------------
Next we'll take a look at the dropsonde data for Hurricane Irma. First, let's take a look at the dropsonde summary:

.. GENERATED FROM PYTHON SOURCE LINES 129-132

.. code-block:: default


    storm.recon.dropsondes





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <tropycal.recon.dropsondes>
    Dataset Summary:
        Storm:                   IRMA 2017
        Missions:                31
        Dropsondes:              458
        Max 500m-avg wind:       170.0 knots
        Max 150m-avg wind:       163.0 knots
        Min sea level pressure:  915.0 hPa
        Source:                  National Hurricane Center (NHC)



.. GENERATED FROM PYTHON SOURCE LINES 133-134

Now use the ``sel`` function to subset to only dropsondes released in the eyewall:

.. GENERATED FROM PYTHON SOURCE LINES 134-137

.. code-block:: default


    storm.recon.dropsondes.sel(location='eyewall')





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <tropycal.recon.dropsondes>
    Dataset Summary:
        Storm:                   IRMA 2017
        Missions:                19
        Dropsondes:              73
        Max 500m-avg wind:       170.0 knots
        Max 150m-avg wind:       163.0 knots
        Min sea level pressure:  926.0 hPa
        Source:                  National Hurricane Center (NHC)



.. GENERATED FROM PYTHON SOURCE LINES 138-139

We can view data from one of the dropsondes by using the ``isel`` method to select a dropsonde number:

.. GENERATED FROM PYTHON SOURCE LINES 139-142

.. code-block:: default


    storm.recon.dropsondes.sel(location='eyewall').isel(23).data





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    [{'lat': 17.6, 'lon': -61.3, 'slp': 948.0, 'TOPlat': 17.61, 'TOPlon': -61.32, 'TOPtime': Timestamp('2017-09-06 03:07:00'), 'BOTTOMlat': 17.49, 'BOTTOMlon': -61.4, 'BOTTOMtime': Timestamp('2017-09-06 03:10:00'), 'MBLdir': nan, 'MBLspd': nan, 'DLMdir': 25.0, 'DLMspd': 139.0, 'WL150dir': 345.0, 'WL150spd': 133.0, 'top': 696.0, 'LSThgt': 13.0, 'software': 'AEV 33304', 'levels':       pres     hgt  temp  dwpt   wdir   wspd
    0   1000.0     NaN   NaN   NaN    NaN    NaN
    1    948.0     0.0  24.6  23.9    NaN    NaN
    2    947.0     NaN   NaN   NaN  340.0  128.0
    3    928.0     NaN   NaN   NaN  355.0  137.0
    4    925.0   220.0  23.8  23.0  355.0  140.0
    5    910.0     NaN   NaN   NaN  355.0  154.0
    6    901.0     NaN   NaN   NaN    0.0  169.0
    7    897.0     NaN   NaN   NaN    0.0  163.0
    8    891.0     NaN   NaN   NaN    0.0  139.0
    9    879.0     NaN   NaN   NaN   10.0  142.0
    10   875.0     NaN   NaN   NaN   10.0  138.0
    11   869.0     NaN   NaN   NaN   15.0  150.0
    12   863.0     NaN   NaN   NaN   20.0  156.0
    13   856.0     NaN   NaN   NaN   20.0  174.0
    14   850.0  1960.0  20.6  19.7   20.0  179.0
    15   832.0     NaN   NaN   NaN   25.0  153.0
    16   812.0     NaN   NaN   NaN   30.0  163.0
    17   804.0     NaN   NaN   NaN   30.0  143.0
    18   756.0     NaN   NaN   NaN   40.0  124.0
    19   732.0     NaN  16.4  14.1    NaN    NaN
    20   700.0  2626.0  11.0   7.6   45.0  145.0
    21   696.0     NaN  10.2   6.6   45.0  146.0, 'mission': '11', 'stormname': 'IRMA', 'obsnum': 21, 'mission_id': 'AF305-1111A-IRMA', 'TOPxdist': -17.607565905503332, 'TOPydist': 18.667114071950596, 'TOPdistance': 25.66101176672228, 'BOTTOMxdist': -24.822341252146206, 'BOTTOMydist': 4.665134677347142, 'BOTTOMdistance': 25.256921958065, 'location': 'EYEWALL', 'octant': 'NW'}]



.. GENERATED FROM PYTHON SOURCE LINES 143-144

Select one of the eyewall dropsondes and plot the Skew-T:

.. GENERATED FROM PYTHON SOURCE LINES 144-147

.. code-block:: default


    storm.recon.dropsondes.sel(location='eyewall').isel(23).plot_skewt()





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <Figure size 1700x1100 with 4 Axes>



.. GENERATED FROM PYTHON SOURCE LINES 148-149

Plot a map of dropsonde points colored by 850mb temperature:

.. GENERATED FROM PYTHON SOURCE LINES 149-152

.. code-block:: default


    storm.recon.dropsondes.plot_points('temp',level=850,prop={'cmap':{1:'dodgerblue',2:'gold',3:'firebrick'},'ms':20})




.. image-sg:: /examples/images/sphx_glr_recon_006.png
   :alt: Hurricane IRMA Recon: 850hPa temperature (deg C), Start ... 21:05 UTC 03 Sep 2017 End ... 21:33 UTC 10 Sep 2017
   :srcset: /examples/images/sphx_glr_recon_006.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <GeoAxes: title={'left': 'Hurricane IRMA\nRecon: 850hPa temperature (deg C)', 'right': 'Start ... 21:05 UTC 03 Sep 2017\nEnd ... 21:33 UTC 10 Sep 2017'}>



.. GENERATED FROM PYTHON SOURCE LINES 153-154

Plot a map of only upper-air dropsondes released at 300mb and above, colored by wind speed at 300mb:

.. GENERATED FROM PYTHON SOURCE LINES 154-157

.. code-block:: default


    storm.recon.dropsondes.sel(top=(None,300)).plot_points('wspd',level=300,prop={'cmap':{1:'dodgerblue',2:'gold',3:'firebrick'},'ms':20})




.. image-sg:: /examples/images/sphx_glr_recon_007.png
   :alt: Hurricane IRMA Recon: 300hPa wind speed (kt), Start ... 17:58 UTC 04 Sep 2017 End ... 01:04 UTC 10 Sep 2017
   :srcset: /examples/images/sphx_glr_recon_007.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <GeoAxes: title={'left': 'Hurricane IRMA\nRecon: 300hPa wind speed (kt)', 'right': 'Start ... 17:58 UTC 04 Sep 2017\nEnd ... 01:04 UTC 10 Sep 2017'}>



.. GENERATED FROM PYTHON SOURCE LINES 158-161

Vortex Data Messages (VDMs)
---------------------------
The last class we'll look at is the VDM class. Let's start off by viewing a summary of VDM data for Hurricane Irma:

.. GENERATED FROM PYTHON SOURCE LINES 161-164

.. code-block:: default


    storm.recon.vdms





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <tropycal.recon.vdms>
    Dataset Summary:
        Storm:                   IRMA 2017
        Missions:                23
        VDMs:                    82
        Min sea level pressure:  915.0 hPa
        Source:                  National Hurricane Center (NHC)



.. GENERATED FROM PYTHON SOURCE LINES 165-166

Let's look at decoded VDMs for a specific pass:

.. GENERATED FROM PYTHON SOURCE LINES 166-169

.. code-block:: default


    storm.recon.vdms.isel(10).data





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    [{'time': datetime.datetime(2017, 9, 4, 23, 7, 5), 'lat': 16.68, 'lon': -54.87, 'Standard Level (hPa)': nan, 'Minimum Height at Standard Level (m)': nan, 'Estimated Maximum Surface Wind Inbound (kt)': 97.0, 'Dropsonde Surface Wind Speed at Center (kt)': 13.0, 'Dropsonde Surface Wind Direction at Center (deg)': 263.0, 'Maximum Flight Level Wind Inbound': '353 deg 110 kt', 'Location of the Maximum Flight Level Wind Inbound': '264 deg 15 nm', 'Minimum Sea Level Pressure (hPa)': 945.0, 'Maximum Flight Level Temp Outside Eye (C)': 13.0, 'Maximum Flight Level Temp Inside Eye (C)': 19.0, 'Dew Point Inside Eye (C)': 18.0, 'Eye character': 'closed', 'Eye Shape': 'circular', 'Eye Diameter (nmi)': 25.0, 'Aircraft': 'NOAA2', 'mission': '05', 'Remarks': ' MAX FL WIND 130 KT 001 / 22 NM 21:45:54Z MAX FL TEMP 19 C 108 / 6 NM FROM FL CNTR CNTR DROPSONDE SFC WIND 130 / 06 KTS  ', 'mission_id': 'NOAA2-0511A-IRMA'}]



.. GENERATED FROM PYTHON SOURCE LINES 170-171

Plot a map of VDM center location, colored by minimum pressure (default):

.. GENERATED FROM PYTHON SOURCE LINES 171-173

.. code-block:: default


    storm.recon.vdms.plot_points(prop={'cmap':{3:'dodgerblue',2:'gold',1:'firebrick'},'ms':40})



.. image-sg:: /examples/images/sphx_glr_recon_008.png
   :alt: Hurricane IRMA Recon: Minimum Sea Level Pressure (hPa), Start ... 21:26 UTC 03 Sep 2017 End ... 20:38 UTC 10 Sep 2017
   :srcset: /examples/images/sphx_glr_recon_008.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <GeoAxes: title={'left': 'Hurricane IRMA\nRecon: Minimum Sea Level Pressure (hPa)', 'right': 'Start ... 21:26 UTC 03 Sep 2017\nEnd ... 20:38 UTC 10 Sep 2017'}>




.. rst-class:: sphx-glr-timing

   **Total running time of the script:** ( 8 minutes  37.605 seconds)


.. _sphx_glr_download_examples_recon.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example




    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: recon.py <recon.py>`

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: recon.ipynb <recon.ipynb>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_