{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "**Brian Blaylock** \n", "*December 18, 2020*\n", "\n", "# 🌎 ABI Geostationary **field-of-view** in Cartopy\n", "\n", "- [GOES-R Series Data Book](https://www.goes-r.gov/downloads/resources/documents/GOES-RSeriesDataBook.pdf)\n", " - Full Disk Coverage: Defined as a 17.4 degree diameter circle (pg 48)\n" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "WARNING! No module named 'dask'\n", "Without dask, you cannot use dask for multiprocessing.\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/p/home/blaylock/BB_python/Carpenter_Workshop/toolbox/cartopy_tools.py:33: UserWarning: No module named 'geopandas' Without geopandas, you cannot subset some NaturalEarthFeatures, like \"Major Highways\" from roads.\n", " warnings.warn(\n" ] } ], "source": [ "from datetime import datetime\n", "import cartopy.crs as ccrs\n", "import matplotlib.pyplot as plt\n", "from shapely.geometry import Point\n", "import numpy as np\n", "\n", "from goes2go.data import goes_latest, goes_nearesttime\n", "\n", "from BB_maps.add_to_cartopy import common_features, pc" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " _______________________________\n", " | Satellite: noaa-goes16 |\n", " | Product: ABI-L2-MCMIPC |\n", " | Domain: C |\n", "📦 Finished downloading [1] files to [/p/cwfs/blaylock/data/noaa-goes16/ABI-L2-MCMIPC]. \n", "📚 Finished reading [1] files into xarray.Dataset. \n", " _______________________________\n", " | Satellite: noaa-goes16 |\n", " | Product: ABI-L2-MCMIPC |\n", " | Domain: C |\n", "📦 Finished downloading [1] files to [/p/cwfs/blaylock/data/noaa-goes16/ABI-L2-MCMIPC]. \n", "📚 Finished reading [1] files into xarray.Dataset. \n", " _______________________________\n", " | Satellite: noaa-goes17 |\n", " | Product: ABI-L2-MCMIPC |\n", " | Domain: C |\n", "📦 Finished downloading [1] files to [/p/cwfs/blaylock/data/noaa-goes17/ABI-L2-MCMIPC]. \n", "📚 Finished reading [1] files into xarray.Dataset. \n" ] } ], "source": [ "# Testing Phase\n", "gcentral = goes_nearesttime(datetime(2017, 10, 5), satellite=\"G16\", product=\"ABI\")\n", "\n", "# GOES East\n", "g16 = goes_latest(satellite=\"G16\", product=\"ABI\")\n", "\n", "# GOES West\n", "g17 = goes_latest(satellite=\"G17\", product=\"ABI\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Get the data's coordinate reference system and field of view with the `goes2go.tools.field_of_view` function\n", "This returns both the field of view as a shapely polygon and the cartopy coordinate reference system." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "
<cartopy.crs.Geostationary object at 0x2b5eb637b400>" ], "text/plain": [ "