X-Git-Url: http://git.kpe.io/?p=ctsim.git;a=blobdiff_plain;f=src%2Fphm2if.cpp;h=573cb7c697525de4790f1143c8fd8b5d72b4adff;hp=3e82eefaadebc2802b5a430352af2f4ef191e13a;hb=2c61ff85796550481227f2fbec53506a6b5bd365;hpb=f4a23943110823118f35756dd41fbd6707f04511 diff --git a/src/phm2if.cpp b/src/phm2if.cpp index 3e82eef..573cb7c 100644 --- a/src/phm2if.cpp +++ b/src/phm2if.cpp @@ -9,7 +9,7 @@ ** This is part of the CTSim program ** Copyright (C) 1983-2000 Kevin Rosenberg ** -** $Id: phm2if.cpp,v 1.9 2000/06/15 19:07:10 kevin Exp $ +** $Id: phm2if.cpp,v 1.10 2000/06/18 10:27:11 kevin Exp $ ** ** This program is free software; you can redistribute it and/or modify ** it under the terms of the GNU General Public License (version 2) as @@ -112,13 +112,13 @@ phm2if_main (int argc, char* argv[]) DomainType opt_filter_domain = D_SPATIAL; char *opt_outfile = NULL; int opt_debug = 0; - char str[256]; - char opt_desc[256], opt_phmfilename[256]; + string opt_desc; + string opt_phmfilename; char *endstr, *endptr; double opt_filter_param = 1; double opt_filter_bw = 1.; int opt_trace = TRACE_NONE; - int opt_verbose = 0; + bool opt_verbose = false; #ifdef HAVE_MPI ImageFile* imLocal = NULL; MPIWorld mpiWorld (argc, argv); @@ -129,8 +129,6 @@ phm2if_main (int argc, char* argv[]) #ifdef HAVE_MPI if (mpiWorld.getRank() == 0) { #endif - strcpy(opt_desc, ""); - strcpy(opt_phmfilename, ""); while (1) { int c = getopt_long(argc, argv, "", my_options, NULL); char *endptr = NULL; @@ -147,8 +145,8 @@ phm2if_main (int argc, char* argv[]) } break; case O_PHMFILE: - strncpy(opt_phmfilename, optarg, sizeof(opt_phmfilename)); - phm.createFromFile(opt_phmfilename); + opt_phmfilename = optarg; + phm.createFromFile(opt_phmfilename.c_str()); #ifdef HAVE_MPI if (mpiWorld.getRank() == 0) cerr << "Can't use phantom from file in MPI mode" << endl; @@ -156,7 +154,7 @@ phm2if_main (int argc, char* argv[]) #endif break; case O_VERBOSE: - opt_verbose = 1; + opt_verbose = true; break; case O_DEBUG: opt_debug = 1; @@ -180,13 +178,13 @@ phm2if_main (int argc, char* argv[]) } break; case O_DESC: - strncpy(opt_desc, optarg, sizeof(opt_desc)); + opt_desc = optarg; break; case O_FILTER_BW: opt_filter_bw = strtod(optarg, &endptr); endstr = optarg + strlen(optarg); if (endptr != endstr) { - fprintf(stderr,"Error setting --filter-bw to %s\n", optarg); + sys_error(ERR_SEVERE,"Error setting --filter-bw to %s\n", optarg); phm2if_usage(argv[0]); return (1); } @@ -195,7 +193,7 @@ phm2if_main (int argc, char* argv[]) opt_filter_param = strtod(optarg, &endptr); endstr = optarg + strlen(optarg); if (endptr != endstr) { - fprintf(stderr,"Error setting --filter-param to %s\n", optarg); + sys_error(ERR_SEVERE,"Error setting --filter-param to %s\n", optarg); phm2if_usage(argv[0]); return (1); } @@ -204,7 +202,7 @@ phm2if_main (int argc, char* argv[]) opt_nsample = strtol(optarg, &endptr, 10); endstr = optarg + strlen(optarg); if (endptr != endstr) { - fprintf(stderr,"Error setting --nsample to %s\n", optarg); + sys_error(ERR_SEVERE,"Error setting --nsample to %s\n", optarg); phm2if_usage(argv[0]); return (1); } @@ -239,38 +237,31 @@ phm2if_main (int argc, char* argv[]) opt_nx = strtol(argv[optind+1], &endptr, 10); endstr = argv[optind+1] + strlen(argv[optind+1]); if (endptr != endstr) { - fprintf(stderr,"Error setting nx to %s\n", argv[optind+1]); + sys_error(ERR_SEVERE,"Error setting nx to %s\n", argv[optind+1]); phm2if_usage(argv[0]); return (1); } opt_ny = strtol(argv[optind+2], &endptr, 10); endstr = argv[optind+2] + strlen(argv[optind+2]); if (endptr != endstr) { - fprintf(stderr,"Error setting ny to %s\n", argv[optind+2]); + sys_error(ERR_SEVERE,"Error setting ny to %s\n", argv[optind+2]); phm2if_usage(argv[0]); return (1); } - snprintf(str, sizeof(str), "nx=%d, ny=%d, nsample=%d, ", opt_nx, opt_ny, opt_nsample); - if (opt_phmfilename[0]) { - strncat(str, "phantom=", sizeof(str)); - strncat(str, opt_phmfilename, sizeof(str)); - } - else if (opt_phmnum != -1) { - strncat(str, "phantom=", sizeof(str)); - strncat(str, name_of_phantom(opt_phmnum), sizeof(str)); - } + ostringstream oss; + oss << "nx=" << opt_nx << ", ny=" << opt_ny << ", nsample=" << opt_nsample << ", "; + if (opt_phmfilename.length()) + oss << "phantom=" << opt_phmfilename; + else if (opt_phmnum != -1) + oss << "phantom=" << name_of_phantom(opt_phmnum); else if (opt_filter != -1) { - strncat(str, "filter=", sizeof(str)); - strncat(str, name_of_filter(opt_filter), sizeof(str)); - strncat(str, " - ", sizeof(str)); - strncat(str, name_of_filter_domain(opt_filter_domain), sizeof(str)); - } - if (opt_desc[0]) { - strncat(str, ": ", sizeof(str)); - strncat(str, opt_desc, sizeof(str)); + oss << "filter=" << name_of_filter(opt_filter); + oss << " - " << name_of_filter_domain(opt_filter_domain); } - strncpy(opt_desc, str, sizeof(opt_desc)); + if (opt_desc.length()) + oss << ": " << opt_desc; + opt_desc = oss.str(); if (opt_verbose) cout << "Rasterize Phantom to Image" << endl << endl; @@ -363,7 +354,7 @@ phm2if_main (int argc, char* argv[]) { imGlobal->arrayDataWrite (); double calctime = timerProgram.timerEnd (); - imGlobal->labelAdd (Array2dFileLabel::L_HISTORY, opt_desc, calctime); + imGlobal->labelAdd (Array2dFileLabel::L_HISTORY, opt_desc.c_str(), calctime); imGlobal->fileClose (); if (opt_verbose) cout << "Time to rasterized phantom: " << calctime << " seconds" << endl;