Take the 2-minute tour ×
Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

I have a database which contains Lat/Long pairs to identify the location of points of interest. I would like to group the points of interest into groups of 10. The group should be geographically local and contain exactly 10 points. Each group should be of a minimum area.

I've looked at various implementations in R but none of them (that I can see) allow you to specify a definite cluster size.

I'd also appreciate a pointer to any beginner level R tutorials as it's all a bit foreign to me at the moment.

(I previously asked this question Grouping map points into fixed cluster sizes but I don't think I was exact enough in my question to get a good answer.)

EDIT:
Geographically local - I think I mean that groups shouldn't significantly overlap. In my application (allocating people to groups for monitoring purposes) it would be ideal if each group was as small as possible in physical area.
Minimum area - again, trying to keep the group area to a minimum. I suppose this could be quantified as keeping the each group's area below a specified threshold (to avoid dozens of small groups and one large one).

share|improve this question
2  
It would help to be a bit more precise about what you're looking for. How would you quantify "geographically local" and "minimum area"? –  whuber Jul 10 '12 at 21:48
    
In some situations, won't "minimum area" and "exactly 10 points" be mutually exclusive? How can you expect both to be used in, for example, you have a 1 mile "minimum area", and 1 feature that has no other features within 10 miles (an outlier, possibly)? –  RyanDalton Jul 12 '12 at 14:26
    
I expect there to be some exceptions, but these can be manually dealt with. I have a fixed size data set and don't mind manually fiddling a few of the groups, but I'd really like an automatic solution for the rest! :) –  Graeme Hilton Jul 13 '12 at 15:26
    
Related: gis.stackexchange.com/questions/15906/…. –  whuber Dec 27 '12 at 22:03
    
The description is still too "hand-wavy" to be solveable in an automated way. Can you make the data public? –  BradHards Jan 18 '13 at 4:23

1 Answer 1

I think you might be looking for a k-nearest neighbor tool. This type of tool can be used to identify the 10 nearest neighbors of all points in your dataset. There seem to be a few different options for this (with some using different algorithms or having slightly different functionalities), and I'm not sure which would be the best option. But here are a few links:

http://stat.ethz.ch/R-manual/R-patched/library/class/html/knn.html http://en.wikibooks.org/wiki/Data_Mining_Algorithms_In_R/Classification/kNN

You may need to combine the results with a clustering algorithm or cluster ensemble tool to identify clusters of points that have similar sets of neighbors to get at classifications of points where there is little-to-no overlap. You may need to do a little bit of manual fiddling with the output, but it should allow you to automate a large portion of the work

Some links: http://jmlr.csail.mit.edu/papers/volume3/strehl02a/strehl02a.pdf http://cran.r-project.org/web/packages/clue/vignettes/clue.pdf

You might also be able to find a k-means clustering tool that will do this all in one step and enforce the 10 point in a cluster rule (just divide the total number of points by 10 and select that as the number of desired clusters for the tool).

share|improve this answer
    
Also, I've found the YaRi tutorial helpful for R: ahandel.myweb.uga.edu/resources.htm –  Jezibelle Jan 23 '13 at 17:18

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.