WcsFit

java.lang.Object
- util.PropertyContainer
- - astro.WcsFit

All Implemented Interfaces:

Cloneable, PropertySupplying
```
public class WcsFit
extends PropertyContainer
```
Tries to fit pixels pairs in a fits image to a WcsProjection that minimizes the differences from the projected pixel pairs to the catalog pairs. The matching, i.e. ensuring that pixel pair at index n correponds to the catalog position a t n, must be performed on the outside, calling setMatches(math.VectorG[], math.VectorG[]) throws an IllegalArgumentException, if the two arrays are not of corresponding size.

Nested Class Summary

Nested Classes
Modifier and Type	Class and Description
`static class`	`WcsFit.Linear` We solve for a wcs, reading the pixel pairs and the matched catalog from an ascii file.
`static class`	`WcsFit.Minimization` Three different minimization functions, increasing in robustness.
`private static class`	`WcsFit.Parser`
`static class`	`WcsFit.Polynomial` We solve for a wcs, reading the pixel pairs and the matched catalog from an ascii file.
`static class`	`WcsFit.PositionExchange` To check if the uncertainty in the WCS determination is from the sextractor positions, we take the IRAF positions as input for the wcs fitter.
`static class`	`WcsFit.Residual` This class tries to fit the residuals on command the command line
`static class`	`WcsFit.SigmaMaps` Takes a valid WCS from a file plus a matched file and outputs the sigma maps.
`static class`	`WcsFit.Stepper` We solve for a wcs, reading the pixel pairs and the matched catalog from an ascii file.
`static class`	`WcsFit.Tnx` We solve for a wcs, reading the pixel pairs and the matched catalog from an ascii file.
`static class`	`WcsFit.Variable` We solve for a wcs, reading the pixel pairs and the matched catalog from an ascii file.
`static class`	`WcsFit.Zpn` We solve for a wcs, reading the pixel pairs and the matched catalog from an ascii file.

Field Summary

Fields
Modifier and Type	Field and Description
`private static String`	`AUX`
`private VectorG[]`	`catalog` The model uses these catalog positions to the stars in the image.
`private static double`	`DEFRESTARTSCALE` Default property value.
`private static double`	`DEFTOLERANCE` Default property value.
`static String`	`KEY_RESTARTSCALE` Muliplier in simplex size on restart.
`static String`	`KEY_TOLERANCE` Muliplier in simplex tolerance.
`private WcsModel`	`model` The wcs model we are using.
`private VectorG[]`	`pixel` The model uses these pixel as the identified location of stars.
`private WcsProjection`	`projection` The model is composed with this projection.
`private VectorG`	`solution` After solving, this is the result.
`private static String`	`WIFSIP`
`private static String`	`WIFSIPOLD`

Fields inherited from class util.PropertyContainer
KEY_LISTSEPARATOR, KEY_MAPKEYVALUECHAR, KEY_MAPSEPARATOR

Fields inherited from interface util.PropertySupplying
CONFIG, KEY_CLASS

Constructor Summary

Constructors
Constructor and Description

WcsFit(Map<String,String> prop)
Constructs a WCS regression analyser.

Constructors
Constructor and Description
`WcsFit(Map<String,String> prop)` Constructs a WCS regression analyser.

Method Summary

All Methods Static Methods Instance Methods Concrete Methods
Modifier and Type	Method and Description
`private VectorG`	`amoebaSolve(VectorG start, VectorG len, WcsFit.Minimization how)` We solve with an Amoeba and return the result.
`private static void`	`around(FitsHeader layout, float[] adu, int x0, int y0, int box, float v)`
`int`	`clip(int maxreject, double sigma)` After solving, this method rejects up to the maximum stated number of measures that are outside the specified RMS-range.
`int`	`clip(int maxreject, VectorG rms)` After solving, this method rejects up to the maximum stated number of measures that are outside the specified RMS-range.
`private static MFits`	`createSigma(WcsFit fit)` Save three sigma plots.
`private static MFits`	`createSigma(WcsFit fit, Dimension size, int comp, int blowup)` Save three sigma plots.
`double`	`getRms()` Returns the rms of the fitted model to the model data.
`WcsModel`	`getWcsModel()` Returns the WCS-model used for fitting the coordinate systems.
`VectorG`	`improve(VectorG start)` We try to solve the WCS with Levenberg Marquardt.
`private static VectorG[]`	`linkMatched(File iraf, File sex, double pix)` Helper method that links an IRAF-output match containing CCD positions, ra and dec, but no magnitude to a daofind-like output provided by sextractor.
`VectorG`	`prepareFit(MHeader h)` From a fits header we try to deduce our starting conditions for fitting the WCS model.
`private static VectorG[][]`	`readIraf(File f)` Helper method that reads an ascii-file into an ra/de vectorG array and a vectorG array of pixel x/y.
`private static VectorG[][]`	`readMatched(File f)` Helper method that reads an ascii-file into an ra/de vectorG array and a vectorG array of pixel x/y.
`private static void`	`saveSigma(WcsFit fit)`
`void`	`setMatches(VectorG[] xy, VectorG[] rade)` Sets the matched pixel position and the corresponding catalog entries.
`void`	`setWcsProjection(WcsProjection wcs)` Sets the WCSprojection we use for fitting our model.
`private MImageHDU`	`sigmaImage(Multidimensional sigma, Dimension box, int compress, int blowup)` We take the residuals to the solution and return an sigma-fits from them.
`VectorG`	`solve(VectorG start)` We try to solve for the WCS with a chi-2 model.
`VectorG`	`solve(VectorG start, WcsFit.Minimization how)` We try to solve for the WCS.

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

- Field Detail
  - KEY_RESTARTSCALE
```
public static final String KEY_RESTARTSCALE
```
    Muliplier in simplex size on restart.
    
    See Also:
    
    Constant Field Values
  - KEY_TOLERANCE
```
public static final String KEY_TOLERANCE
```
    Muliplier in simplex tolerance.
    
    See Also:
    
    Constant Field Values
  - DEFRESTARTSCALE
```
private static final double DEFRESTARTSCALE
```
    Default property value.
    
    See Also:
    
    Constant Field Values
  - DEFTOLERANCE
```
private static final double DEFTOLERANCE
```
    Default property value.
    
    See Also:
    
    Constant Field Values
  - model
```
private WcsModel model
```
    The wcs model we are using.
  - projection
```
private WcsProjection projection
```
    The model is composed with this projection.
  - pixel
```
private VectorG[] pixel
```
    The model uses these pixel as the identified location of stars.
  - catalog
```
private VectorG[] catalog
```
    The model uses these catalog positions to the stars in the image.
  - solution
```
private VectorG solution
```
    After solving, this is the result.
  - WIFSIP
```
private static final String WIFSIP
```
    See Also:
    
    Constant Field Values
  - WIFSIPOLD
```
private static final String WIFSIPOLD
```
    See Also:
    
    Constant Field Values
  - AUX
```
private static final String AUX
```
    See Also:
    
    Constant Field Values
- Constructor Detail
  - WcsFit
```
public WcsFit(Map<String,String> prop)
```
    Constructs a WCS regression analyser.
- Method Detail
  - getWcsModel
```
public WcsModel getWcsModel()
```
    Returns the WCS-model used for fitting the coordinate systems.
  - setWcsProjection
```
public void setWcsProjection(WcsProjection wcs)
```
    Sets the WCSprojection we use for fitting our model. If the pixel and catalogs have been set, we readily construct the WcsModel here. It can be retrieved immediately with getWcsModel().
  - setMatches
```
public void setMatches(VectorG[] xy,
                       VectorG[] rade)
```
    Sets the matched pixel position and the corresponding catalog entries. This method does not try to find matches, it is assumed that the pixel and star arrays have been matched already. Consequentially, we throw an IllegalArgumentException, if the length of the arrays is not identical.
  - prepareFit
```
public VectorG prepareFit(MHeader h)
                   throws HeaderException
```
    From a fits header we try to deduce our starting conditions for fitting the WCS model. This is done by passing the header to the estimateParameter method of the projection.
    
    Throws:
    
    HeaderException
  - improve
```
public VectorG improve(VectorG start)
```
    We try to solve the WCS with Levenberg Marquardt.
  - solve
```
public VectorG solve(VectorG start)
```
    We try to solve for the WCS with a chi-2 model.
  - solve
```
public VectorG solve(VectorG start,
                     WcsFit.Minimization how)
```
    We try to solve for the WCS. If we reach convergence, we re-issue with smaller length scale to the Amoeba solver.
  - amoebaSolve
```
private VectorG amoebaSolve(VectorG start,
                            VectorG len,
                            WcsFit.Minimization how)
```
    We solve with an Amoeba and return the result.
  - getRms
```
public double getRms()
```
    Returns the rms of the fitted model to the model data.
  - clip
```
public int clip(int maxreject,
                double sigma)
```
    After solving, this method rejects up to the maximum stated number of measures that are outside the specified RMS-range. If more than the stated number of measures are outside the specified range, those with the highest residuals are discarded first.
    
    Parameters:
    
    maxreject - Maximum number of measures to reject, zero for no bound
    
    rms - The RMS of the current solution is multiplied with this vector and returns the maximum residual per coordinate
    
    Returns:
    
    The number of rejected measures, zero for no rejection.
  - clip
```
public int clip(int maxreject,
                VectorG rms)
```
    After solving, this method rejects up to the maximum stated number of measures that are outside the specified RMS-range. If more than the stated number of measures are outside the specified range, those with the highest residuals are discarded first.
    
    Parameters:
    
    maxreject - Maximum number of measures to reject, zero for no bound
    
    rms - The RMS of the current solution is multiplied with this vector2D and returns the maximum residual per coordinate
    
    Returns:
    
    The number of rejected measures, zero for no rejection.
  - sigmaImage
```
private MImageHDU sigmaImage(Multidimensional sigma,
                             Dimension box,
                             int compress,
                             int blowup)
                      throws BasicFitsException
```
    We take the residuals to the solution and return an sigma-fits from them.
    
    Throws:
    
    BasicFitsException
  - around
```
private static void around(FitsHeader layout,
                           float[] adu,
                           int x0,
                           int y0,
                           int box,
                           float v)
```
  - linkMatched
```
private static final VectorG[] linkMatched(File iraf,
                                           File sex,
                                           double pix)
```
    Helper method that links an IRAF-output match containing CCD positions, ra and dec, but no magnitude to a daofind-like output provided by sextractor. From the pixel matches we construct a 3-D vector that take the magnitudes from the sextractor file, but the coordinates from the IRAF file. Note that the coordinates differ by one (different origin).
    
    Parameters:
    
    pix - maximum allowed pixel tolerance, >1.
  - readMatched
```
private static final VectorG[][] readMatched(File f)
```
    Helper method that reads an ascii-file into an ra/de vectorG array and a vectorG array of pixel x/y. Ra/de are sexagesimal positions. On the third index, the magnitudes (catalog and instrumental) are provided.
  - readIraf
```
private static final VectorG[][] readIraf(File f)
```
    Helper method that reads an ascii-file into an ra/de vectorG array and a vectorG array of pixel x/y. Ra/de are sexagesimal positions in degrees.
  - createSigma
```
private static MFits createSigma(WcsFit fit,
                                 Dimension size,
                                 int comp,
                                 int blowup)
                          throws BasicFitsException
```
    Save three sigma plots.
    
    Throws:
    
    BasicFitsException
  - createSigma
```
private static MFits createSigma(WcsFit fit)
                          throws BasicFitsException
```
    Save three sigma plots.
    
    Throws:
    
    BasicFitsException
  - saveSigma
```
private static void saveSigma(WcsFit fit)
                       throws BasicFitsException,
                              IOException
```
    Throws:
    
    BasicFitsException
    
    IOException

Class WcsFit

Nested Class Summary

Field Summary

Fields inherited from class util.PropertyContainer

Fields inherited from interface util.PropertySupplying

Constructor Summary

Method Summary

Methods inherited from class util.PropertyContainer

Methods inherited from class java.lang.Object

Field Detail

KEY_RESTARTSCALE

KEY_TOLERANCE

DEFRESTARTSCALE

DEFTOLERANCE

model

projection

pixel

catalog

solution

WIFSIP

WIFSIPOLD

AUX

Constructor Detail

WcsFit

Method Detail

getWcsModel

setWcsProjection

setMatches

prepareFit

improve

solve

solve

amoebaSolve

getRms

clip

clip

sigmaImage

around

linkMatched

readMatched

readIraf

createSigma

createSigma

saveSigma