图像二维信息熵的计算

clear all; clc;
%%读取输入图像                                 
im_T = imread（‘Data/HL/t20.bmp‘）;                                          
%%输入图像RGB转转化为三通道图像
image = im_T;
image_ycbcr = rgb2ycbcr（image）;                                              
image_y = image_ycbcr（: : 1）;
image_cb = image_ycbcr（: : 2）;
image_cr = image_ycbcr（: : 3）;
%初始化
cnt = 0;                                                                   %记录区域块位置
sumend = 0;                                                                %二维信息熵的值
maxcount=0;                                                                %最高频次
center = 0;                                                                %区域块的中心点
average = 0;                                                               %邻域灰度均值
patch_size =5;