aboutsummaryrefslogtreecommitdiff
path: root/pkg/images/immatch/src/linmatch/linmatch.h
diff options
context:
space:
mode:
Diffstat (limited to 'pkg/images/immatch/src/linmatch/linmatch.h')
-rw-r--r--pkg/images/immatch/src/linmatch/linmatch.h298
1 files changed, 298 insertions, 0 deletions
diff --git a/pkg/images/immatch/src/linmatch/linmatch.h b/pkg/images/immatch/src/linmatch/linmatch.h
new file mode 100644
index 00000000..0f776901
--- /dev/null
+++ b/pkg/images/immatch/src/linmatch/linmatch.h
@@ -0,0 +1,298 @@
+# Header file for LINSCALE
+
+define LEN_LSSTRUCT (70 + 12 * SZ_FNAME + 12)
+
+# Quantities that define the current region and the number of regions
+
+define LS_CNREGION Memi[$1] # the current region
+define LS_NREGIONS Memi[$1+1] # the number of regions
+define LS_MAXNREGIONS Memi[$1+2] # the maximum number of regions
+
+# Quantities that are dependent on the number of regions
+
+define LS_RC1 Memi[$1+3] # pointers to first columns of regions
+define LS_RC2 Memi[$1+4] # pointers to last columns of regions
+define LS_RL1 Memi[$1+5] # pointer to first lines of regions
+define LS_RL2 Memi[$1+6] # pointers to last lines of regions
+define LS_RXSTEP Memi[$1+7] # pointers to the x step sizes
+define LS_RYSTEP Memi[$1+8] # pointers to the y step sizes
+define LS_XSHIFT Memr[P2R($1+9)] # the x shift from image to reference
+define LS_YSHIFT Memr[P2R($1+10)] # the y shift from image to reference
+define LS_SXSHIFT Memr[P2R($1+11)] # the x shift from image to reference
+define LS_SYSHIFT Memr[P2R($1+12)] # the y shift from image to reference
+
+define LS_RBUF Memi[$1+14] # pointer to the reference image data
+define LS_RGAIN Memr[P2R($1+15)] # the reference image gain
+define LS_RREADNOISE Memr[P2R($1+16)] # the reference image readout noise
+define LS_RMEAN Memi[$1+17] # pointers to means of ref regions
+define LS_RMEDIAN Memi[$1+18] # pointers to medians of ref regions
+define LS_RMODE Memi[$1+19] # pointers to modes of ref regions
+define LS_RSIGMA Memi[$1+20] # pointers to stdevs of ref regions
+define LS_RSKY Memi[$1+21] # pointers to sky values of ref regions
+define LS_RSKYERR Memi[$1+22] # pointers to sky errors of ref regions
+define LS_RMAG Memi[$1+23] # pointers to magnitudes of ref regions
+define LS_RMAGERR Memi[$1+24] # pointers to mag errors of ref regions
+define LS_RNPTS Memi[$1+25] # pointers to npts of ref regions
+
+define LS_IBUF Memi[$1+27] # pointer to the input image data
+define LS_IGAIN Memr[P2R($1+28)] # the input image gain
+define LS_IREADNOISE Memr[P2R($1+29)] # the input image readout noise
+define LS_IMEAN Memi[$1+30] # pointers to means of image regions
+define LS_IMEDIAN Memi[$1+31] # pointers to medians of image regions
+define LS_IMODE Memi[$1+32] # pointers to modes of image regions
+define LS_ISIGMA Memi[$1+33] # pointers to stdevs of image regions
+define LS_ISKY Memi[$1+34] # pointers to sky values of image regions
+define LS_ISKYERR Memi[$1+35] # pointers to sky errors of image regions
+define LS_IMAG Memi[$1+36] # pointers to magnitudes of image regions
+define LS_IMAGERR Memi[$1+37] # pointers to mag errors of image regions
+define LS_INPTS Memi[$1+38] # pointers to npts of image regions
+
+define LS_RBSCALE Memi[$1+39] # pointers to bscales of regions
+define LS_RBSCALEERR Memi[$1+40] # pointers to bscale errors of regions
+define LS_RBZERO Memi[$1+41] # pointers to bzero errors of regions
+define LS_RBZEROERR Memi[$1+42] # pointers to bzero errors of regions
+define LS_RDELETE Memi[$1+43] # pointer to the delete array
+define LS_RCHI Memi[$1+44] # pointer to the resid array
+
+# Quantities that affect the fitting algorithms
+
+define LS_BSALGORITHM Memi[$1+45] # bscale fitting algorithm
+define LS_BZALGORITHM Memi[$1+46] # bzero fitting algorithm
+define LS_CBZERO Memr[P2R($1+47)] # constant bzero
+define LS_CBSCALE Memr[P2R($1+48)] # constant bscale
+define LS_DNX Memi[$1+49] # x width of data region to extract
+define LS_DNY Memi[$1+50] # y width of data region to extract
+#define LS_PNX Memi[$1+51] # x width of photometry region
+#define LS_PNY Memi[$1+52] # y widht of photometry region
+define LS_DATAMIN Memr[P2R($1+51)] # the minimum good data value
+define LS_DATAMAX Memr[P2R($1+52)] # the maximum good data value
+define LS_MAXITER Memi[$1+53] # maximum number of iterations
+define LS_NREJECT Memi[$1+54] # maximum number of rejections cycles
+define LS_LOREJECT Memr[P2R($1+55)] # low-side sigma rejection criterion
+define LS_HIREJECT Memr[P2R($1+56)] # high-side sigma rejection criterion
+define LS_GAIN Memr[P2R($1+57)] # the constant gain value in e-/adu
+define LS_READNOISE Memr[P2R($1+58)] # the constant readout noise value in e-
+
+# Quantities that define the answers
+
+define LS_TBSCALE Memr[P2R($1+59)] # bzero value
+define LS_TBSCALEERR Memr[P2R($1+60)] # bscale error estimate
+define LS_TBZERO Memr[P2R($1+61)] # bzero value
+define LS_TBZEROERR Memr[P2R($1+62)] # bzero error estimate
+
+# String quantities
+
+define LS_BSSTRING Memc[P2C($1+65)] # bscale string
+define LS_BZSTRING Memc[P2C($1+65+SZ_FNAME+1)] # bzero string
+define LS_CCDGAIN Memc[P2C($1+65+2*SZ_FNAME+2)] # gain keyword
+define LS_CCDREAD Memc[P2C($1+65+3*SZ_FNAME+3)] # readout noise keyword
+define LS_IMAGE Memc[P2C($1+65+4*SZ_FNAME+4)] # input image
+define LS_REFIMAGE Memc[P2C($1+65+5*SZ_FNAME+5)] # reference image
+define LS_REGIONS Memc[P2C($1+65+6*SZ_FNAME+6)] # regions list
+define LS_DATABASE Memc[P2C($1+65+7*SZ_FNAME+7)] # database file
+define LS_OUTIMAGE Memc[P2C($1+65+8*SZ_FNAME+8)] # output image
+define LS_SHIFTSFILE Memc[P2C($1+65+9*SZ_FNAME+9)] # shifts file
+define LS_PHOTFILE Memc[P2C($1+65+10*SZ_FNAME+10)] # shifts file
+define LS_RECORD Memc[P2C($1+65+11*SZ_FNAME+11)] # the record name
+
+
+# Define the bzero and bscale fitting algorithms
+
+define LS_MEAN 1
+define LS_MEDIAN 2
+define LS_MODE 3
+define LS_FIT 4
+define LS_PHOTOMETRY 5
+define LS_FILE 6
+define LS_NUMBER 7
+
+define LS_SCALING "|mean|median|mode|fit|photometry|file|"
+
+# Define the parameters
+
+define CNREGION 1
+define NREGIONS 2
+define MAXNREGIONS 3
+
+define RC1 4
+define RC2 5
+define RL1 6
+define RL2 7
+define RXSTEP 8
+define RYSTEP 9
+define XSHIFT 10
+define YSHIFT 11
+define SXSHIFT 12
+define SYSHIFT 13
+
+define RBUF 14
+define RGAIN 15
+define RREADNOISE 16
+define RMEAN 17
+define RMEDIAN 18
+define RMODE 19
+define RSIGMA 20
+define RSKY 21
+define RSKYERR 22
+define RMAG 23
+define RMAGERR 24
+define RNPTS 25
+
+define IBUF 26
+define IGAIN 27
+define IREADNOISE 28
+define IMEAN 29
+define IMEDIAN 30
+define IMODE 31
+define ISIGMA 32
+define ISKY 33
+define ISKYERR 34
+define IMAG 35
+define IMAGERR 36
+define INPTS 37
+
+define RBSCALE 38
+define RBSCALEERR 39
+define RBZERO 40
+define RBZEROERR 41
+define RDELETE 42
+define RCHI 43
+
+define BZALGORITHM 44
+define BSALGORITHM 45
+define CBZERO 46
+define CBSCALE 47
+define DNX 48
+define DNY 49
+#define PNX 50
+#define PNY 51
+define DATAMIN 50
+define DATAMAX 51
+define MAXITER 52
+
+define NREJECT 53
+define LOREJECT 54
+define HIREJECT 55
+define GAIN 56
+define READNOISE 57
+
+define TBZERO 58
+define TBZEROERR 59
+define TBSCALE 60
+define TBSCALEERR 61
+
+define BSSTRING 62
+define BZSTRING 63
+define CCDGAIN 64
+define CCDREAD 65
+
+define IMAGE 66
+define REFIMAGE 67
+define REGIONS 68
+define DATABASE 69
+define OUTIMAGE 70
+define RECORD 71
+define SHIFTSFILE 72
+define PHOTFILE 73
+
+# Set some default values
+
+define DEF_MAXNREGIONS 100
+define DEF_BZALGORITHM LS_FIT
+define DEF_BSALGORITHM LS_FIT
+define DEF_CBZERO 0.0
+define DEF_CBSCALE 1.0
+define DEF_DNX 31
+define DEF_DNY 31
+define DEF_MAXITER 10
+define DEF_DATAMIN INDEFR
+define DEF_DATAMAX INDEFR
+define DEF_NREJECT 0
+define DEF_LOREJECT INDEFR
+define DEF_HIREJECT INDEFR
+define DEF_GAIN INDEFR
+define DEF_READNOISE INDEFR
+
+# The mode computation parameters.
+
+define LMODE_NMIN 10
+define LMODE_ZRANGE 1.0
+define LMODE_ZBIN 0.1
+define LMODE_ZSTEP 0.01
+define LMODE_HWIDTH 3.0
+
+# The default plot types.
+
+define LS_MMHIST 1
+define LS_MMFIT 2
+define LS_MMRESID 3
+define LS_RIFIT 4
+define LS_RIRESID 5
+define LS_BSZFIT 6
+define LS_BSZRESID 7
+define LS_MAGSKYFIT 8
+define LS_MAGSKYRESID 9
+
+# The bad point deletions code.
+
+define LS_NO 0
+define LS_BADREGION 1
+define LS_BADSIGMA 2
+define LS_DELETED 3
+
+# Commands
+
+define LSCMDS "|input|reference|regions|lintransform|output|photfile|\
+shifts|records|xshift|yshift|dnx|dny|maxnregions|datamin|datamax|\
+maxiter|nreject|loreject|hireject|gain|readnoise|show|markcoords|marksections|"
+
+define LSCMD_IMAGE 1
+define LSCMD_REFIMAGE 2
+define LSCMD_REGIONS 3
+define LSCMD_DATABASE 4
+define LSCMD_OUTIMAGE 5
+define LSCMD_PHOTFILE 6
+define LSCMD_SHIFTSFILE 7
+define LSCMD_RECORD 8
+define LSCMD_XSHIFT 9
+define LSCMD_YSHIFT 10
+define LSCMD_DNX 11
+define LSCMD_DNY 12
+define LSCMD_MAXNREGIONS 13
+define LSCMD_DATAMIN 14
+define LSCMD_DATAMAX 15
+define LSCMD_MAXITER 16
+define LSCMD_NREJECT 17
+define LSCMD_LOREJECT 18
+define LSCMD_HIREJECT 19
+define LSCMD_GAIN 20
+define LSCMD_READNOISE 21
+define LSCMD_SHOW 22
+define LSCMD_MARKCOORDS 23
+define LSCMD_MARKSECTIONS 24
+
+# Keywords
+
+define KY_REFIMAGE "reference"
+define KY_IMAGE "input"
+define KY_REGIONS "regions"
+define KY_DATABASE "lintransform"
+define KY_OUTIMAGE "output"
+define KY_PHOTFILE "photfile"
+define KY_SHIFTSFILE "shifts"
+define KY_RECORD "records"
+define KY_XSHIFT "xshift"
+define KY_YSHIFT "yshift"
+define KY_DNX "dnx"
+define KY_DNY "dny"
+define KY_MAXNREGIONS "maxnregions"
+define KY_DATAMIN "datamin"
+define KY_DATAMAX "datamax"
+define KY_MAXITER "maxiter"
+define KY_NREJECT "nreject"
+define KY_LOREJECT "loreject"
+define KY_HIREJECT "hireject"
+define KY_GAIN "gain"
+define KY_READNOISE "readnoise"
+define KY_NREGIONS "nregions"
+