source: trunk/InputComplete @ 167

Last change on this file since 167 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
RevLine 
[71]1# FULL PARAMETER LIST
2# -------------------
3#
[142]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.
[71]8#
9#
[3]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
[142]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)
[71]17# reconFile -- the location of the FITS file containing the reconstructed image
[3]18
19imageFile       your-file-here
[71]20flagSubsection  0
[3]21subsection      [*,*,*]
[71]22flagReconExists 0
23reconFile       your-reconstructed-file-here
[3]24
25# OUTPUT RELATED
[142]26# flagLog -- log the intermediate results?   
27# logfile-- the file to put that in.
[3]28# outFile -- the final output list.
29# spectraFile -- the postscript file of spectra
[142]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
[17]36# flagMaps -- save postscript versions of the detection and 0th moment maps?
[3]37# detectionMap, momentMap -- the postscript files produced.
38
39flagLog         1
[87]40logFile         duchamp-Logfile.txt
41outFile         duchamp-Results.txt
42spectraFile     duchamp-Spectra.ps
[3]43flagOutputRecon 0
44flagOutputResid 0
45flagVOT         0
[87]46votFile         duchamp-Results.xml
[26]47flagKarma       0
[87]48karmaFile       duchamp-Results.ann
[3]49flagMaps        1
[87]50detectionMap    duchamp-DetectionMap.ps
51momentMap       duchamp-MomentMap.ps
[3]52
53# FIXING UP THE CUBE
54# flagBlankPix -- whether to replace & ignore blank pixels.
[142]55# blankPixValue -- what value they should take. This is overridden by
56#                  FITS header information.
[3]57# flagBaseline -- whether to subtract spectral baselines before searching
[142]58# flagMW -- ignore a range of channels (to remove Milky Way emission)?
59# minMW, maxMW -- the first and last channels to be ignored
[3]60
61flagBlankPix    1
62blankPixValue   -8.00061
[71]63flagMW          0
[3]64minMW           75
65maxMW           112
[87]66flagBaseline    0
[3]67
68# GENERAL DETECTION RELATED
[142]69# flagNegative -- Are the features being searched for negative (set to true)
70#                 or positive (false -- the default)?
[3]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
[63]75flagNegative    0
[71]76snrCut          3.
[3]77flagGrowth      0
[87]78growthCut       2.
[3]79
80# RECONSTRUCTION RELATED
81# flagATrous -- Whether or not to do the reconstruction before searching
[142]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
[3]85# snrRecon -- The threshold used in filtering the wavelet coefficient arrays.
[142]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.
[3]89
90flagATrous      1
[103]91reconDim        3
[3]92scaleMin        1
[71]93snrRecon        4.
[106]94filterCode      1
[3]95
96# FALSE DISCOVERY RATE METHOD
97# flagFDR -- Whether or not to use the false discovery rate method instead of
98#   simple sigma clipping.
[142]99# alphaFDR -- The "alpha" parameter for the FDR method -- what desired
100#             percentage of discoveries will be false. Expressed as a decimal.
[3]101
102flagFDR         0
[71]103alphaFDR        0.01
[3]104
105# MERGING PARAMETERS
[142]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
[3]116
117flagAdjacent    1
118threshSpatial   3
119threshVelocity  7
[71]120minChannels     3
121minPix          2
[3]122
123# OTHER PARAMETERS
[142]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
[3]136verbose         1
137drawBorders     1
[142]138drawBlankEdges  1
[64]139spectralMethod  peak
[103]140spectralUnits   km/s
Note: See TracBrowser for help on using the repository browser.