如何解决R 中的栅格面积计算问题
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
我有一些问题:
-
这些单位是平方米吗?
-
为什么当我尝试(手动)计算总面积(对于两个值)时,结果不同? (例如
31191 (ncell) * (8000(m) * 8000(m)) = 1996224000000
)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。