matlab做的模板匹配，简单有效

clear;clc;

d1=0;d2=0;d3=0;
Imag=imread（‘1.png‘）;%待分割图片
[X Y] = size（Imag）;
% 计算图像直方图
hist = imhist（Imag）;
p = hist/（X*Y）; % 各灰度概率
 
sumP = cumsum（p）;
sumQ = 1-sumP;
 
%将256个灰度作为256个分割阈值，分别计算各阈值下的概率密度函数
c0 = zeros（256256）;
c1 = zeros（256256）;
for i = 1:256
    for j = 1:i
        if sumP（i） > 0
            c0（ij） = p（j）/sumP（i）; %计算各个阈值下的前景概率密度函数
        else
            c0（ij） = 0;
        end
        for k = i+1:256
            if sumQ（i） > 0;
                c1（ik） = p（k）/sumQ（i）; %计算各个阈值下的背景概率密度函数
            else
                c1（ik） = 0;
            end
        end
    end 
end
 
%计算各个阈值下的前景和背景像素的累计熵
H0 = zeros（256256）;
H1 = zeros（256256）;
for i = 1:256
   for j = 1:i
       if c0（ij） ~=0
           H0（ij） =  - c0（ij）.*log10（c0（ij）