R如何將光柵像素的值與這些像素的真實表面相乘？

Question

我想通過對應於這些像素的每R.

下起初的面積乘以一個光柵的每個像素的值沉澱（mm）時，我直接乘以平均像素區域中的光柵（總面積/像素數量）。 但是，如果我們用WGS84投影在全球範圍內進行工作，實際上赤道和高緯度地區的像素區域是非常不同的？ 這就是爲什麼我問這個問題。

預先感謝您使用您的幫助

Answer 1

可以計算出每單位面積（這裏以米爲單位）的每個像素的值：

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