我正在尝试创建一个大的(空间?)数据帧,用于后续多变量统计分析(例如,PCA)中的一组60个光栅文件。
我对使用R中的空间数据有些新意,并且不了解空间坐标数据是如何在数据文件中保留或强制执行的。
我使用以下代码读取光栅文件并从中创建堆栈:
#load libraries
library(raster)
library(sp)
library(rgdal)
library(spatial.tools)
#set wd
setwd("C:/Users/...../data/")
#get raster files from wd
files <- list.raster.files(path = getwd(), pattern = ".tif$", recursive =FALSE, return_rasters = FALSE, return_bbox = FALSE)
#create rasterStack
mystack <- stack(files$raster_files)
#read metadata / summary of rasterStack
mystack
#output has "NA" for coordinate reference system.
为什么R在导入文件时丢失了spaital参考信息? 当我在ArcGIS中查看光栅文件时,他们肯定会有投影信息。
我尝试以这种方式手动应用投影:
#first define the projection using the proj4 syntax
projection <- "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs"
#then apply the projection to the rasterStack
projectRaster(files, projection, method = bilinear, filename = "STACKprj")
这给了我一个错误,即输入数据(&#34; mystack&#34;)有NA作为参考系统,即它没有启动参考系统就无法工作
#I try instead to apply the projection in this way:
STACKprj <- projection(mystack)
这不起作用。
任何提示?谢谢!
答案 0 :(得分:4)
你应该尝试这样的事情:
library(raster)
library(sp)
library(rgdal)
#work wd
wd <- "C:/Users/...../data"
#get raster files from work wd
files <- list.raster.files(path = wd, pattern = ".tif$",
recursive =FALSE, return_rasters = FALSE, return_bbox = FALSE)
#create rasterStack
mystack <- stack(files$raster_files)
proj4string(mystack) <- CRS("+init=epsg:3857") # OSM Mercator projection