灰度共生矩阵_灰度共生_灰度共生矩阵_matlab_matlab灰度共生_分割_

function [Q1]=compute(temp_m,graynum,dist) %%%%%灰度共生矩阵 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%graynum灰度级数 %dist距离 for i=1:graynum for j=1:graynum matH(i,j)=0;matV(i,j)=0;matRD(i,j)=0;matLD(i,j)=0; end end [M,N]=size(temp_m); double total;total=0; temp_m=double(temp_m); for m=1:M for n=1:N if ceil(temp_m(m,n)/(256/graynum))==0 newimage(m,n)=1; else newimage(m,n)=ceil(temp_m(m,n)/(256/graynum));%灰度量化 end end end %%% 0 for m=1:M for n=1:N-dist a=newimage(m,n); b=newimage(m,n+dist); matH(a,b)=matH(a,b)+1; matH(b,a)=matH(b,a)+1; end end % % %%%%%%%%%%%%%%% 90 % for m=1:16-dist % for n=1:16 % matV(newimage(m,n),newimage(m+dist,n))=matV(newimage(m,n),newimage(m+dist,n))+1; % matV(newimage(m+dist,n),newimage(m,n))=matV(newimage(m+dist,n),newimage(m,n))+1; % end % end % % %%%%%%%%%%%%%%%%%%% 135 % for m=1:16-dist % for n=16-dist % matLD(newimage(m,n),newimage(m+dist,n+dist))= matLD(newimage(m,n),newimage(m+dist,n+dist))+1; % matLD(newimage(m+dist,n+dist),newimage(m,n))= matLD(newimage(m+dist,n+dist),newimage(m,n))+1; % end % end % % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 45 % for m=1:M-dist % for n=1:N-dist % a=newimage(m,n); % b=newimage(m+dist,n+dist); % matH(a,b)=matH(a,b)+1; % matH(b,a)=matH(b,a)+1; % % matRD(newimage(m,n),newimage(m-dist+1,n+dist))=matRD(newimage(m,n),newimage(m-dist+1,n+dist))+1; % % matRD(newimage(m-dist+1,n+dist),newimage(m,n))=matRD(newimage(m-dist+1,n+dist),newimage(m,n))+1; % end % end %%%%%%%%%%%%%% for m=1:graynum for n=1:graynum total=total+matH(m,n); end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%5 for m=1:graynum for n=1:graynum mat_temp(m,n)=matH(m,n)/total; mat_temp(m,n)=double(mat_temp(m,n)); end end %%%%%%%%%%计算特征值 Q1=0;Q2=0; for m=1:graynum for n=1:graynum Q1=Q1+mat_temp(m,n)*mat_temp(m,n);%%能量 Q2=Q2+(m-n)*(m-n)*mat_temp(m,n);%%%%对比度 end end