SQL查询基于lat / long - SQLite查询附近的兴趣点

时间:2010-07-02 20:24:08

标签: sqlite geospatial

给定一个包含三个字段的数据库:

纬度 经度 接近

Lat和Long是GPS坐标,而Proximity是(某个单位 - 英尺?秒?分钟?)

并且考虑到用户当前的GPS纬度/长...

我想写一个SQL查询,它将检索用户在这些行的“Proximity”范围内的所有行。

诀窍:这必须在SQLite中工作,它只支持相当原始的数据类型。没有作弊和依赖SQL Server(或其他提供更好地理空间功能的产品)。

有什么建议吗?

3 个答案:

答案 0 :(得分:7)

这是sqlite的基于C的自定义函数[从下面提到的链接复制]。这可以在iOS应用中使用。它假设您有名为纬度和经度的列,并计算它们与您提供的任何纬度/经度坐标之间的差异。优秀的写作,按原样运作。

#define DEG2RAD(degrees) (degrees * 0.01745327) // degrees * pi over 180

static void distanceFunc(sqlite3_context *context, int argc, sqlite3_value **argv)
{
  // check that we have four arguments (lat1, lon1, lat2, lon2)
  assert(argc == 4);
  // check that all four arguments are non-null
  if (sqlite3_value_type(argv[0]) == SQLITE_NULL || sqlite3_value_type(argv[1]) == SQLITE_NULL ||
  sqlite3_value_type(argv[2]) == SQLITE_NULL ||
  sqlite3_value_type(argv[3]) == SQLITE_NULL) {
  sqlite3_result_null(context);
  return;
}

// get the four argument values
double lat1 = sqlite3_value_double(argv[0]);
double lon1 = sqlite3_value_double(argv[1]);
double lat2 = sqlite3_value_double(argv[2]);
double lon2 = sqlite3_value_double(argv[3]);

// convert lat1 and lat2 into radians now, to avoid doing it twice below
double lat1rad = DEG2RAD(lat1);
double lat2rad = DEG2RAD(lat2);

// apply the spherical law of cosines to our latitudes and longitudes, and set the result appropriately

// 6378.1 is the approximate radius of the earth in kilometres
sqlite3_result_double(context, acos(sin(lat1rad) * sin(lat2rad) + cos(lat1rad) * cos(lat2rad) * cos(DEG2RAD(lon2) - DEG2RAD(lon1))) *
6378.1);
}
     

这定义了一个SQL函数距离(Latitude1,Longitude1,   Latitude2,Longitude2),返回距离(以公里为单位)   两点之间。

     

要使用此功能,请添加以上代码...和   然后在调用sqlite3_open后立即添加此行:

sqlite3_create_function(sqliteDatabasePtr, "distance", 4, SQLITE_UTF8, NULL, &distanceFunc, NULL, NULL);
     

...其中sqliteDatabasePtr是您的调用返回的数据库指针   到sqlite3_open。

     

假设您有一个名为Locations的表,其中包含名称   包含值的纬度和经度(均为double类型)   度,您可以在SQL中使用此函数,如下所示:

SELECT * FROM Locations ORDER BY distance(Latitude, Longitude, 51.503357, -0.1199)
     

此示例根据距离对数据库中的位置进行排序   他们来自伦敦眼,现在是51.503357,-0.1199。

编辑:

原始链接http://www.thismuchiknow.co.uk/?p=71已失效,因此评论中提及的人可以使用以下链接:https://web.archive.org/web/20160808122817/http://www.thismuchiknow.co.uk/?p=71获取该网页

答案 1 :(得分:3)

Haversine formula就是你想要的。使用简单的欧几里德距离公式是不够的,因为地球的曲率会影响两点之间的距离。

Creating a Store Locator with PHP, MySQL & Google Maps是一篇关于用MySQL实现此解决方案的文章,但SQLite不支持trig函数。

在Stack Overflow上查看Calculating Great-Circle Distance with SQLite中的答案,了解有关扩展SQLite函数的更多提示,以便解决此问题。

答案 2 :(得分:2)

您还知道有一个名为Spatialite的SQLite插件,它具有与PostGIS和SQLServer相同的功能吗?我假设您正在尝试在iPhone或浏览器中使用SQLLite。如果不是,我会高度推荐空间。