stella.adapter

Class AbstractPinholeModel

java.lang.Object

util.PropertyContainer

stella.adapter.AbstractPinholeShift

stella.adapter.AbstractPinholeModel

All Implemented Interfaces:

Cloneable, PinholeShift, Initializable, PropertySupplying

Direct Known Subclasses:

AveragePinhole, ModelDrift

public abstract class AbstractPinholeModel
extends AbstractPinholeShift

Class that determines the pixel shift between pinhole and star by modelling the star plus pinhole. This class comprises methods to quantify the fit quality.

Nested Class Summary

Nested Classes

Modifier and Type	Class and Description
static class	AbstractPinholeModel.AbstractPinholeFits This is the data model that somehow describes the star plus pinhole.

Nested classes/interfaces inherited from class stella.adapter.AbstractPinholeShift

AbstractPinholeShift.Background, AbstractPinholeShift.Evaluate, AbstractPinholeShift.Reconstruct

Field Summary

Fields

Modifier and Type	Field and Description
private Dimension	box The geometry of the data model to derive the 2-dim back.
private static double	DEFMU Pixel size of detector, mu
private static double	DEFPIXELSCALE Pixel size of detector in arcseconds.
private static boolean	DEFUSESIMPLEX Default we are not using simplex.
protected AbstractPinholeModel.AbstractPinholeFits	fit The data model, where the times are the rectangular grid.
private MImageHDU	hdu We also need the original data.
static String	KEY_MU The pixel size in micron .
static String	KEY_PIXELSCALE The pixel size in micron .
static String	KEY_USESIMPLEX If false, we use Levenberg-Marquardt algorithm, otherwise Simplex.
private VectorG	last The last solution.

Fields inherited from class stella.adapter.AbstractPinholeShift

KEY_ADNOISE, KEY_APPLICABLE, KEY_GAIN, KEY_PINHOLE, KEY_SIGMACLIP, KEY_SKYTHICK

Fields inherited from class util.PropertyContainer

KEY_LISTSEPARATOR, KEY_MAPKEYVALUECHAR, KEY_MAPSEPARATOR

Fields inherited from interface util.PropertySupplying

CONFIG, KEY_CLASS

Constructor Summary

Constructors

Modifier	Constructor and Description
protected	AbstractPinholeModel(Map<String,String> prop) Constructs a model-shift.

Method Summary

Modifier and Type	Method and Description
protected abstract AbstractPinholeModel.AbstractPinholeFits	generatePinholeModel(MImageHDU hdu, Rectangle frame, double back, double backerr, Point2D pinhole) This method intializes the pinhole model used in the fitting process during initFiberView(fits.hdu.MFits).
(package private) Dimension	getBox()
static List<String>	getGnuPlotLines(Dimension size, double[] data) Returns a list of lines for a three-dim GNU plot input.
double[]	getMeasures() Returns the flux of the star as it would have been seen without a pinhole as a two-dim model.
double[]	getModel() Returns the flux of the star as it would have been seen without a pinhole as a two-dim model.
double[]	getModel(VectorG with) Returns the flux of the star as it would have been seen without a pinhole as a two-dim model.
double[]	getRelativeStarModel() We divide the residual by the star model flux at all pixels.
double	getRemainingFlux() Gets the measured ADU of the star.
double[]	getResiduals() Returns the residuals of the model as a two-dimensional array representing the original geometry, x in the first index, y in second.
double	getRms() Returns the RMS of the model, zero if no convergence was obtained
VectorG	getSolution() Returns the final solution.
double[]	getStarModel() Returns the flux of the star as it would have been seen without a pinhole as a two-dim model.
abstract Point2D	getStarPosition() Returns the star position.
Point2D	initFiberView(MFits ccd) Returns the center of gravity shift of the image from the pinhole position.
Point2D	solveFiberView() Finally dolves for the pinhole shift.
static double	sumup(double[] data) Sums a model flux.

Methods inherited from class stella.adapter.AbstractPinholeShift

backgroundFrameMean, backgroundFrameMedian, backgroundMean, backgroundMedian, getAduSignal, getApplicableWindow, getFlux, getFwhm, getPinholePosition, getSignalToNoise, init

Methods inherited from class util.PropertyContainer

augment, augment, augment, augment, augment, clone, defaultBoolean, defaultChar, defaultDouble, defaultFloat, defaultInt, defaultLong, defaultObject, defaultObject, defaultProperties, defaultProperty, getAsBoolean, getAsChar, getAsDouble, getAsEnums, getAsFloat, getAsInt, getAsList, getAsLong, getAsMap, getAsMap, getAsObject, getAsObject, getProperties, getProperty, has, initProperties, isNew, parseObject, parseObject, reload, reload, removeProperty, requires, rescanned, setObject, setProperties, setProperty, stringProperties, toString

Methods inherited from class java.lang.Object

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

Methods inherited from interface stella.adapter.PinholeShift

getBackground

Field Detail

KEY_MU

public static final String KEY_MU

The pixel size in micron .

See Also:

Constant Field Values

KEY_PIXELSCALE

public static final String KEY_PIXELSCALE

The pixel size in micron .

See Also:

Constant Field Values

KEY_USESIMPLEX

public static final String KEY_USESIMPLEX

If false, we use Levenberg-Marquardt algorithm, otherwise Simplex.

See Also:

Constant Field Values

DEFMU

private static final double DEFMU

Pixel size of detector, mu

See Also:

Constant Field Values

DEFPIXELSCALE

private static final double DEFPIXELSCALE

Pixel size of detector in arcseconds.

See Also:

Constant Field Values

DEFUSESIMPLEX

private static final boolean DEFUSESIMPLEX

Default we are not using simplex.

See Also:

Constant Field Values

last

private VectorG last

The last solution.

fit

protected AbstractPinholeModel.AbstractPinholeFits fit

The data model, where the times are the rectangular grid.

box

private Dimension box

The geometry of the data model to derive the 2-dim back.

hdu

private MImageHDU hdu

We also need the original data.

Constructor Detail

AbstractPinholeModel

protected AbstractPinholeModel(Map<String,String> prop)

Constructs a model-shift.

Method Detail

generatePinholeModel

protected abstract AbstractPinholeModel.AbstractPinholeFits generatePinholeModel(MImageHDU hdu,
                                                                                 Rectangle frame,
                                                                                 double back,
                                                                                 double backerr,
                                                                                 Point2D pinhole)
                                                                          throws BasicFitsException

This method intializes the pinhole model used in the fitting process during initFiberView(fits.hdu.MFits). We supply the image data, the applicable image rectangle, the background estimate plus its error and the zero-order pinhole position.

Throws:

BasicFitsException

getBox

Dimension getBox()

Returns:

the box

getRms

public double getRms()

Returns the RMS of the model, zero if no convergence was obtained

getSolution

public VectorG getSolution()

Returns the final solution.

getResiduals

public double[] getResiduals()

Returns the residuals of the model as a two-dimensional array representing the original geometry, x in the first index, y in second.

getRemainingFlux

public double getRemainingFlux()

Gets the measured ADU of the star. First, we call the abstract PinholeShift.getBackground(fits.hdu.MImageHDU) method, then we use sigma-clipping to get all pixel that have a signal. The number of pixels with signal and the ADU sum are returned.

Returns:

The number of pixels with signal in the first index, the ADU-sum in the second.

getStarPosition

public abstract Point2D getStarPosition()

Returns the star position.

getStarModel

public double[] getStarModel()

Returns the flux of the star as it would have been seen without a pinhole as a two-dim model.

getMeasures

public double[] getMeasures()

Returns the flux of the star as it would have been seen without a pinhole as a two-dim model.

getModel

public double[] getModel()

Returns the flux of the star as it would have been seen without a pinhole as a two-dim model.

getModel

public double[] getModel(VectorG with)

Returns the flux of the star as it would have been seen without a pinhole as a two-dim model.

sumup

public static double sumup(double[] data)

Sums a model flux.

getRelativeStarModel

public double[] getRelativeStarModel()

We divide the residual by the star model flux at all pixels.

getGnuPlotLines

public static List<String> getGnuPlotLines(Dimension size,
                                           double[] data)

Returns a list of lines for a three-dim GNU plot input. Plot it with

 gnuplot> splot "data.file" matrix with lines
 

To plot multiple files, use

 gnuplot> splot "data.file" matrix with lines, "data2.file" matrix with
 lines
 

Examples can be found here

initFiberView

public Point2D initFiberView(MFits ccd)

Returns the center of gravity shift of the image from the pinhole position. The true center of gravity is scaled with two functions to return true shifts.

Parameters:

ccd - The image of the pinhole entrance

Returns:

The pixel shift from the pinhole center to the star in a first model.

solveFiberView

public Point2D solveFiberView()

Description copied from interface: PinholeShift

Finally dolves for the pinhole shift. Needs a call to PinholeShift.initFiberView(MFits) before.

Returns:

The pixel shift from the pinhole center to the star in the fitted model