0

我在一个文件夹中有大约 40 个 .tiff 格式的空间栅格。我正在尝试从 R 中的每个栅格生成直方图，并将每个直方图作为 jpeg 保存在单独的文件夹中。我编写了代码来遍历每个光栅，创建一个直方图并使用“jpeg”包保存它。

setwd("G:/Research/MODIS Albedo Oct 08-July 09/Test") 
library(raster)
library(jpeg)
files <- list.files(path="G:/Research/MODIS Albedo Oct 08-July 09/Test", pattern=".tif",all.files=T, full.names=F, no.. = T) #generate a list of rasters in the folder
number<-length(files) #count the number of rasters
for(r in 1:number) #loop over each raster in the folder
{
  x<-raster(files[r], header=F) #load one raster file
  jpeg("G:/Research/MODIS Albedo Oct 08-July 09/test_histplots/r.jpg") #create jpeg using the name 'r' generated by loop 
  hist(x) #generate histogram
  dev.off()  
}

我希望每个生成的 jpeg 都具有不同的名称，最好是原始栅格名称的子集。例如，如果栅格的原始名称为“MODIS101_265”，则 jpeg 的名称应为 265。这里，265 是年份中的儒略日期。我假设这可能涉及在 C 中使用像 %d 这样的格式说明符，但我不确定这在 R 中是如何工作的。

当我运行上面的代码时，我只得到一个直方图。似乎代码在原始栅格上正确循环，但将所有生成的直方图保存到单个 jpeg。

任何建议都会有所帮助！谢谢！

4

1 回答 1

1

使用正则表达式gsub是您从名称中获取数字的朋友。假设您的所有文件都以示例的方式命名（“MODIS101_265.tif”），那么下面的代码将为您工作。

另外，欢迎使用 R ，其中for循环很慢，通常可以用更快的lapply.

saveMyHist <- function(fileName) {
  fileNum <- as.numeric(gsub(".*_(\\d+)\\.tif", "\\1", fileName))
  x <- raster(fileName, header=F)
  jpeg(sprintf("G:/Research/MODIS Albedo Oct 08-July 09/test_histplots/%03d.jpg", 
    fileNum))
  hist(x)
  dev.off()
}

files <- list.files("/Users/home/Documents/Development/Rtesting", 
  pattern=".tif",all.files=T, full.names=F, no.. = T)
lapply(files, saveMyHist)
于 2014-07-09T04:48:32.470 回答