Convertir coordenadas geográficas a UTM

Método para convertir coordenadas geográficas a UTM
Utilizando las fórmulas de Coticchia-Surace.
Ejemplo : Santander, Cantabria, España. Calle Rahela 47
Coordenadas geográficas.
Longitud =3˚48’07.7438’’W
Latitud =43˚28’18.2670’’N
Procedimiento
1.-Convertir a decimales.
Longitud 𝜆= 3+ (48/60)+ (0.734399/3600)= 3.802039972
Latitud 𝜑= 43+ (28/60)+ (18.2670/3600)= 43.4884075

2.- Convertirde decimales a radianes.
Si la longitud está referida al Este del meridiano de Greenwich entonces es positiva.
Si la longitud está referida al Oeste del meridiano de Greenwich entonces es negativa.Si la Latitud está referida al Norte del Ecuador entonces es positiva.
Si la Latitud está referida al Sur del Ecuador entonces es negativa.

Longitud 𝜆= = = -0.0663581158
Latitud 𝜑= = =0.759015897

3.- Cálculo de Huso

Huso= +31
Huso== 30.36 Se truca la parte entera y se toma el 30
Ya que la parte del meridiano central dejaría la parte del huso centrada a laizquierda con coordenadas X negativas, para evitar eso, se suma a todas las coordenadas X la cantidad de 500,000, esto hace que no existan valores negativos para las coordenadas X, de igual forma se realizapara los valores Y, cuyo origen es el Ecuador, si el ecuador es el origen de las Y, toda la parte situada al sur del mismo tendría coordenadas negativas ,para que no ocurra eso , se suma el valor10.000.000 a los valores de Y, pero sólo en el caso de que se trate de coordenadas pertenecientes al hemisferio sur; si las coordenadas pertenecen al hemisferio norte, no se tocan los valores Y.

4.-Cálculo sobre la geometría del elipsoide
Semieje mayor 6378388.0
Semieje menor 6356911.94613
Excentricidad= = =0.08199189
Segunda excentricidad = = =
Cuadrado de la segunda excentricidad==0.00676817

5.- Cálculo de radio polar de curvatura y aplanamiento.
Radio polar de curvatura
= =63999936.608
Aplanamiento
𝛼= = =0.003367003

6.- Distancia angular entre la longitud del...