R语言中RasterLayer赋值后数据丢失的技术问询

嘿，我来帮你搞定这个栅格数据的问题！你遇到的情况是因为用错了赋值方法，咱们一步步拆解清楚：

问题根源 #

你用ras2 = raster(ras1)的时候，这个raster()函数的作用是创建一个和输入栅格拥有相同空间属性（维度、投影、范围）但完全空的RasterLayer——它只会复制栅格的"骨架"信息，并不会拷贝原栅格的像素数据，所以summary(ras2)才会全是NA。

正确的赋值方式 #

根据你的需求，有两种常用的正确操作：

1. 直接引用赋值（内存共享） #

如果不需要修改新变量而不影响原栅格，直接用普通赋值即可，这种方式会让两个变量指向同一份数据，能节省内存：

ras2 <- ras1
summary(ras2)  # 结果会和ras1的统计完全一致

2. 创建独立副本（数据完全复制） #

如果你需要修改ras2但不想影响原ras1，用copy()函数创建一个完全独立的副本：

library(raster)
ras2 <- copy(ras1)
summary(ras2)  # 同样会输出正常的像素统计结果

补充：加载栅格到内存（大栅格慎用） #

你的原栅格是从磁盘读取的Sentinel数据，默认是延迟加载状态（只读取元数据，不把像素加载到内存）。如果想把数据完全加载到内存后再赋值，也可以用readAll()：

ras1 <- readAll(ras1)
ras2 <- ras1

不过这种方式对于Sentinel这类大栅格来说会占用大量内存，非必要不推荐使用。

验证效果 #

运行上面的正确代码后，你再查看ras2的属性和统计值，就会发现数据正常存在了：

ras2  # 会显示完整的栅格属性和数据源信息
summary(ras2)  # 会输出和ras1一致的像素统计结果

内容的提问来源于stack exchange，提问作者Xplorer