计算给定方向列表的区域

时间:2015-01-20 05:24:20

标签: algorithm

让我们说出你的方向列表:

up, up, right, down, right, down, left, left

如果您按照说明操作,您将始终返回起始位置。计算刚刚创建的形状的区域。

上述方向形成的形状如下所示:

 ___
|   |___
|_______|

显然,从图片中可以看出该区域为3。

我尝试使用二维矩阵来追踪方向,但不确定如何从中获取该区域......

例如,在我的2d数组中:

O  O
O  O  O
O  O  O

这可能不是处理这个问题的好方法,任何想法?

4 个答案:

答案 0 :(得分:2)

由于您创建的多边形仅具有轴对齐边,因此您可以计算垂直板的总面积。

假设我们给出了一个顶点列表V。我假设我们已经包含在此列表中,因此我们可以查询V.next(v)每个顶点v in V。对于最后一个,结果是第一个。

首先,尝试找到最左边和最右边的点,以及到达最左边点的顶点(线性时间)。

x = 0                       // current x-position
xMin = inf, xMax = -inf     // leftmost and rightmost point
leftVertex = null           // leftmost vertex
foreach v in V
    x = x + (v is left ? -1 : v is right ? 1 : 0)
    xMax = max(x, xMax)
    if x < xMin
        xMin = x
        leftVertex = V.next(v)

现在我们创建一个简单的数据结构:对于每个垂直平板,我们保留一个最大堆(排序列表也很好,但我们只需要在最后重复获取最大元素)。

width = xMax - xMin
heaps = new MaxHeap[width]

我们现在开始从顶点leftVertex跟踪形状(我们在第一步中找到的最左边的顶点)。我们现在选择这个顶点有x / y位置(0,0),只是因为它很方便。

x = 0, y = 0
v = leftVertex
do
    if v is left
        x = x-1         // use left endpoint for index
        heaps[x].Add(y) // first dec, then store
    if v is right
        heaps[x].Add(y) // use left endpoint for index
        x = x+1         // first store, then inc
    if v is up
        y = y+1
    if v is down
        y = y-1

    v = V.next(v)
until v = leftVertex

您可以在O(n log n)时间内构建此结构,因为添加到堆会花费对数时间。

最后,我们需要从堆中计算区域。对于格式良好的输入,我们需要从堆中获取两个连续的y值并减去它们。

area = 0
foreach heap in heaps
    while heap not empty
        area = heap.PopMax() - heap.PopMax()

return area

同样,这需要O(n log n)时间。


我将算法移植到java实现中(参见Ideone)。两个样本运行:

public static void main (String[] args) {
    //  _
    // | |_
    // |_ _ |
    Direction[] input = { Direction.Up, Direction.Up, 
                          Direction.Right, Direction.Down,
                          Direction.Right, Direction.Down,
                          Direction.Left, Direction.Left };

    System.out.println(computeArea(input));

    //  _
    // |_|_
    //   |_|
    Direction[] input2 = { Direction.Up, Direction.Right, 
                           Direction.Down, Direction.Down,
                           Direction.Right, Direction.Up,
                           Direction.Left, Direction.Left };

    System.out.println(computeArea(input2));
}

返回(按预期):

3
2

答案 1 :(得分:1)

假设某个起点(例如,(0,0))且y方向为正向:

  • left将(-1,0)添加到最后一点。
  • 右(+1,0)加到最后一点。
  • up将(0,+ 1)添加到最后一点。
  • down将(0,-1)添加到最后一点。

然后,一系列方向将产生一个(x,y)顶点坐标列表,从中可以找到How do I calculate the surface area of a 2d polygon?

中所得(隐含的闭合)多边形的面积。

修改

这是Python中的实现和测试。前两个函数来自上面链接的答案:

def segments(p):
    return zip(p, p[1:] + [p[0]])

def area(p):
    return 0.5 * abs(sum(x0*y1 - x1*y0
                         for ((x0, y0), (x1, y1)) in segments(p)))

def mkvertices(pth):
    vert = [(0,0)]
    for (dx,dy) in pth:
        vert.append((vert[-1][0]+dx,vert[-1][1]+dy))
    return vert

left = (-1,0)
right = (+1,0)
up = (0,+1)
down = (0,-1)

#  _
# | |_
# |__|
print (area(mkvertices([up, up, right, down, right, down, left, left])))
#  _
# |_|_
#   |_|
print (area(mkvertices([up, right, down, down, right, up, left, left])))

输出:

3.0
0.0

请注意,对于包含相交线的多边形,此方法会失败,如第二个示例所示。

答案 2 :(得分:1)

对于简单的多边形,可以使用Shoelace公式就地实现。

对于每个段(a, b),我们必须计算(b.x - a.x)*(a.y + b.y)/2。所有段上的总和是多边形的有符号区域。

此外,这里我们只处理长度为1的与轴对齐的线段。由于b.x - a.x = 0,垂直线段可以忽略。 水平段具有a.y + b.y / 2 = a.y = b.yb.x - a.x = +-1。 因此,最后我们只需要跟踪y,添加的区域始终为+-y

这是示例C ++代码:

#include <iostream>
#include <vector>

enum struct Direction
{
    Up, Down, Left, Right
};

int area(const std::vector<Direction>& moves)
{
    int area = 0;
    int y = 0;

    for (auto move : moves)
    {
        switch(move)
        {
            case Direction::Left:
                area += y;
                break;
            case Direction::Right:
                area -= y;
                break;
            case Direction::Up:
                y -= 1;
                break;
            case Direction::Down:
                y += 1;
                break;
        }
    }

    return area < 0 ? -area : area;
}

int main()
{
    std::vector<Direction> moves{{
        Direction::Up, 
        Direction::Up, 
        Direction::Right, 
        Direction::Down, 
        Direction::Right,
        Direction::Down, 
        Direction::Left, 
        Direction::Left
        }};

    std::cout << area(moves);

    return 0;
}

答案 3 :(得分:0)

我假设您正在绘制的形状(轴对齐,多边形图,闭合,非相交线)应该有一些限制,以便能够计算面积。

使用线段表示形状,每个线段由两个点组成,每个点有两个坐标:x和y。

考虑到这些假设,我们可以说任何水平线段都有一个平行线段,其两个点具有相同的x维度,但y维度不同。

这两个线段之间的表面积等于它们之间的高度差。对所有水平线段的区域进行注释,可以得到形状的总表面积。