时间戳

首页 > 解决方案 > 在R中绘制带有大栅格的散点图？

r - 在R中绘制带有大栅格的散点图？

问题描述

我在

library(raster)
library(sp)
library(rgdal)
library(maptools)
library(sf)
library(dplyr)
library(devtools)
library(DGVMTools)
library(Metrics)
library(hydroGOF)
library(sp)
library(grid)
library(latticeExtra)

> Y
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      : layer 
values     : 0, 26.53035  (min, max)

> X
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      : VegH
values     : 0, 17.99169  (min, max)

我正在绘制这两个栅格之间 的plot(x,Y)

散点图 1

然而，由于像素太多，警告来了：

警告消息：在 .local(x, y, ...) 中：绘图使用了 0.7% 的单元格样本。您可以使用“maxpixels”来增加样本）

将两个栅格转换为并绘制后，散点图出现大约需要 1 小时，并且显示一个大的黑色块。 在此处输入图像描述 可重现的栅格：

r1 <- r2 <- raster(nrows=2803, ncols=5303)
values(r1) <- runif(ncell(r1))
values(r2) <- runif(ncell(r2))

我的问题是如何有效地绘制两个大型栅格数据集之间的散点图，这也可以进行视觉诊断？

标签： rrasterscatter-plottiff

解决方案

有这么多的细胞，如果你只是想检查散点图，样本应该是要走的路吗？您可以使用“maxpixels”来增加样本（如您收到的警告消息所示）。例如：

X <- Y <- raster(nrows=2803, ncols=5303)
values(X) <- 10 * runif(ncell(X))
values(Y) <- values(X) * (runif(ncell(Y)) - 0.5)

plot(X, Y, maxpixels=1e5)

但是，这确实需要一段时间。

也许这是一种更好的方法——这也需要一段时间，但可以更好地解释

plot(X, Y, maxpixels=1e5, gridded=TRUE)

推荐阅读