我想将栅格每个像素的降水值（mm）乘以R下每个像素对应的区域。
首先，我直接将光栅乘以平均像素区域（总面积/像素数）。 但是，如果我们使用WGS84投影在全球范围内工作，实际上赤道和高纬度地区的像素区域是非常不同的？ 这就是我提出这个问题的原因。
提前感谢您的帮助
答案 0 :(得分：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)))