座標はフィットを考えてみましょう次のコードを

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 

トリミングされたフレームで座標が異なるのはなぜですか？

---

これは、2つの異なる方法を使用してこれを解決する2つの方法です。

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をコピーします：

prihdr.update(wcs_cropped.to_header()) 

ファイルを保存する前に

一つの方法は、astropy.WCSを使用することです。

cutout_footprintが何であるかわからないので、wcs_croppedを作成するときにスライスインデックスを変更する必要があります。

Cutout2Dという便利な機能があり、デフォルトでWCSが更新されます。

Answer 3

別の答えは、追加のパッケージ必要とするため：astropy.nddata.NDDataに基づいて、特にクラスCCDDataccdprocを：

これは、ファイルを読み取る簡素化：

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