From 94497cd79191eca656eafcd86fe0a79114ed8be9 Mon Sep 17 00:00:00 2001
From: abelvm <abel@cartodb.com>
Date: Wed, 17 Aug 2016 12:22:22 -0400
Subject: [PATCH] fix projection issues

---
 src/pg/sql/09_voronoi.sql | 18 ++++++++++++++----
 1 file changed, 14 insertions(+), 4 deletions(-)

diff --git a/src/pg/sql/09_voronoi.sql b/src/pg/sql/09_voronoi.sql
index 970ec2f..476adb4 100644
--- a/src/pg/sql/09_voronoi.sql
+++ b/src/pg/sql/09_voronoi.sql
@@ -12,6 +12,15 @@ RETURNS geometry  AS $$
 DECLARE
     geomout geometry;
 BEGIN
+    -- we need to make the geometry calculations in (pseudo)meters!!!
+    with a as (
+        SELECT unnest(geomin) as g1
+    ),
+    b as(
+        SELECT st_transform(g1, 3857) g2 from a
+    )
+    SELECT array_agg(g2) INTO geomin from b;
+
     WITH
     convexhull_1 as (
         SELECT
@@ -138,8 +147,9 @@ BEGIN
         FROM voro_cells v
     )
     SELECT
-        st_collect(ST_intersection(c.g, v.g))
-       -- ST_intersection(c.g, v.g)
+        st_collect(
+            ST_Transform(ST_intersection(c.g, v.g), 4326)
+        )
     INTO geomout
     FROM
         voro_set v,
@@ -186,8 +196,8 @@ BEGIN
       RETURN NULL;
    END IF;
    IF ( ST_GeometryType(p_pt1) <> 'ST_Point' OR
-        ST_GeometryType(p_pt1) <> 'ST_Point' OR
-        ST_GeometryType(p_pt1) <> 'ST_Point' ) THEN
+        ST_GeometryType(p_pt2) <> 'ST_Point' OR
+        ST_GeometryType(p_pt3) <> 'ST_Point' ) THEN
       RAISE EXCEPTION 'All supplied geometries must be points.';
       RETURN NULL;
    END IF;