Index: branches/pap/ppSub/src/ppSubArguments.c =================================================================== --- branches/pap/ppSub/src/ppSubArguments.c (revision 23690) +++ branches/pap/ppSub/src/ppSubArguments.c (revision 23704) @@ -189,15 +189,7 @@ } -bool ppSubArgumentsSetup(int argc, char *argv[], pmConfig *config) -{ - // int argnum; // Argument Number +bool ppSubArguments(int argc, char *argv[], pmConfig *config, ppSubData *data) +{ assert(config); - - psMetadata *recipe = psMetadataLookupMetadata(NULL, config->recipes, PPSUB_RECIPE); // Recipe for ppSim - if (!recipe) { - psError(PS_ERR_UNEXPECTED_NULL, false, "Unable to find recipe %s", PPSUB_RECIPE); - return false; - } - psMetadata *arguments = config->arguments; // Command-line arguments @@ -212,45 +204,6 @@ psMetadataAddStr(arguments, PS_LIST_TAIL, "-kernel", 0, "Precalculated kernel to apply", NULL); psMetadataAddStr(arguments, PS_LIST_TAIL, "-stats", 0, "Statistics file", NULL); - psMetadataAddF32(arguments, PS_LIST_TAIL, "-region", 0, "Size of iso-kernel region", NAN); - psMetadataAddS32(arguments, PS_LIST_TAIL, "-size", 0, "Kernel half-size (pixels)", 0); - psMetadataAddS32(arguments, PS_LIST_TAIL, "-order", 0, "Spatial polynomial order", -1); - psMetadataAddStr(arguments, PS_LIST_TAIL, "-type", 0, - "Kernel type (ISIS|POIS|SPAM|FRIES|GUNK|RINGS)", NULL); - psMetadataAddF32(arguments, PS_LIST_TAIL, "-penalty", 0, "Penalty for wideness", NAN); - psMetadataAddStr(arguments, PS_LIST_TAIL, "-isis-widths", 0, - "ISIS Gaussian FWHMs (comma-separated)", NULL); - psMetadataAddStr(arguments, PS_LIST_TAIL, "-isis-orders", 0, - "ISIS polynomial orders (comma-separated)", NULL); - psMetadataAddS32(arguments, PS_LIST_TAIL, "-rings-order", 0, "RINGS polynomial order", -1); - psMetadataAddS32(arguments, PS_LIST_TAIL, "-inner", 0, "SPAM and FRIES inner radius", -1); - psMetadataAddS32(arguments, PS_LIST_TAIL, "-spam-binning", 0, "SPAM kernel binning", 0); - psMetadataAddF32(arguments, PS_LIST_TAIL, "-spacing", 0, "Typical stamp spacing (pixels)", NAN); - psMetadataAddS32(arguments, PS_LIST_TAIL, "-footprint", 0, "Stamp footprint half-size (pixels)", -1); - psMetadataAddF32(arguments, PS_LIST_TAIL, "-source-radius", 0, "Source matching radius (pixels)", NAN); - psMetadataAddS32(arguments, PS_LIST_TAIL, "-stride", 0, "Size of convolution patches (pixels)", -1); - psMetadataAddF32(arguments, PS_LIST_TAIL, "-threshold", 0, "Minimum threshold for stamps (ADU)", NAN); - psMetadataAddS32(arguments, PS_LIST_TAIL, "-iter", 0, "Number of rejection iterations", -1); - psMetadataAddF32(arguments, PS_LIST_TAIL, "-rej", 0, "Rejection thresold (sigma)", NAN); - psMetadataAddF32(arguments, PS_LIST_TAIL, "-sys", 0, "Relative systematic error in kernel", NAN); - psMetadataAddImageMask(arguments, PS_LIST_TAIL, "-mask-bad", 0, "Mask value for bad pixels", 0); - psMetadataAddImageMask(arguments, PS_LIST_TAIL, "-mask-poor", 0, "Mask value for poor pixels", 0); - psMetadataAddF32(arguments, PS_LIST_TAIL, "-poor-frac", 0, "Fraction of variance for poor pixels", NAN); - psMetadataAddF32(arguments, PS_LIST_TAIL, "-badfrac", 0, "Maximum fraction of bad pixels to accept", 1.0); - psMetadataAddBool(arguments, PS_LIST_TAIL, "-reverse", 0, "Reverse sense of subtraction?", false); - psMetadataAddBool(arguments, PS_LIST_TAIL, "-generate-mask", 0, "Generate mask if not supplied?", false); - psMetadataAddStr(arguments, PS_LIST_TAIL, "-stamps", 0, - "Stamps filename; file has x,y on each line", NULL); - psMetadataAddBool(arguments, PS_LIST_TAIL, "-opt", 0, - "Derive optimum parameters for ISIS kernels?", false); - psMetadataAddF32(arguments, PS_LIST_TAIL, "-opt-min", 0, "Minimum value for optimum kernel search", NAN); - psMetadataAddF32(arguments, PS_LIST_TAIL, "-opt-max", 0, "Minimum value for optimum kernel search", NAN); - psMetadataAddF32(arguments, PS_LIST_TAIL, "-opt-step", 0, "Step value for optimum kernel search", NAN); - psMetadataAddF32(arguments, PS_LIST_TAIL, "-opt-tol", 0, "Tolerance for optimum kernel search", NAN); - psMetadataAddS32(arguments, PS_LIST_TAIL, "-opt-order", 0, "Maximum order for optimum kernel search", -1); - psMetadataAddBool(arguments, PS_LIST_TAIL, "-dual", 0, "Dual convolution", false); - psMetadataAddBool(arguments, PS_LIST_TAIL, "-photometry", 0, "Perform photometry?", false); + psMetadataAddStr(arguments, PS_LIST_TAIL, "-stamps", 0, "Stamps filename; x,y on each line", NULL); psMetadataAddS32(arguments, PS_LIST_TAIL, "-threads", 0, "Number of threads", 0); - psMetadataAddS32(arguments, PS_LIST_TAIL, "-bin1", 0, "Binning factor for first level", 0); - psMetadataAddS32(arguments, PS_LIST_TAIL, "-bin2", 0, "Binning factor for second level", 0); psMetadataAddStr(arguments, PS_LIST_TAIL, "-dumpconfig", 0, "file to dump configuration to", NULL); psMetadataAddBool(arguments, PS_LIST_TAIL, "-visual", 0, "Show diagnostic plots", NULL); @@ -302,80 +255,16 @@ } - if (psMetadataLookupBool(NULL, arguments, "-photometry")) { - psMetadataAddBool(recipe, PS_LIST_TAIL, "PHOTOMETRY", PS_META_REPLACE, "Perform photometry?", true); - } - - return true; -} - -bool ppSubArgumentsParse(pmConfig *config) -{ - assert(config); - - psMetadata *arguments = config->arguments; // Command-line arguments - - valueArgStr(arguments, "-stats", "STATS", arguments); - valueArgStr(arguments, "-stamps", "STAMPS", arguments); - - psMetadata *recipe = psMetadataLookupMetadata(NULL, config->recipes, PPSUB_RECIPE); // Recipe for ppSim - if (!recipe) { - psError(PS_ERR_UNEXPECTED_NULL, false, "Unable to find recipe %s", PPSUB_RECIPE); - goto ERROR; - } - - VALUE_ARG_RECIPE_FLOAT("-region", "REGION.SIZE", F32); - VALUE_ARG_RECIPE_INT("-size", "KERNEL.SIZE", S32, 0); - VALUE_ARG_RECIPE_INT("-order", "SPATIAL.ORDER", S32, -1); - VALUE_ARG_RECIPE_FLOAT("-spacing", "STAMP.SPACING", F32); - VALUE_ARG_RECIPE_INT("-rings-order", "RINGS.ORDER", S32, -1); - VALUE_ARG_RECIPE_INT("-inner", "INNER", S32, -1); - VALUE_ARG_RECIPE_INT("-spam-binning", "SPAM.BINNING", S32, 0); - VALUE_ARG_RECIPE_INT("-footprint", "STAMP.FOOTPRINT", S32, -1); - VALUE_ARG_RECIPE_FLOAT("-source-radius", "SOURCE.RADIUS", F32); - VALUE_ARG_RECIPE_INT("-stride", "STRIDE", S32, -1); - VALUE_ARG_RECIPE_FLOAT("-threshold", "STAMP.THRESHOLD", F32); - VALUE_ARG_RECIPE_INT("-iter", "ITER", S32, -1); - VALUE_ARG_RECIPE_FLOAT("-rej", "REJ", F32); - VALUE_ARG_RECIPE_FLOAT("-sys", "SYS", F32); - VALUE_ARG_RECIPE_FLOAT("-badfrac", "BADFRAC", F32); - VALUE_ARG_RECIPE_FLOAT("-penalty", "PENALTY", F32); - VALUE_ARG_RECIPE_FLOAT("-poor-frac", "POOR.FRACTION", F32); - VALUE_ARG_RECIPE_INT("-bin1", "BIN1", S32, 0); - VALUE_ARG_RECIPE_INT("-bin2", "BIN2", S32, 0); - - valueArgRecipeStr(arguments, recipe, "-mask-in", "MASK.IN", recipe); - valueArgRecipeStr(arguments, recipe, "-mask-bad", "MASK.BAD", recipe); - valueArgRecipeStr(arguments, recipe, "-mask-poor", "MASK.POOR", recipe); - - vectorArgRecipe(arguments, "-isis-widths", recipe, "ISIS.WIDTHS", recipe, PS_TYPE_F32); - vectorArgRecipe(arguments, "-isis-orders", recipe, "ISIS.ORDERS", recipe, PS_TYPE_S32); - - psVector *widths = psMetadataLookupPtr(NULL, recipe, "ISIS.WIDTHS"); // ISIS Gaussian widths - psVector *orders = psMetadataLookupPtr(NULL, recipe, "ISIS.ORDERS"); // ISIS Polynomial orders - if (widths->n != orders->n) { - psError(PS_ERR_BAD_PARAMETER_SIZE, true, "Size of vectors for ISIS widths and orders do not match."); - goto ERROR; - } - - if (psMetadataLookupBool(NULL, arguments, "-opt") || psMetadataLookupBool(NULL, recipe, "OPTIMUM")) { - psMetadataAddBool(recipe, PS_LIST_TAIL, "OPTIMUM", PS_META_REPLACE, - "Derive optimum parameters for ISIS kernels?", true); - VALUE_ARG_RECIPE_FLOAT("-opt-min", "OPTIMUM.MIN", F32); - VALUE_ARG_RECIPE_FLOAT("-opt-max", "OPTIMUM.MAX", F32); - VALUE_ARG_RECIPE_FLOAT("-opt-step","OPTIMUM.STEP", F32); - VALUE_ARG_RECIPE_FLOAT("-opt-tol", "OPTIMUM.TOL", F32); - VALUE_ARG_RECIPE_INT("-opt-order", "OPTIMUM.ORDER", S32, -1); - } - - psMetadataAddBool(arguments, PS_LIST_TAIL, "REVERSE", 0, "Reverse sense of subtraction", - psMetadataLookupBool(NULL, arguments, "-reverse")); - psMetadataAddBool(recipe, PS_LIST_TAIL, "MASK.GENERATE", PS_META_REPLACE, "Generate mask if not supplied", - psMetadataLookupBool(NULL, arguments, "-generate-mask")); - psMetadataAddBool(recipe, PS_LIST_TAIL, "DUAL", PS_META_REPLACE, "Dual convolution?", - psMetadataLookupBool(NULL, arguments, "-dual")); - - // Need to update this because it could have been overwritten by the camera's own recipe - if (psMetadataLookupBool(NULL, arguments, "-photometry")) { - psMetadataAddBool(recipe, PS_LIST_TAIL, "PHOTOMETRY", PS_META_REPLACE, "Perform photometry?", true); + data->stamps = psMetadataLookupStr(NULL, arguments, "-stamps"); + + const char *statsName = psMetadataLookupStr(NULL, arguments, "-stats"); // Filename for statistics + if (statsName && strlen(statsName) > 0) { + psString resolved = pmConfigConvertFilename(statsName, config, true, true); // Resolved filename + data->statsFile = fopen(resolved, "w"); + if (!data->statsFile) { + psError(PS_ERR_IO, true, "Unable to open statistics file %s for writing.\n", resolved); + psFree(resolved); + return false; + } + psFree(resolved); } @@ -384,19 +273,5 @@ } - // Translate the kernel type - psString type = psMetadataLookupStr(NULL, arguments, "-type"); // Name of kernel type - if (!type || strlen(type) == 0) { - type = psMetadataLookupStr(NULL, recipe, "KERNEL.TYPE"); - if (!type || strlen(type) == 0) { - psError(PS_ERR_BAD_PARAMETER_VALUE, false, "Unable to find KERNEL.TYPE specified."); - goto ERROR; - } - } - psMetadataAddStr(recipe, PS_LIST_TAIL, "KERNEL.TYPE", PS_META_REPLACE, "Type of kernel", type); - - psTrace("ppSub", 1, "Done reading command-line arguments\n"); - - // XXX move this to ppSubArguments - int threads = psMetadataLookupS32(NULL, config->arguments, "-threads"); // Number of threads + int threads = psMetadataLookupS32(NULL, arguments, "-threads"); // Number of threads if (threads > 0) { if (!psThreadPoolInit(threads)) { @@ -407,6 +282,10 @@ return true; - -ERROR: - return false; -} +} + + + + psTrace("ppSub", 1, "Done reading command-line arguments\n"); + + return true; +}