source: branches/fitshead-branch/InputComplete @ 1441

Last change on this file since 1441 was 96, checked in by Matthew Whiting, 18 years ago

Two sets of large changes:
1) Added reconDim, to select which dimension of reconstruction to use.
2) Changed the way the MW channels are dealt with -- now not set to 0. but

simply ignored for searching purposes.

Summary of changes for each file:
duchamp.hh -- added keyword info for reconDim
param.hh -- Introduced reconDim and flagUsingBlank and isInMW.
param.cc -- Introduced reconDim and flagUsingBlank: initialisation etc commands
InputComplete? -- Added reconDim info
mainDuchamp.cc -- Removed the removeMW call. Changed search function names
docs/Guide.tex -- New code to deal with changes: reconDim, MW removal,

up-to-date output examples, better hints & notes section

Detection/thresholding_functions.cc -- minor typo correction

Cubes/cubes.hh -- added isBlank and removeMW functions in Image class, and

renamed search function prototypes

Cubes/cubes.cc -- added removeMW function for Image class, cleaned up

Cube::removeMW as well (although now not used)

Cubes/outputSpectra.cc -- Improved use of isBlank and isInMW functions: now

show MW channels but don't use them in calculating
the flux range

Cubes/getImage.cc -- added line to indicate whether the Param's blank value

is being used, rather than the FITS header.

Cubes/cubicSearch.cc -- renamed functions: front end is now CubicSearch?, and

searching function is search3DArray.

-- Improved way MW removal is dealt with. Changed

location of stats calculations.

Cubes/saveImage.cc -- added code for saving reconDim info
Cubes/readRecon.cc -- added code for reading reconDim info (and added status

intialisations before all cfitsio commands)

ATrous/ReconSearch.cc -- renamed functions: front end is now ReconSearch?, and

searching function is searchReconArray.
Using reconDim to decide which reconstruction to use.

-- Improved way MW removal is dealt with. Changed

location of stats calculations.

ATrous/atrous_1d_reconstruct.cc -- using Median stats
ATrous/atrous_2d_reconstruct.cc -- made code up to date, to conform with 1- &

3-d code. Removed boundary conditions.

ATrous/atrous_3d_reconstruct.cc -- corrected typo (avGapY /= float(xdim)).

Using median stats.

Cubes/cubicSearchNMerge.cc -- Deleted. (not used)
Cubes/readParams.cc -- Deleted. (not used)

File size: 4.7 KB
Line 
1# FULL PARAMETER LIST
2# -------------------
3#
4# This file contains the full list of user-definable parameters. They are grouped
5#  here according to their use, and the values quoted are the defaults -- if a
6#  given parameter is not included in your parameter file, this is the value it
7#  will take.
8#
9#
10# INPUT RELATED
11# imageFile -- the FITS image.
12# flagSubsection -- whether to get a subsection of that image.
13# subsection -- the subsection to read in, in format [x1:x2,y1:y2,z1:z2],
14#               or * to indicate the full range
15# flagReconExists -- set to true if the reconstructed image exists as a FITS file
16#                    (ie. it has been saved from a previous run of Duchamp)
17# reconFile -- the location of the FITS file containing the reconstructed image
18
19imageFile       your-file-here
20flagSubsection  0
21subsection      [*,*,*]
22flagReconExists 0
23reconFile       your-reconstructed-file-here
24
25# OUTPUT RELATED
26# flagLog -- log the intermediate results?   logfile-- the file to put that in.
27# outFile -- the final output list.
28# spectraFile -- the postscript file of spectra
29# flagOutputRecon/Resid -- whether to save the reconstruction & residual arrays as FITS files
30# flagVOT -- make a VOTable file of the results?   votFile -- the file to put it in
31# flagKarma -- make a Karma annotation file of the results?  karmaFile -- the file to put it in
32# flagMaps -- save postscript versions of the detection and 0th moment maps?
33# detectionMap, momentMap -- the postscript files produced.
34
35flagLog         1
36logFile         duchamp-Logfile.txt
37outFile         duchamp-Results.txt
38spectraFile     duchamp-Spectra.ps
39flagOutputRecon 0
40flagOutputResid 0
41flagVOT         0
42votFile         duchamp-Results.xml
43flagKarma       0
44karmaFile       duchamp-Results.ann
45flagMaps        1
46detectionMap    duchamp-DetectionMap.ps
47momentMap       duchamp-MomentMap.ps
48
49# FIXING UP THE CUBE
50# flagBlankPix -- whether to replace & ignore blank pixels.
51# blankPixValue -- what value they should take. This is overridden by FITS header information.
52# flagBaseline -- whether to subtract spectral baselines before searching
53# flagMW -- whether to excise a range of channels (to remove Milky Way emission)
54# minMW, maxMW -- the first and last channels to be excised
55
56flagBlankPix    1
57blankPixValue   -8.00061
58flagMW          0
59minMW           75
60maxMW           112
61flagBaseline    0
62
63# GENERAL DETECTION RELATED
64# flagNegative -- Are the features being searched for negative (set to true) or positive (false)?
65# snrCut --  How many sigma above mean is a detection when sigma-clipping
66# flagGrowth -- Should the detections be "grown" to a lower significance value?
67# growthCut -- The lower threshold used when growing detections
68
69flagNegative    0
70snrCut          3.
71flagGrowth      0
72growthCut       2.
73
74# RECONSTRUCTION RELATED
75# flagATrous -- Whether or not to do the reconstruction before searching
76# reconDim -- The number of dimensions with which to perform the reconstruction.
77# scaleMin -- The minimum scale (starts at 1) to be included in the reconstruction
78# snrRecon -- The threshold used in filtering the wavelet coefficient arrays.
79# filterCode -- The code number for the choice of filter to be used in the reconstruction:
80# 1 = B3-spline filter, 2 = Triangle function, 3 = Haar wavelet. Other numbers default to 2.
81
82flagATrous      1
83reconDim        3
84scaleMin        1
85snrRecon        4.
86filterCode      2
87
88# FALSE DISCOVERY RATE METHOD
89# flagFDR -- Whether or not to use the false discovery rate method instead of
90#   simple sigma clipping.
91# alphaFDR -- The "alpha" parameter for the FDR method -- what desired percentage
92#   of discoveries will be false. Expressed as a decimal.
93
94flagFDR         0
95alphaFDR        0.01
96
97# MERGING PARAMETERS
98# flagAdjacent -- Whether to use the "adjacent criterion" to judge if objects are to be merged.
99# threshSpatial -- If flagAdjacent=false, this is the maximum spatial separation between objects
100#    for them to be merged.
101# threshVelocity -- The maximum channel separation between objects for them to be merged.
102# minChannels -- The minimum number of consecutive channels an object must have for it to be accepted.
103# minPix -- Minimum number of spatial pixels for a detected object to be counted
104
105flagAdjacent    1
106threshSpatial   3
107threshVelocity  7
108minChannels     3
109minPix          2
110
111# OTHER PARAMETERS
112# verbose -- whether to provide progress indicators on the terminal during execution
113# drawBorders -- whether to draw borders around the detections on the moment maps in the output.
114# spectralMethod -- how to plot the spectra in the output -- the spectrum of the peak pixel
115#      ("peak" -- the default), or integrated over all spatial pixels present ("sum")
116# spectralUnits -- what units you want the spectral axis to be plotted in and values quoted in
117#      (including the units of integrated flux)
118       
119verbose         1
120drawBorders     1
121spectralMethod  peak
122spectralUnits   km/s
Note: See TracBrowser for help on using the repository browser.