可以使用Python中的Astropy库中的photutils模块来找到.fits图像中的星的质心。下面是一个简单的示例代码,演示如何使用photutils模块来找到星的质心:
from astropy.io import fits
from astropy.stats import sigma_clipped_stats
from photutils import centroid_com
# Load the FITS image
image_file = 'star_image.fits'
with fits.open(image_file) as hdul:
image_data = hdul[0].data
# Get the background statistics of the image
mean, median, std = sigma_clipped_stats(image_data)
# Find the centroid of the star
x, y = centroid_com(image_data, mask=None, centroid_func=centroid_com)
print("Centroid of the star: ({}, {})".format(x, y))
这段代码会读取一个名为“star_image.fits”的.fits图像,然后使用sigma_clipped_stats函数来找到图像的平均值、中位数和标准差。之后会使用centroid_com函数来找到星的质心,并使用x和y变量存储结果。最后,代码会打印出星的质心的坐标。
值得注意的是,上面的例子中使用的是centroid_com函数,它通过计算图像的中心位置来找到星的质心。如果你想要使用其他方法来找到质心,你可以在photutils模块中尝试其他提供的方法。