如何匹配R中两个栅格的范围?

时间:2019-08-07 11:32:55

标签: r ggplot2 raster aggregation extent

我正在尝试汇总一张土地覆盖图,分辨率为10kmx10km,而不是25m。我有一个英国的10公里见方的栅格网格,用于汇总土地覆盖图。

当我尝试在“ grainchanger” R程序包中执行此操作时,出现以下错误:

Error in .local(x, y, ...) : extents do not overlap

当我在栅格上运行命令show()时,这就是我得到的: (fine_dat是土地覆盖图,g_raster是10公里的网格)

> show(fine_dat)
class      : RasterLayer 
band       : 1  (of  21  bands)
dimensions : 1300, 700, 910000  (nrow, ncol, ncell)
resolution : 1000, 1000  (x, y)
extent     : 0, 7e+05, 0, 1300000  (xmin, xmax, ymin, ymax)
crs        : +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 
 +y_0=-100000 +ellps=airy 
 +towgs84=446.448,-125.157,542.06,0.15,0.247,0.842,-20.489 +units=m 
+no_defs 
source     : /Users/bexjenks/Documents/MASTERS/RESEARCH 
PROJECT/agri_intensity_wp/LCM2015_GB_1km_percent_cover_target_class.tif 
names      : LCM2015_GB_1km_percent_cover_target_class 
values     : 1, 100  (min, max)



> show(g_raster)
class      : RasterLayer 
band       : 1  (of  3  bands)
dimensions : 295, 162, 47790  (nrow, ncol, ncell)
resolution : 1, 1  (x, y)
extent     : 0, 162, 0, 295  (xmin, xmax, ymin, ymax)
crs        : NA 
source     : /Users/bexjenks/Documents/MASTERS/RESEARCH 
PROJECT/agri_intensity_wp/raster.tif 
names      : raster 
values     : 0, 255  (min, max)

谁能告诉我如何匹配两个栅格的范围,希望可以解决此问题!

0 个答案:

没有答案