Блог

Я пытаюсь обрезать несколько частей из одного изображения и сохранить каждую часть отдельно. Я наткнулся на это, но здесь позиция указана в пикселях, мне нужно указать координаты положения.

 import numpy as np
from astropy.modeling.models import Gaussian2D
from astropy.nddata import Cutout2D
from astropy import units as u
#%matplotlib qt

y, x = np.mgrid[0:500, 0:500]
data = Gaussian2D(1, 50, 100, 10, 5, theta=0.5)(x, y)
import matplotlib.pyplot as plt
#plt.imshow(data, origin='lower')

position = (49.7, 100.1)
size = (41, 51)     # pixels
cutout = Cutout2D(data, position, size)


cutout = Cutout2D(data, position, (41, 51))
plt.imshow(cutout.data, origin='lower')


 

Мне нужно указать позиции из текстового файла, который выглядит так

 RA (Degree)  Dec(Degree)  Size(pixel)
150.2507           2.585      50*50
150.5504           2.585      50*50
150.8508           2.590      50*50
150.9505           2.596      50*50
150.9751           2.599      50*50

 

Заголовок изображения выглядит так

 SIMPLE  =                    T / conforms to FITS standard                      
BITPIX  =                  -64 / array data type                                
NAXIS   =                    3 / number of array dimensions                     
NAXIS1  =                 3307                                                  
NAXIS2  =                 3307                                                  
NAXIS3  =                    1                                                  
COMMENT   FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT   and Astrophysics', volume 376, page 359; bibcode: 2001Aamp;A...376..359H 
LBOUND1 =                -1652 / Pixel origin along axis 1                      
LBOUND2 =                -1652 / Pixel origin along axis 2                      
LBOUND3 =                    1 / Pixel origin along axis 3                      
LABEL   = 'Flux Density'       / Label of the primary array                     
BUNIT   = 'pW      '           / Units of the primary array                     
DATE    = '2017-06-23T01:45:37' / file creation date (YYYY-MM-DDThh:mm:ss UT)   
ORIGIN  = 'Starlink Software'  / Origin of this FITS file                       
HDUCLAS1= 'NDF     '           / Starlink NDF (hierarchical n-dim format)       
HDUCLAS2= 'DATA    '           / Array component subclass                       
HDSTYPE = 'NDF     '           / HDS data type of the component                 
TELESCOP= 'JCMT    '           / Name of Telescope                              
                                                                                
        
 

Пожалуйста, помогите. Спасибо!