时间戳

首页 > 解决方案 > 用纬度和海拔而不是纬度和经度绘制栅格

r - 用纬度和海拔而不是纬度和经度绘制栅格

问题描述

我正在使用 R 并且我有一些问题来解决这个问题：我有 2 个栅格（相同的区域、相同的分辨率、相同的范围、相同的 crs）：

栅格 A：（纬度、经度、值_A）

class      : RasterLayer
dimensions : 832, 541, 450112  (nrow, ncol, ncell)
resolution : 0.008333333, 0.008333333  (x, y)
extent     : -122.2458, -117.7375, 35.0625, 41.99583  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs
names      : values_A
values     : -753.4842, 0  (min, max)

栅格 B：（纬度、经度、海拔）

class      : RasterLayer
dimensions : 832, 541, 450112  (nrow, ncol, ncell)
resolution : 0.008333333, 0.008333333  (x, y)
extent     : -122.2458, -117.7375, 35.0625, 41.99583  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs
names      : dem
values     : 40, 4176  (min, max)

然后，我需要获取作为高程和纬度函数的栅格 A 图（x= 高程，y= 纬度，像素值 = values_A）。

此外，我可能需要聚合具有相同纬度和相同海拔的像素，例如使用均值函数。

你对如何用 R 解决这个问题有什么建议吗？

谢谢！

标签： rcoordinatesrastercoordinate-systems

解决方案

您可以像这样设置示例数据

library(raster)
A <- B <- raster(nrow=83, ncol=54, ext=extent(-122.2458, -117.7375, 35.0625, 41.99583))
values(A) <- rep(1:nrow(A), ncol(A))
values(B) <- 1:ncell(B)

这是一个解决方案

x <- cbind(elevation=values(B), latitude=yFromRow(B, 1:nrow(B)), values(A))
r <- rasterFromXYZ(x)

或者使用一些现有数据

B <- getData("alt", country='CHE')
A <- init(B, "y")
x <- na.omit(cbind(elevation=values(B), latitude=yFromRow(B, 1:nrow(B)), A=values(A)))

你可能想做一些四舍五入

x[,1] <- round(x[,1], -2)
x[,2] <- round(x[,2], 1)

r <- rasterFromXYZ(x)
plot(r, asp=NA)

推荐阅读