mbq*_*mbq 16
这被称为凸壳问题; R内置chull
功能应该做的工作.要计算区域,您可以使用此处的公式.
编辑:更好; splancs
包有areapl
功能.所以解决问题的功能应如下所示:
cha<-function(x,y){
chull(x,y)->i
return(areapl(cbind(x[i],y[i])))
}
Run Code Online (Sandbox Code Playgroud)
例如:
library(splancs);
x<-rnorm(20);rnorm(20)->y;
#Some visualization
i<-chull(x,y);plot(x,y);polygon(x[i],y[i]);
#The area
cha(x,y);
Run Code Online (Sandbox Code Playgroud)