声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1426|回复: 3

[综合讨论] 关于harris算子提取角点的高斯滤波问题

[复制链接]
发表于 2008-6-17 20:05 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
clear;
filename='F:\MATLAB7.0\work\left.jpg';
X= imread(filename);
%
读取图像
Info=imfinfo(filename);
if Info.BitDepth>8

f=rgb2gray(X);

end
%
% fx = [5 0 -5;8 0 -8;5 0 -5];
%
高斯函数一阶微分,x方向(用于改进的Harris角点提取算法)
ori_im=double(f)/255;
%unit8
转化为64为双精度double64
fx = [-2 -1 0 1 2];
% x
方向梯度算子(用于Harris角点提取算法)
Ix = filter2(fx,ori_im);
% x
方向滤波
% fy = [5 8 5;0 0 0;-5 -8 -5];
%
高斯函数一阶微分,y方向(用于改进的Harris角点提取算法)
fy = [-2;-1;0;1;2];
% y
方向梯度算子(用于Harris角点提取算法)
Iy = filter2(fy,ori_im);
% y
方向滤波
Ix2 = Ix.^2;
Iy2 = Iy.^2;
Ixy = Ix.*Iy;
clear Ix;
clear Iy;

h= fspecial('gaussian',[7 7],2);
%
产生7*7的高斯窗函数,sigma=2

Ix2 = filter2(h,Ix2);
Iy2 = filter2(h,Iy2);
Ixy = filter2(h,Ixy);

height = size(ori_im,1);
width = size(ori_im,2);
result = zeros(height,width);
%
纪录角点位置,角点处值为1

R = zeros(height,width);

Rmax = 0;

%
图像中最大的R
for i = 1:height

for j = 1:width


M = [Ix2(i,j) Ixy(i,j);Ixy(i,j) Iy2(i,j)];
%
自相关矩阵

R(i,j) = det(M)-0.06*(trace(M))^2;
%
计算R

if R(i,j) > Rmax


Rmax = R(i,j);


end;


end;

end;

cnt = 0;
for i = 2:height-1

for j = 2:width-1


%
进行非极大抑制,窗口大小3*3

if R(i,j) > 0.01*Rmax && R(i,j) > R(i-1,j-1) && R(i,j) > R(i-1,j) && R(i,j) > R(i-1,j+1) && R(i,j) > R(i,j-1) && R(i,j) > R(i,j+1) && R(i,j) > R(i+1,j-1) && R(i,j) > R(i+1,j) && R(i,j) > R(i+1,j+1)


result(i,j) = 1;


cnt = cnt+1;


end;


end;

end;

[posc, posr] = find(result == 1);
cnt
%
角点个数
imshow(ori_im)
hold on;
plot(posr,posc,'r+');

上面的代码是小弟在论坛里面找到的,能实现Harris角点提取,请问大侠们:x方向梯度算子fx 为什么取[-2 -1 0 1 2]呢?在X,Y方向滤波后为什么还要进行这步处理呢
Ix2 = Ix.^2;
Iy2 = Iy.^2;
Ixy = Ix.*Iy;
最后的 M = [Ix2(i,j) Ixy(i,j);Ixy(i,j) Iy2(i,j)]; 这个矩阵是怎么定义出来的呢?
先谢了

回复
分享到:

使用道具 举报

 楼主| 发表于 2008-6-18 12:35 | 显示全部楼层
版主大侠帮忙解释下啊,怎么没人看呢。。。
 楼主| 发表于 2008-6-18 13:09 | 显示全部楼层

关于滤波的梯度算子问题

fx = [-2 -1 0 1 2];                     % x方向梯度算子

该梯度算子用来实现X方向滤波,请问下这个梯度算子为什么定义为[-2 -1 0 1 2];

fx = [5 0 -5;8 0 -8;5 0 -5];          % x方向高斯函数一阶微分

第一个X方向梯度算子和第二个X方向高斯一阶微分又有什么区别呢?
发表于 2008-6-21 21:54 | 显示全部楼层
Harris角点检测法的原理先弄懂吧,程序没有问题,方向梯度算子可以改的,取决于你所需要的滤波效果

评分

1

查看全部评分

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-28 16:15 , Processed in 0.055397 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表