1 /*****************************************************************************
2 * ** FILE IDENTIFICATION
4 * ** Name: phm2helix.cpp
5 * ** Purpose: Take projections of a phantom object
6 * ** Programmer: Ian Kay
7 * ** Date Started: Aug 2001
9 * ** This is part of the CTSim program
10 * ** Copyright (C) 1983-2000 Kevin Rosenberg
12 * ** $Id: pjHinterp.cpp,v 1.1 2001/09/24 09:40:42 kevin Exp $
14 * ** This program is free software; you can redistribute it and/or modify
15 * ** it under the terms of the GNU General Public License (version 2) as
16 * ** published by the Free Software Foundation.
18 * ** This program is distributed in the hope that it will be useful,
19 * ** but WITHOUT ANY WARRANTY; without even the implied warranty of
20 * ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
21 * ** GNU General Public License for more details.
23 * ** You should have received a copy of the GNU General Public License
24 * ** along with this program; if not, write to the Free Software
25 * ** Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
26 * ******************************************************************************/
32 enum { O_INTERPVIEW, O_VERBOSE, O_TRACE, O_HELP, O_DEBUG, O_VERSION};
35 static struct option my_options[] =
37 {"interpview", 1, 0, O_INTERPVIEW},
38 {"trace", 1, 0, O_TRACE},
39 {"debug", 0, 0, O_DEBUG},
40 {"verbose", 0, 0, O_VERBOSE},
41 {"help", 0, 0, O_HELP},
42 {"version", 0, 0, O_VERSION},
46 static const char* g_szIdStr = "$Id: pjHinterp.cpp,v 1.1 2001/09/24 09:40:42 kevin Exp $";
48 void pjHinterp_usage ( const char *program )
50 std::cout << "usage: " << fileBasename(program) << " raysum-file interp-raysum-file [OPTIONS]" << std::endl;
51 std::cout << "Interpolation of helical data in raw data space" << std::endl;
52 std::cout << " raysum-file Input raysum file" << std::endl;
53 std::cout << " interp-file Output interpolated raysum file " << std::endl;
54 std::cout << " --trace Set tracing to level" << std::endl;
55 std::cout << " none No tracing (default)" << std::endl;
56 std::cout << " console Text level tracing" << std::endl;
57 std::cout << " phantom Trace phantom" << std::endl;
58 std::cout << " proj Trace allrays" << std::endl;
59 std::cout << " plot Trace plotting" << std::endl;
60 std::cout << " clipping Trace clipping" << std::endl;
61 std::cout << " --verbose Turn on verbose mode" << std::endl;
62 std::cout << " --debug Turn on debug mode" << std::endl;
63 std::cout << " --version Print version" << std::endl;
64 std::cout << " --help Print this help message" << std::endl;
69 pjHinterp_main(int argc, char * const argv[])
71 Projections projGlobal;
72 char* pszProjFilename = NULL;
73 char* pszInterpFilename = NULL;
74 bool bOptVerbose = false;
76 int optTrace = Trace::TRACE_NONE;
79 int opt_interpview=-1;
82 int c = getopt_long(argc, argv, "", my_options, NULL);
89 opt_interpview = strtol(optarg, &endptr, 10);
90 endstr = optarg + strlen(optarg);
91 if (endptr != endstr) {
92 std::cerr << "Error setting --interpview to %s" << optarg << std::endl;
93 pjHinterp_usage(argv[0]);
104 if ((optTrace = Trace::convertTraceNameToID(optarg))
105 == Trace::TRACE_INVALID) {
106 pjHinterp_usage(argv[0]);
112 std::cout << "Version " << VERSION << std::endl <<
113 g_szIdStr << std::endl;
115 std::cout << "Unknown version number" << std::endl;
121 pjHinterp_usage(argv[0]);
124 pjHinterp_usage(argv[0]);
129 if (optind + 2 != argc) {
130 pjHinterp_usage(argv[0]);
134 pszProjFilename = argv[optind];
136 pszInterpFilename = argv[optind + 1];
138 Projections projections;
139 if ( projections.read(pszProjFilename) != true ){
140 std::cerr << "Error reading input file " << pszProjFilename << std::endl;
145 projections.printScanInfo(os);
146 std::cout << os.str();
149 int status = projections.Helical180LI(opt_interpview);
150 if ( status != 0 ) return (1);
151 status = projections.HalfScanFeather();
152 if ( status != 0 ) return (1);
153 projections.write( pszInterpFilename );
161 main (int argc, char* argv[])
166 retval = pjHinterp_main(argc, argv);
167 } catch (exception e) {
168 std::cerr << "Exception: " << e.what() << std::endl;
170 std::cerr << "Unknown exception" << std::endl;