技术标签: # MATLAB 笔记 matlab 垂距 曲线压缩
垂距法是指根据中间顶点到其前、后两相邻顶点连线的距离的大小,来确定是否保留该顶点的一种 线要素顶点压缩 算法。当求得的距离大于给定的限差(阈值)时,保留该顶点,否则删除该顶点(如下图所示)。
一般使用所有点到直线距离的 中位数 作为阈值。同时,本人所编写的工具函数不仅可以处理二维数据点,三维甚至更高维度的数据点依旧可以处理。
怕大家找不到工具函数,这里放在最前面啦:
function [newPntSet,vertPnt]=getFeaturePnt(pntSet)
% @author:slandarer
% newPntSet : 特征点
% vertPnt : 原始曲线垂足
[rows,cols]=size(pntSet);
if rows<4||cols<2,error('数据点过少或维度异常');end
innerPntSet=pntSet(2:end-1,:); % 内部点
adjPnt1V=innerPntSet-pntSet(1:end-2,:); % 当前点与前一点向量
adjPnt2V=pntSet(3:end,:)-pntSet(1:end-2,:);% 当前点的两个相邻点向量
adjPnt2V=adjPnt2V./vecnorm(adjPnt2V')'; % 归一化
lVert=sum(adjPnt1V.*adjPnt2V,2);
vertPnt=adjPnt2V.*lVert+pntSet(1:end-2,:); % 获取垂足
vertL=vecnorm((innerPntSet-vertPnt)'); % 计算垂线距离
vertL(isinf(vertL)|isnan(vertL))=0; % 修正/0的情况
innerPntSet(vertL<median(vertL),:)=[]; % 删掉特征性不强的点
newPntSet=[pntSet(1,:);innerPntSet;pntSet(end,:)];
end
这里随机生成一组二维数据(两列),取出特征点后并绘图:
% 随机构造数据
X=linspace(0,25,10)';
Y=randi([0,10],[10,1]);
pntSet=[X,Y];
% 获取特征点
[nPntSet,vertPnt]=getFeaturePnt(pntSet);
% 坐标区域修饰
hold on
ax=gca;
ax.DataAspectRatio=[1,1,1];
% 绘制原始数据曲线
plot(pntSet(:,1),pntSet(:,2),'Color',[0 0.4470 0.7410],'LineWidth',2,'Marker','*');
% 绘制新数据曲线
plot(nPntSet(:,1),nPntSet(:,2),'Color',[0.6350 0.0780 0.1840 .7],'LineWidth',2,'Marker','s');
legend('original-curve','feature-curve')
就把辅助线也画上并各种加属性呗,需要注意的是,垂足并不一定在相邻两点的连线上,有时候是在其延长线上,但因为怕麻烦这里延长线就没有画:
% 随机构造数据
X=linspace(0,25,10)';
Y=randi([0,10],[10,1]);
pntSet=[X,Y];
% 获取特征点
[nPntSet,vertPnt]=getFeaturePnt(pntSet);
% 坐标区域修饰
hold on
ax=gca;
ax.YLim=[0,10];
ax.DataAspectRatio=[1,1,1];
ax.Color=[1,1,1];
ax.XColor=[1,1,1].*.3;
ax.YColor=[1,1,1].*.3;
ax.LineWidth=1.5;
ax.FontName='cambria';
% 绘制原始数据曲线
l1=plot(pntSet(:,1),pntSet(:,2),'Color',[0 0.4470 0.7410],'LineWidth',2,'Marker','*');
% 绘制辅助线及垂线
innerPntSet=pntSet(2:end-1,:);
l2=plot([innerPntSet(:,1),vertPnt(:,1)]',[innerPntSet(:,2),vertPnt(:,2)]','Color',[.3,.3,.3],'LineWidth',1.2,'LineStyle','-.');
plot([pntSet(3:end,1),pntSet(1:end-2,1)]',[pntSet(3:end,2),pntSet(1:end-2,2)]','Color',[.3,.3,.3],'LineWidth',1.2,'LineStyle','-.')
% 绘制新数据曲线
l3=plot(nPntSet(:,1),nPntSet(:,2),'Color',[0.6350 0.0780 0.1840 .7],'LineWidth',2,'Marker','s');
% 增添图例
legend([l1,l2(1),l3],{
'original-curve','auxiliary-line','feature-curve'});
也是几乎完全相同的使用方式,不过绘图用的plot3:
% 随机构造数据
X=linspace(0,25,10)';
Y=randi([0,10],[10,1]);
Z=randi([0,10],[10,1]);
pntSet=[X,Y,Z];
% 获取特征点
[nPntSet,vertPnt]=getFeaturePnt(pntSet);
% 坐标区域修饰
hold on
grid on
ax=gca;
ax.YLim=[0,10];
ax.ZLim=[0,10];
ax.DataAspectRatio=[1,1,1];
ax.Color=[1,1,1];
ax.XColor=[1,1,1].*.3;
ax.YColor=[1,1,1].*.3;
ax.ZColor=[1,1,1].*.3;
ax.LineWidth=1.5;
ax.FontName='cambria';
% 绘制原始数据曲线
l1=plot3(pntSet(:,1),pntSet(:,2),pntSet(:,3),'Color',[0 0.4470 0.7410],'LineWidth',2,'Marker','*');
% 绘制辅助线
innerPntSet=pntSet(2:end-1,:);
l2=plot3([innerPntSet(:,1),vertPnt(:,1)]',[innerPntSet(:,2),vertPnt(:,2)]',...
[innerPntSet(:,3),vertPnt(:,3)]','Color',[.3,.3,.3],'LineWidth',1.2,'LineStyle','-.');
plot3([pntSet(3:end,1),pntSet(1:end-2,1)]',[pntSet(3:end,2),pntSet(1:end-2,2)]',...
[pntSet(3:end,3),pntSet(1:end-2,3)]','Color',[.3,.3,.3],'LineWidth',1.2,'LineStyle','-.')
% 绘制新数据曲线
l3=plot3(nPntSet(:,1),nPntSet(:,2),nPntSet(:,3),'Color',[0.6350 0.0780 0.1840 .7],'LineWidth',2,'Marker','s');
% 增添图例
legend([l1,l2(1),l3],{
'original-curve','auxiliary-line','feature-curve'});
view(3)
文章浏览阅读2.5k次,点赞5次,收藏17次。2018年3月,我与张老师就这么在微信上聊了起来,起初我并没有写书的打算,我们之间只是通过讨论、交流的形式聊聊关于出书的方方面面。最终,敌不过张老师超强的专业能力、细致的解说与盛情相邀,我答应张老师写一本Linux系统运维的图书并由人邮出版。由此,我踏上了漫漫2年多的写书之路。为什么写这本书写书一方面是我对自己所学知识的查漏补缺过程,另一方面也可以向即将进入或已经入行的Linux系统运维同..._linux系统运维指南:从入门到企业实战 pdf
文章浏览阅读2k次,点赞6次,收藏5次。tf.reduce_sum()函数深度解析从矩阵,数组,数据存储的角度 解析axis参数的意义_tf.reduce_sum
文章浏览阅读9.8k次,点赞4次,收藏29次。adb获取app包名的方法_adb获取包名
文章浏览阅读913次,点赞16次,收藏10次。总之,要做好虾皮店铺,不仅需要明确的定位和优质的产品,还需要精心的运营和持续的改进。通过不断优化店铺形象、制定有效的营销策略、提供优质的客户服务以及加强供应链管理等手段,您将能够在激烈的竞争中脱颖而出,实现店铺的长足发展。1.稳定的网络环境是基石,它需要经过技术手段的洗礼,将电脑或手机的底层硬件参数伪装成国外数据,以躲避平台通过IP进行的深度检测。这种真实性高的评价能够帮助商家获得更多的信任和认可,从而提升产品的排名和流量的分配。您可以关注行业动态,学习先进的经营理念和技术,以提高店铺的运营水平。
文章浏览阅读5k次,点赞11次,收藏43次。统计检验_统计测试 cd diagrams
文章浏览阅读332次。购物车_购物车案例请求数据地址
文章浏览阅读2.9k次,点赞3次,收藏37次。1. IIS1. PUT漏洞用户配置不当,exp:https://github.com/hackping/HTTPMLScan.git2. 短文件名猜解IIS的短文件名机制,可以暴力猜解短文件名,访问构造的某个存在的短文件名,会返回404,访问构造的某个不存在的短文件名,返回400。exp:https://github.com/WebBreacher/tilde_enum3.远程代码执行(CVE-2017-7269))**exp**:https://github.com/zcgonv_hrs中间件
文章浏览阅读368次。DB2支持以下两种类型的表空间: 1、 系统管理存储器表空间(SMS-SYSTEM MANAGED STORAGE) 2、 数据库管理存储器表空间(DMS-DATABASE MANAGED STORAGE) SMS、DMS用户表空间的特性对照 特性 ..._db2
文章浏览阅读84次。正在开发中的游戏有个全屏功能--可以在window桌面背景上运行,就像一些视频播放器在桌面背景上播放一样的,花了个上午整了个Demo放出来留个纪念。实现功能:显示图标,双击图标执行相应的程序,右击图标弹出该图标对应得菜单,点击非图标区则弹出桌面菜单。需要完整工程可以点此下载:DesktopWindow.rar。程序效果图如下:在这个程序里,定义了一个XShellItem..._模拟实现windows桌面效果
文章浏览阅读944次。https://www.byhy.net/tut/webdev/django/01/_byhy.net
文章浏览阅读5.8k次,点赞13次,收藏57次。业务场景介绍:H5移动端支持微信支付 [ 微信支付分为微信内支付(JSAPI支付官方API)和微信外支付(H5支付官方API)] && 支付宝支付 [手机网站支付转 APP 支付 官方API ]订单生成逻辑:前端请求后端提交订单,后端去和微信或者支付宝对接生成订单(后续支付都是这个逻辑进行的对接)一、移动端微信支付,vue中如何玩?在移动端微信支付分为微信内支付和微信外支付。1.在订单组件中选择支付方式之后在支付页面先去判断是否是在微信内://判断是否微信 is__移动端支付宝微信支付vue项目怎么写
文章浏览阅读2k次,点赞5次,收藏9次。深度学习编译器主要为解决不同框架下训练的模型部署到指定的某些设备上时所遇到的一系列复杂的问题,即将各种深度学习训练框架的模型部署到各种硬件所面临的问题;_tvm编译器