matlab图像处理,特征提取.docVIP

下载本文档

10
0
约3.42千字
约 4页
2017-08-26 发布于重庆
举报
版权申诉

matlab图像处理,特征提取.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

matlab图像处理,特征提取

MATLAB特征提取代码 for i=1:26??? f=strcat(D:\bishe\,num2str(i));??? image=strcat(f,.jpg);??? PS=imread(image); ??? PS=imresize(PS,[300,300],bilinear);%归一化大小??? PS=rgb2gray(PS);??? [m,n]=size(PS);?????????????????????? %测量图像尺寸参数??? GP=zeros(1,256);???????????????????? %预创建存放灰度出现概率的向量??? for k=0:255??????? GP(k+1)=length(find(PS==k))/(m*n);? %计算每级灰度出现的概率，将其存入GP中相应位置??? end%直方图均衡化??? S1=zeros(1,256);?? for i=1:256?????? for j=1:i?????????? S1(i)=GP(j)+S1(i);????????????? %计算Sk?????? end?? end?? S2=round((S1*256)+0.5);?????????????? %将Sk归到相近级的灰度%图像均衡化?? f=PS;?? for i=0:255?????? f(find(PS==i))=S2(i+1);???????? %将各个像素归一化后的灰度值赋给这个像素?? end?? figure,imshow(f);%边缘检测?? f=edge(f,canny,0.25);?? imshow(f);%二值法锐化图像?? f=double(f);?? [x,y]=gradient(f);?? g=sqrt(x.*x+y.*y);?? i=find(g=0.5);?? g(i)=256;?? j=find(g0.5);?? g(j)=0;?? imshow(g);?? title(二值法锐化图像);%中值滤波g=medfilt2(g); g=dither(g); imshow(g);%提取面积，矩形度，圆形度，拉伸度特征?? %g=im2bw(g);?? [x,y]=size(g);?? BW = bwperim(g,8); % 检测边缘跟踪，用于计算周长 %检测垂直方向连读的周长像素点%?? P1=0;?? P2=0;?? Ny=0; % 记录垂直方向连续周长像素点的个数?? for i=1:x????? for j=1:y????????? if (BW(i,j)0)????????????? P2=j;????????????? if ((P2-P1)==1) % 判断是否为垂直方向连续的周长像素点????????????????? Ny=Ny+1;????????????? end???????????? P1=P2;????????? end????? end?? end%检测水平方向连读的周长像素点?? P1=0;?? P2=0;?? Nx=0; % 记录水平方向连续周长像素点的个数?? for j=1:y?????? for i=1:x?????????? if (BW(i,j)0)?????????????? P2=i;????????????? if ((P2-P1)==1) % 判断是否为水平方向连续的周长像素点?????????????????? Nx=Nx+1;????????????? end????????????? P1=P2;?????????? end?????? end?? end ?? SN=sum(sum(BW)); % 计算周长像素点的总数?? Nd=SN-Nx-Ny; % 计算奇数码的链码数目?? H=max(sum(g)); % 计算目标的高度 ?? W=max(sum(g)); % 图象g经矩阵转置后，计算宽度?? L=sqrt(2)*Nd+Nx+Ny; % 计算周长%====形态特征值计算===%?? A=bwarea(g); % 计算目标的面积?? R=A/(H*W); % 计算矩形度?? E=min(H,W)/max(H,W); % 计算伸长度?? temp1=[A,R,E];%提取不变矩特征?? [M,N]=size(g);?? [x,y]=meshgrid(1:N,1:M);?? x=x(:);?? y=y(:);?? g=g(:);?? m.m00=sum(g);?? if(m.m00==0)????? m.m00=e