I R A F Image Reduction and Analysis Facility PACKAGE = daophot TASK = daofind image = @v_4.stds Input image(s) output = default Output coordinate file(s) (default: image.coo.?) (starmap= ) Output density enhancement image(s) (skymap = ) Output sky image(s) (datapar= ) Data dependent parameters (findpar= ) Object detection parameters (boundar= nearest) Boundary extension (constant|nearest|reflect|wrap) (constan= 0.) Constant for boundary extension (interac= no) Interactive mode ? (icomman= ) Image cursor: [x y wcs] key [cmd] (gcomman= ) Graphics cursor: [x y wcs] key [cmd] (wcsout = )_.wcsout) The output coordinate system (logical,tv,physical) (cache = )_.cache) Cache the image pixels ? (verify = )_.verify) Verify critical daofind parameters ? (update = )_.update) Update critical daofind parameters ? (verbose= )_.verbose) Print daofind messages ? (graphic= )_.graphics) Graphics device (display= )_.display) Display device (mode = ql) |
I R A F Image Reduction and Analysis Facility PACKAGE = daophot TASK = datapars (scale = 0.33) Image scale in units per pixel (fwhmpsf= 1.50) FWHM of the PSF in scale units (emissio= yes) Features are positive ? (sigma = 4.) Standard deviation of background in counts (datamin= -20.) Minimum good data value (datamax= 65500.) Maximum good data value (noise = poisson) Noise model (ccdread= ) CCD readout noise image header keyword (gain = ) CCD gain image header keyword (readnoi= 5.8) CCD readout noise in electrons (epadu = 2.9) Gain in electrons per count (exposur= EXPTIME) Exposure time image header keyword (airmass= AIRMASS) Airmass image header keyword (filter = WFFBAND) Filter image header keyword (obstime= UTSTART) Time of observation image header keyword (itime = INDEF) Exposure time (xairmas= INDEF) Airmass (ifilter= INDEF) Filter (otime = INDEF) Time of observation (mode = ql) |
I R A F Image Reduction and Analysis Facility PACKAGE = daophot TASK = findpars (thresho= 5.) Threshold in sigma for feature detection (nsigma = 1.5) Width of convolution kernel in sigma (ratio = 1.) Ratio of minor to major axis of Gaussian kernel (theta = 0.) Position angle of major axis of Gaussian kernel (sharplo= 0.2) Lower bound on sharpness for feature detection (sharphi= 1.) Upper bound on sharpness for feature detection (roundlo= -1.) Lower bound on roundness for feature detection (roundhi= 2.) Upper bound on roundness for feature detection (mkdetec= yes) Mark detections on the image display ? (mode = ql) |
I R A F Image Reduction and Analysis Facility PACKAGE = astutil TASK = asthedit images = @coo_stds.images Images to be operated upon commands= asthedit_stds.com File of commands (table = ) File of values (colname= ) Column names in table file (prompt = asthedit> ) Prompt for STDIN commands (update = yes) Update image header? (verbose= yes) Verbose output? (oldstyl= no) Use old style format? (mode = ql) |
#example of a calculation ########################################### ###only needed if entered ST instead of UT #$utmidn = 4:00:01 #$stmidns = mst ($utdate, $utmidn, obsdb ("lco", "longitude")) #$stmidn = real($stmidns) #$ut = $st - $stmidn + $utmidn #if($ut<0.) # $ut = $ut + 24. #endif ###################################### $uts = utstart $dobs = @'date-obs' $nepoch = epoch ($dobs,$uts) $nepoch = $nepoch - 2000. newepoch = format ("%12.7f", $nepoch) |
I R A F Image Reduction and Analysis Facility PACKAGE = daophot TASK = phot image = @phot_in Input image(s) coords = @phot_in_2 Input coordinate list(s) (default: image.coo.?) output = default Output photometry file(s) (default: image.mag.?) skyfile = Input sky value file(s) (plotfil= ) Output plot metacode file (datapar= ) Data dependent parameters (centerp= ) Centering parameters (fitskyp= ) Sky fitting parameters (photpar= ) Photometry parameters (interac= no) Interactive mode ? (radplot= no) Plot the radial profiles? (icomman= ) Image cursor: [x y wcs] key [cmd] (gcomman= ) Graphics cursor: [x y wcs] key [cmd] (wcsin = )_.wcsin) The input coordinate system (logical,tv,physical,world) (wcsout = )_.wcsout) The output coordinate system (logical,tv,physical) (cache = )_.cache) Cache the input image pixels in memory ? (verify = )_.verify) Verify critical phot parameters ? (update = )_.update) Update critical phot parameters ? (verbose= )_.verbose) Print phot messages ? (graphic= stdgraph) Graphics device (display= stdimage) Display device (mode = ql) |
I R A F Image Reduction and Analysis Facility PACKAGE = daophot TASK = datapars (scale = 0.33) Image scale in units per pixel (fwhmpsf= 1.50) FWHM of the PSF in scale units (emissio= yes) Features are positive ? (sigma = 4.) Standard deviation of background in counts (datamin= -20.) Minimum good data value (datamax= 65500.) Maximum good data value (noise = poisson) Noise model (ccdread= ) CCD readout noise image header keyword (gain = ) CCD gain image header keyword (readnoi= 5.8) CCD readout noise in electrons (epadu = 2.9) Gain in electrons per count (exposur= EXPTIME) Exposure time image header keyword (airmass= AIRMASS) Airmass image header keyword (filter = WFFBAND) Filter image header keyword (obstime= NEWEPOCH) Time of observation image header keyword (itime = INDEF) Exposure time (xairmas= INDEF) Airmass (ifilter= INDEF) Filter (otime = INDEF) Time of observation (mode = ql) |
I R A F Image Reduction and Analysis Facility PACKAGE = daophot TASK = centerpars (calgori= none) Centering algorithm (cbox = 3.) Centering box width in scale units (cthresh= 0.) Centering threshold in sigma above background (minsnra= 1.) Minimum signal-to-noise ratio for centering algorithim (cmaxite= 10) Maximum iterations for centering algorithm (maxshif= 1.) Maximum center shift in scale units (clean = no) Symmetry clean before centering (rclean = 1.) Cleaning radius in scale units (rclip = 2.) Clipping radius in scale units (kclean = 3.) K-sigma rejection criterion in skysigma (mkcente= no) Mark the computed center (mode = ql) |
I R A F Image Reduction and Analysis Facility PACKAGE = daophot TASK = fitskypars (salgori= median) Sky fitting algorithm (annulus= 10.) Inner radius of sky annulus in scale units (dannulu= 5.) Width of sky annulus in scale units (skyvalu= 0.) User sky value (smaxite= 10) Maximum number of sky fitting iterations (sloclip= 0.) Lower clipping factor in percent (shiclip= 0.) Upper clipping factor in percent (snrejec= 50) Maximum number of sky fitting rejection iterations (sloreje= 3.) Lower K-sigma rejection limit in sky sigma (shireje= 3.) Upper K-sigma rejection limit in sky sigma (khist = 3.) Half width of histogram in sky sigma (binsize= 0.1) Binsize of histogram in sky sigma (smooth = no) Boxcar smooth the histogram (rgrow = 0.) Region growing radius in scale units (mksky = no) Mark sky annuli on the display (mode = ql) |
I R A F Image Reduction and Analysis Facility PACKAGE = daophot TASK = photpars (weighti= constant) Photometric weighting scheme (apertur= 5,6,7,8,10) List of aperture radii in scale units (zmag = 25.) Zero point of magnitude scale (mkapert= no) Draw apertures on the display (mode = ql) |
I R A F Image Reduction and Analysis Facility PACKAGE = daophot TASK = txdump textfile= Input apphot/daophot text database(s) fields = id,xcenter,ycenter,itime,otime,xairmass,ifilter,mag,merr Fields to be extracted expr = PIER[5]=0 Boolean expression for record selection (headers= no) Print the field headers ? (paramet= yes) Print the parameters if headers is yes ? (mode = ql) Mode of task |
textfiles,s,a,"",,,"Input apphot/daophot text database(s)" fields,s,a,"id,xcenter,ycenter,itime,otime,xairmass,ifilter,mag,merr",,,"Fields to be extracted" expr,s,a,"PIER[5]=0",,,"Boolean expression for record selection" headers,b,h,no,,,"Print the field headers ?" parameters,b,h,yes,,,"Print the parameters if headers is yes ?" mode,s,h,"ql",,,"Mode of task" |
textfiles,s,a,,,,"Input apphot/daophot text database(s)" fields,s,h,"id,xcenter,ycenter,itime,otime,xairmass,ifilter,mag,merr",,,"Fields to be extracted" expr,s,h,"yes",,,"Boolean expression for record selection" headers,b,h,no,,,"Print the field headers ?" parameters,b,h,yes,,,"Print the parameters if headers is yes ?" mode,s,h,"ql",,,"Mode of task" |
txdump t421679.fit4.mag.1 > t421679.fit4.txd txdump t421686.fit4.mag.1 > t421686.fit4.txd txdump t421687.fit4.mag.1 > t421687.fit4.txd txdump t421691.fit4.mag.1 > t421691.fit4.txd txdump t421692.fit4.mag.1 > t421692.fit4.txd txdump t421711.fit4.mag.1 > t421711.fit4.txd txdump t421715.fit4.mag.1 > t421715.fit4.txd txdump t421716.fit4.mag.1 > t421716.fit4.txd txdump t421858.fit4.mag.1 > t421858.fit4.txd txdump t421862.fit4.mag.1 > t421862.fit4.txd txdump t421863.fit4.mag.1 > t421863.fit4.txd txdump t421867.fit4.mag.1 > t421867.fit4.txd txdump t421870.fit4.mag.1 > t421870.fit4.txd txdump t421956.fit4.mag.1 > t421956.fit4.txd txdump t421959.fit4.mag.1 > t421959.fit4.txd txdump t421993.fit4.mag.1 > t421993.fit4.txd txdump t422008.fit4.mag.1 > t422008.fit4.txd txdump t422052.fit4.mag.1 > t422052.fit4.txd txdump t422055.fit4.mag.1 > t422055.fit4.txd |
Master Images (V-band) t421686.fit ... SA92 t422254.fit ... SA110 t422681.fit ... SA107 t422503.fit ... SA114 t422225.fit ... PG02 t422905.fit ... SA98 |
I R A F Image Reduction and Analysis Facility PACKAGE = mscred TASK = mscsplit input = t422905.fit List of input MEF files (output = t422905) List of output root names (mefext = .fit) MEF filename extension (delete = no) Delete MEF file after splitting? (verbose= no) Verbose? (fd1 = ) (fd2 = ) (mode = ql) |
I R A F Image Reduction and Analysis Facility PACKAGE = finder TASK = mkgsctab input = Stet_SA114.dat Input coordinate list (RA Dec [ID]) output = t422503_4.fits Output GSC format table (ra_unit= degrees) Input RA units (startid= 1) Starting ID number (if not provided) (region = 0) Region number for new entries (plate = NONE) Plate name for new entries (cdfile = finder$lib/cdfile.gsc) Column definition file (list = Stet_SA114.dat) (mode = ql) |
I R A F Image Reduction and Analysis Facility PACKAGE = finder TASK = tpltsol image = t422503_4.fits Input image name table = t422503_4.fits.tab Input table name database= t422503_4.fits.db Output database name (results= ) Results summary file (imupdat= no) Update the image WCS? (tabupda= no) Update input table? (refitca= no) Recompute XYs for uncentered sources? (verbose= yes) Print verbose progress messages? (dsshead= no) Read plate center from DSS header? (ra_ref = 18:43:13.699) Plate center RA (hours) (dec_ref= +00:32:27.08) Plate center Dec (degrees) (eq_ref = INDEF) Plate center coordinate equinox (inpixsy= logical) Input pixel system (outpixs= logical) Output pixel system (insyste= j2000) Input celestial coordinate system (project= tan) Sky projection geometry (fitgeom= general) Fitting geometry (functio= polynomial) Surface type (xxorder= 4) Order of xi fit in x (xyorder= 4) Order of xi fit in y (xxterms= half) Include cross-terms in xi fit? (yxorder= 4) Order of eta fit in x (yyorder= 4) Order of eta fit in y (yxterms= half) Include cross-terms in eta fit? (reject = 2.) Rejection limit in sigma units (interac= yes) Fit the transformation interactively? (graphic= stdgraph) Default graphics device (cursor = ) Graphics cursor (catpars= ) Catalog description pset (list = ) (mode = ql) |
I R A F Image Reduction and Analysis Facility PACKAGE = finder TASK = tfinder image = t422503_4.fits Image name (rootnam= ) Alternate root name for output files (objects= ) List of program object X,Y coords (scale = 0.33) Plate or image scale ("/pixel) (north = left) Direction of North in the field (east = bottom) Direction of East in the field (ra = 22.694110833333) RA of the reference point (hours) (dec = 1.2038833333333) Dec of the reference point (degrees) (equinox= 2000.) Reference coordinate equinox (xref = 1024.) X coordinate of the reference point (yref = 2048.) Y coordinate of the reference point (date_ob= 2004-09-05) Date of the observation (YYYY-MM-DD) (update = no) Update image header WCS following fit? (interac= yes) Enter interactive image cursor loop? (autocen= no) Center at the catalog coords when entering task? (reselec= yes) Apply selectpars when entering task? (autodis= yes) Redisplay after all-source keystroke command? (verbose= yes) Print a running commentary? (rotate = -0.7) Relative position angle (CCW positive) (boxsize= 9) Centering box full width (edge = 200.) Edge buffer width (pixels) (opaxis = no) Is the reference point on the optical axis? (del_ra = 0.) RA offset of the field center (degrees) (del_dec= 0.) Dec offset of the field center (degrees) (list = ) (mode = ql) |
NORTH.....EAST.....X-REF.....Y-REF.... Chip 1 Left......Bottom...-1089.....2061..... Chip 2 Bottom....Right.....4166......988..... Chip 3 Left......Bottom....3124.....2018..... Chip 4 Left......Bottom....1024.....2048..... |
I R A F Image Reduction and Analysis Facility PACKAGE = imcoords TASK = ccmap input = The input coordinate files database= The output database file (solutio= ) The database plate solution names (images = ) The input images (results= ) The optional results summary files (xcolumn= 1) Column containing the x coordinate (ycolumn= 2) Column containing the y coordinate (lngcolu= 3) Column containing the ra / longitude (latcolu= 4) Column containing the dec / latitude (xmin = INDEF) Minimum logical x pixel value (xmax = INDEF) Maximum logical x pixel value (ymin = INDEF) Minimum logical y pixel value (ymax = INDEF) Maximum logical y pixel value (lngunit= ) Input ra / longitude units (latunit= ) Input dec / latitude units (insyste= j2000) Input celestial coordinate system (refpoin= coords) Source of the reference point definition (lngref = INDEF) Reference point ra / longitude telescope coordinate (latref = INDEF) Reference point dec / latitude telescope coordinate (refsyst= INDEF) Reference point telescope coordinate system (lngrefu= ) Reference point ra / longitude units (latrefu= ) Reference point dec / latitude units (project= tan) Sky projection geometry (fitgeom= general) Fitting geometry (functio= polynomial) Surface type (xxorder= 4) Order of xi fit in x (xyorder= 4) Order of xi fit in y (xxterms= half) Xi fit cross terms type (yxorder= 4) Order of eta fit in x (yyorder= 4) Order of eta fit in y (yxterms= half) Eta fit cross terms type (maxiter= 3) The maximum number of rejection iterations (reject = 3.) Rejection limit in sigma units (update = no) Update the image world coordinate system ? (pixsyst= logical) Input pixel coordinate system (verbose= yes) Print messages about progress of task ? (interac= yes) Fit the transformation interactively ? (graphic= stdgraph) Default graphics device (cursor = ) Graphics cursor (mode = ql) |
I R A F Image Reduction and Analysis Facility PACKAGE = imcoords TASK = cctran input = Stet_SA98.radec The input coordinate files output = Stet_SA98_t422905_4.xy The output coordinate files database= t422905_4.fits.db The input database file solution= t422905_4.fits The input plate solutions (geometr= geometric) Transformation type (linear,geometric) (forward= no) Transform x / y to ra / dec (yes) or vice versa (no) ? (xref = INDEF) The X reference pixel (yref = INDEF) The Y reference pixel (xmag = INDEF) The X axis scale in arcsec per pixel (ymag = INDEF) The Y axis scale in arcsec per pixel (xrotati= INDEF) The X axis rotation angle in degrees (yrotati= INDEF) The Y axis rotation angle in degrees (lngref = INDEF) The ra / longitude reference coordinate in lngunits (latref = INDEF) The dec / latitude reference coordinate in latunits (lngunit= degrees) The input / output ra / longitude reference coordinate units (latunit= degrees) The input / output dec / latitude reference coordinate units (project= tan) The sky projection geometry (xcolumn= 1) Input column containing the x / ra / longitude coordinate (ycolumn= 2) Input column containing the y / dec / latitude coordinate (lngform= ) Output format of the ra / longitude / x coordinate (latform= ) Output format of the dec / latitude / y coordinate (min_sig= 7) Minimum precision of the output coordinates (mode = ql) |
[djo3@fasold INTsept04]$ ./offset_stds Enter file name of: (std star id) <==> (master coord) xref Stet_SA92_t421686_4.xy Enter file of: (star coords on the master frame) t421686.fit4.txm Enter secondary file name ft421685.fit4.txd Enter outputfile name ft421685.fit4.txo enter initial guess p-s 0,0 --- Read std file. t421686.fit4.txm , # Stds = 200 --- --- Read secondary file. ft421685.fit4.txd # Stars = 200 Stet_SA92_t421686_4.xy t421686.fit4.txm ft421685.fit4.txd ft421685.fit4.txo ************************************************ offset, scale factor , rotation in x,y -2.498117 1.000101 2.6654129E-04 1.515641 0.9996654 -5.7293499E-05 ************************************************ sx,sy,sigma,sfac,nrej,n 0.534 0.402 0.473 3.000 4 72 ************************************************ |
as> asttimes.observatory="lapalma" as> asttimes.year=2004 as> asttimes.month=9 as> asttimes.time=16 as> for (i=4; i<15; i+=1) { >>> asttimes (day=i, header=no) >>> } 2004 9 4 SAT 16:00:00.0 16:00:00.0 2004.67671 2453253.1667 13:44:55.3 2004 9 5 SUN 16:00:00.0 16:00:00.0 2004.67944 2453254.1667 13:48:51.9 2004 9 6 MON 16:00:00.0 16:00:00.0 2004.68218 2453255.1667 13:52:48.5 2004 9 7 TUE 16:00:00.0 16:00:00.0 2004.68492 2453256.1667 13:56:45.0 2004 9 8 WED 16:00:00.0 16:00:00.0 2004.68766 2453257.1667 14:00:41.6 2004 9 9 THU 16:00:00.0 16:00:00.0 2004.69039 2453258.1667 14:04:38.1 2004 9 10 FRI 16:00:00.0 16:00:00.0 2004.69313 2453259.1667 14:08:34.7 2004 9 11 SAT 16:00:00.0 16:00:00.0 2004.69587 2453260.1667 14:12:31.2 2004 9 12 SUN 16:00:00.0 16:00:00.0 2004.69861 2453261.1667 14:16:27.8 2004 9 13 MON 16:00:00.0 16:00:00.0 2004.70135 2453262.1667 14:20:24.3 2004 9 14 TUE 16:00:00.0 16:00:00.0 2004.70408 2453263.1667 14:24:20.9 |
L92-S253 394.42 814.40 5.10 4.6865055 1.58 i ft422057.fit4 0.000 16.171 0.019 0 0 0 1 0 0 0 0 0 0 0 |
I R A F Image Reduction and Analysis Facility PACKAGE = photcal TASK = fitparams observat= vtest.obs List of observations files catalogs= Stetson_all.phot List of standard catalog files config = Stet_V_I.con Configuration file paramete= Stet_V_I.out Output parameters file (weighti= equations) Weighting type (uniform,photometric,equations) (addscat= yes) Add a scatter term to the weights ? (toleran= 3.0000000000000E-5) Fit convergence tolerance (maxiter= 15) Maximum number of fit iterations (nreject= 3) Number of rejection iterations (low_rej= 3.) Low sigma rejection factor (high_re= 3.) High sigma rejection factor (grow = 0.) Rejection growing radius (interac= yes) Solve fit interactively ? (logfile= STDOUT) Output log file (log_unm= yes) Log any unmatched stars ? (log_fit= no) Log the fit parameters and statistics ? (log_res= no) Log the results ? (catdir = )_.catdir) The standard star catalog directory (graphic= stdgraph) Output graphics device (cursor = ) Graphics cursor input (mode = ql) |
#This config file is for the Stetson Standard Star list catalog #Tells which column has which data B 2 V 6 R 10 I 14 error(B) 3 error(V) 7 error(R) 11 error(I) 15 #Format from the output of the dan_mkall_noap_sept04.pl script #id 1, x 2, y 3, et 4, ut 5, X 6, filt 7, file 8, apcor 9 m 10 me 11 ... observations xm 6 ut 5 et 4 fi 8 filt 7 apc 9 m 10 error(m) 11 em 11 n1 12 n2 13 n3 14 n4 15 n5 16 n6 17 n7 18 n8 19 n9 20 n10 21 n11 22 #Standard system transformation transformation #a4..a9 not used fit a2 = 0, a3 = .12, a11=0.1, a12=0.1, a13=0.1, a14=0.1, a15=0.1, a16=0.1, a17=0.1, a18=0.1, a19=0.1, a20=0.1, a21=0.1 #const a11 = 0 const a12 = 0 const a13 = 0 const a14 = 0 const a15 = 0 const a16 = 0 #const a17 = 0 #const a18 = 0 const a19 = 0 const a20 = 0 #const a21 = 0 #average extinction coefficients for la palma #U = 0.46, B = 0.22, V = 0.12, R = 0.08, I = 0.04, Z = 0.05, #g' = 0.19, r' = 0.09, i' = 0.05 const a3 = .12 #V extinction default la palma value #const a2 = .01 #correct the counts for the shutter delay of 0.05secs #set p = (25 - mi) / 2.5 #set cts = et * (10**p) #set m = 25 - 2.5 * log10(cts/(et + 0.05)) set VI = V - I set VR = V - R set BV = B - V set mc = m - 2.5*log10(et) EQ1 : m = V + a2*(VI) + a3*(xm) + n1*a11 + n2*a12 + n3*a13 + n4*a14 + n5*a15 + n6*a16 + n7*a17 + n8*a18 + n9*a19 + n10*a20 + n11*a21 #EQ1 : 25 - 2.5*log10(cts) = #- 2.5*log10(et + c) + V + a2*(VI) + a3*(xm) + #n1*a11 + n2*a12 + n3*a13 + #n4*a14 + n5*a15 + n6*a16 + #n7*a17 + n8*a18 + n9*a19 weight(EQ1) = 1.0 / (em + .001)**2 #deriv (EQ1, a1) = 0. #deriv (EQ1, a3) = xm #deriv (EQ1, a2) = VI #deriv (EQ1, a11) = n1 #deriv (EQ1, a12) = n2 #deriv (EQ1, a13) = n3 plot (EQ1) = m - (V + a2*VI + a3*(xm) + n1*a11 + n2*a12 + n3*a13 + n4*a14 + n5*a15 + n6*a16 + n7*a17 + n8*a18 + n9*a19 + n10*a20 + n11*a21 ) , V plot (EQ1) = mc , m - (V + a2*VI + a3*(xm) + n1*a11 + n2*a12 + n3*a13 + n4*a14 + n5*a15 + n6*a16 + n7*a17 + n8*a18 + n9*a19 + n10*a20 + n11*a21) plot (EQ1) = VI , m - (V + a2*VI + a3*(xm) + n1*a11 + n2*a12 + n3*a13 + n4*a14 + n5*a15 + n6*a16 + n7*a17 + n8*a18 + n9*a19 + n10*a20 + n11*a21 ) |
I R A F Image Reduction and Analysis Facility PACKAGE = photcal TASK = chkconfig config = Stet_V_I.con Input configuration file (verbose= no) Verbose output ? (mode = ql) |