если CRPIX1 == NAXIS1/2 и CRPIX2 == NAXIS2/2, то достаточно взять CRVAL1 и CRVAL2.
В других случаях нужно считать.
Понял, буду считать.
Спасибо за питон-скрипт.
Если их нет, провести файл через astrometry.net, получить координаты центра.
тогда нужно поставить соответствующий тэг, чтобы astrometry использовала центральный пиксель в качестве референсного.
Проверил ресолвы astrometry.net, "почти центр" прописан в качестве опорного:
NAXIS1 = 834
NAXIS2 = 624
...
CRPIX1 = 445.763263285 / X reference pixel
CRPIX2 = 321.572100188 / Y reference pixel
То есть считать полюбому придётся.
Интересно, а вывод в stdout от astrometry.net считает координаты центра, как они пишут или этого "почти центра"? Надо проверить.
1: 04:42:11: Run - astrometry: Reading file "/tmp/godzilla.axy"...
1: 04:42:11: Run - astrometry: Only searching for solutions within 1 degrees of RA,Dec (318.262,47.7506)
1: 04:42:11: Run - astrometry: Field 1 did not solve (index index-4210.fits, field objects 21-30).
1: 04:42:11: Run - astrometry: Field 1 did not solve (index index-4209.fits, field objects 21-30).
1: 04:42:12: Run - astrometry: Field 1 did not solve (index index-4208.fits, field objects 21-30).
1: 04:42:12: Run - astrometry: Field 1 did not solve (index index-4207-03.fits, field objects 21-30).
1: 04:42:12: Run - astrometry: log-odds ratio 233.865 (3.68205e+101), 60 match, 0 conflict, 34 distractors, 92 index.
1: 04:42:12: Run - astrometry:
RA,Dec = (318.275,47.7657), pixel scale 3.91511 arcsec/pix.
Я эту последнюю строку регуляркой разбираю, получаю координаты.
Раньше разбирал другую строку:
1: 16:52:24: Run - astrometry: RA,Dec = (318.275,47.7657), pixel scale 3.91511 arcsec/pix.
1: 16:52:24: Run - astrometry: Hit/miss: +++-++++++++--++-+-++++++++--++-++++--+++++-+-+++-++-+-++---+--+-+--+---+-+-+-++-++-+-++-+-+++(best)--++++
1: 16:52:25: Run - astrometry: Field 1: solved with index index-4206-03.fits.
1: 16:52:25: Run - astrometry: Field 1 solved: writing to file /tmp/godzilla.solved to indicate this.
1: 16:52:25: Run - astrometry: Field: /tmp/godzilla.fit
1: 16:52:25: Run - astrometry:
Field center: (RA,Dec) = (318.3, 47.76) deg.
Но сегодня в ночи подумал, что это просто округление цифры выше. Теперь не уверен.