fitsio

printdiff generates a report containing the differences between any two fits files , HDU's their Headers or Data (independently). The output can be generated into a some stream object which can be used for comparison/verification of data.

Basically this can be implemented as a function having overloads based on

1. File path
2. Headers( object can be obtained using get_

Everyone's favorite DS9/IRAF feature!

This feature would allow some way of defining the image scale according to IRAF's z-scale. This method requires some simple statistics and fitting and should return a minimum and maximum color.

Example Usage

using AstroImages
img = load("image.fits")
cmin, cmax = zscale(img, args...)

plot(img, cmin=cmin, cmax=cmax)

Referenc