0

我正在做一个足够简单的操作，将栅格从 WGS84 重新投影到英国国家网格，但我想知道重新投影后的一些结果。结果栅格的总和完全不同；这是由于分辨率和“projectRaster”中双线性插值的工作方式吗？

我从一个 csv 开始，该 csv 具有覆盖纬度和经度（仅整数值）范围为 -180 到 180 度的全局数据，以及一些 z 值。这是子集，制作成投影WGS84的栅格并转换为BNG（下面的子集）：

x <- rep(c(-10:3), times = 10)
n <- 14
y <- rep(48:57, each=n)
z <- rnorm(n=140, mean=20, sd=5)
ind <- which(z %in% sample(z, 45))
z[ind]<-NA
df <- data.frame("x"=x,"y"=y,"value"=z)

bng <- '+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +datum=OSGB36 +units=m +no_defs'
wgs84 <- '+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs'

coordinates(df) = ~x + y
ras = rasterFromXYZ(df, crs=wgs84)
cellStats(ras,sum)
plot(ras)

ras_bng <- projectRaster(ras,crs=bng)
plot(ras_bng)
cellStats(ras_bng,sum)

总和从 1919 到 2625（在我的情况下），相当大。

是否纯粹来自在“边缘”周围创建如此多额外单元的重投影？如果我要重新投影到分辨率小得多（5公里）的栅格，这会大大减少不同的总和吗？

谢谢，小号

> ras
class       : RasterLayer 
dimensions  : 10, 14, 140  (nrow, ncol, ncell)
resolution  : 1, 1  (x, y)
extent      : -10.5, 3.5, 47.5, 57.5  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0 
data source : in memory
names       : layer 
values      : 3.243918, 32.21532  (min, max)

> ras_bng
class       : RasterLayer 
dimensions  : 12, 19, 228  (nrow, ncol, ncell)
resolution  : 67900, 111000  (x, y)
extent      : -375677.8, 914422.2, -343553.2, 988446.8  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +datum=OSGB36 +units=m +no_defs +towgs84=446.448,-125.157,542.060,0.1502,0.2470,0.8421,-20.4894 
data source : in memory
names       : layer 
values      : 8.174968, 28.31331  (min, max)
4

1 回答 1

1

这是因为NA“ras”中的值。

取出时差别很小z[ind]<-NA。这是因为projectRaster隐式使用“na.rm=TRUE”；也许需要一个论据来改变这一点。这不像在其他情况下那么简单，例如，只有一个从 may 插值的单元格是NA，在这种情况下可能应该计算它。

在实践中，很少能找到这么多的 NA 值，通常它们仅限于（土地）边缘。

顺便说WGS84一句，它是一个基准，它不是一个投影。您使用的投影将是“经度/纬度”（也称为其他名称），但这也不是投影，因为投影的整个点都是从这样的角坐标到平面坐标。因此，您正在将一个坐标参考系（'long/lat，WGS84'）转换为另一个（BNG）。

于 2015-06-22T16:41:29.503 回答