У моєму світі за допомогою спеціального SRID (для Карт Google) щось подібне спрацювало:
SELECT * FROM addresses WHERE ST_DWithin(location, ST_SetSRID(ST_MakePoint(longitude, latitude), 3785), radius);
де тип location
- це геометрія (Point, 3785), і longitude
, latitude
і radius
є поплавцями (наприклад, -100, 44, 30 для 100W / 44N / 30 "одиниць" - див. нижче)
Дивіться, який найкращий спосіб знайти усі об’єкти в радіусі іншого об’єкта? у документах Postgis:
ST_DWithin(geometry, geometry, distance)
Функція являє собою зручний спосіб виконання індексованого пошуку відстані. Він працює, створюючи прямокутник пошуку, достатньо великий, щоб укласти радіус відстані, а потім виконувати точний пошук відстані в індексованому підмножині результатів.
ОНОВЛЕННЯ: одиниці не мають миль для SRID 3785 ... вони, здається, або радіани, або градуси, або щось подібне. Але специфікація для мого SRID говорить, що його одиниці є або метрами, або градусами, і це точно не те, принаймні, не без певного перетворення:
alex=# select * from spatial_ref_sys where srid=3785;
srid | auth_name | auth_srid | srtext | proj4text
3785 | EPSG | 3785 | PROJCS["Popular Visualisation CRS / Mercator (deprecated)",GEOGCS["Popular Visualisation CRS",DATUM["Popular_Visualisation_Datum",SPHEROID["Popular Visualisation Sphere",6378137,0,AUTHORITY["EPSG","7059"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6055"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4055"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],PROJECTION["Mercator_1SP"],PARAMETER["central_meridian",0],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs"],AUTHORITY["EPSG","3785"],AXIS["X",EAST],AXIS["Y",NORTH]] | +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs