How to determine blocks in mesh grid?

2 次查看(过去 30 天)
Hello guys,
I have a matrix representing a grid, say
A = [-1 -1; -1 1; 1 -1; 1 1; -1 0; 0 -1; 0 0; 0 1; 1 0];
Every row of A is a Cartesian coordinate. Therefore if you plot a scatter, you'll find that A is a square. My question is: clearly there are four blocks in A: left upper, right upper, left lower, right lower. How can I write a code to generate 4 matrices which each of them represents one block? The result should be something like:
A1 = [-1 -1; 0 -1; -1 0; 0 0],
A2 = [0 -1; 1 -1; 0 0; 1 0],
A3 = [-1 0; 0 0; -1 1; 0 1],
A4 = [0 0; 1 0; 0 1; 1 1].
Thanks!
  2 个评论
Xiaohan Du
Xiaohan Du 2016-5-19
Hi Image Analyst,
Thanks for the reply!
I know what the answer should be like, and I can do it manually. But with the given matrix, what is the algorithm behind it? What if I have a finer grid? I cannot do it manually forever, it needs to be automatic.

请先登录,再进行评论。

采纳的回答

Guillaume
Guillaume 2016-5-19
编辑:Guillaume 2016-5-19
This should work:
A = [-1 -1; -1 1; 1 -1; 1 1; -1 0; 0 -1; 0 0; 0 1; 1 0];
%get the sorted unique x and y coordinates:
xy = cellfun(@unique, num2cell(A, 1), 'UniformOutput', false);
x = xy{1};
y = xy{2};
%create output
Agrid = cell(numel(y) - 1, numel(x) - 1);
%iterate through the grid and fill output
for xidx = 1 : numel(x) - 1
for yidx = 1 : numel(y) - 1
Agrid{yidx, xidx} = [repmat(x(xidx:xidx+1), 2, 1), repelem(y(yidx:yidx+1), 2)];
end
end
  3 个评论
Xiaohan Du
Xiaohan Du 2016-5-19
Hi Guillaume,
I just tested, this works for a finer grid as well! Thank you so much!
Guillaume
Guillaume 2016-5-19
This works for grids of any size, square or non-square and uniform or non-uniform.

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Sources 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by