[FIXED] Get geographical coordinates in a given square

Issue

My task is as follows: knowing the center (starting point), for example - [{'lat': -7.7940023, 'lng': 110.3656535}] and knowing the radius 5km I need to get all the points included in this square in 1 km increments. How do I achieve this?

P.S Using the Haversine formula I can check if a point is in a given square according to the radius

Image

Solution

if you consider a spherical Earth with radius R, the associated angle of a segment with length L (5km in your case) is:

import numpy as np
R = 6378.0 # km
L = 5.0 # km
angle = np.degrees(L/R)

so now, you can easily check if a point is inside your square:

center = {'lat': -7.7940023, 'lng': 110.3656535}
point = {'lat': 'point latitude', 'lng': 'point longitude'} # insert here your values

if (point['lat']-center['lat'] < angle) and (point['lng']-center['lng'] < angle):
    print('Point is inside')
else:
    print('Point is outside')

EDIT: check the one below.

import numpy as np

R = 6378.0  # km
L = 5  # side of the square
center = {'lat': -7.7940023, 'lng': 110.3656535}

square_side = np.linspace(-L/2, L/2, L+1)
angle = np.degrees(square_side/R)
latitude, longitude = np.meshgrid(angle+center['lat'], angle+center['lng'])
points = []
for lat, lng in zip(latitude.flatten(), longitude.flatten()):
    points.append({'lat': lat, 'lng': lng})

Answered By - mauro