RFOC (version 1.0-4)

RectDense: Divide a region into rectangles based on density

Description

Given a set of (x,y) points, partition the field into rectangles each containing a minimum number of points

Usage

RectDense(INx, INy, icut = 1, u = par("usr"), ndivs = 10)

Arguments

INx

x-coordinates

INy

y-coordinates

icut

cut off for number of points

u

user coordinates

ndivs

number of divisions in x-coordinate

Value

• 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

Details

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.

Examples

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')