Pythonでポリゴンの領域を見つけるプログラム

順序付けられたポイントのリストが2D平面上の単純なポリゴンエンドポイントを表すとします。このポリゴンの領域を見つける必要があります。

したがって、入力がpoints =[（0、0）、（0,5）、（3、5）、（3,0）]のような場合、出力は15になります。

これを解決するには、次の手順に従います-

• 関数getInfo（）を定義します。これにはx1、y1、x2、y2が必要です
• return x1 * y2-y1 * x2
• メインの方法から、次の手順を実行します
• N：=ポイントのサイズ
• （firstx、firsty）：=points [0]
• （prevx、prevy）：=（firstx、firsty）
• res：=0
• 1からN-1の範囲のiについては、
  • （nextx、nexty）：=points [i]
  • res：=res + getInfo（prevx、prevy、nextx、nexty）
  • prevx：=nextx
  • prevy：=nexty
• res：=res + getInfo（prevx、prevy、firstx、firsty）
• return | res | / 2.0

例

理解を深めるために、次の実装を見てみましょう-

def getInfo(x1, y1, x2, y2):
   return x1*y2 - y1*x2

def solve(points):
   N = len(points)
   firstx, firsty = points[0]
   prevx, prevy = firstx, firsty
   res = 0

   for i in range(1, N):
      nextx, nexty = points[i]
      res = res + getInfo(prevx,prevy,nextx,nexty)
      prevx = nextx
      prevy = nexty
   res = res + getInfo(prevx,prevy,firstx,firsty)
   return abs(res)/2.0

points = [(0, 0), (0,5), (3, 5), (3,0)]
print(solve(points))

入力

[(0, 0), (0,5), (3, 5), (3,0)]

出力

15.0