77 for (
size_t index=0; index<
itsNdata; index++) {
113 for (
size_t index=0; index<
itsNdata; index++) {
127 for (
size_t index=0; index<
itsNdata; index++) {
185 for (
size_t index=0; index<ndata; index++)
187 pixVal =cellData[index];
191 sigma += pixVal*pixVal;
206 throw Elements::Exception() <<
"Can not compute meaningful stats with statNData=" << statNData <<
"!";
208 mean /= (double)statNData;
209 sigma = sigma/statNData - mean*mean;
210 sigma = sigma>0.0 ?
sqrt(sigma):0.0;
221 for (
size_t index=0; index<ndata; index++)
223 pixVal =cellData[index];
224 if (pixVal>=lcut && pixVal <= hcut)
227 sigma += pixVal*pixVal;
235 throw Elements::Exception() <<
"Can not compute meaningful stats with statNData=" << statNData <<
"!";
238 mean /= (double)statNData;
239 sigma = sigma/statNData - mean*mean;
240 sigma = sigma>0.0 ?
sqrt(sigma):0.0;
262 for (
size_t index=0; index<ndata; index++)
264 pixVal = cellData[index];
265 weightVal = cellWeight[index];
266 if (weightVal < weightThresh && pixVal > -
BIG)
269 sigma += pixVal*pixVal;
270 weightMean += weightVal;
271 weightSigma += weightVal*weightVal;
288 mean /= (double)statNData;
289 sigma = sigma/statNData - mean*mean;
290 sigma = sigma>0.0 ?
sqrt(sigma):0.0;
291 weightMean /= (double)statNData;
292 weightSigma = weightSigma/statNData - weightMean*weightMean;
293 weightSigma = weightSigma>0.0 ?
sqrt(weightSigma):0.0;
310 for (
size_t index=0; index<ndata; index++)
312 pixVal = cellData[index];
313 weightVal = cellWeight[index];
314 if (weightVal < weightThresh && pixVal <=hcut && pixVal>=lcut)
317 sigma += pixVal*pixVal;
319 if (weightVal<=weightHcut && weightVal>= weightLcut)
321 weightMean += weightVal;
322 weightSigma += weightVal*weightVal;
331 throw Elements::Exception() <<
"Can not compute meaningful stats with statNData=" << statNData;
333 mean /= (double)statNData;
334 sigma = sigma/statNData - mean*mean;
335 sigma = sigma>0.0 ?
sqrt(sigma):0.0;
340 throw Elements::Exception() <<
"Can not compute meaningful stats with statNData=" << statNWeight;
342 weightMean /= (double)statNWeight;
343 weightSigma = weightSigma/statNWeight - weightMean*weightMean;
344 weightSigma = weightSigma>0.0 ?
sqrt(weightSigma):0.0;