for j=1:b subplot(a,b,k); imshow(mat2gray(c{i,j})); x(i,j)=std2(c{i,j}); y(i,j)=mean2(c{i,j}); k=k+1; end end function g=statt(f )%统计图像的一维功率谱 [m,n]=size(f); a=floor(m/2); b=floor(n/2); g1=fftshift(fft2(f)); for i=1:m d(i,j)=sqrt((a-i).^2+(b-j).^2); function [c,x,y]=compucell(f,a,b)%对图像进行分块，再对图像进行处理，并显示图像 [m,n]=size(f); r=floor(m/a); t=floor(n/b); p=f(1:a*r,1:b*t); c=mat2cell(p,ones(a,1)*r,ones(b,1)*t); k=1; for i=1:a for j=1:b; end end k=floor(sqrt(a.^2+b.^2)); for i=1:k c=find(di); g(i)=sum(abs(g1(c))); end l=linspace(1,k,k); figure plot(l,g); end