以下是一个 Python 脚本,使用 astropy.coordinatesastropy.units 将 CSV 文件中的 RA 和 Dec 两列数据转换为银河系经度 (Galactic LST),并将转换后的数据保存为新的 CSV 文件。

import csv
from astropy.coordinates import SkyCoord
from astropy import units as u

# 读取 CSV 文件
with open('input.csv', 'r') as f:
    reader = csv.reader(f)
    data = list(reader)

# 转换 RA 和 Dec 数据为 Galactic LST
for row in data:
    ra = float(row[0])
    dec = float(row[1])
    c = SkyCoord(ra=ra*u.deg, dec=dec*u.deg, frame='icrs')
    lst = c.transform_to('galactic').l
    row[0] = str(lst.deg)

# 保存为新的 CSV 文件
with open('output.csv', 'w', newline='') as f:
    writer = csv.writer(f)
    writer.writerows(data)

将上述代码保存为一个 Python 脚本文件,并将输入的 CSV 文件命名为 input.csv,然后运行脚本,即可将转换后的数据保存为新的 CSV 文件 output.csv


原文地址: https://www.cveoy.top/t/topic/o9YF 著作权归作者所有。请勿转载和采集!

免费AI点我,无需注册和登录