Divide a region into rectangles based on density
Given a set of (x,y) points, partition the field into rectangles each containing a minimum number of points
RectDense(INx, INy, icut = 1, u = par("usr"), ndivs = 10)INx |
x-coordinates |
INy |
y-coordinates |
icut |
cut off for number of points |
u |
user coordinates |
ndivs |
number of divisions in x-coordinate |
Based on the user coordinates as returned from par('usr'). Each rectangular region is tested for the number of points that fall within icut or greater.
List:
icorns |
matrix of corners that passed test |
ilens |
vector,number of points in each icorns box |
ipass |
vector, index of the corners that passed icut |
corners |
matrix of all corners |
lens |
vector,number of points for each box |
Jonathan M. Lees<jonathan.lees@unc.edu>
x = rnorm(100)
y = rnorm(100)
plot(x,y)
u = par('usr')
RI = RectDense(x, y, icut=3, u=u, ndivs=10)
rect(RI$icorns[,1],RI$icorns[,2],RI$icorns[,3],RI$icorns[,4], col=NA, border='blue')Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.