darkTubes threshold

src/MATLAB/+RSF/darkTubes.m

@@ -13,13 +13,13 @@ imn(imn>0) = 0;
 imn = abs(imn);
 imp=imf;imp(imp<0)=0;
 
-bw = imbinarize(imn);
+bw = imbinarize(im);
 bw = bwareaopen(bw,1e3);
 bw = imdilate(bw,ones(4*thickness));
 
 imfm = fibermetric(imp);
 % NOTE as with imPreProcess, regional max first, then threshold
-imfm = imfm .* imregionalmax(imfm);
+imfm = imfm .* imregionalmax(imfm,6);
 imfm(~bw) = 0;
 
 imfm = RSF.clipRSF(imfm);