@@ -43,52 +43,22 @@ def __call__(self, transform_xy, x1, y1, x2, y2):
         if hasattr(super(), '_find_transformed_bbox'):
             return super().__call__(transform_xy, x1, y1, x2, y2)
 
-        x, y = np.meshgrid(
-            np.linspace(x1, x2, self.nx), np.linspace(y1, y2, self.ny))
-        lon, lat = transform_xy(np.ravel(x), np.ravel(y))
-
-        with np.errstate(invalid='ignore'):
-            if self.lon_cycle is not None:
-                lon0 = np.nanmin(lon)
-                # Changed from 180 to 360 to be able to span only
-                # 90-270 (left hand side)
-                lon -= 360. * ((lon - lon0) > 360.)
-            if self.lat_cycle is not None:  # pragma: no cover
-                lat0 = np.nanmin(lat)
-                lat -= 360. * ((lat - lat0) > 180.)
-
-        lon_min, lon_max = np.nanmin(lon), np.nanmax(lon)
-        lat_min, lat_max = np.nanmin(lat), np.nanmax(lat)
-
-        lon_min, lon_max, lat_min, lat_max = \
-            self._add_pad(lon_min, lon_max, lat_min, lat_max)
-
-        # check cycle
-        if self.lon_cycle:
-            lon_max = min(lon_max, lon_min + self.lon_cycle)
-        if self.lat_cycle:  # pragma: no cover
-            lat_max = min(lat_max, lat_min + self.lat_cycle)
-
-        if self.lon_minmax is not None:
-            min0 = self.lon_minmax[0]
-            lon_min = max(min0, lon_min)
-            max0 = self.lon_minmax[1]
-            lon_max = min(max0, lon_max)
-
-        if self.lat_minmax is not None:
-            min0 = self.lat_minmax[0]
-            lat_min = max(min0, lat_min)
-            max0 = self.lat_minmax[1]
-            lat_max = min(max0, lat_max)
-
+        transformed_bbox = self._find_extremes(
+            lambda points: np.column_stack(transform_xy(*points.T)),
+            Bbox.from_extents(x1, y1, x2, y2))
+        lon_min, lat_min, lon_max, lat_max = transformed_bbox.extents
         return lon_min, lon_max, lat_min, lat_max
 
     def _find_transformed_bbox(self, trans, bbox):
         """Find grid extremes using the Matplotlib 3.11 interface."""
+        return self._find_extremes(trans.transform, bbox)
+
+    def _find_extremes(self, transform, bbox):
+        """Find grid extremes using a point-array transform."""
         grid = np.reshape(np.meshgrid(
             np.linspace(bbox.x0, bbox.x1, self.nx),
             np.linspace(bbox.y0, bbox.y1, self.ny)), (2, -1)).T
-        lon, lat = trans.transform(grid).T
+        lon, lat = transform(grid).T
 
         with np.errstate(invalid='ignore'):
             if self.lon_cycle is not None: