如何光栅化而不是在r中获取重复信息答案

【问题标题】：How to rasterise and not to get duplicate information in r如何光栅化而不是在r中获取重复信息
【发布时间】：2015-05-30 00:18:34
【问题描述】：

我在使用 R 中的 raster 包对 data.frame 进行光栅化时遇到问题。

我的 data.frame 包含世界海洋的环境数据（温度等）及其坐标（网格 0.5*0.5），十进制经度从 -90 到 90，十进制纬度从 -180 到 180。所以基数包含 90 *2*2 x 180*2*2 = 360*720 = 259200 行，59 个冒号（57 个变量 + 2 个冒号的坐标）。

光栅化后，这是我用plot(r) 得到的：http://postimg.org/image/rqocxcbi3/

所以，重复的图像，方向错误。

我的代码是：

FILE_ENV = read.csv('ENV_DATABASE.csv')
coordinates(FILE_ENV) <- ~LON+LAT
proj4string(FILE_ENV3)=CRS("+init=epsg:4326")
FILE_ENV = spTransform(FILE_ENV,CRS("+init=epsg:4326"))
gridded(FILE_ENV) = TRUE
r = raster(FILE_ENV)

plot(r)

谁能看到我在这里错过/搞砸了什么？

谢谢

编辑： head(FILE_ENV)

LON LAT BAT BAT_CLASSE SLOPE SEDIMENT SST SST_SEAFLOOR SST_SUM SST_WIN SAL_SURF SAL_SEAFLOOR...

1 -179.75 89.75 2804 NA 0.14031838 NA NA NA NA NA NA NA

2 -179.25 89.75 2941 NA 0.12495525 NA NA NA NA NA NA NA

3 -178.75 89.75 3048 NA 0.07784129 NA NA NA NA NA NA NA

4 -178.25 89.75 3093 NA 0.03123910 NA NA NA NA NA NA NA

5 -177.75 89.75 3109 NA 0.01536359 NA NA NA NA NA NA NA

6 -177.25 89.75 3063 NA 0.15619729 NA NA NA NA NA NA NA

【问题讨论】：

听起来你需要在某个地方进行转置，我不知道标准 t() 是否适用于光栅对象，但值得一试。
谢谢！我试过了，但我仍然得到一个重复的情节..
spTransform 用于进行坐标转换——你不需要。我的猜测是 raster 构造函数不理解 FILE_ENV 对象的坐标。将其分解为更明确的部分以供 raster

标签： r environment raster rasterize

【解决方案1】：

由于您没有显示head(FILE_ENV)（我的代码中为d），因此很难说您应该如何做到这一点。但从本质上讲，您应该能够执行以下操作：

library(raster)
d <- read.csv('ENV_DATABASE.csv')
r <- rasterFromXYZ(d, crs="+proj=longlat +datum=WGS84")

【讨论】：

谢谢！我编辑了我的第一条消息以添加head(FILE_ENV)。使用您的代码改进了结果，但我仍然遇到同样的问题（重复和旋转的图像）：s27.postimg.org/6krcblg43/Rplot2.png
我完全重新创建了我的数据库的坐标冒号，它现在可以工作了；我仍然不知道问题出在哪里，但显然是在数据库中

【解决方案2】：

我认为 rasterfromXYZ 会起作用。也就是说，似乎在使用 raster 将事情分解成更小、更明确的块时效果最好。似乎构造函数假设您想要从对象中得到不同的东西，并且它正在丢失矩阵条目的坐标。

这不是最佳做法，但由于您知道您的数据是一个严格的网格，我们可以将条目转换为点，然后将它们栅格化：

library(raster)
coords <- coordinates(FILE_ENV[,c("LON","LAT")])
spdf <- SpatialPointsDataFrame(coords,data=FILE_ENV)
r <- rasterize(spdf,field='value',raster(nrow=360,ncol=720))

现在 r 应该包含您的栅格。

【讨论】：