这个代码我也不是很懂原理，但是效果我觉得还可以。放进来方便我以后使用。
若代码侵权，请联系我。
function bw = multi_detect(I,num)
if isa(I,'double')
I = double(I);
end
range = [min(I(:)) max(I(:))];
smoothValue = 0.01*diff(range).^2;
%imguidedfilter导向滤波函数。方差比较大的区域，保持值不变，在*滑区域，使用临*像素*均（也就退化为普通均值滤波）
J = imguidedfilter(I,'DegreeOfSmoothing',smoothValue);
thresh = multithresh(J,2);
L = imquantize(J,thresh);
L = imfill(L);
bw = ones(size(L));
bw(L ==num) = 0;
放一个效果图：
