R提取栅格功能不起作用

时间:2013-09-12 20:59:44

标签: r extract raster

我在光栅包中的提取功能遇到问题。栅格不是空的,但是当我执行提取功能时,所有值都返回为NA。我可以将光栅上的点绘制得很好,它们会显示在同一位置。

library(dismo)
library(raster)

#Import CSV File
file <- paste("C:/Michelle/Michelle/R/CuskPA_Sept12_2013_1.csv", sep=" ")

#read file
cusk <- read.table(file, header=TRUE, sep=',')
cusks<- subset(cusks, SlongDD!=0 & SlatDD!=0)
cuskll <-cusks[,24:25]

#import raster
salvar <- raster("C:/Michelle/Michelle/R/Variables/salvar")
ext <-extent(-67.6, -64.9, 41.0, 44.7453)
salvarc <-crop(salvar, ext)

#extract from raster   
psalvar <- extract(salvarc, cuskll)

这是我完成提取后得到的结果

head(psalvar)
[1] NA NA NA NA NA NA

summary(psalvar)
Mode    NA's 
logical    3859 

栅格不会显示为空。

> salvar
class       : RasterLayer 
dimensions  : 334, 240, 80160  (nrow, ncol, ncell)
resolution  : 0.01121539, 0.01121539  (x, y)
extent      : -67.59662, -64.90493, 41.00048, 44.74642  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs 
data source : in memory
names       : salvar 
values      : 0.01655413, 2.137321  (min, max)

以下是文件

中纬度/经度的示例
head(cuskll)
SlatDD   SlongDD
1 44.41533 -66.07883
2 44.41167 -66.07733
3 44.44883 -66.07433
4 44.47067 -66.01217
5 44.47050 -66.01183
6 44.47633 -66.01450

1 个答案:

答案 0 :(得分:3)

先经度和纬度第二