Crop FITS verileri PyFITS aracılığıyla küp görüntüsü

Question

Sorunum için bir çözüm bulmaya çalıştım, ancak yapamadım. Bir FITS veri küpüm var ve bunu PyFITS ile kesmem gerekiyor. Komutumla yaptığımda, sonunda 2 boyutlu bir FITS resmim olacak! Birinci boyut enerjidir, ikinci ve üçte ikisi sırasıyla boylam ve enlemdir.

Benim komut aşağıdaki gibidir:

#!/usr/bin/env python 
import pyfits 
import os 
import sys 


def CropFitsFile(src, dst, xs, xe, ys, ye): 
    fh = pyfits.open(src) 
    for eng in range(0,2): 
     img = fh[0].data[eng,ys:ye,xs:xe] 
     header = fh[0].header 
     newfh=pyfits.PrimaryHDU(data=img,header=header) 
     if os.path.exists(dst): 
      os.remove(dst) 
     newfh.writeto(dst) 


if __name__ == "__main__": 
    CropFitsFile(
     src=sys.argv[1], 
     dst=sys.argv[2], 
     xs=int(sys.argv[3]), 
     xe=int(sys.argv[4]), 
     ys=int(sys.argv[5]), 
     ye=int(sys.argv[6]) 
     ) 

Answer 1

doğru anlamak bir 3D dizi dilim ama (sadece boyut 1 bile) üçüncü boyutlarını tutmak istiyorsanız.

Bu, Numpy dizileri hakkında bir sorudur. Bir N boyutlu numpy diziniz olduğunda, bir boyut için bir skaler indeksi geçirdiğinizde, indekslediğiniz eksen boyunca dilimlenmiş bir dizi N-1 döndürür. Örneğin: sebebi ne olursa olsun,

>>> arr[:,0,:] 
array([[ 0, 1, 2], 
     [ 9, 10, 11], 
     [18, 19, 20]]) 

İsterseniz bir N-1 boyutlu dizinin yerine N boyutlu bir dizi dönmek:

>>> arr = np.arange(27).reshape(3, 3, 3) 
>>> arr 
array([[[ 0, 1, 2], 
     [ 3, 4, 5], 
     [ 6, 7, 8]], 

     [[ 9, 10, 11], 
     [12, 13, 14], 
     [15, 16, 17]], 

     [[18, 19, 20], 
     [21, 22, 23], 
     [24, 25, 26]]]) 
>>> arr[0] 
array([[0, 1, 2], 
     [3, 4, 5], 
     [6, 7, 8]]) 
>>> arr[1] 
array([[ 9, 10, 11], 
     [12, 13, 14], 
     [15, 16, 17]]) 

Ayrıca gibi farklı bir eksen boyunca dilim En kolay yol, skaler endeks kullanmak yerine 1 numaralı bir dilim dilimini istemektir. Örneğin:

>>> arr[0:1] 
array([[[0, 1, 2], 
     [3, 4, 5], 
     [6, 7, 8]]]) 

Ben bu gibi diğer sorular verdiğim aynı tavsiyeyi vereceğiz: Bu gerçekten veri bir dosya FITS çıktı gerçeği ötesinde PyFITS hakkında bir soru değil. PyFITS, çoğu bilimsel Python kütüphanesi gibi, verileri numpy dizileri olarak döndürür. Bunlar çoğu bilimsel Python uygulamasında sayısal veriler için kullanılan ana veri yapısıdır, bu yüzden bazı temel bilgileri öğrenmek, Python'da veri analizi yapmak için daha iyi veya daha kötü bir önkoşulun önkoşuludur. Eğer MATLAB kullandıysanız, NumPy dizileri MATLAB'deki dizilere benzer. Sen benim kısa öğretici ile başlayabilir, ancak diğerleri (ve çok muhtemelen daha iyi olanlar :) github.com/embray/notebooks/blob/master/numpy.ipynb vardır tam olarak sorunun ne

Answer 2

from astropy.io import fits 

Ccube = fits.open('Cha_binned_ccube.fits', mode='update') 
Ccube.info() 
Ccube[0].shape 
Ccube[0].data = Ccube[0].data[0:3,0:181,0:402] 
Ccube[0].writeto('Cha_binned_ccube_resize.fits')