Я пытаюсь преобразовать некоторые RA и De c Мне нужно в пиксели, чтобы вырезать подходящее изображение, и я продолжаю получать следующую ошибку:
NoConvergence: «WCS.all_world2pix» не удалось сходиться с требуемой точностью. После 2 итераций решение расходится по крайней мере для одной входной точки.
Это код, который я использую:
coordinates_ID = zip(ID, RA, Dec)
pixRA = []
pixDec = []
for i in coordinates_ID:
pixra, pixdec = w.all_world2pix(i[1], i[2], 1)
pixRA.append(pixra)
pixDec.append(pixdec)
test = (pixRA[35], pixDec[35])
size = (200, 200)
cut = Cutout2D(sciim[0].data, test, size).data
plt.imshow(cut, origin='lower', cmap='viridis')
plt.colorbar()
img_file = "/home/myname/test.fits"
img_hdu = fits.PrimaryHDU(cut, header=sciimheader)
img_hdu.writeto(img_file, overwrite=True)
Я должен указать, что arguments_id равен 3 zip таблицы (как три списка) координат объектов в моем изображении и соответствующие имена. Может кто-нибудь определить, что может заставить меня получить ошибку, которую я получаю?