diff --git a/4_grains/GtThreshold.m b/4_grains/GtThreshold.m index 62a26066845843e6d4d6c92cd530bae1de0d8470..371e3d7d070a45d78d7898d3d5f0714b891757db 100644 --- a/4_grains/GtThreshold.m +++ b/4_grains/GtThreshold.m @@ -249,11 +249,13 @@ classdef GtThreshold < handle stdvol = std( vol(indexes) ); meanvol = mean( vol(indexes) ); - if (stdvol > meanvol) - threshold_val = obj.param.rec.percentile1/100 * meanvol; - else - threshold_val = meanvol - obj.param.rec.percentile2/100 * stdvol; - end + %if (stdvol > meanvol) + % threshold_val = obj.param.rec.percentile1/100 * meanvol; + %else + % threshold_val = meanvol - obj.param.rec.percentile2/100 * stdvol; + %end + % The above cases are not continuous! Use this instead: + threshold_val = meanvol - obj.param.rec.percentile2/100 * stdvol; center(1) = round(size(vol, 1) / 2); center(2) = round(size(vol, 2) / 2);