SpatialPointsDataFrame и вороной - PullRequest
Новая
       7

SpatialPointsDataFrame и вороной

0 голосов
/ 06 мая 2019

Я рисую карты в 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

Когда я прихожу к части вороной части в коде ниже. Я получаю ошибку ниже. Что вызывает это?

Ошибка в deldir :: deldir (xy [, 1], xy [, 2], rw = ext, eps = eps, suppressMsge = 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
...