From 3655ec0ab249b60301244a4f030eeebf20457abb Mon Sep 17 00:00:00 2001
From: Paul Ramsey <pramsey@opengeo.org>
Date: Fri, 31 May 2013 16:12:08 -0700
Subject: [PATCH] Rejigging

---
 lib/pc_api_internal.h       |  8 +++++
 lib/pc_filter.c             | 58 +++++++++++++++++++++++++++++++++++++
 lib/pc_patch_dimensional.c  |  3 +-
 lib/pc_patch_uncompressed.c |  1 +
 lib/pc_stats.c              |  8 +++--
 5 files changed, 75 insertions(+), 3 deletions(-)
 create mode 100644 lib/pc_filter.c

diff --git a/lib/pc_api_internal.h b/lib/pc_api_internal.h
index facfc2e..ebe3eba 100644
--- a/lib/pc_api_internal.h
+++ b/lib/pc_api_internal.h
@@ -83,6 +83,13 @@ enum DIMCOMPRESSIONS {
     PC_DIM_ZLIB = 3
 };
 
+typedef enum {
+    PC_GT,
+    PC_LT,
+    PC_EQUALS,
+    PC_BETWEEN
+} PC_FILTERTYPE;
+
 /** What is the endianness of this system? */
 char machine_endian(void);
 
@@ -222,6 +229,7 @@ uint64_t pc_bytes_sigbits_count_64(const PCBYTES *pcb, uint32_t *nsigbits);
 /** Initialize with very large mins and very small maxes */
 void pc_bounds_init(PCBOUNDS *b);
 PCSTATS* pc_stats_clone(const PCSTATS *stats);
+size_t pc_stats_size(const PCSCHEMA *schema);
 
 
 #endif /* _PC_API_INTERNAL_H */
\ No newline at end of file
diff --git a/lib/pc_filter.c b/lib/pc_filter.c
new file mode 100644
index 0000000..8fdad5f
--- /dev/null
+++ b/lib/pc_filter.c
@@ -0,0 +1,58 @@
+/***********************************************************************
+* pc_filter.c
+*
+*  Pointclound patch filtering.
+*
+*  Copyright (c) 2013 OpenGeo
+*
+***********************************************************************/
+
+#include "pc_api_internal.h"
+
+typedef PC_BITMAP uint8_t;
+
+static PC_BITMAP *
+pc_bitmap_new(uint32_t npoints)
+{
+    if ( npoints == 0 ) 
+        return NULL
+    else
+        return (PC_BITMAP*)pcalloc(sizeof(PC_BITMAP)*npoints);
+}
+static void
+pc_bitmap_free(PCBITMAP *map)
+{
+    pcfree(map);
+}
+
+static PC_BITMAP *
+pc_patch_uncompressed_bitmap(const PCPATCH *pa, PC_FILTERTYPE filter, double val1, double val2)
+{
+    int i;
+    PCPOINT pt;
+    pt.readonly = PC_TRUE;
+    pt.schema = pa->schema;
+    uint8_t *buf = pa->data
+    
+    for ( i = 0; i < pa->npoints; i++ )
+    {
+        pt.data = pa->
+}
+
+PC_BITMAP *
+pc_patch_bitmap(const PCPATCH *pa, PC_FILTERTYPE filter, double val1, double val2)
+{
+    if ( ! pa ) return NULL;
+	switch ( pa->type )
+	{
+		case PC_NONE:
+			return pc_patch_uncompressed_bitmap((PCPATCH_UNCOMPRESSED*)pa, filter, val1, val2);
+		case PC_GHT:
+            // return pc_patch_ght_bitmap((PCPATCH_GHT*)pa, filter, val1, val2);
+		case PC_DIMENSIONAL:
+            // return pc_patch_dimensional_bitmap((PCPATCH_DIMENSIONAL*)pa, filter, val1, val2);
+		default:
+            pc_error("%s: failure", __func__);
+	}
+	return NULL;
+}
diff --git a/lib/pc_patch_dimensional.c b/lib/pc_patch_dimensional.c
index 12cea2e..4f88a9b 100644
--- a/lib/pc_patch_dimensional.c
+++ b/lib/pc_patch_dimensional.c
@@ -285,4 +285,5 @@ pc_patch_dimensional_from_pointlist(const PCPOINTLIST *pdl)
     PCPATCH_DIMENSIONAL *dimpatch = pc_patch_dimensional_from_uncompressed(patch);
     pc_patch_uncompressed_free(patch);
     return dimpatch;
-}
\ No newline at end of file
+}
+
diff --git a/lib/pc_patch_uncompressed.c b/lib/pc_patch_uncompressed.c
index daba9e9..95804f7 100644
--- a/lib/pc_patch_uncompressed.c
+++ b/lib/pc_patch_uncompressed.c
@@ -405,3 +405,4 @@ pc_patch_uncompressed_add_point(PCPATCH_UNCOMPRESSED *c, const PCPOINT *p)
 
 	return PC_SUCCESS;
 }
+
diff --git a/lib/pc_stats.c b/lib/pc_stats.c
index 56d75e3..f316be4 100644
--- a/lib/pc_stats.c
+++ b/lib/pc_stats.c
@@ -1,8 +1,7 @@
 /***********************************************************************
 * pc_stats.c
 *
-*  Pointclound schema handling. Parse and emit the XML format for
-*  representing packed multidimensional point data.
+*  Pointclound patch statistics generation.
 *
 *  Copyright (c) 2013 OpenGeo
 *
@@ -203,6 +202,11 @@ pc_patch_uncompressed_compute_stats(PCPATCH_UNCOMPRESSED *pa)
     return PC_SUCCESS;
 }
 
+size_t
+pc_stats_size(const PCSCHEMA *schema)
+{
+    return 3*schema->size;
+}