R_Calculating具有不规则边界的区域的Thiessen权重

Question

我想计算泰森重量来计算点测量数量的面积降雨量。我正在使用R并且由于同一主题中的一些先前问题，我知道我可以使用 deldir。但问题是我的边界多边形不是矩形;它是一个不规则的多边形（它是使用ArcGIS派生的集水区边界）。但在deldir中，边界只能是一个矩形。是否还有其他包可以计算出由不规则边界覆盖的区域的泰森权重？

以下是我的测量点（meas_points）和（简化的）边界多边形（边界）的坐标

> meas_points
          X      Y
[1,] 415720 432795
[2,] 415513 432834
[3,] 415325 432740
[4,] 415356 432847
[5,] 415374 432858
[6,] 415426 432774
[7,] 415395 432811
[8,] 415626 432762

> boundary
          x      y
[1,] 415491 432947
[2,] 415269 432919
[3,] 415211 432776
[4,] 415247 432657
[5,] 415533 432657
[6,] 415781 432677
[7,] 415795 432836
[8,] 415746 432937

非常感谢任何帮助。感谢。

Answer 1

您可以尝试voronoi的双重性并找到凸包。然后删除所有内部三角形和超出alpha的边。重新插入内部三角形，找到Delaunay三角剖分的对偶。