Back to basic

hvplot/converter.py

@@ -407,6 +407,7 @@ def __init__(
         self.dynamic = dynamic
         self.geo = any([geo, crs, global_extent, projection, project, coastline, features])
         self.crs = self._process_crs(data, crs) if self.geo else None
+        self.output_projection = self.crs
         self.project = project
         self.coastline = coastline
         self.features = features
@@ -584,6 +585,8 @@ def __init__(
 
         if self.crs and global_extent:
             plot_opts['global_extent'] = global_extent
+        if self.crs != self.output_projection:
+            plot_opts['projection'] = self.output_projection
         title = title if title is not None else getattr(self, '_title', None)
         if title is not None:
             plot_opts['title'] = title
@@ -1260,12 +1263,10 @@ def method_wrapper(ds, x, y):
                 obj = method(x, y)
                 obj._dataset = dataset
 
-        if self.geo and self.crs != self.output_projection:
+        if self.geo and self.project:
+            # Apply projection before rasterizing
             import geoviews as gv
-            if isinstance(obj, gv.element.geo._Element):
-                obj.opts(projection=self.output_projection)
-            else:
-                obj = gv.project(obj, projection=self.output_projection)
+            obj = gv.project(obj, projection=self.output_projection)
 
         if not (self.datashade or self.rasterize or self.downsample):
             layers = self._apply_layers(obj)

hvplot/tests/testgeo.py

@@ -140,6 +140,7 @@ def test_geo_with_rasterize(self):
             rasterize=True,
             dynamic=False,
             aggregator="max",
+            project=True,
         )
 
         p1 = gv.Points(ds.isel(time=0), kdims=["lon", "lat"], crs=ccrs.PlateCarree())