R：填充轮廓中的封闭区域

Question

我在R中绘制一个knn-classification图，使用绘图绘制样本和轮廓以绘制分类平面的线。

这是我的代码：

k<-1
datax<-rbind(matrix(rnorm(30,-1,5.25),15,2),matrix(rnorm(36,1,5.25),18,2))
datay<-rbind(matrix(1,15,1),matrix(0,18,1))
plot(datax[,1], datax[,2],pch = datay+1,axes=FALSE,ann=FALSE)
box()
n <- 1000
xp <- seq(length=n, from = min(datax[,1]), to = max(datax[,1]))
yp <- seq(length=n,from = min(datax[,2]) ,to = max(datax[,2]))
gr <- expand.grid(xp, yp)
library(class)
z <- as.numeric(knn(datax, gr, datay,k))-1
zM <- matrix(z, n, n, byrow = FALSE)
contour(xp, yp, zM, xlab="x",ylab="",nlevels = 1 ,lwd=2, add=TRUE, drawlabels =FALSE)

我的问题是：如何为情节中的封闭区域着色？我试过fill.contour，但没有添加参数。我只想要分类器为= 0白色的区域和蓝色分类= 1的区域。我该怎么办？ 感谢

Answer 1

您可以使用contour来保留轮廓线边缘的坐标，而不是contourLines，而不是polygon。

plot(datax[,1], datax[,2],axes=FALSE,ann=FALSE, type="n")
box()
cL <- contourLines(xp, yp, zM,nlevels = 1)
lapply(cL,function(x)polygon(x$x,x$y,col="red"))
points(datax[,1], datax[,2],pch = datay+1)

然而，到达绘图边缘的轮廓线并不完美（请参见第二个图的左下角），因此需要进行一些手工调整：

编辑：在嵌套轮廓线的情况下，我认为没有一种简单的方法来处理它，但这是一种方法：

library(splancs)
ord <- sapply(lapply(cL,function(x)datay[inout(datax,cbind(x$x,x$y))]),
              median) #Check what values are present in the polygon and 
                      #take the most common one
plot(datax[,1], datax[,2],axes=FALSE,ann=FALSE, type="n")
box()
lapply(cL[ord==1],function(x)polygon(x$x,x$y,col="blue"))
lapply(cL[ord==0],function(x)polygon(x$x,x$y,col="white"))
points(datax[,1], datax[,2],pch = datay+1)

第二次编辑：当然也有可能在你的情况下使用函数image：

image(xp, yp, zM, col=c("transparent","blue"))
points(datax[,1], datax[,2],pch = datay+1)