![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/231061/image.png)
How to repeatedly subdive vector rectangle into equal area?
2 次查看(过去 30 天)
显示 更早的评论
I have x,y point data that I am trying to analyze in "swaths" or rectangles along the x-axis. I need to take the rectangle and subdivide it into equal areas, and count the number of areas that contain points.
Consider this example with a very small x,y dataset:
% make dummy data
x = 10*rand(30,1);
y = 25*rand(30,1);
% set up basis for building swaths or rectangles
step = 1;
xmax = step*(ceil(max(x)/step)); % rounds up to nearest step size
centers = 0:step/2:xmax;
centers = centers';
ymax = step*(ceil(max(y)/step));
%% will ultimately loop over all centers, but just consider one for now
jj = 2;
xr = [centers(jj-1), centers(jj+1), centers(jj+1), centers(jj-1), centers(jj-1)];
yr = [ymax+step, ymax+step, 0-step, 0-step, ymax+step];
% plot for discussion
scatter(x,y,'.')
box on
grid on
hold on
axis equal
plot(xr,yr,'blue')
% this gives point counts inside big rectangle as a sanity check
% but I need to do this for all the equal areas that I need to build
[in,on] = inpolygon(x,y,xr,yr);
num_points_example = sum(in);
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/230993/image.jpeg)
So in this example, my blue rectangle is 26x1 units. To make equal areas, I need to keep the x-axis length the same (=1) and break the rectangle down into 26 square along the y-axis. My question is, how do I programmatically build a grid or mesh for this 26x1? Note that in the next step, I'll be subdividing it further (52x2), etc., always maintaining equal area. I only have the basic Matlab, no fancy toolboxes, R2018a.
0 个评论
采纳的回答
Akira Agata
2019-7-24
编辑:Akira Agata
2019-7-24
The following is an example.
% make dummy data
x = 10*rand(30,1);
y = 25*rand(30,1);
% min. and max. of x and y
xlim = [-2 12];
ylim = [-5 30];
% Number of rectangle units (N x M)
N = 12;
M = 2;
% Count (x,y) points in each unit
xEdges = linspace(xlim(1),xlim(2),N+1);
yEdges = linspace(ylim(1),ylim(2),M+1);
figure
scatter(x,y,'.')
hold on
h = histogram2(x,y,xEdges,yEdges,...
'DisplayStyle', 'tile',...
'ShowEmptyBins','on',...
'FaceAlpha', 0.2);
for kk1 = 1:N
for kk2 = 1:M
text(mean(xEdges(kk1:kk1+1)),mean(yEdges(kk2:kk2+1)),...
num2str(h.Values(kk1,kk2)),...
'FontSize',12)
end
end
colorbar
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/231061/image.png)
3 个评论
Akira Agata
2019-7-24
编辑:Akira Agata
2019-7-24
Thank you for your comment, and sorry for causing some confusion.
This is because 'FaceAlpha' option was set to 0.2 to clearly show both scattered data points and numbers.
To keep consistency, please delete 'FaceAlpha' option, and move scatter(x,y,'.') after histogram2 (but in that case, the plot would be non-user-friendly...)
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Annotations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!