1 /*****************************************************************************
5 ** Purpose: Display linogram sampling
6 ** Programmer: Kevin Rosenberg
7 ** Date Started: April 2003
9 ** This is part of the CTSim program
10 ** Copyright (C) 1983-2009 Kevin Rosenberg
12 ** This program is free software; you can redistribute it and/or modify
13 ** it under the terms of the GNU General Public License (version 2) as
14 ** published by the Free Software Foundation.
16 ** This program is distributed in the hope that it will be useful,
17 ** but WITHOUT ANY WARRANTY; without even the implied warranty of
18 ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
19 ** GNU General Public License for more details.
21 ** You should have received a copy of the GNU General Public License
22 ** along with this program; if not, write to the Free Software
23 ** Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
24 ******************************************************************************/
28 enum { O_XY, O_POLAR_RT, O_VERBOSE, O_HELP, O_VERSION, O_DEBUG };
30 static struct option my_options[] =
33 {"polar-rt", 0, 0, O_POLAR_RT},
34 {"debug", 0, 0, O_DEBUG},
35 {"verbose", 0, 0, O_VERBOSE},
36 {"help", 0, 0, O_HELP},
37 {"version", 0, 0, O_VERSION},
41 static const char* g_szIdStr = "$Id$";
45 linogram_usage (const char *program)
47 std::cout << "usage: " << fileBasename(program) << " n d [OPTIONS]\n";
48 std::cout << "Imagefile information\n";
49 std::cout << std::endl;
50 std::cout << " n Linogram N\n";
51 std::cout << " d Max detector spacing\n";
52 std::cout << " --xy Output x,y pairs\n";
53 std::cout << " --polar-rt Output r,t pairs\n";
54 std::cout << " --debug Debug mode\n";
55 std::cout << " --verbose Verbose mode\n";
56 std::cout << " --version Print version\n";
57 std::cout << " --help Print this help message\n";
61 linogram_main (int argc, char *const argv[])
70 int c = getopt_long (argc, argv, "", my_options, NULL);
91 std::cout << "Version " << VERSION << std::endl << g_szIdStr << std::endl;
93 std::cout << "Unknown version number\n";
98 linogram_usage(argv[0]);
101 linogram_usage(argv[0]);
106 if (optind + 2 != argc) {
107 linogram_usage (argv[0]);
111 const char* in_n = argv[optind];
112 const char* in_d = argv[optind+1];
115 double d = atof (in_d);
116 int size = 4 * n + 3;
119 double theta_base = PI/4;
122 double theta_vec [size];
123 for (int i = 0; i < size; i++) {
124 int m = i - (2 * n + 1);
125 theta_vec[i] = atan (static_cast<double>(2 * m) / size);
130 for (m = 0; m < size; m++) {
131 double step = d * cos(theta_vec[m]);
132 for (int id = min; id <= max; id++) {
133 double r = id * step;
134 double x = r * cos(theta_vec[m] + theta_base);
135 double y = r * sin(theta_vec[m] + theta_base);
136 printf ("%lf,%lf ", x, y);
141 for (m = 0; m < size; m++) {
142 double step = d * cos(theta_vec[m]);
143 for (int id = min; id <= max; id++) {
144 double r = id * step;
145 double x = r * cos(theta_vec[m] + PI/2 + theta_base);
146 double y = r * sin(theta_vec[m] + PI/2 + theta_base);
147 printf ("%lf,%lf ", x, y);
153 for (m = 0; m < size; m++) {
155 printf ("%lf: ", theta_vec[m] + theta_base);
156 double step = d * cos(theta_vec[m]);
157 for (int id = min; id <= max; id++) {
159 printf ("%lf,", theta_vec[m] + theta_base);
160 printf ("%lf ", id * step);
165 for (m = 0; m < size; m++) {
167 printf ("%lf: ", theta_vec[m] + PI/2 + theta_base);
168 double step = d * cos(theta_vec[m]);
169 for (int id = min; id <= max; id++) {
171 printf ("%lf,", theta_vec[m] + PI/2 + theta_base);
172 printf ("%lf ", id * step);
183 main (int argc, char *const argv[])
188 retval = linogram_main(argc, argv);
189 } catch (exception e) {
190 std::cerr << "Exception: " << e.what() << std::endl;
192 std::cerr << "Unknown exception\n";