2

我有三个SpatialPointsDataFrame对象，实际上每个对象只有一个点。我的目的是为它们中的每一个制作一个包含点的栅格，这样除了点之外的所有单元格都是“NA”，所以我可以使用distance()包中的函数raster生成一个栅格层，其中z 值是到 z 不是“NA”的唯一像元的距离。

我的代码对三个对象中的第一个对象没有问题，但其他两个对象出现以下错误：

    error in seq.default(zrng[1], zrng[2], length.out = cuts + 2) : 
    'from' cannot be NA, NaN or infinite
    In addition: Warning messages:
    1: In asMethod(object) :
      complete map seems to be NA's -- no selection was made
    2: In min(x) : no non-missing arguments to min; returning Inf
    3: In max(x) : no non-missing arguments to max; returning -Inf

我已经两次和三次检查了我的点是否包含在栅格范围内，但我真的无法确定问题所在

这是我的代码：

library(raster)

TIM <- data.frame()
TIM[1,1] <- -13.8309
TIM[1,2] <- 28.9942

VEN <- data.frame()
VEN[1,1] <- -15.7886
VEN[1,2] <- 27.8444

MCL <- data.frame()
MCL[1,1] <- -13.5325
MCL[1,2] <- 29.2914


coordinates(TIM) <- ~V1+V2
coordinates(VEN) <- ~V1+V2
coordinates(MCL) <- ~V1+V2

bb2 <- matrix(c(-20, -9.5, 20.5, 31.5), nrow = 2, ncol = 2, byrow = T)
bb2 <- extent(bb2)

r <- raster(nrows = 1217, ncols = 1047)
r <- setExtent(r, bb2, keepres=F)

rMCL <- rasterize(MCL, r)
spplot(rMCL)

#so far so good, but from now on it doesn't work

rVEN <- rasterize(VEN, r)
spplot(rVEN)


rTIM <- rasterize(TIM, r)
spplot(rTIM)

编辑：我尝试将其转换为 SpatialGridDataFrame 并绘制它，但我的观点不在栅格范围内，即绘图为空。代码：

rr <- as(rTIM, "SpatialGridDataFrame")
spplot(rr)
#this produces an empty plot

我还尝试在没有预定列数和行数的栅格中绘制它，它可以工作：

r <- raster()
r <- setExtent(r, bb2, keepres=F)
rTIM <- rasterize(TIM, r)
spplot(rTIM)
# this produces a raster containing my point

问题是，我真的需要设置绘图的分辨率，以便栅格的每个单元代表大约 1 平方千米，这就是我之前使用的行数和列数。有任何想法吗？

4

1 回答 1

0

我可以通过将所有三组坐标添加到同一个数据框，然后在创建栅格时使用 count 函数来使其工作：

library(raster)

# Add all 3 sets of coordinates to the same dataframe

df <- data.frame()
df[1,1] <- -13.8309
df[1,2] <- 28.9942
df[2,1] <- -15.7886
df[2,2] <- 27.8444
df[3,1] <- -13.5325
df[3,2] <- 29.2914

# Create new column in dataframe (we will use this for the count function)

df$x <- c(1,1,1)

# Convert to spatial points dataframe

df.sp <- df
coordinates(df.sp) <- ~ V1+V2

# Make raster

bb2 <- matrix(c(-20, -9.5, 20.5, 31.5), nrow = 2, ncol = 2, byrow = T)
bb2 <- extent(bb2)

r <- raster(nrows = 1217, ncols = 1047)
r <- setExtent(r, bb2, keepres=F)

# Rasterise using the count function

raster <- rasterize(df.sp, r, field= "x", fun="count")

# The table shows there are 3 cells with a value of 1 so it seems to have worked

table(values(raster))

1 
3 

spplot(raster)

raster

class       : RasterLayer 
dimensions  : 1217, 1047, 1274199  (nrow, ncol, ncell)
resolution  : 0.01002865, 0.00903862  (x, y)
extent      : -20, -9.5, 20.5, 31.5  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
data source : in memory
names       : layer 
values      : 1, 1  (min, max)

这个情节对我来说没有多大意义，但我认为那是因为你的栅格中有很多单元格，所以你什么都看不到。栅格中的值绝对是 3 个值为 1 的像元，其余的都是 NA，所以我认为这是您想要的。

于 2016-06-09T22:04:43.403 回答