clear?all;
close?all;
clc;

%img=double（imread（‘1.jpg‘））;
img1=imread（‘123.jpg‘）;
figure（1）
imshow（img1）;
img1?=?rgb2gray（img1）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%??灰度圖像取反??%%%%%%%%%%%
%Igray3?=?img1;
%d?=?find（img1）;
%img1（d（1:end））?=255;
%img?=?abs（img1-Igray3）;
%figure（2）
%imshow（img）;

img=double（img1）;

[m?n]=size（img）;
lastimg?=?img?-?img;
compimg?=?lastimg;
%w=4;
%sigma?=?10;
%court=0;
%for?sigma=[03510];
????%court=court+1;
%[x?y]=meshgrid（-w:w-w:w）;???%%%%%%%%返回的x，y是兩個矩陣，且兩個矩陣的行數和列數分別相等


H=cell（mn）;
sigmaseq?=?[pi/4pi/23*pi/4pi5*pi/46*pi/47*pi/420];
%test?=?0;
H1=cell（18）;
for?seq?=?1?:?8

????sigma?=?sigmaseq（seq）;
????r=-sqrt（2*pi*sigma）;
????w=4;
????[x?y]=meshgrid（-w:w-w:w）;
????Dxx?=?1/（2*pi*sigma^4）*（x.^2/sigma^2-1）*exp（-（y.^2

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-12-16?22:26??gaijin_hessian\
?????文件???????53848??2015-05-08?08:59??gaijin_hessian\123.jpg
?????文件????????2466??2015-05-12?11:42??gaijin_hessian\gaijin.m