def scan2colfil(x, y, x0, y0, scale, tipo=0):
    """
    Transforma de x/y en proyeccion geoestacionaria a
    En base a 5.2.8.2 de PUG3
    Parameters
    ----------
    x : float, float arr
       coordenada vertical, en radianes
    x : float
       coordenada vertical del primer punto en radianes
    x0 : float
        posición del x[0] en radianes
    y0 : float
        coordenada horizontal del primer punto, en radianes.
        Paralelo al eje terrestre
    scale : float
        tamaño del pixel en radianes
    tipo : TYPE, optional
        tipo de salida. The default is 0 para float, 1 para int.
    Returns
    -------
    col : columna
    fil : fila
    """
    col = (x - x0) / scale
    fil = -(y - y0) / scale
    if tipo == 0:
        return col, fil
    elif tipo == 1:
        return round(col), round(fil)
    else:
        raise TypeError("Type must be 0 (float) or 1 (int)")