MATLAB实例：散点密度图

MATLAB绘制用颜色表示数据密度的散点图

1. demo.m

% 用颜色表示数据密度的散点图

data_load=dlmread('E:\scanplot\gauss.txt');

X=data_load(:,1:2);

scatplot(X(:,1),X(:,2),'circles', sqrt((range(X(:, 1))/30)^2 + (range(X(:,2))/30)^2), 100, 5, 1, 8);

% colormap jet

print(gcf,'-dpng','散点密度图.png');

2. scatplot.m

function out = scatplot(x,y,method,radius,N,n,po,ms)

% Scatter plot with color indicating data density

% https://www.mathworks.com/matlabcentral/fileexchange/8577-scatplot

% USAGE:

% out = scatplot(x,y,method,radius,N,n,po,ms)

% out = scatplot(x,y,dd)

%

% DESCRIPTION:

% Draws a scatter plot with a colorscale

% representing the data density computed

% using three methods

%

% INPUT VARIABLES:

% x,y - are the data points

% method - is the method used to calculate data densities:

% 'circles' - uses circles with a determined area

% centered at each data point

% 'squares' - uses squares with a determined area

% centered at each data point

% 'voronoi' - uses voronoi cells to determin data densities

% default method is 'voronoi'

% radius - is the radius used for the circles or squares

% used to calculate the data densities if

% (Note: only used in methods 'circles' and 'squares'

% default radius is sqrt((range(x)/30)^2 + (range(y)/30)^2)

% N - is the size of the square mesh (N x N) used to

% filter and calculate contours

% default is 100

% n - is the number of coeficients used in the 2-D

% running mean filter

% default is 5

% (Note: if n is length(2), n(2) is tjhe number of

% of times the filter is applied)

% po - plot options:

% 0 - No plot

% 1 - plots only colored data points (filtered)

% 2 - plots colored data points and contours (filtered)

% 3 - plots only colored data points (unfiltered)

% 4 - plots colored data points and contours (unfiltered)

% default is 1

% ms - uses this marker size for filled circles

% default is 4

%

% OUTPUT VARIABLE:

% out - structure array that contains the following fields:

% dd - unfiltered data densities at (x,y)

% ddf - filtered data densities at (x,y)

% radius - area used in 'circles' and 'squares'

% methods to calculate densities

% xi - x coordenates for zi matrix

% yi - y coordenates for zi matrix

% zi - unfiltered data densities at (xi,yi)

% zif - filtered data densities at (xi,yi)

% [c,h] = contour matrix C as described in

% CONTOURC and a handle H to a contourgroup object

% hs = scatter points handles

%

%Copy-Left, Alejandro Sanchez-Barba, 2005

if nargin==0

scatplotdemo

return

end

if nargin<3 | isempty(method)

method = 'vo';

end

if isnumeric(method)

gsp(x,y,method,2)

return

else

method = method(1:2);

end

if nargin<4 | isempty(n)

n = 5; %number of filter coefficients

end

if nargin<5 | isempty(radius)

radius = sqrt((range(x)/30)^2 + (range(y)/30)^2);

end

if nargin<6 | isempty(po)

po = 1; %plot option

end

if nargin<7 | isempty(ms)

ms = 7; %markersize

end

if nargin<8 | isempty(N)

N = 100; %length of grid

end

%Correct data if necessary

x = x(:);

y = y(:);

%Asuming x and y match

idat = isfinite(x);

x = x(idat);

y = y(idat);

holdstate = ishold;

if holdstate==0

cla

end

hold on

%--------- Caclulate data density ---------

dd = datadensity(x,y,method,radius);

%------------- Gridding -------------------

xi = repmat(linspace(min(x),max(x),N),N,1);

yi = repmat(linspace(min(y),max(y),N)',1,N);

zi = griddata(x,y,dd,xi,yi);

%----- Bidimensional running mean filter -----

zi(isnan(zi)) = 0;

coef = ones(n(1),1)/n(1);

zif = conv2(coef,coef,zi,'same');

if length(n)>1

for k=1:n(2)

zif = conv2(coef,coef,zif,'same');

end

end

%-------- New Filtered data densities --------

ddf = griddata(xi,yi,zif,x,y);

%----------- Plotting --------------------

switch po

case {1,2}

if po==2

[c,h] = contour(xi,yi,zif);

out.c = c;

out.h = h;

end %if

hs = gsp(x,y,ddf,ms);

out.hs = hs;

colorbar

case {3,4}

if po>3

[c,h] = contour(xi,yi,zi);

out.c = c;

end %if

hs = gsp(x,y,dd,ms);

out.hs = hs;

colorbar

end %switch

%------Relocate variables and place NaN's ----------

dd(idat) = dd;

dd(~idat) = NaN;

ddf(idat) = ddf;

ddf(~idat) = NaN;

%--------- Collect variables ----------------

out.dd = dd;

out.ddf = ddf;

out.radius = radius;

out.xi = xi;

out.yi = yi;

out.zi = zi;

out.zif = zif;

if ~holdstate

hold off

end

return

%~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

function scatplotdemo

po = 2;

method = 'squares';

radius = [];

N = [];

n = [];

ms = 5;

x = randn(1000,1);

y = randn(1000,1);

out = scatplot(x,y,method,radius,N,n,po,ms)

return

%~~~~~~~~~~ Data Density ~~~~~~~~~~~~~~

function dd = datadensity(x,y,method,r)

%Computes the data density (points/area) of scattered points

%Striped Down version

%

% USAGE:

% dd = datadensity(x,y,method,radius)

%

% INPUT:

% (x,y) - coordinates of points

% method - either 'squares','circles', or 'voronoi'

% default = 'voronoi'

% radius - Equal to the circle radius or half the square width

Ld = length(x);

dd = zeros(Ld,1);

switch method %Calculate Data Density

case 'sq' %---- Using squares ----

for k=1:Ld

dd(k) = sum( x>(x(k)-r) & x(y(k)-r) & y

end %for

area = (2*r)^2;

dd = dd/area;

case 'ci'

for k=1:Ld

dd(k) = sum( sqrt((x-x(k)).^2 + (y-y(k)).^2) < r );

end

area = pi*r^2;

dd = dd/area;

case 'vo' %----- Using voronoi cells ------

[v,c] = voronoin([x,y]);

for k=1:length(c)

%If at least one of the indices is 1,

%then it is an open region, its area

%is infinity and the data density is 0

if all(c{k}>1)

a = polyarea(v(c{k},1),v(c{k},2));

dd(k) = 1/a;

end %if

end %for

end %switch

return

%~~~~~~~~~~ Graf Scatter Plot ~~~~~~~~~~~

function varargout = gsp(x,y,c,ms)

%Graphs scattered poits

map = colormap;

ind = fix((c-min(c))/(max(c)-min(c))*(size(map,1)-1))+1;

h = [];

%much more efficient than matlab's scatter plot

for k=1:size(map,1)

if any(ind==k)

h(end+1) = line('Xdata',x(ind==k),'Ydata',y(ind==k), ...

'LineStyle','none','Color',map(k,:), ...

'Marker','.','MarkerSize',ms);

end

end

if nargout==1

varargout{1} = h;

end

return

3. 结果

MATLAB实例：绘制折线图

MATLAB实例:绘制折线图 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 条形图的绘制见:MATLAB实例:绘制条形图 用MATLAB将几组不同的数 ...

MATLAB实例：构造网络连接图&lpar;Network Connection&rpar;及计算图的代数连通度&lpar;Algebraic Connectivity&rpar;

MATLAB实例:构造网络连接图(Network Connection)及计算图的代数连通度(Algebraic Connectivity) 作者:凯鲁嘎吉 - 博客园 http://www.cnbl ...

Matlab plotyy画双纵坐标图实例

Matlab plotyy画双纵坐标图实例 x = 0:0.01:20;y1 = 200*exp(-0.05*x).*sin(x);y2 = 0.8*exp(-0.5*x).*sin(10*x);[A ...

MATLAB实例：求相关系数、绘制热图并找到强相关对

MATLAB实例:求相关系数.绘制热图并找到强相关对 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 用MATLAB编程,求给定数据不同维度之间的相关系 ...

MATLAB实例：聚类网络连接图

MATLAB实例:聚类网络连接图 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 本文给出一个简单实例,先生成2维高斯数据,得到数据之后,用模糊C均值( ...

MATLAB实例：二元高斯分布图

MATLAB实例:二元高斯分布图 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 1. MATLAB程序 %% demo Multivariate No ...

Python图表数据可视化Seaborn：1&period; 风格&vert; 分布数据可视化-直方图&vert; 密度图&vert; 散点图

conda  install seaborn  是安装到jupyter那个环境的 1. 整体风格设置 对图表整体颜色.比例等进行风格设置,包括颜色色板等调用系统风格进行数据可视化 set() / se ...

MATLAB实例：绘制条形图

MATLAB实例:绘制条形图 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 用MATLAB绘制条形图,自定义条形图的颜色.图例位置.横坐标名称.显示条 ...

MATLAB实例：将批量的图片保存为&period;mat文件

MATLAB实例:将批量的图片保存为.mat文件 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 一.彩色图片 图片数据:horse.rar 1. MA ...

随机推荐

体验报告：微信小程序在安卓机和苹果机上的区别

很多人可能会问:微信小程序和在微信里面浏览一个网页有什么区别? 首先,小程序的运行是全屏的,界面跟进入了一个APP很像,更为沉浸跟在微信里面访问h5不一样:其次,它的浏览体验更为稳定. 不过,这还不够 ...

POJ1014 解题报告&lpar;DFS&rpar;

题目在此:http://poj.org/problem?id=1014 要看清题意呢,题中要求输入的是价值分别为1,2,3,4,5,6的大理石的个数,而不是6块价值为输入数字的大理石!选这个题主要想练 ...

51nod 1412 AVL树的种类&lpar;dp&rpar;

题目链接:51nod 1412 AVL树的种类 开始做的时候把深度开得过小了结果一直WA,是我天真了.. #include #include ...

Regex&period;Replace的基本用法

Regex构造函数Regex(string pattern)Regex(string pattern,RegexOptions options)参数说明pattern:要匹配的正则表达式模式optio ...

使用Redis的理由

Redis是一个远程内存数据库,它不仅性能强劲,而且还具有复制特性以及为解决问题而生的独一无二的数据模型.Redis提供了5种不同类型的数据结构,各式各样的问题都可以很自然地映射到这些数据结构上:Re ...

Could not load type System&period;ServiceModel&period;Activation&period;HttpModule解决办法

等注册完成后网站就可以打开了. win2008下提示未能从程序集“System.ServiceModel, Version=3.0.0.0问题解决 在Windows Server 2008中的IIS服 ...

DailyTick 开发实录 —— UI 设计

上次的文章中描述了 DailyTick 的设计理念.经过两周左右的设计和开发,现在 DailyTick 的主要 UI 已经完成了原型的设计和初步的实现.既然是原型,当然看起来就有点粗糙. 主 UI 主 ...

C&num; 读取系统日志

.NET框架类库提供了EventLog类和EventLogEntry类与系统日志进行交互二者属于System.Diagnostics命名空间 EventLog 类的属性主要有 Entris返回一个Ev ...

docker dead but pid file exists

CentOS 6安装docker 报docker dead but pid file exists 执行 yum install epel-release yum install docker-io ...

SpringBoot系列&colon; 集成MyBatis

本文主要修改自下面博客:http://www.ityouknow.com/springboot/2016/11/06/spring-boo-mybatis.htmlhttp://tengj.top/2 ...