reproject.mosaicking.find_optimal_celestial_wcs(input_data, hdu_in=None, frame=None, auto_rotate=False, projection='TAN', resolution=None, reference=None)[source]

Given one or more images, return an optimal WCS projection object and shape.

This currently only works with 2-d images with celestial WCS.


One or more input data specifications to include in the calculation of the final WCS. This should be an iterable containing one entry for each specification, where a single data specification is one of:

If only one input data needs to be provided, it is also possible to pass it in without including it in an iterable.

hdu_inint or str, optional

If input_data is a FITS file or an HDUList instance, specifies the HDU to use.

framestr or BaseCoordinateFrame

The coordinate system for the final image (defaults to the frame of the first image specified).


Whether to rotate the header to minimize the final image area (if True, requires shapely>=1.6 to be installed).


Three-letter code for the WCS projection.


The resolution of the final image. If not specified, this is the smallest resolution of the input images.


The reference coordinate for the final header. If not specified, this is determined automatically from the input images.


The optimal WCS determined from the input images.


The optimal shape required to cover all the output.