Вопрос по FITS формату.
Есть массив изображений в данном формате. Написал утилиту чтения изображений. Тут проблем нет. Но в данном формате для определения реальных мировых координат из координат текущего пиксела в изображении есть ряд полей.
CRPIX1, CRPIX2 - это координаты опорной точки в пикселах
CRVAL1, CRVAL2 - это её мировые координаты
CD1_1, CD1_2, CD2_1, CD2_2 - это матрица поворота, чтобы преобразовать пиксельные координаты в мировые.
Нахожу разницу в пиксельных координатах относительно интересующей точки и опорной (DX, DY).
Потом получаю координаты
CD1_2*DX+CD2_2*DY - это для DEC - вычисляется верно, проверял сторонней программой FV.
CD1_1*DX+СD2_1*DY - это для RA - вычисляется неверно, проверял так же сторонней программой FV.;
RA всё время сдвинута куда-то.
Комбинировал с расположением коэффициентов CD в данной матрице. Получить аналогичное FV значение для RA не получается.
На сколько я понимаю, нужно матрицу 2х2 которая есть эти CDij умножить на матрицу 1х2 - которое есть разница в пиксельных координатах между текущей и опорной точкой.
Что может быть не так? Ведь если DEC вычисляется верно, то значит и RA должно верно вычисляться.
Так же вопрос, кто-то работал с cfitsio.dll и портировал ли хидерный файл в Delphi?
