机器学习初探（手写数字识别）HOG图片

这里我们讲一下使用HOG的方法进行手写数字识别：

首先把

代码分享出来：

hog1.m

function B = hog1(A)
%A是28*28的
B=[];
[x,y] = size(A);
%外圈补0
A(:,y+) = ;
A(x+,:) = ;
for i=:x
    deltax(:,i)=A(:,i+)-A(:,i);
end
for i=:y
    deltay(i,:)=A(i+,:)-A(i,:);
end
for i = :
    for j = :
        Px=deltax(i*-:i*+,j*-:j*+);
        Py=deltay(i*-:i*+,j*-:j*+);
        B=[B;hog2(Py,Px)];
    end
end

hog2.m

function res = hog2(Px,Py);
%这里Px和Py是x方向和Y方向的偏导矩阵
%返回是一个7*1的矩阵
res=zeros(,);
[pxx,pxy]=size(Px);
for k=:pxx*pxy
    len = sqrt(Px(k)^+Py(k)^);
    if atan2(Py(k),Px(k))>=-pi && atan2(Py(k),Px(k))<*pi/-pi
        res()=res()+len;
    elseif atan2(Py(k),Px(k))>=*pi/-pi && atan2(Py(k),Px(k))<*pi/-pi
        res()=res()+len;
    elseif atan2(Py(k),Px(k))>=*pi/-pi && atan2(Py(k),Px(k))<*pi/-pi
        res()=res()+len;
    elseif atan2(Py(k),Px(k))>=*pi/-pi && atan2(Py(k),Px(k))<*pi/-pi
        res()=res()+len;
    elseif atan2(Py(k),Px(k))>=*pi/-pi && atan2(Py(k),Px(k))<*pi/-pi
        res()=res()+len;
    elseif atan2(Py(k),Px(k))>=*pi/-pi && atan2(Py(k),Px(k))<*pi/-pi
        res()=res()+len;
    elseif atan2(Py(k),Px(k))>=*pi/-pi && atan2(Py(k),Px(k))<*pi/-pi
        res()=res()+len;
    end
end

hog1得到的是一个1183*1的矩阵。

http://blog.csdn.net/hqh45/article/details/44228715