stella.asiva

Class SkyQualityModel

java.lang.Object

math.AbstractDataModel

stella.asiva.SkyQualityModel

All Implemented Interfaces:

Clipping, DataModel, QuickDesignMatrix

public class SkyQualityModel
extends AbstractDataModel
implements QuickDesignMatrix

A sky quality model has a single dependent variable, the sky quality, as a function of azimuth and altitude. The model is based on real-valued spherical harmonics orthogonal on the upper hemisphere.

Author:

wooz

Nested Class Summary

Nested Classes

Modifier and Type	Class and Description
static class	SkyQualityModel.QuickDM Class to pre-calculate the design matrix for a stated number of parameters and height validity zone.

Nested classes/interfaces inherited from class math.AbstractDataModel

AbstractDataModel.ChiSquare, AbstractDataModel.Lorentzian, AbstractDataModel.Robust

Field Summary

Fields

Modifier and Type	Field and Description
private int	box Averaging box size, only needed for QuickDesignMatrix.
private static String	DESIGN
private static String	EXT
private int	lmax Maximum order of zonal harmonics.
private static String	NAME
private double	xmax Sine of maximum height.
private double	xmin Sine of minimum height.

Constructor Summary

Constructors

Modifier	Constructor and Description
protected	SkyQualityModel()
protected	SkyQualityModel(double hmin, double hmax)
protected	SkyQualityModel(double hmin, double hmax, int b)
protected	SkyQualityModel(VectorG[] azalt, double[] sky, double[] err) Use factory methods to construct a sky quality model.

Method Summary

Modifier and Type	Method and Description
double	evaluateModel(VectorG p, VectorG azalt) The ordering of the parameters must be such that at index 0 we have Y00, then Y1-1, Y10, Y11; Y2-2, etc.
MImageHDU	evaluateSkyImage(VectorG params)
MImageHDU	evaluateSkyImage(VectorG params, AltAzFromFits aa) Using the float data provided by evaluateSkyModel(math.VectorG) we return an ImageHDU ready to be writtn to a fits file.
MImageHDU	evaluateSkyImageResiduals(MImageHDU orgfits, VectorG params, AltAzFromFits aa)
float[]	evaluateSkyModel(VectorG params)
float[]	evaluateSkyModel(VectorG params, AltAzFromFits aa) Given the probably fitted model parameters (the parameters that are the solution from the sky quality fit) we re-construct the original sky-quality fits data in a rectangular area as returned as the fits size in the provided projection (640x480 per default).
float[]	evaluateSkyResiduals(float[] sqadu, VectorG params, AltAzFromFits aa) Given the probably fitted model parameters (the parameters that are the solution from the sky quality fit) we re-construct the original sky-quality fits data in a rectangular area as returned as the fits size in the provided projection (640x480 per default).
double[][]	getDesignMatrix() We try to locate a resource file
double	getHeightMax() Sets the maximum height the fit should operate on, normally 90 degrees.
double	getHeightMin() Minimum height the fit should be applied for.
private String	getId() Returns the maximum l, the height zone and the averaging box size as a string-part for an resource name.
int	getLmax()
String	getName() Returns a resource name that identifies the maximum l, the height zone and the averaging box size.
int	getParameterCount() The number of parameters is given by the maximum degree lmax.
String	getPrecalculatedResourceName() Returns a resource name that could provide a serialized version of the design matrix.
static SkyQualityModel	getSkyQualityModel(float[] adu, LayoutImpl size, int box, float validmin, double hmin, double hmax, AltAzFromFits azalt)
static SkyQualityModel	getSkyQualityModel(MFits avisa, int box, float validmin) Takes a fits file of the sky quality, an averaging box size, and a minimum sky quality value to be considered valid and return a sky model.
static SkyQualityModel	getSkyQualityModel(MImageHDU avisa, int box, float validmin)
static SkyQualityModel	getSkyQualityModel(MImageHDU avisa, int box, float validmin, double hmin, double hmax, AltAzFromFits azalt)
private boolean	inside(VectorG azalt)
protected void	setHeightMax(double hmax) Sets the maximum zenith distance the fit should operate on, normally 90 degrees.
protected void	setHeightMin(double hmin) Sets the minimum height the fit should be applied for.
void	setLmax(int lmax) Sets the degree of the highest spherical harmonics (zonal for the upper hemisphere) used in the model.

Methods inherited from class math.AbstractDataModel

clip, getChi2, getChi2, getChiSquareModel, getChiSquareModel, getComponent, getFixParameters, getFixParameters, getLorentzianModel, getLorentzianSolving, getMeasureCount, getMeasurementErrors, getMeasurementTimes, getMeasures, getModel, getModel, getResiduals, getResiduals, getRms, getRms, getRobustModel, getRobustSolving, getSimplexSolving, getTimeOffset, getTimes, getTotalErrors, setMeasures, setTimes, setTotalErrors

Methods inherited from class java.lang.Object

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

Field Detail

EXT

private static final String EXT

See Also:

Constant Field Values

DESIGN

private static final String DESIGN

See Also:

Constant Field Values

NAME

private static final String NAME

See Also:

Constant Field Values

lmax

private int lmax

Maximum order of zonal harmonics.

xmin

private double xmin

Sine of minimum height.

xmax

private double xmax

Sine of maximum height.

box

private int box

Averaging box size, only needed for QuickDesignMatrix.

Constructor Detail

SkyQualityModel

protected SkyQualityModel(VectorG[] azalt,
                          double[] sky,
                          double[] err)

Use factory methods to construct a sky quality model.

Parameters:

azalt - All azimuth and altitude, where sky quality exists

sky - All sky-quality pixel values.

err - Currently null as no error information provided

SkyQualityModel

protected SkyQualityModel()

SkyQualityModel

protected SkyQualityModel(double hmin,
                          double hmax)

SkyQualityModel

protected SkyQualityModel(double hmin,
                          double hmax,
                          int b)

Method Detail

getPrecalculatedResourceName

public String getPrecalculatedResourceName()

Returns a resource name that could provide a serialized version of the design matrix. Depends on the maximum order l, the height zone and the box size.

getName

public String getName()

Returns a resource name that identifies the maximum l, the height zone and the averaging box size.

getId

private String getId()

Returns the maximum l, the height zone and the averaging box size as a string-part for an resource name.

getDesignMatrix

public double[][] getDesignMatrix()

We try to locate a resource file

Specified by:

getDesignMatrix in interface QuickDesignMatrix

getLmax

public int getLmax()

Returns:

the lmax, the highest degree to which Ylm are used

setLmax

public void setLmax(int lmax)

Sets the degree of the highest spherical harmonics (zonal for the upper hemisphere) used in the model. By convention, all orders from -l to +l are used.

Parameters:

lmax - the the highest degree to which Ylm are used.

getHeightMax

public double getHeightMax()

Sets the maximum height the fit should operate on, normally 90 degrees.

Returns:

the max height in degrees

setHeightMax

protected void setHeightMax(double hmax)

Sets the maximum zenith distance the fit should operate on, normally 90 degrees.

Parameters:

zmax - the zmax to set in degrees

getHeightMin

public double getHeightMin()

Minimum height the fit should be applied for. 15 degrees is reasonable.

Returns:

the zmin in degrees

setHeightMin

protected void setHeightMin(double hmin)

Sets the minimum height the fit should be applied for. 15 degrees is reasonable.

Parameters:

zmin - the zmin to set in degrees

evaluateModel

public double evaluateModel(VectorG p,
                            VectorG azalt)

The ordering of the parameters must be such that at index 0 we have Y00, then Y1-1, Y10, Y11; Y2-2, etc. l^2+2l+1 parameters have to be stated.

Specified by:

evaluateModel in interface DataModel

Returns:

A double with the sky quality model.

getParameterCount

public int getParameterCount()

The number of parameters is given by the maximum degree lmax.

Specified by:

getParameterCount in interface DataModel

evaluateSkyImage

public MImageHDU evaluateSkyImage(VectorG params)
                           throws BasicFitsException

Throws:

BasicFitsException

evaluateSkyImage

public MImageHDU evaluateSkyImage(VectorG params,
                                  AltAzFromFits aa)
                           throws BasicFitsException

Using the float data provided by evaluateSkyModel(math.VectorG) we return an ImageHDU ready to be writtn to a fits file.

Parameters:

params -

aa -

Returns:

Throws:

BasicFitsException

FitsException

evaluateSkyImageResiduals

public MImageHDU evaluateSkyImageResiduals(MImageHDU orgfits,
                                           VectorG params,
                                           AltAzFromFits aa)
                                    throws BasicFitsException

Throws:

BasicFitsException

evaluateSkyModel

public float[] evaluateSkyModel(VectorG params)

evaluateSkyModel

public float[] evaluateSkyModel(VectorG params,
                                AltAzFromFits aa)

Given the probably fitted model parameters (the parameters that are the solution from the sky quality fit) we re-construct the original sky-quality fits data in a rectangular area as returned as the fits size in the provided projection (640x480 per default). Non-defined pixel return zero.

Parameters:

params -

aa -

Returns:

An adu array, first index y, second x (FITS standard)

evaluateSkyResiduals

public float[] evaluateSkyResiduals(float[] sqadu,
                                    VectorG params,
                                    AltAzFromFits aa)

Given the probably fitted model parameters (the parameters that are the solution from the sky quality fit) we re-construct the original sky-quality fits data in a rectangular area as returned as the fits size in the provided projection (640x480 per default). Non-defined pixel return zero.

Parameters:

params -

aa -

Returns:

An adu array, first index y, second x (FITS standard)

inside

private boolean inside(VectorG azalt)

getSkyQualityModel

public static SkyQualityModel getSkyQualityModel(MFits avisa,
                                                 int box,
                                                 float validmin)
                                          throws BasicFitsException,
                                                 IOException

Takes a fits file of the sky quality, an averaging box size, and a minimum sky quality value to be considered valid and return a sky model.

Parameters:

avisa -

box -

validmin -

Returns:

Throws:

FitsException

IOException

BasicFitsException

getSkyQualityModel

public static SkyQualityModel getSkyQualityModel(MImageHDU avisa,
                                                 int box,
                                                 float validmin)
                                          throws BasicFitsException

Throws:

BasicFitsException

getSkyQualityModel

public static SkyQualityModel getSkyQualityModel(MImageHDU avisa,
                                                 int box,
                                                 float validmin,
                                                 double hmin,
                                                 double hmax,
                                                 AltAzFromFits azalt)
                                          throws BasicFitsException

Throws:

BasicFitsException

getSkyQualityModel

public static SkyQualityModel getSkyQualityModel(float[] adu,
                                                 LayoutImpl size,
                                                 int box,
                                                 float validmin,
                                                 double hmin,
                                                 double hmax,
                                                 AltAzFromFits azalt)