Otsu Threshold algorithm function for blockproc

2 次查看（过去 30 天）
Issara Laosuwan
Issara Laosuwan 2016-5-25
评论： Issara Laosuwan 2016-5-25
Hello! I try to write an Otsu Threshold algorithm function to be used later along with the blockproc. Here is my code:
clc; clear; close all;
function out = graythresh(im)
if isstruct(im)
im = im.data;
end
if ~isempty(im)
I = im2uint8(im(:));
num_bins = 256;
counts = imhist(I,num_bins);
P = counts/sum(counts);
omega = cumsum(p);
mu = cusum(p.*(1:num_bins)');
mu_t = mu(end);
previous_state = warning('off','MATLAB:divideBYZero');
sigma_b_squared = (mu_t*omega-mu).^2./(omega.*(1-omega));
warning(previous_state);
maxval = max(sigma_b_squared);
isfinite_maxval = isfinite(maxval);
if isfinite_maxval
idx = mean(find(sigma_b_squared == maxval));
level = (idx-1)/(num_bins-1);
else
level = 0.0;
end
if nargout > 1
if isfinite_maxval
em = maxval/sum(p.*((1:num_bins).^2)')-mu_t^2);
else
em = 0;
end
end
out = (im>level*max(im(:)));
disp(['Threshold =',num2str(round(255*level))])
However, this code cannot be run successfully. The command window state that
Error: File: OTSUBLOCK.m Line: 3 Column: 1
Function definitions are not permitted in this context.
Please kindly give me advices
Thankyou!

请先登录，再回答此问题。

采纳的回答

Walter Roberson
Walter Roberson 2016-5-25
编辑：Walter Roberson 2016-5-25
You must delete the first line, the one that says
clc; clear; close all;
Also, you must either rename the file from OTSUBLOCK.m to graythresh.m or else you must change the "function out = graythresh(im)" to
function out = OUTSUBLOCK(im)

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Image Thresholding 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by