how can i solve problem in my code for adaptive median filter
2 次查看(过去 30 天)
显示 更早的评论
function f= dd(r,Smax)
if (Smax <= 1) | (Smax/2 == round(Smax/2)) | (Smax ~= round(Smax))
error('SMAX must be an odd integer > 1.')
end
[M, N] = size(r);
% Initial setup.
f = r;
f(:) = 0;
alreadyProcessed = false(size(r));
%Begin filtering.
for k = 3:2:Smax
zmin = ordfilt2(r, 1, ones(k, k), 'symmetric');
zmax = ordfilt2(r, k * k, ones(k, k), 'symmetric');
zmed = medfilt2(r, [k k], 'symmetric');
processUsingLevelB = (zmed > zmin) & (zmax > zmed) & ...
~alreadyProcessed;
zB = (r > zmin) & (zmax > r);
outputZxy = processUsingLevelB & zB;
outputZmed = processUsingLevelB & ~zB;
f(outputZxy) = r(outputZxy);
f(outputZmed) = zmed(outputZmed);
alreadyProcessed = alreadyProcessed | processUsingLevelB;
if all(alreadyProcessed(:))
break;
end
end
% Output zmed for any remaining unprocessed pixels. Note that this
% zmed was computed using a window of size Smax-by-Smax, which is
% the final value of k in the loop.
f(~alreadyProcessed) = zmed(~alreadyProcessed);
采纳的回答
Walter Roberson
2017-10-28
You appear to be attempting to pass in characters for the first input to your function, the one named "r" here. That input must be numeric for this code to work.
0 个评论
更多回答(1 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Numerical Integration and Differential Equations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!