Esta respuesta calcula la distancia que se mueve cada barco, sin considerar el tiempo transcurrido.
Para cada barco calcula:
- Posición promedio en base a su lista de coordenadas.
- Posición extrema, la coordenada del punto más alejado de la posición promedio.
- La distancia entre ambos.
Dados estos tres datos, es posible reducir su lista de coordenadas a sólo dos puntos: la posición promedio y la posición extrema.
Para calcular la distancia se ocupa la formula euclidiana ajustada según la latitud del barco.
Los resultados se pueden verificar usando calculadoras en línea como Latitude/Longitude Distance Calculator. Los resultados pueden variar, dado que el sitio ocupa la formula más precisa y general (Haversine).
No me queda claro como quieres la salida, por lo que no produce ningún archivo. Se puede fácilmente producir uno con los barcos inmóviles y sus dos coordenadas, pero no tengo claro que pasa con los que si se mueven. Espero tus comentarios.
Código
import csv
from math import sqrt, cos
miles_per_degree = 59.71 # En el ecuador
def small_distance(lat, lng, lat0, lng0):
"""
Calcula distancia entre dos coordenadas para pequeñas distancias.
Se aplica la formula euclidiana para distancie entre dos puntos,
corregidas por la latitud.
:param lat: Latitud punto partida
:param lng: Longitud punto de partida.
:param lat0: Latitud punto llegada
:param lng0: Longitud punto llegada.
:return: Distancia en millas nauticas.
"""
x = lat - lat0
y = (lng - lng0) * cos(lat0)
return miles_per_degree * sqrt(x * x + y * y)
def max_distance(list_coord):
"""
Calcula la dispersion de un conjunto de coordenadas.
Calcula el promedio de una nube de coordenadas y luego
la distancia de cada punto respecto al promedio.
:param list_coord: Lista de coordenadas (lon, lat)
:return: La maxima distancia en millas naúticas, la
coord promedio y la coord más lejana
"""
lon_average = sum([x[0] for x in list_coord]) / len(list_coord)
lat_average = sum([x[1] for x in list_coord]) / len(list_coord)
print(" coord average =", (lon_average, lat_average))
max_dist = 0
max_coord = None
for (lon, lat) in list_coord:
distance = small_distance(lat, lon, lat_average, lon_average)
if distance >= max_dist:
max_dist = distance
max_coord = (lon, lat)
print(" coord fartest =", max_coord)
return max_dist, (lon_average, lat_average), max_coord
def producir_resultados(mmsi, list_coord):
"""
Calcula e imprime para un barco en particular.
:param mmsi: El barco.
:param list_coord: La lista de coordenadas
:return: Maxima distancia recorrida, coord promedio y coord más lejana
"""
print("MMSI", mmsi, "(", len(list_coord), " puntos)")
if len(list_coord) > 1:
max_dist, average, fartest = max_distance(list_coord)
else:
max_dist = 0.0
average = fartest = list_coord[0]
print(" Maxima distancia =", max_dist)
print()
return max_dist, average, fartest
def calcula_distancia(csv_file):
"""
Calcula la distancia para los barcos en un archivos.
:param csv_file: El archivo con coord.
"""
with open(csv_file) as csv_file:
csv_reader = csv.reader(csv_file, delimiter=',')
mmsi_act = None
list_coord = []
header = next(csv_reader, None)
for mmsi, datetime, lon, lat in csv_reader:
coord = (float(lon), float(lat))
if mmsi_act is None:
mmsi_act = mmsi
if mmsi_act != mmsi:
producir_resultados(mmsi_act, list_coord)
mmsi_act = mmsi
list_coord = []
list_coord.append(coord)
producir_resultados(mmsi_act, list_coord)
#
# Revisar distancias usando
# https://www.meridianoutpost.com/resources/etools/calculators/calculator-latitude-longitude-distance.php?
#
# Ejemplo 1: Valores correctos son 49.79 miles / 80.13 kilometers / 43.24 nautical miles
#
print("Test para MMSI 209698000")
mmsi_209698000 = [(-1.89616501331329, 36.9659996032715),
(-0.959831655025482, 37.571834564209),
(-0.959831655025482, 37.571834564209),
(-0.959831655025482, 37.571834564209)]
max_dist, pos_average, pos_fartest = max_distance(mmsi_209698000)
print(" Distancia = ", max_dist, " millas naúticas.")
#
# Ejemplo 2: El barco no se mueve de su posición.
#
print("Test para MMSI 100000000")
mmsi_100000000 = [(-5.35431146621704, 36.1574516296387),
(-5.35431146621704, 36.1574516296387),
(-5.35431146621704, 36.1574516296387),
(-5.35431146621704, 36.1574516296387),
(-5.35431146621704, 36.1574516296387)]
max_dist, pos_average, pos_fartest = max_distance(mmsi_100000000)
print(" Distancia = ", max_dist, " millas naúticas.")
#
# Ejemplo 3: Pequeña distancia
# Valores correctos: 0.00 miles / 0.00 kilometers / 0.00 nautical miles
#
print("Test para MMSI 205834010")
mmsi_205834010 = [(-2.60475993156433 , 36.7570495605469),
(-2.60475826263428 , 36.7570571899414),
(-2.6047465801239 , 36.7570381164551),
(-2.60472822189331 , 36.7570648193359),
(-2.60472321510315 , 36.7570724487305),
(-2.60469508171082 , 36.7570648193359),
(-2.60473656654358 , 36.7570724487305),
(-2.60474991798401 , 36.7570610046387),
(-2.60470509529114 , 36.7570533752441)]
max_dist, pos_average, pos_fartest = max_distance(mmsi_205834010)
print(" Distancia = ", max_dist, " millas naúticas.")
print()
print("---- Proceso archivo ----")
calcula_distancia('barcos.csv')
Produce
Test para MMSI 209698000
coord average = (-1.193914994597434, 37.420375823974624)
coord fartest = (-1.89616501331329, 36.9659996032715)
Distancia = 48.59234211469109 millas naúticas.
Test para MMSI 100000000
coord average = (-5.35431146621704, 36.1574516296387)
coord fartest = (-5.35431146621704, 36.1574516296387)
Distancia = 0.0 millas naúticas.
Test para MMSI 205834010
coord average = (-2.6047336525387244, 36.757059309217674)
coord fartest = (-2.60469508171082, 36.7570648193359)
Distancia = 0.0013938850605571295 millas naúticas.
---- Proceso archivo ----
MMSI 17620 ( 2418 puntos)
coord average = (-6.180518489518966, 36.497049667699166)
coord fartest = (-6.18065643310547, 36.4969139099121)
Maxima distancia = 0.008632633913399173
MMSI 20233 ( 40 puntos)
coord average = (-6.051580440998078, 35.45204601287842)
coord fartest = (-6.05156850814819, 35.4520416259766)
Maxima distancia = 0.0005171909496390382
MMSI 2424103 ( 1 puntos)
Maxima distancia = 0.0
MMSI 2424105 ( 227 puntos)
coord average = (-4.661880942693366, 35.21154432254741)
coord fartest = (-4.66185998916626, 35.2115249633789)
Maxima distancia = 0.0015238578744766443
...
MMSI 620069000 ( 8 puntos)
coord average = (-0.30466312170028675, 35.85506343841553)
coord fartest = (-0.304684996604919, 35.8550300598145)
Maxima distancia = 0.002023965931953086
MMSI 636012425 ( 1 puntos)
Maxima distancia = 0.0
MMSI 671222100 ( 62 puntos)
coord average = (-2.466488011421696, 36.83163489064863)
coord fartest = (-2.46646165847778, 36.8316688537598)
Maxima distancia = 0.0022689686308574636
mmsi
es la identificación del barco?