示例如下，自行延伸： 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), 来自: 180.218.154.163 ※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1446217622.A.D0E.html