clear all;

clc;

aa=xlsread('E:\chilun05188.xlsx');

N=length(aa);

ang=0:(0.4*2*205*pi*0.000001/60):(0.4*2*205*pi*0.000001/60)*(N-1);

ang=ang';

aa=aa+116.6479;

%基于三点说构成的夹角提取特征点

cc=1;

s=20;  %间隔s个点进行一次计算

t=0.12;     %用于调整特征点的位置

for k=1:s:N-2*s

A=[ang(k) aa(k)];

B=[ang(k+s) aa(k+s)];

C=[ang(k+2*s) aa(k+2*s)];

a1=B(1)-A(1)+j*(B(2)-A(2));

a2=C(1)-B(1)+j*(C(2)-B(2));

alpha(1,cc)=(angle(a2)-angle(a1))*180/pi;

cc=cc+1;

end

alpha=alpha';

tw=find(abs(alpha)>t);   %满足要求的齿顶齿根数据点在原数组中的索引位置

l_tw=length(tw);

tw_diff=diff(tw);

w_tw=find(tw_diff>20);

w1_feature_points1=tw(w_tw)*s;

w2_feature_points1=tw(w_tw+1)*s;

tw_one=tw(1)*s

tw_last=tw(l_tw)*s

w1_feature_points=[w1_feature_points1' tw_last];    %齿形分离特征点在原数组中的索引位置

w2_feature_points=[tw_one w2_feature_points1'];     %齿形分离特征点在原数组中的索引位置

l_w1_feature_points=length(w1_feature_points');

w_feature_points=sort([w1_feature_points w2_feature_points]);

t1=1;

t2=1;

for k1=1:2:l_w1_feature_points

w_one_1(1,t1)=w_feature_points(2*k1-1);

w_one_2(1,t2)=w_feature_points(2*k1);

t1=t1+1;

t2=t2+1;

end

w_one=sort([w_one_1 w_one_2]);       %齿顶特征点在原数组中的索引位置

t3=1;

t4=1;

for k2=2:2:l_w1_feature_points

w_other_1(1,t3)=w_feature_points(2*k2-1);

w_other_2(1,t4)=w_feature_points(2*k2);

t3=t3+1;

t4=t4+1;

end

w_other=sort([w_other_1 w_other_2]);     %齿根特征点在原数组中的索引位置

figure(1)

polar(ang,aa,'k')

hold on

polar(ang(w_one),aa(w_one),'.r')

hold on

polar(ang(w_other),aa(w_other),'.g')

figure(2)

plot(ang,aa)

hold on

plot(ang(w_one),aa(w_one),'.r','MarkerSize',20)

hold on

plot(ang(w_other),aa(w_other),'.g','MarkerSize',20)

基于三点说构成的夹角提取特征点相关推荐

  1. java kinect 人体识别_基于三个kinect的人体建模

    基于三个kinect的人体重建,首先要对三个kinect相机进行标定,得到相机内外参数.这样每个相机看到空间一点的坐标是一致的,三个相机坐标系得以统一.标定方法是两两标定,即取某一kinect相机(记 ...

  2. python实现dem输出三维模型_资源三号卫星影像立体像对如何提取DEM数据的方法

    原标题:资源三号卫星影像立体像对如何提取DEM数据的方法 OrthoMapping是ArcGIS 10.5推出的基于无人机.大飞机.卫星拍摄的原始影像获取专业级别信息产品的生产能力.使用OrthoMa ...

  3. 基于三个kinect的人体建模

    单个kinect的人体重建,在Kinect SDK 1.8中,Kinect Fusion的效果已经很不错了.其缺点显而易见,一是扫描时间长,重建对象也需要长时间保持静态:二是需要人体或者kinect转 ...

  4. matlab点云聚类,基于区域聚类分割的点云特征线提取

    王晓辉 , 吴禄慎 , 陈华伟 , 胡赟 , 石雅莹 . . 基于区域聚类分割的点云特征线提取. 光学学报, 2018, 38(11): 1110001-. Wang Xiaohui , Wu Lus ...

  5. 基于K-Means聚类算法的主颜色提取

    点击上方"小白学视觉",选择加"星标"或"置顶" 重磅干货,第一时间送达 01.简介 本期我们将一起实现基于K-Means聚类算法的主色提取 ...

  6. 零基础基于U-Net网络实战眼底图像血管提取

    文章目录 1 前言 2 血管提取任务概述 3 U-Net架构简介 4 眼底图像血管分割代码 5 结果评估可视化(ROC曲线) 6 改进U-Net网络完成眼底图像血管提取任务思路 1 前言 本文基于U- ...

  7. 图像五值化与基于三值图像的车牌识别(1)

    庆祝2011年新年 张忠义 海口愚佬会教育科技有限公司 全文下载 www.yulaohui.com/color5_2/ 摘 要:      图像五值化是通过公式Y =(R+G)/2,将RGB空间映射到 ...

  8. ML之SVM(三种):基于三种SVM(linearSVR、polySVR、RBFSVR)对Boston(波士顿房价)数据集(506,13+1)进行价格回归预测并对比各自性能

    ML之SVM(三种):基于三种SVM(linearSVR.polySVR.RBFSVR)对Boston(波士顿房价)数据集(506,13+1)进行价格回归预测并对比各自性能 目录 输出结果 设计思路 ...

  9. EL之DTRFGBT:基于三种算法(DT、RF、GBT)对泰坦尼克号乘客数据集进行二分类(是否获救)预测并对比各自性能

    EL之DT&RF&GBT:基于三种算法(DT.RF.GBT)对泰坦尼克号乘客数据集进行二分类(是否获救)预测并对比各自性能 目录 输出结果 ​设计思路 核心代码 输出结果 设计思路 核 ...

  10. 基于三值图像的车牌识别测试程序发布

    基于三值图像的车牌识别测试程序正式发布,初始版本是0.9.0.由于还有匹配模板没有完全导入,因此,识别精度还低.但是,目的是演示识别过程.其意义是,1.基于彩色图像直接进行三值化(即彩色二值化)并基于 ...

最新文章

  1. python文字教程-Python在图片中添加文字的两种方法
  2. python sin(x)/x 图像
  3. 模块化数据机房具备的五大优势
  4. vba 单元格 一系例操作
  5. # 20175213 2018-2019-2 《Java程序设计》第1周学习总结
  6. js深度复制一个对象使用JSON.stringify是最好的办法吗
  7. MySQL5.5.15_linux下mysql-5.5.15安装详细步骤
  8. WCF与AJAX编程开发实践(1):AJAX基础概念和纯AJAX示例
  9. C++编写简单的俄罗斯方块游戏
  10. 事态升级是什么意思_俗语“头不顶桑脚不踩槐;生不睡柳死不睡杨”啥意思 有道理吗...
  11. 【大数据-第二期】java基础第五天作业
  12. vue中引用swiper轮播插件
  13. xp 游戏计算机没有了怎么办,我的xp系统没有自带游戏怎么办
  14. 【POJ P3311】【状压DP】Hie with the Pie
  15. 异常处理(六)--------SpringBoot+Maven项目运行异常:Unable to find a single main class from the following candidat
  16. 【解决方案 十二】一文彻底解决文件格式判别问题
  17. 13.JVM-垃圾回收相关概念
  18. Python之基础详解(十二):进阶,用pyecharts绘制世界地图、玫瑰图、雷达图
  19. 编程之路上的感悟-编程语言该怎么选
  20. 空调噪音测试软件,检测空调噪音的四种方法

热门文章

  1. 大学计算机课程教学建议,计算机应用基础课程教学的建议
  2. Uniapp或H5之ORC识别与自定义照相机
  3. Python代码编辑器jupyter的使用
  4. 网站建设中百度快照劫持是什么?劫持百度快照是怎么回事?
  5. 并行计算,网格计算与分布式计算的…
  6. ubuntu 下启用ipv6
  7. java实时监控系统_JavaMonitor
  8. 流控大师 panabit
  9. HDU-5855 Less Time, More profit(最大权闭合图+二分)
  10. java短信接口开发完整项目_java项目接入第三方短信接口