码迷,mamicode.com
首页 > 其他好文 > 详细

Matlab程序学习(散点边界线/包络线)

时间:2017-06-26 18:55:43      阅读:3403      评论:0      收藏:0      [点我收藏+]

标签:tle   数据   mat   led   color   不规则   sort   length   主程   

针对散点边界线,Matlab重的convhull为凸包络线,往往不能满足实际工作需要。

作散点的紧致包络线,需要考虑到凹凸点,为减少包络线与散点数据的过度紧密(穿过全散点),并不过度包络(类似convhull),从以下思路来实现:

1.散点数据按行走方向排序,并判断点的凹凸;

2.针对凹点,添加阈值,通过夹角判断,去除部分凹点,然后继续此种循环判断,直到所有凹点满足要求;

程序文件:


function y2=ScatterHull(Data,ang_2)
% Thinking from Dsp Tian 
% http://www.cnblogs.com/tiandsp/p/4006575.html
% y1(:,1)=x;y1(:,2)=y;ang_2:concave min angle
% ang_2=180:conhull;ang_2=0:all scatters
y1=unique(Data,rows);% Remove duplicate rows 
%%
% sort of scatters
    cen=mean(y1);
    ang=atan2(y1(:,1)-cen(1),y1(:,2)-cen(2)); %每个点到坐标中心极角
    y1=[y1,ang];
    y1=sortrows(y1,3);    %按极角排序
    y1=y1(:,1:2);
%%
y2=y1;
im=1;
[n m]=size(y1);
%%
while length(im)>0
 
im=[];
%向量计算,判断凹点
for i=1:n    
    if i==1  %处理第一个点
        v1=y1(n,:)-y1(1,:);
        v2=y1(2,:)-y1(1,:); 
    elseif i==n %最后一个点
        v1=y1(n-1,:)-y1(n,:);    
        v2=y1(1,:)-y1(n,:);        
    else 
        v1=y1(i-1,:)-y1(i,:);     
        v2=y1(i+1,:)-y1(i,:);
    end
    r=det([v1;v2]);  %叉乘后向量方向
    c=dot(v1,v2)/(norm(v1)*norm(v2));
    ang=rad2deg(acos(c));
        
    if r<=0&(ang<ang_2|isnan(ang)==1)%concave
            ang(ang<ang_2)=[];
            im=[im;i];
    end
end
y2(im,:)=[];
y1=y2;
[n m]=size(y1);
end
%%
y2=[y2;y2(1,:)];%results curve is closed
end


 

 

主程序:

function samplef
%UNTITLED5 Summary of this function goes here
%   Detailed explanation goes here
close all
clear all
clc
t=load(data2.txt);
t(t(:,1)==0|t(:,2)==0,:)=[];
x=t(:,1);
y=t(:,2);
figure(1)
ang=[0 60 120 180];
for i=1:4
t2=ScatterHull(t,ang(i));
t2=[t2;t2(1,:)];%画图使封闭
subplot(2,2,i)
hold on
plot(t(:,1),t(:,2),ro)
plot(t2(:,1),t2(:,2),g-)
legend(scatter,HullCurve,Location,NorthWest)
title(num2str(ang(i)))
end
end

 

结果:

 技术分享

注意:部分散点极度不规则的情况,效果并不理想,例如散点图明显分区,具有几个明显的聚点中心

 

 

 

Matlab程序学习(散点边界线/包络线)

标签:tle   数据   mat   led   color   不规则   sort   length   主程   

原文地址:http://www.cnblogs.com/ClownZzb/p/7081830.html

(0)
(1)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!