Transformacion de coordenadas MLP

tcoord_mlp.py

import math as m


def mlp2sirgas(x_in, y_in):
    K   = 0.9993861340849
    W   = 0.0025759302920
    T_x = 299836.7973
    T_y = 6399665.955
    teta= W * m.pi / 180.0

    x = K * (  x_in * m.cos(teta) + y_in * m.sin(teta) ) + T_x
    y = K * ( -x_in * m.sin(teta) + y_in * m.cos(teta) ) + T_y

    return (x, y)


def sirgas2mlp(x_in, x_in):

    K   = 1.0006142429036
    W   = -0.0025759302936
    T_x = -299733.0738
    T_y = -6403610.386
    teta= W * m.pi / 180.0

    x = K * (  x_in * m.cos(teta) + y_in * m.sin(teta) ) + T_x
    y = K * ( -x_in * m.sin(teta) + y_in * m.cos(teta) ) + T_y

    return (x, y)