Skip to content

Instantly share code, notes, and snippets.

@rutj3

rutj3/Satpy_OLCI.ipynb

Created March 7, 2018 10:54
Show Gist options
  • Save rutj3/ee45afe4daa5e08df10a2675cf7841c2 to your computer and use it in GitHub Desktop.
Save rutj3/ee45afe4daa5e08df10a2675cf7841c2 to your computer and use it in GitHub Desktop.
Download ZIP
Satpy OLCI
Display the source blob
Display the rendered blob
Raw
Satpy_OLCI.ipynb
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"from satpy.utils import debug_on; debug_on()\n",
"from satpy import Scene, find_files_and_readers\n",
"\n",
"from datetime import datetime"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"[DEBUG: 2018-03-07 11:54:20 : satpy.readers] Reading ['C:\\\\Miniconda3\\\\envs\\\\satpy\\\\lib\\\\site-packages\\\\satpy\\\\etc\\\\readers\\\\nc_olci_l1b.yaml', 'C:\\\\Miniconda3\\\\envs\\\\satpy\\\\lib\\\\site-packages\\\\satpy\\\\etc\\\\readers\\\\nc_olci_l1b.yaml', 'C:\\\\Miniconda3\\\\envs\\\\satpy\\\\lib\\\\site-packages\\\\satpy\\\\etc\\\\readers\\\\nc_olci_l1b.yaml']\n",
"C:\\Miniconda3\\envs\\satpy\\lib\\site-packages\\h5py\\__init__.py:36: FutureWarning: Conversion of the second argument of issubdtype from `float` to `np.floating` is deprecated. In future, it will be treated as `np.float64 == np.dtype(float).type`.\n",
" from ._conv import register_converters as _register_converters\n",
"[DEBUG: 2018-03-07 11:54:20 : satpy.scene] Setting 'PPP_CONFIG_DIR' to 'C:\\Miniconda3\\envs\\satpy\\lib\\site-packages\\satpy\\etc'\n",
"[DEBUG: 2018-03-07 11:54:20 : satpy.readers] Reading ['C:\\\\Miniconda3\\\\envs\\\\satpy\\\\lib\\\\site-packages\\\\satpy\\\\etc\\\\readers\\\\nc_olci_l1b.yaml', 'C:\\\\Miniconda3\\\\envs\\\\satpy\\\\lib\\\\site-packages\\\\satpy\\\\etc\\\\readers\\\\nc_olci_l1b.yaml', 'C:\\\\Miniconda3\\\\envs\\\\satpy\\\\lib\\\\site-packages\\\\satpy\\\\etc\\\\readers\\\\nc_olci_l1b.yaml']\n",
"[DEBUG: 2018-03-07 11:54:21 : satpy.readers.yaml_reader] Assigning to nc_olci_l1b: ['D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa06_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa09_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa16_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa19_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa11_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa17_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa15_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa10_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa18_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa02_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa04_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa15_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa05_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa13_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa21_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa08_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa14_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa05_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa09_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa17_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa18_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa16_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa20_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa12_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa13_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa19_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa03_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa10_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa21_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa08_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa11_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa20_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa12_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa07_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa03_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa07_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa02_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa14_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa06_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa01_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\Oa04_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\Oa01_radiance.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\tie_geometries.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\tie_geometries.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124529_20171011T124829_20171012T172726_0179_023_152_2340_LN1_O_NT_002.SEN3\\\\geo_coordinates.nc', 'D:\\\\Data\\\\Sentinel\\\\Sentinel3\\\\OLCI\\\\S3A_OL_1_EFR____20171011T124829_20171011T125129_20171012T172749_0179_023_152_2520_LN1_O_NT_002.SEN3\\\\geo_coordinates.nc']\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"[DEBUG: 2018-03-07 11:54:21 : satpy.composites] Looking for composites config file olci.yaml\n",
"[DEBUG: 2018-03-07 11:54:21 : satpy.composites] Looking for composites config file visir.yaml\n"
]
}
],
"source": [
"base_dir = r'D:\\Data\\Sentinel\\Sentinel3\\OLCI'\n",
"\n",
"files = find_files_and_readers(sensor=\"olci\",\n",
" start_time=datetime(2017, 10, 11, 12, 0),\n",
" end_time=datetime(2017, 10, 11, 12, 59),\n",
" base_dir=base_dir,\n",
" reader='nc_olci_l1b')\n",
"\n",
"\n",
"scn = Scene(filenames=files)"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"scrolled": false
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"[DEBUG: 2018-03-07 11:54:25 : satpy.readers.yaml_reader] No coordinates found for DatasetID(name='longitude', wavelength=None, resolution=300, polarization=None, calibration=None, modifiers=())\n",
"[DEBUG: 2018-03-07 11:54:25 : satpy.readers.nc_olci] Reading longitude.\n",
"[DEBUG: 2018-03-07 11:54:26 : satpy.readers.nc_olci] Reading longitude.\n",
"[DEBUG: 2018-03-07 11:54:26 : satpy.readers.yaml_reader] No coordinates found for DatasetID(name='latitude', wavelength=None, resolution=300, polarization=None, calibration=None, modifiers=())\n",
"[DEBUG: 2018-03-07 11:54:26 : satpy.readers.nc_olci] Reading latitude.\n",
"[DEBUG: 2018-03-07 11:54:27 : satpy.readers.nc_olci] Reading latitude.\n",
"[DEBUG: 2018-03-07 11:54:28 : satpy.readers.nc_olci] Reading Oa08.\n",
"[DEBUG: 2018-03-07 11:54:29 : satpy.readers.nc_olci] Reading Oa08.\n",
"[DEBUG: 2018-03-07 11:54:30 : satpy.readers.nc_olci] Reading Oa03.\n",
"[DEBUG: 2018-03-07 11:54:31 : satpy.readers.nc_olci] Reading Oa03.\n",
"[DEBUG: 2018-03-07 11:54:32 : satpy.readers.nc_olci] Reading solar_zenith_angle.\n",
"[DEBUG: 2018-03-07 11:54:36 : satpy.readers.nc_olci] Reading solar_zenith_angle.\n",
"[DEBUG: 2018-03-07 11:54:40 : satpy.readers.nc_olci] Reading satellite_azimuth_angle.\n",
"[DEBUG: 2018-03-07 11:54:44 : satpy.readers.nc_olci] Reading satellite_azimuth_angle.\n",
"[DEBUG: 2018-03-07 11:54:48 : satpy.readers.nc_olci] Reading satellite_zenith_angle.\n",
"[DEBUG: 2018-03-07 11:54:48 : satpy.readers.nc_olci] Reading satellite_zenith_angle.\n",
"[DEBUG: 2018-03-07 11:54:48 : satpy.readers.nc_olci] Reading Oa06.\n",
"[DEBUG: 2018-03-07 11:54:49 : satpy.readers.nc_olci] Reading Oa06.\n",
"[DEBUG: 2018-03-07 11:54:50 : satpy.readers.nc_olci] Reading solar_azimuth_angle.\n",
"[DEBUG: 2018-03-07 11:54:50 : satpy.readers.nc_olci] Reading solar_azimuth_angle.\n",
"[DEBUG: 2018-03-07 11:54:51 : satpy.composites] Applying sun zen correction\n",
"[DEBUG: 2018-03-07 11:54:51 : satpy.composites] Computing sun zenith angles.\n",
"[DEBUG: 2018-03-07 11:54:54 : satpy.composites] Apply the effective solar atmospheric path length correction method by Li and Shibata\n",
"[DEBUG: 2018-03-07 11:55:00 : satpy.composites] Sun-zenith correction applied. Computation time: 9.1 (sec)\n",
"[DEBUG: 2018-03-07 11:55:00 : satpy.composites] Applying sun zen correction\n",
"[DEBUG: 2018-03-07 11:55:00 : satpy.composites] Apply the standard sun-zenith correction [1/cos(sunz)]\n",
"[DEBUG: 2018-03-07 11:55:03 : satpy.composites] Sun-zenith correction applied. Computation time: 3.0 (sec)\n",
"[INFO: 2018-03-07 11:55:03 : satpy.composites] Removing Rayleigh scattering and aerosol absorption\n",
"[INFO: 2018-03-07 11:55:05 : pyspectral.rayleigh] Atmosphere chosen: us-standard\n",
"[DEBUG: 2018-03-07 11:55:05 : pyspectral.rayleigh] LUT filename: C:\\Users\\Rutger\\AppData\\Local\\pytroll\\pyspectral\\marine_clean_aerosol\\rayleigh_lut_us-standard.h5\n",
"[DEBUG: 2018-03-07 11:55:05 : pyspectral.rsr_reader] Filename: C:\\Users\\Rutger\\AppData\\Local\\pytroll\\pyspectral\\rsr_olci_Sentinel-3A.h5\n",
"[DEBUG: 2018-03-07 11:55:05 : pyspectral.rsr_reader] Filename: C:\\Users\\Rutger\\AppData\\Local\\pytroll\\pyspectral\\rsr_olci_Sentinel-3A.h5\n",
"[DEBUG: 2018-03-07 11:55:05 : pyspectral.rsr_reader] No detectors found - assume only one...\n",
"[WARNING: 2018-03-07 11:55:05 : satpy.composites] Could not get the reflectance correction using band name: Oa08\n",
"[WARNING: 2018-03-07 11:55:05 : satpy.composites] Will try use the wavelength, however, this may be ambiguous!\n",
"[DEBUG: 2018-03-07 11:55:05 : pyspectral.rsr_reader] Filename: C:\\Users\\Rutger\\AppData\\Local\\pytroll\\pyspectral\\rsr_olci_Sentinel-3A.h5\n",
"[DEBUG: 2018-03-07 11:55:05 : pyspectral.rsr_reader] Filename: C:\\Users\\Rutger\\AppData\\Local\\pytroll\\pyspectral\\rsr_olci_Sentinel-3A.h5\n",
"[DEBUG: 2018-03-07 11:55:05 : pyspectral.rsr_reader] No detectors found - assume only one...\n"
]
},
{
"ename": "AttributeError",
"evalue": "More than one band found with that wavelength! ['band=ch7', 'band=ch8', 'band=ch9', 'band=ch10', 'band=ch11', 'band=ch12', 'band=ch13', 'band=ch14']",
"output_type": "error",
"traceback": [
"\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
"\u001b[1;31mKeyError\u001b[0m Traceback (most recent call last)",
"\u001b[1;32mC:\\Miniconda3\\envs\\satpy\\lib\\site-packages\\satpy\\composites\\__init__.py\u001b[0m in \u001b[0;36m__call__\u001b[1;34m(self, projectables, optional_datasets, **info)\u001b[0m\n\u001b[0;32m 375\u001b[0m \u001b[1;32mtry\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 376\u001b[1;33m \u001b[0mrefl_cor_band\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mcorrector\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mget_reflectance\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0msunz\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msatz\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mssadiff\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mvis\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mid\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mname\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mred\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 377\u001b[0m \u001b[1;32mexcept\u001b[0m \u001b[0mKeyError\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32mC:\\Miniconda3\\envs\\satpy\\lib\\site-packages\\pyspectral\\rayleigh.py\u001b[0m in \u001b[0;36mget_reflectance\u001b[1;34m(self, sun_zenith, sat_zenith, azidiff, bandname, redband)\u001b[0m\n\u001b[0;32m 159\u001b[0m \u001b[1;31m# Get wavelength in nm for band:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 160\u001b[1;33m \u001b[0mwvl\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mget_effective_wavelength\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mbandname\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m*\u001b[0m \u001b[1;36m1000.0\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 161\u001b[0m \u001b[0mrayl\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mwvl_coord\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mazid_coord\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msatz_sec_coord\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msunz_sec_coord\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mget_reflectance_lut\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32mC:\\Miniconda3\\envs\\satpy\\lib\\site-packages\\pyspectral\\rayleigh.py\u001b[0m in \u001b[0;36mget_effective_wavelength\u001b[1;34m(self, bandname)\u001b[0m\n\u001b[0;32m 140\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 141\u001b[1;33m wvl, resp = rsr.rsr[bandname][\n\u001b[0m\u001b[0;32m 142\u001b[0m 'det-1']['wavelength'], rsr.rsr[bandname]['det-1']['response']\n",
"\u001b[1;31mKeyError\u001b[0m: 'ch8'",
"\nDuring handling of the above exception, another exception occurred:\n",
"\u001b[1;31mAttributeError\u001b[0m Traceback (most recent call last)",
"\u001b[1;32m<ipython-input-3-159cb4b8618b>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[0;32m 1\u001b[0m \u001b[0mcomposite\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;34m'true_color'\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 2\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 3\u001b[1;33m \u001b[0mscn\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mload\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mcomposite\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 4\u001b[0m \u001b[0mscn\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mshow\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mcomposite\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32mC:\\Miniconda3\\envs\\satpy\\lib\\site-packages\\satpy\\scene.py\u001b[0m in \u001b[0;36mload\u001b[1;34m(self, wishlist, calibration, resolution, polarization, compute, unload, **kwargs)\u001b[0m\n\u001b[0;32m 565\u001b[0m \u001b[0mkeepables\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;32mNone\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 566\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mcompute\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 567\u001b[1;33m \u001b[0mkeepables\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcompute\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 568\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmissing_datasets\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 569\u001b[0m \u001b[1;31m# copy the set of missing datasets because they won't be valid\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32mC:\\Miniconda3\\envs\\satpy\\lib\\site-packages\\satpy\\scene.py\u001b[0m in \u001b[0;36mcompute\u001b[1;34m(self, nodes)\u001b[0m\n\u001b[0;32m 512\u001b[0m \u001b[0mnodes\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mset\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdep_tree\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mtrunk\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mnodes\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mrequired_nodes\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m-\u001b[0m\u001b[0;31m \u001b[0m\u001b[0;31m\\\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 513\u001b[0m \u001b[0mset\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdatasets\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mkeys\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 514\u001b[1;33m \u001b[1;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mread_composites\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mnodes\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 515\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 516\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0m_remove_failed_datasets\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mkeepables\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32mC:\\Miniconda3\\envs\\satpy\\lib\\site-packages\\satpy\\scene.py\u001b[0m in \u001b[0;36mread_composites\u001b[1;34m(self, compositor_nodes)\u001b[0m\n\u001b[0;32m 486\u001b[0m \u001b[0mkeepables\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mset\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 487\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mitem\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mcompositor_nodes\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 488\u001b[1;33m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_generate_composite\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mitem\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mkeepables\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 489\u001b[0m \u001b[1;32mreturn\u001b[0m \u001b[0mkeepables\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 490\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32mC:\\Miniconda3\\envs\\satpy\\lib\\site-packages\\satpy\\scene.py\u001b[0m in \u001b[0;36m_generate_composite\u001b[1;34m(self, comp_node, keepables)\u001b[0m\n\u001b[0;32m 461\u001b[0m composite = compositor(prereq_datasets,\n\u001b[0;32m 462\u001b[0m \u001b[0moptional_datasets\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0moptional_datasets\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 463\u001b[1;33m **self.info)\n\u001b[0m\u001b[0;32m 464\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdatasets\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mcomposite\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mid\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mcomposite\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 465\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mcomp_node\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mname\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mwishlist\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32mC:\\Miniconda3\\envs\\satpy\\lib\\site-packages\\satpy\\composites\\__init__.py\u001b[0m in \u001b[0;36m__call__\u001b[1;34m(self, projectables, optional_datasets, **info)\u001b[0m\n\u001b[0;32m 379\u001b[0m \u001b[0mLOG\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mwarning\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"Will try use the wavelength, however, this may be ambiguous!\"\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 380\u001b[0m refl_cor_band = corrector.get_reflectance(sunz, satz, ssadiff,\n\u001b[1;32m--> 381\u001b[1;33m vis.id.wavelength[1], red)\n\u001b[0m\u001b[0;32m 382\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 383\u001b[0m proj = Dataset(vis - refl_cor_band,\n",
"\u001b[1;32mC:\\Miniconda3\\envs\\satpy\\lib\\site-packages\\pyspectral\\rayleigh.py\u001b[0m in \u001b[0;36mget_reflectance\u001b[1;34m(self, sun_zenith, sat_zenith, azidiff, bandname, redband)\u001b[0m\n\u001b[0;32m 158\u001b[0m \u001b[1;34m\"\"\"Get the reflectance from the three sun-sat angles.\"\"\"\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 159\u001b[0m \u001b[1;31m# Get wavelength in nm for band:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 160\u001b[1;33m \u001b[0mwvl\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mget_effective_wavelength\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mbandname\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m*\u001b[0m \u001b[1;36m1000.0\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 161\u001b[0m \u001b[0mrayl\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mwvl_coord\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mazid_coord\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msatz_sec_coord\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msunz_sec_coord\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mget_reflectance_lut\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 162\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32mC:\\Miniconda3\\envs\\satpy\\lib\\site-packages\\pyspectral\\rayleigh.py\u001b[0m in \u001b[0;36mget_effective_wavelength\u001b[1;34m(self, bandname)\u001b[0m\n\u001b[0;32m 137\u001b[0m raise BandFrequencyOutOfRange(\n\u001b[0;32m 138\u001b[0m 'Requested band frequency should be between 0.4 and 0.8 microns!')\n\u001b[1;32m--> 139\u001b[1;33m \u001b[0mbandname\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mget_bandname_from_wavelength\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msensor\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mbandname\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrsr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mrsr\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 140\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 141\u001b[0m wvl, resp = rsr.rsr[bandname][\n",
"\u001b[1;32mC:\\Miniconda3\\envs\\satpy\\lib\\site-packages\\pyspectral\\utils.py\u001b[0m in \u001b[0;36mget_bandname_from_wavelength\u001b[1;34m(sensor, wavelength, rsr, epsilon, multiple_bands)\u001b[0m\n\u001b[0;32m 283\u001b[0m \u001b[0mbstrlist\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[1;34m'band={}'\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mformat\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mb\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mb\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mchfound\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 284\u001b[0m \u001b[1;32mif\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mmultiple_bands\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 285\u001b[1;33m \u001b[1;32mraise\u001b[0m \u001b[0mAttributeError\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"More than one band found with that wavelength! {}\"\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mformat\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mstr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mbstrlist\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 286\u001b[0m \u001b[1;32melse\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 287\u001b[0m \u001b[0mLOG\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdebug\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"More than one band found with requested wavelength: %s\"\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mstr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mbstrlist\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;31mAttributeError\u001b[0m: More than one band found with that wavelength! ['band=ch7', 'band=ch8', 'band=ch9', 'band=ch10', 'band=ch11', 'band=ch12', 'band=ch13', 'band=ch14']"
]
}
],
"source": [
"composite = 'true_color'\n",
"\n",
"scn.load([composite])\n",
"scn.show(composite)"
]
},
{
"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.6.4"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment