Re: [postgis-users] ST_Intersects differing results

2016-05-25 Thread Paul Ramsey
If your data model expects strict containment (all parcels must be in only one zone) then the simplest thing is to proxy the parcel polygon with a point, then you get a guaranteed 1:1 relationship between parcel and zone. Using ST_PointOnSurface() is required since ST_Centroid() is not guaranteed t

Re: [postgis-users] ST_Intersects differing results

2016-05-25 Thread Rémi Cura
Hey, if you want to do it the robust way, you can use a shared surface distance : for instance, ST_area(ST_Intersection(A,B)) / ST_Area(A) where A is you polygon, and B your zoning information. If everything is snapped, you could also use topology. Cheers, Rémi-C 2016-05-25 3:45 GMT+02:00 Paul