notebooks/00_geosat_explaining_steps.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<img src=\"images/logos/ncar_pythia.png\" alt=\"pythia ncar\" />\n",
    "\n",
    "**This notebook is developed during the pythia cook-off at NCAR Mesa-Lab Boulder Colorado, June 12-14, 2024**\n",
    "\n",
    "Participants in the workshop event have the chance to practice collaborative problem-solving and hands-on learning in the field of Python programming.\n",
    "\n",
    "This notebook is part of the **Breakout Topic: Geostationary on AWS**, lead by **Jorge Humberto Bravo Mendez** jbravo2@stevens.edu, from Stevens Institute of Technology "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Advanced Baseline Imager (ABI) data with Satpy\n",
    "\n",
    "\n",
    "Using Satpy to read and  Advanced Baseline Imager (ABI) data from GOES-R satellites. Here's a step-by-step guide:"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Overview\n",
    "It takes a methodical approach to efficiently access, evaluate, and visualize satellite imagery in order to create a Jupyter Notebook to read Advanced Baseline Imager (ABI) data with Satpy. The notebook starts by configuring the environment and importing the modules before gaining access to ABI data from AWS. A Satpy Scene object is initialized by using the Python library Satpy, which is designed for earth-observing satellite instruments. It specifies the ABI data file(s) and the relevant reader. The desired ABI channels or products are then loaded into the Scene object. Then, the loaded data is accessed and visualized using the built-in visualization tools of Matplotlib or Satpy. To ensure that the analysis is clear and repeatable, the notebook places a strong emphasis on comprehensive documentation and explanation. Users can learn more about the analysis of satellite imagery with this methodical approach, opening up a world of possibilities for everything from environmental research to weather monitoring."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Prerequisites\n",
    "\n",
    "<div class=\"admonition alert alert-info\">\n",
    "    <p class=\"admonition-title\" style=\"font-weight:bold\">Info</p>\n",
    "    Some prerequisites to understanding Satpy effectively:\n",
    "</div>\n",
    "\n",
    "\n",
    "1. **Basic Python:** Since Satpy is a Python library, knowledge of the Python programming language is required. It is essential to comprehend ideas such as variables, data types, loops, functions, and libraries/modules.\n",
    "\n",
    "1. **Knowledge of Satellite Data:** Although not required, it can be helpful to have some knowledge of satellite data and its uses. This entails being familiar with satellite instruments, data formats, and the different kinds of data that satellites gather.\n",
    "\n",
    "1. **Concepts of Data Processing:** It is beneficial to have a basic understanding of concepts related to reading, modifying, and visualizing data. Understanding these ideas will help in making the most of Satpy, as it offers features for reading, processing, and analyzing satellite data.\n",
    "\n",
    "1. **Python Libraries:** It can be helpful to have some knowledge of other Python libraries that are frequently used in scientific computing, such as NumPy, Matplotlib, and Pandas. These libraries are frequently integrated with Satpy for tasks involving data manipulation and visualization."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<div class=\"admonition alert alert-warning\">\n",
    "    <p class=\"admonition-title\" style=\"font-weight:bold\">Warning</p>\n",
    "    Before starting, we must obtain data\n",
    "</div>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "### List of files to download for testing\n",
    "urls2dwn = [ \n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C01_G16_s20210992350171_e20210992359479_c20210992359525.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C02_G16_s20210992350171_e20210992359479_c20210992359522.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C03_G16_s20210992350171_e20210992359479_c20210992359530.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C04_G16_s20210992350171_e20210992359479_c20210992359514.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C05_G16_s20210992350171_e20210992359479_c20210992359533.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C06_G16_s20210992350171_e20210992359484_c20210992359528.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C07_G16_s20210992350171_e20210992359490_c20210992359539.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C08_G16_s20210992350171_e20210992359479_c20210992359535.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C09_G16_s20210992350171_e20210992359484_c20210992359551.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C10_G16_s20210992350171_e20210992359490_c20210992359536.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C11_G16_s20210992350171_e20210992359479_c20210992359545.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C12_G16_s20210992350171_e20210992359484_c20210992359539.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C13_G16_s20210992350171_e20210992359490_c20210992359555.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C14_G16_s20210992350171_e20210992359479_c20210992359547.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C15_G16_s20210992350171_e20210992359484_c20210992359557.nc',\n",
    "  'https://noaa-goes16.s3.amazonaws.com/ABI-L1b-RadF/2021/099/23/OR_ABI-L1b-RadF-M6C16_G16_s20210992350171_e20210992359490_c20210992359542.nc']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "File 'OR_ABI-L1b-RadF-M6C01_G16_s20210992350171_e20210992359479_c20210992359525.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C02_G16_s20210992350171_e20210992359479_c20210992359522.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C03_G16_s20210992350171_e20210992359479_c20210992359530.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C04_G16_s20210992350171_e20210992359479_c20210992359514.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C05_G16_s20210992350171_e20210992359479_c20210992359533.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C06_G16_s20210992350171_e20210992359484_c20210992359528.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C07_G16_s20210992350171_e20210992359490_c20210992359539.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C08_G16_s20210992350171_e20210992359479_c20210992359535.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C09_G16_s20210992350171_e20210992359484_c20210992359551.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C10_G16_s20210992350171_e20210992359490_c20210992359536.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C11_G16_s20210992350171_e20210992359479_c20210992359545.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C12_G16_s20210992350171_e20210992359484_c20210992359539.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C13_G16_s20210992350171_e20210992359490_c20210992359555.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C14_G16_s20210992350171_e20210992359479_c20210992359547.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C15_G16_s20210992350171_e20210992359484_c20210992359557.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n",
      "File 'OR_ABI-L1b-RadF-M6C16_G16_s20210992350171_e20210992359490_c20210992359542.nc' downloaded and saved to 'Output_data/ABI-L1b-RadF/s20210992350171'.\n"
     ]
    }
   ],
   "source": [
    "import requests\n",
    "import os\n",
    "\n",
    "# Specify the local directory where you want to save the files.\n",
    "# local_directory = input(\"Enter the path to the download folder: \")\n",
    "local_directory = \"Output_data/ABI-L1b-RadF/s20210992350171\"\n",
    "# Ensure that the local directory exists; create it if it doesn't.\n",
    "os.makedirs(local_directory, exist_ok=True)\n",
    "\n",
    "# Iterate through the URLs and download files.\n",
    "for urld in urls2dwn:\n",
    "    # Extract the filename from the URL.\n",
    "    ntw = urld.split('/')[-1]\n",
    "    \n",
    "    # Construct the complete path to save the file in the local directory.\n",
    "    file_path = os.path.join(local_directory, ntw)\n",
    "    \n",
    "    # Send an HTTP GET request to the URL.\n",
    "    resp = requests.get(urld)\n",
    "    \n",
    "    # Check if the response is successful (status code 200).\n",
    "    if resp.status_code == 200:\n",
    "        # Write the content to the file in binary mode.\n",
    "        with open(file_path, \"wb\") as file:\n",
    "            file.write(resp.content)\n",
    "        print(f\"File '{ntw}' downloaded and saved to '{local_directory}'.\")\n",
    "    else:\n",
    "        print(f\"Failed to download '{ntw}' from the URL: {urld}\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Imports"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# To stop extraneous warning messages from clogging the notebook output, \n",
    "# import the warnings module and set it to ignore all warnings.\n",
    "import warnings\n",
    "warnings.filterwarnings('ignore')\n",
    "\n",
    "# Scene is a representation of satellite data that enables functions such as data reading, \n",
    "# data resampling, data compositing, and data saving.\n",
    "from satpy.scene import Scene\n",
    "# This line imports the debug_on function from the satpy.utils module, which is used to enable detailed \n",
    "# debug logging in Satpy. This can be helpful for troubleshooting and understanding the internal workings of Satpy processes.\n",
    "from satpy.utils import debug_on\n",
    "\n",
    "# Importing datetime for date and time operations. Processing time-stamped satellite data \n",
    "# is necessary for accurate time-based operations and analysis.\n",
    "from datetime import datetime\n",
    "\n",
    "# Using glob allows for efficient and flexible file handling, especially when dealing with \n",
    "# large datasets typically generated by geostationary satellites.\n",
    "from glob import glob"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# The magic command '%matplotlib inline' configures the Jupyter Notebook to display plots directly below the code cells.\n",
    "# This setting is essential for interactive data visualization, especially when plotting data from geostationary satellites.\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Starting to create satpy scenes"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<div class=\"admonition alert alert-danger\">\n",
    "    <p class=\"admonition-title\" style=\"font-weight:bold\">Attention</p>\n",
    "    SatPy requires the original file names at all times!<p>\n",
    "    Therefore, don't alter them when you save the data locally. SatPy won't be able to open the files if this is the case.\n",
    "</div>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['./Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C08_G16_s20210992350171_e20210992359479_c20210992359535.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C14_G16_s20210992350171_e20210992359479_c20210992359547.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C01_G16_s20210992350171_e20210992359479_c20210992359525.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C07_G16_s20210992350171_e20210992359490_c20210992359539.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C06_G16_s20210992350171_e20210992359484_c20210992359528.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C02_G16_s20210992350171_e20210992359479_c20210992359522.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C03_G16_s20210992350171_e20210992359479_c20210992359530.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C16_G16_s20210992350171_e20210992359490_c20210992359542.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C15_G16_s20210992350171_e20210992359484_c20210992359557.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C04_G16_s20210992350171_e20210992359479_c20210992359514.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C10_G16_s20210992350171_e20210992359490_c20210992359536.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C05_G16_s20210992350171_e20210992359479_c20210992359533.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C12_G16_s20210992350171_e20210992359484_c20210992359539.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C11_G16_s20210992350171_e20210992359479_c20210992359545.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C09_G16_s20210992350171_e20210992359484_c20210992359551.nc',\n",
       " './Output_data/ABI-L1b-RadF/s20210992350171/OR_ABI-L1b-RadF-M6C13_G16_s20210992350171_e20210992359490_c20210992359555.nc']"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sat_files = glob(\"./Output_data/ABI-L1b-RadF/s20210992350171/*\")\n",
    "sat_files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "['C01', 'C02', 'C03', 'C04', 'C05', 'C06', 'C07', 'C08', 'C09', 'C10', 'C11', 'C12', 'C13', 'C14', 'C15', 'C16']\n"
     ]
    }
   ],
   "source": [
    "# Create a Scene object named 'scn' by providing a list of filenames obtained from 'fGRl1b'.\n",
    "# The Scene object is a core part of the Satpy library, which organizes and manages multiple data sources typically from satellite observations.\n",
    "scn = Scene(filenames = sat_files, reader = 'abi_l1b')\n",
    "\n",
    "# Retrieve the names of all datasets available in the 'scn' object, which represent various satellite data channels.\n",
    "# This information is crucial for understanding what types of data are available for analysis, such as different spectral bands and derived products.\n",
    "dataset_names = scn.all_dataset_names()\n",
    "\n",
    "# Output the list of dataset names. This is useful for educational purposes to show students the variety of data that can be processed\n",
    "# and to select specific datasets for further analysis in practical exercises.\n",
    "print(dataset_names)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Load multiple datasets using a list comprehension to generate dataset names.\n",
    "# List comprehensions provide a concise way to create lists based on existing lists or ranges.\n",
    "# In this case, we generate names for the datasets 'C01' to 'C16', which are typical channel identifiers in satellite data.\n",
    "scn.load([f'C{x:02d}' for x in range(1, 17)])\n",
    "\n",
    "# Explanation of the list comprehension:\n",
    "# [f'C{x:02d}' for x in range(1, 17)] creates a list of strings from 'C01' to 'C16'.\n",
    "# 'f' before the string starts an f-string, allowing us to insert variables directly into the string.\n",
    "# '{x:02d}' formats the number 'x' as a two-digit decimal, padding with zeros if necessary.\n",
    "\n",
    "# The 'scn.load' function is then used to load these specific datasets into the Scene object.\n",
    "# Loading multiple channels like this is common in the analysis of satellite imagery,\n",
    "# where each channel can represent different spectral bands and contain different types of environmental information."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "['24h_microphysics', 'airmass', 'ash', 'cimss_cloud_type', 'cimss_cloud_type_raw', 'cimss_green', 'cimss_green_sunz', 'cimss_green_sunz_rayleigh', 'cimss_true_color', 'cimss_true_color_sunz', 'cimss_true_color_sunz_rayleigh', 'cira_day_convection', 'cira_fire_temperature', 'cloud_phase', 'cloud_phase_distinction', 'cloud_phase_distinction_raw', 'cloud_phase_raw', 'cloudtop', 'color_infrared', 'colorized_ir_clouds', 'convection', 'day_microphysics', 'day_microphysics_abi', 'day_microphysics_eum', 'dust', 'fire_temperature_awips', 'fog', 'geo_color', 'geo_color_background_with_low_clouds', 'geo_color_high_clouds', 'geo_color_low_clouds', 'geo_color_night', 'green', 'green_crefl', 'green_nocorr', 'green_raw', 'green_snow', 'highlight_C14', 'ir108_3d', 'ir_cloud_day', 'land_cloud', 'land_cloud_fire', 'natural_color', 'natural_color_nocorr', 'natural_color_raw', 'natural_color_raw_with_night_ir', 'night_fog', 'night_ir_alpha', 'night_ir_with_background', 'night_ir_with_background_hires', 'night_microphysics', 'night_microphysics_eum', 'overview', 'overview_raw', 'rocket_plume_day', 'rocket_plume_night', 'snow', 'snow_fog', 'so2', 'tropical_airmass', 'true_color', 'true_color_crefl', 'true_color_nocorr', 'true_color_raw', 'true_color_reproduction', 'true_color_reproduction_corr', 'true_color_reproduction_uncorr', 'true_color_with_night_fires', 'true_color_with_night_fires_nocorr', 'true_color_with_night_ir', 'true_color_with_night_ir_hires', 'water_vapors1', 'water_vapors2']\n"
     ]
    }
   ],
   "source": [
    "# The method 'available_composite_names' is called on the 'scn' object.\n",
    "# This method retrieves a list of all the composite images that can be created using the loaded data channels.\n",
    "# Composite images are made by combining multiple data channels to enhance the visualization and interpretation of satellite data.\n",
    "# This feature is particularly useful in the study of geostationary satellites, as it allows for more detailed and informative visual representations of atmospheric phenomena.\n",
    "\n",
    "# Retrieve and print the list of available composite names, providing a crucial insight into the data visualization capabilities.\n",
    "print(scn.available_composite_names())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Assign the dataset name 'airmass' to the variable 'rgb_im'.\n",
    "# This variable naming provides clarity when referencing the dataset in multiple places, \n",
    "# ensuring consistency and reducing the likelihood of errors in dataset identification.\n",
    "rgb_im = 'airmass'\n",
    "\n",
    "# Load the dataset named 'airmass' using the 'scn.load' method.\n",
    "# The 'airmass' composite is particularly useful in meteorology as it combines several spectral bands to highlight features like dust, ash, and water vapor,\n",
    "# making it easier to analyze atmospheric conditions from geostationary satellite data.\n",
    "scn.load([rgb_im])\n",
    "\n",
    "# Display the loaded 'airmass' dataset using 'scn.show'.\n",
    "# This method visualizes the specified dataset, allowing students to see the practical application of satellite data composites \n",
    "# and understand their relevance in real-world atmospheric monitoring and analysis.\n",
    "scn.show(rgb_im)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body[data-theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block !important;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-index-preview {\n",
       "  grid-column: 2 / 5;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data,\n",
       ".xr-index-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data,\n",
       ".xr-index-data-in:checked ~ .xr-index-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-index-name div,\n",
       ".xr-index-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data,\n",
       ".xr-index-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2,\n",
       ".xr-no-icon {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.DataArray &#x27;where-303720c4e2e42379319fb14493c19145&#x27; (bands: 3, y: 5424,\n",
       "                                                            x: 5424)&gt;\n",
       "dask.array&lt;where, shape=(3, 5424, 5424), dtype=float64, chunksize=(1, 4096, 4096), chunktype=numpy.ndarray&gt;\n",
       "Coordinates:\n",
       "    crs      object PROJCRS[&quot;unknown&quot;,BASEGEOGCRS[&quot;unknown&quot;,DATUM[&quot;unknown&quot;,E...\n",
       "  * y        (y) float64 5.434e+06 5.432e+06 5.43e+06 ... -5.432e+06 -5.434e+06\n",
       "  * x        (x) float64 -5.434e+06 -5.432e+06 -5.43e+06 ... 5.432e+06 5.434e+06\n",
       "  * bands    (bands) &lt;U1 &#x27;R&#x27; &#x27;G&#x27; &#x27;B&#x27;\n",
       "Attributes: (12/29)\n",
       "    long_name:               Brightness Temperature\n",
       "    end_time:                2021-04-09 23:59:48.350000\n",
       "    area:                    Area ID: GOES-East\\nDescription: 2km at nadir\\nP...\n",
       "    orbital_parameters:      {&#x27;projection_longitude&#x27;: -75.0, &#x27;projection_lati...\n",
       "    observation_type:        Rad\n",
       "    reader:                  abi_l1b\n",
       "    ...                      ...\n",
       "    optional_datasets:       []\n",
       "    name:                    airmass\n",
       "    _satpy_id:               DataID(name=&#x27;airmass&#x27;, resolution=2000)\n",
       "    prerequisites:           [DataQuery(name=&#x27;_airmass_dep_0&#x27;), DataQuery(nam...\n",
       "    optional_prerequisites:  []\n",
       "    mode:                    RGB</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'>'where-303720c4e2e42379319fb14493c19145'</div><ul class='xr-dim-list'><li><span class='xr-has-index'>bands</span>: 3</li><li><span class='xr-has-index'>y</span>: 5424</li><li><span class='xr-has-index'>x</span>: 5424</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-c8618b10-4ed5-4c9e-8cd8-86b547d3a53d' class='xr-array-in' type='checkbox' checked><label for='section-c8618b10-4ed5-4c9e-8cd8-86b547d3a53d' title='Show/hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-array-preview xr-preview'><span>dask.array&lt;chunksize=(1, 4096, 4096), meta=np.ndarray&gt;</span></div><div class='xr-array-data'><table>\n",
       "    <tr>\n",
       "        <td>\n",
       "            <table style=\"border-collapse: collapse;\">\n",
       "                <thead>\n",
       "                    <tr>\n",
       "                        <td> </td>\n",
       "                        <th> Array </th>\n",
       "                        <th> Chunk </th>\n",
       "                    </tr>\n",
       "                </thead>\n",
       "                <tbody>\n",
       "                    \n",
       "                    <tr>\n",
       "                        <th> Bytes </th>\n",
       "                        <td> 673.37 MiB </td>\n",
       "                        <td> 128.00 MiB </td>\n",
       "                    </tr>\n",
       "                    \n",
       "                    <tr>\n",
       "                        <th> Shape </th>\n",
       "                        <td> (3, 5424, 5424) </td>\n",
       "                        <td> (1, 4096, 4096) </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <th> Dask graph </th>\n",
       "                        <td colspan=\"2\"> 12 chunks in 67 graph layers </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <th> Data type </th>\n",
       "                        <td colspan=\"2\"> float64 numpy.ndarray </td>\n",
       "                    </tr>\n",
       "                </tbody>\n",
       "            </table>\n",
       "        </td>\n",
       "        <td>\n",
       "        <svg width=\"194\" height=\"184\" style=\"stroke:rgb(0,0,0);stroke-width:1\" >\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"10\" y1=\"0\" x2=\"24\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"10\" y1=\"90\" x2=\"24\" y2=\"105\" />\n",
       "  <line x1=\"10\" y1=\"120\" x2=\"24\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"10\" y1=\"0\" x2=\"10\" y2=\"120\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"14\" y1=\"4\" x2=\"14\" y2=\"124\" />\n",
       "  <line x1=\"19\" y1=\"9\" x2=\"19\" y2=\"129\" />\n",
       "  <line x1=\"24\" y1=\"14\" x2=\"24\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"10.0,0.0 24.9485979497544,14.948597949754403 24.9485979497544,134.9485979497544 10.0,120.0\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"10\" y1=\"0\" x2=\"130\" y2=\"0\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"14\" y1=\"4\" x2=\"134\" y2=\"4\" />\n",
       "  <line x1=\"19\" y1=\"9\" x2=\"139\" y2=\"9\" />\n",
       "  <line x1=\"24\" y1=\"14\" x2=\"144\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"10\" y1=\"0\" x2=\"24\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"100\" y1=\"0\" x2=\"115\" y2=\"14\" />\n",
       "  <line x1=\"130\" y1=\"0\" x2=\"144\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"10.0,0.0 130.0,0.0 144.9485979497544,14.948597949754403 24.9485979497544,14.948597949754403\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"24\" y1=\"14\" x2=\"144\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"24\" y1=\"105\" x2=\"144\" y2=\"105\" />\n",
       "  <line x1=\"24\" y1=\"134\" x2=\"144\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"24\" y1=\"14\" x2=\"24\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"115\" y1=\"14\" x2=\"115\" y2=\"134\" />\n",
       "  <line x1=\"144\" y1=\"14\" x2=\"144\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"24.9485979497544,14.948597949754403 144.9485979497544,14.948597949754403 144.9485979497544,134.9485979497544 24.9485979497544,134.9485979497544\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Text -->\n",
       "  <text x=\"84.948598\" y=\"154.948598\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" >5424</text>\n",
       "  <text x=\"164.948598\" y=\"74.948598\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(-90,164.948598,74.948598)\">5424</text>\n",
       "  <text x=\"7.474299\" y=\"147.474299\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(45,7.474299,147.474299)\">3</text>\n",
       "</svg>\n",
       "        </td>\n",
       "    </tr>\n",
       "</table></div></div></li><li class='xr-section-item'><input id='section-b5500fc0-2d99-4de3-b0a7-c814ade9f536' class='xr-section-summary-in' type='checkbox'  checked><label for='section-b5500fc0-2d99-4de3-b0a7-c814ade9f536' class='xr-section-summary' >Coordinates: <span>(4)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span>crs</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>PROJCRS[&quot;unknown&quot;,BASEGEOGCRS[&quot;u...</div><input id='attrs-aa4b9711-0a51-44b7-9a1d-226f11a22994' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-aa4b9711-0a51-44b7-9a1d-226f11a22994' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e4be1acd-d1e9-450b-93d9-c5279cdeb0e6' class='xr-var-data-in' type='checkbox'><label for='data-e4be1acd-d1e9-450b-93d9-c5279cdeb0e6' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array(&lt;Derived Projected CRS: PROJCRS[&quot;unknown&quot;,BASEGEOGCRS[&quot;unknown&quot;,DATUM[&quot;unk ...&gt;\n",
       "Name: unknown\n",
       "Axis Info [cartesian]:\n",
       "- E[east]: Easting (metre)\n",
       "- N[north]: Northing (metre)\n",
       "Area of Use:\n",
       "- undefined\n",
       "Coordinate Operation:\n",
       "- name: unknown\n",
       "- method: Geostationary Satellite (Sweep X)\n",
       "Datum: unknown\n",
       "- Ellipsoid: GRS 1980\n",
       "- Prime Meridian: Greenwich\n",
       ", dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>y</span></div><div class='xr-var-dims'>(y)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>5.434e+06 5.432e+06 ... -5.434e+06</div><input id='attrs-34181645-af48-4881-952b-416bc23d9365' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-34181645-af48-4881-952b-416bc23d9365' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-017f98aa-bbb2-471c-b0b4-6da61c4907cb' class='xr-var-data-in' type='checkbox'><label for='data-017f98aa-bbb2-471c-b0b4-6da61c4907cb' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>meter</dd></dl></div><div class='xr-var-data'><pre>array([ 5433892.876412,  5431888.859124,  5429884.841836, ..., -5429884.841836,\n",
       "       -5431888.859124, -5433892.876412])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>x</span></div><div class='xr-var-dims'>(x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-5.434e+06 -5.432e+06 ... 5.434e+06</div><input id='attrs-d27d1b2d-daa5-43ff-9751-b01c5a3006f4' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-d27d1b2d-daa5-43ff-9751-b01c5a3006f4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3d85d148-a42f-41be-8406-da215baaa29c' class='xr-var-data-in' type='checkbox'><label for='data-3d85d148-a42f-41be-8406-da215baaa29c' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>meter</dd></dl></div><div class='xr-var-data'><pre>array([-5433892.876412, -5431888.859124, -5429884.841836, ...,  5429884.841836,\n",
       "        5431888.859124,  5433892.876412])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>bands</span></div><div class='xr-var-dims'>(bands)</div><div class='xr-var-dtype'>&lt;U1</div><div class='xr-var-preview xr-preview'>&#x27;R&#x27; &#x27;G&#x27; &#x27;B&#x27;</div><input id='attrs-4f5fe7c7-8ea4-4fe3-9b78-60fcc914edca' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-4f5fe7c7-8ea4-4fe3-9b78-60fcc914edca' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-1c8cafa7-81bf-418a-90f1-bb2bbb292f72' class='xr-var-data-in' type='checkbox'><label for='data-1c8cafa7-81bf-418a-90f1-bb2bbb292f72' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;R&#x27;, &#x27;G&#x27;, &#x27;B&#x27;], dtype=&#x27;&lt;U1&#x27;)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-2df10e2d-1092-4d7f-b50c-5cd07e20264b' class='xr-section-summary-in' type='checkbox'  ><label for='section-2df10e2d-1092-4d7f-b50c-5cd07e20264b' class='xr-section-summary' >Indexes: <span>(3)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-index-name'><div>y</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-27c5043b-2696-49db-9f6d-62fee39ef593' class='xr-index-data-in' type='checkbox'/><label for='index-27c5043b-2696-49db-9f6d-62fee39ef593' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Float64Index([5433892.8764120005,     5431888.859124,  5429884.841836001,\n",
       "              5427880.8245480005,      5425876.80726,  5423872.789972001,\n",
       "                  5421868.772684,     5419864.755396,  5417860.738108001,\n",
       "                   5415856.72082,\n",
       "              ...\n",
       "                  -5415856.72082,    -5417860.738108, -5419864.755395999,\n",
       "                 -5421868.772684,    -5423872.789972, -5425876.807259999,\n",
       "              -5427880.824547999,    -5429884.841836, -5431888.859123999,\n",
       "              -5433892.876411999],\n",
       "             dtype=&#x27;float64&#x27;, name=&#x27;y&#x27;, length=5424))</pre></div></li><li class='xr-var-item'><div class='xr-index-name'><div>x</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-cf43562e-3c94-4d3d-bae7-8d9a6bf9683b' class='xr-index-data-in' type='checkbox'/><label for='index-cf43562e-3c94-4d3d-bae7-8d9a6bf9683b' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Float64Index([-5433892.8764120005,     -5431888.859124,  -5429884.841836001,\n",
       "              -5427880.8245480005,      -5425876.80726,  -5423872.789972001,\n",
       "                  -5421868.772684,     -5419864.755396,  -5417860.738108001,\n",
       "                   -5415856.72082,\n",
       "              ...\n",
       "                    5415856.72082,      5417860.738108,   5419864.755395999,\n",
       "                   5421868.772684,      5423872.789972,   5425876.807259999,\n",
       "                5427880.824547999,      5429884.841836,   5431888.859123999,\n",
       "                5433892.876411999],\n",
       "             dtype=&#x27;float64&#x27;, name=&#x27;x&#x27;, length=5424))</pre></div></li><li class='xr-var-item'><div class='xr-index-name'><div>bands</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-cbd603a8-35d7-463a-8a93-9e6314cf5d56' class='xr-index-data-in' type='checkbox'/><label for='index-cbd603a8-35d7-463a-8a93-9e6314cf5d56' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([&#x27;R&#x27;, &#x27;G&#x27;, &#x27;B&#x27;], dtype=&#x27;object&#x27;, name=&#x27;bands&#x27;))</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-74a783f2-650e-4dee-9edb-803bba44ab41' class='xr-section-summary-in' type='checkbox'  ><label for='section-74a783f2-650e-4dee-9edb-803bba44ab41' class='xr-section-summary' >Attributes: <span>(29)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>long_name :</span></dt><dd>Brightness Temperature</dd><dt><span>end_time :</span></dt><dd>2021-04-09 23:59:48.350000</dd><dt><span>area :</span></dt><dd>Area ID: GOES-East\n",
       "Description: 2km at nadir\n",
       "Projection ID: abi_fixed_grid\n",
       "Projection: {&#x27;ellps&#x27;: &#x27;GRS80&#x27;, &#x27;h&#x27;: &#x27;35786023&#x27;, &#x27;lon_0&#x27;: &#x27;-75&#x27;, &#x27;no_defs&#x27;: &#x27;None&#x27;, &#x27;proj&#x27;: &#x27;geos&#x27;, &#x27;sweep&#x27;: &#x27;x&#x27;, &#x27;type&#x27;: &#x27;crs&#x27;, &#x27;units&#x27;: &#x27;m&#x27;, &#x27;x_0&#x27;: &#x27;0&#x27;, &#x27;y_0&#x27;: &#x27;0&#x27;}\n",
       "Number of columns: 5424\n",
       "Number of rows: 5424\n",
       "Area extent: (-5434894.8851, -5434894.8851, 5434894.8851, 5434894.8851)</dd><dt><span>orbital_parameters :</span></dt><dd>{&#x27;projection_longitude&#x27;: -75.0, &#x27;projection_latitude&#x27;: 0.0, &#x27;projection_altitude&#x27;: 35786023.0, &#x27;satellite_nominal_latitude&#x27;: 0.0, &#x27;satellite_nominal_longitude&#x27;: -75.19999694824219, &#x27;satellite_nominal_altitude&#x27;: 35786023.4375, &#x27;yaw_flip&#x27;: False}</dd><dt><span>observation_type :</span></dt><dd>Rad</dd><dt><span>reader :</span></dt><dd>abi_l1b</dd><dt><span>start_time :</span></dt><dd>2021-04-09 23:50:17.100000</dd><dt><span>platform_shortname :</span></dt><dd>G16</dd><dt><span>platform_name :</span></dt><dd>GOES-16</dd><dt><span>grid_mapping :</span></dt><dd>goes_imager_projection</dd><dt><span>sensor :</span></dt><dd>abi</dd><dt><span>production_site :</span></dt><dd>WCDAS</dd><dt><span>standard_name :</span></dt><dd>airmass</dd><dt><span>orbital_slot :</span></dt><dd>GOES-East</dd><dt><span>timeline_ID :</span></dt><dd>None</dd><dt><span>scene_abbr :</span></dt><dd>F</dd><dt><span>resolution :</span></dt><dd>2000</dd><dt><span>instrument_ID :</span></dt><dd>FM1</dd><dt><span>scan_mode :</span></dt><dd>M6</dd><dt><span>cell_methods :</span></dt><dd>t: point area: point</dd><dt><span>scene_id :</span></dt><dd>Full Disk</dd><dt><span>ancillary_variables :</span></dt><dd>[]</dd><dt><span>wavelength :</span></dt><dd>None</dd><dt><span>optional_datasets :</span></dt><dd>[]</dd><dt><span>name :</span></dt><dd>airmass</dd><dt><span>_satpy_id :</span></dt><dd>DataID(name=&#x27;airmass&#x27;, resolution=2000)</dd><dt><span>prerequisites :</span></dt><dd>[DataQuery(name=&#x27;_airmass_dep_0&#x27;), DataQuery(name=&#x27;_airmass_dep_1&#x27;), DataQuery(name=&#x27;C08&#x27;)]</dd><dt><span>optional_prerequisites :</span></dt><dd>[]</dd><dt><span>mode :</span></dt><dd>RGB</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray 'where-303720c4e2e42379319fb14493c19145' (bands: 3, y: 5424,\n",
       "                                                            x: 5424)>\n",
       "dask.array<where, shape=(3, 5424, 5424), dtype=float64, chunksize=(1, 4096, 4096), chunktype=numpy.ndarray>\n",
       "Coordinates:\n",
       "    crs      object PROJCRS[\"unknown\",BASEGEOGCRS[\"unknown\",DATUM[\"unknown\",E...\n",
       "  * y        (y) float64 5.434e+06 5.432e+06 5.43e+06 ... -5.432e+06 -5.434e+06\n",
       "  * x        (x) float64 -5.434e+06 -5.432e+06 -5.43e+06 ... 5.432e+06 5.434e+06\n",
       "  * bands    (bands) <U1 'R' 'G' 'B'\n",
       "Attributes: (12/29)\n",
       "    long_name:               Brightness Temperature\n",
       "    end_time:                2021-04-09 23:59:48.350000\n",
       "    area:                    Area ID: GOES-East\\nDescription: 2km at nadir\\nP...\n",
       "    orbital_parameters:      {'projection_longitude': -75.0, 'projection_lati...\n",
       "    observation_type:        Rad\n",
       "    reader:                  abi_l1b\n",
       "    ...                      ...\n",
       "    optional_datasets:       []\n",
       "    name:                    airmass\n",
       "    _satpy_id:               DataID(name='airmass', resolution=2000)\n",
       "    prerequisites:           [DataQuery(name='_airmass_dep_0'), DataQuery(nam...\n",
       "    optional_prerequisites:  []\n",
       "    mode:                    RGB"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Access the dataset named 'airmass' from the 'scn' object using the previously defined variable 'rgb_im'.\n",
    "# The variable 'rgb_im' holds the string 'airmass', which acts as a key to retrieve the corresponding dataset from the Scene object.\n",
    "# This operation is critical in satellite data processing as it allows for direct manipulation and analysis of specific datasets.\n",
    "result = scn[rgb_im]\n",
    "\n",
    "# The retrieved dataset can then be used for further analysis, visualization, or processing.\n",
    "# It's important for students to understand how to efficiently access and work with specific datasets within a complex data structure.\n",
    "result"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=()),\n",
       " DataID(name='C02', wavelength=WavelengthRange(min=0.59, central=0.64, max=0.69, unit='µm'), resolution=500, calibration=<calibration.reflectance>, modifiers=()),\n",
       " DataID(name='C03', wavelength=WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=()),\n",
       " DataID(name='C04', wavelength=WavelengthRange(min=1.3705, central=1.378, max=1.3855, unit='µm'), resolution=2000, calibration=<calibration.reflectance>, modifiers=()),\n",
       " DataID(name='C05', wavelength=WavelengthRange(min=1.58, central=1.61, max=1.64, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=()),\n",
       " DataID(name='C06', wavelength=WavelengthRange(min=2.225, central=2.25, max=2.275, unit='µm'), resolution=2000, calibration=<calibration.reflectance>, modifiers=()),\n",
       " DataID(name='C07', wavelength=WavelengthRange(min=3.8, central=3.9, max=4.0, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
       " DataID(name='C08', wavelength=WavelengthRange(min=5.77, central=6.185, max=6.6, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
       " DataID(name='C09', wavelength=WavelengthRange(min=6.75, central=6.95, max=7.15, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
       " DataID(name='C10', wavelength=WavelengthRange(min=7.24, central=7.34, max=7.44, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
       " DataID(name='C11', wavelength=WavelengthRange(min=8.3, central=8.5, max=8.7, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
       " DataID(name='C12', wavelength=WavelengthRange(min=9.42, central=9.61, max=9.8, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
       " DataID(name='C13', wavelength=WavelengthRange(min=10.1, central=10.35, max=10.6, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
       " DataID(name='C14', wavelength=WavelengthRange(min=10.8, central=11.2, max=11.6, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
       " DataID(name='C15', wavelength=WavelengthRange(min=11.8, central=12.3, max=12.8, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
       " DataID(name='C16', wavelength=WavelengthRange(min=13.0, central=13.3, max=13.6, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
       " DataID(name='airmass', resolution=2000)]"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Retrieve the keys (dataset names) available in the 'scn' object\n",
    "# This function lists all dataset identifiers stored in the Scene object, each represented as a DataID object.\n",
    "keys = scn.keys()\n",
    "\n",
    "# The output is a list of DataID objects, each encapsulating the metadata for a different satellite data channel or composite.\n",
    "# These DataIDs include crucial information such as:\n",
    "# - name: The identifier of the dataset, like 'C01', 'C02', ..., 'C16', 'airmass'.\n",
    "# - wavelength: A WavelengthRange object indicating the spectral range each channel covers, important for identifying the type of observations (e.g., visible, infrared).\n",
    "# - resolution: The spatial resolution of the data in meters, which affects the detail level visible in the imagery.\n",
    "# - calibration: The type of calibration applied to the data, which could be reflectance or brightness temperature, impacting how the data should be interpreted.\n",
    "# - modifiers: Any additional processing applied to the data channel.\n",
    "\n",
    "# Print the keys to show the available datasets and their properties, aiding in understanding the capabilities and focus areas of the satellite instrument.\n",
    "keys"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Area ID: GOES-East\n",
       "Description: 2km at nadir\n",
       "Projection ID: abi_fixed_grid\n",
       "Projection: {'ellps': 'GRS80', 'h': '35786023', 'lon_0': '-75', 'no_defs': 'None', 'proj': 'geos', 'sweep': 'x', 'type': 'crs', 'units': 'm', 'x_0': '0', 'y_0': '0'}\n",
       "Number of columns: 5424\n",
       "Number of rows: 5424\n",
       "Area extent: (-5434894.8851, -5434894.8851, 5434894.8851, 5434894.8851)"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Access the area information associated with the 'C13' dataset in the 'scn' object.\n",
    "# The area property of a dataset within a Scene object provides geographical and geometric details about the satellite data.\n",
    "# This includes information such as the projection, extent, resolution, and size of the dataset, which are crucial for spatial analysis.\n",
    "\n",
    "# For the 'C13' dataset, which typically represents an infrared channel on geostationary satellites, \n",
    "# understanding its area is essential for accurately interpreting spatial phenomena observed in the imagery, such as cloud formations or surface temperatures.\n",
    "area_info = scn[\"C13\"].area\n",
    "\n",
    "# Print the area information to provide insights into the spatial characteristics of the 'C13' dataset.\n",
    "# This information supports tasks such as mapping, data integration with other geospatial datasets, and more precise environmental monitoring.\n",
    "area_info"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Area ID: GOES-East\n",
       "Description: 1km at nadir\n",
       "Projection ID: abi_fixed_grid\n",
       "Projection: {'ellps': 'GRS80', 'h': '35786023', 'lon_0': '-75', 'no_defs': 'None', 'proj': 'geos', 'sweep': 'x', 'type': 'crs', 'units': 'm', 'x_0': '0', 'y_0': '0'}\n",
       "Number of columns: 10848\n",
       "Number of rows: 10848\n",
       "Area extent: (-5434894.8851, -5434894.8851, 5434894.8851, 5434894.8851)"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Access the area definition information for the 'C01' dataset in the 'scn' object.\n",
    "# The 'area' property of a dataset provides detailed geographic and geometric information about where the satellite data was collected.\n",
    "# This includes details such as projection type, coordinate reference system, image extent, and pixel resolution.\n",
    "\n",
    "# The 'C01' channel, which often captures data in a visible light wavelength (approximately 0.47 µm),\n",
    "# provides high-resolution imagery useful for detailed visual inspections of cloud cover, surface features, and atmospheric conditions.\n",
    "area_info = scn[\"C01\"].area\n",
    "\n",
    "# Print the area information to give insights into the geographic scope and detail level of the 'C01' dataset.\n",
    "# This is crucial for applications such as mapping, tracking environmental changes, and integrating satellite data with other geographic information systems (GIS).\n",
    "area_info"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Area ID: GOES-East\n",
       "Description: 0.5km at nadir\n",
       "Projection ID: abi_fixed_grid\n",
       "Projection: {'ellps': 'GRS80', 'h': '35786023', 'lon_0': '-75', 'no_defs': 'None', 'proj': 'geos', 'sweep': 'x', 'type': 'crs', 'units': 'm', 'x_0': '0', 'y_0': '0'}\n",
       "Number of columns: 21696\n",
       "Number of rows: 21696\n",
       "Area extent: (-5434894.8851, -5434894.8851, 5434894.8851, 5434894.8851)"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Access the area definition information for the 'C02' dataset in the 'scn' object.\n",
    "# The 'area' property of a dataset provides crucial geographic and geometric information about the data's capture region.\n",
    "# This includes the projection type, coordinate system, extent of the image, and pixel resolution.\n",
    "\n",
    "# The 'C02' channel is typically in the visible spectrum (centered around 0.64 µm), offering detailed imagery suitable for analyzing surface features,\n",
    "# cloud formations, and atmospheric phenomena. This channel is instrumental in meteorological analysis and environmental monitoring.\n",
    "area_info = scn[\"C02\"].area\n",
    "\n",
    "# Print the area information to provide insights into the geographic scope and resolution of the 'C02' dataset.\n",
    "# Understanding this information is vital for accurate mapping, environmental monitoring, and integration with other geospatial data sources.\n",
    "area_info"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING:satpy.scene:The following datasets were not created and may require resampling to be generated: DataID(name='natural_color')\n"
     ]
    }
   ],
   "source": [
    "# Load the \"natural_color\" dataset using the 'scn.load' method.\n",
    "# The \"natural_color\" composite is a popular visual representation that combines multiple spectral bands\n",
    "# to produce an image that approximates what the human eye would see from space. \n",
    "# This is particularly useful in earth observation for visualizing land cover, water bodies, and other natural features.\n",
    "\n",
    "# The composite typically leverages red, green, and blue spectral bands to enhance the clarity and detail of the image,\n",
    "# which makes it ideal for presentations, educational purposes, and initial visual assessments in environmental studies.\n",
    "scn.load([\"natural_color\"])\n",
    "\n",
    "# This step is crucial for subsequent visualization and analysis tasks, as it prepares the data for easy interpretation and application."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "INFO:satpy.resample:Using default KDTree resampler\n",
      "INFO:satpy.resample:Using default KDTree resampler\n",
      "INFO:satpy.resample:Using default KDTree resampler\n"
     ]
    }
   ],
   "source": [
    "# Get the area definition of the \"C13\" dataset from the 'scn' object.\n",
    "# The 'area' property contains critical geographic and geometric information, including projection, extent, and resolution,\n",
    "# which is essential for accurate geographic referencing in satellite imagery analysis.\n",
    "rs = scn[\"C13\"].area\n",
    "\n",
    "# Resample the scene to the specified area definition.\n",
    "# Resampling is a critical process in satellite data handling, allowing datasets to be standardized to a common spatial grid.\n",
    "# This standardization is necessary for accurate comparison and integration of different data types and sources.\n",
    "\n",
    "# Here, 'scn.resample(rs)' adjusts all data in the scene to match the area definition of the 'C13' dataset.\n",
    "# This is particularly useful when preparing data for detailed analysis or visualization,\n",
    "# ensuring consistency across different datasets within the same scene.\n",
    "lscn = scn.resample(rs)\n",
    "\n",
    "# The resulting 'lscn' is a new Scene object containing all the original data,\n",
    "# but now aligned to the same geographic grid as the 'C13' dataset.\n",
    "# This uniformity is crucial for subsequent processing steps, such as creating composites or conducting multi-temporal analysis."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Load the \"natural_color\" dataset for the resampled scene.\n",
    "# The 'natural_color' composite is designed to mimic the colors visible to the human eye by combining specific spectral bands.\n",
    "# Loading this dataset into the resampled scene ('lscn') ensures that the data is properly aligned and standardized across the same geographic area.\n",
    "\n",
    "lscn.load([\"natural_color\"])\n",
    "\n",
    "# Display the \"natural_color\" dataset.\n",
    "# This step involves visualizing the loaded dataset using the 'show' method, which renders the satellite image according to the data's true or natural color representation.\n",
    "# Visualizing data in this way is particularly useful for presentations and educational purposes, providing a realistic view of earth's features such as vegetation, water bodies, and urban areas.\n",
    "\n",
    "# lscn.show(\"natural_color\")\n",
    "\n",
    "# This visualization step is crucial for analyzing environmental and atmospheric conditions,\n",
    "# as it allows observers to easily identify and assess visible features without the need for specialized image interpretation skills."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Crop the resampled scene (lscn) to a specific geographic region defined by latitude and longitude bounds.\n",
    "# The bounds are given as a tuple in the format (lon_min, lat_min, lon_max, lat_max).\n",
    "# This operation is useful for focusing the analysis on a particular area of interest, reducing data volume and enhancing processing efficiency.\n",
    "scn_c1 = lscn.crop(ll_bbox=(-65.7, 10.7, -55.9, 20.1))\n",
    "\n",
    "# Display the \"C13\" dataset from the cropped scene (scn_c1).\n",
    "# The \"C13\" channel typically represents an infrared wavelength used for observing cloud structure and surface temperature,\n",
    "# crucial for meteorological studies and weather forecasting.\n",
    "# Displaying this dataset allows for detailed observation of atmospheric conditions within the specified geographic area.\n",
    "\n",
    "### Uncomment to show it\n",
    "#scn_c1.show(\"C13\")\n",
    "\n",
    "# Visualizing specific channels like \"C13\" in defined geographic regions helps in targeted analysis,\n",
    "# such as monitoring storm development or evaluating climate patterns in detail.\n",
    "# This capability is particularly valuable in educational settings for demonstrating real-world applications of satellite data analysis."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Load the \"ash\" dataset from the cropped scene (scn_c1).\n",
    "# The \"ash\" composite is specifically designed to detect and visualize volcanic ash in the atmosphere using satellite data.\n",
    "# This dataset is crucial for monitoring volcanic activity and assessing the distribution and movement of ash clouds,\n",
    "# which can have significant impacts on air quality and aviation safety.\n",
    "scn_c1.load(['ash'])\n",
    "\n",
    "# Display the \"ash\" dataset from the cropped scene (scn_c1).\n",
    "# Displaying this dataset allows for visual assessment of ash presence within the specified geographic area.\n",
    "# The visualization is particularly useful in educational and operational settings for demonstrating how satellite data can be applied to real-world environmental challenges.\n",
    "\n",
    "### Uncomment to show it\n",
    "# scn_c1.show('ash')\n",
    "\n",
    "# This step not only aids in the educational demonstration of satellite capabilities but also provides practical insights into the management of natural disasters.\n",
    "# Such visualizations are essential tools in emergency response planning and environmental monitoring."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Load the \"so2\" dataset from the cropped scene (scn_c1).\n",
    "# The \"so2\" dataset is designed to detect sulfur dioxide (SO2) concentrations in the atmosphere using satellite data.\n",
    "# SO2 is a significant volcanic gas and industrial pollutant, making this dataset crucial for monitoring air quality and volcanic activity.\n",
    "scn_c1.load(['so2'])\n",
    "\n",
    "# Display the \"so2\" dataset from the cropped scene (scn_c1).\n",
    "# Displaying this dataset allows for a visual assessment of SO2 distribution within the specified geographic area.\n",
    "# The visualization helps in understanding the spatial extent and concentration of sulfur dioxide, which is essential for both environmental and public health assessments.\n",
    "\n",
    "### Uncomment to show it\n",
    "# scn_c1.show('so2')\n",
    "\n",
    "# This step not only aids in the educational demonstration of how satellite imagery can be utilized to monitor environmental pollutants but also provides practical insights into the management of air quality.\n",
    "# Such visualizations are valuable tools for researchers, policymakers, and educators in understanding and addressing atmospheric pollution."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.9"
  },
  "nbdime-conflicts": {
   "local_diff": [
    {
     "diff": [
      {
       "diff": [
        {
         "key": 0,
         "op": "addrange",
         "valuelist": [
          "Python 3"
         ]
        },
        {
         "key": 0,
         "length": 1,
         "op": "removerange"
        }
       ],
       "key": "display_name",
       "op": "patch"
      }
     ],
     "key": "kernelspec",
     "op": "patch"
    }
   ],
   "remote_diff": [
    {
     "diff": [
      {
       "diff": [
        {
         "key": 0,
         "op": "addrange",
         "valuelist": [
          "Python3"
         ]
        },
        {
         "key": 0,
         "length": 1,
         "op": "removerange"
        }
       ],
       "key": "display_name",
       "op": "patch"
      }
     ],
     "key": "kernelspec",
     "op": "patch"
    }
   ]
  },
  "toc-autonumbering": false
 },
 "nbformat": 4,
 "nbformat_minor": 4
}