如何在R中读取超出内存的大TIFF文件指定区域？

当然可以！你完全不用把整个超大TIFF文件塞进内存，R里有好几种方法能实现仅读取指定范围的内容，和C#的思路类似——都是基于文件连接来按需读取数据。下面给你两个实用的方案：

方案1：用terra包（推荐，高效简洁） #

terra是现在处理栅格数据的主流包，替代了旧的raster包，对大型文件的支持更友好，操作也更直观：

# 先安装并加载包
install.packages("terra")
library(terra)

# 打开TIFF文件——这一步只是建立文件连接，不会把全量数据载入内存
large_tiff <- rast("你的大型TIFF文件路径.tif")

# 方法A：通过地理/投影坐标范围读取
# 比如定义x从1000到5000，y从2000到6000（根据你的文件坐标系调整）
target_extent <- ext(1000, 5000, 2000, 6000)
subset_data <- crop(large_tiff, target_extent)

# 方法B：直接通过行列索引读取（如果你知道具体的行列范围）
# 比如读取第100到500列，第200到600行的区域
subset_data <- large_tiff[100:500, 200:600]

# 查看读取的结果
plot(subset_data)

方案2：用rgdal包（兼容旧项目） #

如果你需要兼容一些旧代码，rgdal包也能实现这个需求，它直接调用GDAL库，和C#里用GDAL的逻辑很像：

install.packages("rgdal")
library(rgdal)

# 打开TIFF文件的GDAL连接
tiff_connection <- GDAL.open("你的大型TIFF文件路径.tif")

# 定义要读取的窗口参数：
# xoff = 起始列, yoff = 起始行, xsize = 要读取的列数, ysize = 要读取的行数
read_window <- c(xoff = 100, yoff = 200, xsize = 401, ysize = 401)
# 注意：xsize是结束列 - 起始列 + 1，比如从100到500就是401列

# 读取指定窗口的数据
subset_matrix <- GDAL.read(tiff_connection, 
                           offset = c(read_window["xoff"], read_window["yoff"]),
                           region.dim = c(read_window["xsize"], read_window["ysize"]))

# 记得关闭文件连接
GDAL.close(tiff_connection)

# 可选：把读取到的矩阵转成栅格对象方便后续分析
library(raster)
subset_raster <- raster(subset_matrix, 
                        xmn = 1000, xmx = 5000, # 对应你读取区域的x范围
                        ymn = 2000, ymx = 6000) # 对应y范围

注意事项 #

• 不管用哪个包，关键都是先建立文件连接，再按需读取指定区域，避免全量加载。
• 如果你的TIFF是带地理坐标系的，要确保定义的范围和文件的坐标系一致，不然会读取到错误的区域。
• 对于多波段的TIFF，这两个方法也都支持指定读取某几个波段，比如large_tiff[[1:3]]（terra）或者在GDAL.read里指定band参数。

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