X-Git-Url: http://git.kpe.io/?p=ctsim.git;a=blobdiff_plain;f=src%2Fphm2if.cpp;h=a6a7e260aa50ac75fe7ae9a00204c839c5d0ae3d;hp=3e82eefaadebc2802b5a430352af2f4ef191e13a;hb=2f3d6e2580db607105bb072b13e4aff453ae4495;hpb=b5857e74e5735455b5ef11cbea5044ae7b2e8a0d diff --git a/src/phm2if.cpp b/src/phm2if.cpp index 3e82eef..a6a7e26 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.15 2000/06/26 21:15:24 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 @@ -61,6 +61,7 @@ phm2if_usage (const char *program) cout << " ny Number of pixels Y-axis" << endl; cout << " --phantom Phantom to use for projection" << endl; cout << " herman Herman head phantom" << endl; + cout << " bherman Bordered Herman head phantom" << endl; cout << " rowland Rowland head phantom" << endl; cout << " browland Bordered Rowland head phantom" << endl; cout << " unitpulse Unit pulse phantom" << endl; @@ -107,18 +108,18 @@ phm2if_main (int argc, char* argv[]) Phantom phm; int opt_nx = 0, opt_ny = 0; int opt_nsample = 1; - int opt_phmnum = -1; - FilterType opt_filter = static_cast(-1); - DomainType opt_filter_domain = D_SPATIAL; + string optPhmName = Phantom::PHM_HERMAN_STR; + string optFilterName; + string optDomainName = SignalFilter::DOMAIN_SPATIAL_STR; 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 +130,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; @@ -141,22 +140,13 @@ phm2if_main (int argc, char* argv[]) switch (c) { case O_PHANTOM: - if ((opt_phmnum = opt_set_phantom(optarg)) < 0) { - phm2if_usage(argv[0]); - return (1); - } + optPhmName = optarg; break; case O_PHMFILE: - strncpy(opt_phmfilename, optarg, sizeof(opt_phmfilename)); - phm.createFromFile(opt_phmfilename); -#ifdef HAVE_MPI - if (mpiWorld.getRank() == 0) - cerr << "Can't use phantom from file in MPI mode" << endl; - return (1); -#endif + opt_phmFileName = optarg; break; case O_VERBOSE: - opt_verbose = 1; + opt_verbose = true; break; case O_DEBUG: opt_debug = 1; @@ -168,25 +158,19 @@ phm2if_main (int argc, char* argv[]) } break; case O_FILTER: - if ((opt_filter = opt_set_filter(optarg)) < 0) { - phm2if_usage (argv[0]); - return (1); - } + optFilterName = optarg; break; case O_FILTER_DOMAIN: - if ((opt_filter_domain = opt_set_filter_domain(optarg)) < 0) { - phm2if_usage (argv[0]); - return (1); - } + optDomainName = optarg; 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 +179,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 +188,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); } @@ -225,8 +209,8 @@ phm2if_main (int argc, char* argv[]) } } - if (phm.nPElem() == 0 && opt_phmnum == -1 && opt_filter == -1) { - cerr << "No phantom defined" << endl; + if (optPhmName == "" && optFilterName == "" && opt_phmFileName == "") { + cerr << "No phantom defined" << endl << endl; phm2if_usage(argv[0]); return (1); } @@ -239,39 +223,49 @@ 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)); + ostringstream oss; + oss << "phm2if: nx=" << opt_nx << ", ny=" << opt_ny << ", nsample=" << opt_nsample << ", "; + if (opt_phmFileName != "") + oss << "phantom=" << opt_phmFileName; + else if (optPhmName != "") + oss << "phantom=" << optPhmName; + else if (optFilterName != "") { + oss << "filter=" << optFilterName << " - " << optDomainName; } - else if (opt_phmnum != -1) { - strncat(str, "phantom=", sizeof(str)); - strncat(str, name_of_phantom(opt_phmnum), sizeof(str)); - } - 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 != "") + oss << ": " << opt_desc; + opt_desc = oss.str(); + + if (optPhmName != "") { + phm.createFromPhantom (optPhmName.c_str()); + if (phm.fail()) { + cout << phm.failMessage() << endl << endl; + phm2if_usage(argv[0]); + return (1); + } } - if (opt_desc[0]) { - strncat(str, ": ", sizeof(str)); - strncat(str, opt_desc, sizeof(str)); + + if (opt_phmFileName != "") { + phm.createFromFile(opt_phmFileName.c_str()); +#ifdef HAVE_MPI + if (mpiWorld.getRank() == 0) + cerr << "Can't use phantom from file in MPI mode" << endl; + return (1); +#endif } - strncpy(opt_desc, str, sizeof(opt_desc)); - + if (opt_verbose) cout << "Rasterize Phantom to Image" << endl << endl; #ifdef HAVE_MPI @@ -279,43 +273,34 @@ phm2if_main (int argc, char* argv[]) #endif #ifdef HAVE_MPI - TimerMPI timerBcast (mpiWorld.getComm()); + TimerCollectiveMPI timerBcast (mpiWorld.getComm()); + mpiWorld.BcastString (optPhmName); mpiWorld.getComm().Bcast (&opt_verbose, 1, MPI::INT, 0); mpiWorld.getComm().Bcast (&opt_debug, 1, MPI::INT, 0); mpiWorld.getComm().Bcast (&opt_trace, 1, MPI::INT, 0); mpiWorld.getComm().Bcast (&opt_nx, 1, MPI::INT, 0); mpiWorld.getComm().Bcast (&opt_ny, 1, MPI::INT, 0); mpiWorld.getComm().Bcast (&opt_nsample, 1, MPI::INT, 0); - mpiWorld.getComm().Bcast (&opt_phmnum, 1, MPI::INT, 0); - mpiWorld.getComm().Bcast (&opt_filter, 1, MPI::INT, 0); - mpiWorld.getComm().Bcast (&opt_filter_domain, 1, MPI::INT, 0); mpiWorld.getComm().Bcast (&opt_filter_param, 1, MPI::DOUBLE, 0); mpiWorld.getComm().Bcast (&opt_filter_bw, 1, MPI::DOUBLE, 0); + mpiWorld.BcastString (optFilterName); + mpiWorld.BcastString (optDomainName); + if (opt_verbose) timerBcast.timerEndAndReport ("Time to broadcast variables"); mpiWorld.setTotalWorkUnits (opt_nx); + if (mpiWorld.getRank() > 0 && optPhmName != "") + phm.createFromPhantom (optPhmName.c_str()); + if (mpiWorld.getRank() == 0) { - imGlobal = new ImageFile (opt_outfile, opt_nx, opt_ny); - imGlobal->fileCreate(); + imGlobal = new ImageFile (opt_nx, opt_ny); } imLocal = new ImageFile (opt_nx, opt_ny); #else - imGlobal = new ImageFile (opt_outfile, opt_nx, opt_ny); - imGlobal->fileCreate (); -#endif - - if (opt_phmnum >= 0) - phm.create (opt_phmnum); - -#ifdef HAVE_MPI - else { - if (mpiWorld.getRank() == 0) - cerr << "phmnum < 0" << endl; - exit(1); - } + imGlobal = new ImageFile (opt_nx, opt_ny); #endif ImageFileArray v; @@ -327,17 +312,17 @@ phm2if_main (int argc, char* argv[]) if (mpiWorld.getRank() == 0) { v[opt_nx/2][opt_ny/2] = 1.; } - } else if (opt_filter != -1) { + } else if (optFilterName != "") { if (mpiWorld.getRank() == 0) { - image_filter_response (*imGlobal, opt_filter_domain, opt_filter_bw, opt_filter, opt_filter_param, opt_trace); + imGlobal->filterResponse (optDomainName.c_str(), opt_filter_bw, optFilterName.c_str(), opt_filter_param); } } else { - TimerMPI timerRasterize (mpiWorld.getComm()); - phm_to_imagefile (phm, *imLocal, mpiWorld.getMyStartWorkUnit(), mpiWorld.getMyLocalWorkUnits(), opt_nsample, opt_trace); + TimerCollectiveMPI timerRasterize (mpiWorld.getComm()); + phm.convertToImagefile (*imLocal, opt_nsample, opt_trace, mpiWorld.getMyStartWorkUnit(), mpiWorld.getMyLocalWorkUnits()); if (opt_verbose) timerRasterize.timerEndAndReport ("Time to rasterize phantom"); - TimerMPI timerGather (mpiWorld.getComm()); + TimerCollectiveMPI timerGather (mpiWorld.getComm()); mpi_gather_image (mpiWorld, imGlobal, imLocal, opt_debug); if (opt_verbose) timerGather.timerEndAndReport ("Time to gather image"); @@ -346,14 +331,14 @@ phm2if_main (int argc, char* argv[]) v = imGlobal->getArray (); if (phm.getComposition() == P_UNIT_PULSE) { v[opt_nx/2][opt_ny/2] = 1.; - } else if (opt_filter != -1) { - image_filter_response (*imGlobal, opt_filter_domain, opt_filter_bw, opt_filter, opt_filter_param, opt_trace); + } else if (optFilterName != "") { + imGlobal->filterResponse (optDomainName.c_str(), opt_filter_bw, optFilterName.c_str(), opt_filter_param); } else { #if HAVE_SGP if (opt_trace >= TRACE_PHM) phm.show(); #endif - phm_to_imagefile (phm, *imGlobal, 0, opt_nx, opt_nsample, opt_trace); + phm.convertToImagefile (*imGlobal, opt_nsample, opt_trace); } #endif @@ -361,10 +346,9 @@ phm2if_main (int argc, char* argv[]) if (mpiWorld.getRank() == 0) #endif { - imGlobal->arrayDataWrite (); double calctime = timerProgram.timerEnd (); - imGlobal->labelAdd (Array2dFileLabel::L_HISTORY, opt_desc, calctime); - imGlobal->fileClose (); + imGlobal->labelAdd (Array2dFileLabel::L_HISTORY, opt_desc.c_str(), calctime); + imGlobal->fileWrite (opt_outfile); if (opt_verbose) cout << "Time to rasterized phantom: " << calctime << " seconds" << endl; @@ -376,7 +360,7 @@ phm2if_main (int argc, char* argv[]) scanf ("%d", &nscale); printf ("Enter minimum and maximum densities (min, max): "); scanf ("%lf %lf", &dmin, &dmax); - image_display_scale (*imGlobal, nscale, dmin, dmax); + imGlobal->displayScaling (nscale, dmin, dmax); } } @@ -414,6 +398,16 @@ void mpi_gather_image (MPIWorld& mpiWorld, ImageFile* imGlobal, ImageFile* imLoc int main (int argc, char* argv[]) { - return (phm2if_main(argc, argv)); + int retval = 1; + + try { + retval = phm2if_main(argc, argv); + } catch (exception e) { + cerr << "Exception: " << e.what() << endl; + } catch (...) { + cerr << "Unknown exception" << endl; + } + + return (retval); } #endif