利用灰度共生矩阵对图像分类

%**************************************************************************
%                   图像检索——纹理特征
%基于共生矩阵纹理特征提取，d=1θ=0°45°90°135°共四个矩阵
%所用图像灰度级均为256
%参考《基于颜色空间和纹理特征的图像检索》
%function : T=Texture（Image） 
%Image    : 输入图像数据
%T        : 返回八维纹理特征行向量
%**************************************************************************
clc;
clear;
Image =imread（‘D:\毕业设计\纹理图\mosaic_02_train_03.bmp‘）;
figure（1）
subplot（221）
imshow（Image）
Gray=graycomatrix（Image）; 
[MN]=size（Gray）;


%--------------------------------------------------------------------------
%2.为了减少计算量，对原始图像灰度级压缩，将Gray量化成16级
%--------------------------------------------------------------------------
for i = 1:M
    for j = 1:N
        for n = 1:256/16
            if （n-1）*16==Gray（ij）&Gray（ij）==（n-1）*