voronoi2dCellColour.m File Reference

Go to the source code of this file.

Functions
and colour the % cells with the value of the orientational correlations defined at each %	site (col=1)
and colour the % cells with the value of the orientational correlations defined at each % or with the area of the	cell (col=0). % x
id	Returns:Testcase: (id variance, [Testcase] -1, 1, 20 linspace)
%	voronoi2dCellColour (kron(x', ones(length(y), 1)), kron(ones(length(x), 1), y'), 0.75, zeros(length(x)^ 2, 1), 1, 0)
% If at least one of the indices % then it is an open region and we can t patch that	area (jj)
	p6cp6 (jj)
% use area as color end if	edgeCol (jj)< 0 h.EdgeColor
	get (h)
	plot (x, y, 'r *')
	plot (x, y, 'bo')
colorbar	set (gca, 'FontName', 'Latin Modern Roman', 'FontSize', 22)
	set (gca, 'TickLabelInterpreter', 'latex')
	ylabel (' $y$(m)', 'Interpreter', 'latex')
	xlabel (' $x$(m)', 'Interpreter', 'latex')

Variables
	function [p6cp6, area, avg_area, num_edges, var]
and colour the % cells with the value of the orientational correlations defined at each % or with the area of the	y =x
and colour the % cells with the value of the orientational correlations defined at each % or with the area of the	edgeCol ==0 is black
and colour the % cells with the value of the orientational correlations defined at each % or with the area of the edgeCol is white %	dx = 1
%	radius = 0.5
	q = ones(length(x)*length(y),1)
	DT = delaunayTriangulation([x y])
if	colScheme
end	num_edges =zeros(length(c),1)
	avg_area = 0
	area =zeros(length(c),1)
	clf
	R = DT.Points
% Credit to MathWorks support team for neighboringVertices m	vTriAtt = vertexAttachments(DT)
for	ii
% Find all the unique vertices and remove the current vertex	neighboursOfInternal {ii} = setdiff(unique(verticesOfTI), ii)
end % Calculate the Voronoi diagram for	jj
% If at least one of the indices	is
	h = patch(v(c{jj},1),v(c{jj},2),p6cp6(jj))
h	LineWidth = '2'
h	LineStyle = '-'
	var = sum((area(area>1) - avg_area).^2)
hold	on
axis	square

Function Documentation

area()

% If at least one of the indices % then it is an open region and we can t patch that area

(

jj

)

cell()

and colour the % cells with the value of the orientational correlations defined at each % or with the area of the cell

(

col

= 0

)

edgeCol()

% use area as color end if edgeCol

(

jj

)

get()

get

(

h

)

p6cp6()

p6cp6

(

jj

)

plot() [1/2]

plot	(	x	,
		y	,
		'r *'
	)

plot() [2/2]

plot	(	x	,
		y	,
		'bo'
	)

Returns:Testcase:()

id Returns:Testcase: ( id  variance, [Testcase] -  1, 1  , 20  linspace  )

virtual

set() [1/2]

colorbar set	(	gca	,
		'FontName'	,
		'Latin Modern Roman'	,
		'FontSize'	,
		22
	)

set() [2/2]

set	(	gca	,
		'TickLabelInterpreter'	,
		'latex'
	)

site()

and colour the % cells with the value of the orientational correlations defined at each % site

(

col

= 1

)

voronoi2dCellColour()

% voronoi2dCellColour	(	kron(x', ones(length(y), 1))	,
		kron(ones(length(x), 1), y')	,
		0.	75,
		zeros(length(x)^ 2, 1)	,
		1	,
		0
	)

xlabel()

xlabel	(	' $x$(m)'	,
		'Interpreter'	,
		'latex'
	)

ylabel()

ylabel	(	' $y$(m)'	,
		'Interpreter'	,
		'latex'
	)

Variable Documentation

area

area =zeros(length(c),1)

Definition at line 32 of file voronoi2dCellColour.m.

avg_area

end end avg_area = 0

Definition at line 31 of file voronoi2dCellColour.m.

clf

clf

Definition at line 34 of file voronoi2dCellColour.m.

colScheme

if colScheme

Initial value:

==1
    p6cp6=zeros(length(c),1)

Definition at line 27 of file voronoi2dCellColour.m.

DT

DT = delaunayTriangulation([x y])

Definition at line 25 of file voronoi2dCellColour.m.

dx

dx = 1

Definition at line 10 of file voronoi2dCellColour.m.

edgeCol

and colour the % cells with the value of the orientational correlations defined at each % or with the area of the edgeCol ==0 is black

Definition at line 10 of file voronoi2dCellColour.m.

function

function[p6cp6, area, avg_area, num_edges, var]

Initial value:

=voronoi2dCellColour(x,y,radius,edgeCol,dx,colScheme)
%vorCellColour Determine the Voronoi diagram of the input data

Definition at line 1 of file voronoi2dCellColour.m.

h

% use g6 as color else h = patch(v(c{jj},1),v(c{jj},2),p6cp6(jj))

Definition at line 58 of file voronoi2dCellColour.m.

ii

for ii

Initial value:

= 1:size(R,1)
    % 2. Use the connectivity list to get the vertex indices of all these
    % triangles
    verticesOfTI         = DT.ConnectivityList(vTriAtt{ii},:)

Definition at line 39 of file voronoi2dCellColour.m.

is

% If at least one of the indices is

Definition at line 50 of file voronoi2dCellColour.m.

jj

end % Calculate the Voronoi diagram for jj

Initial value:

= 1:length(c)
    %Credit to MathWorks team for voronoi examples. https:
    if (all(c{jj}~=1) && all(sqrt((x(jj)).^2 + (y(jj)).^2) < radius))

Definition at line 48 of file voronoi2dCellColour.m.

LineStyle

h LineStyle = '-'

Definition at line 65 of file voronoi2dCellColour.m.

LineWidth

: LineWidth = '2'

Definition at line 64 of file voronoi2dCellColour.m.

neighboursOfInternal

% Find all the unique vertices and remove the current vertex neighboursOfInternal {ii} = setdiff(unique(verticesOfTI), ii)

Definition at line 44 of file voronoi2dCellColour.m.

num_edges

end num_edges =zeros(length(c),1)

Definition at line 30 of file voronoi2dCellColour.m.

on

hold on

Definition at line 82 of file voronoi2dCellColour.m.

q

q = ones(length(x)*length(y),1)

Definition at line 22 of file voronoi2dCellColour.m.

R

R = DT.Points

Definition at line 35 of file voronoi2dCellColour.m.

radius

% radius = 0.5

Definition at line 22 of file voronoi2dCellColour.m.

square

axis square

Definition at line 85 of file voronoi2dCellColour.m.

var

var = sum((area(area>1) - avg_area).^2)

Definition at line 79 of file voronoi2dCellColour.m.

vTriAtt

% Credit to MathWorks support team for neighboringVertices m vTriAtt = vertexAttachments(DT)

Definition at line 38 of file voronoi2dCellColour.m.

y

y =x

Definition at line 10 of file voronoi2dCellColour.m.