source: trunk/InputComplete @ 278

Last change on this file since 278 was 275, checked in by Matthew Whiting, 17 years ago

Many changes -- almost a week's worth!

  • Major change is to enable spatial smoothing. This incorporated:
    • A new class GaussSmooth? (in eponymous files) to implement this.
    • New parameters concerning setup of Gaussian.
    • New Cube functions to do the smoothing and searching: smoothCube() and smoothNSearch() (the latter not used, but still in the file).
  • Some changes to the statistical calculations:
    • Cube::setCubeStats() now re-worked to be clearer. Smoothed data dealt with directly. setupFDR() now has a switching function to go before it.
    • Improved the way the stats are written to the results file, with new functions.
    • new functions to StatsContainer?: scaleNoise and sigmaToMADFM
    • a new findAllStats() function to return all four statistical parameters, which is made use of by the StatsContainer? calculate() functions.

Other changes include:

  • Lutz_detect now does not check for the minsize criterion.
  • New function Param::makeBlankMask() to make a mask array.
  • FitsHeader? now has a bpaKeyword member
  • Fixed drawContours() function so that it can deal with data outside box.
  • Added new files to Makefile.in
  • Hanning class now has a define(int) function. Definition moved out of the constructor.
  • In the map plots, each object is now marked by the peak location, not the centre.
  • Documentation somewhat updated (bit still to do for smoothing, especially statistics)
  • Other minor changes, including comments.
File size: 6.7 KB
Line 
1# FULL PARAMETER LIST
2# -------------------
3#
4# This file contains the full list of user-definable parameters. They are
5#  grouped here according to their use, and the values quoted are the defaults:
6#  if a given parameter is not included in your parameter file, this is the
7#  value it 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 cube exists as a FITS
16#                    file (ie. it has been saved from a previous Duchamp run)
17# reconFile -- the location of the FITS file containing the reconstructed cube
18# flagSmoothExists -- set to true if the smoothed cube exists as a FITS
19#                    file (ie. it has been saved from a previous Duchamp run)
20# smoothFile -- the location of the FITS file containing the smoothed cube
21
22imageFile       your-file-here
23flagSubsection  0
24subsection      [*,*,*]
25flagReconExists 0
26reconFile       your-reconstructed-file-here
27flagReconExists 1
28reconFile       your-smoothed-file-here
29
30# OUTPUT RELATED
31# flagLog -- log the intermediate results?   
32# logfile-- the file to put that in.
33# outFile -- the final output list.
34# spectraFile -- the postscript file of spectra
35# flagOutputSmooth -- whether to save the smoothed arrays as a FITS file
36# flagOutputRecon/Resid -- whether to save the reconstruction & residual arrays
37#                          as FITS files
38# flagVOT -- make a VOTable file of the results? 
39# votFile -- the file to put it in
40# flagKarma -- make a Karma annotation file of the results? 
41# karmaFile -- the file to put it in
42# flagMaps -- save postscript versions of the detection and 0th moment maps?
43# detectionMap, momentMap -- the postscript files produced.
44# flagXOutput -- display the moment map in a pgplot window
45
46flagLog         1
47logFile         duchamp-Logfile.txt
48outFile         duchamp-Results.txt
49spectraFile     duchamp-Spectra.ps
50flagOutputSmooth 0
51flagOutputRecon 0
52flagOutputResid 0
53flagVOT         0
54votFile         duchamp-Results.xml
55flagKarma       0
56karmaFile       duchamp-Results.ann
57flagMaps        1
58detectionMap    duchamp-DetectionMap.ps
59momentMap       duchamp-MomentMap.ps
60flagXOutput     1
61
62# FIXING UP THE CUBE
63# flagBlankPix -- whether to replace & ignore blank pixels.
64# blankPixValue -- what value they should take. This is overridden by
65#                  FITS header information.
66# flagBaseline -- whether to subtract spectral baselines before searching
67# flagMW -- ignore a range of channels (to remove Milky Way emission)?
68# minMW, maxMW -- the first and last channels to be ignored
69
70flagBlankPix    1
71blankPixValue   -8.00061
72flagMW          0
73minMW           75
74maxMW           112
75flagBaseline    0
76
77# GENERAL DETECTION RELATED
78# flagStatSec -- Whether to only use a subsection of the cube to calculate
79#                the statistics.
80# StatSec -- The subsection used for statistics calculations. It has the
81#            same format as the pixel subsection.
82# flagNegative -- Are the features being searched for negative (set to true)
83#                 or positive (false -- the default)?
84# snrCut --  How many sigma above mean is a detection when sigma-clipping
85# threshold -- The threshold flux dividing source and non-source. Used instead
86#              of calculating it from the cube's statistics. If not specified,
87#              it will be calculated.
88# flagGrowth -- Should the detections be "grown" to a lower significance value?
89# growthCut -- The lower threshold used when growing detections
90
91flagStatSec     0
92StatSec         [*,*,*]
93flagNegative    0
94snrCut          3.
95threshold       0.
96flagGrowth      0
97growthCut       2.
98
99# RECONSTRUCTION RELATED
100# flagATrous -- Whether or not to do the reconstruction before searching
101# reconDim -- The number of dimensions in which to perform the reconstruction.
102# scaleMin -- The minimum scale (starts at 1) to be included in the
103#              reconstruction
104# snrRecon -- The threshold used in filtering the wavelet coefficient arrays.
105# filterCode -- The code number for the choice of filter to be used in the
106#               reconstruction:  1 = B3-spline filter, 2 = Triangle function,
107#                                3 = Haar wavelet. Other numbers default to 2.
108
109flagATrous      1
110reconDim        1
111scaleMin        1
112snrRecon        4.
113filterCode      1
114
115# SMOOTHING
116# flagSmooth -- whether to smooth the cube.
117# smoothType -- either "spectral" or "spatial"
118# hanningWidth -- the width parameter of the Hanning (spectral smoothing)
119#                 function
120# kernMaj -- the FWHM of the major axis of the 2D spatial smoothing gaussian
121# kernMin -- the FWHM of the minor axis of the 2D spatial smoothing gaussian
122# kernPA  -- the position angle of the major axis of the 2D spatial
123#            smoothing gaussian
124
125flagSmooth      false
126smoothType      spectral
127hanningWidth    5
128kernMaj         3.
129kernMin         3.
130kernPA          0.
131
132# FALSE DISCOVERY RATE METHOD
133# flagFDR -- Whether or not to use the false discovery rate method instead of
134#   simple sigma clipping.
135# alphaFDR -- The "alpha" parameter for the FDR method -- what desired
136#             percentage of discoveries will be false. Expressed as a decimal.
137
138flagFDR         0
139alphaFDR        0.01
140
141# MERGING PARAMETERS
142# flagAdjacent -- Whether to use the "adjacent criterion" to judge if objects
143#                  are to be merged.
144# threshSpatial -- If flagAdjacent=false, this is the maximum spatial
145#                   separation between objects for them to be merged.
146# threshVelocity -- The maximum channel separation between objects for them
147#                   to be merged.
148# minChannels -- The minimum number of consecutive channels an object must have
149#                 for it to be accepted.
150# minPix -- Minimum number of spatial pixels a detected object must have for
151#            it to be counted
152
153flagAdjacent    1
154threshSpatial   3
155threshVelocity  7
156minChannels     3
157minPix          2
158
159# OTHER PARAMETERS
160# verbose -- whether to provide progress indicators on the terminal during
161#             execution
162# drawBorders -- whether to draw borders around the detections on the moment
163#                 maps in the output.
164# drawBlankEdges -- whether to draw an outline around any BLANK regions in the
165#                    moment maps and image cutouts.
166# spectralMethod -- how to plot the spectra in the output -- the spectrum of
167#                    the peak pixel ("peak" -- the default), or integrated over
168#                    all spatial pixels present ("sum")
169# spectralUnits -- what units you want the spectral axis to be plotted in and
170#                   values quoted in (including the units of integrated flux)
171
172verbose         1
173drawBorders     1
174drawBlankEdges  1
175spectralMethod  peak
176spectralUnits   km/s
Note: See TracBrowser for help on using the repository browser.