SQlite获取最近的位置(经度和纬度)

Jur*_*ure 82 java sqlite android geolocation

我的SQLite数据库中存储了纬度和经度数据,我希望得到最接近我输入参数的位置(例如我的当前位置 - lat/lng等).

我知道这在MySQL中是可能的,我已经做了很多研究,SQLite需要为Haversine公式定义外部函数(计算球体上的距离),但我还没有找到任何用Java编写的工作.

此外,如果我想添加自定义函数,我需要org.sqlite.jar(for org.sqlite.Function),这会给应用程序增加不必要的大小.

另一方面,我需要SQL中的Order by函数,因为单独显示距离并不是一个大问题 - 我已经在我的自定义SimpleCursorAdapter中做了,但是我无法对数据进行排序,因为我我的数据库中没有距离列.这意味着每次更改位置时都会更新数据库,这会浪费电池和性能.因此,如果有人知道使用不在数据库中的列对游标进行排序,我将不胜感激!

我知道有很多Android应用程序使用这个功能,但有人可以解释一下魔术.

顺便说一句,我找到了另一种方法:查询在SQLite中获取基于Radius的记录?

它建议为lat和lng的cos和sin值创建4个新列,但是还有其他任何不那么冗余的方法吗?

Bob*_*obs 107

1)首先使用良好的近似值过滤SQLite数据,并减少需要在java代码中评估的数据量.为此,请使用以下过程:

要获得确定性阈值和更准确的数据过滤,最好在java代码中计算中心点的北,西,东和南米的4个位置,然后通过小于和大于SQL运算符(>,<)确定数据库中的点是否在该矩形中.radius

该方法calculateDerivedPosition(...)为您计算这些点(图中的p1,p2,p3,p4).

在此输入图像描述

/**
* Calculates the end-point from a given source at a given range (meters)
* and bearing (degrees). This methods uses simple geometry equations to
* calculate the end-point.
* 
* @param point
*           Point of origin
* @param range
*           Range in meters
* @param bearing
*           Bearing in degrees
* @return End-point from the source given the desired range and bearing.
*/
public static PointF calculateDerivedPosition(PointF point,
            double range, double bearing)
    {
        double EarthRadius = 6371000; // m

        double latA = Math.toRadians(point.x);
        double lonA = Math.toRadians(point.y);
        double angularDistance = range / EarthRadius;
        double trueCourse = Math.toRadians(bearing);

        double lat = Math.asin(
                Math.sin(latA) * Math.cos(angularDistance) +
                        Math.cos(latA) * Math.sin(angularDistance)
                        * Math.cos(trueCourse));

        double dlon = Math.atan2(
                Math.sin(trueCourse) * Math.sin(angularDistance)
                        * Math.cos(latA),
                Math.cos(angularDistance) - Math.sin(latA) * Math.sin(lat));

        double lon = ((lonA + dlon + Math.PI) % (Math.PI * 2)) - Math.PI;

        lat = Math.toDegrees(lat);
        lon = Math.toDegrees(lon);

        PointF newPoint = new PointF((float) lat, (float) lon);

        return newPoint;

    }
Run Code Online (Sandbox Code Playgroud)

现在创建您的查询:

PointF center = new PointF(x, y);
final double mult = 1; // mult = 1.1; is more reliable
PointF p1 = calculateDerivedPosition(center, mult * radius, 0);
PointF p2 = calculateDerivedPosition(center, mult * radius, 90);
PointF p3 = calculateDerivedPosition(center, mult * radius, 180);
PointF p4 = calculateDerivedPosition(center, mult * radius, 270);

strWhere =  " WHERE "
        + COL_X + " > " + String.valueOf(p3.x) + " AND "
        + COL_X + " < " + String.valueOf(p1.x) + " AND "
        + COL_Y + " < " + String.valueOf(p2.y) + " AND "
        + COL_Y + " > " + String.valueOf(p4.y);
Run Code Online (Sandbox Code Playgroud)

COL_X是数据库中存储纬度值的列的名称,COL_Y用于经度.

因此,您可以在中心点附近获得一些数据,并且具有良好的近似值.

2)现在,您可以循环使用这些过滤后的数据,并使用以下方法确定它们是否真正靠近您的点(在圆圈中):

public static boolean pointIsInCircle(PointF pointForCheck, PointF center,
            double radius) {
        if (getDistanceBetweenTwoPoints(pointForCheck, center) <= radius)
            return true;
        else
            return false;
    }

public static double getDistanceBetweenTwoPoints(PointF p1, PointF p2) {
        double R = 6371000; // m
        double dLat = Math.toRadians(p2.x - p1.x);
        double dLon = Math.toRadians(p2.y - p1.y);
        double lat1 = Math.toRadians(p1.x);
        double lat2 = Math.toRadians(p2.x);

        double a = Math.sin(dLat / 2) * Math.sin(dLat / 2) + Math.sin(dLon / 2)
                * Math.sin(dLon / 2) * Math.cos(lat1) * Math.cos(lat2);
        double c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a));
        double d = R * c;

        return d;
    }
Run Code Online (Sandbox Code Playgroud)

请享用!

我使用并定制了这个参考并完成了它.


小智 68

Chris的答案非常有用(谢谢!),但只有在使用直线坐标(例如UTM或OS网格参考)时才能使用.如果使用lat/lng的度数(例如WGS84),则上述仅适用于赤道.在其他纬度,您需要减少经度对排序顺序的影响.(想象一下,你接近北极......纬度仍然与任何地方相同,但经度可能只有几英尺.这意味着排序顺序不正确).

如果您不在赤道,请根据您当前的纬度预先计算软糖因子:

<fudge> = Math.pow(Math.cos(Math.toRadians(<lat>)),2);
Run Code Online (Sandbox Code Playgroud)

然后订购:

((<lat> - LAT_COLUMN) * (<lat> - LAT_COLUMN) + (<lng> - LNG_COLUMN) * (<lng> - LNG_COLUMN) * <fudge>)

它仍然只是一个近似值,但比第一个好得多,因此排序顺序的不准确性将更加罕见.

  • 这是一个非常有趣的观点,关于纵向线在两极会聚并使结果偏离你越近.很好的修复. (3认同)

Chr*_*son 67

我知道这已经得到了回答和接受,但我想我会添加我的经验和解决方案.

虽然我很高兴在设备上执行半正功能来计算用户当前位置与任何特定目标位置之间的准确距离,但仍需要按距离的顺序对查询结果进行排序和限制.

不太令人满意的解决方案是返回批次并在事后进行排序和过滤,但这会导致第二个游标并返回并丢弃许多不必要的结果.

我首选的解决方案是传递long和lats的平方delta值的排序顺序:

((<lat> - LAT_COLUMN) * (<lat> - LAT_COLUMN) +
 (<lng> - LNG_COLUMN) * (<lng> - LNG_COLUMN))
Run Code Online (Sandbox Code Playgroud)

没有必要只为排序顺序执行完整的hasrsine,并且不需要将结果平方根,因此SQLite可以处理计算.

编辑:

这个答案仍在接受爱.它在大多数情况下都能正常工作,但是如果你需要更高的准确度,请查看@Teasel下面的答案,它增加了一个"软糖"因子,可以解决随着纬度接近90而增加的不准确性.

  • 在我的应用程序BostonBusMap中,我使用此解决方案显示最接近当前位置的停靠点.但是你需要通过"cos(纬度)"来缩放经度距离,使纬度和经度大致相等.见http://en.wikipedia.org/wiki/Geographical_distance#Spherical_Earth_projected_to_a_plane (6认同)
  • @iMatoria - 这只是毕达哥拉斯着名定理的缩减版.给定两组坐标,两个X值之间的差异表示直角三角形的一侧,Y值之间的差异是另一侧.要获得斜边(以及点之间的距离),可以将这两个值的平方加在一起,然后将结果平方根.在我们的例子中,我们不做最后一点(平方根),因为我们不能.幸运的是,这不是排序顺序所必需的. (3认同)