JPH08167034A - ２楕円の共通接線作成方法及び図形処理装置 - Google Patents
２楕円の共通接線作成方法及び図形処理装置Info
- Publication number
- JPH08167034A JPH08167034A JP30756794A JP30756794A JPH08167034A JP H08167034 A JPH08167034 A JP H08167034A JP 30756794 A JP30756794 A JP 30756794A JP 30756794 A JP30756794 A JP 30756794A JP H08167034 A JPH08167034 A JP H08167034A
- Authority
- JP
- Japan
- Prior art keywords
- ellipse
- ellipses
- point
- tangent
- straight line
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Image Generation (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】
【目的】 安定して矛盾なく２楕円の共通接線を作成す
る２楕円の共通接線作成方法及び図形処理装置を提供す
る。 【構成】 計算機を使用して２楕円の共通接線を作成す
る２楕円の共通接線作成方法であって、一方の楕円の中
心点が原点と一致し、且つ長軸ベクトルがＸ軸と一致す
るような座標系に、２楕円の座標変換を行うステップＳ
１，Ｓ２と、前記変換後の２楕円の交点を算出して、求
まった交点の数と２楕円の位置とから２楕円の位相関係
を求めるステップＳ３，Ｓ４と、前記求めた位相関係に
より共通接線の本数を定めるステップＳ５と、前記座標
変換された２楕円から、所定のパラメトリック表現の評
価関数に基づいて、前記共通接線の本数だけのパラメト
リック表現の解を求めるステップＳ６と、前記求まった
パラメトリック表現の解を実空間の座標値に変換するス
テップＳ７とを備える。
る２楕円の共通接線作成方法及び図形処理装置を提供す
る。 【構成】 計算機を使用して２楕円の共通接線を作成す
る２楕円の共通接線作成方法であって、一方の楕円の中
心点が原点と一致し、且つ長軸ベクトルがＸ軸と一致す
るような座標系に、２楕円の座標変換を行うステップＳ
１，Ｓ２と、前記変換後の２楕円の交点を算出して、求
まった交点の数と２楕円の位置とから２楕円の位相関係
を求めるステップＳ３，Ｓ４と、前記求めた位相関係に
より共通接線の本数を定めるステップＳ５と、前記座標
変換された２楕円から、所定のパラメトリック表現の評
価関数に基づいて、前記共通接線の本数だけのパラメト
リック表現の解を求めるステップＳ６と、前記求まった
パラメトリック表現の解を実空間の座標値に変換するス
テップＳ７とを備える。
Description
【０００１】
【産業上の利用分野】本発明はＣＡＤ等における２楕円
の共通接線作成方法及び図形処理装置に関するものであ
る。
の共通接線作成方法及び図形処理装置に関するものであ
る。
【０００２】
【従来の技術】従来から、ＣＡＤ等で計算機を使って２
楕円の共通接線を算出するる方法として、４次の評価関
数を作成してその評価関数を代数的に解く方法が知られ
ている。
楕円の共通接線を算出するる方法として、４次の評価関
数を作成してその評価関数を代数的に解く方法が知られ
ている。
【０００３】
【発明が解決しようとする課題】しかしながら、従来の
方法だと２楕円の位相関係を考慮せずに処理を行ってい
るので、解の共通接線の個数に位相関係との矛盾が生じ
たり、かならずしも安定的には解が求まらないという問
題点があった。本発明は、前記従来の欠点を除去し、安
定して矛盾なく２楕円の共通接線を作成する２楕円の共
通接線作成方法及び図形処理装置を提供する。
方法だと２楕円の位相関係を考慮せずに処理を行ってい
るので、解の共通接線の個数に位相関係との矛盾が生じ
たり、かならずしも安定的には解が求まらないという問
題点があった。本発明は、前記従来の欠点を除去し、安
定して矛盾なく２楕円の共通接線を作成する２楕円の共
通接線作成方法及び図形処理装置を提供する。
【０００４】
【課題を解決するための手段】この課題を解決するため
に、本発明の２楕円の共通接線作成方法は、計算機を使
用して２楕円の共通接線を作成する２楕円の共通接線作
成方法であって、一方の楕円の中心点が原点と一致し、
且つ長軸ベクトルがＸ軸と一致するような座標系に、２
楕円の座標変換を行うステップと、前記変換後の２楕円
の交点を算出して、求まった交点の数と２楕円の位置と
から２楕円の位相関係を求めるステップと、前記求めた
位相関係により共通接線の本数を定めるステップと、前
記座標変換された２楕円から、所定のパラメトリック表
現の評価関数に基づいて、前記共通接線の本数だけのパ
ラメトリック表現の解を求めるステップと、前記求まっ
たパラメトリック表現の解を実空間の座標値に変換する
ステップとを備えることを特徴とする。
に、本発明の２楕円の共通接線作成方法は、計算機を使
用して２楕円の共通接線を作成する２楕円の共通接線作
成方法であって、一方の楕円の中心点が原点と一致し、
且つ長軸ベクトルがＸ軸と一致するような座標系に、２
楕円の座標変換を行うステップと、前記変換後の２楕円
の交点を算出して、求まった交点の数と２楕円の位置と
から２楕円の位相関係を求めるステップと、前記求めた
位相関係により共通接線の本数を定めるステップと、前
記座標変換された２楕円から、所定のパラメトリック表
現の評価関数に基づいて、前記共通接線の本数だけのパ
ラメトリック表現の解を求めるステップと、前記求まっ
たパラメトリック表現の解を実空間の座標値に変換する
ステップとを備えることを特徴とする。
【０００５】ここで、前記２楕円の座標変換を行うステ
ップは、一方の楕円の中心点が原点と一致し、且つ長軸
ベクトルがＸ軸と一致するような変換を行う変換マトリ
ックスを生成するステップと、前記変換マトリックスに
より２楕円の座標変換を行うステップとを有する。ま
た、前記解を求めるステップでは、中心点が原点となっ
た一方の楕円上の一点Ｐ１から他方の楕円に引かれた２
本の接線Ｌ１，Ｌ２が、前記一方の楕円の接線となる条
件を満たした場合の前記評価関数の変数を解とする。ま
た、前記評価関数は、前記点Ｐ１の原点回りの角度をθ
とした場合に、 ｆ( θ)=(PY2-b・sinθ)a・sinθ+b・cosθ(PX2-a・cosθ) …(1) ｇ( θ)=(PY3-b・sinθ)a・sinθ+b・cosθ(PX3-a・cosθ) …(2) ｆ( θ) は直線Ｌ１に対する評価関数であり、ｇ( θ)
は直線Ｌ２に対する評価関数である。
ップは、一方の楕円の中心点が原点と一致し、且つ長軸
ベクトルがＸ軸と一致するような変換を行う変換マトリ
ックスを生成するステップと、前記変換マトリックスに
より２楕円の座標変換を行うステップとを有する。ま
た、前記解を求めるステップでは、中心点が原点となっ
た一方の楕円上の一点Ｐ１から他方の楕円に引かれた２
本の接線Ｌ１，Ｌ２が、前記一方の楕円の接線となる条
件を満たした場合の前記評価関数の変数を解とする。ま
た、前記評価関数は、前記点Ｐ１の原点回りの角度をθ
とした場合に、 ｆ( θ)=(PY2-b・sinθ)a・sinθ+b・cosθ(PX2-a・cosθ) …(1) ｇ( θ)=(PY3-b・sinθ)a・sinθ+b・cosθ(PX3-a・cosθ) …(2) ｆ( θ) は直線Ｌ１に対する評価関数であり、ｇ( θ)
は直線Ｌ２に対する評価関数である。
【０００６】なおここで、ａ＝標準位置にある楕円Ｅ１
の長軸ベクトルの長さ、ｂ＝標準位置にある楕円Ｅ１の
短軸ベクトルの長さ、PX2 ＝パラメータθに対する標準
位置の楕円Ｅ１上の一点Ｐ１から他方の楕円Ｅ２への接
線上の一点のｘ成分（直線Ｌ１）、PY2 ＝パラメータθ
に対する標準位置の楕円Ｅ１上の一点Ｐ１から他方の楕
円Ｅ２への接線上の一点のｙ成分（直線Ｌ１）、PX3 ＝
パラメータθに対する標準位置の楕円Ｅ１上の一点Ｐ１
から他方の楕円Ｅ２への接線上の一点のｘ成分（直線Ｌ
２）、PY3 ＝パラメータθに対する標準位置の楕円Ｅ１
上の一点Ｐ１から他方の楕円Ｅ２への接線上の一点のｙ
成分（直線Ｌ２）である。
の長軸ベクトルの長さ、ｂ＝標準位置にある楕円Ｅ１の
短軸ベクトルの長さ、PX2 ＝パラメータθに対する標準
位置の楕円Ｅ１上の一点Ｐ１から他方の楕円Ｅ２への接
線上の一点のｘ成分（直線Ｌ１）、PY2 ＝パラメータθ
に対する標準位置の楕円Ｅ１上の一点Ｐ１から他方の楕
円Ｅ２への接線上の一点のｙ成分（直線Ｌ１）、PX3 ＝
パラメータθに対する標準位置の楕円Ｅ１上の一点Ｐ１
から他方の楕円Ｅ２への接線上の一点のｘ成分（直線Ｌ
２）、PY3 ＝パラメータθに対する標準位置の楕円Ｅ１
上の一点Ｐ１から他方の楕円Ｅ２への接線上の一点のｙ
成分（直線Ｌ２）である。
【０００７】また、前記２楕円の共通接線の本数を定め
るステップの後に、解の存在範囲を確定するステップを
更に有する。また、前記接線の数を確定するステップの
後に、接点での接線が２楕円の共通接線となる場合にお
いて、そこでの接線を直接算出するステップを更に有す
る。又、本発明の図形処理装置は、予め作成された図形
に関連して新たな図形を作成する図形処理装置におい
て、２楕円が与えられた場合に、一方の楕円の中心点が
原点と一致し、且つ長軸ベクトルがＸ軸と一致するよう
な座標系に、２楕円の座標変換を行う座標変換手段と、
前記変換後の２楕円の交点を算出して、求まった交点の
数と２楕円の位置とから２楕円の位相関係を求め、前記
求めた位相関係により共通接線の本数を定める接線本数
決定手段と、前記座標変換された２楕円から、所定のパ
ラメトリック表現の評価関数に基づいて、前記共通接線
の本数だけのパラメトリック表現の解を求める直線評価
手段と、前記求まったパラメトリック表現の解を実空間
の座標値に変換して２楕円の共通接線として出力する共
通接線作成手段とを備えることを特徴とする。
るステップの後に、解の存在範囲を確定するステップを
更に有する。また、前記接線の数を確定するステップの
後に、接点での接線が２楕円の共通接線となる場合にお
いて、そこでの接線を直接算出するステップを更に有す
る。又、本発明の図形処理装置は、予め作成された図形
に関連して新たな図形を作成する図形処理装置におい
て、２楕円が与えられた場合に、一方の楕円の中心点が
原点と一致し、且つ長軸ベクトルがＸ軸と一致するよう
な座標系に、２楕円の座標変換を行う座標変換手段と、
前記変換後の２楕円の交点を算出して、求まった交点の
数と２楕円の位置とから２楕円の位相関係を求め、前記
求めた位相関係により共通接線の本数を定める接線本数
決定手段と、前記座標変換された２楕円から、所定のパ
ラメトリック表現の評価関数に基づいて、前記共通接線
の本数だけのパラメトリック表現の解を求める直線評価
手段と、前記求まったパラメトリック表現の解を実空間
の座標値に変換して２楕円の共通接線として出力する共
通接線作成手段とを備えることを特徴とする。
【０００８】ここで、前記直線評価手段は、中心点が原
点となった一方の楕円上の一点Ｐ１から他方の楕円に引
かれた２本の接線Ｌ１，Ｌ２が、前記一方の楕円の接線
となる条件を満たした場合の前記評価関数の変数を解と
し、前記評価関数は、前記点Ｐ１の原点回りの角度をθ
とした場合に、 ｆ( θ)=(PY2-b・sinθ)a・sinθ+b・cosθ(PX2-a・cosθ) …(1) ｇ( θ)=(PY3-b・sinθ)a・sinθ+b・cosθ(PX3-a・cosθ) …(2) ｆ( θ) は直線Ｌ１に対する評価関数であり、ｇ( θ)
は直線Ｌ２に対する評価関数である。
点となった一方の楕円上の一点Ｐ１から他方の楕円に引
かれた２本の接線Ｌ１，Ｌ２が、前記一方の楕円の接線
となる条件を満たした場合の前記評価関数の変数を解と
し、前記評価関数は、前記点Ｐ１の原点回りの角度をθ
とした場合に、 ｆ( θ)=(PY2-b・sinθ)a・sinθ+b・cosθ(PX2-a・cosθ) …(1) ｇ( θ)=(PY3-b・sinθ)a・sinθ+b・cosθ(PX3-a・cosθ) …(2) ｆ( θ) は直線Ｌ１に対する評価関数であり、ｇ( θ)
は直線Ｌ２に対する評価関数である。
【０００９】なおここで、ａ＝標準位置にある楕円Ｅ１
の長軸ベクトルの長さ、ｂ＝標準位置にある楕円Ｅ１の
短軸ベクトルの長さ、PX2 ＝パラメータθに対する標準
位置の楕円Ｅ１上の一点Ｐ１から他方の楕円Ｅ２への接
線上の一点のｘ成分（直線Ｌ１）、PY2 ＝パラメータθ
に対する標準位置の楕円Ｅ１上の一点Ｐ１から他方の楕
円Ｅ２への接線上の一点のｙ成分（直線Ｌ１）、PX3 ＝
パラメータθに対する標準位置の楕円Ｅ１上の一点Ｐ１
から他方の楕円Ｅ２への接線上の一点のｘ成分（直線Ｌ
２）、PY3 ＝パラメータθに対する標準位置の楕円Ｅ１
上の一点Ｐ１から他方の楕円Ｅ２への接線上の一点のｙ
成分（直線Ｌ２）である。
の長軸ベクトルの長さ、ｂ＝標準位置にある楕円Ｅ１の
短軸ベクトルの長さ、PX2 ＝パラメータθに対する標準
位置の楕円Ｅ１上の一点Ｐ１から他方の楕円Ｅ２への接
線上の一点のｘ成分（直線Ｌ１）、PY2 ＝パラメータθ
に対する標準位置の楕円Ｅ１上の一点Ｐ１から他方の楕
円Ｅ２への接線上の一点のｙ成分（直線Ｌ１）、PX3 ＝
パラメータθに対する標準位置の楕円Ｅ１上の一点Ｐ１
から他方の楕円Ｅ２への接線上の一点のｘ成分（直線Ｌ
２）、PY3 ＝パラメータθに対する標準位置の楕円Ｅ１
上の一点Ｐ１から他方の楕円Ｅ２への接線上の一点のｙ
成分（直線Ｌ２）である。
【００１０】
【実施例】図１は、本発明の共通接線作成方法を実現す
る図形処理装置のブロック図である。バス１（制御線、
データ線及びアドレス線を含む）には、中央処理装置
（ＣＰＵ）２、リード・オンリ・メモリ（ＲＯＭ）３、
ランダム・アクセス・メモリ（ＲＡＭ）４、入力インタ
ーフェイス５を介して入力装置６、ＣＲＴインターフェ
イス７を介してＣＲＴ８、外部記憶装置インターフェイ
ス９を介して、磁気ディスク、磁気テープ等の外部記憶
装置１０が接続されている。ＣＰＵ２は、ＲＯＭ３に記
憶されたプログラムに従って、ＲＡＭ４を一時記憶装置
として種々の処理及び制御、例えば図形入力制御，図形
表示処理，ピック処理，交点算出などを行う。
る図形処理装置のブロック図である。バス１（制御線、
データ線及びアドレス線を含む）には、中央処理装置
（ＣＰＵ）２、リード・オンリ・メモリ（ＲＯＭ）３、
ランダム・アクセス・メモリ（ＲＡＭ）４、入力インタ
ーフェイス５を介して入力装置６、ＣＲＴインターフェ
イス７を介してＣＲＴ８、外部記憶装置インターフェイ
ス９を介して、磁気ディスク、磁気テープ等の外部記憶
装置１０が接続されている。ＣＰＵ２は、ＲＯＭ３に記
憶されたプログラムに従って、ＲＡＭ４を一時記憶装置
として種々の処理及び制御、例えば図形入力制御，図形
表示処理，ピック処理，交点算出などを行う。
【００１１】本例では、ＲＯＭ３に、２楕円の共通接線
を作成する接線作成プログラム３ａと作成された線が接
線か否かを評価する評価関数３ｂが、ＲＡＭ４に、一方
の楕円の中心を原点とし長軸をＸ軸とするする座標への
変換マトリックス４ａ，２楕円の交点数４ｂ，楕円及び
接線の実空間座標値４ｃ，上記変換マトリックス４ａに
より変換された楕円及び接線の変換座標値４ｄが格納さ
れる。
を作成する接線作成プログラム３ａと作成された線が接
線か否かを評価する評価関数３ｂが、ＲＡＭ４に、一方
の楕円の中心を原点とし長軸をＸ軸とするする座標への
変換マトリックス４ａ，２楕円の交点数４ｂ，楕円及び
接線の実空間座標値４ｃ，上記変換マトリックス４ａに
より変換された楕円及び接線の変換座標値４ｄが格納さ
れる。
【００１２】入力装置６はキーボード，タブレット，マ
ウス等であり、図形データの入力を行なうが、この図形
データはホストコンピュータから受けてもよい。ＣＲＴ
８は必要に応じて複数のビットマップ・プレーン等を含
んでおり図形を表示する。以下、図１〜図４を参照して
本発明の好適な実施例を説明するが、以下に述べる処理
及び判断はＣＰＵ２がＲＯＭ３のプログラムに従って行
なう。
ウス等であり、図形データの入力を行なうが、この図形
データはホストコンピュータから受けてもよい。ＣＲＴ
８は必要に応じて複数のビットマップ・プレーン等を含
んでおり図形を表示する。以下、図１〜図４を参照して
本発明の好適な実施例を説明するが、以下に述べる処理
及び判断はＣＰＵ２がＲＯＭ３のプログラムに従って行
なう。
【００１３】図２の本実施例の流れを示すフローチャー
トに従って、２楕円の共通接線作成処理の説明を行な
う。 （ステップＳ１：変換マトリックスの算出）本実施例で
は、まず対象２楕円のうち、一方の楕円の中心点が原点
に一致し、その楕円の長軸ベクトルがＸ軸に一致するよ
うな変換マトリックスを求め、ＲＡＭ４または外部記憶
インターフェイス９を介して外部記憶装置１０に登録す
る。
トに従って、２楕円の共通接線作成処理の説明を行な
う。 （ステップＳ１：変換マトリックスの算出）本実施例で
は、まず対象２楕円のうち、一方の楕円の中心点が原点
に一致し、その楕円の長軸ベクトルがＸ軸に一致するよ
うな変換マトリックスを求め、ＲＡＭ４または外部記憶
インターフェイス９を介して外部記憶装置１０に登録す
る。
【００１４】（ステップＳ２：２楕円に変換を行なう）
ステップＳ１で求めた変換マトリックスによる変換を対
象２楕円に対し行なう。このステップによって一方の楕
円は中心点が原点にあり、長軸ベクトルがＸ軸と一致す
るような標準位置に移動し、他方の楕円は標準位置より
元の２楕円の相対的な位置関係分離れたところに移動す
る（図４参照）。
ステップＳ１で求めた変換マトリックスによる変換を対
象２楕円に対し行なう。このステップによって一方の楕
円は中心点が原点にあり、長軸ベクトルがＸ軸と一致す
るような標準位置に移動し、他方の楕円は標準位置より
元の２楕円の相対的な位置関係分離れたところに移動す
る（図４参照）。
【００１５】（ステップＳ３：２楕円の交点を算出す
る）２楕円の交点を算出する。なおこの時接点は交点と
みなしてよい。そしてその交点数をＲＡＭ４又は外部記
憶インターフェイス９を介して外部記憶装置１０に登録
する。 （ステップＳ４：２楕円の位相関係を算出する）２楕円
の位相関係を、ＲＡＭ４又は外部記憶装置１０に登録し
てある交点数を基に算出する。例えば、ステップＳ３で
算出された交点数が“０”の場合、可能な位相関係は図
３に示すa-1 ないしa-2 になる。この２つを分別するに
は、一方の楕円が他方の楕円に含まれるか含まれないか
の判定を行えばよい。
る）２楕円の交点を算出する。なおこの時接点は交点と
みなしてよい。そしてその交点数をＲＡＭ４又は外部記
憶インターフェイス９を介して外部記憶装置１０に登録
する。 （ステップＳ４：２楕円の位相関係を算出する）２楕円
の位相関係を、ＲＡＭ４又は外部記憶装置１０に登録し
てある交点数を基に算出する。例えば、ステップＳ３で
算出された交点数が“０”の場合、可能な位相関係は図
３に示すa-1 ないしa-2 になる。この２つを分別するに
は、一方の楕円が他方の楕円に含まれるか含まれないか
の判定を行えばよい。
【００１６】この判定は以下のように行なう。２楕円の
中心点を結ぶ線分を作成し、その線分の長さがもし幾何
学的許容誤差ε以内であれば２楕円の中心点は一致して
いるとみなす。この場合、すでに２楕円は交点をもたな
いということがわかっているので、２楕円の位相関係は
図３に示すa-2 になる。また、２楕円の中心点を結ぶ線
分の長さが幾何学的許容誤差ε以上であれば、その線分
と２楕円のうち大きい長軸を持つ楕円との交点算出を行
なう。もし交点が存在するならば２楕円の位相関係は図
３に示すa-1 になり、もし交点が存在しなければa-2 と
する。
中心点を結ぶ線分を作成し、その線分の長さがもし幾何
学的許容誤差ε以内であれば２楕円の中心点は一致して
いるとみなす。この場合、すでに２楕円は交点をもたな
いということがわかっているので、２楕円の位相関係は
図３に示すa-2 になる。また、２楕円の中心点を結ぶ線
分の長さが幾何学的許容誤差ε以上であれば、その線分
と２楕円のうち大きい長軸を持つ楕円との交点算出を行
なう。もし交点が存在するならば２楕円の位相関係は図
３に示すa-1 になり、もし交点が存在しなければa-2 と
する。
【００１７】ステップＳ３で算出された交点数が“１”
の場合、可能な位相関係は図３に示すb-1 ないしb-2 に
なる。この２つを分別するには、一方の楕円が他方の楕
円に含まれるか含まれないかの判定を行えば良い。この
判定方法は、上記した交点数が“０”の時と同じように
行える。ステップＳ３で算出された交点数が“２”の場
合、可能な位相関係は図３に示すc-1 になる。
の場合、可能な位相関係は図３に示すb-1 ないしb-2 に
なる。この２つを分別するには、一方の楕円が他方の楕
円に含まれるか含まれないかの判定を行えば良い。この
判定方法は、上記した交点数が“０”の時と同じように
行える。ステップＳ３で算出された交点数が“２”の場
合、可能な位相関係は図３に示すc-1 になる。
【００１８】ステップＳ３で算出された交点数が“３”
の場合、可能な位相関係は図３に示すd-1 になる。ステ
ップＳ３で算出された交点数が“４”の場合、可能な位
相関係は図３に示すe-1 になる。ステップＳ３で２楕円
が同じだと判定された場合、可能な位相関係は図３に示
すｆ-1になる。
の場合、可能な位相関係は図３に示すd-1 になる。ステ
ップＳ３で算出された交点数が“４”の場合、可能な位
相関係は図３に示すe-1 になる。ステップＳ３で２楕円
が同じだと判定された場合、可能な位相関係は図３に示
すｆ-1になる。
【００１９】（ステップＳ５：共通接線の本数を定め
る）ステップＳ４で算出された２楕円の位相関係を基
に、求めるべき共通接点の本数を定める。 ケース a-1 a-2 b-1 b-2 c-1 d-1 e-1 f-1 本数 ４ ０ ３ １ ２ ３ ４ ０ （ステップＳ６：パラメータ解を算出する）ステップＳ
５で定まった接線の本数分、下記の評価関数(1) 及び評
価関数(2)が“０”となるθの解を割線法を用いて算出
する。以下、割線法による解の算出を、図４の例を参照
にして説明する。
る）ステップＳ４で算出された２楕円の位相関係を基
に、求めるべき共通接点の本数を定める。 ケース a-1 a-2 b-1 b-2 c-1 d-1 e-1 f-1 本数 ４ ０ ３ １ ２ ３ ４ ０ （ステップＳ６：パラメータ解を算出する）ステップＳ
５で定まった接線の本数分、下記の評価関数(1) 及び評
価関数(2)が“０”となるθの解を割線法を用いて算出
する。以下、割線法による解の算出を、図４の例を参照
にして説明する。
【００２０】 ｆ( θ)=(PY2-b・sinθ)a・sinθ+b・cosθ(PX2-a・cosθ) …(1) ｇ( θ)=(PY3-b・sinθ)a・sinθ+b・cosθ(PX3-a・cosθ) …(2) ｆ( θ) は直線Ｌ１に対する評価関数であり、ｇ( θ)
は直線Ｌ２に対する評価関数である。なおここで、 ａ＝標準位置にある楕円Ｅ１の長軸ベクトルの長さ ｂ＝標準位置にある楕円Ｅ１の短軸ベクトルの長さ PX2 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｘ成分
（直線Ｌ１） PY2 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｙ成分
（直線Ｌ１） PX3 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｘ成分
（直線Ｌ２） PY3 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｙ成分
（直線Ｌ２） 従って、評価関数(1) 及び評価関数(2) は、楕円Ｅ１上
の点Ｐ１から他方の楕円Ｅ２に引かれた接線ベクトルと
点Ｐ１での法線ベクトルとの内積であり、共通接線の場
合は、ｆ( θ) ，ｇ( θ) は“０”となる。
は直線Ｌ２に対する評価関数である。なおここで、 ａ＝標準位置にある楕円Ｅ１の長軸ベクトルの長さ ｂ＝標準位置にある楕円Ｅ１の短軸ベクトルの長さ PX2 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｘ成分
（直線Ｌ１） PY2 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｙ成分
（直線Ｌ１） PX3 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｘ成分
（直線Ｌ２） PY3 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｙ成分
（直線Ｌ２） 従って、評価関数(1) 及び評価関数(2) は、楕円Ｅ１上
の点Ｐ１から他方の楕円Ｅ２に引かれた接線ベクトルと
点Ｐ１での法線ベクトルとの内積であり、共通接線の場
合は、ｆ( θ) ，ｇ( θ) は“０”となる。
【００２１】（ステップＳ７：パラメータ解を実空間の
座標値に変換する）ステップＳ６で算出した、標準位置
にある楕円Ｅ１上のパラメータθにより移動する点Ｐ１
が、その楕円Ｅ１の接線である場合が解、すなわち両楕
円の共通接線である。上記パラメータは標準位置に移動
する前の楕円に関しても有効であるため、そのパラメー
タθで表わされた接線を実空間の座標値に変換する。
座標値に変換する）ステップＳ６で算出した、標準位置
にある楕円Ｅ１上のパラメータθにより移動する点Ｐ１
が、その楕円Ｅ１の接線である場合が解、すなわち両楕
円の共通接線である。上記パラメータは標準位置に移動
する前の楕円に関しても有効であるため、そのパラメー
タθで表わされた接線を実空間の座標値に変換する。
【００２２】算出された実空間の接線座標値を、ＲＡＭ
４または外部記憶インターフェイス９を介して外部記憶
装置１０に登録する。尚、図２には図示しなかったが、
算出された実空間の接線座標値に基づいて、ＣＲＴ８の
画面上に２楕円の共通接線が描画される。 ＜本実施例の変形例＞上記実施例の図２のステップＳ５
とＳ６との間に、以下のステップＳ５−１及び／又はＳ
５−２を追加することで、算出すべき解の数を減らすこ
とが可能である。
４または外部記憶インターフェイス９を介して外部記憶
装置１０に登録する。尚、図２には図示しなかったが、
算出された実空間の接線座標値に基づいて、ＣＲＴ８の
画面上に２楕円の共通接線が描画される。 ＜本実施例の変形例＞上記実施例の図２のステップＳ５
とＳ６との間に、以下のステップＳ５−１及び／又はＳ
５−２を追加することで、算出すべき解の数を減らすこ
とが可能である。
【００２３】（ステップＳ５−２：解の存在範囲を確定
する）ステップＳ５で定まった接線の本数分、評価関数
(1) 及び評価関数(2) を計算し、ステップＳ６の前処理
として２分法により解の存在範囲の確定を予め行なう。 （ステップＳ５−１：接点での接線の算出）ステップＳ
５において接線の数を確定したが、さらにb-1,b-2 及び
d-1 の場合において接点での接線が２楕円の共通接線と
なるためそこでの接線を直接(3) 式によって求める。そ
して収束演算で算出すべき解の数を減らす。
する）ステップＳ５で定まった接線の本数分、評価関数
(1) 及び評価関数(2) を計算し、ステップＳ６の前処理
として２分法により解の存在範囲の確定を予め行なう。 （ステップＳ５−１：接点での接線の算出）ステップＳ
５において接線の数を確定したが、さらにb-1,b-2 及び
d-1 の場合において接点での接線が２楕円の共通接線と
なるためそこでの接線を直接(3) 式によって求める。そ
して収束演算で算出すべき解の数を減らす。
【００２４】楕円上の点Ｐ１(a・cosθ,b・sinθ) におけ
る接線 x・cosθ/a + y・sinθ/b = 1 …(3) 尚、ステップＳ４における位相関係の判定で、２楕円の
交点が１点あるいは０点の場合において、一方の楕円の
中心点より半無限直線を引き、その半無限直線が他方の
楕円と干渉するかしないかによって、一方の楕円が他方
の楕円の中に含まれるかの判定を行なうことができる。
もしその半無限直線と他方の楕円との交点が一点存在す
れば（接点は交点とみなさない）、a-1/b-1 ではなくa-
2/b-2 と判定する。
る接線 x・cosθ/a + y・sinθ/b = 1 …(3) 尚、ステップＳ４における位相関係の判定で、２楕円の
交点が１点あるいは０点の場合において、一方の楕円の
中心点より半無限直線を引き、その半無限直線が他方の
楕円と干渉するかしないかによって、一方の楕円が他方
の楕円の中に含まれるかの判定を行なうことができる。
もしその半無限直線と他方の楕円との交点が一点存在す
れば（接点は交点とみなさない）、a-1/b-1 ではなくa-
2/b-2 と判定する。
【００２５】又、ステップＳ６において割線法を用いて
解を算出したが、他の収束演算を用いても同様の処理を
行なうことができる。尚、本発明は、複数の機器から構
成されるシステムに適用しても、１つの機器から成る装
置に適用しても良い。また、本発明はシステム或は装置
にプログラムを供給することによって達成される場合に
も適用できることはいうまでもない。
解を算出したが、他の収束演算を用いても同様の処理を
行なうことができる。尚、本発明は、複数の機器から構
成されるシステムに適用しても、１つの機器から成る装
置に適用しても良い。また、本発明はシステム或は装置
にプログラムを供給することによって達成される場合に
も適用できることはいうまでもない。
【００２６】
【発明の効果】以上説明したような２楕円の共通接線作
成を行なうことにより、共通接線の解の個数の位相関係
との矛盾を防ぐことができ、安定的に共通接線が求まる
という効果がある。
成を行なうことにより、共通接線の解の個数の位相関係
との矛盾を防ぐことができ、安定的に共通接線が求まる
という効果がある。
【図１】本発明を実現する図形処理装置のブロック図で
ある。
ある。
【図２】本発明の好適な実施例を示す流れ図である。
【図３】位相関係によって種別した２楕円の位置関係図
である。
である。
【図４】一方の楕円を標準位置に移動し、あるパラメー
タθより他方の楕円Ｅ２への接線を引いた例を示す図で
ある。
タθより他方の楕円Ｅ２への接線を引いた例を示す図で
ある。
【図５】本発明の好適な他の実施例を示す流れ図であ
る。
る。
Claims (8)
- 【請求項１】 計算機を使用して２楕円の共通接線を作
成する２楕円の共通接線作成方法であって、 一方の楕円の中心点が原点と一致し、且つ長軸ベクトル
がＸ軸と一致するような座標系に、２楕円の座標変換を
行うステップと、 前記変換後の２楕円の交点を算出して、求まった交点の
数と２楕円の位置とから２楕円の位相関係を求めるステ
ップと、 前記求めた位相関係により共通接線の本数を定めるステ
ップと、 前記座標変換された２楕円から、所定のパラメトリック
表現の評価関数に基づいて、前記共通接線の本数だけの
パラメトリック表現の解を求めるステップと、 前記求まったパラメトリック表現の解を実空間の座標値
に変換するステップとを備えることを特徴とする２楕円
の共通接線作成方法。 - 【請求項２】 前記２楕円の座標変換を行うステップ
は、 一方の楕円の中心点が原点と一致し、且つ長軸ベクトル
がＸ軸と一致するような変換を行う変換マトリックスを
生成するステップと、 前記変換マトリックスにより２楕円の座標変換を行うス
テップとを有することを特徴とする請求項１記載の２楕
円の共通接線作成方法。 - 【請求項３】 前記解を求めるステップでは、中心点が
原点となった一方の楕円上の一点Ｐ１から他方の楕円に
引かれた２本の接線Ｌ１，Ｌ２が、前記一方の楕円の接
線となる条件を満たした場合の前記評価関数の変数を解
とすることを特徴とする請求項１記載の２楕円の共通接
線作成方法。 - 【請求項４】 前記評価関数は、前記点Ｐ１の原点回り
の角度をθとした場合に、 ｆ( θ)=(PY2-b・sinθ)a・sinθ+b・cosθ(PX2-a・cosθ) …(1) ｇ( θ)=(PY3-b・sinθ)a・sinθ+b・cosθ(PX3-a・cosθ) …(2) ｆ( θ) は直線Ｌ１に対する評価関数であり、ｇ( θ)
は直線Ｌ２に対する評価関数であることを特徴とする請
求項３記載の２楕円の共通接線作成方法。なおここで、 ａ＝標準位置にある楕円Ｅ１の長軸ベクトルの長さ、 ｂ＝標準位置にある楕円Ｅ１の短軸ベクトルの長さ、 PX2 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｘ成分
（直線Ｌ１）、 PY2 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｙ成分
（直線Ｌ１）、 PX3 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｘ成分
（直線Ｌ２）、 PY3 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｙ成分
（直線Ｌ２）である。 - 【請求項５】 前記２楕円の共通接線の本数を定めるス
テップの後に、解の存在範囲を確定するステップを更に
有することを特徴とする請求項１記載の２楕円の共通接
線作成方法。 - 【請求項６】 前記接線の数を確定するステップの後
に、接点での接線が２楕円の共通接線となる場合におい
て、そこでの接線を直接算出するステップを更に有する
ことを特徴とする請求項１または５記載の２楕円の共通
接線作成方法。 - 【請求項７】 予め作成された図形に関連して新たな図
形を作成する図形処理装置において、 ２楕円が与えられた場合に、一方の楕円の中心点が原点
と一致し、且つ長軸ベクトルがＸ軸と一致するような座
標系に、２楕円の座標変換を行う座標変換手段と、 前記変換後の２楕円の交点を算出して、求まった交点の
数と２楕円の位置とから２楕円の位相関係を求め、前記
求めた位相関係により共通接線の本数を定める接線本数
決定手段と、 前記座標変換された２楕円から、所定のパラメトリック
表現の評価関数に基づいて、前記共通接線の本数だけの
パラメトリック表現の解を求める直線評価手段と、 前記求まったパラメトリック表現の解を実空間の座標値
に変換して２楕円の共通接線として出力する共通接線作
成手段とを備えることを特徴とする図形処理装置。 - 【請求項８】 前記直線評価手段は、中心点が原点とな
った一方の楕円上の一点Ｐ１から他方の楕円に引かれた
２本の接線Ｌ１，Ｌ２が、前記一方の楕円の接線となる
条件を満たした場合の前記評価関数の変数を解とし、 前記評価関数は、前記点Ｐ１の原点回りの角度をθとし
た場合に、 ｆ( θ)=(PY2-b・sinθ)a・sinθ+b・cosθ(PX2-a・cosθ) …(1) ｇ( θ)=(PY3-b・sinθ)a・sinθ+b・cosθ(PX3-a・cosθ) …(2) ｆ( θ) は直線Ｌ１に対する評価関数であり、ｇ( θ)
は直線Ｌ２に対する評価関数であることを特徴とする請
求項７記載の図形処理装置。なおここで、 ａ＝標準位置にある楕円Ｅ１の長軸ベクトルの長さ、 ｂ＝標準位置にある楕円Ｅ１の短軸ベクトルの長さ、 PX2 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｘ成分
（直線Ｌ１）、 PY2 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｙ成分
（直線Ｌ１）、 PX3 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｘ成分
（直線Ｌ２）、 PY3 ＝パラメータθに対する標準位置の楕円Ｅ１上の一
点Ｐ１から他方の楕円Ｅ２への接線上の一点のｙ成分
（直線Ｌ２）である。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP30756794A JP3332626B2 (ja) | 1994-12-12 | 1994-12-12 | ２楕円の共通接線作成方法及び図形処理装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP30756794A JP3332626B2 (ja) | 1994-12-12 | 1994-12-12 | ２楕円の共通接線作成方法及び図形処理装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08167034A true JPH08167034A (ja) | 1996-06-25 |
| JP3332626B2 JP3332626B2 (ja) | 2002-10-07 |
Family
ID=17970638
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP30756794A Expired - Fee Related JP3332626B2 (ja) | 1994-12-12 | 1994-12-12 | ２楕円の共通接線作成方法及び図形処理装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3332626B2 (ja) |
-
1994
- 1994-12-12 JP JP30756794A patent/JP3332626B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP3332626B2 (ja) | 2002-10-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6966588B2 (ja) | 交差点の描画方法、装置、サーバ及び記憶媒体 | |
| US20020040256A1 (en) | Method of designing forms of cable clamp and cables using three-dimensional CAD system, and computer readable storage medium storing relevant processes | |
| JP2577510B2 (ja) | 図形オブジェクト処理装置及び方法 | |
| US5471569A (en) | Retrieval space determination system for three-dimensionally shaped parts in a CAD system | |
| JP2809762B2 (ja) | 図形整形装置 | |
| JPH08167034A (ja) | ２楕円の共通接線作成方法及び図形処理装置 | |
| CN100501667C (zh) | 三次元量测程序转换系统及方法 | |
| US6598185B1 (en) | Pattern data inspection method and storage medium | |
| JP3413733B2 (ja) | 形状相互間の干渉チェック方法 | |
| US6957176B2 (en) | Reduction processing method and computer readable storage medium having program stored thereon for causing computer to execute the method | |
| JP2014063406A (ja) | 検証支援プログラム、検証支援方法、および検証支援装置 | |
| JP2662856B2 (ja) | 形状特徴計測装置および方法 | |
| JPS59182348A (ja) | 菊池パタ−ン又は擬菊池パタ−ンの解析方法 | |
| JP3535598B2 (ja) | 図形データ整合方法 | |
| JP3156488B2 (ja) | モデル形式変換装置 | |
| JP3679436B2 (ja) | ３次元上領域における領域内外の点の抽出方法及び装置、及び同一曲線上の点の並び順判定方法及び装置 | |
| JPH01193973A (ja) | ３次元形状の干渉チェック方法 | |
| JPS61158395A (ja) | 高速円弧発生装置 | |
| JPH11305833A (ja) | ランドマーク情報の生成装置 | |
| EP0537685A2 (en) | Method for processing detection of candidate of branched direction | |
| JPH08329268A (ja) | 図形処理方法とその装置 | |
| JPH0895627A (ja) | 形状計測手順の教示方法 | |
| JP2004030135A (ja) | Ｃａｄシステムにおける円の一部の描画方法及びプログラム | |
| JP2765312B2 (ja) | 印刷配線板設計装置 | |
| CN116543011A (zh) | 基于单个轮胎的航向检测方法、装置、设备及存储介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20020628 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080726 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080726 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090726 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090726 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100726 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100726 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110726 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120726 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120726 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130726 Year of fee payment: 11 |
|
| LAPS | Cancellation because of no payment of annual fees |