1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
|
# FLATCOMBINE -- Process and combine flat field CCD images.
procedure flatcombine (input)
string input {prompt="List of flat field images to combine"}
file output="Flat" {prompt="Output flat field root name"}
string combine="average" {prompt="Type of combine operation",
enum="average|median"}
string reject="avsigclip" {prompt="Type of rejection",
enum="none|minmax|ccdclip|crreject|sigclip|avsigclip|pclip"}
string ccdtype="flat" {prompt="CCD image type to combine"}
bool process=yes {prompt="Process images before combining?"}
bool subsets=yes {prompt="Combine images by subset parameter?"}
bool delete=no {prompt="Delete input images after combining?"}
bool clobber=no {prompt="Clobber existing output image?"}
string scale="mode" {prompt="Image scaling",
enum="none|mode|median|mean|exposure"}
string statsec="" {prompt="Image section for computing statistics"}
int nlow=1 {prompt="minmax: Number of low pixels to reject"}
int nhigh=1 {prompt="minmax: Number of high pixels to reject"}
int nkeep=1 {prompt="Minimum to keep (pos) or maximum to reject (neg)"}
bool mclip=yes {prompt="Use median in sigma clipping algorithms?"}
real lsigma=3. {prompt="Lower sigma clipping factor"}
real hsigma=3. {prompt="Upper sigma clipping factor"}
string rdnoise="0." {prompt="ccdclip: CCD readout noise (electrons)"}
string gain="1." {prompt="ccdclip: CCD gain (electrons/DN)"}
string snoise="0." {prompt="ccdclip: Sensitivity noise (fraction)"}
real pclip=-0.5 {prompt="pclip: Percentile clipping parameter"}
real blank=1. {prompt="Value if there are no pixels"}
begin
string ims
ims = input
# Process images first if desired.
if (process == YES)
ccdproc (ims, output="", ccdtype=ccdtype, noproc=no)
# Combine the flat field images.
combine (ims, output=output, plfile="", sigma="", combine=combine,
reject=reject, ccdtype=ccdtype, subsets=subsets, delete=delete,
clobber=clobber, project=no, outtype="real", offsets="none",
masktype="none", blank=blank, scale=scale, zero="none", weight=no,
statsec=statsec, lthreshold=INDEF, hthreshold=INDEF, nlow=nlow,
nhigh=nhigh, nkeep=nkeep, mclip=mclip, lsigma=lsigma, hsigma=hsigma,
rdnoise=rdnoise, gain=gain, snoise=snoise, sigscale=0.1,
pclip=pclip, grow=0)
end
|
