我正在尝试从包含列中物种名称、经度和纬度的矩阵创建栅格图层。
数据名称为“base”,当我尝试使用base
命令创建范围为 的栅格图层时
r<-raster(base)
我收到一个错误:
Error en .local(x, ...) : list has no "x"
你能帮我吗?谢谢
你没有给我们太多继续,但你可能想试试这个(下面的例子假设你的坐标数据是使用 WGS84 数据的经纬度格式。如果它不是你的点在 2- D空间会出错,所以请尝试查找您的投影信息...
require( sp )
require( raster )
require( maptools )
#Some dummy data in the format you have described
base <- data.frame( longitude = sample( c(-20:20) , 10 ) , latitude = sample( c(-30:30) , 10 ) , species = sample( letters[1:3] , 10 , repl = TRUE ) )
head( base )
# longitude latitude species
# 1 -10 7 a
# 2 3 14 b
# 3 -9 29 b
# 4 -8 -25 a
# 5 0 -19 b
# 6 -19 19 c
#Make spatial points dataframe, which will have a spatial extent
spdf <- SpatialPointsDataFrame( base[ c("longitude" , "latitude") ] , data = data.frame( base$species ) , proj4string = CRS("+proj=longlat +datum=WGS84") )
#make raster based on the extent of your data
r <- raster( extent( spdf ) )
r
# class : RasterLayer
# dimensions : 10, 10, 100 (nrow, ncol, ncell)
# resolution : 3.3, 5.8 (x, y)
# extent : -20, 13, -29, 29 (xmin, xmax, ymin, ymax)
# coord. ref. : NA
# values : none