使用matlab进行DBscan聚类_dbscan聚类分析图用什么软件-程序员宅基地

技术标签: matlab  聚类  计算机视觉  

使用matlab对输入数据进行DBscan聚类。算法的思想来自基于密度方法的聚类
可直接运行代码以及数据文件可从此下载

在这里插入图片描述
在这里插入图片描述

clear;
close all;
clc;
k = 3;
Eps = 2;

%% 生成模拟数据
n = 200;
a = linspace(0,8*pi,n/2);
u = [5*cos(a)+5 10*cos(a)+5]'+1*rand(n,1);
v = [5*sin(a)+5 10*sin(a)+5]'+1*rand(n,1);
mu1 = [20 20];
S1 = [10 0;0 10];
data1 = mvnrnd(mu1,S1,100);
data = [u v;data1];

% image = imread('data.png');
% image = image(:,:,1);
% [x,y]=find(image == 0);
% data=[x,y];

%% 准备变量,输出原始结果
[m,n] = size(data);
data=[(1:m)',data];
n = n + 1;
type = zeros(1,m);
cluster_No = 1;
visited = zeros(m,1);
class = zeros(1,m)-2;

figure(2);
plot(data(:,2),data(:,3),'k.');
grid on
daspect([1 1 1]);
xlabel('x');ylabel('y');
title('原始输入点');
hold on;

%% DBscan
Kdtree = KDTreeSearcher(data(:,2:3));

for i = 1:m
    % 抽取一个未访问点
    if visited(i)==0
        % 标为访问
        visited(i) = 1;
        point_now = data(i,:);
        Idx_range = rangesearch(Kdtree, point_now(2:3), Eps);
        index = Idx_range{
    1};
        if length(index) > k
            class(i) = cluster_No;
            while index
                if visited(index(1)) == 0
                    visited(index(1)) = 1;
                    if class(index(1)) <= 0
                        class(index(1)) = cluster_No;
                    end
                    point_now = data(index(1),:);
                    Idx_range = rangesearch(Kdtree, point_now(2:3), Eps);
                    index_temp = Idx_range{
    1};
                    index(1) = [];
                    if length(index_temp) > k
                        index = [index, index_temp];
                    end
                else
                    index(1) = [];
                end
            end
            cluster_No = cluster_No + 1;
        end
    end
end

%% DBscan聚类结果
figure;
for i = 1: cluster_No
    color = [rand(),rand(),rand()];
    data_class = data(find(class==i),:);
    plot(data_class(:,2),data_class(:,3),'.','Color',color,'MarkerFaceColor',color);
    hold on
end
data_class = data(find(class<=0),:);
plot(data_class(:,2),data_class(:,3),'k*');
hold on
grid on
daspect([1 1 1]);
xlabel('x');ylabel('y');
title('DBscan聚类结果');
版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/john_xia/article/details/107662116

智能推荐

上传计算机桌面文件图标不见,关于桌面上图标都不见了这类问题的解决方法-程序员宅基地

文章浏览阅读8.9k次。关于桌面上图标都不见了这类问题的解决方法1、在桌面空白处右击鼠标-->排列图标-->勾选显示桌面图标。2、如果问题还没解决,那么打开任务管理器(同时按“Ctrl+Alt+Del”即可打开),点击“文件”→“新建任务”,在打开的“创建新任务”对话框中输入“explorer”,单击“确定”按钮后,稍等一下就可以见到桌面图标了。3、问题还没解决,按Windows键+R(或者点开始-->..._上传文件时候怎么找不到桌面图标

LINUX 虚拟网卡tun例子——修改_怎么设置tun的接收缓冲-程序员宅基地

文章浏览阅读1.5k次。参考:http://blog.csdn.net/zahuopuboss/article/details/9259283 #include #include #include #include #include #include #include #include #include #include #include #include _怎么设置tun的接收缓冲

UITextView 评论输入框 高度自适应-程序员宅基地

文章浏览阅读741次。创建一个inputView继承于UIView- (instancetype)initWithFrame:(CGRect)frame{ self = [superinitWithFrame:frame]; if (self) { self.backgroundColor = [UIColorcolorWithRed:0.13gre

赛效:怎么一键人像抠图-程序员宅基地

文章浏览阅读149次。2:人像图上传成功后自动进入抠图页面,我们点击下方的“下载图片”,可将抠图后的人像图下载到电脑本地。3:图片下载成功后,效果如下图所示。如果对抠图要求较高,可以自行修改细节。1:在人像抠图页面,点击下方的“本地上传”,将本地人像图片添加上去。”官网输入关键词,搜索查找相关教程。如果你想了解更多好用的。

4.html+css+javascript网页设计实例/“成都“旅游主题介绍/web前端期末大作业/-程序员宅基地

文章浏览阅读1k次,点赞23次,收藏21次。本实例为html+css+javascript代码。该实例里面有设置导航栏效果、动态轮播特效、DIV、Banner、表格、视频、注册、登录页面等实例比较全面,有助于各位大小朋友学习,本文以“成都”旅游为主题来设计实践,介绍如何通过从头开始设计旅游网站并将其转换为代码的过程。

PerfDog-移动端性能测试-基本使用_perfdog在跑性能的过程中能录屏吗-程序员宅基地

文章浏览阅读941次。常见的腾讯性能测试工具:腾讯gt、腾讯wetest、腾讯perfdog腾讯perfdog: https://perfdog.qq.com/一、介绍:移动全平台iOS/Android性能测试、分析工具平台。快速定位分析性能问题,提升APP应用及游戏性能和品质。手机无需ROOT/越狱,手机硬件、游戏及应用APP也无需做任何修改,极简化即插即用。Windows & Mac OS X平台PerfDog桌面应用程序版本都支持对iOS和Android设备进行测试。PC上PerfDog可多开,单PC可._perfdog在跑性能的过程中能录屏吗

随便推点

element-ui switch开关打开和关闭时的文字设置样式-程序员宅基地

文章浏览阅读3.3k次,点赞2次,收藏2次。element switch开关文字显示element中switch开关把on-text 和 off-text 属性改为 active-text 和 inactive-text 属性.怎么把文字描述显示在开关上?下面就是实现方法: 1 <el-table-column label="状态"> 2 <template slot-scope="scope">..._el-switch 不同状态显示不同字

HttpRequestUtil方法get、post、JsonToPost_httprequestutil.httpget-程序员宅基地

文章浏览阅读785次。java后台发起请求使用的工具类package com.cennavi.utils;import org.apache.http.Header;import org.apache.http.HttpResponse;import org.apache.http.HttpStatus;import org.apache.http.client.HttpClient;import org.apache.http.client.methods.HttpPost;import org.apach_httprequestutil.httpget

App-V轻量级应用程序虚拟化之三客户端测试-程序员宅基地

文章浏览阅读137次。在前两节我们部署了App-V Server并且序列化了相应的软件,现在可谓是万事俱备,只欠东风。在这篇博客里面主要介绍一下如何部署客户端并实现应用程序的虚拟化。在这里先简要的说一下应用虚拟化的工作原理吧!App-V Streaming 就是利用templateServer序列化出一个软件运行的虚拟环境,然后上传到app-v Server上,最后客户..._app-v 客户端

实时视频传输方案汇总-java_eclipse视频传输设计-程序员宅基地

文章浏览阅读4.6k次。实时视频传输方案汇总-java目录libstreamingh264jrtsp-h264-clientOpenFlowihmc-video-codecs目录libstreaminghttps://github.com/fyhertz/libstreamingIntroductionWhat it doeslibstreaming is an API that allows you, wi..._eclipse视频传输设计

unbantu上python安装步骤_如何在Ubuntu中安装Python 3.6?-程序员宅基地

文章浏览阅读230次。Python是增长最快的主要通用编程语言。原因有很多,比如它的可读性和灵活性,易于学习和使用,可靠和高效。有两个主要的Python版本被使用- 2和3 (Python的现在和未来);前者将看不到新的主要版本,而后者正在积极开发中,在过去几年已经看到了许多稳定的版本。Python 3的最新稳定版本是3.6版。Ubuntu 18.04和Ubuntu 17.10都预装了Python 3.6,这与老版本的..._ubuntu 14.04 安装python 3.6

NXP NFC Reader Library 移植思路_nxpnfcreader-程序员宅基地

文章浏览阅读1.1k次。记录 NXP 的 SDK 移植到 STM32 的过程,没有完成移植,中途暂停了。 _nxpnfcreader

推荐文章

热门文章

相关标签