小波矩特征提取代码.doc

《小波矩特征提取代码.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、小波矩特征提取matlab代码分类：ImageRecognitionMatlab2010-12-1010:00122人阅读评论(0)收藏举报这是我上研究生时写的小波矩特征提取代码： %新归一化方法小波矩特征提取----------------------------------------------------------F=imread('a1.bmp');F=im2bw(F);F=imresize(F,[128128]);%求取最上点fori=1:128   forj=1:128       if(F(i,j)==1)           ytop=i;        

2、   xtop=j;           break;       end   end   if(F(i,j)==1)       break;   endend%求取最下点fori=1:128   forj=1:128       if(F(129-i,j)==1)           ybottom=129-i;           xbottom=j;           break;       end   end   if(F(129-i,j)==1)       break;   endend%求取最左点fori=1:128   forj=1:128       

3、if(F(j,i)==1)           yleft=j;           xleft=i;           break;       end   end   if(F(j,i)==1)       break;   endend%求取最右点fori=1:128   forj=1:128       if(F(j,129-i)==1)           yright=j;           xright=129-i;           break;       end   end   if(F(j,129-i)==1)       break;   end

4、end%求取中心点x0=(xright-xleft)/2+xleft;y0=(ybottom-ytop)/2+ytop;x0=round(x0);y0=round(y0);%图像平移F=double(F);[M,N]=size(F);F1=zeros(M,N);M0=M/2;N0=N/2;fori=1:M forj=1:N       ifF(i,j)==1          F1(i+M0-y0,j+N0-x0)=1;       end   endend%figure,imshow(F1);%图像缩放max=0;fori=1:128   forj=1:128       

5、if(F(i,j)==1)           d=sqrt((i-y0)^2+(j-x0)^2);           if(max

6、=round((500-M)/2);x1=round((500-N)/2);fori=1:M   forj=1:N       if(F2(i,j)==1)           f3(y1+i,x1+j)=1;       end   endend%figure,imshow(f3);%图像从笛卡儿坐标转换为极坐标------------------------------------------------------------%角度间隔为2*pi/(128*128),128个像素长设为图像的单位圆半径%f1矩阵里面放着对应极半径和角度的值f1=zeros(128,163

7、84);%直角坐标与极坐标建立起一一对应的关系          fori=1:128   forj=1:16384       a=j*2*pi/16384.0;       r=i;       y=round(r*sin(a));       x=round(r*cos(a));      if(f3(250+x,250+y)==1)           f1(i,j)=1;       end   endendF3=zeros(512,512);%将极坐标转换后的图像显示出来fori=1:128