微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

R 中的栅格面积计算问题

如何解决R 中的栅格面积计算问题

我有一个包含以下信息的 RasterLayer

dimensions : 111,281,31191  (nrow,ncol,ncell)
resolution : 0.08333333,0.08333333  (x,y)
crs        : +proj=longlat +datum=wgs84 +no_defs +ellps=wgs84 +towgs84=0,0 

我的栅格的值是 1 和 6,我想要的是计算每个值所具有的以平方公里为单位的面积。

为此,我正在使用

tapply(area(my_raster),my_raster[],sum,na.rm=TRUE,weight=TRUE)

这给了我以下内容

1        6 
87299.97 19834.27

我有一些问题:

  1. 这些单位是平方米吗?

  2. 为什么当我尝试(手动)计算总面积(对于两个值)时,结果不同? (例如31191 (ncell) * (8000(m) * 8000(m)) = 1996224000000

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。