FIX: add support for imshow extent to have units

pratimugale · pratimugale · commit b0e758df7741 · 2022-10-15T22:14:39.000+05:30
diff --git a/doc/users/next_whats_new/imshow_extent_units.rst b/doc/users/next_whats_new/imshow_extent_units.rst
@@ -0,0 +1,26 @@
+The *extent* of ``imshow`` can now be expressed with units
+----------------------------------------------------------
+The *extent* parameter of `~.axes.Axes.imshow` and  `~.AxesImage.set_extent`
+can now be expressed with units.
+
+.. plot::
+    :include-source: true
+
+    import matplotlib.pyplot as plt
+    import numpy as np
+    from matplotlib.dates import HourLocator, DateFormatter
+    from matplotlib.ticker import AutoMinorLocator
+
+    fig, ax = plt.subplots()
+    dateFirst = np.datetime64('2020-01-01', 'D')
+    dateLast = np.datetime64('2020-01-11', 'D')
+
+    arr = [[i+j for i in range(10)] for j in range(10)]
+
+    ax.imshow(arr, origin='lower', extent=[1, 11, dateFirst, dateLast])
+
+    ax.yaxis.set_major_formatter(DateFormatter('%d/%m/%y:- %H00hours'))
+    ax.yaxis.set_major_locator(HourLocator(byhour=[0, 6, 12, 18, 24]))
+    ax.yaxis.set_minor_locator(AutoMinorLocator())
+
+    plt.show()
diff --git a/lib/matplotlib/axes/_axes.py b/lib/matplotlib/axes/_axes.py
@@ -5523,6 +5523,7 @@ def imshow(self, X, cmap=None, norm=None, aspect=None,
 
         extent : floats (left, right, bottom, top), optional
             The bounding box in data coordinates that the image will fill.
+            These values may be unitful and match the units of the Axes.
             The image is stretched individually along x and y to fill the box.
 
             The default extent is determined by the following conditions.
@@ -5615,7 +5616,7 @@ def imshow(self, X, cmap=None, norm=None, aspect=None,
 
         # update ax.dataLim, and, if autoscaling, set viewLim
         # to tightly fit the image, regardless of dataLim.
-        im.set_extent(im.get_extent())
+        im.set_extent(im.get_extent(), **kwargs)
 
         self.add_image(im)
         return im
diff --git a/lib/matplotlib/image.py b/lib/matplotlib/image.py
@@ -953,7 +953,7 @@ def _check_unsampled_image(self):
         """Return whether the image would be better drawn unsampled."""
         return self.get_interpolation() == "none"
 
-    def set_extent(self, extent):
+    def set_extent(self, extent, **kwargs):
         """
         Set the image extent.
 
@@ -962,6 +962,10 @@ def set_extent(self, extent):
         extent : 4-tuple of float
             The position and size of the image as tuple
             ``(left, right, bottom, top)`` in data coordinates.
+        **kwargs
+            Other parameters from which unit info (i.e., the *xunits*,
+            *yunits*, *zunits* (for 3D axes), *runits* and *thetaunits* (for
+            polar axes) entries are applied, if present.
 
         Notes
         -----
@@ -970,7 +974,21 @@ def set_extent(self, extent):
         state is not changed, so following this with ``ax.autoscale_view()``
         will redo the autoscaling in accord with ``dataLim``.
         """
-        self._extent = xmin, xmax, ymin, ymax = extent
+        (xmin, xmax), (ymin, ymax) = self.axes._process_unit_info(
+            [("x", [extent[0], extent[1]]),
+             ("y", [extent[2], extent[3]])],
+            kwargs)
+        xmin = self.axes._validate_converted_limits(
+            xmin, self.convert_xunits)
+        xmax = self.axes._validate_converted_limits(
+            xmax, self.convert_xunits)
+        ymin = self.axes._validate_converted_limits(
+            ymin, self.convert_yunits)
+        ymax = self.axes._validate_converted_limits(
+            ymax, self.convert_yunits)
+        extent = [xmin, xmax, ymin, ymax]
+
+        self._extent = extent
         corners = (xmin, ymin), (xmax, ymax)
         self.axes.update_datalim(corners)
         self.sticky_edges.x[:] = [xmin, xmax]
diff --git a/lib/matplotlib/tests/baseline_images/test_axes/extent_units.pdf b/lib/matplotlib/tests/baseline_images/test_axes/extent_units.pdf
diff --git a/lib/matplotlib/tests/test_axes.py b/lib/matplotlib/tests/test_axes.py
@@ -8195,3 +8195,15 @@ def test_bar_leading_nan():
         for b in rest:
             assert np.isfinite(b.xy).all()
             assert np.isfinite(b.get_width())
+
+
+@image_comparison(["extent_units.pdf"], style="mpl20")
+def test_extent_units():
+    mpl.style.use("mpl20")
+    _, ax = plt.subplots()
+    dateFirst = np.datetime64('2020-01-01', 'D')
+    dateLast = np.datetime64('2020-01-11', 'D')
+
+    arr = [[i+j for i in range(10)] for j in range(10)]
+    ax.imshow(arr, origin='lower',
+              extent=[1, 11, dateFirst, dateLast])
