def colfil2scan(col, fil, x0, y0, scale):
    """
    Transforma filas/columnas de la imagen a x/y en proyeccion geoestacionaria
    En base a 5.2.8.2 de PUG3
    Parameters
    ----------
    col : int, int arr
        columna
    fil : int, int arr
        fila
    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
    Returns
    -------
    x : float, float arr
       coordenada horizontal, en radianes.
    y : float, float arr
       coordenada vertical, en radianes. Paralelo al eje terrestre
    """
    x = col * scale + x0
    y = -fil * scale + y0
    return x, y