Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Cannot load datasets of multiple SEVIRI native files #1656

Closed
harry9713 opened this issue Apr 30, 2021 · 11 comments · Fixed by #1663
Closed

Cannot load datasets of multiple SEVIRI native files #1656

harry9713 opened this issue Apr 30, 2021 · 11 comments · Fixed by #1663

Comments

@harry9713
Copy link

harry9713 commented Apr 30, 2021 •
edited

Describe the bug

I was trying to load a Multiscene with list of SEVIRI data in native format(reader : "seviri_l1b_native") and failed to generate the dataset error pops up everytime. (I have edited the original issue since it was actually a wrong usage of Scene with list of files as pointed out by @djhoese in the comment . Error, Code and everything is now about Multiscene.)

To Reproduce
Creating a Multiscene object with a list of files in native format that downloaded from EUMETSAT data store and then trying to load data from different channels. I can read data from each file seperately.

from satpy import MultiScene,Scene
from pyresample.geometry import AreaDefinition
from satpy.multiscene import timeseries
import glob
from satpy.utils import debug_on
debug_on()
pdir ='/DataShared/DataFiles/TB data/Met7_SEVIRI_20180812_20180819/native files/MSG1-SEVI-MSG15-0100-NA-20180812*.nat'
fs = glob.glob(pdir)
msc = MultiScene.from_files(fs,reader='seviri_l1b_native')
msc.load(['IR_108'])

new_msc = msc.resample('eurol')
blend_msc = new_msc.blend(blend_function=timeseries)

Expected behavior
Load the data from all the input file nicely in ascending order of time.

Actual results

[DEBUG: 2021-05-03 02:01:21 : satpy.readers.yaml_reader] Reading ('/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/etc/readers/seviri_l1b_native.yaml',)
[DEBUG: 2021-05-03 02:01:21 : satpy.multiscene] Forcing iteration of generator-like object of Scenes
[DEBUG: 2021-05-03 02:01:21 : satpy.readers.yaml_reader] Reading ('/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/etc/readers/seviri_l1b_native.yaml',)
[DEBUG: 2021-05-03 02:01:21 : satpy.readers.yaml_reader] Assigning to seviri_l1b_native: ['/DataShared/DataFiles/TB data/Met7_SEVIRI_20180812_20180819/native files/MSG1-SEVI-MSG15-0100-NA-20180812002741.561000000Z-NA.nat']
[DEBUG: 2021-05-03 02:01:21 : satpy.composites.config_loader] Looking for composites config file seviri.yaml
[DEBUG: 2021-05-03 02:01:21 : satpy.composites.config_loader] Looking for composites config file visir.yaml
[DEBUG: 2021-05-03 02:01:21 : native_msg] Calibration time 0:00:00.012176
[DEBUG: 2021-05-03 02:01:23 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:23 : satpy.readers.yaml_reader] Requested orientation for Dataset None is 'native' (default). No flipping is applied.
[DEBUG: 2021-05-03 02:01:23 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:23 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:23 : satpy.scene] Resampling DataID(name='IR_108', wavelength=WavelengthRange(min=9.8, central=10.8, max=11.8, unit='µm'), resolution=3000.403165817, calibration=<calibration.brightness_temperature>, modifiers=())
[DEBUG: 2021-05-03 02:01:23 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:23 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:23 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:23 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:23 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/geometry.py:1918: DeprecationWarning: 'get_xy_from_lonlat' is deprecated, please use 'get_array_indices_from_lonlat' instead.
  warnings.warn("'get_xy_from_lonlat' is deprecated, please use "
[DEBUG: 2021-05-03 02:01:25 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:25 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[INFO: 2021-05-03 02:01:25 : satpy.resample] Using default KDTree resampler
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/numpy/lib/function_base.py:1281: RuntimeWarning: invalid value encountered in subtract
  a = op(a[slice1], a[slice2])
[DEBUG: 2021-05-03 02:01:25 : satpy.resample] Check if ./resample_lut-e9afa54463fbbae510911c25098512db47a93609.npz exists
[DEBUG: 2021-05-03 02:01:25 : satpy.resample] Computing kd-tree parameters
[DEBUG: 2021-05-03 02:01:25 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:25 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1025: DeprecationWarning: `np.float` is a deprecated alias for the builtin `float`. To silence this warning, use `float` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.float64` here.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  input_coords = input_coords.astype(np.float)
[DEBUG: 2021-05-03 02:01:25 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:25 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1047: DeprecationWarning: `np.int` is a deprecated alias for the builtin `int`. To silence this warning, use `int` by itself. Doing this will not modify any behavior and is safe. When replacing `np.int`, you may wish to use e.g. `np.int64` or `np.int32` to specify the precision. If you wish to review your current use, check the release note link for additional information.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  radius=self.radius_of_influence, dtype=np.int,
[DEBUG: 2021-05-03 02:01:25 : satpy.resample] Resampling reshape-4a774fa2bc0bd6de2fedca401a9a84d9
[DEBUG: 2021-05-03 02:01:25 : satpy.readers.yaml_reader] Reading ('/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/etc/readers/seviri_l1b_native.yaml',)
[DEBUG: 2021-05-03 02:01:25 : satpy.readers.yaml_reader] Assigning to seviri_l1b_native: ['/DataShared/DataFiles/TB data/Met7_SEVIRI_20180812_20180819/native files/MSG1-SEVI-MSG15-0100-NA-20180812004241.640000000Z-NA.nat']
[DEBUG: 2021-05-03 02:01:26 : satpy.composites.config_loader] Looking for composites config file seviri.yaml
[DEBUG: 2021-05-03 02:01:26 : satpy.composites.config_loader] Looking for composites config file visir.yaml
[DEBUG: 2021-05-03 02:01:26 : native_msg] Calibration time 0:00:00.010361
[DEBUG: 2021-05-03 02:01:28 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:28 : satpy.readers.yaml_reader] Requested orientation for Dataset None is 'native' (default). No flipping is applied.
[DEBUG: 2021-05-03 02:01:28 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:28 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:28 : satpy.scene] Resampling DataID(name='IR_108', wavelength=WavelengthRange(min=9.8, central=10.8, max=11.8, unit='µm'), resolution=3000.403165817, calibration=<calibration.brightness_temperature>, modifiers=())
[DEBUG: 2021-05-03 02:01:28 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:28 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:28 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:28 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:28 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/geometry.py:1918: DeprecationWarning: 'get_xy_from_lonlat' is deprecated, please use 'get_array_indices_from_lonlat' instead.
  warnings.warn("'get_xy_from_lonlat' is deprecated, please use "
[DEBUG: 2021-05-03 02:01:30 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:30 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[INFO: 2021-05-03 02:01:30 : satpy.resample] Using default KDTree resampler
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
[DEBUG: 2021-05-03 02:01:30 : satpy.resample] Check if ./resample_lut-e9afa54463fbbae510911c25098512db47a93609.npz exists
[DEBUG: 2021-05-03 02:01:30 : satpy.resample] Computing kd-tree parameters
[DEBUG: 2021-05-03 02:01:30 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:30 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1025: DeprecationWarning: `np.float` is a deprecated alias for the builtin `float`. To silence this warning, use `float` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.float64` here.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  input_coords = input_coords.astype(np.float)
[DEBUG: 2021-05-03 02:01:30 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:30 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1047: DeprecationWarning: `np.int` is a deprecated alias for the builtin `int`. To silence this warning, use `int` by itself. Doing this will not modify any behavior and is safe. When replacing `np.int`, you may wish to use e.g. `np.int64` or `np.int32` to specify the precision. If you wish to review your current use, check the release note link for additional information.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  radius=self.radius_of_influence, dtype=np.int,
[DEBUG: 2021-05-03 02:01:30 : satpy.resample] Resampling reshape-486fca7282a054331c92e7a9c41a6e04
[DEBUG: 2021-05-03 02:01:30 : satpy.readers.yaml_reader] Reading ('/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/etc/readers/seviri_l1b_native.yaml',)
[DEBUG: 2021-05-03 02:01:30 : satpy.readers.yaml_reader] Assigning to seviri_l1b_native: ['/DataShared/DataFiles/TB data/Met7_SEVIRI_20180812_20180819/native files/MSG1-SEVI-MSG15-0100-NA-20180812005740.517000000Z-NA.nat']
[DEBUG: 2021-05-03 02:01:30 : satpy.composites.config_loader] Looking for composites config file seviri.yaml
[DEBUG: 2021-05-03 02:01:30 : satpy.composites.config_loader] Looking for composites config file visir.yaml
[DEBUG: 2021-05-03 02:01:30 : native_msg] Calibration time 0:00:00.009589
[DEBUG: 2021-05-03 02:01:32 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:32 : satpy.readers.yaml_reader] Requested orientation for Dataset None is 'native' (default). No flipping is applied.
[DEBUG: 2021-05-03 02:01:33 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:33 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:33 : satpy.scene] Resampling DataID(name='IR_108', wavelength=WavelengthRange(min=9.8, central=10.8, max=11.8, unit='µm'), resolution=3000.403165817, calibration=<calibration.brightness_temperature>, modifiers=())
[DEBUG: 2021-05-03 02:01:33 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:33 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:33 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:33 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:33 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/geometry.py:1918: DeprecationWarning: 'get_xy_from_lonlat' is deprecated, please use 'get_array_indices_from_lonlat' instead.
  warnings.warn("'get_xy_from_lonlat' is deprecated, please use "
[DEBUG: 2021-05-03 02:01:35 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:35 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[INFO: 2021-05-03 02:01:35 : satpy.resample] Using default KDTree resampler
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
[DEBUG: 2021-05-03 02:01:35 : satpy.resample] Check if ./resample_lut-e9afa54463fbbae510911c25098512db47a93609.npz exists
[DEBUG: 2021-05-03 02:01:35 : satpy.resample] Computing kd-tree parameters
[DEBUG: 2021-05-03 02:01:35 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:35 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1025: DeprecationWarning: `np.float` is a deprecated alias for the builtin `float`. To silence this warning, use `float` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.float64` here.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  input_coords = input_coords.astype(np.float)
[DEBUG: 2021-05-03 02:01:35 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:35 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1047: DeprecationWarning: `np.int` is a deprecated alias for the builtin `int`. To silence this warning, use `int` by itself. Doing this will not modify any behavior and is safe. When replacing `np.int`, you may wish to use e.g. `np.int64` or `np.int32` to specify the precision. If you wish to review your current use, check the release note link for additional information.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  radius=self.radius_of_influence, dtype=np.int,
[DEBUG: 2021-05-03 02:01:35 : satpy.resample] Resampling reshape-364c754d3d1d5d117db8942f52064c3b
[DEBUG: 2021-05-03 02:01:35 : satpy.readers.yaml_reader] Reading ('/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/etc/readers/seviri_l1b_native.yaml',)
[DEBUG: 2021-05-03 02:01:35 : satpy.readers.yaml_reader] Assigning to seviri_l1b_native: ['/DataShared/DataFiles/TB data/Met7_SEVIRI_20180812_20180819/native files/MSG1-SEVI-MSG15-0100-NA-20180812011240.595000000Z-NA.nat']
[DEBUG: 2021-05-03 02:01:35 : satpy.composites.config_loader] Looking for composites config file seviri.yaml
[DEBUG: 2021-05-03 02:01:35 : satpy.composites.config_loader] Looking for composites config file visir.yaml
[DEBUG: 2021-05-03 02:01:35 : native_msg] Calibration time 0:00:00.010877
[DEBUG: 2021-05-03 02:01:37 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:37 : satpy.readers.yaml_reader] Requested orientation for Dataset None is 'native' (default). No flipping is applied.
[DEBUG: 2021-05-03 02:01:37 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:37 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:37 : satpy.scene] Resampling DataID(name='IR_108', wavelength=WavelengthRange(min=9.8, central=10.8, max=11.8, unit='µm'), resolution=3000.403165817, calibration=<calibration.brightness_temperature>, modifiers=())
[DEBUG: 2021-05-03 02:01:37 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:37 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:37 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:37 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:37 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/geometry.py:1918: DeprecationWarning: 'get_xy_from_lonlat' is deprecated, please use 'get_array_indices_from_lonlat' instead.
  warnings.warn("'get_xy_from_lonlat' is deprecated, please use "
[DEBUG: 2021-05-03 02:01:39 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:39 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[INFO: 2021-05-03 02:01:39 : satpy.resample] Using default KDTree resampler
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
[DEBUG: 2021-05-03 02:01:39 : satpy.resample] Check if ./resample_lut-e9afa54463fbbae510911c25098512db47a93609.npz exists
[DEBUG: 2021-05-03 02:01:39 : satpy.resample] Computing kd-tree parameters
[DEBUG: 2021-05-03 02:01:39 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:39 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1025: DeprecationWarning: `np.float` is a deprecated alias for the builtin `float`. To silence this warning, use `float` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.float64` here.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  input_coords = input_coords.astype(np.float)
[DEBUG: 2021-05-03 02:01:39 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:39 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1047: DeprecationWarning: `np.int` is a deprecated alias for the builtin `int`. To silence this warning, use `int` by itself. Doing this will not modify any behavior and is safe. When replacing `np.int`, you may wish to use e.g. `np.int64` or `np.int32` to specify the precision. If you wish to review your current use, check the release note link for additional information.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  radius=self.radius_of_influence, dtype=np.int,
[DEBUG: 2021-05-03 02:01:39 : satpy.resample] Resampling reshape-f69e0f8c2d2b5f5cc816708afb58d3f0
[DEBUG: 2021-05-03 02:01:39 : satpy.readers.yaml_reader] Reading ('/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/etc/readers/seviri_l1b_native.yaml',)
[DEBUG: 2021-05-03 02:01:39 : satpy.readers.yaml_reader] Assigning to seviri_l1b_native: ['/DataShared/DataFiles/TB data/Met7_SEVIRI_20180812_20180819/native files/MSG1-SEVI-MSG15-0100-NA-20180812012740.674000000Z-NA.nat']
[DEBUG: 2021-05-03 02:01:39 : satpy.composites.config_loader] Looking for composites config file seviri.yaml
[DEBUG: 2021-05-03 02:01:40 : satpy.composites.config_loader] Looking for composites config file visir.yaml
[DEBUG: 2021-05-03 02:01:40 : native_msg] Calibration time 0:00:00.010977
[DEBUG: 2021-05-03 02:01:42 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:42 : satpy.readers.yaml_reader] Requested orientation for Dataset None is 'native' (default). No flipping is applied.
[DEBUG: 2021-05-03 02:01:42 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:42 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:42 : satpy.scene] Resampling DataID(name='IR_108', wavelength=WavelengthRange(min=9.8, central=10.8, max=11.8, unit='µm'), resolution=3000.403165817, calibration=<calibration.brightness_temperature>, modifiers=())
[DEBUG: 2021-05-03 02:01:42 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:42 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:42 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:42 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:42 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/geometry.py:1918: DeprecationWarning: 'get_xy_from_lonlat' is deprecated, please use 'get_array_indices_from_lonlat' instead.
  warnings.warn("'get_xy_from_lonlat' is deprecated, please use "
[DEBUG: 2021-05-03 02:01:44 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:44 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[INFO: 2021-05-03 02:01:44 : satpy.resample] Using default KDTree resampler
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
[DEBUG: 2021-05-03 02:01:44 : satpy.resample] Check if ./resample_lut-e9afa54463fbbae510911c25098512db47a93609.npz exists
[DEBUG: 2021-05-03 02:01:44 : satpy.resample] Computing kd-tree parameters
[DEBUG: 2021-05-03 02:01:44 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:44 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1025: DeprecationWarning: `np.float` is a deprecated alias for the builtin `float`. To silence this warning, use `float` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.float64` here.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  input_coords = input_coords.astype(np.float)
[DEBUG: 2021-05-03 02:01:44 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:44 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1047: DeprecationWarning: `np.int` is a deprecated alias for the builtin `int`. To silence this warning, use `int` by itself. Doing this will not modify any behavior and is safe. When replacing `np.int`, you may wish to use e.g. `np.int64` or `np.int32` to specify the precision. If you wish to review your current use, check the release note link for additional information.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  radius=self.radius_of_influence, dtype=np.int,
[DEBUG: 2021-05-03 02:01:44 : satpy.resample] Resampling reshape-c1507ddb7c2422c6b641288fce24bfbb
[DEBUG: 2021-05-03 02:01:44 : satpy.readers.yaml_reader] Reading ('/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/etc/readers/seviri_l1b_native.yaml',)
[DEBUG: 2021-05-03 02:01:44 : satpy.readers.yaml_reader] Assigning to seviri_l1b_native: ['/DataShared/DataFiles/TB data/Met7_SEVIRI_20180812_20180819/native files/MSG1-SEVI-MSG15-0100-NA-20180812014240.752000000Z-NA.nat']
[DEBUG: 2021-05-03 02:01:44 : satpy.composites.config_loader] Looking for composites config file seviri.yaml
[DEBUG: 2021-05-03 02:01:44 : satpy.composites.config_loader] Looking for composites config file visir.yaml
[DEBUG: 2021-05-03 02:01:44 : native_msg] Calibration time 0:00:00.009976
[DEBUG: 2021-05-03 02:01:46 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:46 : satpy.readers.yaml_reader] Requested orientation for Dataset None is 'native' (default). No flipping is applied.
[DEBUG: 2021-05-03 02:01:46 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:46 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:47 : satpy.scene] Resampling DataID(name='IR_108', wavelength=WavelengthRange(min=9.8, central=10.8, max=11.8, unit='µm'), resolution=3000.403165817, calibration=<calibration.brightness_temperature>, modifiers=())
[DEBUG: 2021-05-03 02:01:47 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:47 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:47 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:47 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:47 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/geometry.py:1918: DeprecationWarning: 'get_xy_from_lonlat' is deprecated, please use 'get_array_indices_from_lonlat' instead.
  warnings.warn("'get_xy_from_lonlat' is deprecated, please use "
[DEBUG: 2021-05-03 02:01:48 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:48 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[INFO: 2021-05-03 02:01:48 : satpy.resample] Using default KDTree resampler
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
[DEBUG: 2021-05-03 02:01:48 : satpy.resample] Check if ./resample_lut-e9afa54463fbbae510911c25098512db47a93609.npz exists
[DEBUG: 2021-05-03 02:01:48 : satpy.resample] Computing kd-tree parameters
[DEBUG: 2021-05-03 02:01:48 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:48 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1025: DeprecationWarning: `np.float` is a deprecated alias for the builtin `float`. To silence this warning, use `float` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.float64` here.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  input_coords = input_coords.astype(np.float)
[DEBUG: 2021-05-03 02:01:48 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:48 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1047: DeprecationWarning: `np.int` is a deprecated alias for the builtin `int`. To silence this warning, use `int` by itself. Doing this will not modify any behavior and is safe. When replacing `np.int`, you may wish to use e.g. `np.int64` or `np.int32` to specify the precision. If you wish to review your current use, check the release note link for additional information.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  radius=self.radius_of_influence, dtype=np.int,
[DEBUG: 2021-05-03 02:01:48 : satpy.resample] Resampling reshape-ce4cb1af385eb4f3a2656a40436fab79
[DEBUG: 2021-05-03 02:01:48 : satpy.readers.yaml_reader] Reading ('/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/etc/readers/seviri_l1b_native.yaml',)
[DEBUG: 2021-05-03 02:01:48 : satpy.readers.yaml_reader] Assigning to seviri_l1b_native: ['/DataShared/DataFiles/TB data/Met7_SEVIRI_20180812_20180819/native files/MSG1-SEVI-MSG15-0100-NA-20180812015740.831000000Z-NA.nat']
[DEBUG: 2021-05-03 02:01:49 : satpy.composites.config_loader] Looking for composites config file seviri.yaml
[DEBUG: 2021-05-03 02:01:49 : satpy.composites.config_loader] Looking for composites config file visir.yaml
[DEBUG: 2021-05-03 02:01:49 : native_msg] Calibration time 0:00:00.009599
[DEBUG: 2021-05-03 02:01:51 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:51 : satpy.readers.yaml_reader] Requested orientation for Dataset None is 'native' (default). No flipping is applied.
[DEBUG: 2021-05-03 02:01:51 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:51 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:51 : satpy.scene] Resampling DataID(name='IR_108', wavelength=WavelengthRange(min=9.8, central=10.8, max=11.8, unit='µm'), resolution=3000.403165817, calibration=<calibration.brightness_temperature>, modifiers=())
[DEBUG: 2021-05-03 02:01:51 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:51 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:51 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:51 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:51 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/geometry.py:1918: DeprecationWarning: 'get_xy_from_lonlat' is deprecated, please use 'get_array_indices_from_lonlat' instead.
  warnings.warn("'get_xy_from_lonlat' is deprecated, please use "
[DEBUG: 2021-05-03 02:01:53 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:53 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[INFO: 2021-05-03 02:01:53 : satpy.resample] Using default KDTree resampler
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
[DEBUG: 2021-05-03 02:01:53 : satpy.resample] Check if ./resample_lut-e9afa54463fbbae510911c25098512db47a93609.npz exists
[DEBUG: 2021-05-03 02:01:53 : satpy.resample] Computing kd-tree parameters
[DEBUG: 2021-05-03 02:01:53 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:53 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1025: DeprecationWarning: `np.float` is a deprecated alias for the builtin `float`. To silence this warning, use `float` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.float64` here.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  input_coords = input_coords.astype(np.float)
[DEBUG: 2021-05-03 02:01:53 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:53 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1047: DeprecationWarning: `np.int` is a deprecated alias for the builtin `int`. To silence this warning, use `int` by itself. Doing this will not modify any behavior and is safe. When replacing `np.int`, you may wish to use e.g. `np.int64` or `np.int32` to specify the precision. If you wish to review your current use, check the release note link for additional information.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  radius=self.radius_of_influence, dtype=np.int,
[DEBUG: 2021-05-03 02:01:53 : satpy.resample] Resampling reshape-a34556adee3f74e3e67b3a7fcd5190af
[DEBUG: 2021-05-03 02:01:53 : satpy.readers.yaml_reader] Reading ('/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/etc/readers/seviri_l1b_native.yaml',)
[DEBUG: 2021-05-03 02:01:53 : satpy.readers.yaml_reader] Assigning to seviri_l1b_native: ['/DataShared/DataFiles/TB data/Met7_SEVIRI_20180812_20180819/native files/MSG1-SEVI-MSG15-0100-NA-20180812021240.910000000Z-NA.nat']
[DEBUG: 2021-05-03 02:01:53 : satpy.composites.config_loader] Looking for composites config file seviri.yaml
[DEBUG: 2021-05-03 02:01:53 : satpy.composites.config_loader] Looking for composites config file visir.yaml
[DEBUG: 2021-05-03 02:01:53 : native_msg] Calibration time 0:00:00.010009
[DEBUG: 2021-05-03 02:01:55 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:55 : satpy.readers.yaml_reader] Requested orientation for Dataset None is 'native' (default). No flipping is applied.
[DEBUG: 2021-05-03 02:01:56 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:56 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:56 : satpy.scene] Resampling DataID(name='IR_108', wavelength=WavelengthRange(min=9.8, central=10.8, max=11.8, unit='µm'), resolution=3000.403165817, calibration=<calibration.brightness_temperature>, modifiers=())
[DEBUG: 2021-05-03 02:01:56 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:56 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:56 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:56 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:56 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/geometry.py:1918: DeprecationWarning: 'get_xy_from_lonlat' is deprecated, please use 'get_array_indices_from_lonlat' instead.
  warnings.warn("'get_xy_from_lonlat' is deprecated, please use "
[DEBUG: 2021-05-03 02:01:58 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:58 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[INFO: 2021-05-03 02:01:58 : satpy.resample] Using default KDTree resampler
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
[DEBUG: 2021-05-03 02:01:58 : satpy.resample] Check if ./resample_lut-e9afa54463fbbae510911c25098512db47a93609.npz exists
[DEBUG: 2021-05-03 02:01:58 : satpy.resample] Computing kd-tree parameters
[DEBUG: 2021-05-03 02:01:58 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:58 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1025: DeprecationWarning: `np.float` is a deprecated alias for the builtin `float`. To silence this warning, use `float` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.float64` here.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  input_coords = input_coords.astype(np.float)
[DEBUG: 2021-05-03 02:01:58 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
[DEBUG: 2021-05-03 02:01:58 : pyproj] PROJ_ERROR: proj_crs_get_sub_crs: Object is not a CompoundCRS
/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/pyresample/kd_tree.py:1047: DeprecationWarning: `np.int` is a deprecated alias for the builtin `int`. To silence this warning, use `int` by itself. Doing this will not modify any behavior and is safe. When replacing `np.int`, you may wish to use e.g. `np.int64` or `np.int32` to specify the precision. If you wish to review your current use, check the release note link for additional information.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
  radius=self.radius_of_influence, dtype=np.int,
[DEBUG: 2021-05-03 02:01:58 : satpy.resample] Resampling reshape-5238f226a0a2058555d988dec5d0ebd7
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
Traceback (most recent call last):
  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/dataset/metadata.py", line 179, in _all_values_equal
    return _pairwise_all(nan_allclose, values)
  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/dataset/metadata.py", line 144, in _pairwise_all
    if not func(values[0], value):
  File "<__array_function__ internals>", line 5, in allclose
  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/numpy/core/numeric.py", line 2256, in allclose
    res = all(isclose(a, b, rtol=rtol, atol=atol, equal_nan=equal_nan))
  File "<__array_function__ internals>", line 5, in isclose
  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/numpy/core/numeric.py", line 2362, in isclose
    xfin = isfinite(x)
TypeError: ufunc 'isfinite' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule ''safe''

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/home/edith005/Project2020/ProjectCode/Metsat_SEVIRI/error_multiscene.py", line 13, in <module>
    blend_msc = new_msc.blend(blend_function=timeseries)
  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/multiscene.py", line 306, in blend
    new_scn[ds_id] = blend_function(datasets)
  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/multiscene.py", line 67, in timeseries
    res.attrs = combine_metadata(*[x.attrs for x in expanded_ds])
  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/dataset/metadata.py", line 53, in combine_metadata
    return _combine_shared_info(shared_keys, info_dicts, average_times)
  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/dataset/metadata.py", line 81, in _combine_shared_info
    elif _are_values_combinable(values):
  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/dataset/metadata.py", line 116, in _are_values_combinable
    return _all_values_equal(values)
  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/dataset/metadata.py", line 181, in _all_values_equal
    return _pairwise_all(eq, values)
  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/dataset/metadata.py", line 144, in _pairwise_all
    if not func(values[0], value):
**ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()**

Screenshots
If applicable, add screenshots to help explain your problem.

Environment Info:

  • OS: Linux Mint 20 x86_64
  • Satpy Version: 0.27.1.dev60+gcec1f0e1
  • PyResample Version: 1.19.0
  • Readers and writers dependencies (when relevant): [run from satpy.utils import check_satpy; check_satpy()]

Additional context
Add any other context about the problem here.
@djhoese
Copy link
Member

djhoese commented Apr 30, 2021

I think this was fixed in #1643. Is there any chance you could try installing satpy from the master branch on github and try running again?

@djhoese
Copy link
Member

djhoese commented Apr 30, 2021

oh wait, are you providing more than one "time step" to the Scene? Try only providing one (ex. one full disk set of files for a specific time).

@harry9713
Copy link
Author

oh wait, are you providing more than one "time step" to the Scene? Try only providing one (ex. one full disk set of files for a specific time).

Yes that works perfectly and I've mentioned it in my issue. I can read each files seperately. The issue is when i try to create a scene object witha a bunch of files simultaneously and load a channel.
So I guess there is no way to load that as of now, right? :)

@djhoese
Copy link
Member

djhoese commented Apr 30, 2021

@harry9713 Sorry for missing that. You can use the MultiScene object. The Scene object is not meant to handle more than one time step (although this error is particularly unhelpful) and will likely never be changed to do that.

https://satpy.readthedocs.io/en/stable/multiscene.html

@harry9713
Copy link
Author

@harry9713 Sorry for missing that. You can use the MultiScene object. The Scene object is not meant to handle more than one time step (although this error is particularly unhelpful) and will likely never be changed to do that.

https://satpy.readthedocs.io/en/stable/multiscene.html

Thank you for your quick response. I did try Multiscene module to load the data but encountered the below error. Infact, I used the given example of Multiscene to blend, but still the error pops up. I think there is something I am not understanding here.

from satpy.multiscene import timeseries

msc.load(['IR_108'])
new_msc = msc.resample('euro4',resampler='nearest')
blend_msc = new_msc.blend(blend_function=timeseries)

gives the error

  a = op(a[slice1], a[slice2])
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
/home/edith005/anaconda3/envs/natsat/lib/python3.9/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
  proj_string = self.to_proj4()
Traceback (most recent call last):

  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/dataset/metadata.py", line 179, in _all_values_equal
    return _pairwise_all(nan_allclose, values)

  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/dataset/metadata.py", line 144, in _pairwise_all
    if not func(values[0], value):

  File "<__array_function__ internals>", line 5, in allclose

  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/numpy/core/numeric.py", line 2256, in allclose
    res = all(isclose(a, b, rtol=rtol, atol=atol, equal_nan=equal_nan))

  File "<__array_function__ internals>", line 5, in isclose

  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/numpy/core/numeric.py", line 2362, in isclose
    xfin = isfinite(x)

TypeError: ufunc 'isfinite' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule ''safe''


During handling of the above exception, another exception occurred:

Traceback (most recent call last):

  File "/home/edith005/Project2020/ProjectCode/Metsat_SEVIRI/multiscene_try.py", line 24, in <module>
    blend_msc = new_msc.blend(blend_function=timeseries)

  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/multiscene.py", line 306, in blend
    new_scn[ds_id] = blend_function(datasets)

  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/multiscene.py", line 67, in timeseries
    res.attrs = combine_metadata(*[x.attrs for x in expanded_ds])

  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/dataset/metadata.py", line 53, in combine_metadata
    return _combine_shared_info(shared_keys, info_dicts, average_times)

  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/dataset/metadata.py", line 81, in _combine_shared_info
    elif _are_values_combinable(values):

  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/dataset/metadata.py", line 116, in _are_values_combinable
    return _all_values_equal(values)

  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/dataset/metadata.py", line 181, in _all_values_equal
    return _pairwise_all(eq, values)

  File "/home/edith005/anaconda3/envs/natsat/lib/python3.9/site-packages/satpy/dataset/metadata.py", line 144, in _pairwise_all
    if not func(values[0], value):

ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()

@djhoese
Copy link
Member

djhoese commented May 1, 2021

How are you creating the MultiScene? It may just be that although you are doing things correctly this is still a bug that had to be fixed (and is fixed in the master branch...hopefully).

@harry9713
Copy link
Author

harry9713 commented May 1, 2021 •
edited

How are you creating the MultiScene? It may just be that although you are doing things correctly this is still a bug that had to be fixed (and is fixed in the master branch...hopefully).

I did actually switch to master as per your instruction and using it from then. Still the issue persists. In the example code given , the multiscene object from a list of files is created using mscn = MultiScene.from_files(glob('/data/abi/day_1/*C0[12]*.nc'), reader='abi_l1b') , am doing the same by changing the reader.

pdir ='/DataShared/DataFiles/TB data/Met7_SEVIRI_20180812_20180819/native files/MSG1-SEVI-MSG15-0100-NA-20180812*.nat'
fs = glob.glob(pdir)
msc = MultiScene.from_files(fs,reader='seviri_l1b_native')
msc.load(['IR_108'])

Then I used the below example code for blending the multiscene object as timeseries using

>>> mscn.load(['I04'])
>>> new_mscn = mscn.resample(my_area)
>>> blended_scene = new_mscn.blend(blend_function=timeseries)

the last line gives me the errors I posted in the previous comment. I can still work with single dataset using
mscn.Scenes[2]['IR_108'] which gives me the dataset in the third timestep without any problem.

@simonrp84
Copy link
Member

mscn.load(['I04']) is not a SEVIRI channel. Change it to: mscn.load(['IR_108'])

@harry9713
Copy link
Author

mscn.load(['I04']) is not a SEVIRI channel. Change it to: mscn.load(['IR_108'])

I was explaining about the example code I've used. I've already made required changes with channels and reader. :)

@djhoese
Copy link
Member

djhoese commented May 2, 2021

I don't have that much SEVIRI data to test this. Maybe one of the euro-based satpy devs can look at this a little more. @mraspaud @pnuu @gerritholl ?

I'm a little confused at what has changed that is causing these new metadata related issues lately.

@simonrp84
Copy link
Member

Just to note, I get the same error as @harry9713 when running similar code for SEVIRI native. Something must've changed in the metadata.

@sfinkens sfinkens mentioned this issue May 6, 2021
Merged
2 tasks
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

Successfully merging a pull request may close this issue.

Ignore raw metadata when combining metadata sfinkens/satpy
3 participants
@djhoese @simonrp84 @harry9713