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

为什么 Eckert IV 投影在 ArcGIS 和 R 中的处理方式不同?

如何解决为什么 Eckert IV 投影在 ArcGIS 和 R 中的处理方式不同?

我在 ArcGIS Pro(2.6 版)和 R(4.0.3 版)中将一些多边形数据投影到 EckertIV 时遇到了不同的问题。我正在尝试计算此多边形内的栅格单元数(新西兰是专属经济区)。数据可在此处获得:https://drive.google.com/drive/folders/1OE_q_W97KHw7iC3l0v9y5EMrYMSRLW_u?usp=sharing。在 R 中,我加载数据,使用 proj4 符号投影到 Eckert IV,并使用 raster::extract() 计算多边形内像素的总和。

代码
图书馆(光栅)
图书馆(sf)
poly poly 单元格 ex

使用 R,我得到 557 个单元格的值。

在 ArcGIS 中,我将数据添加到项目中,然后将多边形数据投影到 EckertIV 中,并使用分区统计作为表格。这导致 841 个单元格。如果我在 ArcGIS 中使用 R 中制作的投影多边形并使用分区统计,我得到的值为 557,因此表明投影是问题所在。

我可以在 ArcGIS 中直观地看到,与 ArcGIS 中投影的多边形相比,R 中投影的多边形部分缺失。谁能解释为什么会有这种差异以及我可以做些什么来确保我在 R 中获得准确的结果? The portion of the polygon that is being missed is shown in pink

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