生成灰度共生矩阵 chengcu.doc

《生成灰度共生矩阵 chengcu.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、functionout=cooccurrence(input,dir,dist,symmetric);input=round(input);[rc]=size(input);min_intensity=min(min(input));max_intensity=max(max(input));out=zeros(max_intensity-min_intensity+1);if(dir==0)dir_x=0;dir_y=1;endif(dir==1)dir_x=1;dir_y=1;endif(dir==2)dir_x=1;dir_y=0;endif(d

2、ir==3)dir_x=1;dir_y=-1;endif(dir==4)dir_x=0;dir_y=-1;endif(dir==5)dir_x=-1;dir_y=-1;endif(dir==6)dir_x=-1;dir_y=0;endif(dir==7)dir_x=-1;dir_y=1;enddir_x=dir_x*dist;dir_y=dir_y*dist;out_ind_x=0;out_ind_y=0;forintensity1=min_intensity:max_intensityout_ind_x=out_ind_x+1;out_ind_y=0

3、;[ind_x1ind_y1]=find(input==intensity1);ind_x1=ind_x1+dir_x;ind_y1=ind_y1+dir_y;forintensity2=min_intensity:max_intensityout_ind_y=out_ind_y+1;[ind_x2ind_y2]=find(input==intensity2);count=0;fori=1:size(ind_x1,1)forj=1:size(ind_x2,1)if((ind_x1(i)==ind_x2(j))&&(ind_y1(i)==ind_y2(j

4、)))count=count+1;endendendout(out_ind_x,out_ind_y)=count;endendif(symmetric)if(dir<4)dir=dir+4;elsedir=mod(dir,4);endout=out+cooccurrence(input,dir,dist,0);end