rtree: prepare for future modifications

git-svn-id: https://svn.osgeo.org/grass/grass/trunk@52519 15284696-431f-4ddb-bdfa-cd5b030d7da7

lib/vector/rtree/index.c

@@ -207,9 +207,15 @@ void RTreeFreeIndex(struct RTree *t)
  * Search in an index tree for all data retangles that
  * overlap or touch the argument rectangle.
  * Return the number of qualifying data rects.
+ * 
+ * add option to select operator to select rectangles ?
+ * current: overlap
+ * possible alternatives: 
+ *  - select all rectangles that are fully contained in r
+ *  - select all rectangles that fully contain r
  */
-int RTreeSearch(struct RTree *t, struct RTree_Rect *r, SearchHitCallback *shcb,
-		void *cbarg)
+int RTreeSearch(struct RTree *t, struct RTree_Rect *r,
+                SearchHitCallback *shcb, void *cbarg)
 {
     assert(r && t);
 

lib/vector/rtree/index.h

@@ -59,15 +59,15 @@ void RTreeReInsertNode(struct RTree_Node *, struct RTree_ListNode **);
 void RTreeFreeListBranch(struct RTree_ListBranch *);
 
 /* indexm.c */
-int RTreeSearchM(struct RTree *, struct RTree_Rect *, SearchHitCallback *,
-		       void *);
+int RTreeSearchM(struct RTree *, struct RTree_Rect *,
+                 SearchHitCallback *, void *);
 int RTreeInsertRectM(struct RTree_Rect *, union RTree_Child, int, struct RTree *);
 int RTreeDeleteRectM(struct RTree_Rect *, union RTree_Child, struct RTree *);
 int RTreeValidChildM(union RTree_Child *child);
 
 /* indexf.c */
-int RTreeSearchF(struct RTree *, struct RTree_Rect *, SearchHitCallback *,
-		       void *);
+int RTreeSearchF(struct RTree *, struct RTree_Rect *,
+                 SearchHitCallback *, void *);
 int RTreeInsertRectF(struct RTree_Rect *, union RTree_Child, int, struct RTree *);
 int RTreeDeleteRectF(struct RTree_Rect *, union RTree_Child, struct RTree *);
 int RTreeValidChildF(union RTree_Child *);

lib/vector/rtree/rect.c

@@ -436,11 +436,38 @@ int RTreeOverlap(struct RTree_Rect *r, struct RTree_Rect *s, struct RTree *t)
 
 
 /*-----------------------------------------------------------------------------
-| Decide whether rectangle r is contained in rectangle s.
+| Decide whether rectangle s is contained in rectangle r.
 -----------------------------------------------------------------------------*/
 int RTreeContained(struct RTree_Rect *r, struct RTree_Rect *s, struct RTree *t)
 {
-    register int i, j, result;
+    register int i, j;
+
+    /* assert(r && s); */
+
+    /* undefined rect is contained in any other */
+    if (Undefined(r, t))
+	return TRUE;
+
+    /* no rect (except an undefined one) is contained in an undef rect */
+    if (Undefined(s, t))
+	return FALSE;
+
+    for (i = 0; i < t->ndims; i++) {
+	j = i + t->ndims_alloc;	/* index for high sides */
+	if (s->boundary[i] < r->boundary[i] ||
+	    s->boundary[j] > r->boundary[j])
+	    return FALSE;
+    }
+    return TRUE;
+}
+
+
+/*-----------------------------------------------------------------------------
+| Decide whether rectangle s fully contains rectangle r.
+-----------------------------------------------------------------------------*/
+int RTreeContains(struct RTree_Rect *r, struct RTree_Rect *s, struct RTree *t)
+{
+    register int i, j;
 
     /* assert(r && s); */
 
@@ -452,11 +479,11 @@ int RTreeContained(struct RTree_Rect *r, struct RTree_Rect *s, struct RTree *t)
     if (Undefined(s, t))
 	return FALSE;
 
-    result = TRUE;
     for (i = 0; i < t->ndims; i++) {
 	j = i + t->ndims_alloc;	/* index for high sides */
-	result = result && r->boundary[i] >= s->boundary[i]
-	    && r->boundary[j] <= s->boundary[j];
+	if (s->boundary[i] > r->boundary[i] ||
+	    s->boundary[j] < r->boundary[j])
+	    return FALSE;
     }
-    return result;
+    return TRUE;
 }

lib/vector/rtree/rtree.h

@@ -195,13 +195,15 @@ struct RTree
 };
 
 /* RTree main functions */
-int RTreeSearch(struct RTree *, struct RTree_Rect *, SearchHitCallback *,
-		       void *);
+int RTreeSearch(struct RTree *, struct RTree_Rect *,
+                SearchHitCallback *, void *);
 int RTreeInsertRect(struct RTree_Rect *, int, struct RTree *);
 int RTreeDeleteRect(struct RTree_Rect *, int, struct RTree *);
 struct RTree *RTreeNewIndex(int, off_t, int);
 void RTreeFreeIndex(struct RTree *);
 int RTreeOverlap(struct RTree_Rect *, struct RTree_Rect *, struct RTree *);
+int RTreeContained(struct RTree_Rect *, struct RTree_Rect *, struct RTree *);
+int RTreeContains(struct RTree_Rect *, struct RTree_Rect *, struct RTree *);
 
 /* RTree node management */
 struct RTree_Node *RTreeNewNode(struct RTree *, int);