R语言使用另一个栅格扩展栅格的空间范围

whhtz7ly 于 2023-06-27 发布在其他

关注(0)|答案(1)|浏览(163)

我想知道是否有人知道我如何扩展（物理和维度）一个小光栅的空间范围，以匹配一个大光栅。我在terra中使用过extend函数，但它返回的只是相同的大小。我已经尝试了同样的事情在ArcGIS Pro与光栅计算器，但没有运气。我需要扩展的部分是小图像的最北部分。扩展部分将填充零。
数据在这里：https://drive.google.com/drive/folders/1c5AZS9onKUOzNoIVw3py3QFIRxITCjZ_?usp=sharing
以下是栅格信息

>large 

class       : SpatRaster 
dimensions  : 762, 762, 1  (nrow, ncol, nlyr)
resolution  : 3100, 3100  (x, y)
extent      : -3034878, -672677.7, 2022791, 4384991  (xmin, xmax, ymin, ymax)
coord. ref. : Canada_Lambert_Conformal_Conic 
source      : Yukon_poly_raster_project.tif 
name        : Yukon_poly_raster_project 
min value   :                         0 
max value   :                       127

>small 

class       : SpatRaster 
dimensions  : 426, 400, 1  (nrow, ncol, nlyr)
resolution  : 3096.673, 3098.812  (x, y)
extent      : -2284678, -1046008, 2551588, 3871682  (xmin, xmax, ymin, ymax)
coord. ref. : Canada_Lambert_Conformal_Conic 
source      : FRT_Canada_yukon_ras_1.tif 
name        : FRT_Canada_yukon_ras_1 
min value   :                      3 
max value   :                     11

以下是我的R代码：

library(terra)
r = extend(small, large)

r

来源：https://stackoverflow.com/questions/76518898/expand-spatial-exent-of-a-raster-using-another-raster

1条答案

按热度按时间

vulvrdjw1#

你的R代码有什么问题？我觉得extend的结果很好。

library(terra)
#terra 1.7.38

large <- rast(nrow=762, ncol=762, ext=c(-3034878, -672677.7, 2022791, 4384991))
small <- rast(nrow=426, ncol=400, ext=c(-2284678, -1046008, 2551588, 3871682))

e <- extend(small, large)
e
#class       : SpatRaster 
#dimensions  : 763, 763, 1  (nrow, ncol, nlyr)
#resolution  : 3096.675, 3098.812  (x, y)
#extent      : -3034073, -671310.3, 2021691, 4386085  (xmin, xmax, ymin, ymax)
#coord. ref. :  

plot(ext(e), lwd=3)
lines(ext(small), col="red", lwd=2)

但我认为你的目的是，在这种情况下：

z <- resample(small, large)

赞(0）回复(0）举报 2023-06-27

我来回答

R语言使用另一个栅格扩展栅格的空间范围

1条答案

相关问题

热门标签

最新问答

R语言 使用另一个栅格扩展栅格的空间范围

1条答案

相关问题

热门标签

最新问答

R语言使用另一个栅格扩展栅格的空间范围