Lesen einer .tif-Datei in R

Question

Ich lese eine .tif-Datei in R und bekomme 4 Warnungen unten aufgelistet. Wenn ich der Anweisung in der 4. Meldung folge, bleiben die ersten 3 Warnungen erhalten, aber die aus der Datei gelesenen Werte ändern sich drastisch bei jedem Pixel. Bitte helfen Sie mir, die Daten korrekt aus den .tif-Dateien zu lesen. Beispieldatei finden Sie auf den Link zu finden: ftp://ftp.ntsg.umt.edu/pub/MODIS/NTSG_Products/MOD16/MOD16A2_MONTHLY.MERRA_GMAO_1kmALB/GEOTIFF_0.05degree/

mein Code:

remove(list=ls()) 

library(tiff) 

library(raster) 

str_name<-'MOD16A2_ET_0.05deg_GEO_2008M01.tif' 

read_file<-readTIFF(str_name) 

Warnmeldungen:

1: In readTIFF(str_name) : 
    TIFFReadDirectory: Unknown field with tag 33550 (0x830e) encountered 
2: In readTIFF(str_name) : 
    TIFFReadDirectory: Unknown field with tag 33922 (0x8482) encountered 
3: In readTIFF(str_name) : 
    TIFFReadDirectory: Unknown field with tag 34735 (0x87af) encountered 
4: In readTIFF(str_name) : 
    tiff package currently only supports unsigned integer or float sample formats in direct mode, but the image contains signed integer format - it will be treated as unsigned (use native=TRUE or convert=TRUE to avoid this issue) 

Bitte helfen Sie mir mit diesem Thema richtig tif-Dateien zu lesen. Danke im Voraus.

Answer 1

Haben Sie einfach die Raster-Paket-Rasterfunktion versucht (oder stapeln Sie wenn mehrere Ebenen tif)? Das Raster Paket wurde hergestellt mit georeferenzierten Raster-Dataset zu behandeln:

library(raster) 
str_name<-'MOD16A2_ET_0.05deg_GEO_2008M01.tif' 
imported_raster=raster(str_name) 

Der einfache Code über Werke und ergibt ein Rasterobjekt mit den folgenden Eigenschaften:

class  : RasterLayer 
dimensions : 2800, 7200, 20160000 (nrow, ncol, ncell) 
resolution : 0.05, 0.05 (x, y) 
extent  : -180, 180, -60, 80 (xmin, xmax, ymin, ymax) 
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : C:\Users\lfortini\Downloads\MOD16A2_ET_0.05deg_GEO_2000M01.tif 
names  : MOD16A2_ET_0.05deg_GEO_2000M01 
values  : -32768, 32767 (min, max) 

Answer 2

einfach den Pixel als unsigned lesen und wandeln sie in unterzeichnet:

t = readTIFF("MOD16A2_ET_0.05deg_GEO_2008M01.tif", as.is=TRUE) 
t[t >= 32738L] = -65536L + t[t >= 32738L] 

im Bild suchen, können Sie wollen auch -32.768-NA konvertieren, wie, dass die Verwendung in der Datei zu sein scheint:

t[t == -32768L] = NA 

Wenn Sie jetzt die ganzen Zahlen auf [-1,1] reellen Zahlen konvertieren wollen, tun nur

t = t/32768 

Die ersten drei Warnungen nur Sie sagen, dass es zusätzliche benutzerdefinierte Tags in der Datei .