本文整理汇总了C++中ModelObject::SetupModelImage方法的典型用法代码示例。如果您正苦于以下问题:C++ ModelObject::SetupModelImage方法的具体用法?C++ ModelObject::SetupModelImage怎么用?C++ ModelObject::SetupModelImage使用的例子?那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类ModelObject
的用法示例。
在下文中一共展示了ModelObject::SetupModelImage方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: main
//.........这里部分代码省略.........
nColumns = options.nColumns;
nRows = options.nRows;
}
nPixels_tot = nColumns * nRows;
/* Read in PSF image, if supplied */
if (options.psfImagePresent) {
printf("Reading PSF image (\"%s\") ...\n", options.psfFileName.c_str());
psfPixels = ReadImageAsVector(options.psfFileName, &nColumns_psf, &nRows_psf);
if (psfPixels == NULL) {
fprintf(stderr, "\n*** ERROR: Unable to read PSF image file \"%s\"!\n\n",
options.psfFileName.c_str());
exit(-1);
}
nPixels_psf = nColumns_psf * nRows_psf;
printf("naxis1 [# pixels/row] = %d, naxis2 [# pixels/col] = %d; nPixels_tot = %d\n",
nColumns_psf, nRows_psf, nPixels_psf);
}
else
printf("* No PSF image supplied -- no image convolution will be done!\n");
if (! options.subsamplingFlag)
printf("* Pixel subsampling has been turned OFF.\n");
/* Set up the model object */
theModel = new ModelObject();
theModel->SetDebugLevel(options.debugLevel);
/* Add functions to the model object; also tells model object where function
sets start */
status = AddFunctions(theModel, functionList, functionBlockIndices, options.subsamplingFlag);
if (status < 0) {
fprintf(stderr, "*** ERROR: Failure in AddFunctions!\n\n");
exit(-1);
}
// Add PSF image vector, if present (needs to be added prior to image data, so that
// ModelObject can figure out proper internal model-image size
if (options.psfImagePresent) {
status = theModel->AddPSFVector(nPixels_psf, nColumns_psf, nRows_psf, psfPixels);
if (status < 0) {
fprintf(stderr, "*** ERROR: Failure in ModelObject::AddPSFVector!\n\n");
exit(-1);
}
}
/* Define the size of the requested model image */
theModel->SetupModelImage(nColumns, nRows);
theModel->PrintDescription();
// Set up parameter vector(s), now that we know how many total parameters
// there will be
nParamsTot = theModel->GetNParams();
printf("%d total parameters\n", nParamsTot);
if (nParamsTot != (int)parameterList.size()) {
fprintf(stderr, "*** ERROR: number of input parameters (%d) does not equal",
(int)parameterList.size());
fprintf(stderr, " required number of parameters for specified functions (%d)!\n\n",
nParamsTot);
exit(-1);
}
/* Copy parameters into C array and generate the model image */
paramsVect = (double *) calloc(nParamsTot, sizeof(double));
for (int i = 0; i < nParamsTot; i++)
paramsVect[i] = parameterList[i];
// Generate the image (including convolution, if requested), repeatedly
gettimeofday(&timer_start, NULL);
for (int ii = 0; ii < options.nIterations; ii++)
theModel->CreateModelImage(paramsVect);
gettimeofday(&timer_end, NULL);
microsecs = timer_end.tv_usec - timer_start.tv_usec;
time_elapsed = timer_end.tv_sec - timer_start.tv_sec + microsecs/1e6;
printf("\nELAPSED TIME FOR %d ITERATIONS: %.6f sec\n", options.nIterations, time_elapsed);
printf("Mean time per iteration = %.7f\n", time_elapsed/options.nIterations);
/* Save model image: */
printf("\nSaving output model image (\"%s\") ...\n", options.outputImageName.c_str());
vector<string> commentStrings;
SaveVectorAsImage(theModel->GetModelImageVector(), options.outputImageName,
nColumns, nRows, commentStrings);
printf("Done!\n\n");
// Free up memory
if (options.psfImagePresent)
free(psfPixels);
free(paramsVect);
delete theModel;
return 0;
}