Skip to content

Instantly share code, notes, and snippets.

@raybellwaves

raybellwaves/satpy_test_all_composites.ipynb

Created September 24, 2020 02:38
Show Gist options
  • Save raybellwaves/c9ef25bbb9819b3b88c1201d11340b00 to your computer and use it in GitHub Desktop.
Save raybellwaves/c9ef25bbb9819b3b88c1201d11340b00 to your computer and use it in GitHub Desktop.
Download ZIP
Display the source blob
Display the rendered blob
Raw
satpy_test_all_composites.ipynb
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"# Create a new env and install satpy and s3fs\n",
"# conda create -n test_env python=3.8\n",
"# conda activate test_env\n",
"# conda install -c conda-forge satpy s3fs jupyter\n",
"#\n",
"# Install\n",
"# satpy==0.23.0\n",
"# pyresample==1.16.0"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"import s3fs\n",
"fs = s3fs.S3FileSystem(anon=True)"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"['noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C01_G16_s20202680000157_e20202680009465_c20202680009527.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C02_G16_s20202680000157_e20202680009465_c20202680009507.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C03_G16_s20202680000157_e20202680009465_c20202680009524.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C04_G16_s20202680000157_e20202680009465_c20202680009495.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C05_G16_s20202680000157_e20202680009465_c20202680009519.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C06_G16_s20202680000157_e20202680009471_c20202680009515.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C07_G16_s20202680000157_e20202680009477_c20202680009531.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C08_G16_s20202680000157_e20202680009465_c20202680009535.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C09_G16_s20202680000157_e20202680009471_c20202680009544.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C10_G16_s20202680000157_e20202680009477_c20202680009532.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C11_G16_s20202680000157_e20202680009465_c20202680009527.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C12_G16_s20202680000157_e20202680009471_c20202680009554.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C13_G16_s20202680000157_e20202680009477_c20202680009536.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C14_G16_s20202680000157_e20202680009465_c20202680009547.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C15_G16_s20202680000157_e20202680009471_c20202680009558.nc',\n",
" 'noaa-goes16/ABI-L1b-RadF/2020/268/00/OR_ABI-L1b-RadF-M6C16_G16_s20202680000157_e20202680009477_c20202680009541.nc']"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"fs.glob('noaa-goes16/ABI-L1b-RadF/2020/268/00/*s20202680000157*')"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[None,\n",
" None,\n",
" None,\n",
" None,\n",
" None,\n",
" None,\n",
" None,\n",
" None,\n",
" None,\n",
" None,\n",
" None,\n",
" None,\n",
" None,\n",
" None,\n",
" None,\n",
" None]"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"fs.get('noaa-goes16/ABI-L1b-RadF/2020/268/00/*s20202680000157*', 'data')"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"import glob"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"filenames = glob.glob('data/*')"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"16"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"len(filenames)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/resample.py:150: UserWarning: Gradient search resampler not available, upgrade Pyresample.\n",
" warnings.warn('Gradient search resampler not available, upgrade Pyresample.')\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"Readers\n",
"=======\n",
"abi_l1b: ok\n",
"abi_l1b_scmi: ok\n",
"abi_l2_nc: ok\n",
"acspo: ok\n",
"agri_l1: ok\n",
"ahi_hrit: ok\n",
"ahi_hsd: ok\n",
"ahi_l1b_gridded_bin: ok\n",
"ami_l1b: ok\n",
"amsr2_l1b: ok\n",
"amsr2_l2: ok\n",
"avhrr_l1b_aapp: ok\n",
"avhrr_l1b_eps: ok\n",
"avhrr_l1b_gaclac: cannot find module 'satpy.readers.avhrr_l1b_gaclac' (No module named 'pygac')\n",
"avhrr_l1b_hrpt: cannot find module 'satpy.readers.hrpt' (No module named 'pygac')\n",
"caliop_l2_cloud: cannot find module 'satpy.readers.caliop_l2_cloud' (No module named 'pyhdf')\n",
"clavrx: cannot find module 'satpy.readers.clavrx' (No module named 'pyhdf')\n",
"cmsaf-claas2_l2_nc: ok\n",
"electrol_hrit: ok\n",
"fci_l1c_fdhsi: ok\n",
"fci_l2_nc: ok\n",
"generic_image: ok\n",
"geocat: ok\n",
"ghrsst_l3c_sst: cannot find module 'satpy.readers.ghrsst_l3c_sst' (cannot import name 'Dataset' from 'satpy.dataset' (/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/dataset/__init__.py))\n",
"glm_l2: ok\n",
"goes-imager_hrit: ok\n",
"goes-imager_nc: ok\n",
"gpm_imerg: ok\n",
"grib: cannot find module 'satpy.readers.grib' (No module named 'pygrib')\n",
"hsaf_grib: cannot find module 'satpy.readers.hsaf_grib' (No module named 'pygrib')\n",
"hy2_scat_l2b_h5: ok\n",
"iasi_l2: ok\n",
"iasi_l2_so2_bufr: cannot find module 'satpy.readers.iasi_l2_so2_bufr' (('Missing eccodes-python and/or eccodes C-library installation. Use conda to install eccodes.\\n Error: ', ModuleNotFoundError(\"No module named 'eccodes'\")))\n",
"jami_hrit: ok\n",
"li_l2: cannot find module 'satpy.readers.li_l2' (No module named 'h5netcdf')\n",
"maia: ok\n",
"mersi2_l1b: ok\n",
"mimicTPW2_comp: ok\n",
"modis_l1b: cannot find module 'satpy.readers.modis_l1b' (No module named 'pyhdf')\n",
"modis_l2: cannot find module 'satpy.readers.modis_l2' (No module named 'pyhdf')\n",
"msi_safe: cannot find module 'satpy.readers.msi_safe' (No module named 'glymur')\n",
"mtsat2-imager_hrit: ok\n",
"nucaps: ok\n",
"nwcsaf-geo: ok\n",
"nwcsaf-msg2013-hdf5: ok\n",
"nwcsaf-pps_nc: ok\n",
"olci_l1b: ok\n",
"olci_l2: ok\n",
"omps_edr: ok\n",
"safe_sar_l2_ocn: ok\n",
"sar-c_safe: ok\n",
"scatsat1_l2b: cannot find module 'satpy.readers.scatsat1_l2b' (cannot import name 'Dataset' from 'satpy.dataset' (/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/dataset/__init__.py))\n",
"seviri_l1b_hrit: ok\n",
"seviri_l1b_icare: cannot find module 'satpy.readers.seviri_l1b_icare' (No module named 'pyhdf')\n",
"seviri_l1b_native: ok\n",
"seviri_l1b_nc: ok\n",
"seviri_l2_bufr: cannot find module 'satpy.readers.seviri_l2_bufr' (Missing eccodes-python and/or eccodes C-library installation. Use conda to install eccodes)\n",
"seviri_l2_grib: cannot find module 'satpy.readers.seviri_l2_grib' (Missing eccodes-python and/or eccodes C-library installation. Use conda to install eccodes)\n",
"slstr_l1b: ok\n",
"slstr_l2: ok\n",
"smos_l2_wind: ok\n",
"tropomi_l2: ok\n",
"vaisala_gld360: ok\n",
"vii_l1b_nc: ok\n",
"vii_l2_nc: ok\n",
"viirs_compact: ok\n",
"viirs_edr_active_fires: ok\n",
"viirs_edr_flood: cannot find module 'satpy.readers.viirs_edr_flood' (No module named 'pyhdf')\n",
"viirs_l1b: ok\n",
"viirs_sdr: ok\n",
"virr_l1b: ok\n",
"\n",
"Writers\n",
"=======\n",
"cf: ok\n",
"geotiff: ok\n",
"mitiff: ok\n",
"ninjotiff: cannot find module 'satpy.writers.ninjotiff' (No module named 'pyninjotiff')\n",
"scmi: ok\n",
"simple_image: ok\n",
"\n",
"Extras\n",
"======\n",
"cartopy: No module named 'cartopy'\n",
"geoviews: No module named 'geoviews'\n",
"\n"
]
}
],
"source": [
"from satpy.config import check_satpy; check_satpy()"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [],
"source": [
"from satpy import Scene\n",
"from satpy.utils import debug_on; debug_on()"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"[DEBUG: 2020-09-23 22:35:43 : satpy.scene] Setting 'PPP_CONFIG_DIR' to '/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/etc'\n",
"[DEBUG: 2020-09-23 22:35:43 : satpy.readers.yaml_reader] Reading ('/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/etc/readers/abi_l1b.yaml',)\n",
"[DEBUG: 2020-09-23 22:35:43 : satpy.readers.yaml_reader] Assigning to abi_l1b: ['data/OR_ABI-L1b-RadF-M6C01_G16_s20202680000157_e20202680009465_c20202680009527.nc', 'data/OR_ABI-L1b-RadF-M6C02_G16_s20202680000157_e20202680009465_c20202680009507.nc', 'data/OR_ABI-L1b-RadF-M6C03_G16_s20202680000157_e20202680009465_c20202680009524.nc', 'data/OR_ABI-L1b-RadF-M6C04_G16_s20202680000157_e20202680009465_c20202680009495.nc', 'data/OR_ABI-L1b-RadF-M6C05_G16_s20202680000157_e20202680009465_c20202680009519.nc', 'data/OR_ABI-L1b-RadF-M6C06_G16_s20202680000157_e20202680009471_c20202680009515.nc', 'data/OR_ABI-L1b-RadF-M6C07_G16_s20202680000157_e20202680009477_c20202680009531.nc', 'data/OR_ABI-L1b-RadF-M6C08_G16_s20202680000157_e20202680009465_c20202680009535.nc', 'data/OR_ABI-L1b-RadF-M6C09_G16_s20202680000157_e20202680009471_c20202680009544.nc', 'data/OR_ABI-L1b-RadF-M6C10_G16_s20202680000157_e20202680009477_c20202680009532.nc', 'data/OR_ABI-L1b-RadF-M6C11_G16_s20202680000157_e20202680009465_c20202680009527.nc', 'data/OR_ABI-L1b-RadF-M6C12_G16_s20202680000157_e20202680009471_c20202680009554.nc', 'data/OR_ABI-L1b-RadF-M6C13_G16_s20202680000157_e20202680009477_c20202680009536.nc', 'data/OR_ABI-L1b-RadF-M6C14_G16_s20202680000157_e20202680009465_c20202680009547.nc', 'data/OR_ABI-L1b-RadF-M6C15_G16_s20202680000157_e20202680009471_c20202680009558.nc', 'data/OR_ABI-L1b-RadF-M6C16_G16_s20202680000157_e20202680009477_c20202680009541.nc']\n",
"[DEBUG: 2020-09-23 22:35:43 : satpy.composites] Looking for composites config file abi.yaml\n",
"[DEBUG: 2020-09-23 22:35:43 : satpy.composites] Looking for composites config file visir.yaml\n"
]
}
],
"source": [
"scn = Scene(reader='abi_l1b', filenames=filenames)"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='brightness_temperature'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='brightness_temperature')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='satellite_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='satellite_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='satellite_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='satellite_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='satellite_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='satellite_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='satellite_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='satellite_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='satellite_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='satellite_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='satellite_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='satellite_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:35:53 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n"
]
},
{
"data": {
"text/plain": [
"['C01',\n",
" 'C02',\n",
" 'C03',\n",
" 'C04',\n",
" 'C05',\n",
" 'C06',\n",
" 'C07',\n",
" 'C08',\n",
" 'C09',\n",
" 'C10',\n",
" 'C11',\n",
" 'C12',\n",
" 'C13',\n",
" 'C14',\n",
" 'C15',\n",
" 'C16',\n",
" 'airmass',\n",
" 'ash',\n",
" 'cimss_cloud_type',\n",
" 'cimss_green',\n",
" 'cimss_green_sunz',\n",
" 'cimss_green_sunz_rayleigh',\n",
" 'cimss_true_color',\n",
" 'cimss_true_color_sunz',\n",
" 'cimss_true_color_sunz_rayleigh',\n",
" 'cira_day_convection',\n",
" 'cira_fire_temperature',\n",
" 'cloud_phase',\n",
" 'cloud_phase_distinction',\n",
" 'cloud_phase_distinction_raw',\n",
" 'cloud_phase_raw',\n",
" 'cloudtop',\n",
" 'color_infrared',\n",
" 'colorized_ir_clouds',\n",
" 'convection',\n",
" 'day_microphysics',\n",
" 'day_microphysics_abi',\n",
" 'day_microphysics_eum',\n",
" 'dust',\n",
" 'fire_temperature_awips',\n",
" 'fog',\n",
" 'green',\n",
" 'green_crefl',\n",
" 'green_nocorr',\n",
" 'green_raw',\n",
" 'green_snow',\n",
" 'ir108_3d',\n",
" 'ir_cloud_day',\n",
" 'land_cloud',\n",
" 'land_cloud_fire',\n",
" 'natural_color',\n",
" 'natural_color_nocorr',\n",
" 'natural_color_raw',\n",
" 'night_fog',\n",
" 'night_ir_alpha',\n",
" 'night_ir_with_background',\n",
" 'night_ir_with_background_hires',\n",
" 'night_microphysics',\n",
" 'night_microphysics_abi',\n",
" 'overview',\n",
" 'overview_raw',\n",
" 'snow',\n",
" 'snow_fog',\n",
" 'so2',\n",
" 'tropical_airmass',\n",
" 'true_color',\n",
" 'true_color_crefl',\n",
" 'true_color_nocorr',\n",
" 'true_color_raw',\n",
" 'true_color_with_night_ir',\n",
" 'true_color_with_night_ir_hires',\n",
" 'water_vapors1',\n",
" 'water_vapors2']"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"datasets = scn.available_dataset_names(composites=True)\n",
"datasets"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='satellite_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='satellite_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='solar_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='satellite_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='satellite_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='solar_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='satellite_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='satellite_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='solar_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:01 : satpy.node] Skipping optional DataQuery(name='satellite_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.node] Skipping optional DataQuery(name='satellite_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.node] Skipping optional DataQuery(name='solar_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.node] Skipping optional DataQuery(name='satellite_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.node] Skipping optional DataQuery(name='satellite_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.node] Skipping optional DataQuery(name='solar_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='brightness_temperature'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='brightness_temperature')\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.node] Skipping optional DataQuery(name='satellite_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.node] Skipping optional DataQuery(name='satellite_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='satellite_zenith_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.node] Skipping optional DataQuery(name='solar_azimuth_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_azimuth_angle', calibration='reflectance')\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.node] Skipping optional DataQuery(name='solar_zenith_angle', calibration='reflectance'): Unknown dataset DataQuery(name='solar_zenith_angle', calibration='reflectance')\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C15.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to brightness temperatures\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C16.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to brightness temperatures\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C14.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to brightness temperatures\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C04.\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to reflectances\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C12.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to brightness temperatures\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C10.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to brightness temperatures\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C02.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to reflectances\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C03.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to reflectances\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C06.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to reflectances\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C08.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to brightness temperatures\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C11.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to brightness temperatures\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C01.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to reflectances\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C09.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to brightness temperatures\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C05.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to reflectances\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C07.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to brightness temperatures\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Reading in get_dataset C13.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.readers.abi_l1b] Calibrating to brightness temperatures\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.composites] Applying sun zen correction\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.composites] Computing sun zenith angles.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.composites] Apply the standard sun-zenith correction [1/cos(sunz)]\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.composites] Sun-zenith correction applied. Computation time: 0.0 (sec)\n",
"[INFO: 2020-09-23 22:36:02 : satpy.composites] Getting reflective part of C07\n",
"[DEBUG: 2020-09-23 22:36:02 : pyspectral.rsr_reader] Filename: /home/ray/.local/share/pyspectral/rsr_abi_GOES-16.h5\n",
"[DEBUG: 2020-09-23 22:36:02 : pyspectral.rsr_reader] Filename: /home/ray/.local/share/pyspectral/rsr_abi_GOES-16.h5\n",
"[DEBUG: 2020-09-23 22:36:02 : pyspectral.rsr_reader] No detectors found - assume only one...\n",
"[DEBUG: 2020-09-23 22:36:02 : pyspectral.solar] Begin and end wavelength/wavenumber: 3.694672 4.076807 \n",
"[INFO: 2020-09-23 22:36:02 : pyspectral.near_infrared_reflectance] No lut filename available in config file. Will generate filename automatically\n",
"[INFO: 2020-09-23 22:36:02 : pyspectral.near_infrared_reflectance] lut filename: /tmp/tb2rad_lut_goes-16_abi_ch7.npz\n",
"[INFO: 2020-09-23 22:36:02 : pyspectral.near_infrared_reflectance] File was there and has been read!\n",
"[DEBUG: 2020-09-23 22:36:02 : pyspectral.near_infrared_reflectance] Apply sun-zenith angle clipping between 0 and 85.00\n",
"[DEBUG: 2020-09-23 22:36:02 : pyspectral.near_infrared_reflectance] Derive the 3.9 micron radiance CO2 correction coefficent\n",
"[INFO: 2020-09-23 22:36:02 : pyspectral.near_infrared_reflectance] CO2 correction applied...\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.scene] Delaying generation of DataID(name='day_microphysics') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.composites] Applying sun zen correction\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.composites] Computing sun zenith angles.\n",
"[DEBUG: 2020-09-23 22:36:02 : satpy.composites] Apply the standard sun-zenith correction [1/cos(sunz)]\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"[DEBUG: 2020-09-23 22:36:02 : satpy.composites] Sun-zenith correction applied. Computation time: 0.0 (sec)\n",
"[INFO: 2020-09-23 22:36:02 : satpy.composites] Removing Rayleigh scattering with atmosphere 'us-standard' and aerosol type 'rayleigh_only' for 'C02'\n",
"[INFO: 2020-09-23 22:36:03 : pyspectral.rayleigh] Atmosphere chosen: us-standard\n",
"[DEBUG: 2020-09-23 22:36:03 : pyspectral.rayleigh] LUT filename: /home/ray/.local/share/pyspectral/rayleigh_only/rayleigh_lut_us-standard.h5\n",
"[DEBUG: 2020-09-23 22:36:03 : pyspectral.rsr_reader] Filename: /home/ray/.local/share/pyspectral/rsr_abi_GOES-16.h5\n",
"[DEBUG: 2020-09-23 22:36:03 : pyspectral.rsr_reader] Filename: /home/ray/.local/share/pyspectral/rsr_abi_GOES-16.h5\n",
"[DEBUG: 2020-09-23 22:36:03 : pyspectral.rsr_reader] No detectors found - assume only one...\n",
"[DEBUG: 2020-09-23 22:36:03 : pyspectral.rayleigh] Band name: ch2 Effective wavelength: 0.635429\n",
"[DEBUG: 2020-09-23 22:36:03 : pyspectral.rayleigh] Time - Interpolation: 0.028588\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.composites] Applying sun zen correction\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.composites] Computing sun zenith angles.\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.composites] Apply the standard sun-zenith correction [1/cos(sunz)]\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.composites] Sun-zenith correction applied. Computation time: 0.0 (sec)\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected')) because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='green') because of dependency's delayed generation: DataID(name='C03', wavelength=WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delayed prerequisite for 'DataID(name='green')': 'DataID(name='C03', wavelength=WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))'\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='true_color') because of dependency's delayed generation: DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected'))\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delayed prerequisite for 'DataID(name='true_color')': 'DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected'))'\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Setting 'PPP_CONFIG_DIR' to '/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/etc'\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.readers.yaml_reader] Reading ('/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/etc/readers/generic_image.yaml',)\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.readers.yaml_reader] Assigning to generic_image: ['BlackMarble_2016_3km_geo.tif']\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] Entering env context: <rasterio.env.Env object at 0x7f346c074f70>\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] Starting outermost env\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] No GDAL environment exists\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] New GDAL environment <rasterio._env.GDALEnv object at 0x7f346c23d610> created\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._env] GDAL_DATA found in environment: '/home/ray/local/bin/anaconda3/envs/test_env/share/gdal'.\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._env] PROJ_LIB found in environment: '/home/ray/local/bin/anaconda3/envs/test_env/share/proj'.\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._env] Started GDALEnv <rasterio._env.GDALEnv object at 0x7f346c23d610>.\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] Entered env context: <rasterio.env.Env object at 0x7f346c074f70>\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._base] Sharing flag: 0\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._base] Nodata success: 0, Nodata value: -10000000000.000000\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._base] Nodata success: 0, Nodata value: -10000000000.000000\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._base] Nodata success: 0, Nodata value: -10000000000.000000\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._base] Dataset <open DatasetReader name='BlackMarble_2016_3km_geo.tif' mode='r'> is started.\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] Exiting env context: <rasterio.env.Env object at 0x7f346c074f70>\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] Cleared existing <rasterio._env.GDALEnv object at 0x7f346c23d610> options\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._env] Stopped GDALEnv <rasterio._env.GDALEnv object at 0x7f346c23d610>.\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] Exiting outermost env\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] Exited env context: <rasterio.env.Env object at 0x7f346c074f70>\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:53: FutureWarning: '+init=<authority>:<code>' syntax is deprecated. '<authority>:<code>' is the preferred initialization method. When making the change, be mindful of axis order changes: https://pyproj4.github.io/pyproj/stable/gotchas.html#axis-order-changes-in-proj-6\n",
" return _prepare_from_string(\" \".join(pjargs))\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] Entering env context: <rasterio.env.Env object at 0x7f346c09b8b0>\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] Starting outermost env\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] No GDAL environment exists\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] New GDAL environment <rasterio._env.GDALEnv object at 0x7f346c26ba30> created\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._env] GDAL_DATA found in environment: '/home/ray/local/bin/anaconda3/envs/test_env/share/gdal'.\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._env] PROJ_LIB found in environment: '/home/ray/local/bin/anaconda3/envs/test_env/share/proj'.\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._env] Started GDALEnv <rasterio._env.GDALEnv object at 0x7f346c26ba30>.\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] Entered env context: <rasterio.env.Env object at 0x7f346c09b8b0>\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._base] Sharing flag: 0\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._base] Nodata success: 0, Nodata value: -10000000000.000000\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._base] Nodata success: 0, Nodata value: -10000000000.000000\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._base] Nodata success: 0, Nodata value: -10000000000.000000\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._base] Dataset <open DatasetReader name='BlackMarble_2016_3km_geo.tif' mode='r'> is started.\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] Exiting env context: <rasterio.env.Env object at 0x7f346c09b8b0>\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] Cleared existing <rasterio._env.GDALEnv object at 0x7f346c26ba30> options\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio._env] Stopped GDALEnv <rasterio._env.GDALEnv object at 0x7f346c26ba30>.\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] Exiting outermost env\n",
"[DEBUG: 2020-09-23 22:36:03 : rasterio.env] Exited env context: <rasterio.env.Env object at 0x7f346c09b8b0>\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.composites] Looking for composites config file images.yaml\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.composites] No composite config found called images.yaml\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.readers.generic_image] Reading DataID(name='image', modifiers=()).\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='night_ir_with_background_hires') because of incompatible areas\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='true_color_with_night_ir_hires') because of dependency's delayed generation: DataID(name='night_ir_with_background_hires')\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delayed prerequisite for 'DataID(name='true_color_with_night_ir_hires')': 'DataID(name='night_ir_with_background_hires')'\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.composites.crefl_utils] No average elevation information provided in CREFL\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.composites.crefl_utils] Using ABI CREFL algorithm\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/dask/array/core.py:1377: FutureWarning: The `numpy.polyval` function is not implemented by Dask array. You may want to use the da.map_blocks function or something similar to silence this warning. Your code may stop working in a future release.\n",
" warnings.warn(\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='_convection_dep_2') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='C03', wavelength=WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected')) because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='_cira_day_convection_dep_2') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='cira_day_convection') because of dependency's delayed generation: DataID(name='_cira_day_convection_dep_2')\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delayed prerequisite for 'DataID(name='cira_day_convection')': 'DataID(name='_cira_day_convection_dep_2')'\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='cimss_green') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='cimss_true_color') because of dependency's delayed generation: DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delayed prerequisite for 'DataID(name='cimss_true_color')': 'DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=())'\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='color_infrared') because of dependency's delayed generation: DataID(name='green')\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delayed prerequisite for 'DataID(name='color_infrared')': 'DataID(name='green')'\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='cimss_green_sunz') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='cimss_green_sunz_rayleigh') because of dependency's delayed generation: DataID(name='C03', wavelength=WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delayed prerequisite for 'DataID(name='cimss_green_sunz_rayleigh')': 'DataID(name='C03', wavelength=WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))'\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='cimss_true_color_sunz_rayleigh') because of dependency's delayed generation: DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected'))\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delayed prerequisite for 'DataID(name='cimss_true_color_sunz_rayleigh')': 'DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected'))'\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='cimss_true_color_sunz') because of dependency's delayed generation: DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delayed prerequisite for 'DataID(name='cimss_true_color_sunz')': 'DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))'\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='cimss_green') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.composites.crefl_utils] No average elevation information provided in CREFL\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.composites.crefl_utils] Using ABI CREFL algorithm\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/dask/array/core.py:1377: FutureWarning: The `numpy.polyval` function is not implemented by Dask array. You may want to use the da.map_blocks function or something similar to silence this warning. Your code may stop working in a future release.\n",
" warnings.warn(\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='green_nocorr') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='true_color_nocorr') because of dependency's delayed generation: DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delayed prerequisite for 'DataID(name='true_color_nocorr')': 'DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=())'\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='cloud_phase_distinction_raw') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='_cira_day_convection_dep_2') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.composites.crefl_utils] No average elevation information provided in CREFL\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.composites.crefl_utils] Using ABI CREFL algorithm\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/dask/array/core.py:1377: FutureWarning: The `numpy.polyval` function is not implemented by Dask array. You may want to use the da.map_blocks function or something similar to silence this warning. Your code may stop working in a future release.\n",
" warnings.warn(\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='natural_color_raw') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='green') because of dependency's delayed generation: DataID(name='C03', wavelength=WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delayed prerequisite for 'DataID(name='green')': 'DataID(name='C03', wavelength=WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))'\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='green_raw') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delaying generation of DataID(name='true_color_raw') because of dependency's delayed generation: DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))\n",
"[DEBUG: 2020-09-23 22:36:03 : satpy.scene] Delayed prerequisite for 'DataID(name='true_color_raw')': 'DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))'\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='land_cloud_fire') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Setting 'PPP_CONFIG_DIR' to '/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/etc'\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.readers.yaml_reader] Reading ('/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/etc/readers/generic_image.yaml',)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.readers.yaml_reader] Assigning to generic_image: ['BlackMarble_2016_01deg_geo.tif']\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] Entering env context: <rasterio.env.Env object at 0x7f3468598ca0>\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] Starting outermost env\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] No GDAL environment exists\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] New GDAL environment <rasterio._env.GDALEnv object at 0x7f34685f97c0> created\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._env] GDAL_DATA found in environment: '/home/ray/local/bin/anaconda3/envs/test_env/share/gdal'.\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._env] PROJ_LIB found in environment: '/home/ray/local/bin/anaconda3/envs/test_env/share/proj'.\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._env] Started GDALEnv <rasterio._env.GDALEnv object at 0x7f34685f97c0>.\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] Entered env context: <rasterio.env.Env object at 0x7f3468598ca0>\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._base] Sharing flag: 0\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._base] Nodata success: 0, Nodata value: -10000000000.000000\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._base] Nodata success: 0, Nodata value: -10000000000.000000\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._base] Nodata success: 0, Nodata value: -10000000000.000000\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._base] Dataset <open DatasetReader name='BlackMarble_2016_01deg_geo.tif' mode='r'> is started.\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] Exiting env context: <rasterio.env.Env object at 0x7f3468598ca0>\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] Cleared existing <rasterio._env.GDALEnv object at 0x7f34685f97c0> options\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._env] Stopped GDALEnv <rasterio._env.GDALEnv object at 0x7f34685f97c0>.\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] Exiting outermost env\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] Exited env context: <rasterio.env.Env object at 0x7f3468598ca0>\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:53: FutureWarning: '+init=<authority>:<code>' syntax is deprecated. '<authority>:<code>' is the preferred initialization method. When making the change, be mindful of axis order changes: https://pyproj4.github.io/pyproj/stable/gotchas.html#axis-order-changes-in-proj-6\n",
" return _prepare_from_string(\" \".join(pjargs))\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] Entering env context: <rasterio.env.Env object at 0x7f34685dc160>\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] Starting outermost env\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] No GDAL environment exists\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] New GDAL environment <rasterio._env.GDALEnv object at 0x7f3468598490> created\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._env] GDAL_DATA found in environment: '/home/ray/local/bin/anaconda3/envs/test_env/share/gdal'.\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._env] PROJ_LIB found in environment: '/home/ray/local/bin/anaconda3/envs/test_env/share/proj'.\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._env] Started GDALEnv <rasterio._env.GDALEnv object at 0x7f3468598490>.\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] Entered env context: <rasterio.env.Env object at 0x7f34685dc160>\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._base] Sharing flag: 0\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._base] Nodata success: 0, Nodata value: -10000000000.000000\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._base] Nodata success: 0, Nodata value: -10000000000.000000\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._base] Nodata success: 0, Nodata value: -10000000000.000000\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._base] Dataset <open DatasetReader name='BlackMarble_2016_01deg_geo.tif' mode='r'> is started.\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] Exiting env context: <rasterio.env.Env object at 0x7f34685dc160>\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] Cleared existing <rasterio._env.GDALEnv object at 0x7f3468598490> options\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio._env] Stopped GDALEnv <rasterio._env.GDALEnv object at 0x7f3468598490>.\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] Exiting outermost env\n",
"[DEBUG: 2020-09-23 22:36:04 : rasterio.env] Exited env context: <rasterio.env.Env object at 0x7f34685dc160>\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.composites] Looking for composites config file images.yaml\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.composites] No composite config found called images.yaml\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.readers.generic_image] Reading DataID(name='image', modifiers=()).\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='night_ir_with_background') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='green_crefl') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='green_snow') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='day_microphysics_abi') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='overview_raw') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='cira_fire_temperature') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='cloud_phase_raw') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.composites] Applying sun zen correction\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.composites] Computing sun zenith angles.\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.composites] Apply the standard sun-zenith correction [1/cos(sunz)]\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.composites] Sun-zenith correction applied. Computation time: 0.0 (sec)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='snow_fog') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='green_raw') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='snow') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.composites] Applying sun zen correction\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.composites] Computing sun zenith angles.\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.composites] Apply the standard sun-zenith correction [1/cos(sunz)]\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.composites] Sun-zenith correction applied. Computation time: 0.0 (sec)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='convection') because of dependency's delayed generation: DataID(name='_convection_dep_2')\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delayed prerequisite for 'DataID(name='convection')': 'DataID(name='_convection_dep_2')'\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='true_color') because of dependency's delayed generation: DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected'))\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delayed prerequisite for 'DataID(name='true_color')': 'DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected'))'\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='day_microphysics_eum') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='green_nocorr') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='cimss_cloud_type') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='true_color_with_night_ir') because of dependency's delayed generation: DataID(name='night_ir_with_background')\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delayed prerequisite for 'DataID(name='true_color_with_night_ir')': 'DataID(name='night_ir_with_background')'\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='true_color_crefl') because of dependency's delayed generation: DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected_crefl'))\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delayed prerequisite for 'DataID(name='true_color_crefl')': 'DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected_crefl'))'\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='natural_color_nocorr') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='cloud_phase') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='fire_temperature_awips') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected')) because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='natural_color') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='land_cloud') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='cloud_phase_distinction') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='cimss_green_sunz_rayleigh') because of dependency's delayed generation: DataID(name='C03', wavelength=WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delayed prerequisite for 'DataID(name='cimss_green_sunz_rayleigh')': 'DataID(name='C03', wavelength=WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))'\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='overview') because of incompatible areas\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Delaying generation of DataID(name='night_ir_with_background_hires') because of incompatible areas\n",
"[WARNING: 2020-09-23 22:36:04 : satpy.scene] The following datasets were not created and may require resampling to be generated: DataID(name='day_microphysics'), DataID(name='day_microphysics_abi'), DataID(name='overview_raw'), DataID(name='true_color_with_night_ir_hires'), DataID(name='cira_fire_temperature'), DataID(name='cloud_phase_raw'), DataID(name='snow_fog'), DataID(name='cira_day_convection'), DataID(name='green_raw'), DataID(name='snow'), DataID(name='cimss_true_color'), DataID(name='color_infrared'), DataID(name='cimss_green_sunz'), DataID(name='cimss_true_color_sunz_rayleigh'), DataID(name='cimss_true_color_sunz'), DataID(name='convection'), DataID(name='true_color'), DataID(name='cimss_green'), DataID(name='day_microphysics_eum'), DataID(name='green_nocorr'), DataID(name='cimss_cloud_type'), DataID(name='true_color_nocorr'), DataID(name='true_color_with_night_ir'), DataID(name='cloud_phase_distinction_raw'), DataID(name='true_color_crefl'), DataID(name='natural_color_nocorr'), DataID(name='cloud_phase'), DataID(name='fire_temperature_awips'), DataID(name='natural_color_raw'), DataID(name='true_color_raw'), DataID(name='green'), DataID(name='natural_color'), DataID(name='land_cloud'), DataID(name='land_cloud_fire'), DataID(name='cloud_phase_distinction'), DataID(name='night_ir_with_background'), DataID(name='green_crefl'), DataID(name='cimss_green_sunz_rayleigh'), DataID(name='overview'), DataID(name='night_ir_with_background_hires'), DataID(name='green_snow')\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_night_microphysics_dep_0', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_night_fog_dep_0', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_night_fog_dep_1', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_dust_dep_1', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_so2_dep_1', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_airmass_dep_0', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_ash_dep_1', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_dust_dep_0', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_water_vapors2_dep_0', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_tropical_airmass_dep_0', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_tropical_airmass_dep_1', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_night_microphysics_abi_dep_1', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_ash_dep_0', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_night_microphysics_dep_1', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_fog_dep_1', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_airmass_dep_1', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_so2_dep_0', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_night_microphysics_abi_dep_0', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:04 : satpy.scene] Unloading dataset: DataID(name='_fog_dep_0', resolution=2000)\n"
]
}
],
"source": [
"scn.load(datasets)"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Bidirectional Reflectance',\n",
" 'standard_name': 'toa_bidirectional_reflectance',\n",
" 'sensor_band_bit_depth': 10,\n",
" 'units': '%',\n",
" 'resolution': 1000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C01',\n",
" 'wavelength': WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'),\n",
" 'calibration': 'reflectance',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 46, 500000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C01'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Bidirectional Reflectance',\n",
" 'standard_name': 'toa_bidirectional_reflectance',\n",
" 'sensor_band_bit_depth': 12,\n",
" 'units': '%',\n",
" 'resolution': 500,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C02',\n",
" 'wavelength': WavelengthRange(min=0.59, central=0.64, max=0.69, unit='µm'),\n",
" 'calibration': 'reflectance',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 46, 500000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C02', wavelength=WavelengthRange(min=0.59, central=0.64, max=0.69, unit='µm'), resolution=500, calibration=<calibration.reflectance>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C02'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Bidirectional Reflectance',\n",
" 'standard_name': 'toa_bidirectional_reflectance',\n",
" 'sensor_band_bit_depth': 10,\n",
" 'units': '%',\n",
" 'resolution': 1000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C03',\n",
" 'wavelength': WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'),\n",
" 'calibration': 'reflectance',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 46, 500000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C03', wavelength=WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 14,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C03'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Bidirectional Reflectance',\n",
" 'standard_name': 'toa_bidirectional_reflectance',\n",
" 'sensor_band_bit_depth': 11,\n",
" 'units': '%',\n",
" 'resolution': 2000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C04',\n",
" 'wavelength': WavelengthRange(min=1.3705, central=1.378, max=1.3855, unit='µm'),\n",
" 'calibration': 'reflectance',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 46, 500000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C04', wavelength=WavelengthRange(min=1.3705, central=1.378, max=1.3855, unit='µm'), resolution=2000, calibration=<calibration.reflectance>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 15,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C04'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Bidirectional Reflectance',\n",
" 'standard_name': 'toa_bidirectional_reflectance',\n",
" 'sensor_band_bit_depth': 10,\n",
" 'units': '%',\n",
" 'resolution': 1000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C05',\n",
" 'wavelength': WavelengthRange(min=1.58, central=1.61, max=1.64, unit='µm'),\n",
" 'calibration': 'reflectance',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 46, 500000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C05', wavelength=WavelengthRange(min=1.58, central=1.61, max=1.64, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 16,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C05'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Bidirectional Reflectance',\n",
" 'standard_name': 'toa_bidirectional_reflectance',\n",
" 'sensor_band_bit_depth': 10,\n",
" 'units': '%',\n",
" 'resolution': 2000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C06',\n",
" 'wavelength': WavelengthRange(min=2.225, central=2.25, max=2.275, unit='µm'),\n",
" 'calibration': 'reflectance',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 47, 100000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C06', wavelength=WavelengthRange(min=2.225, central=2.25, max=2.275, unit='µm'), resolution=2000, calibration=<calibration.reflectance>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 17,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C06'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Brightness Temperature',\n",
" 'standard_name': 'toa_brightness_temperature',\n",
" 'sensor_band_bit_depth': 14,\n",
" 'units': 'K',\n",
" 'resolution': 2000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C07',\n",
" 'wavelength': WavelengthRange(min=3.8, central=3.9, max=4.0, unit='µm'),\n",
" 'calibration': 'brightness_temperature',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 47, 700000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C07', wavelength=WavelengthRange(min=3.8, central=3.9, max=4.0, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 18,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C07'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Brightness Temperature',\n",
" 'standard_name': 'toa_brightness_temperature',\n",
" 'sensor_band_bit_depth': 12,\n",
" 'units': 'K',\n",
" 'resolution': 2000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C08',\n",
" 'wavelength': WavelengthRange(min=5.77, central=6.185, max=6.6, unit='µm'),\n",
" 'calibration': 'brightness_temperature',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 46, 500000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C08', wavelength=WavelengthRange(min=5.77, central=6.185, max=6.6, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 19,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C08'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Brightness Temperature',\n",
" 'standard_name': 'toa_brightness_temperature',\n",
" 'sensor_band_bit_depth': 11,\n",
" 'units': 'K',\n",
" 'resolution': 2000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C09',\n",
" 'wavelength': WavelengthRange(min=6.75, central=6.95, max=7.15, unit='µm'),\n",
" 'calibration': 'brightness_temperature',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 47, 100000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C09', wavelength=WavelengthRange(min=6.75, central=6.95, max=7.15, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 20,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C09'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Brightness Temperature',\n",
" 'standard_name': 'toa_brightness_temperature',\n",
" 'sensor_band_bit_depth': 12,\n",
" 'units': 'K',\n",
" 'resolution': 2000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C10',\n",
" 'wavelength': WavelengthRange(min=7.24, central=7.34, max=7.44, unit='µm'),\n",
" 'calibration': 'brightness_temperature',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 47, 700000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C10', wavelength=WavelengthRange(min=7.24, central=7.34, max=7.44, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 21,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C10'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 22,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Brightness Temperature',\n",
" 'standard_name': 'toa_brightness_temperature',\n",
" 'sensor_band_bit_depth': 12,\n",
" 'units': 'K',\n",
" 'resolution': 2000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C11',\n",
" 'wavelength': WavelengthRange(min=8.3, central=8.5, max=8.7, unit='µm'),\n",
" 'calibration': 'brightness_temperature',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 46, 500000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C11', wavelength=WavelengthRange(min=8.3, central=8.5, max=8.7, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 22,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C11'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 23,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Brightness Temperature',\n",
" 'standard_name': 'toa_brightness_temperature',\n",
" 'sensor_band_bit_depth': 11,\n",
" 'units': 'K',\n",
" 'resolution': 2000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C12',\n",
" 'wavelength': WavelengthRange(min=9.42, central=9.61, max=9.8, unit='µm'),\n",
" 'calibration': 'brightness_temperature',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 47, 100000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C12', wavelength=WavelengthRange(min=9.42, central=9.61, max=9.8, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 23,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C12'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 24,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Brightness Temperature',\n",
" 'standard_name': 'toa_brightness_temperature',\n",
" 'sensor_band_bit_depth': 12,\n",
" 'units': 'K',\n",
" 'resolution': 2000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C13',\n",
" 'wavelength': WavelengthRange(min=10.1, central=10.35, max=10.6, unit='µm'),\n",
" 'calibration': 'brightness_temperature',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 47, 700000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C13', wavelength=WavelengthRange(min=10.1, central=10.35, max=10.6, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 24,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C13'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Brightness Temperature',\n",
" 'standard_name': 'toa_brightness_temperature',\n",
" 'sensor_band_bit_depth': 12,\n",
" 'units': 'K',\n",
" 'resolution': 2000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C14',\n",
" 'wavelength': WavelengthRange(min=10.8, central=11.2, max=11.6, unit='µm'),\n",
" 'calibration': 'brightness_temperature',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 46, 500000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C14', wavelength=WavelengthRange(min=10.8, central=11.2, max=11.6, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 25,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C14'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 26,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Brightness Temperature',\n",
" 'standard_name': 'toa_brightness_temperature',\n",
" 'sensor_band_bit_depth': 12,\n",
" 'units': 'K',\n",
" 'resolution': 2000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C15',\n",
" 'wavelength': WavelengthRange(min=11.8, central=12.3, max=12.8, unit='µm'),\n",
" 'calibration': 'brightness_temperature',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 47, 100000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C15', wavelength=WavelengthRange(min=11.8, central=12.3, max=12.8, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 26,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C15'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 27,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"{'orbital_parameters': {'projection_longitude': -75.0,\n",
" 'projection_latitude': 0.0,\n",
" 'projection_altitude': 35786023.0,\n",
" 'satellite_nominal_latitude': 0.0,\n",
" 'satellite_nominal_longitude': -75.19999694824219,\n",
" 'satellite_nominal_altitude': 35786023.4375,\n",
" 'yaw_flip': False},\n",
" 'long_name': 'Brightness Temperature',\n",
" 'standard_name': 'toa_brightness_temperature',\n",
" 'sensor_band_bit_depth': 10,\n",
" 'units': 'K',\n",
" 'resolution': 2000,\n",
" 'grid_mapping': 'goes_imager_projection',\n",
" 'cell_methods': 't: point area: point',\n",
" 'platform_name': 'GOES-16',\n",
" 'sensor': 'abi',\n",
" 'name': 'C16',\n",
" 'wavelength': WavelengthRange(min=13.0, central=13.3, max=13.6, unit='µm'),\n",
" 'calibration': 'brightness_temperature',\n",
" 'modifiers': (),\n",
" 'observation_type': 'Rad',\n",
" 'scene_abbr': 'F',\n",
" 'scan_mode': 'M6',\n",
" 'platform_shortname': 'G16',\n",
" 'scene_id': 'Full Disk',\n",
" 'orbital_slot': 'GOES-East',\n",
" 'instrument_ID': 'FM1',\n",
" 'production_site': 'RBU',\n",
" 'timeline_ID': None,\n",
" 'start_time': datetime.datetime(2020, 9, 24, 0, 0, 15, 700000),\n",
" 'end_time': datetime.datetime(2020, 9, 24, 0, 9, 47, 700000),\n",
" 'area': 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),\n",
" '_satpy_id': DataID(name='C16', wavelength=WavelengthRange(min=13.0, central=13.3, max=13.6, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=()),\n",
" 'ancillary_variables': []}"
]
},
"execution_count": 27,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scn['C16'].attrs"
]
},
{
"cell_type": "code",
"execution_count": 28,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Setting 'PPP_CONFIG_DIR' to '/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/etc'\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C15', wavelength=WavelengthRange(min=11.8, central=12.3, max=12.8, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : pyresample.geometry] Projections for data and slice areas are identical: PROJCRS[\"unknown\",BASEGEOGCRS[\"unknown\",DATUM[\"unknown\",ELLIPSOID[\"GRS 1980\",6378137,298.257222096042,LENGTHUNIT[\"metre\",1,ID[\"EPSG\",9001]]]],PRIMEM[\"Greenwich\",0,ANGLEUNIT[\"degree\",0.0174532925199433],ID[\"EPSG\",8901]]],CONVERSION[\"unknown\",METHOD[\"Geostationary Satellite (Sweep X)\"],PARAMETER[\"Longitude of natural origin\",-75,ANGLEUNIT[\"degree\",0.0174532925199433],ID[\"EPSG\",8802]],PARAMETER[\"Satellite Height\",35786023,LENGTHUNIT[\"metre\",1,ID[\"EPSG\",9001]]],PARAMETER[\"False easting\",0,LENGTHUNIT[\"metre\",1],ID[\"EPSG\",8806]],PARAMETER[\"False northing\",0,LENGTHUNIT[\"metre\",1],ID[\"EPSG\",8807]]],CS[Cartesian,2],AXIS[\"(E)\",east,ORDER[1],LENGTHUNIT[\"metre\",1,ID[\"EPSG\",9001]]],AXIS[\"(N)\",north,ORDER[2],LENGTHUNIT[\"metre\",1,ID[\"EPSG\",9001]]]]\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[INFO: 2020-09-23 22:36:19 : satpy.resample] Using default KDTree resampler\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/numpy/lib/function_base.py:1280: RuntimeWarning: invalid value encountered in subtract\n",
" a = op(a[slice1], a[slice2])\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/numpy/lib/function_base.py:1280: RuntimeWarning: invalid value encountered in subtract\n",
" a = op(a[slice1], a[slice2])\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C16', wavelength=WavelengthRange(min=13.0, central=13.3, max=13.6, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C14', wavelength=WavelengthRange(min=10.8, central=11.2, max=11.6, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C04', wavelength=WavelengthRange(min=1.3705, central=1.378, max=1.3855, unit='µm'), resolution=2000, calibration=<calibration.reflectance>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C12', wavelength=WavelengthRange(min=9.42, central=9.61, max=9.8, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C10', wavelength=WavelengthRange(min=7.24, central=7.34, max=7.44, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C02', wavelength=WavelengthRange(min=0.59, central=0.64, max=0.69, unit='µm'), resolution=500, calibration=<calibration.reflectance>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : pyresample.geometry] Projections for data and slice areas are identical: PROJCRS[\"unknown\",BASEGEOGCRS[\"unknown\",DATUM[\"unknown\",ELLIPSOID[\"GRS 1980\",6378137,298.257222096042,LENGTHUNIT[\"metre\",1,ID[\"EPSG\",9001]]]],PRIMEM[\"Greenwich\",0,ANGLEUNIT[\"degree\",0.0174532925199433],ID[\"EPSG\",8901]]],CONVERSION[\"unknown\",METHOD[\"Geostationary Satellite (Sweep X)\"],PARAMETER[\"Longitude of natural origin\",-75,ANGLEUNIT[\"degree\",0.0174532925199433],ID[\"EPSG\",8802]],PARAMETER[\"Satellite Height\",35786023,LENGTHUNIT[\"metre\",1,ID[\"EPSG\",9001]]],PARAMETER[\"False easting\",0,LENGTHUNIT[\"metre\",1],ID[\"EPSG\",8806]],PARAMETER[\"False northing\",0,LENGTHUNIT[\"metre\",1],ID[\"EPSG\",8807]]],CS[Cartesian,2],AXIS[\"(E)\",east,ORDER[1],LENGTHUNIT[\"metre\",1,ID[\"EPSG\",9001]]],AXIS[\"(N)\",north,ORDER[2],LENGTHUNIT[\"metre\",1,ID[\"EPSG\",9001]]]]\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[INFO: 2020-09-23 22:36:19 : satpy.resample] Using default KDTree resampler\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/numpy/lib/function_base.py:1280: RuntimeWarning: invalid value encountered in subtract\n",
" a = op(a[slice1], a[slice2])\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/numpy/lib/function_base.py:1280: RuntimeWarning: invalid value encountered in subtract\n",
" a = op(a[slice1], a[slice2])\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-b8a59214b6cb72587530e735a4b8938ac275ec22.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C03', wavelength=WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : pyresample.geometry] Projections for data and slice areas are identical: PROJCRS[\"unknown\",BASEGEOGCRS[\"unknown\",DATUM[\"unknown\",ELLIPSOID[\"GRS 1980\",6378137,298.257222096042,LENGTHUNIT[\"metre\",1,ID[\"EPSG\",9001]]]],PRIMEM[\"Greenwich\",0,ANGLEUNIT[\"degree\",0.0174532925199433],ID[\"EPSG\",8901]]],CONVERSION[\"unknown\",METHOD[\"Geostationary Satellite (Sweep X)\"],PARAMETER[\"Longitude of natural origin\",-75,ANGLEUNIT[\"degree\",0.0174532925199433],ID[\"EPSG\",8802]],PARAMETER[\"Satellite Height\",35786023,LENGTHUNIT[\"metre\",1,ID[\"EPSG\",9001]]],PARAMETER[\"False easting\",0,LENGTHUNIT[\"metre\",1],ID[\"EPSG\",8806]],PARAMETER[\"False northing\",0,LENGTHUNIT[\"metre\",1],ID[\"EPSG\",8807]]],CS[Cartesian,2],AXIS[\"(E)\",east,ORDER[1],LENGTHUNIT[\"metre\",1,ID[\"EPSG\",9001]]],AXIS[\"(N)\",north,ORDER[2],LENGTHUNIT[\"metre\",1,ID[\"EPSG\",9001]]]]\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/pyproj/crs/crs.py:543: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems\n",
" proj_string = self.to_proj4()\n",
"[INFO: 2020-09-23 22:36:19 : satpy.resample] Using default KDTree resampler\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/numpy/lib/function_base.py:1280: RuntimeWarning: invalid value encountered in subtract\n",
" a = op(a[slice1], a[slice2])\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/numpy/lib/function_base.py:1280: RuntimeWarning: invalid value encountered in subtract\n",
" a = op(a[slice1], a[slice2])\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-d0da6cd853a15c6916288d0d2455e9abbd0ec6aa.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C06', wavelength=WavelengthRange(min=2.225, central=2.25, max=2.275, unit='µm'), resolution=2000, calibration=<calibration.reflectance>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C08', wavelength=WavelengthRange(min=5.77, central=6.185, max=6.6, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C11', wavelength=WavelengthRange(min=8.3, central=8.5, max=8.7, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-d0da6cd853a15c6916288d0d2455e9abbd0ec6aa.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C09', wavelength=WavelengthRange(min=6.75, central=6.95, max=7.15, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C05', wavelength=WavelengthRange(min=1.58, central=1.61, max=1.64, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-d0da6cd853a15c6916288d0d2455e9abbd0ec6aa.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C07', wavelength=WavelengthRange(min=3.8, central=3.9, max=4.0, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C13', wavelength=WavelengthRange(min=10.1, central=10.35, max=10.6, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=())\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.scene] Resampling DataID(name='C03', wavelength=WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Check if ./resample_lut-d0da6cd853a15c6916288d0d2455e9abbd0ec6aa.npz exists\n",
"[DEBUG: 2020-09-23 22:36:19 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='C07', wavelength=WavelengthRange(min=3.8, central=3.9, max=4.0, unit='µm'), resolution=2000, calibration=<calibration.brightness_temperature>, modifiers=('nir_reflectance',))\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling mul-f9c25de35a85e6f5dd6419bdad781ac1\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='C02', wavelength=WavelengthRange(min=0.59, central=0.64, max=0.69, unit='µm'), resolution=500, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-b8a59214b6cb72587530e735a4b8938ac275ec22.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='C02', wavelength=WavelengthRange(min=0.59, central=0.64, max=0.69, unit='µm'), resolution=500, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected'))\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-b8a59214b6cb72587530e735a4b8938ac275ec22.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-d0da6cd853a15c6916288d0d2455e9abbd0ec6aa.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='night_ir_alpha', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling concatenate-b1b5a83d0132c13318f5721ff4166a27\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='_night_background_hires')\n",
"[INFO: 2020-09-23 22:36:20 : satpy.scene] Not reducing data before resampling.\n",
"[INFO: 2020-09-23 22:36:20 : satpy.resample] Using default KDTree resampler\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/numpy/lib/function_base.py:1280: RuntimeWarning: invalid value encountered in subtract\n",
" a = op(a[slice1], a[slice2])\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-29bca24069d8c72778e6c10048b9d130bb1eeced.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling None\n",
"[WARNING: 2020-09-23 22:36:20 : pyresample.kd_tree] Fill value incompatible with integer data using 255 instead.\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='_cira_day_convection_dep_0', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='C03', wavelength=WavelengthRange(min=0.8455, central=0.865, max=0.8845, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected_crefl'))\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-d0da6cd853a15c6916288d0d2455e9abbd0ec6aa.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='night_fog', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling where-7933d2ba747ce384fe491459b3f91182\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='_cira_day_convection_dep_1', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='ir_cloud_day', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling concatenate-73c873f1e2f3daa8a3dc93cb76e973ad\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='_convection_dep_1', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='C02', wavelength=WavelengthRange(min=0.59, central=0.64, max=0.69, unit='µm'), resolution=500, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected_crefl'))\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-b8a59214b6cb72587530e735a4b8938ac275ec22.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='dust', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling where-19337e31300e8829d14270b04dd57988\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='water_vapors2', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling where-b00248997362f4caf5286b1dec6c3b07\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='ir108_3d', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling truediv-bdb003f9ab4c1d1c755390beead41da9\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='C01', wavelength=WavelengthRange(min=0.45, central=0.47, max=0.49, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected', 'rayleigh_corrected_crefl'))\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-d0da6cd853a15c6916288d0d2455e9abbd0ec6aa.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='_night_background')\n",
"[INFO: 2020-09-23 22:36:20 : satpy.scene] Not reducing data before resampling.\n",
"[INFO: 2020-09-23 22:36:20 : satpy.resample] Using default KDTree resampler\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/numpy/lib/function_base.py:1280: RuntimeWarning: invalid value encountered in subtract\n",
" a = op(a[slice1], a[slice2])\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-5cec9b560e7c4efdc84e7888af57aeb351420a84.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling None\n",
"[WARNING: 2020-09-23 22:36:20 : pyresample.kd_tree] Fill value incompatible with integer data using 255 instead.\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='ash', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling where-19337e31300e8829d14270b04dd57988\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='night_microphysics', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling where-7933d2ba747ce384fe491459b3f91182\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='colorized_ir_clouds', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling truediv-fe9ea93b1eeeede061fc2e690dc32e6c\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.scene] Resampling DataID(name='_snow_fog_dep_0', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:20 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.scene] Resampling DataID(name='C05', wavelength=WavelengthRange(min=1.58, central=1.61, max=1.64, unit='µm'), resolution=1000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Check if ./resample_lut-d0da6cd853a15c6916288d0d2455e9abbd0ec6aa.npz exists\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.scene] Resampling DataID(name='tropical_airmass', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Resampling where-a8812975960d0e3a76cc2b0b3a64f62d\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.scene] Resampling DataID(name='airmass', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Resampling where-442264f292323db323bec1cc7b51ed8b\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"[DEBUG: 2020-09-23 22:36:21 : satpy.scene] Resampling DataID(name='C06', wavelength=WavelengthRange(min=2.225, central=2.25, max=2.275, unit='µm'), resolution=2000, calibration=<calibration.reflectance>, modifiers=('sunz_corrected',))\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Resampling None\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.scene] Resampling DataID(name='so2', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Resampling where-00635709537d455ce9d904675b30cea5\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.scene] Resampling DataID(name='_convection_dep_0', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Resampling Rad\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.scene] Resampling DataID(name='cloudtop', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Resampling where-4280cbb20ca6602c4a61daf2970ef208\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.scene] Resampling DataID(name='night_microphysics_abi', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Resampling where-e439e0c2fd5dd4f867684e6bfed8059b\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.scene] Resampling DataID(name='water_vapors1', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Resampling where-d65e77a653eda425a6ca9a692c84d4df\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.scene] Resampling DataID(name='fog', resolution=2000)\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Check if ./resample_lut-a2170994bfac1c6e8cae83fa7084cbe4f7109a7c.npz exists\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Computing kd-tree parameters\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.resample] Resampling where-2355d8e9dad8ddf9547f69836209b20a\n",
"/home/ray/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/node.py:105: UserWarning: Attribute access to DataIDs is deprecated, use key access instead.\n",
" return self.name == other.name\n",
"[INFO: 2020-09-23 22:36:21 : satpy.composites] Removing Rayleigh scattering with atmosphere 'us-standard' and aerosol type 'rayleigh_only' for 'C01'\n",
"[INFO: 2020-09-23 22:36:21 : pyspectral.rayleigh] Atmosphere chosen: us-standard\n",
"[DEBUG: 2020-09-23 22:36:21 : pyspectral.rayleigh] LUT filename: /home/ray/.local/share/pyspectral/rayleigh_only/rayleigh_lut_us-standard.h5\n",
"[DEBUG: 2020-09-23 22:36:21 : pyspectral.rsr_reader] Filename: /home/ray/.local/share/pyspectral/rsr_abi_GOES-16.h5\n",
"[DEBUG: 2020-09-23 22:36:21 : pyspectral.rsr_reader] Filename: /home/ray/.local/share/pyspectral/rsr_abi_GOES-16.h5\n",
"[DEBUG: 2020-09-23 22:36:21 : pyspectral.rsr_reader] No detectors found - assume only one...\n",
"[DEBUG: 2020-09-23 22:36:21 : pyspectral.rayleigh] Band name: ch1 Effective wavelength: 0.469946\n",
"[DEBUG: 2020-09-23 22:36:21 : pyspectral.rayleigh] Time - Interpolation: 0.047816\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.composites] Sharpening image with high resolution red band\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.writers] Enhancement configuration options: [{'name': 'stretch', 'method': <function stretch at 0x7f34685f7430>, 'kwargs': {'stretch': 'linear', 'cutoffs': [0.02, 0.02]}}, {'name': 'inverse', 'method': <function invert at 0x7f34685f7af0>, 'args': [[True, True, True, True]]}]\n",
"[DEBUG: 2020-09-23 22:36:21 : trollimage.xrimage] Applying stretch linear with parameters {'cutoffs': [0.02, 0.02]}\n",
"[DEBUG: 2020-09-23 22:36:21 : trollimage.xrimage] Perform a linear contrast stretch.\n",
"[DEBUG: 2020-09-23 22:36:21 : trollimage.xrimage] Calculate the histogram quantiles: \n",
"[DEBUG: 2020-09-23 22:36:21 : trollimage.xrimage] Left and right quantiles: 0.02 0.02\n",
"[DEBUG: 2020-09-23 22:36:21 : trollimage.xrimage] Applying invert with parameters [True, True, True, True]\n",
"[DEBUG: 2020-09-23 22:36:21 : satpy.writers] Enhancement configuration options: [{'name': 'stretch', 'method': <function stretch at 0x7f34685f7430>, 'kwargs': {'stretch': 'linear'}}]\n",
"[DEBUG: 2020-09-23 22:36:21 : trollimage.xrimage] Applying stretch linear with parameters {}\n",
"[DEBUG: 2020-09-23 22:36:21 : trollimage.xrimage] Perform a linear contrast stretch.\n",
"[DEBUG: 2020-09-23 22:36:21 : trollimage.xrimage] Calculate the histogram quantiles: \n",
"[DEBUG: 2020-09-23 22:36:21 : trollimage.xrimage] Left and right quantiles: 0.005 0.005\n"
]
},
{
"ename": "ValueError",
"evalue": "operands could not be broadcast together with shapes (4,) (0,) ",
"output_type": "error",
"traceback": [
"\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
"\u001b[0;31mValueError\u001b[0m Traceback (most recent call last)",
"\u001b[0;32m<ipython-input-28-d99dafe80c78>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mscn_reampled\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mscn\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mresample\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mscn\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'C01'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mattrs\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'area'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/scene.py\u001b[0m in \u001b[0;36mresample\u001b[0;34m(self, destination, datasets, generate, unload, resampler, reduce_data, **resample_kwargs)\u001b[0m\n\u001b[1;32m 1126\u001b[0m \u001b[0;31m# multiple resolutions, etc.)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 1127\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mgenerate\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1128\u001b[0;31m \u001b[0mkeepables\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mnew_scn\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_composites\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 1129\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 1130\u001b[0m \u001b[0;31m# don't lose datasets that we may need later for generating\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/scene.py\u001b[0m in \u001b[0;36m_generate_composites\u001b[0;34m(self, nodes)\u001b[0m\n\u001b[1;32m 895\u001b[0m \u001b[0mrequired_nodes\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_wishlist\u001b[0m \u001b[0;34m-\u001b[0m \u001b[0mset\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_datasets\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mkeys\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 896\u001b[0m \u001b[0mnodes\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mset\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_dependency_tree\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mtrunk\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mnodes\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mrequired_nodes\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m-\u001b[0m \u001b[0mset\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_datasets\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mkeys\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 897\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_read_composites\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mnodes\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 898\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 899\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0m_remove_failed_datasets\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkeepables\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/scene.py\u001b[0m in \u001b[0;36m_read_composites\u001b[0;34m(self, compositor_nodes)\u001b[0m\n\u001b[1;32m 871\u001b[0m \u001b[0mkeepables\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mset\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 872\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mitem\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mcompositor_nodes\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 873\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_composite\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mitem\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkeepables\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 874\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mkeepables\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 875\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/scene.py\u001b[0m in \u001b[0;36m_generate_composite\u001b[0;34m(self, comp_node, keepables)\u001b[0m\n\u001b[1;32m 811\u001b[0m \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 812\u001b[0m \u001b[0mdelayed_prereq\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;32mFalse\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 813\u001b[0;31m prereq_datasets = self._get_prereq_datasets(\n\u001b[0m\u001b[1;32m 814\u001b[0m \u001b[0mcomp_node\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 815\u001b[0m \u001b[0mprereqs\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/scene.py\u001b[0m in \u001b[0;36m_get_prereq_datasets\u001b[0;34m(self, comp_id, prereq_nodes, keepables, skip)\u001b[0m\n\u001b[1;32m 761\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mprereq_id\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_datasets\u001b[0m \u001b[0;32mand\u001b[0m \u001b[0mprereq_id\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mkeepables\u001b[0m\u001b[0;31m \u001b[0m\u001b[0;31m\\\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 762\u001b[0m \u001b[0;32mand\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0mprereq_node\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mis_leaf\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 763\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_composite\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mprereq_node\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkeepables\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 764\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 765\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mprereq_node\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_dependency_tree\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mempty_node\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/scene.py\u001b[0m in \u001b[0;36m_generate_composite\u001b[0;34m(self, comp_node, keepables)\u001b[0m\n\u001b[1;32m 845\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 846\u001b[0m \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 847\u001b[0;31m composite = compositor(prereq_datasets,\n\u001b[0m\u001b[1;32m 848\u001b[0m \u001b[0moptional_datasets\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0moptional_datasets\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 849\u001b[0m **self.attrs)\n",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/composites/__init__.py\u001b[0m in \u001b[0;36m__call__\u001b[0;34m(self, projectables, *args, **kwargs)\u001b[0m\n\u001b[1;32m 1654\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 1655\u001b[0m \u001b[0;31m# Get merged metadata\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1656\u001b[0;31m \u001b[0mattrs\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mcombine_metadata\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mforeground\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mbackground\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 1657\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mattrs\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'sensor'\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 1658\u001b[0m \u001b[0;31m# sensor can be a set\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/dataset/metadata.py\u001b[0m in \u001b[0;36mcombine_metadata\u001b[0;34m(average_times, *metadata_objects)\u001b[0m\n\u001b[1;32m 50\u001b[0m \u001b[0mshared_keys\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0m_shared_keys\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0minfo_dicts\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 51\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 52\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0m_combine_shared_info\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mshared_keys\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0minfo_dicts\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0maverage_times\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 53\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 54\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/dataset/metadata.py\u001b[0m in \u001b[0;36m_combine_shared_info\u001b[0;34m(shared_keys, info_dicts, average_times)\u001b[0m\n\u001b[1;32m 78\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0;34m'time'\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mkey\u001b[0m \u001b[0;32mand\u001b[0m \u001b[0misinstance\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mvalues\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdatetime\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32mand\u001b[0m \u001b[0maverage_times\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 79\u001b[0m \u001b[0mshared_info\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mkey\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0maverage_datetimes\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mvalues\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 80\u001b[0;31m \u001b[0;32melif\u001b[0m \u001b[0m_are_values_combinable\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mvalues\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 81\u001b[0m \u001b[0mshared_info\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mkey\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mvalues\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 82\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mshared_info\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/dataset/metadata.py\u001b[0m in \u001b[0;36m_are_values_combinable\u001b[0;34m(values)\u001b[0m\n\u001b[1;32m 113\u001b[0m \u001b[0;31m# checked any_arrays so this false positive should have no impact\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 114\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0m_all_list_of_arrays_equal\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mvalues\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 115\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0m_all_values_equal\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mvalues\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 116\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 117\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/dataset/metadata.py\u001b[0m in \u001b[0;36m_all_values_equal\u001b[0;34m(values)\u001b[0m\n\u001b[1;32m 176\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0m_all_values_equal\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mvalues\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 177\u001b[0m \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 178\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0m_pairwise_all\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mnan_allclose\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvalues\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 179\u001b[0m \u001b[0;32mexcept\u001b[0m \u001b[0mTypeError\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 180\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0m_pairwise_all\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0meq\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvalues\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/satpy/dataset/metadata.py\u001b[0m in \u001b[0;36m_pairwise_all\u001b[0;34m(func, values)\u001b[0m\n\u001b[1;32m 141\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0m_pairwise_all\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfunc\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvalues\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 142\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mvalue\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mvalues\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m1\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 143\u001b[0;31m \u001b[0;32mif\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0mfunc\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mvalues\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvalue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 144\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0;32mFalse\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 145\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0;32mTrue\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m<__array_function__ internals>\u001b[0m in \u001b[0;36mallclose\u001b[0;34m(*args, **kwargs)\u001b[0m\n",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/numpy/core/numeric.py\u001b[0m in \u001b[0;36mallclose\u001b[0;34m(a, b, rtol, atol, equal_nan)\u001b[0m\n\u001b[1;32m 2187\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 2188\u001b[0m \"\"\"\n\u001b[0;32m-> 2189\u001b[0;31m \u001b[0mres\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mall\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0misclose\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0ma\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mb\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mrtol\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mrtol\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0matol\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0matol\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mequal_nan\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mequal_nan\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 2190\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mbool\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mres\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 2191\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m<__array_function__ internals>\u001b[0m in \u001b[0;36misclose\u001b[0;34m(*args, **kwargs)\u001b[0m\n",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/numpy/core/numeric.py\u001b[0m in \u001b[0;36misclose\u001b[0;34m(a, b, rtol, atol, equal_nan)\u001b[0m\n\u001b[1;32m 2288\u001b[0m \u001b[0myfin\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0misfinite\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0my\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 2289\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mxfin\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32mand\u001b[0m \u001b[0mall\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0myfin\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 2290\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0mwithin_tol\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0my\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0matol\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mrtol\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 2291\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 2292\u001b[0m \u001b[0mfinite\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mxfin\u001b[0m \u001b[0;34m&\u001b[0m \u001b[0myfin\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m~/local/bin/anaconda3/envs/test_env/lib/python3.8/site-packages/numpy/core/numeric.py\u001b[0m in \u001b[0;36mwithin_tol\u001b[0;34m(x, y, atol, rtol)\u001b[0m\n\u001b[1;32m 2274\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mwithin_tol\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0my\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0matol\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mrtol\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 2275\u001b[0m \u001b[0;32mwith\u001b[0m \u001b[0merrstate\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0minvalid\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'ignore'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 2276\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0mless_equal\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mabs\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m-\u001b[0m\u001b[0my\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0matol\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0mrtol\u001b[0m \u001b[0;34m*\u001b[0m \u001b[0mabs\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0my\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 2277\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 2278\u001b[0m \u001b[0mx\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0masanyarray\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0ma\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;31mValueError\u001b[0m: operands could not be broadcast together with shapes (4,) (0,) "
]
}
],
"source": [
"scn_reampled = scn.resample(scn['C01'].attrs['area'])"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"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.8.3"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment