이미지 받기 SDSS DR10의 파일에 맞춤

Question

기본적으로 나는 astropy에 통합 된 피치를 사용하여 here이라는 r-waveband 적합 파일에서 이미지 데이터를 얻고 싶습니다. 내 시도는 다음과 같습니다.

from astropy.io import fits 
import matplotlib.pyplot as plt 

hdulist = fits.open("frame-r-004646-1-0019.fits") 
hdulist.info() 

imgplot = plt.imshow(hdulist[0].data) 
imgplot.set_cmap('binary_r') 

plt.show() 
hdulist.close() 

현재는 jpeg 이미지에서 보이는 밝은 별 대신 작은 점 만 볼 수 있습니다. 또한 ImageHDU에는 이미지가 없습니다 (1D 배열이므로 불가능) 이것이 어떤 이유입니까?

마지막으로,이 이미지에서 내 목표 은하를 (ra, dec)로 식별하려고 시도하고 있습니다. 피트 파일에는 (ra, dec)가 있고 픽셀 당 0.396 arcsec의 픽셀 스케일이 있습니다. 그래서 통합 된 (ra, dec) 코너에 대해 알고 싶습니다. 그래서 이미지 내의 타겟 위치를 외삽 할 수 있습니다!

Answer 1

당신은이 같은 몇 가지 괜찮은 컬러 정상화 할 것입니다 : IM = hdulist가 [0] WCS를 렌더링

을 .DATA

imshow(im,vmin=np.percentile(im,5),vmax=np.percentile(im,95),cmap='gray') 

을, 나는 당신이 경우, aplpy를 사용하는 것이 좋습니다 ds9와 같은 실제 뷰어에서 벗어나지 마십시오.

gc = aplpy.FITSFigure("frame-r-004646-1-0019.fits") 
gc.show_grayscale(vmin=np.percentile(im,5),vmax=np.percentile(im,95))