function [Tp1 Tp2 Tp3 Tp4 Tp5 Tp6] = Track_Point2(Start_x,Start_y,I,BW,row,column)
%%%%%code by zeng hui 边缘点跟踪
% I=wiener2(I,[3 3]);%%%%%去噪声
% I=intrans(I,'stretch',mean2(im2double(I)),8);
h = fspecial('gaussian', 3);
It = imfilter(I, h);
Gxm = [-1 -2 -1; 0 0 0; 1 2 1];
Gym = [-1 0 1; -2 0 2; -1 0 1];
Start_x0 = Start_x; Start_y0 = Start_y;
%%%%%%%%%%%%%%%%%%%%从初始点向右上跟踪
%%%%若连续出现15个右上,则截止。
Fb = 0;k = 0;Tp1 = [ ]; kb = 0;
while Fb ==0
k =k+1;
Np_c = [Start_x-1 Start_y;Start_x-1 Start_y+1;Start_x Start_y+1];
[ff,fg] = FTY(Np_c,BW);
% plot(Np_c(:,2),Np_c(:,1),'.r')
if ff ==1
Start_x = Np_c(fg,1);Start_y = Np_c(fg,2);
if fg ==1
kb = kb+1;
end
elseif ff==0
break;
else
for i=1:3
[Gxp_c(i),Gyp_c(i),Gp_c(i)] = Point_Filter_p(Np_c(i,1),Np_c(i,2),It,Gxm,Gym);
end
[ac as] =max(Gp_c);
if as ==1
kb = kb+1;
end
Start_x = Np_c(as,1);Start_y = Np_c(as,2);
end
if kb >5
Fb =1;
end
if Start_x<3
Tp1(k,:) = [Start_x,Start_y];
break;
end
Tp1(k,:) = [Start_x,Start_y];
end
%%%%%%%%%%%%%%%%%%%%从初始点向右下跟踪
for k=1:50
Np_c = [Start_x+1 Start_y;Start_x+1 Start_y+1;Start_x Start_y+1];
[ff,fg] = FTY(Np_c,BW);
% plot(Np_c(:,2),Np_c(:,1),'.r')
if ff ==1
Start_x = Np_c(fg,1);Start_y = Np_c(fg,2);
elseif ff==0
break;
else
for i=1:3
[Gxp_c(i),Gyp_c(i),Gp_c(i)] = Point_Filter_p(Np_c(i,1),Np_c(i,2),It,Gxm,Gym);
end
[ac as] =max(Gp_c);
Start_x = Np_c(as,1);Start_y = Np_c(as,2);
end
Tp2(k,:) = [Start_x,Start_y];
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Start_x = Start_x0; Start_y = Start_y0;
%%%%%%%%%%%%%%%%%%%%从初始点向左下跟踪
%%%%若出现一个像素点数目大于5的竖线,则截止。
Fb = 0;k = 0;Tp3 = [ ]; kb = 0;
while Fb ==0
k = k+1;
Np_c = [Start_x Start_y-1;Start_x+1 Start_y-1;Start_x+1 Start_y];
[ff,fg] = FTY(Np_c,BW);
% plot(Np_c(:,2),Np_c(:,1),'.r')
if ff ==1
Start_x = Np_c(fg,1);Start_y = Np_c(fg,2);
if fg ==3
kb = kb+1;
else
kb = 0;
end
elseif ff==0
break;
else
for i=1:3
[Gxp_c(i),Gyp_c(i),Gp_c(i)] = Point_Filter_p(Np_c(i,1),Np_c(i,2),It,Gxm,Gym);
end
[ac as] =max(Gp_c);
if as ==3
kb = kb+1;
else
kb =0;
end
Start_x = Np_c(as,1);Start_y = Np_c(as,2);
end
if kb >4
Fb =1;
end
if Start_y<3
Tp3(k,:) = [Start_x,Start_y];
break;
end
Tp3(k,:) = [Start_x,Start_y];
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%从向下跟踪
%%%%若连续出现40个右下,则截止。连续出现2个左下,截止
Fb = 0;k = 0;Tp4 = [ ]; kb = 0; kp =0;
while Fb ==0
k = k+1;
Np_c = [Start_x+1 Start_y-1;Start_x+1 Start_y;Start_x+1 Start_y+1];
[ff,fg] = FTY(Np_c,BW);
% plot(Np_c(:,2),Np_c(:,1),'.r')
if ff ==1
Start_x = Np_c(fg,1);Start_y = Np_c(fg,2);
if fg ==3
kb = kb+1;
end
if fg ==1
kp = kp+1;
else
kp =0;
end
else
for i=1:3
[Gxp_c(i),Gyp_c(i),Gp_c(i)] = Point_Filter_p(Np_c(i,1),Np_c(i,2),It,Gxm,Gym);
end
[ac as] =max(Gp_c);
if as ==3
kb = kb+1;
end
if as ==1
kp = kp+1;
else
kp =0;
end
Start_x = Np_c(as,1);Start_y = Np_c(as,2);
end
if kb >40
Fb =1;
end
if kp >2
Fb =1;
end
Tp4(k,:) = [Start_x,Start_y];
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%从初始点向右下跟踪
%%%若出现一个像素点数目大于5的横线,则截止
%%%%%在出现像素数目大于3的横线后,若出现4个以上的非边缘点,截止。
Fb = 0;k = 0;Tp5 = [ ]; kb = 0; kb2 = 0; kp =1;
while Fb ==0
k =k+1;
Np_c = [Start_x+1 Start_y;Start_x+1 Start_y+1;Start_x Start_y+1];
[ff,fg] = FTY(Np_c,BW);
% plot(Np_c(:,2),Np_c(:,1),'.r')
if ff==1
Start_x = Np_c(fg,1);Start_y = Np_c(fg,2);
if fg==3
kb=kb+1;
else
kb=0;
end
else
for i=1:3
[Gxp_c(i),Gyp_c(i),Gp_c(i)] = Point_Filter_p(Np_c(i,1),Np_c(i,2),It,Gxm,Gym);
end
[ac as] =max(Gp_c);
if as==3
kb=kb+1;
else
kb=0;
end
Start_x = Np_c(as,1);Start_y = Np_c(as,2);
end
if kb>7
Fb=1;
end
if row-Start_x<2
Tp5(k,:) = [Start_x,Start_y];
break;
end
Tp5(k,:) = [Start_x,Start_y];
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%从初始点向右上跟踪
%%%%若连续出现15个右上,则截止。
Fb = 0;k = 0;Tp6 = [ ]; kb = 0;
while Fb ==0
k =k+1;
Np_c = [Start_x-1 Start_y;Start_x-1 Start_y+1;Start_x Start_y+1];
[ff,fg] = FTY(Np_c,BW);
% plot(Np_c(:,2),Np_c(:,1),'.r')
if ff ==1
Start_x = Np_c(fg,1);Start_y = Np_c(fg,2);
elseif ff==0
break;
else
for i=1:3
[Gxp_c(i),Gyp_c(i),Gp_c(i)] = Point_Filter_p(Np_c(i,1),Np_c(i,2),It,Gxm,Gym);
end
[ac as] =max(Gp_c);
Start_x = Np_c(as,1);Start_y = Np_c(as,2);
end
if column-Start_y<2
Tp6(k,:) = [Start_x,Start_y];
break;
end
Tp6(k,:) = [Start_x,Start_y];
end
没有合适的资源?快使用搜索试试~ 我知道了~
本科毕业设计——人耳识别(SIFT&OIFT)
共408个文件
bmp:336个
m:46个
p:6个
4星 · 超过85%的资源 需积分: 21 87 下载量 82 浏览量
2011-05-02
23:08:26
上传
评论 6
收藏 8.29MB RAR 举报
温馨提示
人耳识别技术是一种新的生物特征识别技术,由于人耳独特的生理特征结构和生理位置,以及其不受外界环境(刺激)和内心活动对生物特征体影响的特点,逐渐引起同领域研究人员的广泛关注。本文主要从人耳图像的特征提取技术方面进行探索与研究。 人耳图像的特征提取是整个识别过程的核心环节。本文的工作主要是围绕如何提取人耳的局部特征和几何特征及如何将这些特征应用于人耳识别等问题展开的,主要工作有: 1)针对人耳图像存在大量相似纹理,直接应用SIFT描述子进行特征点匹配会产生大量误匹配的情况,提出了利用基于全局上下文信息的SIFT描述子进行图像匹配的方法。结合全局上下文信息的SIFT描述子不仅保持了SIFT描述子对图像尺度、旋转、光照变化和图像噪声的良好性能,而且可以较好地对相似纹理区域进行辨别分析。实验结果表明,基于结合全局上下文信息的SIFT描述子的图像匹配算法可有效避免相似纹理区域特征点之间的误匹配,提高了人耳图像匹配的效率。 2)针对已有的人耳几何特征提取方法受姿态变化影响都较为严重的情况,本文提出了一种基于射影不变量的人耳特征提取方法。首先,检测人耳长轴上的五个边缘特征点,然后,任取其中的4个边缘特征点并计算它们之间的交比,可得到5个交比值。最后,将计算出的5个交比值联合起来,构造成人耳特征向量。实验结果表明,与利用特征点长度比值作为人耳几何特征的识别方法相比,基于射影不变量的人耳识别方法对人耳姿态变化具有更好的鲁棒性。 3)为了进一步提高识别率,本文将基于结合全局上下文信息SIFT描述子的图像匹配方法与基于射影不变量的人耳几何特征提取方法结合起来用于人耳识别,并通过大量的实验验证了该方法的有效性。
资源推荐
资源详情
资源评论
收起资源包目录
本科毕业设计——人耳识别(SIFT&OIFT) (408个子文件)
Track_Point2.asv 6KB
Track_Point.asv 3KB
test.asv 2KB
Point_Filter_p.asv 2KB
Search_Start.asv 550B
Track_Point_Candidate_you.asv 416B
29-1.bmp 84KB
04-4.bmp 82KB
55-2.bmp 80KB
32-3.bmp 78KB
69-3.BMP 77KB
27-4.bmp 75KB
51-3.BMP 75KB
17-4.bmp 74KB
17-1.bmp 73KB
77-3.bmp 73KB
77-1.bmp 73KB
69-1.BMP 70KB
41-3.BMP 69KB
74-3.BMP 67KB
test.bmp 67KB
35-4.bmp 67KB
43-4.BMP 66KB
21-3.BMP 66KB
50-3.BMP 66KB
46-1.BMP 66KB
46-3.BMP 65KB
77-4.bmp 65KB
07-2.bmp 65KB
09-3.BMP 65KB
58-3.BMP 64KB
41-4.BMP 64KB
55-4.BMP 64KB
54-4.BMP 64KB
28-3.bmp 64KB
43-3.BMP 64KB
54-1.BMP 64KB
69-4.BMP 64KB
46-4.BMP 63KB
51-4.BMP 62KB
50-1.BMP 62KB
05-3.BMP 62KB
41-1.BMP 62KB
30-3.bmp 62KB
43-1.BMP 62KB
28-4.bmp 61KB
55-3.BMP 61KB
62-3.BMP 61KB
76-4.bmp 61KB
52-3.BMP 61KB
76-2.bmp 61KB
09-2.bmp 60KB
12-4.BMP 60KB
54-3.BMP 60KB
17-3.BMP 60KB
76-3.bmp 60KB
53-3.BMP 59KB
04-3.BMP 59KB
55-1.BMP 59KB
57-1.BMP 59KB
23-3.BMP 59KB
12-2.BMP 59KB
23-3.BMP 59KB
76-1.bmp 59KB
16-3.BMP 59KB
21-4.BMP 59KB
20-3.BMP 58KB
50-4.BMP 58KB
05-1.BMP 58KB
05-1.BMP 58KB
54-2.BMP 58KB
57-4.BMP 57KB
37-4.BMP 57KB
40-1.BMP 57KB
48-1.BMP 57KB
13-3.BMP 57KB
48-3.BMP 57KB
40-4.BMP 57KB
18-4.bmp 57KB
01-2.bmp 57KB
01-2.bmp 57KB
01-2.bmp 57KB
60-1.BMP 57KB
63-3.BMP 57KB
47-3.BMP 56KB
72-3.BMP 56KB
29-4.BMP 56KB
16-1.BMP 56KB
36-1.bmp 56KB
53-1.BMP 56KB
58-1.BMP 56KB
05-4.BMP 56KB
07-3.BMP 56KB
29-3.BMP 56KB
37-1.BMP 56KB
28-1.bmp 56KB
52-1.BMP 55KB
23-1.BMP 55KB
51-2.BMP 55KB
40-3.BMP 55KB
共 408 条
- 1
- 2
- 3
- 4
- 5
帅羊羊
- 粉丝: 2
- 资源: 6
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功
- 1
- 2
- 3
前往页