Cálculo de Coordenadas Geodésicas, Distancias y Tiempos Astronómicos en MATLAB
Cálculo de Coordenadas Geodésicas, Distancias y Tiempos Astronómicos
El siguiente código en MATLAB implementa un conjunto de cálculos de geodesia y astronomía de posición: distancias entre estaciones, transformaciones de coordenadas entre marcos de referencia y cálculos de tiempos astronómicos.
1. Longitud de arco paralelo entre estaciones
Se convierte la latitud de la primera estación a radianes y se calcula el radio de curvatura en el primer vertical N1. La diferencia de longitudes se obtiene como lonT = lon1 - lon2.
lat1 = deg2rad(dms2degrees(latY));
lon1 y2 ... lonT = lon1 - lon2;
N1 = a / sqrt(1 - e2*sin(lat1)^2);
Distancia del paralelo: pdist1 = N1*cos(lat1).*lonT;
DifLongP = abs(pdist1 - pdist2);