我正在尝试汇总一张土地覆盖图,分辨率为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)
谁能告诉我如何匹配两个栅格的范围,希望可以解决此问题!