问题描述
我在 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 (将#修改为@）