r语言 - 栅格重投影到拉特龙改变值

Question

我从以下栅格开始:

class      : RasterLayer 
dimensions : 900, 900, 810000  (nrow, ncol, ncell)
resolution : 1, 1  (x, y)
extent     : -523.4622, 376.5378, -4658.645, -3758.645  (xmin, xmax, ymin, ymax)
crs        : +proj=stere +lat_0=90 +lon_0=10 +k=0.93301270189 +x_0=0 +y_0=0 +ellps=WGS84 +units=km +no_defs 
source     : memory
names      : layer 
values     : 0, 2500  (min, max)

如果我用:

raster_projection <- projectRaster(rb, crs = "+proj=longlat +datum=WGS84 +no_defs")

它创建了一个重投影，在数据上有显著的变化:

class      : RasterLayer 
dimensions : 944, 1014, 957216  (nrow, ncol, ncell)
resolution : 0.0136, 0.00857  (x, y)
extent     : 2.00348, 15.79388, 46.99609, 55.08617  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
source     : memory
names      : layer 
values     : -0.6159579, 3583.708  (min, max)

我认为数据点是这样说的:"太高了"。

这是意料之中的吗?

如何创建一个更真实的投影?

谢谢!

Answer 1

新的极值可能在边缘;所以你的问题可能是由几个细胞引起的。更合适的比较方法是看箱线图或分位数。你也可以用terra::project代替。我认为这种方法不太容易出现这种情况。