WCS Plotter (#73)

* Catch ValueError in 1D animation update

* Catch error setting clim in 2D plots

* Changelog

* Check for NaNs instead of catching ValueError

* And again

* Well obviously that only works for dask arrays

* Use nanmin/nanmax in case some but not all of the data is nan

* [pre-commit.ci] auto fixes from pre-commit.com hooks

for more information, see https://pre-commit.ci

* Not sure how this became this PR's problem, but fix the docs build

* Update tox.ini

* Update pyproject.toml

* Update tox.ini

---------

Co-authored-by: Drew Leonard <[email protected]>
Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com>
Co-authored-by: Stuart Mumford <[email protected]>

Author: 4 people

changelog/73.trivial.rst

@@ -0,0 +1 @@
+Catch errors when animating plots with missing data and replace them with warnings.

mpl_animators/wcs.py

@@ -1,3 +1,4 @@
+import warnings
 from functools import partial
 
 import numpy as np
@@ -277,8 +278,13 @@ def update_plot_1d(self, val, line, slider):
 
         # If we are not setting ylim globally then we set it per frame.
         if self.ylim == 'dynamic':
-            self.axes.set_ylim(float(self.data[self.frame_index].min()),
-                               float(self.data[self.frame_index].max()))
+            vmin = float(np.nanmin(self.data[self.frame_index]))
+            vmax = float(np.nanmax(self.data[self.frame_index]))
+            if np.isnan(vmin) or np.isnan(vmax):
+                warnings.warn(UserWarning(f"No data found for data slice {self.frame_index} - cannot set ylim"))
+                return
+
+            self.axes.set_ylim(vmin, vmax)
         slider.cval = val
 
     def plot_start_image_2d(self, ax):
@@ -315,7 +321,12 @@ def _get_2d_plot_limits(self):
         Get vmin, vmax of a data slice when clip_interval is specified.
         """
         percent_limits = self.clip_interval.to('%').value
-        vmin, vmax = AsymmetricPercentileInterval(*percent_limits).get_limits(self.data_transposed)
+        if np.isnan(self.data_transposed).all():
+            warnings.warn(UserWarning(f"No data found for data slice {self.frame_index} - cannot set vmin, vmax"))
+            vmin, vmax = 0, 0
+        else:
+            vmin, vmax = AsymmetricPercentileInterval(*percent_limits).get_limits(self.data_transposed)
+
         return vmin, vmax
 
     def update_plot_2d(self, val, im, slider):

pyproject.toml

@@ -39,7 +39,7 @@ docs = [
   "sunpy-sphinx-theme",
   "packaging",
   "sphinx_gallery",
-  "sunpy",
+  "sunpy[map]",
   "scipy",
 ]
 all = [