关于MATLAB图像相关性分析的问题2008-10-2013:532829人阅读评论(2)收藏举报对于下边的程序,它的功能是对读入的图像随机的选取1000对像素点,然后对像素点进行相关性分析:clearall;a=imread('lenac.bmp');fork=1:1000i=fix(rand*250);j=fix(rand*250);x(k)=a(i+1,j+1);y(k)=a(i+2,j+1);endsum=0.0;fori=1:1000sum=sum+x(i);endex=sum/1000;sum=0.0;fori=1:1000sum=sum+y(i);endey=sum/1000;sum=0.0;fori=1:1000sum=sum+x(i)*x(i);endex2=sum/1000;sum=0.0;fori=1:1000sum=sum+y(i)*y(i);endey2=sum/1000;dx=ex2-ex*ex;dy=ey2-ey*ey;sum=0.0;fori=1:1000sum=sum+x(i)*y(i);endexy=sum/1000;cov=exy-ex*ey;xg=cov/(sqrt(dx)*sqrt(dy))这里有两个错误:首先系统会提示:sqrt使用错误,原因是dx和dy是unit8类型的。还有一个就是计算错误,你会发现不管你怎么加sum后边的值都一直会是255。上边两处错误的原因其实都是出在unit8类型上,默认情况下,系统读入图像后储存的类型就是unit8。此时我们直接对储存图像的数据进行操作,则后边程序中出现的所有的变量就都是unit8类型的,由于unit8类型的最大值就是255,所以sum会一直停留在255不变,不管你加多少次。要改正这种情况有以下两种方式:1、将存储读入图像的矩阵强制改成double型;2、存储矩阵不变,但是在后边的计算中将每个变量都强制改为double型;如下程序就是用第二种方式改正后的程序:clearall;a=imread('lenac.bmp');fork=1:1000i=fix(rand*250);j=fix(rand*250);x(k)=a(i+1,j+1);y(k)=a(i+2,j+1);endsum=0.0;sum=double(sum);fori=1:1000sum=sum+double(x(i));endex=sum/1000;sum=0.0;sum=double(sum);fori=1:1000sum=sum+double(y(i));endey=sum/1000;sum=0.0;sum=double(sum);fori=1:1000sum=sum+double(x(i))*double(x(i));endex2=sum/1000;sum=0.0;sum=double(sum);fori=1:1000sum=sum+double(y(i))*double(y(i));endey2=sum/1000;dx=ex2-ex*ex;dy=ey2-ey*ey;sum=0.0;sum=double(sum);fori=1:1000sum=sum+double(x(i))*double(y(i));endexy=sum/1000;cov=exy-ex*ey;dx=double(dx);dy=double(dy);xg=cov/(sqrt(dx)*sqrt(dy))分享到: