int
unrrdu_gammaMain(int argc, const char **argv, const char *me,
hestParm *hparm) {
hestOpt *opt = NULL;
char *out, *err;
Nrrd *nin, *nout;
double min, max, gamma;
airArray *mop;
int pret, blind8BitRange;
NrrdRange *range;
hestOptAdd(&opt, "g,gamma", "gamma", airTypeDouble, 1, 1, &gamma, NULL,
"gamma > 1.0 brightens; gamma < 1.0 darkens. "
"Negative gammas invert values (like in xv). ");
hestOptAdd(&opt, "min,minimum", "value", airTypeDouble, 1, 1, &min, "nan",
"Value to implicitly map to 0.0 prior to calling pow(). "
"Defaults to lowest value found in input nrrd.");
hestOptAdd(&opt, "max,maximum", "value", airTypeDouble, 1, 1, &max, "nan",
"Value to implicitly map to 1.0 prior to calling pow(). "
"Defaults to highest value found in input nrrd.");
hestOptAdd(&opt, "blind8", "bool", airTypeBool, 1, 1, &blind8BitRange,
nrrdStateBlind8BitRange ? "true" : "false",
"Whether to know the range of 8-bit data blindly "
"(uchar is always [0,255], signed char is [-128,127]).");
OPT_ADD_NIN(nin, "input nrrd");
OPT_ADD_NOUT(out, "output nrrd");
mop = airMopNew();
airMopAdd(mop, opt, (airMopper)hestOptFree, airMopAlways);
USAGE(_unrrdu_gammaInfoL);
PARSE();
airMopAdd(mop, opt, (airMopper)hestParseFree, airMopAlways);
nout = nrrdNew();
airMopAdd(mop, nout, (airMopper)nrrdNuke, airMopAlways);
range = nrrdRangeNew(min, max);
airMopAdd(mop, range, (airMopper)nrrdRangeNix, airMopAlways);
nrrdRangeSafeSet(range, nin, blind8BitRange);
if (nrrdArithGamma(nout, nin, range, gamma)) {
airMopAdd(mop, err = biffGetDone(NRRD), airFree, airMopAlways);
fprintf(stderr, "%s: error doing gamma:\n%s", me, err);
airMopError(mop);
return 1;
}
SAVE(out, nout, NULL);
airMopOkay(mop);
return 0;
}
int
baneGkms_scatMain(int argc, char **argv, char *me, hestParm *hparm) {
hestOpt *opt = NULL;
char *out[2], *perr, err[BIFF_STRLEN];
Nrrd *hvol, *nvgRaw, *nvhRaw, *nvgQuant, *nvhQuant;
NrrdRange *vgRange, *vhRange;
airArray *mop;
int pret, E;
double gamma;
hestOptAdd(&opt, "g", "gamma", airTypeDouble, 1, 1, &gamma, "1.0",
"gamma used to brighten/darken scatterplots. "
"gamma > 1.0 brightens; gamma < 1.0 darkens. "
"Negative gammas invert values (like in xv). ");
hestOptAdd(&opt, "i", "hvolIn", airTypeOther, 1, 1, &hvol, NULL,
"input histogram volume (from \"gkms hvol\")",
NULL, NULL, nrrdHestNrrd);
hestOptAdd(&opt, "o", "vgOut vhOut", airTypeString, 2, 2, out, NULL,
"Filenames to use for two output scatterplots, (gradient "
"magnitude versus value, and 2nd derivative versus value); "
"can use PGM or PNG format");
mop = airMopNew();
airMopAdd(mop, opt, (airMopper)hestOptFree, airMopAlways);
USAGE(_baneGkms_scatInfoL);
PARSE();
airMopAdd(mop, opt, (airMopper)hestParseFree, airMopAlways);
nvgRaw = nrrdNew();
nvhRaw = nrrdNew();
nvgQuant = nrrdNew();
nvhQuant = nrrdNew();
airMopAdd(mop, nvgRaw, (airMopper)nrrdNuke, airMopAlways);
airMopAdd(mop, nvhRaw, (airMopper)nrrdNuke, airMopAlways);
airMopAdd(mop, nvgQuant, (airMopper)nrrdNuke, airMopAlways);
airMopAdd(mop, nvhQuant, (airMopper)nrrdNuke, airMopAlways);
if (baneRawScatterplots(nvgRaw, nvhRaw, hvol, AIR_TRUE)) {
sprintf(err, "%s: trouble creating raw scatterplots", me);
biffAdd(BANE, err); airMopError(mop); return 1;
}
vgRange = nrrdRangeNewSet(nvgRaw, nrrdBlind8BitRangeFalse);
vhRange = nrrdRangeNewSet(nvhRaw, nrrdBlind8BitRangeFalse);
airMopAdd(mop, vgRange, (airMopper)nrrdRangeNix, airMopAlways);
airMopAdd(mop, vhRange, (airMopper)nrrdRangeNix, airMopAlways);
E = 0;
if (!E) E |= nrrdArithGamma(nvgRaw, nvgRaw, vgRange, gamma);
if (!E) E |= nrrdArithGamma(nvhRaw, nvhRaw, vhRange, gamma);
if (!E) E |= nrrdQuantize(nvgQuant, nvgRaw, vgRange, 8);
if (!E) E |= nrrdQuantize(nvhQuant, nvhRaw, vhRange, 8);
if (E) {
sprintf(err, "%s: trouble doing gamma or quantization", me);
biffMove(BANE, err, NRRD); airMopError(mop); return 1;
}
if (!E) E |= nrrdSave(out[0], nvgQuant, NULL);
if (!E) E |= nrrdSave(out[1], nvhQuant, NULL);
if (E) {
sprintf(err, "%s: trouble saving scatterplot images", me);
biffMove(BANE, err, NRRD); airMopError(mop); return 1;
}
airMopOkay(mop);
return 0;
}