Question

我有两个栅格：

raster1
    class      : RasterLayer 
band       : 6  (of  6  bands)
dimensions : 2803, 5303, 14864309  (nrow, ncol, ncell)
resolution : 0.008333333, 0.008333333  (x, y)
extent     : 60.85, 105.0417, 15.95833, 39.31667  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
source     : XXXX
names      : GPP
values     : 0, 435.5 (min, max)

raster2
class      : RasterLayer 
dimensions : 2803, 5303, 14864309  (nrow, ncol, ncell)
resolution : 0.008333333, 0.008333333  (x, y)
extent     : 60.85, 105.0417, 15.95833, 39.31667  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
source     : memory
names      : GPP
values     : 0, 506.4  (min, max)

我想根据土地覆盖情况在这两个栅格之间绘制散点图。我的土地覆盖物栅格文件是：

landcover
class      : RasterLayer 
dimensions : 2803, 5303, 14864309  (nrow, ncol, ncell)
resolution : 0.008333333, 0.008333333  (x, y)
extent     : 60.85, 105.0417, 15.95833, 39.31667  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
source     : XXXX 
names      : landuse
values     : 1, 12  (min, max)

对于变量GPP，我如何在raster1和raster2之间绘制散点图，对于可以说是常绿森林的土地覆盖，在“ landcover”的栅格文件中用值1表示。

Answer 1

第一步是根据您确定的土地用途对每个栅格进行分组。然后，您可以直接绘制它们：

plot(raster1[landcover==1],raster2[landcover==1])

如何根据土地覆盖率计算R中的散点图？

1 个答案: