Fixed the OGR driver to use point or line spatial filter geometries i…

…n degenerated cases (#4420)

HISTORY.TXT

@@ -16,6 +16,8 @@ Current Version (git master, 6.3-dev, future 6.4):
 
 - Add notable changes here. Less important changes should go only in the
   commit message and not appear here.
+
+- Fixed the OGR driver to use point or line spatial filter geometries in degenerated cases (#4420)
 
 Version 6.2 (beta1: 20120629): 
 -------------------------------------------------

mapogr.cpp

@@ -1225,20 +1225,44 @@ static int msOGRFileWhichShapes(layerObj *layer, rectObj rect,
    * ------------------------------------------------------------------ */
   ACQUIRE_OGR_LOCK;
 
-  OGRGeometryH hSpatialFilterPolygon = OGR_G_CreateGeometry( wkbPolygon );
-  OGRGeometryH hRing = OGR_G_CreateGeometry( wkbLinearRing );
+  if (rect.minx == rect.maxx && rect.miny == rect.maxy)
+  {
+      OGRGeometryH hSpatialFilterPoint = OGR_G_CreateGeometry( wkbPoint );
 
-  OGR_G_AddPoint_2D( hRing, rect.minx, rect.miny);
-  OGR_G_AddPoint_2D( hRing, rect.maxx, rect.miny);
-  OGR_G_AddPoint_2D( hRing, rect.maxx, rect.maxy);
-  OGR_G_AddPoint_2D( hRing, rect.minx, rect.maxy);
-  OGR_G_AddPoint_2D( hRing, rect.minx, rect.miny);
+      OGR_G_SetPoint_2D( hSpatialFilterPoint, 0, rect.minx, rect.miny );
+
+      OGR_L_SetSpatialFilter( psInfo->hLayer, hSpatialFilterPoint );
 
-  OGR_G_AddGeometryDirectly( hSpatialFilterPolygon, hRing );
+      OGR_G_DestroyGeometry( hSpatialFilterPoint );
+  }
+  else if (rect.minx == rect.maxx || rect.miny == rect.maxy)
+  {
+      OGRGeometryH hSpatialFilterLine = OGR_G_CreateGeometry( wkbLineString );
+
+      OGR_G_AddPoint_2D( hSpatialFilterLine, rect.minx, rect.miny );
+      OGR_G_AddPoint_2D( hSpatialFilterLine, rect.maxx, rect.maxy );
+
+      OGR_L_SetSpatialFilter( psInfo->hLayer, hSpatialFilterLine );
 
-  OGR_L_SetSpatialFilter( psInfo->hLayer, hSpatialFilterPolygon );
+      OGR_G_DestroyGeometry( hSpatialFilterLine );
+  }
+  else
+  {
+      OGRGeometryH hSpatialFilterPolygon = OGR_G_CreateGeometry( wkbPolygon );
+      OGRGeometryH hRing = OGR_G_CreateGeometry( wkbLinearRing );
 
-  OGR_G_DestroyGeometry( hSpatialFilterPolygon );
+      OGR_G_AddPoint_2D( hRing, rect.minx, rect.miny);
+      OGR_G_AddPoint_2D( hRing, rect.maxx, rect.miny);
+      OGR_G_AddPoint_2D( hRing, rect.maxx, rect.maxy);
+      OGR_G_AddPoint_2D( hRing, rect.minx, rect.maxy);
+      OGR_G_AddPoint_2D( hRing, rect.minx, rect.miny);
+
+      OGR_G_AddGeometryDirectly( hSpatialFilterPolygon, hRing );
+
+      OGR_L_SetSpatialFilter( psInfo->hLayer, hSpatialFilterPolygon );
+
+      OGR_G_DestroyGeometry( hSpatialFilterPolygon );
+  }
 
   psInfo->rect = rect;