我有这个查询，在MySQL中可以正常工作。
SELECT ((ACOS(SIN(12.345 * PI() / 180) * SIN(lat * PI() / 180) +
COS(12.345 * PI() / 180) * COS(lat * PI() / 180) * COS((67.89 - lon) *
PI() / 180)) * 180 / PI()) * 60 * 1.1515 * 1.609344) AS distance, poi.*
FROM poi
WHERE lang='eng'
HAVING distance<='30'
距离以公里为单位，输入参数为lat=12.345和lon=67.89
SQLite版本为3，在Android上无法运行自定义函数。因为它不是标准SQLite的一部分，所以我也没有acos()等函数可用。
如何在SQLite中执行上述查询？
这是一个用Java实现的在Android设备上构建基于位置查询的例子。
这个想法来自于KennyTM（见被接受的回答），它需要在你的表中添加4列，以存储纬度和经度的正弦和余弦值。
以下是为“商店”表准备插入数据的代码：
public static void injectLocationValues(ContentValues values, double latitude, double longitude) {
values.put(LocationColumns.LATITUDE, latitude);
values.put(LocationColumns.LONGITUDE, longitude);
values.put(LocationColumns.COSLAT, Math.cos(MathUtil.deg2rad(latitude)));
values.put(LocationColumns.SINLAT, Math.sin(MathUtil.deg2rad(latitude)));
values.put(LocationColumns.COSLNG, Math.cos(MathUtil.deg2rad(longitude)));
values.put(LocationColumns.SINLNG, Math.sin(MathUtil.deg2rad(longitude)));
}
public static double deg2rad(double deg) {
return (deg * Math.PI / 180.0);
}
你可以使用以下函数来构建你的投影：
/**
* Build query based on distance using spherical law of cosinus
*
* d = acos(sin(lat1).sin(lat2)+cos(lat1).cos(lat2).cos(long2−long1)).R
* where R=6371 and latitudes and longitudes expressed in radians
*
* In Sqlite we do not have access to acos() sin() and lat() functions.
* Knowing that cos(A-B) = cos(A).cos(B) + sin(A).sin(B)
* We can determine a distance stub as:
* d = sin(lat1).sin(lat2)+cos(lat1).cos(lat2).(cos(long2).cos(long1)+sin(long2).sin(long1))
*
* First comparison point being fixed, sin(lat1) cos(lat1) sin(long1) and cos(long1)
* can be replaced by constants.
*
* Location aware table must therefore have the following columns to build the equation:
* sinlat => sin(radians(lat))
* coslat => cos(radians(lat))
* coslng => cos(radians(lng))
* sinlng => sin(radians(lng))
*
* Function will return a real between -1 and 1 which can be used to order the query.
* Distance in km is after expressed from R.acos(result)
*
* @param latitude, latitude of search
* @param longitude, longitude of search
* @return selection query to compute the distance
*/
public static String buildDistanceQuery(double latitude, double longitude) {
final double coslat = Math.cos(MathUtil.deg2rad(latitude));
final double sinlat = Math.sin(MathUtil.deg2rad(latitude));
final double coslng = Math.cos(MathUtil.deg2rad(longitude));
final double sinlng = Math.sin(MathUtil.deg2rad(longitude));
//@formatter:off
return "(" + coslat + "*" + LocationColumns.COSLAT
+ "*(" + LocationColumns.COSLNG + "*" + coslng
+ "+" + LocationColumns.SINLNG + "*" + sinlng
+ ")+" + sinlat + "*" + LocationColumns.SINLAT
+ ")";
//@formatter:on
}
public static double convertPartialDistanceToKm(double result) {
return Math.acos(result) * 6371;
}
如果你想使用偏移距离（partial distance）排序查询结果，那么你需要按照降序（DESC）而不是升序（ASC）的方式排序。
在使用iOS的sqlite3时，我遇到了同样的问题，在尝试了一些公式后，我找到了一种不需要使用sql端函数的方法（伪代码）：
为存储在数据库中的每个项目预先计算这些值（并存储它们）：
cos_lat = cos(lat * PI / 180)
sin_lat = sin(lat * PI / 180)
cos_lng = cos(lng * PI / 180)
sin_lng = sin(lng * PI / 180)
在搜索时预先计算这些值（针对给定的位置cur_lat、cur_lng）
CUR_cos_lat = cos(cur_lat * PI / 180)
CUR_sin_lat = sin(cur_lat * PI / 180)
CUR_cos_lng = cos(cur_lng * PI / 180)
CUR_sin_lng = sin(cur_lng * PI / 180)
cos_allowed_distance = cos(2.0 / 6371) # This is 2km
您的SQL查询将类似于以下内容（将CUR_*替换为刚刚计算的值）
SELECT * FROM position WHERE CUR_sin_lat * sin_lat + CUR_cos_lat * cos_lat * (cos_lng* CUR_cos_lng + sin_lng * CUR_sin_lng) > cos_allowed_distance;
lat和lon的正弦和余弦值。由于cos(a+b) = cos a cos b - sin a sin b，以及其他出现在程序中的像SIN(12.345 * PI() / 180)这样的sin和cos值可以在运行查询之前在程序中计算，所以大的"distance"表达式可以简化为形如P * SIN_LAT + Q * COS_LAT + ...的形式，这些可以由SQLite3处理。org.sqlite不在Android中，并且不清楚是否存在一个org.sqlite实现可以与Android的SQLite环境协同工作。 - Pentium10ACOS()？使用 sqlite3_create_function() 创建自定义 SQL 函数？在 Android 中该怎么做？ - Afriza N. Arief