matlab特征提取代码（Matlab feature extraction code）.docVIP

matlab特征提取代码（Matlab feature extraction code） 我= 1:26 F = strcat（想：\ bishe \”，num2str（我））； 图像= strcat（F，“.jpg”）； PS = imread（图像）； PS = imresize（PS，[ 300300 ]，bilinear’）；%归一化大小 PS = rgb2gray（PS）； [ M ]，N =大小（PS）；%测量图像尺寸参数 GP =零（1256）；%预创建存放灰度出现概率的向量 K = 0:255 GP（K + 1）=长度（找到（PS = = k））/（m×n）；%计算每级灰度出现的概率，将其存入GP中相应位置 结束 %直方图均衡化 S1 =零（1256）； 我为1:256 J = 1：我 S1（我）= GP（J）+ S1（我）；%计算SK 结束 结束 S2 =圆（（S1×256）+ 0.5）；%将SK归到相近级的灰度 %图像均衡化 F＝PS； 我= 0:255 F（找到（PS = =我））= S2（i + 1）；%将各个像素归一化后的灰度值赋给这个像素 结束 图，imshow（f）； %边缘检测 F =边缘（F，canny ，0.25）； imshow（f）； %二值法锐化图像 F =双（F）； 【X，Y ] =梯度（F）； G = SQRT（X * X + Y Y）； 我发现（G＞＝0.5）； G（我）= 256； J =找到（g＜0.5）； G（j）= 0； imshow（G）； 标题（“二值法锐化图像”）； %中值滤波 G = medfilt2（G）； G =抖动（G）； imshow（G）； %提取面积，矩形度，圆形度，拉伸度特征 % G = im2bw（G）； 【X，Y ] =大小（G）； BW = bwperim（G，8）；%检测边缘跟踪，用于计算周长 检测垂直方向连读的周长像素点% % P1 = 0； P2 = 0； 纽约= 0；%记录垂直方向连续周长像素点的个数 我= 1 x J = 1：Y 如果（BW（i，j） 0） P2 = J； 如果（（p2-p1）= = 1）%判断是否为垂直方向连续的周长像素点 纽约纽约+ 1； 结束 P1 = P2； 结束 结束 结束 %检测水平方向连读的周长像素点 P1 = 0； P2 = 0； NX = 0；%记录水平方向连续周长像素点的个数 J = 1：Y 我= 1 x 如果（BW（i，j） 0） P2 =我； 如果（（p2-p1）= = 1）%判断是否为水平方向连续的周长像素点 NX = NX + 1； 结束 P1 = P2； 结束 结束 结束 Sn =金额（金额（BW））；%计算周长像素点的总数 ND = Sn NX纽约；%计算奇数码的链码数目 H=max (sum (g));% calculated the target height W=max (sum (G)); g% by image matrix transpose, calculating width L=sqrt (2) *Nd+Nx+Ny%; calculate the perimeter %==== form eigenvalue calculation ===% A=bwarea (g);% calculated the target area R=A/ (H*W);% rectangular degree E=min (H, W) /max (H, W);% calculated length Temp1=[A, R, E]; % extraction of invariant feature [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=eps; End M.m10=sum (x.*g); M.m01=sum (y.*g); M.m11=sum (x.*y.*g); M.m20=sum (x.^2.*g); M.m02=sum (y.^2.*g); M.m30=sum (x.^3.*g); M.m03=sum (y.^3.*g); M.m12=sum (x.*y.^2.*g); M.m21=sum (x.^2.*y.*g); Xbar=m.m10/m.m00; Ybar=m.m01/m.m00; E.eta11= (m.m11-ybar*m.m10) /m.m00^2; E.eta20= (m.m20-xbar*m.m10) /m.m00^2; E.eta02= (m.m02-ybar*m.m01) /m.m00^2;