在R中導入，處理，編輯和導出.tif文件

Question

我正在處理從紅杉鸚鵡傳感器捕獲的.tif圖像。我想要進行輻射校準並以相同格式（.tif）導出生成的圖像。

我將圖像作爲柵格導入，然後使用一些算法進行處理，最後嘗試導出爲.tif文件，但無法打開。生成的文件爲7 MB，但無法查看圖像。

這裏是我的腳本：

setwd("/where the images are/") 
rlist=list.files(getwd(), pattern="TIF$", full.names=F) 
options(digits=20) 

for(i in rlist){ 
    data <- raster(i) 

meta <- exifr(i, recursive = FALSE, quiet = TRUE, exiftoolargs = NULL) 
SM <- meta$SensorModel 
SM <- strsplit(SM, ",")[[1]] 
A <- as.numeric(SM[1]) 
B <- as.numeric(SM[2]) 
C <- as.numeric(gsub("[^0-9\\.]", "", SM[3])) 

Ep <- meta$ExposureTime ## Epsilon 
f <- meta$FNumber ## Focus Number 
ys <- meta$ISO ##ISO 

I <- f^2*(data-B)/(A*Ep*ys+C) 
I <- flip(I,"x") 
I <- flip(I,"y") 

Answer 1

是原來的劇本以及縮進？用您的示例圖像，一切正常。

我爲更好地瞭解小的變化：在``先計算I`

library(raster) 
library(exifr) 

setwd("/where the images are/") 
rlist=list.files(getwd(), pattern="TIF$", full.names=F) 
options(digits=20) 

for(i in seq_along(rlist)){ # small change 

    data <- raster(rlist[i]) # small change 

    meta <- exifr(rlist[i], recursive = FALSE, quiet = TRUE, exiftoolargs = NULL) # small change 
    SM <- meta$SensorModel 
    SM <- strsplit(SM, ",")[[1]] 
    A <- as.numeric(SM[1]) 
    B <- as.numeric(SM[2]) 
    C <- as.numeric(gsub("[^0-9\\.]", "", SM[3])) 

    Ep <- meta$ExposureTime ## Epsilon 
    f <- meta$FNumber ## Focus Number 
    ys <- meta$ISO ##ISO 

    I <- calc(data,fun = function(x){f^2*(x-B)/(A*Ep*ys+C)}) # small change 
    I <- flip(I,"x") 
    I <- flip(I,"y") 
    print(plot(I)) 
}