Algorithme pour trouver tous Latitude Longitude endroits à l'intérieur d'une certaine distance à partir d'un Lat de Gnl emplacement

Commençons par Comparer la distance entre les latitudes. Chaque degré de latitude est d'environ 69 km (111 kilomètres) de distance. La gamme varie (en raison de la terre légèrement la forme ellipsoïde) à partir de 68.703 miles (110.567 km) au niveau de l'équateur à 69.407 (111.699 km) dans les pôles. La distance entre les deux endroits sera égale ou plus grande que la distance entre leurs latitudes.

Noter que ce n'est pas vrai pour les longitudes - la longueur de chaque degré de longitude est, dépend de la latitude. Toutefois, si vos données est limité à une certaine zone (un seul pays, par exemple) - vous pouvez calculer une valeur minimale et maximale des valeurs limites pour les longitudes ainsi.

---

Continuer une faible précision, rapide calcul de la distance que suppose sphérique de la terre:

Le grand cercle de la distance d entre deux points de coordonnées {lat1,lon1} et {lat2,lon2} est donnée par:

d = acos(sin(lat1)*sin(lat2)+cos(lat1)*cos(lat2)*cos(lon1-lon2))

Une mathématiquement équivalent de la formule, ce qui est moins sujet à erreur d'arrondi pour les courtes distances est:

d = 2*asin(sqrt((sin((lat1-lat2)/2))^2 +
    cos(lat1)*cos(lat2)*(sin((lon1-lon2)/2))^2))

d est la distance (en radians

distance_km ≈ radius_km * distance_radians ≈ 6371 * d

(6371 km à l'est de la rayon moyen de la terre)

Cette méthode de calcul des exigences sont mimimal. Cependant, le résultat est très précis pour les petites distances.

---

Puis, si elle est à une distance donnée, plus ou moins, utiliser une méthode plus précise.

GeographicLib est la mise en œuvre plus précise, je sais, mais Vincenty formule d'inversion peut être utilisé en tant que bien.

---

Si vous utilisez un SGBDR, définir la latitude de la clé primaire et la longitude comme une clé secondaire. Requête pour un intervalle de latitude, ou pour une latitude/longitude), comme il est décrit ci-dessus, puis calculer les distances exactes de l'ensemble des résultats.

Noter que les versions modernes de tous les principaux Sgbdr support géographique, les types de données et des requêtes en mode natif.

PostgreSQL SIG extensions pourraient être utiles - comme dans, il peut déjà mettre en œuvre la plupart des fonctionnalités que vous pensez de la mise en œuvre.

Réservoirs Yogihosting

J'ai dans ma base de données, un ensembles de tables à partir d'Ouvrir Streep Cartes et j'ai testé succès.

Travail à Distance fine en mètres.

SET @orig_lat=-8.116137;
SET @orig_lon=-34.897488;
SET @dist=1000;

SELECT *,(((acos(sin((@orig_lat*pi()/180)) * sin((dest.latitude*pi()/180))+cos((@orig_lat*pi()/180))*cos((dest.latitude*pi()/180))*cos(((@orig_lon-dest.longitude)*pi()/180))))*180/pi())*60*1.1515*1609.344) as distance FROM nodes AS dest HAVING distance < @dist ORDER BY distance ASC LIMIT 100;

Comme biziclop mentionné, une sorte d'espace métrique arbre serait probablement votre meilleure option. J'ai de l'expérience à l'aide de kd-tree et quad arbres pour faire ces sortes de requêtes de plage et ils sont étonnamment rapide; ils ne sont pas non plus difficile à écrire. Je vous suggère de regarder dans l'un de ces structures, comme ils vous permettent également de répondre à d'autres questions intéressantes comme "quel est le point le plus proche de mon jeu de données à ce point?"

Vous pouvez convertir latitude-longitude au format UTM qui est format métrique qui peuvent vous aider à calculer les distances. Ensuite, vous pouvez facilement décider si le point tombe dans l'emplacement spécifique.

vous pouvez vérifier cette équation
je pense que ça va aider

SELECT id, ( 3959 * acos( cos( radians(37) ) * cos( radians( lat ) ) * cos( radians( lng ) - radians(-122) ) + sin( radians(37) ) * sin( radians( lat ) ) ) ) AS distance FROM markers HAVING distance < 25 ORDER BY distance LIMIT 0 , 20;