좌표를 맞는이 고려 다음 코드를

Question

파일 (다운로드 test.fits) 자르기 후 보존되지 않습니다

from astropy.io import fits 
from photutils.utils import cutout_footprint 

# Read fits file. 
hdulist = fits.open('test.fits') 
hdu_data = hdulist[0].data 
hdulist.close() 

# Some center and box to crop 
xc, yc, xbox, ybox = 267., 280., 50., 100. 
# Crop image. 
hdu_crop = cutout_footprint(hdu_data, (xc, yc), (ybox, xbox))[0] 
# Add comment to header 
prihdr = hdulist[0].header 
prihdr['COMMENT'] = "= Cropped fits file") 
# Write cropped frame to new fits file. 
fits.writeto('crop.fits', hdu_crop, prihdr) 

을 원본 (왼쪽)와 (오른쪽) 이미지는 다음과 같이립니다 :

프레임 중앙의 별에 대한 (ra, dec) equatorial coordinates은 다음과 같습니다.

Original frame: 12:10:32 +39:24:17 
Cropped frame: 12:12:07 +39:06:50 

자른 프레임에서 좌표가 다른 이유는 무엇입니까?

---

두 가지 방법을 사용하여이를 해결할 수있는 두 가지 방법이 있습니다.

from astropy.io import fits 
from photutils.utils import cutout_footprint 
from astropy.wcs import WCS 
from astropy.nddata.utils import Cutout2D 
import datetime 

# Read fits file. 
hdulist = fits.open('test.fits') 
hdu = hdulist[0].data 
# Header 
hdr = hdulist[0].header 
hdulist.close() 

# Some center and box to crop 
xc, yc, xbox, ybox = 267., 280., 50., 100. 

# First method using cutout_footprint 
# Crop image. 
hdu_crop = cutout_footprint(hdu, (xc, yc), (ybox, xbox))[0] 
# Read original WCS 
wcs = WCS(hdr) 
# Cropped WCS 
wcs_cropped = wcs[yc - ybox // 2:yc + ybox // 2, xc - xbox // 2:xc + xbox // 2] 
# Update WCS in header 
hdr.update(wcs_cropped.to_header()) 
# Add comment to header 
hdr['COMMENT'] = "= Cropped fits file ({}).".format(datetime.date.today()) 
# Write cropped frame to new fits file. 
fits.writeto('crop.fits', hdu_crop, hdr) 

# Second method using Cutout2D 
# Crop image 
hdu_crop = Cutout2D(hdu, (xc, yc), (xbox, ybox), wcs=WCS(hdr)) 
# Cropped WCS 
wcs_cropped = hdu_crop.wcs 
# Update WCS in header 
hdr.update(wcs_cropped.to_header()) 
# Add comment to header 
hdr['COMMENT'] = "= Cropped fits file ({}).".format(datetime.date.today()) 
# Write cropped frame to new fits file. 
fits.writeto('crop.fits', hdu_crop.data, hdr) 

Answer 1

photutils.utils.cutout_footprint은 픽셀 만 잘라내고 픽셀 좌표와 세계 좌표를 변환하는 데 사용되는 WCS는 업데이트하지 않습니다.

대신 astropy.nddata.utils.Cutout2D을 사용하십시오.

Answer 2

이미지를 슬라이스했지만 WCS 정보 (특히 참조 픽셀 값)를 변경하지 않았기 때문에 좌표가 변경되었습니다.

from astropy.wcs import WCS 
wcs = WCS(hdulist[0].header) 
wcs_cropped = wcs[280-50 : 280+50 , 267-25 : 267+25] 

다음 헤더에이 업데이트 WCS 복사 :

한 가지 방법은 astropy.WCS을 사용하는 것입니다

prihdr.update(wcs_cropped.to_header()) 

파일을 저장하기 전에.

cutout_footprint이 무엇인지 잘 모르겠습니다. 따라서 wcs_cropped을 만들 때 조각 인덱스를 변경해야 할 수도 있습니다.

Cutout2D이라는 편의 기능이 WCS으로 기본적으로 업데이트됩니다.

Answer 3

가 추가 패키지를 필요로하기 때문에 또 다른 답 : ccdproc 특히 astropy.nddata.NDData에 따라 클래스 CCDData :

그것은 파일을 읽는 단순화가 :

from ccdproc import CCDData 
ccd = CCDData.read(filename, unit='erg/cm**2/s/AA') 

단위는 헤더 장치 때문에 지정해야합니다을 astropy.units 모듈과 호환되지 않습니다.

CCDData에 대한 중요한 것은 (대부분) 당신이 그 (것)들을, 그들이 속성으로 저장되고 가장 기본적인 작업 보존 (업데이트) 직접 units의 관리, wcs, header 및 masks을 할 필요가 없다는 것입니다 따라서.

ccd_cropped.write('cropped.fits') 

을 그리고 그것은 제대로 업데이트 좌표를 가져야한다 : 예를 들어 슬라이스를 들어 : 당신은 단순히 (다시 저장 등)을 계속 처리 할 수 ​​있도록

xc, yc, xbox, ybox = 267., 280., 50., 100. 
ccd_cropped = ccd[int(yc - ybox // 2) : int(yc + ybox // 2), 
        int(xc - xbox // 2) : int(xc + xbox // 2)] 

이, ccd_cropped 이미 WCS 정보를 업데이트 슬라이스. 슬라이싱 부분은 실제로 astropy.nddata.NDDataRef을 사용하여 가능하며 read 및 write 부분 만 명시 적으로 구현됩니다. ccdproc.CCDData