fixed corner case centroid=PIA

8 years ago · 47251daa5f
parent d8604f3c9b
commit 47251daa5f
1 changed files with 1 additions and 0 deletions
--- a/src/pg/sql/13_PIA.sql
+++ b/src/pg/sql/13_PIA.sql
@ -46,6 +46,7 @@ BEGIN
    SELECT array_agg(c) INTO cells FROM c1;

    -- 1st guess: centroid
+    best_c := polygon;
    best_d := cdb_crankshaft._Signed_Dist(polygon, ST_Centroid(Polygon));

    -- looping the loop