astro.fits

Class MFitsStatistic

java.lang.Object

astro.fits.MFitsStatistic

public class MFitsStatistic
extends Object

Author:

wooz

Nested Class Summary

Nested Classes

Modifier and Type	Class and Description
static class	MFitsStatistic.Scale

Field Summary

Fields

Modifier and Type	Field and Description
static String	BIASAV
private static String	BIASAVREM
static String	BIASRMS
private static String	BIASRMSREM
static String	BIASSEM
private static String	BIASSEMREM
static String	FAVREJ
private static String	FAVREJREM
static String	FCOMBINE
private static String	FCOMBINEREM
static String	FLATAV
private static String	FLATAVREM
static String	FLATRMS
private static String	FLATRMSREM
static String	FLATSEM
private static String	FLATSEMREM
static String	FMAXREJ
private static String	FMAXREJREM
static String	FRMSMAX
private static String	FRMSMAXREM
static String	FXMAX
private static String	FXMAXREM
static String	FXRMS
private static String	FXRMSREM
static String	FYMAX
private static String	FYMAXREM
static String	FYRMS
private static String	FYRMSREM
static String	MMAV
private static String	MMAVREM
static String	MMCOMB
private static String	MMCOMBREM
static String	MMRMS
static String	MMRMSMAX
private static String	MMRMSMAXREM
private static String	MMRMSREM
static String	MMSEM
private static String	MMSEMREM
static String	MMXMAX
private static String	MMXMAXREM
static String	MMYMAX
private static String	MMYMAXREM
static String	RGAINMAX
private static String	RGAINMAXREM
static String	RGAINMIN
private static String	RGAINMINREM
static String	ROWGAIN
private static String	ROWGAINREM
static String	ROWMAX
private static String	ROWMAXREM
static String	ROWMIN
private static String	ROWMINREM
static String	ZAVREJ
private static String	ZAVREJREM
static String	ZCOMBINE
private static String	ZCOMBINEREM
static String	ZMAXREJ
private static String	ZMAXREJREM
static String	ZRMSMAX
private static String	ZRMSMAXREM
static String	ZXMAX
private static String	ZXMAXREM
static String	ZXRMS
private static String	ZXRMSREM
static String	ZYMAX
private static String	ZYMAXREM
static String	ZYRMS
private static String	ZYRMSREM

Constructor Summary

Constructors

Constructor and Description
MFitsStatistic()

Method Summary

Modifier and Type	Method and Description
private static boolean	allPrimary(List<MImageHDU> img) If all image hdus are primary, return true.
static MImageHDU	artificialFloatBias(double level, double readnoise, int... axes) Creates an artificial bias image with the average bias level and a read-noise in ADUs.
static MImageHDU	artificialFloatFlat(double bias, double readnoise, double light, double gain, int... axes) Creates an artificial flat image with the stated bias, read-noise, gain and illumination properties.
static MImageHDU	artificialShortBias(double level, double readnoise, int... axes) Creates an artificial bias image with the average bias level and a read-noise in ADUs.
static MImageHDU	artificialShortFlat(double bias, double readnoise, double light, double gain, int... axes) Creates an artificial flat image with the stated bias, read-noise, gain and illumination properties.
static MImageHDU	avSigClipCombine(List<MImageHDU> flat, MFitsStatistic.Scale mode, boolean combav, double siglo, double sighi, int keep, boolean useav, int nmax) Combines images with the avsigclip algorithm, which is suitable for photon-noise limited data (or no read-out noise data).
static int	fillAduCovSums(ArrayLayout al, double[] fcor1, double[] fcor2, Map<Statistic.Cov,Double> cov, Map<Statistic.Imo,Double> img1, Map<Statistic.Imo,Double> img2, Map<Statistic.ImCov,Double> imcov) Calculates all statistical covariance sums over the two (flat) ADUs.
static int	fillAduCovSums(ArrayLayout al, float[] fcor1, float[] fcor2, Map<Statistic.Cov,Double> cov, Map<Statistic.Imo,Double> img1, Map<Statistic.Imo,Double> img2, Map<Statistic.ImCov,Double> imcov) Calculates all statistical covariance sums over the two (flat) ADUs.
static Vector2D	getFwhmAndIntensityM(double[] adu, double back, double thresh) From the raw adu counts and an background estimation, we calculate peak intensity and FWHM of a gaussian distribution with same average and variance.
static MImageHDU	minMaxCombine(List<MImageHDU> dark, boolean dumphi, boolean dumplo) Combines images by rejecting just the maximum abd/or minimum of the combined ADUs.
private static float[][]	sanityCheck(List<MImageHDU> bias, ArrayLayout al) Takes a list of images and returns a two-dimensional float array, first index is the image number, then element number.
private static FitsHeader	sanitySize(List<MImageHDU> bias) Takes a list of images and returns their layout, if equal
static MImageHDU	zeroCombine(List<MImageHDU> bias, boolean combav, double rn, double gain, double siglo, double sighi, int keep, boolean useav, int nmax) Takes a list of input image hdus and performs a zero-combination on these.

Methods inherited from class java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Field Detail

FLATRMS

public static final String FLATRMS

See Also:

Constant Field Values

FLATSEM

public static final String FLATSEM

See Also:

Constant Field Values

FLATAV

public static final String FLATAV

See Also:

Constant Field Values

FCOMBINE

public static final String FCOMBINE

See Also:

Constant Field Values

FRMSMAX

public static final String FRMSMAX

See Also:

Constant Field Values

FXRMS

public static final String FXRMS

See Also:

Constant Field Values

FYRMS

public static final String FYRMS

See Also:

Constant Field Values

FMAXREJ

public static final String FMAXREJ

See Also:

Constant Field Values

FXMAX

public static final String FXMAX

See Also:

Constant Field Values

FYMAX

public static final String FYMAX

See Also:

Constant Field Values

FAVREJ

public static final String FAVREJ

See Also:

Constant Field Values

ROWGAIN

public static final String ROWGAIN

See Also:

Constant Field Values

RGAINMAX

public static final String RGAINMAX

See Also:

Constant Field Values

ROWMAX

public static final String ROWMAX

See Also:

Constant Field Values

RGAINMIN

public static final String RGAINMIN

See Also:

Constant Field Values

ROWMIN

public static final String ROWMIN

See Also:

Constant Field Values

MMRMS

public static final String MMRMS

See Also:

Constant Field Values

MMSEM

public static final String MMSEM

See Also:

Constant Field Values

MMAV

public static final String MMAV

See Also:

Constant Field Values

MMCOMB

public static final String MMCOMB

See Also:

Constant Field Values

MMRMSMAX

public static final String MMRMSMAX

See Also:

Constant Field Values

MMXMAX

public static final String MMXMAX

See Also:

Constant Field Values

MMYMAX

public static final String MMYMAX

See Also:

Constant Field Values

BIASRMS

public static final String BIASRMS

See Also:

Constant Field Values

BIASSEM

public static final String BIASSEM

See Also:

Constant Field Values

BIASAV

public static final String BIASAV

See Also:

Constant Field Values

ZRMSMAX

public static final String ZRMSMAX

See Also:

Constant Field Values

ZXRMS

public static final String ZXRMS

See Also:

Constant Field Values

ZYRMS

public static final String ZYRMS

See Also:

Constant Field Values

ZCOMBINE

public static final String ZCOMBINE

See Also:

Constant Field Values

ZMAXREJ

public static final String ZMAXREJ

See Also:

Constant Field Values

ZXMAX

public static final String ZXMAX

See Also:

Constant Field Values

ZYMAX

public static final String ZYMAX

See Also:

Constant Field Values

ZAVREJ

public static final String ZAVREJ

See Also:

Constant Field Values

FLATRMSREM

private static final String FLATRMSREM

See Also:

Constant Field Values

FLATSEMREM

private static final String FLATSEMREM

See Also:

Constant Field Values

FLATAVREM

private static final String FLATAVREM

See Also:

Constant Field Values

FRMSMAXREM

private static final String FRMSMAXREM

See Also:

Constant Field Values

FXRMSREM

private static final String FXRMSREM

See Also:

Constant Field Values

FYRMSREM

private static final String FYRMSREM

See Also:

Constant Field Values

FCOMBINEREM

private static final String FCOMBINEREM

See Also:

Constant Field Values

FMAXREJREM

private static final String FMAXREJREM

See Also:

Constant Field Values

FXMAXREM

private static final String FXMAXREM

See Also:

Constant Field Values

FYMAXREM

private static final String FYMAXREM

See Also:

Constant Field Values

FAVREJREM

private static final String FAVREJREM

See Also:

Constant Field Values

ROWGAINREM

private static final String ROWGAINREM

See Also:

Constant Field Values

RGAINMAXREM

private static final String RGAINMAXREM

See Also:

Constant Field Values

ROWMAXREM

private static final String ROWMAXREM

See Also:

Constant Field Values

RGAINMINREM

private static final String RGAINMINREM

See Also:

Constant Field Values

ROWMINREM

private static final String ROWMINREM

See Also:

Constant Field Values

MMRMSREM

private static final String MMRMSREM

See Also:

Constant Field Values

MMSEMREM

private static final String MMSEMREM

See Also:

Constant Field Values

MMAVREM

private static final String MMAVREM

See Also:

Constant Field Values

MMCOMBREM

private static final String MMCOMBREM

See Also:

Constant Field Values

MMRMSMAXREM

private static final String MMRMSMAXREM

See Also:

Constant Field Values

MMXMAXREM

private static final String MMXMAXREM

See Also:

Constant Field Values

MMYMAXREM

private static final String MMYMAXREM

See Also:

Constant Field Values

BIASRMSREM

private static final String BIASRMSREM

See Also:

Constant Field Values

BIASSEMREM

private static final String BIASSEMREM

See Also:

Constant Field Values

BIASAVREM

private static final String BIASAVREM

See Also:

Constant Field Values

ZRMSMAXREM

private static final String ZRMSMAXREM

See Also:

Constant Field Values

ZXRMSREM

private static final String ZXRMSREM

See Also:

Constant Field Values

ZYRMSREM

private static final String ZYRMSREM

See Also:

Constant Field Values

ZCOMBINEREM

private static final String ZCOMBINEREM

See Also:

Constant Field Values

ZMAXREJREM

private static final String ZMAXREJREM

See Also:

Constant Field Values

ZXMAXREM

private static final String ZXMAXREM

See Also:

Constant Field Values

ZYMAXREM

private static final String ZYMAXREM

See Also:

Constant Field Values

ZAVREJREM

private static final String ZAVREJREM

See Also:

Constant Field Values

Constructor Detail

MFitsStatistic

public MFitsStatistic()

Method Detail

getFwhmAndIntensityM

public static Vector2D getFwhmAndIntensityM(double[] adu,
                                            double back,
                                            double thresh)

From the raw adu counts and an background estimation, we calculate peak intensity and FWHM of a gaussian distribution with same average and variance.

Parameters:

adu - array of the image ADUs.

back - background estimation in ADUs.

thresh - If negative, use all pixel, otherwise only those at least thresh above back.

Returns:

At index 0 the peak intensity, at index 1 the fwhm.

artificialFloatBias

public static MImageHDU artificialFloatBias(double level,
                                            double readnoise,
                                            int... axes)
                                     throws BasicFitsException

Creates an artificial bias image with the average bias level and a read-noise in ADUs.

Parameters:

level - The bias level

readnoise - The read-noise in ADUs (no gain)

axes - Dimension of image

Returns:

Throws:

BasicFitsException

artificialShortBias

public static MImageHDU artificialShortBias(double level,
                                            double readnoise,
                                            int... axes)
                                     throws BasicFitsException

Creates an artificial bias image with the average bias level and a read-noise in ADUs.

Parameters:

level - The bias level

readnoise - The read-noise in ADUs (no gain)

axes - Dimension of image

Returns:

Throws:

BasicFitsException

artificialShortFlat

public static MImageHDU artificialShortFlat(double bias,
                                            double readnoise,
                                            double light,
                                            double gain,
                                            int... axes)
                                     throws BasicFitsException

Creates an artificial flat image with the stated bias, read-noise, gain and illumination properties.

Parameters:

bias - The bias level in ADUs

readnoise - in electrons

light - the illumination level in ADUs

gain - The gain of the CCD, ADU*gain=e-

axes - Dimension of image

Returns:

Throws:

BasicFitsException

artificialFloatFlat

public static MImageHDU artificialFloatFlat(double bias,
                                            double readnoise,
                                            double light,
                                            double gain,
                                            int... axes)
                                     throws BasicFitsException

Creates an artificial flat image with the stated bias, read-noise, gain and illumination properties.

Parameters:

bias - The bias level in ADUs

readnoise - in electrons

light - the illumination level in ADUs

gain - The gain of the CCD, ADU*gain=e-

axes - Dimension of image

Returns:

Throws:

BasicFitsException

fillAduCovSums

public static final int fillAduCovSums(ArrayLayout al,
                                       double[] fcor1,
                                       double[] fcor2,
                                       Map<Statistic.Cov,Double> cov,
                                       Map<Statistic.Imo,Double> img1,
                                       Map<Statistic.Imo,Double> img2,
                                       Map<Statistic.ImCov,Double> imcov)

Calculates all statistical covariance sums over the two (flat) ADUs. Mostly used in flat-field context to check for similarities. This sums can than be used to calculate gradients, scaling, etc.

Parameters:

fcor1 - First (flat) field

fcor2 - Second (flat) field.

Returns:

The number of elements summed up. The argument maps are filled with the available numbers if non-null.

fillAduCovSums

public static final int fillAduCovSums(ArrayLayout al,
                                       float[] fcor1,
                                       float[] fcor2,
                                       Map<Statistic.Cov,Double> cov,
                                       Map<Statistic.Imo,Double> img1,
                                       Map<Statistic.Imo,Double> img2,
                                       Map<Statistic.ImCov,Double> imcov)

Calculates all statistical covariance sums over the two (flat) ADUs. Mostly used in flat-field context to check for similarities. This sums can than be used to calculate gradients, scaling, etc.

Parameters:

fcor1 - First (flat) field

fcor2 - Second (flat) field.

Returns:

The number of elements summed up. The argument maps are filled with the available numbers if non-null.

zeroCombine

public static MImageHDU zeroCombine(List<MImageHDU> bias,
                                    boolean combav,
                                    double rn,
                                    double gain,
                                    double siglo,
                                    double sighi,
                                    int keep,
                                    boolean useav,
                                    int nmax)
                             throws BasicFitsException

Takes a list of input image hdus and performs a zero-combination on these. The images should be overscan-corrected and clipped to the useful area. No header entries are retained. On return, an average RMS value is added to the returned header with the key.

Throws:

BasicFitsException

minMaxCombine

public static MImageHDU minMaxCombine(List<MImageHDU> dark,
                                      boolean dumphi,
                                      boolean dumplo)
                               throws BasicFitsException

Combines images by rejecting just the maximum abd/or minimum of the combined ADUs.

Throws:

BasicFitsException

avSigClipCombine

public static MImageHDU avSigClipCombine(List<MImageHDU> flat,
                                         MFitsStatistic.Scale mode,
                                         boolean combav,
                                         double siglo,
                                         double sighi,
                                         int keep,
                                         boolean useav,
                                         int nmax)
                                  throws BasicFitsException

Combines images with the avsigclip algorithm, which is suitable for photon-noise limited data (or no read-out noise data). No sigma-adjustment is done with the scaling, so images should have comparable intensity levels. On calculating the statistic, no sampling is done, thus we might loose some performance.

Throws:

BasicFitsException

allPrimary

private static boolean allPrimary(List<MImageHDU> img)

If all image hdus are primary, return true.

sanitySize

private static FitsHeader sanitySize(List<MImageHDU> bias)
                              throws BasicFitsException

Takes a list of images and returns their layout, if equal

Parameters:

bias -

Throws:

BasicFitsException - If list is null, empty or not all images have identical size.

sanityCheck

private static float[][] sanityCheck(List<MImageHDU> bias,
                                     ArrayLayout al)
                              throws BasicFitsException

Takes a list of images and returns a two-dimensional float array, first index is the image number, then element number.

Parameters:

bias -

Returns:

null, if only a single image is found

Throws:

BasicFitsException - If list is null, empty or not all images have identical size.