clear?all;



%{

[filename?pathname?filter]=uigetfile（‘*.*‘‘圖像文件‘）;
if??filter==0
????return?
end?
img=fullfile（pathnamefilename）;
a=imread（img）;

%}
a=[0255;?0255];?
[hw]=size（a）
%[hwc]=size（a）;

%水平方向相關性分析??
for?k=1:1000
????i=fix（rand*（w-1））+1;%i的取值范圍為[0w-1]
????j=fix（rand*w）+1;???%j的取值為[0w]
????x（k）=a（ij）;%
????y（k）=a（i+1j）;
end
sum=0.0;
sum=double（sum）;
for?i=1:1000
????sum=sum+double（x（i））;
end
ex=sum/1000;%求均值E（x）
sum=0.0;
sum=double（sum）;
for?i=1:1000
????sum=sum+double（y（i））;
end
ey=sum/1000;%求均值E（y）
sum=0.0;
sum=double（sum）;
for?i=1:1000
????sum=sum+（double（x（i））-double（ex））*（double（x（i））-double（ex））;
end
dx=double（sum/1000）;
sum=0.0
sum=double（sum）;
for?i=1:1000
????sum=sum+（double（y（i））-double（ey））*（double（y（i））-double（ey））;
end
dy=double（sum/1000）;
sum=0.0
sum=double（sum）;
for?i=1:1000
????sum=sum+double（double（x（i））-double（ex））*double（double（y（i））-double（ey））
end
cov=double（sum/1000）;
rela=cov/（sqrt（dx）*sqrt（dy））;
rela

%繪制相關性分析圖
figure
plot（xy‘.‘‘MarkerSize‘4）;%marker是圖上畫上點的地方標上符號，size指標號的大小
xlabel（‘像素值（xy）‘）;
ylabel（‘像素值（x+1y）‘）;

%{
?%垂直方向相關性分析
for?k=1:1000
????i=fix（rand*h）+1;
????j=fix（rand*（h-1））+1;
????x（k）=a（ij）;
????y（k）=a（ij+1）;
end
sum=0.0;
sum=double（sum）;
for?i=1:1000
????sum=sum+double（x（i））;
end

ex=sum/1000;
sum=0.0;
sum=double（sum）;
for?i=1:1000
????sum=sum+double（y（i））;
end
ey=sum/1000;
sum=0.0;
sum