library(sp) data(meuse.grid) coordinates(meuse.grid) <- ~x+y plot(meuse.grid) XY = coordinates(meuse.grid) xy = XY[sample(nrow(XY),10),] library(magrittr) ExtractbyCircle = function(xy, XY, con){ d = as.matrix(dist(rbind(xy, XY)))[-1,1] %>% set_names(1:nrow(XY)) s = d[d<=con] return(as.integer(names(s))) } SS = list() for(i in 1:nrow(xy)){ SS[[i]] = ExtractbyCircle(xy=xy[i,], XY=XY, con = 100) } str(SS) XY[SS[[1]],] ※ 引述《celestialgod (天)》之铭言： : 示例如下，自行延伸： : library(pracma) : # circle : circleFunc = function(r, x0 = 0, y0 = 0, n.split = 1000){ : cbind(r*cos(seq(0, 2*pi, length.out = n.split)) + x0, : r*sin(seq(0, 2*pi, length.out = n.split)) + y0) : } : squreArea = matrix(runif(1000), 100) : circlePoints = circleFunc(0.5, 0.5, 0.5) : inCicle = inpolygon(squreArea[,1], squreArea[,2], : circlePoints[,1], circlePoints[,2]) : plot(squreArea[,1], squreArea[,2], : xlim = c(0,1), ylim = c(0,1), xlab = "x", ylab = "y") : lines(circlePoints[,1], circlePoints[,2], col = "blue") : points(squreArea[inCicle,1], squreArea[inCicle,2], col = "red") : http://i.imgur.com/H1vyxBa.png : 我之前是365万个点，是用Rcpp逐点处理的 (距离矩阵要开满大的满麻烦的) : 如果点不多的话，用inpolygon很快XD : ※ 引述《allen63521 (GoGoPadres)》之铭言： : : [问题类型]: : : 程式谘询(我想用R 做某件事情，但是我不知道要怎麽用R 写出来) : : [软体熟悉度]: : : 新手(没写过程式，R 是我的第一次) : : [问题叙述]: 想要分析位在某些目标附近的资料 : : 手上有一笔关於不动产成交的资料，内容包含成交价格、坪数、该房屋位置的横、纵座标 : : 我想要用散布图把所有data的横、纵座标在平面座标上排出来，并且把几个目标定位在座标上，去讨论这几个目标附近（可能是一个圆或是矩形范围内）的房屋成交价 : : 目前的想法是： : : 如果想讨论的范围是圆，可以对"该目标跟所有data"去算点和点之间的距离 : : 再从所有算出来的距离中挑比较近（抓某个比例）的那些房子出来观察 : : 但这个方法在目标不只有一个的时候可能就会变得比较麻烦，因为要把"每个目标"和"每栋房子"的组合都算过 : : 想要请问有没有可能做到： : : 以每个目标为圆心，设定一个圆（先不考虑怎麽决定半径大小），让系统去判定有哪些房子是落在这些范围里面的，就可以直接把他们抓出来看 : : [关键字]: : : 范围内的资料 --

※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 1.164.12.185 ※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1446220015.A.909.html ※ 编辑: Edster (111.249.20.219), 10/31/2015 08:14:06