Como calcular a distância entre dois pontos na terra em Python

# vamos importar o módulo Math
import math

# função que recebe dois pontos na terra e retorna a distância
# entre eles em quilômetros
def calcularDistancia(lat1, lat2, lon1, lon2):
  raio_terra = 6378.137 # raio da terra em quilômetros
    
  # o primeiro passo é converter as latitudes e longitudes
  # para radianos
  lon1 = math.radians(lon1)
  lon2 = math.radians(lon2)
  lat1 = math.radians(lat1)
  lat2 = math.radians(lat2)
 
  # agora aplicamos a Fórmula de Haversine
  dlon = lon2 - lon1
  dlat = lat2 - lat1
  a = math.pow(math.sin(dlat / 2), 2) + math.cos(lat1) * math.cos(lat2) \
    * math.pow(math.sin(dlon / 2),2)
             
  c = 2 * math.asin(math.sqrt(a))
 
  # e retornamos a distância    
  return(c * raio_terra)

# método principal
def main():
  # vamos solicitar a latitude e longitude das duas localizações
  lat1 = float(input("Informe a primeira latitude: "))
  lon1 = float(input("Informe a primeira longitude: "))
  lat2 = float(input("Informe a segunda latitude: "))
  lon2 = float(input("Informe a segunda longitude: "))

  # vamos calcular a distância entre os dois pontos em Kms
  distancia = calcularDistancia(lat1, lat2, lon1, lon2)
    
  # mostramos o resultado
  print("A distância entre os dois pontos é: {0} kms".format(distancia))

if __name__== "__main__":
  main()