图像去噪基于matlab自适应中值滤波图像去噪【含Matlab 1844期】.zip

function f = adaptive_median_filter (g, Smax) %% 判断邻域是否合理 if (Smax <= 1) || (Smax/2 == round(Smax/2)) || (Smax ~= round(Smax)) error ('SMAX must be an odd integer > 1.') end %% f = g; f(:) = 0; %% 标记是否已处理过 alreadyProcessed = false (size(g)); %% 开始自适应滤波 for k = 3:2:Smax zmin = ordfilt2(g, 1, ones(k, k),'symmetric'); zmax = ordfilt2(g, k * k, ones(k, k), 'symmetric'); zmed = medfilt2(g, [k k], 'symmetric'); %% 判断是否进入进程B processUsingLevelB = (zmed > zmin) & (zmax > zmed) & ~alreadyProcessed; %% 若g不是脉冲，保留原值 zB = (g > zmin) & (zmax > g); outputZxy = processUsingLevelB & zB; %% 若是脉冲，用Zmed替换 outputZmed = processUsingLevelB & ~zB; f (outputZxy) = g(outputZxy); f (outputZmed) = zmed(outputZmed); %% 已处理记录 alreadyProcessed = alreadyProcessed | processUsingLevelB; %% 是否退出 if all (alreadyProcessed (:)) break; end end %% 大于窗口尺寸后，Zxy替换成Zmed输出 f (~alreadyProcessed) = zmed (~alreadyProcessed); end