source: tags/release-1.0.5/InputComplete

Last change on this file was 142, checked in by Matthew Whiting, 18 years ago

Large range of changes -- related to image cutouts and spectral units mostly.

duchamp.hh -- changed name of spectral type, and introduced one for frequency.
Cubes/getImage.cc -- changed the way the spectral type and spectral units are

looked at. Now is able to use frequency units (defaults
to MHz) when there isn't a good velocity transformation
possible (ie. if there is no restfreq defined).

Cubes/drawMomentCutout.cc -- Range of changes here:

  • improved the way unwanted pixels (those BLANK and outside image coundaries) are dealt with in the image array in drawMomentCutout
  • removed the way the blank pixel information was changed. Now behaves in a consistent way whether or not there are blank pixels
  • improved call to plotBlankEdges()
  • added call to drawFieldEdge() -- a new function that draws a yellow line around the boundary of the pixels of the image.
  • improved the tick mark that shows the angular scale of the cutout. Now adaptable to any pixel scale.
  • added calls to cpgtest() to all functions

Also these files were changed in relation to these edits:
Utils/drawBlankEdges.cc -- improved execution, with "blank" array.
Cubes/cubes.cc -- added call to Param::drawBlankEdge in Cube::plotBlankEdges()
Cubes/outputSpectra.cc -- moved spectra away from left and right axes.
param.cc -- added necessary calls for the new parameter. Other minor changes

to formatting. Added a missed call to stringize.

mainDuchamp.cc -- added #include <param.hh>
param.hh -- Added a new parameter: blankEdge
Cubes/cubes.hh -- improved appearance of comments
Cubes/plots.hh -- improved appearance of comments
InputComplete? -- added new parameter.
docs/Guide.tex -- added text about blank edge plotting.
All six images were changed as well.

CHANGES -- some of these changes -- not up to date yet!

File size: 5.2 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 image 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 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?   
27# logfile-- the file to put that in.
28# outFile -- the final output list.
29# spectraFile -- the postscript file of spectra
30# flagOutputRecon/Resid -- whether to save the reconstruction & residual arrays
31#                          as FITS files
32# flagVOT -- make a VOTable file of the results? 
33# votFile -- the file to put it in
34# flagKarma -- make a Karma annotation file of the results? 
35# karmaFile -- the file to put it in
36# flagMaps -- save postscript versions of the detection and 0th moment maps?
37# detectionMap, momentMap -- the postscript files produced.
38
39flagLog         1
40logFile         duchamp-Logfile.txt
41outFile         duchamp-Results.txt
42spectraFile     duchamp-Spectra.ps
43flagOutputRecon 0
44flagOutputResid 0
45flagVOT         0
46votFile         duchamp-Results.xml
47flagKarma       0
48karmaFile       duchamp-Results.ann
49flagMaps        1
50detectionMap    duchamp-DetectionMap.ps
51momentMap       duchamp-MomentMap.ps
52
53# FIXING UP THE CUBE
54# flagBlankPix -- whether to replace & ignore blank pixels.
55# blankPixValue -- what value they should take. This is overridden by
56#                  FITS header information.
57# flagBaseline -- whether to subtract spectral baselines before searching
58# flagMW -- ignore a range of channels (to remove Milky Way emission)?
59# minMW, maxMW -- the first and last channels to be ignored
60
61flagBlankPix    1
62blankPixValue   -8.00061
63flagMW          0
64minMW           75
65maxMW           112
66flagBaseline    0
67
68# GENERAL DETECTION RELATED
69# flagNegative -- Are the features being searched for negative (set to true)
70#                 or positive (false -- the default)?
71# snrCut --  How many sigma above mean is a detection when sigma-clipping
72# flagGrowth -- Should the detections be "grown" to a lower significance value?
73# growthCut -- The lower threshold used when growing detections
74
75flagNegative    0
76snrCut          3.
77flagGrowth      0
78growthCut       2.
79
80# RECONSTRUCTION RELATED
81# flagATrous -- Whether or not to do the reconstruction before searching
82# reconDim -- The number of dimensions in which to perform the reconstruction.
83# scaleMin -- The minimum scale (starts at 1) to be included in the
84#              reconstruction
85# snrRecon -- The threshold used in filtering the wavelet coefficient arrays.
86# filterCode -- The code number for the choice of filter to be used in the
87#               reconstruction:  1 = B3-spline filter, 2 = Triangle function,
88#                                3 = Haar wavelet. Other numbers default to 2.
89
90flagATrous      1
91reconDim        3
92scaleMin        1
93snrRecon        4.
94filterCode      1
95
96# FALSE DISCOVERY RATE METHOD
97# flagFDR -- Whether or not to use the false discovery rate method instead of
98#   simple sigma clipping.
99# alphaFDR -- The "alpha" parameter for the FDR method -- what desired
100#             percentage of discoveries will be false. Expressed as a decimal.
101
102flagFDR         0
103alphaFDR        0.01
104
105# MERGING PARAMETERS
106# flagAdjacent -- Whether to use the "adjacent criterion" to judge if objects
107#                  are to be merged.
108# threshSpatial -- If flagAdjacent=false, this is the maximum spatial
109#                   separation between objects for them to be merged.
110# threshVelocity -- The maximum channel separation between objects for them
111#                   to be merged.
112# minChannels -- The minimum number of consecutive channels an object must have
113#                 for it to be accepted.
114# minPix -- Minimum number of spatial pixels a detected object must have for
115#            it to be counted
116
117flagAdjacent    1
118threshSpatial   3
119threshVelocity  7
120minChannels     3
121minPix          2
122
123# OTHER PARAMETERS
124# verbose -- whether to provide progress indicators on the terminal during
125#             execution
126# drawBorders -- whether to draw borders around the detections on the moment
127#                 maps in the output.
128# drawBlankEdges -- whether to draw an outline around any BLANK regions in the
129#                    moment maps and image cutouts.
130# spectralMethod -- how to plot the spectra in the output -- the spectrum of
131#                    the peak pixel ("peak" -- the default), or integrated over
132#                    all spatial pixels present ("sum")
133# spectralUnits -- what units you want the spectral axis to be plotted in and
134#                   values quoted in (including the units of integrated flux)
135
136verbose         1
137drawBorders     1
138drawBlankEdges  1
139spectralMethod  peak
140spectralUnits   km/s
Note: See TracBrowser for help on using the repository browser.