Changes between Version 1 and Version 2 of UsersWikiExamplesFindNearbyLatLon

Timestamp:

08/11/09 10:45:42 (15 years ago)

Author:

pramsey

Comment:

—

UsersWikiExamplesFindNearbyLatLon

@@ -20 +20 @@
         ' LANGUAGE 'SQL' IMMUTABLE;
 }}}
+If you are willing to accept some slop in the final distance calculation, you can replace it with an approximation based on the latitude of one of the arguments. This will allow the function to accept arguments that are lines and polygons as well as points.
+{{{
+#!sql
+CREATE OR REPLACE FUNCTION ST_DWithin_Sphere(geometry, geometry, float8)
+        RETURNS boolean
+        AS '
+            SELECT $1 && ST_Expand($2,$3 * 1.79866403673916e-05) 
+               AND $2 && ST_Expand($1,$3 * 1.79866403673916e-05) 
+               AND ST_Distance($1, $2) < 
+                        $3 / ST_Distance_Sphere(
+                                   ST_MakePoint(0, ST_Y(ST_Centroid($1))), 
+                                   ST_MakePoint(1, ST_Y(ST_Centroid($1))))
+        ' LANGUAGE 'SQL' IMMUTABLE;
+}}}