diff --git a/xpublish_wms/grids/regular.py b/xpublish_wms/grids/regular.py
index b594137..6e7fff9 100644
--- a/xpublish_wms/grids/regular.py
+++ b/xpublish_wms/grids/regular.py
@@ -20,7 +20,7 @@ def name(self) -> str:
 
     @property
     def render_method(self) -> RenderMethod:
-        return RenderMethod.Raster
+        return RenderMethod.Quad
 
     @property
     def crs(self) -> str:
diff --git a/xpublish_wms/wms/get_map.py b/xpublish_wms/wms/get_map.py
index 2b089d7..162ba0c 100644
--- a/xpublish_wms/wms/get_map.py
+++ b/xpublish_wms/wms/get_map.py
@@ -352,11 +352,17 @@ def render(
                 da,
             )
         elif ds.gridded.render_method == RenderMethod.Quad:
-            mesh = cvs.quadmesh(
-                da,
-                x="x",
-                y="y",
-            )
+            try:
+                mesh = cvs.quadmesh(
+                    da,
+                    x="x",
+                    y="y",
+                )
+            except Exception as e:
+                logger.warning(f"Error rendering quadmesh: {e}, falling back to raster")
+                mesh = cvs.raster(
+                    da,
+                )
         elif ds.gridded.render_method == RenderMethod.Triangle:
             triangles = ds.gridded.tessellate(da)
             if x is not None: