前回、３点からなる三角形の面積を外積を用いて求めました。
これを多角形へ応用したいと思います。
まず、外積のおさらいから。
Ｚ成分が０（ゼロ）の２つのベクトル
の外積は
となり、Ｚ成分の大きさが２つのベクトルのなす平行四辺形の面積となり、三角形の面積はこの半分（1/2）となります。
さらに、ベクトルａ から ベクトルｂ への向きが反時計方向の場合、
ベクトルａ と ベクトルｂ の外積のＺ成分の値は正となり、
逆に時計方向の場合、Ｚ成分は負となります。
これを踏まえて、３点からなる三角形の面積を求めるの時は三角形の辺上にベクトルを取りましたが、今回は原点と多角形の頂点の座標とで成すベクトルとします。
ここで、多角形の頂点の座標をＰ１～Ｐ３のように反時計方向に定義します。
ただし、Ｚ座標は０（ゼロ）とします。
ベクトルＰ０→Ｐ１ と ベクトルＰ０→Ｐ２ の外積のＺ成分の値は時計方向なので、負となります。
同様に、ベクトルＰ０→Ｐ２ と ベクトルＰ０→Ｐ３ の外積のＺ成分の値は反時計方向なので、正となります。
ベクトルＰ０→Ｐ３ と ベクトルＰ０→Ｐ１ の外積のＺ成分の値も反時計方向なので、正となります。
これらの外積の結果のＺ成分を足して1/2にすると、求めたい三角形Ｐ１Ｐ２Ｐ３の面積が求まります。
この事をｎ点からなる多角形へ応用すると、下図のような図形の場合、
点Ｐ１～Ｐ３までは時計方向となるので、外積のＺ成分は負となります。
点Ｐ３～Ｐ７、Ｐ１までは反時計方向となるので、外積のＺ成分は正となります。
これら全ての外積のＺ成分を足し、1/2にすると多角形の面積が求まります。
この事を一般式で書くと、頂点の座標をPi (xi, yi) とすると
となります。
i は i = 1, 2, 3・・・nのインデックス番号、
| | は絶対値です。
ただし、i = n のとき、n+1 = 1 とします。
また、絶対値を取っているのは、頂点の座標が時計方向へ割り振られた場合にも対応できるようにしています。