使用 R 只读一段 tif答案

【问题标题】：Read only a section of tif using R使用 R 只读一段 tif
【发布时间】：2018-03-29 08:34:25
【问题描述】：

我有一个不适合内存的大 tif 文件。我想知道是否可以根据范围仅读取此 tif 文件的一部分。

我正在使用 R，到目前为止还没有找到解决方案。我知道它可以使用 C# 来完成，所以我想 R 中应该有办法。

谢谢！

【问题讨论】：

降低分辨率试试

标签： r raster tiff

【解决方案1】：

查看rtiff 包。 readTiff 函数可以减小图片大小。

【讨论】：

【解决方案2】：

如果您在 R 中使用栅格，您可能会使用 raster 包。

你可以这样做：

library(raster)

# testraster

r <- raster(system.file("external/test.grd", package="raster"))

# not in memory

inMemory(r)
# [1] FALSE

如果我们查看栅格对象，我们可以看到栅格有 115 行和 80 列：

r
# class       : RasterLayer 
# dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
# resolution  : 40, 40  (x, y)
# extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
# coord. ref. : +init=epsg:28992 +towgs84=565.237,50.0087,465.658,-0.406857,0.350733,-1.87035,4.0812 +proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs 
# data source : /usr/local/lib/R/site-library/raster/external/test.grd 
# names       : test 
# values      : 128.434, 1805.78  (min, max)

所以让我们从中心加载一个 10 x 10 的块：

arr <- getValuesBlock(r,row=58,nrows=10,col=40, ncols=10,format='matrix')

这将返回一个矩阵，您可以对其进行进一步处理、绘图等。

plot(raster(arr))

【讨论】：