冈萨雷斯 《数字图像处理matlab版》.源代码含修订版+源图像

function f = adpmedian（g Smax）
%ADPMEDIAN Perform adaptive median filtering.
%   F = ADPMEDIAN（G SMAX） performs adaptive median filtering of
%   image G.  The median filter starts at size 3-by-3 and iterates up
%   to size SMAX-by-SMAX. SMAX must be an odd integer greater than 1.

%   Copyright 2002-2004 R. C. Gonzalez R. E. Woods & S. L. Eddins
%   Digital Image Processing Using MATLAB Prentice-Hall 2004
%   $Revision: 1.5 $  $Date: 2003/11/21 14:19:05 $

% SMAX must be an odd positive integer greater than 1.
if （Smax <= 1） | （Smax/2 == round（Smax/2）） | （Smax ~= round（Smax））
   error（‘SMAX must be an odd integer > 1.‘）
end
[M N] = size（g）;

% Initial setup.
f = g;
f（:） = 0;
alreadyProcessed = false（size（g））;

% Begin filtering.
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‘）;
   
   processUsingLevelB = （zmed > zmin） & （zmax > zmed） & ...
       ~alreadyProcessed; 
   zB = （g > zmin） & （zmax > g）;
   outputZxy  = processUsingLevelB & zB;
   outputZmed = processUsingLevelB & ~zB;
   f（outputZxy） = g（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）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       5120  2004-12-16 13:36  dipum_1.1.3_mcodes_Read_Me_First.txt

     文件       1572  2004-12-16 13:36  dipum_1.1.3_mcodesadpmedian.m

     文件        651  2004-12-16 13:36  dipum_1.1.3_mcodesaverage.m

     文件       3235  2004-06-06 16:16  dipum_1.1.3_mcodesayesgauss （2）.m

     文件       3034  2004-12-16 13:36  dipum_1.1.3_mcodesayesgauss.m

     文件       3063  2004-12-16 13:36  dipum_1.1.3_mcodesound2eight.m

     文件       2477  2004-12-16 13:36  dipum_1.1.3_mcodesound2four.m

     文件       2224  2004-12-16 13:36  dipum_1.1.3_mcodesound2im.m

     文件       5776  2004-12-16 13:36  dipum_1.1.3_mcodesoundaries.m

     文件       3966  2004-12-16 14:36  dipum_1.1.3_mcodessubsamp （2）.m

     文件       3966  2004-12-16 13:36  dipum_1.1.3_mcodessubsamp.m

     文件        592  2004-12-16 13:36  dipum_1.1.3_mcodeschangeclass.m

     文件       2353  2004-12-16 13:36  dipum_1.1.3_mcodescolorgrad.m

     文件       2428  2004-12-16 13:36  dipum_1.1.3_mcodescolorseg.m

     文件       1177  2004-12-16 13:36  dipum_1.1.3_mcodescompare.m

     文件       1029  2010-05-26 16:33  dipum_1.1.3_mcodesconnectpoly.asv

     文件       1040  2004-12-16 13:36  dipum_1.1.3_mcodesconnectpoly.m

     文件       6128  2004-12-16 13:36  dipum_1.1.3_mcodesContents.m

     文件        799  2004-12-16 13:36  dipum_1.1.3_mcodesconwaylaws.m

     文件        977  2004-12-16 13:36  dipum_1.1.3_mcodescovmatrix.m

     文件        687  2004-12-16 13:36  dipum_1.1.3_mcodesdftcorr.m

     文件        882  2004-12-16 13:36  dipum_1.1.3_mcodesdftfilt.m

     文件        736  2004-12-16 13:36  dipum_1.1.3_mcodesdftuv.m

     文件       6622  2004-12-16 13:36  dipum_1.1.3_mcodesdiameter.m

     文件        874  2004-12-16 13:36  dipum_1.1.3_mcodesendpoints.m

     文件       1028  2004-12-16 13:36  dipum_1.1.3_mcodesentropy.m

     文件       8758  2004-12-16 13:36  dipum_1.1.3_mcodesfchcode.m

     文件       1454  2004-12-16 13:36  dipum_1.1.3_mcodesfrdescp.m

     文件        992  2004-12-16 13:36  dipum_1.1.3_mcodesfwtcompare.m

     文件        447  2004-12-16 13:36  dipum_1.1.3_mcodesgmean.m

............此处省略84个文件信息