这个代码我也不是很懂原理,但是效果我觉得还可以。放进来方便我以后使用。
若代码侵权,请联系我。
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;
放一个效果图:

浙公网安备 33010602011771号