python: une façon élégante de trouver les coordonnées GPS d'un cercle autour d'un certain emplacement GPS

Question

J'ai un ensemble de coordonnées GPS en notation décimale, et je cherche un moyen de trouver les coordonnées dans un cercle avec un rayon variable autour chaque emplacement.

Here is an example de ce dont j'ai besoin. C'est un cercle avec 1km rayon autour de la coordonnée 47,11. Ce dont j'ai besoin est l'algorithme pour trouver les coordonnées du cercle, donc je peux l'utiliser dans mon fichier kml en utilisant un polygone. Idéalement pour python.

Des idées?

Answer 1

Utilisez la formule pour "distance donnée du point de destination et portant du point de départ" ici:

http://www.movable-type.co.uk/scripts/latlong.html

avec votre point central comme point de départ, votre rayon comme la distance, et boucle sur un certain nombre de roulements de 0 degrés à 360 degrés. Cela vous donnera les points sur un cercle, et travaillera aux pôles parce qu'il utilise de grands cercles partout.

Answer 2

voir aussi Adding distance to a GPS coordinate pour les relations simples entre les distances lat/lon et courte distance.

cela fonctionne:

import math 

# inputs 
radius = 1000.0 # m - the following code is an approximation that stays reasonably accurate for distances < 100km 
centerLat = 30.0 # latitude of circle center, decimal degrees 
centerLon = -100.0 # Longitude of circle center, decimal degrees 

# parameters 
N = 10 # number of discrete sample points to be generated along the circle 

# generate points 
circlePoints = [] 
for k in xrange(N): 
    # compute 
    angle = math.pi*2*k/N 
    dx = radius*math.cos(angle) 
    dy = radius*math.sin(angle) 
    point = {} 
    point['lat']=centerLat + (180/math.pi)*(dy/6378137) 
    point['lon']=centerLon + (180/math.pi)*(dx/6378137)/math.cos(centerLat*math.pi/180) 
    # add to list 
    circlePoints.append(point) 

print circlePoints 

Answer 3

C'est un problème de trigonométrie simple.

Définissez votre système de coordonnées XOY au centre de votre cercle. Commencez par y = 0 et trouvez votre valeur x avec x = r. Ensuite, il suffit de faire pivoter votre rayon autour de l'origine par l'angle a (en radians). Vous pouvez trouver les coordonnées de votre prochain point sur le cercle avec Xi = r * cos(a), Yi = r * sin(a). Répétez le dernier 2 * Pi/a fois.

C'est tout.

MISE À JOUR

Prendre le commentaire de @poolie en compte, le problème peut être résolu de la manière suivante (en supposant que la Terre étant la sphère droite). Considérons une section transversale de la Terre avec son plus grand diamètre D à travers notre point (appelez-le L). Le diamètre de 1 km de notre cercle devient alors un accord (appelez le AB) du cercle de la section transversale de la Terre. Ainsi, la longueur de l'arc AB devient (AB) = D * Theta, où Theta = 2 * sin(|AB|/2). En outre, il est facile de trouver toutes les autres dimensions.