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

【问题标题】：How to compute scatterplots in R according to land cover?如何根据土地覆盖计算R中的散点图？
【发布时间】：2020-03-08 17:44:56
【问题描述】：

我有两个栅格：

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 之间的散点图，例如常绿森林的土地覆盖，它由“土地覆盖”的光栅文件中的值 1 表示。

【问题讨论】：

标签： r raster

【解决方案1】：

第一步是根据您确定的土地利用对每个栅格进行子集化。然后你可以直接绘制它们：

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

【讨论】：

这段代码似乎为不同的土地覆盖类别给出了相同的散点图。您能否在“按您确定的土地用途对每个栅格进行子集”中详细说明？ raster1[landcover==1] 意思是你在上面的代码中写的吗？
这就是您所要求的“我如何为变量 GPP 绘制 raster1 和 raster2 之间的散点图，例如常绿森林的土地覆盖，它由光栅文件中的值 1 表示的“土地覆盖”。”