【人脸识别】 HOG特征提取人脸识别【含Matlab源码 641期】.zip

function featureVec= hog(img) img=double(img); step=8; %8*8个像素作为一个cell [m1 n1]=size(img); %改变图像尺寸为step的最近整数倍 img=imresize(img,[floor(m1/step)*step,floor(n1/step)*step],'nearest'); [m n]=size(img); % 2、%伽马校正 img=sqrt(img); % 3、求梯度和方向 fy=[-1 0 1]; %定义竖直模板 fx=fy'; %定义水平模板 Iy=imfilter(img,fy,'replicate'); %竖直梯度 Ix=imfilter(img,fx,'replicate'); %水平梯度 Ied=sqrt(Ix.^2+Iy.^2); %梯度幅值 Iphase=Iy./Ix; %边缘斜率，有些为inf,-inf,nan，其中nan需要再处理一下 the=atan(Iphase)*180/3.14159; %求梯度角度 for i=1:m for j=1:n if(Ix(i,j)>=0&&Iy(i,j)>=0) %第一象限 the(i,j)=the(i,j); elseif(Ix(i,j)<=0&&Iy(i,j)>=0) %第二象限 the(i,j)=the(i,j)+180; elseif(Ix(i,j)<=0&&Iy(i,j)<=0) %第三象限 the(i,j)=the(i,j)+180; elseif(Ix(i,j)>=0&&Iy(i,j)<=0) %第四象限 the(i,j)=the(i,j)+360; end if isnan(the(i,j))==1 %0/0会得到nan，如果像素是nan，重设为0 the(i,j)=0; end end end the=the+0.000001; %防止角度为 % 4、划分cell，求cell的直方图( 1 cell = 8*8 pixel ) %下面是求cell step=8; %step*step个像素作为一个cell orient=9; %方向直方图的方向个数 jiao=360/orient; %每个方向包含的角度数 Cell=cell(1,1); %所有的角度直方图,cell是可以动态增加的，所以先设了一个 ii=1; jj=1; for i=1:step:m ii=1; for j=1:step:n Hist1(1:orient)=0; for p=1:step for q=1:step %梯度方向直方图 Hist1(ceil(the(i+p-1,j+q-1)/jiao))=Hist1(ceil(the(i+p-1,j+q-1)/jiao))+Ied(i+p-1,j+q-1); end end Cell{ii,jj}=Hist1; %放入Cell中 ii=ii+1; end jj=jj+1; end % 5、划分block，求block的特征值,使用重叠方式( 1 block = 2*2 cell ) [m n]=size(Cell); feature=cell(1,(m-1)*(n-1)); for i=1:m-1 for j=1:n-1 block=[]; block=[Cell{i,j}(:)' Cell{i,j+1}(:)' Cell{i+1,j}(:)' Cell{i+1,j+1}(:)']; block=block./sum(block); %归一化 feature{(i-1)*(n-1)+j}=block; end end % 6、图像的HOG特征值 [m n]=size(feature); l=2*2*orient; featureVec=zeros(1,n*l); for i=1:n featureVec((i-1)*l+1:i*l)=feature{i}(:); end % [m n]=size(img); % img=sqrt(img); %伽马校正 % %下面是求边缘 % fy=[-1 0 1]; %定义竖直模板 % fx=fy'; %定义水平模板 % Iy=imfilter(img,fy,'replicate'); %竖直边缘 % Ix=imfilter(img,fx,'replicate'); %水平边缘 % Ied=sqrt(Ix.^2+Iy.^2); %边缘强度 % Iphase=Iy./Ix; %边缘斜率，有些为inf,-inf,nan，其中nan需要再处理一下 % % % %下面是求cell % step=16; %step*step个像素作为一个单元 % orient=9; %方向直方图的方向个数 % jiao=360/orient; %每个方向包含的角度数 % Cell=cell(1,1); %所有的角度直方图,cell是可以动态增加的，所以先设了一个 % ii=1; % jj=1; % for i=1:step:m %如果处理的m/step不是整数，最好是i=1:step:m-step % ii=1; % for j=1:step:n %注释同上 % tmpx=Ix(i:i+step-1,j:j+step-1); % tmped=Ied(i:i+step-1,j:j+step-1); % tmped=tmped/sum(sum(tmped)); %局部边缘强度归一化 % tmpphase=Iphase(i:i+step-1,j:j+step-1); % Hist=zeros(1,orient); %当前step*step像素块统计角度直方图,就是cell % for p=1:step % for q=1:step % if isnan(tmpphase(p,q))==1 %0/0会得到nan，如果像素是nan，重设为0 % tmpphase(p,q)=0; % end % ang=atan(tmpphase(p,q)); %atan求的是[-90 90]度之间 % ang=mod(ang*180/pi,360); %全部变正，-90变270 % if tmpx(p,q)<0 %根据x方向确定真正的角度 % if ang<90 %如果是第一象限 % ang=ang+180; %移到第三象限 % end % if ang>270 %如果是第四象限 % ang=ang-180; %移到第二象限 % end % end % ang=ang+0.0000001; %防止ang为0 % Hist(ceil(ang/jiao))=Hist(ceil(ang/jiao))+tmped(p,q); %ceil向上取整，使用边缘强度加权 % end % end % Hist=Hist/sum(Hist); %方向直方图归一化 % Cell{ii,jj}=Hist; %放入Cell中 % ii=ii+1; %针对Cell的y坐标循环变量 % end % jj=jj+1; %针对Cell的x坐标循环变量 % end % % %下面是求feature,2*2个cell合成一个block,没有显式的求block % [m n]=size(Cell); % feature=cell(1,(m-1)*(n-1)); % for i=1:m-1 % for j=1:n-1 % f=[]; % f=[f Cell{i,j}(:)' Cell{i,j+1}(:)' Cell{i+1,j}(:)' Cell{i+1,j+1}(:)']; % feature{(i-1)*(n-1)+j}=f; % end % end % % % % 到此结束，feature即为所求 % % 下面是为了显示而写的 % l=length(feature); % f=[]; % for i=1:l % f=[f;feature{i}(:)']; % end % % figure % % mesh(f)