42 addUsageLine(
"Eliminate images whose variance is extremely large");
47 addParamsLine(
"[--confidence <conf=0.99>] : Remove an image if its variance is outside this confidence beyond sigma^2_0");
49 addParamsLine(
"[--minSigma2 <sigma2=0.01>] : Minimum variance");
60 <<
"Sigma2: " <<
sigma20 << std::endl
61 <<
"MinSigma2: " <<
minSigma2 << std::endl;
71 I().computeAvgStdev(avg,stddev);
72 double sigma2=stddev*stddev;
76 if (z>zalpha || sigma2<
minSigma2 || !std::isfinite(stddev) || fabs(avg)>
sigma20/9.0)
double getDoubleParam(const char *param, int arg=0)
double icdf_gauss(double p)
void processImage(const FileName &fnImg, const FileName &fnImgOut, const MDRow &rowIn, MDRow &rowOut)
Process image or volume.
int verbose
Verbosity level.
void setValue(MDLabel label, const T &d, bool addLabel=true)
int read(const FileName &name, DataMode datamode=DATA, size_t select_img=ALL_IMAGES, bool mapData=false, int mode=WRITE_READONLY)
void addUsageLine(const char *line, bool verbatim=false)
void finishProcessing()
Finish processing.
void addParamsLine(const String &line)