R Wie multipliziert man die Werte der Rasterpixel mit den wahren Flächen dieser Pixel?

Question

Ich möchte die Niederschlagswerte (mm) jedes Pixels eines Rasters mit der Fläche multiplizieren, die jedem dieser Pixel unter R.

entspricht Zuerst multiplizierte ich das Raster direkt mit der durchschnittlichen Pixelfläche (Gesamtfläche/Anzahl der Pixel). Aber wenn wir mit der WGS84-Projektion auf globaler Ebene arbeiten, ist die Pixelfläche am Äquator und in den Regionen mit hoher geografischer Breite sehr unterschiedlich? Deshalb stelle ich diese Frage.

Vielen Dank im Voraus für Ihre Hilfe

Answer 1

Sie können den Wert jedes Pixels pro Flächeneinheit berechnet werden (hier in Metern) mit:

rows = init(raster , v='y') # latitude bands 
PI = 3.14159265359 # PI number 
earth_radius = 6371007.181 # radius of earth in meter 
raster_resolution = xres(raster) # resolution of the raster layer 

new_raster <- raster * ((raster_resolution * (PI/180) * earth_radius) * (raster_resolution * (3.14159265359/180) * earth_radius * cos(rows/180*PI)))