30.293400N,30.848787E

python
import sys, os, math, json, pkgutil, glob
mods = ['geopandas','shapely','geopy','reverse_geocoder','timezonefinder','cartopy','fiona','pyproj']
{k: pkgutil.find_loader(k) is not None for k in mods}


{'geopandas': True,
 'shapely': True,
 'geopy': True,
 'reverse_geocoder': False,
 'timezonefinder': False,
 'cartopy': False,
 'fiona': True,
 'pyproj': True}

python
import geopandas as gpd, os, glob, pkgutil, pathlib, json, shapely
try:
    print(gpd.datasets.available)
except Exception as e: print("datasets error", e)
try:
    path = gpd.datasets.get_path('naturalearth_lowres')
    print(path)
except Exception as e: print("get_path error", e)

python
from math import radians, sin, cos, sqrt, atan2
lat1, lon1 = 30.2934, 30.848787
places = {
"Cairo": (30.0444,31.2357),
"Giza": (30.0131,31.2089),
"Al Khatatbah Air Base": (30.3111,30.8764),
"Al-Khatatba village": (30.3595172,30.8176582),
"Khatatba station": (30.35,30.8167),
"Sheikh Zayed": (30.047,30.976),
"6th of October City": (29.9285,30.9188),
"Sadat City": (30.367,30.504),
"Ashmoun": (30.297,30.976)
}
def hav(a,b,c,d):
    R=6371
    phi1,phi2=radians(a),radians(c)
    dphi=radians(c-a); dl=radians(d-b)
    x=sin(dphi/2)**2+cos(phi1)*cos(phi2)*sin(dl/2)**2
    return 2*R*atan2(sqrt(x),sqrt(1-x))
for k,v in places.items():
    print(k, hav(lat1,lon1,*v))
# DMS convert
def dms(deg):
    d=int(abs(deg))
    m=int((abs