Поиск городов в километрах "X" (или миль)

Ответ 1

Это запрос MySQL, который будет делать именно то, что вы хотите. Имейте в виду, что такие вещи, как правило, являются приближениями, поскольку земля не является абсолютно сферической, и это не учитывает горы, холмы, долины и т.д. Мы используем этот код на AcademicHomes.com с PHP и MySQL, он возвращает записи в пределах $radius miles $latitude, $longitude.

$res = mysql_query("SELECT
    * 
FROM
    your_table
WHERE
    (
        (69.1 * (latitude - " . $latitude . ")) * 
        (69.1 * (latitude - " . $latitude . "))
    ) + ( 
        (69.1 * (longitude - " . $longitude . ") * COS(" . $latitude . " / 57.3)) * 
        (69.1 * (longitude - " . $longitude . ") * COS(" . $latitude . " / 57.3))
    ) < " . pow($radius, 2) . " 
ORDER BY 
    (
        (69.1 * (latitude - " . $latitude . ")) * 
        (69.1 * (latitude - " . $latitude . "))
    ) + ( 
        (69.1 * (longitude - " . $longitude . ") * COS(" . $latitude . " / 57.3)) * 
        (69.1 * (longitude - " . $longitude . ") * COS(" . $latitude . " / 57.3))
    ) ASC");

Ответ 2

EDIT: Если у вас есть где-то список всех городов мира вместе с их лат. и долго. значения, вы можете выполнить поиск. В этом случае см. Мою первую ссылку ниже для формулы для расчета ширины одной продольной степени на широте :

Честно говоря, осложнения, стоящие за этой проблемой, таковы, что вам будет гораздо лучше использовать такую ​​услугу, как Google Карты, для получения ваших данных. В частности, Земля не является идеальной сферой, и расстояние между двумя градусами изменяется, поскольку вы ближе или ближе от экватора.

См. http://en.wikipedia.org/wiki/Geographic_coordinate_system для примеров того, что я имею в виду, и проверьте API Карт Google.

Ответ 3

В зависимости от того, сколько городов вы включаете, вы можете предварительно скопировать список. Мы делаем это здесь для внутреннего приложения, где погрешность +100 м слишком велика для нашей настройки. Он работает, имея две ключевые таблицы location1, location2, distance. Мы можем быстро отвести местоположение x от места1 очень быстро.

Кроме того, поскольку вычисления могут выполняться в автономном режиме, это не влияет на работу системы. Пользователи также получают более быстрые результаты.

Ответ 4

Я пробовал использовать вышеприведенный код, и ответы были слишком сильными, когда расстояние между точками находилось в диапазоне 20-30 миль, и я в порядке с несколькими минутами ошибки. Поговорил с моим приятелем, и мы придумали это. Код - это python, но вы можете легко его перевести. Чтобы избежать постоянной конверсии в радианы, я переделал свою базу данных, переведя точки lat/lng из градусов в радианы. Приятная часть этого состоит в том, что большая часть математики в основном выполняется один раз.

ra = 3963.1906 # radius @ equator in miles, change to km  if you want distance in km
rb = 3949.90275  # radius @ poles in miles, change to km  if you want distance in km
ra2 = ra * ra
rb2 = rb * rb

phi = self.lat

big_ol_constant = (math.pow(ra2*math.cos(phi), 2) + pow(rb2*math.sin(phi), 2))/ (pow(ra*math.cos(phi), 2) + pow(rb*math.sin(phi), 2))

sqlWhere = "%(distance)g > sqrt((power(lat - %(lat)g,2) + power(lng-%(lng)g,2)) * %(big_ol_constant)g)" % {
    'big_ol_constant': big_ol_constant, 'lat': self.lat, 'lng': self.lng, 'distance': distance}

# This is the Django portion of it, where the ORM kicks in.  sqlWhere is what you would put after the WHERE part of your SQL Query.
qs = ZipData.objects.extra(where=[sqlWhere]);

Кажется, очень точная, когда расстояние между ними невелико и в пределах 10 миль или около того, когда расстояние увеличивается до 200 миль (конечно, к тому времени у вас возникают проблемы с "как вороньи мухи" и "дороги с твердым покрытием" ).

Вот модель ZipData, о которой я упоминал выше.

class ZipData(models.Model):
    zipcode = ZipCodeField(null=False, blank=False, verbose_name="ZipCode", primary_key=True)
    city = models.CharField(max_length=32, null=False, blank=False)
    state = models.CharField(max_length=2)
    lat = models.FloatField(null=False, blank=False)
    lng = models.FloatField(null=False, blank=False)

Дополнительным примечанием является то, что вы можете получить БОЛЬШИЕ геоданные, связанные с почтовыми кодами, в GeoNames.org, и у них даже есть некоторый веб-сервис API, которые вы также можете использовать.

Ответ 5

Существует много (плохих опций)

• Рассчитать расстояние, используя математическую формулу (рассматривать X1-X2 и Y1-Y2) в качестве векторов.

• Создайте таблицу поиска заранее со всеми комбинациями и сохраните расстояния.

• Рассмотрим использование GIS-специфического расширения MySQL. Вот одна статья. Об этом я нашел.

Ответ 6

lessthandot.com фактически имеет 3 разных способа сделать это. вам придется немного прокручивать блоги, но они есть. http://blogs.lessthandot.com/

Ответ 7

Ниже приведена функция nerddinner (пример приложения ASP.NET MVC доступно в codeplex) database (MSSQL).

ALTER FUNCTION [dbo].[DistanceBetween] (@Lat1 as real,
                @Long1 as real, @Lat2 as real, @Long2 as real)
RETURNS real
AS
BEGIN

DECLARE @dLat1InRad as float(53);
SET @dLat1InRad = @Lat1 * (PI()/180.0);
DECLARE @dLong1InRad as float(53);
SET @dLong1InRad = @Long1 * (PI()/180.0);
DECLARE @dLat2InRad as float(53);
SET @dLat2InRad = @Lat2 * (PI()/180.0);
DECLARE @dLong2InRad as float(53);
SET @dLong2InRad = @Long2 * (PI()/180.0);

DECLARE @dLongitude as float(53);
SET @dLongitude = @dLong2InRad - @dLong1InRad;
DECLARE @dLatitude as float(53);
SET @dLatitude = @dLat2InRad - @dLat1InRad;
/* Intermediate result a. */
DECLARE @a as float(53);
SET @a = SQUARE (SIN (@dLatitude / 2.0)) + COS (@dLat1InRad)
                 * COS (@dLat2InRad)
                 * SQUARE(SIN (@dLongitude / 2.0));
/* Intermediate result c (great circle distance in Radians). */
DECLARE @c as real;
SET @c = 2.0 * ATN2 (SQRT (@a), SQRT (1.0 - @a));
DECLARE @kEarthRadius as real;
/* SET kEarthRadius = 3956.0 miles */
SET @kEarthRadius = 6376.5;        /* kms */

DECLARE @dDistance as real;
SET @dDistance = @kEarthRadius * @c;
return (@dDistance);
END

Я предполагаю, что это может быть полезно.

Ответ 8

Вы можете использовать теорему Пифагора для вычисления близости двух пар точек lat/lon.

Если у вас есть два местоположения (Alpha и Beta), вы можете рассчитать их расстояние на расстоянии:

SQRT( POW(Alpha_lat - Beta_lat,2) + POW(Alpha_lon - Beta_lon,2) )

Ответ 9

Используя настройку со следующего URL, Ive построил запрос ниже. (Обратите внимание, что я использую codeIgnitor для запроса базы данных)

http://howto-use-mysql-spatial-ext.blogspot.com/2007/11/using-circular-area-selection.html

function getRadius($point="POINT(-29.8368 30.9096)", $radius=2)
{
    $km = 0.009;
    $center = "GeomFromText('$point')";
    $radius = $radius*$km;
    $bbox = "CONCAT('POLYGON((',
        X($center) - $radius, ' ', Y($center) - $radius, ',',
        X($center) + $radius, ' ', Y($center) - $radius, ',',
        X($center) + $radius, ' ', Y($center) + $radius, ',',
        X($center) - $radius, ' ', Y($center) + $radius, ',',
        X($center) - $radius, ' ', Y($center) - $radius, '
    ))')";

    $query = $this->db->query("
    SELECT id, AsText(latLng) AS latLng, (SQRT(POW( ABS( X(latLng) - X({$center})), 2) + POW( ABS(Y(latLng) - Y({$center})), 2 )))/0.009 AS distance
    FROM crime_listing
    WHERE Intersects( latLng, GeomFromText($bbox) )
    AND SQRT(POW( ABS( X(latLng) - X({$center})), 2) + POW( ABS(Y(latLng) - Y({$center})), 2 )) < $radius
    ORDER BY distance
        ");

    if($query->num_rows()>0){
        return($query->result());
    }else{
        return false;
    }
}

Ответ 10

Не изобретайте велосипед. Это пространственный запрос. Используйте MySQL встроенные пространственные расширения для хранения данных координат широты-долготы в родной геометрии MySQL тип столбца. Затем используйте функцию Distance для запроса точек, находящихся на заданном расстоянии друг от друга.

Отказ от ответственности: это основано на чтении документации, я сам этого не пробовал.