Importieren, Bearbeiten, Bearbeiten und Exportieren von .tif-Dateien in R

Question

Ich arbeite mit einem .tif-Bild, das von einem Sequoia Parrot-Sensor aufgenommen wurde. Ich möchte eine radiometrische Kalibrierung durchführen und das resultierende Bild im selben Format (.tif) exportieren.

Ich importiere das Bild als Raster, dann mit einigen Algorithmen und schließlich versuchen, als eine .tif-Datei zu exportieren, aber es ist unmöglich zu öffnen. Die resultierende Datei ist 7 MB, aber das Bild kann nicht angezeigt werden.

Hier ist mein Skript:

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

Original-Skript ist gut gegliedert? Mit Ihrem Beispielbild funktioniert alles richtig.

Ich habe kleine Änderungen zum besseren Verständnis:

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)) 
}