Question

我有两个栅格，我想要一个到另一个的空间范围。然后将其另存为新栅格。我使用以下代码。但是，我无法将具有新空间范围的2013图像另存为新栅格。非常感谢任何指导。

raster_2013 <- raster("avgt2013.tif")
extent(raster_2013)
class       : Extent 
xmin        : 112.91 
xmax        : 153.64 
ymin        : -43.75 
ymax        : -9 
> res(raster_2013)
[1] 0.01 0.01
> 
> raster_2015 <- raster("avgt2015.tif")
> extent(raster_2015)
class       : Extent 
xmin        : 112 
xmax        : 154 
ymin        : -44 
ymax        : -9 
> res(raster_2015)
[1] 0.01 0.01
> 
> e <- extent(112, 154, -44, -9)
> 
> ex = extent(raster_2015)
> r2 = crop(raster_2013, ex)
> 
> 
> new_2013 <- alignExtent(e, raster_2013, snap='near')
> str(new_2013)
Formal class 'Extent' [package "raster"] with 4 slots
  ..@ xmin: num 112
  ..@ xmax: num 154
  ..@ ymin: num -44
  ..@ ymax: num -9
> 
> rc <- crop(raster_2013, e, snap='near')
> extent(rc)
class       : Extent 
xmin        : 112.91 
xmax        : 153.64 
ymin        : -43.75 
ymax        : -9

Answer 1

首先，请make a simple reproducible example to ask a question。

library(raster)

set.seed(11)
raster_2013 = raster(ext=extent(112.91, 153.64, -43.75, -9), res=c(0.01, 0.01))
raster_2013[] = rnorm(ncell(raster_2013))
raster_2015 = raster(ext=extent(112, 154, -44, -9), res=c(0.01, 0.01))
raster_2015[] = rnorm(ncell(raster_2015))

然后，您的代码存在一些问题。在您的情况下，alignExtent是无用的，因为两个栅格具有相同的分辨率，并且它们的范围与此分辨率相对应。

如果你的目标是将raster_2015的范围赋予raster_2013，你需要意识到extent(raster_2015)相对于xmin更短（更小），但在其他地方更大或相等。因此crop ping只会影响raster_2013的xmin。您首先需要extend，然后才需要crop才能获得完全相同的范围：

new_2013 <- crop(extend(raster_2013, raster_2015), raster_2015)
all.equal(extent(raster_2015), extent(new_2013))
#[1] TRUE

正如@ Geo-sp所提到的，你也可以resample raster_2013，但如果rwo栅格不对齐，你通常会使用它（并且要知道它会在case，由于插值导致修改数据）。在这里，由于它们是，它将给出与crop(extend())相同的结果，但它会更慢，更耗费资源：

system.time(new_2013 <- crop(extend(raster_2013, raster_2015), raster_2015))
#  user  system elapsed 
# 0.676   0.036   0.712 
system.time(new_2013_res <- resample(raster_2013, raster_2015))
#   user  system elapsed 
# 10.324   0.536  10.869
all.equal(new_2013, new_2013_res)
# [1] TRUE

最后，为了保存它，你可以使用writeRaster，因为阅读文档会引导你; - ）

writeRaster(new_2013, "raster_2013_extent2015.grd")

在R中更改栅格空间范围

1 个答案: