Saving mosaic image as a fits file with appropriate header.

Hi.

I wonder if there is a way to save the optimum wcs and the array as a fits image without copying the header of one of the fits file images?

We have data which contains 8 different ccds that we want to mosaic together. We can plot the outputted array in matplotlib which seems to work fine. However when trying to save this as a fits file then the documentation seems to want to do something along the lines of taking one of the headers and using that when writing the fits file. But in this case that would result in chip specific information being used for the entire mosaic image.

Does reproject have a method to write a combined header object in some way which includes the optimum wcs? If not, would it be a feasible addition? And would the current correct method be to create a new header object from scratch?

The issue names no repository files or tests. Start by reading the FITS-writing documentation and the current reprojection output API, then determine whether a combined header is already supported or would need a new entry point. Done means establishing a documented way to save the mosaic array with an appropriate WCS header without reusing chip-specific metadata.

Written by the indexing model from the issue text.