从一个FITS文件复制空间中的标题从一个FITS文件到新创建的FITS文件

Question

我有一个程序处理现有的FITS文件（xbulge-w1.fits），并将处理后的图像保存到新的FITS（w1_resampled.fits）。我想将原始标题复制到新的标题中，以便它们处于相同的坐标（即星系）。
我想用下面的代码来做到这一点：

# Open the FITS files as input image and mask 
# Process the images 
# Rescale image to galactic coordinates and display 
# Plot and save median filtered images as png, rescaled to galactic coordinates 
# Save as FITS files and close 
# Edit FITS headers to recenter images at galactic center 
header = fits.getdata('xbulge-w1.fits', header=True) 

header['COMMENT'] = 'Resampled with median filtered pixels' 
header['IMAGEW'] = 877 
header['IMAGEH'] = 901 
header['WCSAXES'] = (2, 'Number of coordinate axes') 
header['CRPIX1'] = 438.5 
header['CRPIX2'] = 450.5             
header['PC1_1'] = (-0.0333333333333, 'Coordinate transformation matrix element') 
header['PC2_2'] = (0.0333333333333, 'Coordinate transformation matrix element') 
header['CDELT1'] = (1., '[deg] Coordinate increment at reference point') 
header['CDELT2'] = (1., '[deg] Coordinate increment at reference point') 
header['CUNIT1'] = ('deg  ', 'Units of coordinate increment and value') 
header['CUNIT2'] = ('deg  ', 'Units of coordinate increment and value') 
header['CTYPE1'] = 'GLON-AIT'     
header['CTYPE2'] = 'GLAT-AIT'                
header['CRVAL1'] = (0., '[deg] Coordinate value at reference point') 
header['CRVAL2'] = (0., '[deg] Coordinate value at reference point') 
header['LONPOLE'] = (0., '[deg] Native longitude of celestial pole') 
header['LATPOLE'] = (90., '[deg] Native latitude of celestial pole') 
header['RADESYS'] = ('ICRS ', 'Equatorial coordinate system') 

fits.writeto('w1_resampled.fits', header, overwrite=True) 
fits.writeto('w2_resampled.fits', header, overwrite=True) 

hdulist.close() 
hdulist2.close() 
hdulist3.close() 

第5条评论列出工作职能的程序，它只是我遇到问题的标题。当我在DS9打开w1_resampled.fits，头部包含

SIMPLE =     T/conforms to FITS standard      
BITPIX =     -32/array data type         
NAXIS =     2/number of array dimensions      
NAXIS1 =     877             
NAXIS2 =     901             
EXTEND =     T             
END 

，而不是从xbulge-w1.fits头复制的数据。
如何将数据从一个标题复制到另一个标题？

Answer 1

哦，我明白了。我太过于复杂了。这就是我需要做的：在[文件]

# Edit FITS headers to recenter images at galactic center 
w1_resampled_header = w1header 
w2_resampled_header = w2header 

w1_resampled_header['CRPIX1'] = w1header['CRPIX1'] 
w1_resampled_header['CRPIX2'] = w1header['CRPIX2'] 
w2_resampled_header['CRPIX1'] = w2header['CRPIX1'] 
w2_resampled_header['CRPIX2'] = w2header['CRPIX2']             

# Save as FITS files and close 
fits.writeto('w1_resampled.fits', 
      w1_resampled, 
      w1_resampled_header, 
      overwrite = True) 
fits.writeto('w2_resampled.fits', 
      w2_resampled, 
      w2_resampled_header, 
      overwrite = True)