matlab-多尺度分割

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;