我正在R Statistics中绘制地图。 我有一个像下面的SpatialPointsDataFrame'dsp'。
class : SpatialPointsDataFrame
features : 905
extent : 3605539, 4471451, 1886482, 2428298 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0
variables : 4
names : res_id, alert, Longitude, Latitude
min values : 3, 1, -122.299294, -5.675994
max values : 4095, 8, 44.474445, 68.8162
和区域定义“ grootepeel”
class : SpatialPolygonsDataFrame
features : 13
extent : 3605539, 4471451, 1886482, 2428298 (xmin, xmax, ymin, ymax)
coord. ref. : NA
variables : 1
names : Id
min values : 0
max values : 0
当我进入下面代码中的voronoi部分时。我得到下面的错误。是什么原因造成的？
deldir :: deldir（xy [，1]，xy [，2]，rw = ext，eps = eps，preventMsge = TRUE）中的错误： 给定的矩形窗口内没有点 且未指定虚拟点。因此有 没有三角剖分/细分的点。
grootepeel <- readShapeSpatial("Alpine_Convention_Perimeter_2018.shp")
proj4string(grootepeel) <- "+proj=longlat +datum=WGS84" # specify projection
plot(grootepeel)
bbpeel <- bbox(grootepeel)
dsp <- SpatialPoints(dfpeel[,2:3], proj4string = CRS("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"),bbox = bbpeel)
d <- dfloc
d <- as.data.frame(d)
dsp <- SpatialPointsDataFrame(dsp,d)
v <- voronoi(dsp,ext=extent(grootepeel))
plot(v)
d包含
res_id alert Longitude Latitude
1 3 2 10.665000 44.145599
2 4 5 13.345600 47.563801