改进后的海森矩阵算法（Hessian ）代码

一种改进后的海森矩阵算法代码，在医学图像处理中效果比较好，可以识别较细的血管。

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