clear?all;?close?all;?clc;

img=double（imread（‘lena.jpg‘））;
imshow（img[]）;
[m?n]=size（img）;

img=sqrt（img）;??????%伽馬校正

%下面是求邊緣
fy=[-1?0?1];????????%定義豎直模板
fx=fy‘;?????????????%定義水平模板
Iy=imfilter（imgfy‘replicate‘）;????%豎直邊緣
Ix=imfilter（imgfx‘replicate‘）;????%水平邊緣
Ied=sqrt（Ix.^2+Iy.^2）;??????????????%邊緣強(qiáng)度
Iphase=Iy./Ix;??????????????%邊緣斜率，有些為inf-infnan，其中nan需要再處理一下


%下面是求cell
step=16;????????????????%step*step個(gè)像素作為一個(gè)單元
orient=9;???????????????%方向直方圖的方向個(gè)數(shù)
jiao=360/orient;????????%每個(gè)方向包含的角度數(shù)
Cell=cell（11）;??????????????%所有的角度直方圖cell是可以動(dòng)態(tài)增加的，所以先設(shè)了一個(gè)
ii=1;??????????????????????
jj=1;
for?i=1:step:m??????????%如果處理的m/step不是整數(shù)，最好是i=1:step:m-step
????ii=1;
????for?j=1:step:n??????%注釋同上
????????tmpx=Ix（i:i+step-1j:j+step-1）;
????????tmped=Ied（i:i+step-1j:j+ste

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????目錄???????????0??2014-02-17?23:48??hog?feature\
?????文件????????2695??2014-02-17?18:08??hog?feature\hog.m
?????文件???????10555??2014-02-17?18:05??hog?feature\lena.jpg