2016-03-01 1 views

ответ

3
library(raster) 
library(rgdal) 

files <- list.files(path=".", pattern="tif$", full.names=TRUE) # select all the tiff files in the directory 
s <- stack(files) # stack all of them using r raster library 
shp <- shapefile("points.shp") # read the point shape file 
ext <- extract(s, shp) # extract the values from raster stack. 
+0

@RobertH Я не хочу складывать слои, мне нужно пройти через каждый из них и извлечь значения, используя файл формы точки. – jmutua

+0

@jmutua вы не получите разные результаты с циклом for. вы можете использовать: 'mylist <- list() для (i в файлах) {r <- растровый (i) ext <- extract (r, loc) dat <- data.frame (as.vector (ext)) mylist [[i]] <- dat } ' –

+0

@ Geo-sp - это не обязательно так. Вы предполагаете, что все растры имеют одинаковую степень и разрешение (т. Е. Они стекируются) – jbaums

 Смежные вопросы

  • Нет связанных вопросов^_^