00001 #ifndef BEZIER_H
00002 #define BEZIER_H
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031 #include <qfits.h>
00032 #include "spectrum_ops.h"
00033
00034
00035
00036
00037
00038 #define cubeNONEIGHBOR -1000000
00039
00040 #define cubePT_BADPIXEL 0
00041 #define cubePT_USE 1
00042 #define cubePT_FIND 2
00043 #define cubePT_NOINFO 3
00044
00045 typedef struct dim_ {
00046 unsigned short x;
00047 unsigned short y;
00048 unsigned short z;
00049 } Dim;
00050
00051 typedef struct xyzw_ {
00052 double x;
00053 double y;
00054 double z;
00055 double w;
00056 } XYZW;
00057
00058
00059 typedef struct lookup_ {
00060 OneImage *X;
00061 OneImage *hX;
00062 OneImage *Y;
00063 OneImage *Z;
00064 OneCube *id;
00065 OneCube *jd;
00066 } Lookup;
00067
00068
00069
00070 #define i2(im,X,Y) X+Y*im->lx
00071 #define i3(im,X,Y,Z) X+Y*im->lx+Z*im->lx*im->ly
00072
00073
00074
00075
00076
00077
00078 OneImage * cbezierInterpolateImage( OneImage *im, OneImage *mask, Lookup *look, short rx, short ry, short rz,
00079 int max_rad , float ** slit_edges );
00080
00081
00082 OneImage * cbezierFindBad( OneImage *im, OneImage *mask, short rx, short ry, short rz,
00083 short lowerI, short highI, short lowerJ, short highJ, float factor );
00084
00085
00086
00087
00088 float cbezierCorrectPixel( int ipos ,int jpos,
00089 OneImage *im,
00090 OneImage *mask,
00091 OneCube *sc_im,
00092 OneCube *drs_sc_mask,
00093 Lookup *look, short rx, short ry, short rz );
00094 float cbezierCorrectPixel2D( int ipos ,int jpos,
00095 OneImage *im,
00096 OneImage *mask,
00097 OneCube *sc_im,
00098 OneCube *drs_sc_mask,
00099 short rx, short ry, short rz , double *med, double *stdev, float factor );
00100
00101 float cbezierInterpol( OneCube * im, OneCube * action );
00102 int Bezier( XYZW *p, int n, double mu, double munk, XYZW *res );
00103
00104 float cbezierCorrectCosmic( int ipos ,int jpos,
00105 OneImage *im,
00106 OneImage *mask,
00107 OneCube *sc_im,
00108 OneCube *drs_sc_mask,
00109 Lookup *look, short rx, short ry,
00110 short rz , double *med, double *stdev, float factor );
00111
00112
00113 OneImage * cbezierFindCosmic( OneImage *im, OneImage *mask, Lookup *look,
00114 short rx, short ry, short rz,short lowerI,
00115 short highI, short lowerJ, short highJ, float factor );
00116
00117
00118
00119 int cCreateXYZ( Lookup *l );
00120
00121
00122
00123 Lookup * newLookup( void ) ;
00124
00125
00126
00127 void destroyLookup( Lookup *l ) ;
00128
00129
00130
00131 int changeMask( OneImage * mask, OneImage * im ) ;
00132
00133
00134
00135 #endif
00137