绘制一个大的光栅文件会产生一个白色的框架
Plotting a big raster file results in a white frame
我目前正在使用 256x256 像素的 ASCII 矩阵。我将它正确地导入到 R 中,对其进行栅格化,并且值是我所期望的(即,正确的 x 和 y 边界以及最小和最大“z”值)。然而,在绘制它时我得到一个空白光栅,就像矩阵中的每个值都是零一样。
我尝试创建另一个文件作为 5x5 矩阵,我没有遇到任何问题。我错过了什么吗?
文件及截图如下:
my 256x256 raster
https://gofile.io/d/JGApXI ascii 矩阵 link
您的栅格几乎是空的,因为它只有 2% 的值 !=0
。但是,如果您导出栅格并在 GIS 软件(如 Qgis 或 ArcMap)中将其可视化,通过为 0 值设置 100% 透明度,您可以看到剩余的值:
举个例子:
library(raster)
x <- read.table("D:/muon sideways0000.txt")
x <- as.matrix(x)
r <- raster(x)
writeRaster(r,"D:/r.tif")
z <- apply(x, 1, function(x)sum(x!=0))
sum(z)/ncell(r)*100
为了帮助可视化,您可以这样做
library(terra)
x <- read.table("muon sideways0000.txt")
x <- as.matrix(x)
r <- rast(x)
plot(r > 1)
或其他一些转换,例如
rr <- clamp(r, 0, 100)
plot(rr)
我目前正在使用 256x256 像素的 ASCII 矩阵。我将它正确地导入到 R 中,对其进行栅格化,并且值是我所期望的(即,正确的 x 和 y 边界以及最小和最大“z”值)。然而,在绘制它时我得到一个空白光栅,就像矩阵中的每个值都是零一样。
我尝试创建另一个文件作为 5x5 矩阵,我没有遇到任何问题。我错过了什么吗?
文件及截图如下: my 256x256 raster
https://gofile.io/d/JGApXI ascii 矩阵 link
您的栅格几乎是空的,因为它只有 2% 的值 !=0
。但是,如果您导出栅格并在 GIS 软件(如 Qgis 或 ArcMap)中将其可视化,通过为 0 值设置 100% 透明度,您可以看到剩余的值:
举个例子:
library(raster)
x <- read.table("D:/muon sideways0000.txt")
x <- as.matrix(x)
r <- raster(x)
writeRaster(r,"D:/r.tif")
z <- apply(x, 1, function(x)sum(x!=0))
sum(z)/ncell(r)*100
为了帮助可视化,您可以这样做
library(terra)
x <- read.table("muon sideways0000.txt")
x <- as.matrix(x)
r <- rast(x)
plot(r > 1)
或其他一些转换,例如
rr <- clamp(r, 0, 100)
plot(rr)