cta_transient_analysis.py 文件源码

def loadFile(input_file):
    df = pd.read_csv(input_file).dropna()

    x = df['stereo:estimated_direction:x']
    y = df['stereo:estimated_direction:y']
    z = df['stereo:estimated_direction:z']

    _, lat, lon = cartesian_to_spherical(
        x.values * u.m, y.values * u.m, z.values * u.m)

    alt = Angle(90 * u.deg - lat).degree
    az = Angle(lon).wrap_at(180 * u.deg).degree
    df['alt'] = alt
    df['az'] = az
    return df