how can i solve problem in my code for adaptive median filter
Afficher commentaires plus anciens
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);
1 commentaire
Daleel Ahmed
le 28 Oct 2017
Modifié(e) : Walter Roberson
le 28 Oct 2017
Réponse acceptée
Plus de réponses (1)
Zeeshan Salam
le 24 Mar 2019
0 votes
can u send me working of this algorithm and flowchart
Catégories
En savoir plus sur Filter Analysis dans Centre d'aide et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!