R中的空间数据：多类SVM的图决策区域

Question

我有data.frame个x和y列，以及class列，它给出了现有多类SVM模型下每个点的分类。这是一些示例代码：

library(rgdal)
library(rgeos)
library(e1071)  # for svm()
library(sp)
library(raster)
library(maptools)
library(plyr)

## Create a mask of the data region, as a data frame of x/y points.
covering <- function(data, xlen=150, ylen=150) {
    # Convex hulls of each class's data points:
    polys <- dlply(data, .(class), function(x) Polygon(x[chull(x[-3]), -3]))
    # Union of the hulls:
    bbs <- unionSpatialPolygons(SpatialPolygons(list(Polygons(polys, 1))), 1)

    # Pixels that are inside the union polygon:
    grid <- expand.grid(x=seq(min(data$x), max(data$x), length.out=xlen),
                        y=seq(min(data$y), max(data$y), length.out=ylen))
    grid[!is.na(over(SpatialPoints(grid), bbs)), ]
}

set.seed(123)
data <- rbind(data.frame(x=rnorm(1000, 5*runif(1)), y=rnorm(1000, 5*runif(1)), class='A'),
              data.frame(x=rnorm(1000, 5*runif(1)), y=rnorm(1000, 5*runif(1)), class='B'),
              data.frame(x=rnorm(1000, 5*runif(1)), y=rnorm(1000, 5*runif(1)), class='C'))

m <- svm(class ~ x+y, data)
grid <- covering(data)

par(mfrow=c(1,2))
plot(y ~ x, data, col=data$class)
plot(y ~ x, grid, col=predict(m, grid), pch=20)

我接下来要做的是将此结果转换为某种SpatialPolygons对象，每个因子级别有一个Polygons个对象，以便导出到GeoJSON，以便可以在映射中使用应用。这样做的好方法是什么？我是否需要自己编写例程来跟踪图像（作为矩阵）并找到区域之间的边界？

我查看了rasterToPolygons()的文档，但我无法弄清楚如何将它应用到我的情况中，所以我欢迎一些帮助。

最后，我的数据 将成为具有真实纬度/经度信息的地理空间，但我想先尝试这个更简单的案例。

Answer 1

您需要将第二张图片转换为栅格（请参阅here）。

sp.grid <- cbind(grid, col = predict(m, grid))
coordinates(sp.grid) <- ~ x + y
gridded(sp.grid) <- TRUE
sp.grid <- raster(sp.grid)

然后你的尝试有效。

grid.pols <- rasterToPolygons(sp.grid, n = 16, dissolve = TRUE)
plot(grid.pols)

class       : SpatialPolygonsDataFrame 
features    : 3 
extent      : -1.842044, 7.259169, -2.298892, 5.512507  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
variables   : 1
names       : col 
min values  :   1 
max values  :   3