x=cp_image([4 84]);     x  =x(1:128,128+(1:256)); x  =round(8*(x+128)/256);
x_1=cp_image([4 84],1); x_1=x_1(1:128,1:128);     x_1=round(8*(x_1+128)/256);
x_2=cp_image([4 84],2); x_2=x_2(1:128,1:128);     x_2=round(8*(x_2+128)/256);
whos
[asm,con,cor,ent]    =cooc(x        ,[1 -1],[16 16]);
[asm0,con0,cor0,ent0]=cooc([x_1 x_2],[1 -1],[16 16]);
[asm1,con1,cor1,ent1]=cooc(x_1      ,[1 -1],[16 16]);
[asm2,con2,cor2,ent2]=cooc(x_2      ,[1 -1],[16 16]);
save cooc.mat

x=cp_image([4 84]);     x  =x(1:128,128+(1:256)); x  =round(7*x  /255);
x_1=cp_image([4 84],1); x_1=x_1(1:128,1:128);     x_1=round(7*x_1/255);
x_2=cp_image([4 84],2); x_2=x_2(1:128,1:128);     x_2=round(7*x_2/255);
whos
[asm,con,cor,ent]    =cooc(x        ,[1 -1],[16 16]);
[asm0,con0,cor0,ent0]=cooc([x_1 x_2],[1 -1],[16 16]);
[asm1,con1,cor1,ent1]=cooc(x_1      ,[1 -1],[16 16]);
[asm2,con2,cor2,ent2]=cooc(x_2      ,[1 -1],[16 16]);
save cooc2.mat
exit
