centerLat = 30.0
centerLon = -100.0
N = 10
circlePoints = []
for k in range(N):
    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)
circlePoints.append(point)
print(circlePoints)
python: elegant way of finding the GPS coodinates of a circle around a certain GPS location