source: branches/fitshead-branch/duchamp.hh @ 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: 2.1 KB
Line 
1#include <string.h>
2
3// how to convey whether a function has worked
4enum OUTCOME {SUCCESS, FAILURE};
5
6
7const std::string ERR_USAGE_MSG =
8"Usage:: Duchamp -p [parameter file]\n\
9  Other options:\n\
10    -f <file>  Use default parameters with imageFile=<file>\n\
11    -v         Version number\n\
12    -h         This help information\n";
13
14const std::string ERR_USAGE_MSG_SHORT = "Usage:: Duchamp -p [parameter file]\n";
15
16const std::string PROGNAME = "Duchamp";
17
18const std::string VERSION = "1.0";
19
20// The spectral type that we want the wcsprm structs to be in
21const char duchampSpectralType[9] = "VELO-F2V";
22
23// The following are the FITS Header Keywords corresponding to the
24//  parameters related to the atrous reconstruction.
25const std::string keyword_flagMW     = "DU_MWFLG";
26const std::string keyword_minMW      = "DU_MWMIN";
27const std::string keyword_maxMW      = "DU_MWMAX";
28const std::string keyword_scaleMin   = "DU_MINSC";
29const std::string keyword_snrRecon   = "DU_ATCUT";
30const std::string keyword_reconDim   = "DU_ATDIM";
31const std::string keyword_filterCode = "DU_FILTR";
32const std::string keyword_ReconResid = "DU_RECON"; // is this file the reconstruction or the residual?
33// And these are the comments corresponding to the relevant keywords
34const std::string comment_flagMW     = "Duchamp parameter flagMW";
35const std::string comment_minMW      = "Duchamp parameter minMW";
36const std::string comment_maxMW      = "Duchamp parameter maxMW";
37const std::string comment_scaleMin   = "Duchamp parameter scaleMin";
38const std::string comment_snrRecon   = "Duchamp parameter snrRecon";
39const std::string comment_reconDim   = "Duchamp parameter reconDim";
40const std::string comment_filterCode = "Duchamp parameter filterCode";
41const std::string comment_ReconResid = "Is this the reconstruction or residual?";
42// Descriptive Headers:
43const std::string header_history1 = "Reconstructed with a trous wavelet technique";
44const std::string header_history2 = "Reconstruction by Duchamp v." + VERSION;
45const std::string header_comment  = "The following are the Duchamp parameters used in reconstruction";
46const std::string flagMW_comment  = "Channels containing Galactic emission have been removed";
Note: See TracBrowser for help on using the repository browser.