我想将栅格中每个像素的降雨值(mm)乘以R下每个像素对应的面积。
首先,我直接将栅格乘以平均像素面积(总面积/像素数)。但是如果我们在全球范围内使用WGS84投影,实际上赤道和高纬度地区的像素面积是非常不同的?这就是我问这个问题的原因。
预先感谢您的帮助
发布于 2017-07-19 20:43:58
您可以使用以下命令计算每单位面积的每个像素的值(此处以米为单位):
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)))
https://stackoverflow.com/questions/45076550
复制