空間上の二点を結ぶ直線上に任意の点が存在するかどうかの関数

Cの初心者です。

空間上に存在する2点間を結ぶ直線上に任意の点が存在するかどうかの
関数を作りたいのですがどのような公式を用いて評価すればいいのか分かりません。
どなたかご教示ください。

引数
　　　始点　A(x,y,z)
　　　終点　B(x,y,z)
　　　直線上に存在するであろう任意の点
　　　　　　C(x,y,z)

　関数のイメージ
　　　boolean isOnLine (Ax, Ay, Az, Bx, By, Bz, Cx, Cy, Cz); 配列でもよいです。

返り値
　　　True、　False　　（ 0 or 1 ）

よろしくお願いします。

No.2 ベストアンサー 回答者： chie65535 回答日時： 2009/08/25 13:24

直線ＡＢの式を求め、その式に任意の点Ｃの座標を代入し、式が成り立つか判定すればＯＫ。

(Cx,Cy,Cz)=(Ax,Ay,Az)+t(Bx-Ax,By-Ay,Bz-Az)
この式が、点Ｃを通る直線ＡＢの式。

この式をt=の形に直せば
t=(Cx-Ax)/(Bx-Ax)=(Cy-Ay)/(By-Ay)=(Cz-Az)/(Bz-Az)
となる。

２つの点｛(Cx-Ax)、(Cy-Ay)、(Cz-Az)｝、｛(Bx-Ax)、(By-Ay)、(Bz-Az)｝は「点Ａが原点になるように、直線ＡＢと点Ｃを並行移動した時の、点Ｂと点Ｃの座標」を表わしている。

そして、原点と点Ｂの延長線上に点Ｃがあるならば各x、y、zの比は等しくなるので
(Cx-Ax):(Cy-Ay):(Cz-Az)=(Bx-Ax):(By-Ay):(Bz-Az)
が成り立つ。

つまり
(Cx-Ax)/(Bx-Ax)=(Cy-Ay)/(By-Ay)=(Cz-Az)/(Bz-Az)
が成り立つ。

プログラムでは
(Bx-Ax)が０、かつ、(By-Ay)が０、かつ、(Bz-Az)が０
(Bx-Ax)が０、かつ、(By-Ay)が０
(Bx-Ax)が０、かつ、(Bz-Az)が０
(By-Ay)が０、かつ、(Bz-Az)が０
(Bx-Ax)が０
(By-Ay)が０
(Bz-Az)が０
と言う７つの例外条件を取り除いた上
(Cx-Ax)/(Bx-Ax)と(Cy-Ay)/(By-Ay)と(Cz-Az)/(Bz-Az)が等しいかどうかを判定すれば良い。

７つの例外条件では、それぞれ

(Bx-Ax)が０、かつ、(By-Ay)が０、かつ、(Bz-Az)が０
の時
(Cx-Ax)が０、かつ、(Cy-Ay)が０、かつ、(Cz-Az)が０
ならば真（３つの点がすべて同じ座標にある）

(Bx-Ax)が０、かつ、(By-Ay)が０
の時
(Cx-Ax)が０、かつ、(Cy-Ay)が０
ならば真（直線ＡＢがＺ軸と並行）

(Bx-Ax)が０、かつ、(Bz-Az)が０
の時
(Cx-Ax)が０、かつ、(Cz-Az)が０
ならば真（直線ＡＢがＹ軸と並行）

(By-Ay)が０、かつ、(Bz-Az)が０
の時
(Cy-Ay)が０、かつ、(Cz-Az)が０
ならば真（直線ＡＢがＸ軸と並行）

(Bx-Ax)が０
の時
(Cx-Ax)が０、かつ、(Cy-Ay)/(By-Ay)と(Cz-Az)/(Bz-Az)が等しいならば真（直線ＡＢが「平面(Bx-Ax)=0」上にある）

(By-Ay)が０
の時
(Cy-Ay)が０、かつ、(Cx-Ax)/(Bx-Ax)と(Cz-Az)/(Bz-Az)が等しいならば真（直線ＡＢが「平面(By-Ay)=0」上にある）

(Bz-Az)が０
の時
(Cz-Az)が０、かつ、(Cx-Ax)/(Bx-Ax)と(Cy-Ay)/(By-Ay)が等しいならば真（直線ＡＢが「平面(Bz-Az)=0」上にある）

と判定すればよい。

なお、この判定では「点Ｃが直線ＡＢの延長線上にあり、点Ａと点Ｂの中間にない場合（線上にあるが外側にある場合）」にも真となる。

点Ｃが点Ａと点Ｂの中間にあるかどうかは「Ax、Bx、Cxの大小関係」で判定可能なので、上記の判定の前に判定してしまうと良い（Ay、By、Cy、Az、Bz、Czの大小関係は無視して良い）

No.8 回答者： Tacosan 回答日時： 2009/08/25 15:57

えぇと....

ベクトルの取り方を確認して下さい＞#7. #5 では「AB と AC」なのでなす角は 0 です. もちろんこれだと「C が B より A に近い」という条件を組み込まなければなりません. 一方, #3 のように「CA と CB」でとればなす角が π というだけですみます.
いずれにしても「角度を求める」という方針だと「C と A または B が一致する」可能性を念頭に置く必要がありますね.
ああ, 線分だったら Line より LineSegment (もしくは Segment) の方がいいかもしれない.
というか, この辺はほとんど数学の問題であって, 全く C には関係ない気がする.

No.7 回答者： qbr2 回答日時： 2009/08/25 15:44

No.5の方の内容に補足させてください。

ABの線分上にCが存在するかどうかなので、
θ=180度で、cosθ = -1　の場合　じゃないかと思います。

ABの直線上にCがあるかだと、cosθ = -1 or 1 ですね。

No.6 回答者： titokani 回答日時： 2009/08/25 15:42

計算方法についてはみなさんがおっしゃるとおりですが、関数の使われ方によっては誤差の処理が必要でしょうね。

No.5 回答者： episteme 回答日時： 2009/08/25 15:18

> ところでその"二つのベクトル AB と AC の内積から求められる"とは

> 具体的にどのようにどのような公式になるのでしょうか？

ベクトルABを x, ベクトルACを y とします。
ｘ・y = |x||y|cosθ なので
cosθ= (ｘ・y) / (|x||y|)
二つのベクトルが重なるときθ=0だから、cosθ = 1 すなわち
(ｘ・y) = |x||y|

No.4 回答者： Tacosan 回答日時： 2009/08/25 14:40

一応確認しておきたいのですが, 判定したいのは

・AB を通る直線上に C が存在するかどうか
・AB を端点とする線分上に C が存在するかどうか
のどちらでしょうか?
後者の方がちょっとだけ面倒なので.

この回答への補足

ご確認ありがとうございます。

当方が行いたいのは後者の方です。
よろしくおねがいします。

補足日時：2009/08/25 14:45

No.3 回答者： qbr2 回答日時： 2009/08/25 13:59

答えがひとつ出ているようですが、内積を使う場合の参考に。

http://w3e.kanazawa-it.ac.jp/math/category/vecto …
「ベクトルのなす角」に関する説明ページです。

Cx, Cy, Czを原点として、Ax,Ay,Azに向かうベクトル、
Cx, Cy, Czを原点として、Bx,By,Bzに向かうベクトル　のなす角が、
180度であれば、直線上にCx,Cy,Czがあると判定できます。

No.1 回答者： episteme 回答日時： 2009/08/25 12:44

二つのベクトル AB と AC の内積から求められるのでは?

この回答への補足

ご回答ありがとうございます。

ところでその"二つのベクトル AB と AC の内積から求められる"とは
具体的にどのようにどのような公式になるのでしょうか？

中学高校の数学を忘れてしまって全然思い出せません。
ご教示頂けませんか？

補足日時：2009/08/25 13:08