数字图像处理实验——滤波，皮肤毛发去除，傅里叶变换

数字图像处理实验——滤波，皮肤毛发去除，傅里叶变换

function thresh = otsu（datapix_num）
%计算分割阈值

img_var = zeros（2561）;
for i=1:256
    w0 = sum（sum（data<=i-1））./pix_num;
    w1 = 1-w0;
    u0 = sum（sum（data.*double（data<=i-1）））./（w0*pix_num）;
    u1 = sum（sum（data.*double（data>i-1）））./（w1*pix_num）;
    img_var（i） = w0.*w1.*（（u0-u1）.^2）;
end

[~I] = max（img_var）;
thresh = I-1;

end

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

     文件     958374  2008-04-17 16:53  dip experiment\datapicture\1.bmp

     文件     214326  2008-04-17 16:53  dip experiment\datapicture\2.bmp

     文件     889074  2008-04-17 16:53  dip experiment\datapicture\3.bmp

     文件    1243280  2009-05-02 09:23  dip experiment\datapicture\CHEN2-7.BMP

     文件    1203344  2009-05-02 11:02  dip experiment\datapicture\CHEN2-8.BMP

     文件      66616  2009-06-03 12:36  dip experiment\datapicture\google-256.bmp

     文件      66614  2010-05-07 10:35  dip experiment\datapicture\grid-fft.bmp

     文件      66616  2001-11-09 01:40  dip experiment\datapicture\grid.bmp

     文件      66616  2006-11-05 14:04  dip experiment\datapicture\lena.bmp

     文件      66616  2009-03-28 17:47  dip experiment\datapicture\rect-45度.bmp

     文件      66616  2009-03-28 17:46  dip experiment\datapicture\rect.bmp

     文件      74520  2009-06-03 12:37  dip experiment\datapicture\参考图.bmp

     文件      17464  2009-06-03 12:37  dip experiment\datapicture\基准图.bmp

     文件     940856  2007-05-21 20:37  dip experiment\datapicture\皮肤镜图像.bmp

     文件     940854  2009-06-11 11:40  dip experiment\datapicture\结果.bmp

     文件        341  2012-05-15 22:30  dip experiment\otsu.m

     文件       1940  2012-05-16 10:19  dip experiment\task1.m

     文件       2450  2012-05-16 10:19  dip experiment\task2.m

     文件       1353  2012-05-16 10:15  dip experiment\task3.m

     文件        392  2012-05-15 09:30  dip experiment\test.m

     目录          0  2012-05-15 17:41  dip experiment\datapicture

     目录          0  2012-05-16 10:13  dip experiment

----------- ---------  ---------- -----  ----

              6888262                    22