• 大小: 25.13MB
    文件类型: .zip
    金币: 1
    下载: 0 次
    发布日期: 2023-07-04
  • 语言: 其他
  • 标签: 纹理分析  图片  

资源简介

本图片对需要进行纹理分析的朋友很有帮助,希望大家喜欢

资源截图

代码片段和文件信息

%计算Tamura纹理特征
close all
clear all
clc
t0=cputime;
I = imread(‘D1.gif‘);
[NxNy] = size(I);
Ng=256;
G=double(I);
%计算粗糙度(coarseness)
Sbest=zeros(NxNy);
E0h=zeros(NxNy);
E0v=zeros(NxNy);
E1h=zeros(NxNy);
E1v=zeros(NxNy);
E2h=zeros(NxNy);
E2v=zeros(NxNy);
E3h=zeros(NxNy);
E3v=zeros(NxNy);
E4h=zeros(NxNy);
E4v=zeros(NxNy);
E5h=zeros(NxNy);
E5v=zeros(NxNy);
flag=0;

for i=1:Nx
    for j=2:Ny
        E0h(ij)=G(ij)-G(ij-1);
    end
end
E0h=E0h/2;

for i=1:Nx-1
    for j=1:Ny
        E0v(ij)=G(ij)-G(i+1j);
    end
end
E0v=E0v/2;

%图片大小必须大于4*4才能计算E1h、E1v
if (Nx<4||Ny<4)
    flag=1;
end
if(flag==0)
    for i=1:Nx-1
        for j=3:Ny-1
            E1h(ij)=sum(sum(G(i:i+1j:j+1)))-sum(sum(G(i:i+1j-2:j-1)));
        end
    end
    for i=2:Nx-2
        for j=2:Ny
            E1v(ij)=sum(sum(G(i-1:ij-1:j)))-sum(sum(G(i+1:i+2j-1:j)));
        end
    end
    E1h=E1h/4;
    E1v=E1v/4;
end


%图片大小必须大于8*8才能计算E2h、E2v
if (Nx<8||Ny<8)
    flag=1;
end
if(flag==0)
    for i=2:Nx-2
        for j=5:Ny-3
            E2h(ij)=sum(sum(G(i-1:i+2j:j+3)))-sum(sum(G(i-1:i+2j-4:j-1)));
        end
    end
    for i=4:Nx-4
        for j=3:Ny-1
            E2v(ij)=sum(sum(G(i-3:ij-2:j+1)))-sum(sum(G(i+1:i+4j-2:j+1)));
        end
    end
    E2h=E2h/16;
    E2v=E2v/16;
end


%图片大小必须大于16*16才能计算E3h、E3v
if (Nx<16||Ny<16)
    flag=1
end
if(flag==0)
    for i=4:Nx-4
        for j=9:Ny-7
            E3h(ij)=sum(sum(G(i-3:i+4j:j+7)))-sum(sum(G(i-3:i+4j-8:j-1)));
        end
    end
    for i=8:Nx-8
        for j=5:Ny-3
            E3v(ij)=sum(sum(G(i-7:ij-4:j+3)))-sum(sum(G(i+1:i+8j-4:j+3)));
        end
    end
    E3h=E3h/64;
    E3v=E3v/64;
end


%图片大小必须大于32*32才能计算E4h、E4v 
if (Nx<32||Ny<32)
    flag=1;
end
if(flag==0)
    for i=8:Nx-8
        for j=17:Ny-15
            E4h(ij)=sum(sum(G(i-7:i+8j:j+15)))-sum(sum(G(i-7:i+8j-16:j-1)));
        end
    end
    for i=16:Nx-16
        for j=9:Ny-7
            E4v(ij)=sum(sum(G(i-15:ij-8:j+7)))-sum(sum(G(i+1:i+16j-8:j+7)));
        end
    end
    E4h=E4h/256;
    E4v=E4v/256;
end


%图片大小必须大于64*64才能计算E5h、E5v
if (Nx<64||Ny<64)
    flag=1;
end
if(flag==0)
    for i=16:Nx-16
        for j=33:Ny-31
            E5h(ij)=sum(sum(G(i-15:i+16j:j+31)))-sum(sum(G(i-15:i+16j-32:j-31)));
        end
    end
    for i=32:Nx-32
        for j=17:Ny-15
            E5v(ij)=sum(sum(G(i-31:ij-16:j+15)))-sum(sum(G(i+1:i+32j-16:j+15)));
        end
    end
    E5h=E5h/1024;
    E5v=E5v/1024;
end

for i=1:Nx
    for j=1:Ny
        [maxvindex]=max([E0h(ij)E0v(ij)E1h(ij)E1v(ij)E2h(ij)E2v(ij)E3h(ij)E3v(ij)E4h(ij)E4v(ij)E5h(ij)E5v(ij)]);
        k=floor((index+1)/2);
        Sbest(ij)=2.^k;
    end
end
Fcoarseness=sum(sum(Sbest))/(Nx*Ny);
%计算对比度
[countsgraylevels]=imhist(I);
PI=counts/(Nx*Ny);
averagevalue=sum(graylevels.*PI);
u4=sum((graylevels-repmat(averagevalue[2561])).^4.*PI);
standarddeviation=sum((graylevels-repmat(averagevalue[2561])).^2.*PI);
alpha4=u4/standarddeviation^2;
Fcontrast=

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-10-20 21:07  Brodatz纹理图像库\
     目录           0  2012-10-20 21:07  Brodatz纹理图像库\Brodatz纹理图像库1.1\
     文件      282020  2006-11-18 12:14  Brodatz纹理图像库\D1.gif
     文件      275134  2006-11-18 12:25  Brodatz纹理图像库\D10.gif
     文件      409259  2006-11-18 14:39  Brodatz纹理图像库\D100.gif
     文件      196183  2006-11-18 14:42  Brodatz纹理图像库\D101.gif
     文件      196649  2006-11-18 14:41  Brodatz纹理图像库\D102.gif
     文件      306666  2006-11-18 14:43  Brodatz纹理图像库\D103.gif
     文件      448402  2006-11-18 14:46  Brodatz纹理图像库\D105.gif
     文件      435449  2006-11-18 15:04  Brodatz纹理图像库\D106.gif
     文件      293287  2006-11-18 15:03  Brodatz纹理图像库\D107.gif
     文件      269763  2006-11-18 15:06  Brodatz纹理图像库\D108.gif
     文件      290104  2006-11-18 15:09  Brodatz纹理图像库\D109.gif
     文件      350058  2006-11-18 15:12  Brodatz纹理图像库\D110.gif
     文件      369494  2006-11-18 15:12  Brodatz纹理图像库\D111.gif
     文件      333496  2006-11-18 15:14  Brodatz纹理图像库\D112.gif
     文件      372412  2006-11-18 12:17  Brodatz纹理图像库\D2.gif
     文件      335990  2006-11-18 12:17  Brodatz纹理图像库\D3.gif
     文件      415253  2006-11-18 12:21  Brodatz纹理图像库\D4.gif
     文件      329992  2006-11-18 12:52  Brodatz纹理图像库\D41.gif
     文件      350969  2006-11-18 12:52  Brodatz纹理图像库\D42.gif
     文件      187161  2006-11-18 12:54  Brodatz纹理图像库\D43.gif
     文件      187154  2006-11-18 12:54  Brodatz纹理图像库\D44.gif
     文件      215492  2006-11-18 12:55  Brodatz纹理图像库\D45.gif
     文件      301623  2006-11-18 12:56  Brodatz纹理图像库\D46.gif
     文件      273741  2006-11-18 12:56  Brodatz纹理图像库\D47.gif
     文件      249355  2006-11-18 12:57  Brodatz纹理图像库\D48.gif
     文件      231847  2006-11-18 12:57  Brodatz纹理图像库\D49.gif
     文件      351599  2006-11-18 12:20  Brodatz纹理图像库\D5.gif
     文件      355489  2006-11-18 12:59  Brodatz纹理图像库\D50.gif
     文件      339311  2006-11-18 12:58  Brodatz纹理图像库\D51.gif
............此处省略55个文件信息

评论

共有 条评论