Python 脚本:将 RA 和 Dec 坐标转换为银河系经度 (Galactic LST)
以下是一个 Python 脚本,使用 astropy.coordinates 和 astropy.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 著作权归作者所有。请勿转载和采集!