python - 对齐天文图像

Question

我在不同夜晚拍摄的不同滤镜（U、B、V、R、I）中的天文来源的 .fits 文件（图像）很少。现在我已经对这些图像进行了天体测量并获得了 wcs 信息。如何对齐这些图像，以便我可以在每个过滤器中添加多个图像并增加信噪比。它应该考虑平移和旋转。一些 python 包或代码或任何可用于 IRAF 的任务会更好。

我已经尝试过 IRAF 中可用的 geomap 和 geotran 任务。