自适应阈值分割的源程序

国外网上找的自适应阈值分割的源程序，测试过，好用。 可以下载

function f_threshold（file_name）
% Evan Ruzanski and Chiou-Wei Tsai
% UCSB 11/3/2003

clc

% Initialize variables
thresh = 77;
gray_range = 256;

% Read and display test image as given in Lee p. 712
im = imread（file_name）;
figure
imshow（im）; title（‘Original 2D image‘）;

% Perform binary spatial lifting
s = size（im）;
imz（1:s（1）1:s（2）1:gray_range） = 0; 
im = double（im）; % Typecast for arithmetic operations

for i = 1:s（1）
    for j = 1:s（2）
        if im（ij） > 1
            imz（ij1:im（ij）） = 1; % Make 3rd Dim binary
        end
    end
end


% 3D representation of lifting the image
figure（2）
x = 1:s（1）;
y = 1:s（2）;
[XY] = meshgrid（xy）;
colormap（gray）
rotate3d
mesh（XYim）; axis tight
title（‘Lifted representation of 2D image‘）


% Take 3D FFT o

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

     文件       3072  2003-10-30 13:05  f_threshold\278Aprojectimage.jpg

     文件       1212  2004-12-14 01:47  f_threshold\f_threshold.m

     目录          0  2008-06-12 17:24  f_threshold

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

                 4284                    3