Запрос на получение записей на основе Radius в SQLite?

Ответ 1

Вы можете создать 4 новых столбца: sin и cos lat и lon. Поскольку cos(a+b) = cos a cos b - sin a sin b, а другие появления sin и cos, такие как SIN(12.345 * PI() / 180), могут быть вычислены в программе перед запуском запроса, большое выражение "distance" сводится к чему-то вроде формы P * SIN_LAT + Q * COS_LAT + ..., которую может обрабатывать SQLite3.

BTW, см. также Sqlite на Android: как создать функцию sqlite dist db - для использования в приложении для расчета расстояния с использованием lat, long.

Ответ 2

Вот реализация в 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;
}

Если вы хотите заказать свой запрос с использованием частичного расстояния, вам необходимо заказать DESC, а не ASC.

Ответ 3

Если во время работы над sqlite3 для ios такая же проблема, если немного поиграть с формулой, это способ сделать это без использования функции со стороны 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;