SIP-совместимый пакет решений Python для планшетов

Question

Мне нужно выполнить точные преобразования координат пиксель-мир в файлах FITS, которые были изначально созданы с использованием Maxim DL. Максим использует Pinpoint для решения пластины, которое генерирует коэффициенты искажения TRi_j. Они несовместимы с функциями преобразования координат astropy.wcs, которые я предлагал использовать, так как они предполагают коэффициенты искажения SIP.

Поэтому я ищу варианты повторного объединения файлов FITS для генерации коэффициентов SIP.

Пока все, что я нашел, это astrometry.net, но это онлайн-сервис. Я действительно ищу офлайн-листинг (желательно с локальной копией GSC), который я могу выполнять синхронно как часть рабочего процесса моего приложения.

Существуют ли какие-либо связанные с Astropy (или другие) пакеты Python, которые выполняют совместимые с SIP планшеты против GSC?

В качестве альтернативы, есть ли какие-нибудь эквиваленты для wcs.all_pix2world, которые могут использовать коэффициенты искажения TRi_j, чтобы я мог работать с данными Maxim DL?

Большое спасибо

Nigel

Answer 1

Вы можете сделать:

from astropy.io import fits
from astropy.wcs import WCS

hdul = fits.open(fitsfilename)[0]

wcs = WCS(hdul.header)

ax = fig.gca()
ax.scatter([34], [3.2], transform=ax.get_transform('world'))

( На основании этого Q .)

Answer 2

В дополнение к коэффициентам SIP методы astropy.wcs будут работать с коэффициентами искажения TPV . Это означает, что вы можете использовать вывод SCAMP astrometric solver напрямую с astropy.wcs. Если вы хотите преобразовать коэффициенты TPV в форму SIP, вы можете использовать пакет sip_tpv , для которого я являюсь ведущим участником. Я не знаю, как будет упакован пакет Python для SCAMP - я обернул его для конвейера Zwicky Transient Facility , но этот код не является общедоступным.