可以使用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模块中尝试其他提供的方法。