MATLAB编程经验.doc

MATLAB编程经验 由imread读进来的彩色图像一般是RGB格式，是一个三维矩阵，其每个元素都uint8类型，而matlab的基本矩阵运算和傅立叶变换都不支持uint8类型。因此必须将矩阵转化为double类型才能参与运算。 用到的函数： rgb2gray：RGB转256阶灰度 double()：强制类型转换 Matlab的imshow只能显示像素值为uint8的彩色或黑白图像，如果用double类型会显示出一片白色，如果用uint16类型会显示出一片黑色 用到的函数： uint8()：强制类型转换 利用傅立叶变换通过快速傅立叶变换进行矩阵运算 计算各子矩阵与模板矩阵的元素乘积和（即互相关）： abw=abs(ifft2(fft2(bw) .* fft2(rot90(a,2),m,n))); 所得矩阵abw的各元素abw(i,j)的值为以bw(i,j)为右下角元素的子矩阵与模板矩阵的元素乘积和 计算各子矩阵的元素平方和： bw2=abs(ifft2(fft2(bw.^2) .* fft2(ones(size(a)),m,n))); bwsqrt=sqrt(bw2); 所得矩阵bw2的各元素bw2(i,j)的值为以bw(i,j)为右下角元素的子矩阵的各元素的平??和的平方根 计算模板矩阵的元素平方和的平方根： asqrt=sqrt(sumsqr(a)); 计算各子矩阵与模板矩阵的相关系数： T=abw./bwsqrt./asqrt; 所得矩阵T的各元素T(i,j) 的值为以bw(i,j)为右下角元素的子矩阵与模板矩阵的相关系数 找出T中元素的最大值所在的位置(i,j)即为模板在大图中的右下角的位置 Tmax=max(max(T)); [r,c]=find(T=Tmax); 计算矩阵中以(i,j)为右下角，大小为N*N的子矩阵中各元素的和 T=real(ifft2(fft2(result(:,:,dimension)) .* fft2(ones(N,N),m,n))); 计算矩阵中以(i,j)为“中心”，大小为N*N的子矩阵中各元素的和 T=real(ifft2(fft2(result(:,:,dimension)) .* fft2(ones(N,N),m,n))); for i=1:1:m for j=1:1:n %result(i,j)=floor((T(i,j)+T(mod(i+N,m),j)+T(i,mod(j+N,n))+T(mod(i+N,m),mod(j+N,n)+1))/N^2/4); irange=i+N; jrange=j+N; if irangem irange=mod(irange,m); end if jrangen jrange=mod(jrange,n); end result(i,j,dimension)=floor((T(i,j)+T(irange,j)+T(i,jrange)+T(irange,jrange))/N^2/4); end end