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）;??????????????%邊緣強度
Iphase=Iy./Ix;??????????????%邊緣斜率，有些為inf-infnan，其中nan需要再處理一下


%下面是求cell
step=16;????????????????%step*step個像素作為一個單元
orient=9;???????????????%方向直方圖的方向個數
jiao=360/orient;????????%每個方向包含的角度數
Cell=cell（11）;??????????????%所有的角度直方圖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-1j:j+step-1）;
????????tmped=Ied（i:i+step-1j:j+ste

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????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