JP2010183995A - 長軸部先端からの光軸を有する物体における光軸位置と該物体の位置姿勢を定義する手段との3次元相対関係測定方法、装置およびシステム - Google Patents
長軸部先端からの光軸を有する物体における光軸位置と該物体の位置姿勢を定義する手段との3次元相対関係測定方法、装置およびシステム Download PDFInfo
- Publication number
- JP2010183995A JP2010183995A JP2009029216A JP2009029216A JP2010183995A JP 2010183995 A JP2010183995 A JP 2010183995A JP 2009029216 A JP2009029216 A JP 2009029216A JP 2009029216 A JP2009029216 A JP 2009029216A JP 2010183995 A JP2010183995 A JP 2010183995A
- Authority
- JP
- Japan
- Prior art keywords
- optical axis
- tip
- dimensional
- measured
- marker
- 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
Images
Landscapes
- Endoscopes (AREA)
Abstract
【解決手段】 被測定物体10の長軸部11先端からの光軸と第1標識部12との3次元相対関係を測定する3次元相対関係測定装置であって、光軸が当たるターゲット23と3次元形状測定装置により測定可能な第2標識部22とを有する第1の較正用物体20と、3次元形状測定装置により測定可能な第3標識部を有する第2の較正用物体30と、第1の較正用物体を被測定物体10と相対的に移動させる移動手段とを有し、移動手段により第1の較正用物体20を移動させてターゲット23に長軸部11先端からの光軸が当たるようにした後、3次元形状測定装置で第1乃至第3標識部の3次元相対位置を測定することにより被測定物体10の光軸と第1標識部12との3次元相対関係を測定する。
【選択図】 図1
Description
そして、長軸部先端からの光軸を有する物体が硬性内視鏡である場合は、光軸位置を遠方まで精度よく測定することで高い精度の手術ナビゲーションを行うことができるようにすることにある。
3次元形状測定装置と、
長軸部と、長軸部先端からの光軸と、前記3次元形状測定装置により位置および姿勢が測定可能な第1標識部とを有する被測定物体と、
前記被測定物体を固定する固定手段と、
前記被測定物体の長軸部先端からの光軸が当たるターゲットと、前記3次元形状測定装置により位置および立体形状が測定可能な第2標識部とを有し、前記ターゲットの中心座標と前記第2標識部との3次元相対関係が予めわかっているか、または前記3次元形状測定装置により前記ターゲットの中心座標と前記第2標識部との3次元相対関係が測定可能な第1の較正用物体と、
前記3次元形状測定装置により位置および立体形状が測定可能な第3標識部を有し、前記固定手段により前記被測定物体が固定された状態において、前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標と前記第3標識部との3次元相対関係が予めわかっているか、または前記3次元形状測定装置により前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸上に含まれる座標と前記第3標識部との3次元相対関係が測定可能な第2の較正用物体と、
を有する3次元形状測定システムにおける、前記被測定物体の長軸部先端からの光軸と前記第1標識部との3次元相対関係を測定する3次元相対関係測定方法であって、
前記被測定物体を前記固定手段に固定し、前記第1の較正用物体のターゲットの中心に前記被測定物体の長軸部先端からの光軸が当たるようにするステップと、
前記3次元形状測定装置により、前記被測定物体、前記第1の較正用物体および前記第2の較正用物体の立体形状データ群を取得するステップと、
前記立体形状データ群から、前記第1標識部の位置および姿勢と前記第2標識部の位置および立体形状と前記第3標識部の位置および立体形状とを算出するステップと、
前記第2標識部の位置および立体形状と、前記第3標識部の位置および立体形状と、前記ターゲットの中心座標と前記第2標識部との3次元相対関係と、前記長軸部先端の座標または前記長軸部先端からの光軸上に含まれる座標と前記第3標識部との3次元相対関係とに基づいて、前記長軸部先端からの光軸の位置を算出するステップと、
前記第1標識部の位置および姿勢と、前記長軸部先端からの光軸の位置とに基づいて、前記被定物体の長軸部先端からの光軸と前記第1標識部との3次元相対関係を算出するステップと、
を有する、3次元相対関係測定方法。
長軸部と、長軸部先端からの光軸と、前記3次元形状測定装置により位置および姿勢が測定可能な第1標識部とを有する被測定物体と、
前記被測定物体を固定する固定手段と、
前記被測定物体の長軸部先端からの光軸が当たるターゲットと、前記3次元形状測定装置により位置および立体形状が測定可能な第2標識部とを有し、前記ターゲットの中心座標と前記第2標識部との3次元相対関係が予めわかっているか、または前記3次元形状測定装置により前記ターゲットの中心座標と前記第2標識部との3次元相対関係が測定可能な第1の較正用物体と、
前記3次元形状測定装置により位置および立体形状が測定可能な第3標識部を有し、前記固定手段により前記被測定物体が固定された状態において、前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標と前記第3標識部との3次元相対関係が予めわかっているか、または前記3次元形状測定装置により前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸上に含まれる座標と前記第3標識部との3次元相対関係が測定可能な第2の較正用物体と、
を有する3次元形状測定システムにおける、前記被測定物体の長軸部先端からの光軸と前記第1標識部との3次元相対関係を測定する3次元相対関係測定プログラムであって、
前記被測定物体を前記固定手段に固定し、前記第1の較正用物体のターゲットの中心に前記被測定物体の長軸部先端からの光軸を一致させた状態で、前記3次元形状測定装置により、前記被測定物体、前記第1の較正用物体および前記第2の較正用物体の立体形状データ群を取得するステップと、
前記立体形状データ群から、前記第1標識部の位置および姿勢と前記第2標識部の位置および立体形状と前記第3標識部の位置および立体形状とを算出するステップと、
前記第2標識部の位置および立体形状と、前記第3標識部の位置および立体形状と、前記ターゲットの中心座標と前記第2標識部との3次元相対関係と、前記長軸部先端の座標または前記長軸部先端からの光軸上に含まれる座標と前記第3標識部との3次元相対関係とに基づいて、前記長軸部先端からの光軸の位置を算出するステップと、
前記第1標識部の位置および姿勢と、前記長軸部先端からの光軸の位置とに基づいて、前記被定物体の長軸部先端からの光軸と前記第1標識部との3次元相対関係を算出するステップと、
を有する、3次元相対関係測定プログラム。
長軸部と、長軸部先端からの光軸と、前記3次元形状測定装置により位置および姿勢が測定可能な第1標識部とを有する被測定物体と、
前記被測定物体を固定する固定手段と、
前記被測定物体の長軸部先端からの光軸が当たるターゲットと、前記3次元形状測定装置により位置および立体形状が測定可能な第2標識部とを有し、前記ターゲットの中心座標と前記第2標識部との3次元相対関係が予めわかっているか、または前記3次元形状測定装置により前記ターゲットの中心座標と前記第2標識部との3次元相対関係が測定可能な第1の較正用物体と、
前記3次元形状測定装置により位置および立体形状が測定可能な第3標識部を有し、前記固定手段により前記被測定物体が固定された状態において、前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標と前記第3標識部との3次元相対関係が予めわかっているか、または前記3次元形状測定装置により前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸上に含まれる座標と前記第3標識部との3次元相対関係が測定可能な第2の較正用物体と、を有し、
前記被測定物体を前記固定手段に固定し、前記第1の較正用物体のターゲットの中心に前記被測定物体の長軸部先端からの光軸を一致させた状態で、前記3次元形状測定装置により、前記被測定物体、前記第1の較正用物体および前記第2の較正用物体の立体形状データ群を取得する手段と、
前記立体形状データ群から、前記第1標識部の位置および姿勢と前記第2標識部の位置および立体形状と前記第3標識部の位置および立体形状とを算出する手段と、
前記第2標識部の位置および立体形状と、前記第3標識部の位置および立体形状と、前記ターゲットの中心座標と前記第2標識部との3次元相対関係と、前記長軸部先端の座標または前記長軸部先端からの光軸上に含まれる座標と前記第3標識部との3次元相対関係とに基づいて、前記長軸部先端からの光軸の位置を算出する手段と、
前記第1標識部の位置および姿勢と、前記長軸部先端からの光軸の位置とに基づいて、前記被定物体の長軸部先端からの光軸と前記第1標識部との3次元相対関係を算出する手段と、
を有する、3次元相対関係測定システム。
前記被測定物体を固定可能な固定手段と、
前記被測定物体の長軸部先端からの光軸が当たるターゲットと、3次元形状測定装置により位置および立体形状が測定可能な第2標識部とを有し、前記ターゲットの中心座標と前記第2標識部との3次元相対関係が予めわかっているか、または3次元形状測定装置により前記ターゲットの中心座標と前記第2標識部との3次元相対関係が測定可能な第1の較正用物体と、
3次元形状測定装置により位置および立体形状が測定可能な第3標識部を有し、前記固定手段により前記被測定物体が固定された状態において、前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標と前記第3標識部との3次元相対関係が予めわかっているか、または3次元形状測定装置により前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸上に含まれる座標と前記第3標識部との3次元相対関係が測定可能な第2の較正用物体と、
前記第1の較正用物体を前記被測定物体と相対的に移動させる移動手段と、を有し、
前記移動手段により前記第1の較正用物体を移動させて前記ターゲットに前記長軸部先端からの光軸が当たるようにした後、3次元形状測定装置で前記第1乃至第3標識部の3次元相対位置を測定することにより、前記被測定物体の前記光軸と前記第1標識部との3次元相対関係を測定可能な3次元相対関係測定装置。
前記第1の較正用物体の第2標識部は、予め立体形状がわかっており、
前記ターゲットの中心座標と前記第2標識部との3次元相対関係は、前記第2標識部の立体形状に関連するベクトルまたは定点位置と、前記ターゲットの中心座標との関係により定義される。
前記第1の較正用物体のターゲットの中心座標と前記第2標識部の定点位置とは等しいか、または前記第1の較正用物体のターゲットの中心座標は前記第2標識部における2つ以上の定点を用いて定義される点である。
前記第2の較正用物体の第3標識部は、予め立体形状がわかっており、
前記長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標と前記第3標識部との3次元相対関係は、前記第3標識部の立体形状に関連するベクトルまたは定点位置と、前記長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標との関係により定義される。
前記長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標と前記第3標識部の定点位置とは等しいか、あるいは前記長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標は前記第3標識部における2つ以上の定点を用いて定義される点である。
前記長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標は、前記第3標識部における2つの定点を結んだ直線または前記第3標識部における定点とベクトルから定まる直線上において、定点から前記第3標識部の立体形状に関連する長さにより定まる点である。
前記長軸部先端からの光軸に含まれる座標は、光を透過する平板に形成されたターゲットの中心座標として定義される。
前記第1の較正用物体および前記第2の較正用物体は、球体、多面体、円柱、円錐のいずれかの形状を一部または全体に有している。
前記被測定物体は体内挿入器具であって、前記長軸部は体内挿入部であり、前記長軸部先端からの光軸は長軸部先端からの撮像方向であり、前記第1標識部は体内に挿入されない部分に設けられている。
前記第2標識部および第3標識部は3次元形状測定装置により立体形状を測定できるものであれば何でも良い。前記第1の較正用物体および第2の較正用物体の立体形状そのものを第2標識部および第3標識部としても良いし、前記第1の較正用物体および第2の較正用物体の一部に球体などを取り付けて第2標識部および第3標識部としても良い。なお演算処理のスピードなどを考慮すると、前記第1の較正用物体においては較正用物体の一部に球体等の第2標識部を取り付ける方が好ましい。
これにより、硬性内視鏡で撮像している部位を画像上に表示しながら手術ナビゲーションを行う手術支援システムにおいて、画面に表示される撮像部位(硬性内視鏡の光軸と体内腔の交点)を精度よく行うことができ、高い精度の手術ナビゲーションを行うことができる。
・少なくとも2つの定点と1つのベクトル
・少なくとも1つの定点と2つのベクトル
が定義可能であってもよく。また定点やベクトルが定義可能な物体としては、球体以外に多面体、円柱、円錐等がある。よって、第1標識部には様々な形態が可能である。
図2のように第1の較正用物体20のターゲット23の中心にピンを取り付けた球体26を刺し込む。球体26は刺し込んだ箇所における直方体21の前面の垂線が球体26の中心を通るように形成されており、球体26の径は、球体22A,22B,22Cとは異なっている。
この状態で3次元形状測定を行い、点群データを取得する。3次元形状測定の方向は、4つの球体22A,22B,22C,26および直方体21の前面を測定できる方向であればどのような方向でもよい。
また3次元形状測定装置は非接触で物体表面の点群データを取得することができるものであれば、どのような測定原理のものでもよい。
取得した点群データから、4つの球体22A,22B,22C,26の点群データを抽出する。4つの球体22A,22B,22C,26の径はコンピュータに記憶されており、球体という情報および径の数値を用いて抽出を行う。抽出方法は特開2004−333371号公報に詳細に記載されているのでそちらを参照する。
抽出した点群データから球体の式(x−a)2+(y−b)2+(z−c)2
= d2を最小2乗法により求める。(a,b,c)が球体の中心座標になる。
取得した点群データから4つの球体22A,22B,22C,26の点群データを除いた残りの点群データの中で球体26の中心座標の近傍にある一部の点群データを抽出し、平面の式a・x+b・y+c・z+1=0に代入して最小2乗法により係数a,b,cを計算する。次にすべての点群データから、この平面からの距離が所定の範囲内にある点群データを抽出する。このとき抽出された点群データが所定数以上あれば、最初に抽出した点群データはすべて直方体21の前面のものと判定して、抽出した点群データすべてを、平面の式a・x+b・y+c・z+1=0に代入して最小2乗法により係数a,b,cを計算する。
もし、抽出された点群データが所定数未満であれば、最初に抽出した点群データは直方体21の前面以外のものを含むと判定して、最初に抽出した点群データからより球体26の中心座標の近傍にあるデータを抽出して同一の処理を行う。この処理を係数a,b,cが計算されるまで行う。
直方体21の前面の法線ベクトルは(a,b,c),(−a,−b,−c)の2つがあり、直方体21の前面から球体26側に向かう法線ベクトルを選定する。それには、3次元形状測定装置の原点座標から球体26の中心座標に向かうベクトルA(ベクトル成分は球体26の中心座標に等しい)と(a,b,c),(−a,−b,−c)の法線ベクトルをベクトルBとしたときの以下の内積の式から角度Θを計算し、角度Θが90度より大きい方の法線ベクトルを選定する。
A・B=|A|・|B|・COSΘ
球体26の中心座標(xd,yd,zd)を通り、直方体21の前面の法線ベクトル(a,
b,c)に平行な直線の方程式(x−xd)/a = (y−yd)/b = (z−zd)/cと、直方体21の前面の平面方程式a・x+b・y+c・z+1=0からなる連立方程式を解いてx,y,zを算出する。このx,y,zが、直線と直方体21の前面の交点座標であり、ターゲット23の中心座標(xt,yt,zt)である。
ここまでの処理により3つの球体22A,22B,22Cの中心座標(xaA,yaA,zaA),(xbA,ybA,zbA),(xcA,ycA,zcA)およびターゲット23の中心座標(xtA,ytA,ztA)を同一座標系で得ることができる。この座標系を座標系Aとする。
なお(1)〜(4)による座標系Aによる3つの球体22A,22B,22Cの中心座標とターゲット23の中心座標との取得は、1度取得したデータを記憶しておけば、第1の較正用物体20を交換しない限り2回目からは実施しなくてよい。
硬性内視鏡10を第2の較正用物体(球体)30に挿入して、ターゲット23に硬性内視鏡10の先端が当たる箇所で硬性内視鏡10を固定する。次に直方体21を移動させ図3に示すように硬性内視鏡10の先端がターゲット23の中心に一致するようにする。
この状態で3次元形状測定を行い、点群データを取得する。3次元形状測定の方向は、図1に示すように光軸位置測定装置1に硬性内視鏡10を固定したとき、4つの球体22A,22B,22C,30および硬性内視鏡10の複数の球体12を測定できる方向であればどのような方向でもよいが、この後、硬性内視鏡10の光軸位置を測定する時と同一にする必要があるので、このとき設定した3次元形状測定装置の位置は動かさない。
上記(2)と同様の方法により3つの球体22A,22B,22Cと第2の較正用物体(球体)30の中心座標(xaB,yaB,zaB),(xbB,ybB,zbB),(xcB,ycB,zcB),(xeB,yeB,zeB)を計算する。このときの座標系を座標系Bとする。
座標系Aによる座標値を座標系Bによる座標値に変換する座標変換関数FABを計算する。
座標変換係数Fは、以下の式における回転成分Mと移動成分(α,β,γ)である。座標変換係数Fを計算するとは、回転成分Mのg11〜g33と移動成分α,β,γを計算することである。
上記式のx’,y’,z’に(xaB,yaB,zaB)を代入し、x,y,zに(xaA,yaA,zaA)を代入した式と、上記式のx’,y’,z’に(xbB,ybB,zbB)を代入し、x,y,zに(xbA,ybA,zbA)を代入した式と、上記式のx’,y’,z’に(xcB,ycB,zcB)を代入し、x,y,zに(xcA,ycA,zcA)を代入した式により9つの式が成立する。
上記(4)で計算された座標系Aによるターゲット23の中心座標(xtA,ytA,ztA)を、(7)で計算された座標変換関数FABにより座標系Bによる座標(xtB,ytB,ztB)に座標変換する。硬性内視鏡10の先端はターゲット23の中心に一致させているので、ターゲット23の中心座標(xtB,ytB,ztB)は、硬性内視鏡10の先端座標である。
上記(6)で計算された第2の較正用物体(球体)30の中心座標(xeB,yeB,zeB)から(8)で計算された硬性内視鏡10の先端座標(xtB,ytB,ztB)に向かうベクトルの成分(xtB−xeB,ytB−yeB,ztB−zeB)を計算する。前述のように第2の較正用物体(球体)30の中心から硬性内視鏡10の先端までの距離は予め測定されているのでこの距離をLとすると、図1に示すように光軸位置測定装置1に硬性内視鏡10を固定したときの、座標系Bによる硬性内視鏡10の先端座標(xsB,ysB,zsB)は、
で計算される。
ここまでの処理により硬性内視鏡10の光軸位置を検出する準備が完了する。
なお、この実施形態においては第1の較正用物体の球体22A,22B,22Cと第2の較正用物体(球体)30を用いて硬性内視鏡10の先端座標を取得しており、請求項における第2の較正用物体の第3標識部は第1の較正用物体の第2標識部も含んでいるとみなす。
図1に示すように光軸位置測定装置1に硬性内視鏡10を固定し直し、硬性内視鏡10の撮像する画面を見ながら第1の較正用物体20を2方向に移動させてターゲット23の中心が撮像画像の中心に来るようにする。これにより硬性内視鏡10の光軸がターゲット23の中心を通るようになる。なお、光軸位置測定装置1に硬性内視鏡10を固定し直すとき光軸位置測定装置1は動かないようにする。
3次元形状測定を行い、点群データを取得する。3次元形状測定装置の位置は、(5)で設定した位置のままとする。即ち、座標系Bによる座標値を取得する。
上記(2)と同様の方法により3つの球体22A,22B,22Cと第2の較正用物体(球体)30の中心座標(xaB’,yaB’,zaB’),(xbB’,ybB’,zbB’),(xcB’,ycB’,zcB’),(xeB’,yeB’,zeB’)と硬性内視鏡10の複数の球体12の中心座標(x1B,y1B,z1B),(x2B,y2B,z2B),(x3B,y3B,z3B),(x4B,y4B,z4B)を計算する。
上記(7)と同様に、座標変換関数FAB’を計算する。
上記(8)と同様の方法で座標変換を行い、座標系Bによるターゲット23の中心座標(xtB’,ytB’,ztB’)を計算する。この座標は硬性内視鏡10の光軸が通る座標である。
上記(9)で計算された硬性内視鏡10の先端座標(xsB,ysB,zsB)から上記(14)で計算されたターゲット23の中心座標(xtB’,ytB’,ztB’)に向かうベクトルの成分(xtB’−xsB,ytB’−ysB,ztB’−zsB)を計算し、硬性内視鏡10の先端座標(xsB,ysB,zsB)と共に記憶する。
これにより硬性内視鏡10の光軸位置が光軸の起点の座標(xsB,ysB,zsB)と光軸方向のベクトルの成分(xtB’−xsB,ytB’−ysB,ztB’−zsB)という形で取得されたことになる。
また、上記(12)で計算された硬性内視鏡10の複数の球体12の中心座標(x1B,y1B,z1B),(x2B,y2B,z2B),(x3B,y3B,z3B),(x4B,y4B,z4B)も記憶する。
これにより座標系Bで(即ち同一座標系で)硬性内視鏡10の光軸位置と硬性内視鏡10の位置および姿勢を定義する手段である複数の球体12の中心座標とを取得したことになる。
これらの座標およびベクトル成分を記憶しておけば、手術の期間中の3次元形状測定から検出した硬性内視鏡10の複数の球体12の中心座標と記憶している複数の球体12の中心座標とから座標変換関数を算出し、算出した座標変換関数を用いて記憶している硬性内視鏡10の起点の座標と光軸方向のベクトルの成分を座標変換すれば、同一座標系で人体の3次元形状データと硬性内視鏡10の光軸位置とを得ることができる。
上記実施形態における第1の較正用物体20,第2の較正用物体30は様々な変形が可能である。
上記実施形態では第1の較正用物体20として3つの大きさが異なる球体が取り付けられた直方体を有する物体を用いたが、第1の較正用物体20を以下の1)〜3)のようにし、第1の較正用物体20における座標系A,座標系Bでの定点およびベクトル検出と座標変換係数FABの計算を以下のようにすることもできる。なお、いずれの場合も座標系Aによるターゲット23の中心座標の取得方法は上記実施形態と同一である。
上記実施形態の(2)と同じ方法により2つの球体の中心座標を座標系A,座標系Bで取得し、上記実施形態の(3)と同じ方法により直方体前面の法線ベクトルを座標系A,座標系Bで取得する。そして座標変換係数FABを2つの球体の中心座標と1つの法線ベクトルを用いて計算する。計算は以下のように行う。
座標変換係数を計算するとは、上記実施形態の(7)で示したように、回転成分Mのg11〜g33と移動成分α,β,γを計算することである。座標系Aでの2つの球体の中心座標を(xaA,yaA,zaA),(xbA,ybA,zbA)、法線ベクトルを(aA,bA,cA)、座標系Bでの2つの球体の中心座標を(xaB,yaB,zaB),(xbB,ybB,zbB)、法線ベクトルを(aB,bB,cB)とし、(xaA,yaA,zaA)と(xaB,yaB,zaB)および(xbA,ybA,zbA)と(xbB,ybB,zbB)を上記実施形態の(7)で示した式に代入し、(aA,bA,cA)と(aB,bB,cB)を上記実施形態の(7)で示した式の移動成分(α,β,γ)を除いた式へ代入すると9つの式が成立する。
次に(xaB,yaB,zaB)から(xbB,ybB,zbB)へ向かうベクトルと法線ベクトル(aB,bB,cB)の外積によるベクトルの成分(dB,eB,fB)を計算し、(xaA,yaA,zaA)から(xbA,ybA,zbA)へ向かうベクトルと法線ベクトル(aA,bA,cA)の外積によるベクトルの成分(dA,eA,fA)を計算する。そして上記実施形態の(7)で示した式の移動成分(α,β,γ)を除いた式へ代入すると3つの式が成立し、合わせて12の式が成立する。この12の式からg11〜g33が同じものを集めて連立方程式を作成し、この連立方程式を解くことによりg11〜g33およびα,β,γを計算することができる。詳細は、特開2005−249402号公報の段落0104〜0118を参照する。
上記実施形態の(2)と同じ方法により1つの球体22Aの中心座標を座標系A,座標系Bで取得する。そして、直方体21の前面と直方体21に隣接する回転平板24の表面の2つの法線ベクトルのベクトル成分の検出を以下の方法で行う。
・座標系A(第1の較正用物体20に球体26を取り付けての3次元測定)
上記実施形態の(3)と同じ方法により直方体21の前面の平面式および法線ベクトルを算出する。
次に直方体21の前面以外の点群データを用いて上記実施形態の(3)に書かれた平面の点群データ抽出と平面式の計算方法によりいずれかの平面の平面式を算出する。この平面と球体26の中心座標との距離および平面の点群データの数から平面が直方体21の側面か上面か直方体に隣接する回転平板24の表面か光軸位置測定装置1の土台の表面かを判定し、直方体21に隣接する回転平板24の表面の平面式が算出されるまで上記処理を繰り返す。平面式の係数a,b,cから法線ベクトル(a,b,c),(−a,−b,−c)が出るので、回転平板24の表面から回転平板24の外側に向かう法線ベクトルの選定を上記実施形態の(3)と同じ方法により行う。
・座標系B(光軸位置測定装置1に硬性内視鏡10を固定しての3次元形状測定)
上記実施形態の(3)に書かれた平面の点群データ抽出と平面式の計算方法によりいずれかの平面の平面式を算出する。この平面の法線ベクトルと硬性内視鏡10の複数の球体12の1つ球体における中心座標から球体22Aの中心座標へ向かうベクトルとがなす角度により平面式が直方体21の前面の平面か否かを判定する。この処理を直方体21の前面の平面式と判定されるまで行う。
次に球体22Aを球体26の代わりに用いて上記の座標系Aの場合と同じ方法により回転平板24の表面の平面式を算出する。直方体21の前面および回転平板24の表面とも平面式の係数a,b,cから法線ベクトル(a,b,c),(−a,−b,−c)が出るので、回転平板24の表面から回転平板24の外側に向かう法線ベクトルの選定を上記実施形態の(3)と同じ方法により行う。この場合は球体26はないので、3次元形状測定装置の原点座標から球体22Aの中心座標に向かうベクトルをベクトルAとする。
座標系A,座標系Bでの球体22Aの中心座標と2つの平面の法線ベクトルとから座標変換係数を計算する。計算は以下のように行う。
座標変換係数を計算するとは、上記実施形態の(7)で示したように、回転成分Mのg11〜g33と移動成分α,β,γを計算することである。座標系Aでの球体の中心座標を(xaA,yaA,zaA),2つの法線ベクトルを(a1A,b1A,c1A),(a2A,b2A,c2A)とし、座標系Bでの球体の中心座標を(xaB,yaB,zaB),2つの法線ベクトルを(a1B,b1B,c1B),(a2B,b2B,c2B)とし、これらの座標を上記実施形態の(7)で示した式に代入し、ベクトル成分を上記実施形態の(7)で示した式から移動成分(α,β,γ)を除いた式に代入すると9つの式が成立する。
次に(a1A,b1A,c1A)と(a2A,b2A,c2A)との外積によるベクトルの成分(dA,eA,fA)を計算し、(a1B,b1B,c1B)と(a2B,b2B,c2B)との外積によるベクトルの成分(dB,eB,fB)を計算する。そして上記実施形態の(7)で示した式の移動成分(α,β,γ)を除いた式のx’,y’,z’に(dB,eB,fB)を代入し、x,y,zに(dA,eA,fA)を代入した式により3つの式が成立し、合わせて12の式が成立する。この12の式からg11〜g33が同じものを集めて連立方程式を作成し、この連立方程式を解くことによりg11〜g33およびα,β,γを計算することができる。詳細は、特開2005−249402号公報の式(数19)における座標の座標変換の1つをベクトルの座標変換に変えるのみで段落0104〜0118に示される計算の仕方と同じである。
直方体21の前面と側面と上面でできる角の座標の算出と直方体21の前面,側面,上面または直方体21に隣接する回転平板24の表面の2つ以上の法線ベクトルのベクトル成分の算出を以下の方法で行う。
上記実施形態の(3)に示された平面の点群データ抽出と平面式の計算方法により、複数の平面の平面式を可能な限り算出する。算出される平面式は直方体21の前面、側面、上面、直方体21に隣接する回転平板24の表面、光軸位置測定装置1の土台の表面の平面式5つである。
法線ベクトルのなす角度を内積の式から算出し、互いに90度近い角度になる3つの平面の平面式でなる連立方程式を解いて3つの平面式でできる角の座標を計算する。計算される角の座標は直方体21の前面と側面と上面でできる角と、直方体21の前面と上面と回転平板24の表面でできる角と、直方体21の前面と側面と光軸位置測定装置1の土台の表面でできる角と、直方体21の前面と回転平板24の表面と光軸位置測定装置1の土台の表面でできる角の4つである。
次に法線ベクトルのなす角度を内積の式から算出し、互いに90度近い角度になる2つの平面の平面式で直線の式を算出する。これにより8つの直線の式が計算される。
直線の式に点群データを代入して、直線の近傍にある点群データを抽出する。そして抽出した点群データの内、直線内にある角の座標からの距離が最も大きいものを抽出する。
直方体21の3辺の長さと回転平板24の横方向直線部分の長さおよび光軸位置測定装置1の土台から直方体21までの距離が既知であれば、この長さと距離から8つの直線の中から直方体21の辺に相当する直線を判別できる。8つの直線の式すべてにおいてこの距離を計算して直方体21の3辺の直線の式を判別し、この直線の式を算出した平面式がいずれであるかより直方体21の前面と側面と上面の3つの平面式を判別する。
3つの平面式による連立方程式を解くことで直方体21の前面と側面と上面でできる角の座標が取得され、平面式から直方体21の前面,側面,上面,回転平板24の表面の2つ以上の法線ベクトルのベクトル成分が算出される。
なお座標系Aの場合は球体26と平面との距離を平面式を判別する条件に入れてもよい。
平面式の係数a,b,cから出る法線ベクトル(a,b,c),(−a,−b,−c)の内、直方体21および回転平板24の平面から直方体21および回転平板24の外側に向かう法線ベクトルの選定は上記の2)球体を1つ取り付ける場合と同じである。なお座標系Bの場合は球体26の中心座標の代わりに直方体21
の前面と側面と上面でできる角の座標を使う。座標変換関数の計算は上記の2)球体を1つ取り付ける場合と同じである。
この変形例においては、角の座標および2つ以上のベクトル成分の取得までのデータ処理に時間を要するため、直方体21の反射率を回転平板24および光軸位置測定装置1の土台の反射率と大きく異ならせたり、直方体21のそれぞれの平面の反射率を異ならせ、3次元形状測定の際、同時に反射光量を測定して反射光量により点群データを分類すれば、データ処理の時間を大幅に短縮することができる。
上記実施形態では(1)〜(15)の手順により、第1の較正用物体20の定点座標とターゲット23の中心座標とを座標系Aにより求め、硬性内視鏡10を固定したときの第1の較正用物体10の定点座標を座標系Bにより求め、座標系Aによる座標を座標系Bによる座標に座標変換する座標変換関数FABを求め、座標系Bによるターゲット23の中心座標を、座標系Aによるターゲット23の中心座標を座標系Bの座標に座標変換することで求めたが、第1の較正用物体20を以下のようにすれば、第1の較正用物体20に定義される定点座標を用いてターゲット23の中心座標を求めることができる。すなわち、上記実施形態の(1)〜(4)の作業が不要になり、最初から光軸位置測定装置1に硬性内視鏡10を固定し、硬性内視鏡10の光軸位置を検出することができる。これによれば、第1の較正用物体20の加工を精度よく行う必要があるが、硬性内視鏡10の光軸位置と硬性内視鏡10の位置姿勢を定義する手段である複数の球体12の中心座標との関係を同一座標系で検出する作業が簡単になるという効果がある。
この場合はターゲット23の中心座標を次のように求める。光軸位置測定装置1に硬性内視鏡10を固定して3次元形状測定し、取得した点群データを用いて上記第1の較正用物体20の変形例1の3)に記載された平面式の算出と平面式の判別から、直方体21の前面、側面、上面と回転平板24の表面を判別し、直方体21の前面、側面、上面でなる角P1と直方体21の前面、上面、回転平板24の表面でなる角P2の座標を求める、次に既知である直方体21の3辺の長さと直方体21の前面、側面、上面の判別の過程で求めた直線の式から、角P3,P4の座標を求める。
次に角P1と角P4、角P2と角P3を結ぶ直線の式を計算し、直線の式からなる連立方程式を解いて直線の交点座標を計算する。計算された交点座標がターゲット23の中心座標になる。
なお、この場合も上記第1の較正用物体20の変形例1の3)と同様、直方体21の反射率を回転平板24および光軸位置測定装置1の土台の反射率と大きく異ならせたり、直方体21のそれぞれの平面の反射率を異ならせ、3次元形状測定の際、同時に反射光量を測定して反射光量により点群データを分類すれば、データ処理の時間を大幅に短縮することができる。
この場合はターゲット43の中心座標を次のように求める。光軸位置測定装置1に硬性内視鏡10を固定して3次元形状測定し、取得した点群データを用いて上記実施形態の(3)に書かれた平面の点群データ抽出と平面式の計算方法によりいずれかの平面の平面式を算出する。この平面を法線ベクトル方向に所定距離移動した式を算出し、その平面から所定距離範囲内にある点群データを抽出して、抽出した点群データが円の方程式に当てはまるか否かを判定する。当てはまれば平面式は円柱41の上面のものであると判定し、当てはまらなければ別の平面と判定する。この処理を円柱41の上面の平面式が算出されるまで行う。
次に円柱41の上面である平面をその法線ベクトル方向に所定間隔ごとに移動したときの平面式を複数算出し、それぞれの平面において平面から所定距離範囲内にある点群データを抽出して、抽出した点群データから最小2乗法により円の方程式を算出し、円の中心座標を算出する。
次に複数の円の中心座標から最小2乗法により直線の方程式を算出し、直線の方程式と円柱41の上面の平面式からなる連立方程式を解いて円柱41の中心線と円柱41の上面が交わる点(即ち円である上面の中心点)の座標を算出する。
算出された交点座標がターゲット43の中心座標になる。
この場合はターゲット53の中心座標を次のように求める。光軸位置測定装置1に硬性内視鏡10を固定して3次元形状測定し、取得した点群データを用いて上記実施形態の(3)に書かれた平面の点群データ抽出と平面式の計算方法によりいずれかの平面の平面式を算出する。この平面を法線ベクトル方向に所定距離移動した式を算出し、その平面から所定距離範囲内にある点群データを抽出して、抽出した点群データが円の方程式に当てはまるか否かを判定する。当てはまれば平面式は直方体51の前面のものであると判定し、当てはまらなければそれ以外の平面と判定する。この処理を直方体51の前面の平面式が算出されるまで行う。
次に直方体51の前面である平面をその法線ベクトル方向に所定間隔ごとに移動したときの平面式を複数算出し、それぞれの平面において平面から所定距離範囲内にある点群データを抽出して、抽出した点群データから最小2乗法により円の方程式を算出し、円の半径と円の中心座標を算出する。次に複数の円の中心座標から最小2乗法により直線の方程式を算出し、直線の方程式と直方体51の前面の平面式からなる連立方程式を解いて円錐52の中心線と直方体51の前面が交わる点の座標Cを算出する。そして直方体51の前面からの所定間隔ごとにおけるそれぞれの円の半径が減少する割合から半径が0になる(即ち、円錐52の頂点に到達する)直方体51の前面からの距離Lを算出し、直方体51の前面の法線ベクトルで大きさがLのベクトルDの成分を計算する。
交点座標CにベクトルDの成分を加算した座標(即ち、円錐52の頂点座標)がターゲット53の中心座標になる。
上記実施形態では第2の較正用物体30は全体を1つの球体にし、この球体の中心座標と第1の較正用物体20のターゲット23の座標と予め測定されている球体の中心座標から硬性内視鏡10の先端までの距離により、硬性内視鏡10の先端座標(光軸の起点座標)を計算したが、第2の較正用物体30を以下のようにすれば、最初から硬性内視鏡10を図1に示される正規の位置にセットして、第2の較正用物体30に定義される定点座標およびベクトルから、硬性内視鏡10の先端座標(光軸の起点座標)を取得することができる。すなわち、上記実施形態の(5)〜(9)の作業が不要になり、硬性内視鏡10の光軸位置と硬性内視鏡10の位置姿勢を定義する手段である複数の球体12の中心座標との関係を同一座標系で検出する作業が簡単になるという効果がある。
この場合は硬性内視鏡10の先端座標(光軸の起点座標)を次のように求める。硬性内視鏡10を2つの球体61,62の穴に挿入させて固定し、3次元形状測定して取得した点群データから上記実施形態の(2)と同じ方法により2つの球体61,62の中心座標を取得する。そして、球体62の中心座標から球体61の中心座標に向かうベクトルの成分を計算し、このベクトル成分と球体61の中心座標と予め測定されている球体61の中心座標から硬性内視鏡10の先端までの距離Lにより、上記実施形態(9)に示されている計算方法で、硬性内視鏡10の先端座標を計算する。
この場合は硬性内視鏡10の先端座標(光軸の起点座標)を次のように求める。硬性内視鏡10を円柱70の穴に挿入させて固定し、3次元形状測定して取得した点群データから上記第1の較正用物体20の変形例2の(2)に記載された方法により円柱70の中心線の式と円柱70の上面の中心座標を取得する。
次に、円柱70の中心線の式から得られるこの直線方向のベクトルの成分と円柱70の上面の中心座標と、予め測定されている円柱70の上面から硬性内視鏡10の先端までの距離Lにより、上記実施形態(9)に示されている計算方法で、硬性内視鏡10の先端座標を計算する。
この場合は硬性内視鏡10の先端座標(光軸の起点座標)を次のように求める。硬性内視鏡10を直方体80の穴に挿入させて固定し、3次元形状測定して取得した点群データから上記第1の較正用物体20の変形例1の3)に記載された方法により直方体80前面の法線ベクトルの成分、前面と2つの側面とからなる角の座標、前面と2つの側面とからなる3辺のベクトル成分(直線の式)を計算する。
次に上記第1の較正用物体20の変形例2の1)に記載された方法により直方体80前面の中心座標を算出し、この中心座標と直方体前面の法線ベクトルの成分と、予め測定されている直方体80の前面から硬性内視鏡10の先端までの距離Lにより、上記実施形態(9)に示されている計算方法で、硬性内視鏡10の先端座標を計算する。
上記実施形態および第2の較正用物体30の変形例2では、硬性内視鏡10の長軸部11は直線であることを前提にしていた。しかし、硬性内視鏡10の長軸部11が直線でない場合は、上記実施形態および第2の較正用物体30の変形例1による方法では、長軸部11の先端座標(光軸の起点座標)を得ることはできない。
この場合は第2の較正用物体30を以下のようにすれば、長軸部11が直線でない場合でも硬性内視鏡10の先端座標(光軸の起点座標)を取得することができる。
この場合は硬性内視鏡10の先端座標(光軸の起点座標)を次のように求める。硬性内視鏡10を取り付ける前に上記実施形態1)〜4)と同じ方法により、座標系Cによる球体の中心座標と平板91の第1の較正用物体20側の表面における穴の中心座標とを取得する。穴の中心座標を取得するには図2のピン付け球体の代わりに穴の径に合う円柱が付けられた球体を穴に挿入すればよい。
次に硬性内視鏡10の先端を平板91の穴に挿入して光軸位置測定装置1に取り付け、3次元形状測定を行う。このとき先端位置が平板91の第1の較正用物体20側の表面に合うようにする。そして取得した点群データから上記実施形態(2)で示されている方法により座標系Bによる球体の中心座標を取得し、上記実施形態(7)で示されている方法により座標変換関数FCBを計算する。得られた座標変換関数FCBにより座標系Cによる平板91の第1の較正用物体20側の表面における穴の中心座標を座標系Bによる座標に変換する。これにより、硬性内視鏡10の先端座標(光軸の起点座標)と硬性内視鏡10の位置姿勢を定義する手段である複数の球体12の中心座標とを同じ座標系Bで取得することができる。
この場合の第2の較正用物体90は第1の較正用物体20の変形例1の1)〜3)と同じ変形が可能である。
この場合は硬性内視鏡10の先端座標(光軸の起点座標)の取得方法は上記第2の較正用物体30の変形例1の2)および3)と同一である。
上記実施形態および第2の較正用物体30の変形例1および変形例2では、硬性内視鏡10の先端の座標(厳密には、先端の断面と長軸部11の中心線が交差する点の座標)を光軸の起点の座標とした。しかしながら、硬性内視鏡10によっては図13に示すように先端の座標と光軸の起点座標とが一致していない場合がある。
この場合は、第2の較正用物体30を以下のようにすれば光軸位置および光軸の起点座標を取得することができる。
この場合は硬性内視鏡10の先端付近の光軸上の座標は、上記第2の較正用物体30の変形例2の1)と同じ方法により取得することができる。なお、ピン付の球体をターゲット113の中心に取り付けて3次元形状測定する場合、薄厚の透明板はそのままでは点群データは得られないので、表面に反射率の高いシートを貼るようにする。
そして、第1の較正用物体20のターゲット23の中心座標が別に取得されるので、これと硬性内視鏡10の先端付近の光軸上の座標とから光軸位置(光軸の直線の式)が取得される。
次に硬性内視鏡10の先端断面に合うように平板を当て、同じ測定位置で(即ち座標系Bで)3次元形状測定を行い、上記実施形態(3)に示された方法で平面の式を得る。この平面の式と先に得られた光軸の直線の式との連立方程式を解くことで、平面と光軸が交差する点である光軸の起点座標を取得することができる。
この場合の第2の較正用物体110は第1の較正用物体20の変形例1の1)〜3)と同じ変形が可能である。
この場合は硬性内視鏡10の先端付近の光軸上の座標は、上記第1の較正用物体20の変形例2の2)と同じ方法により取得することができる。これ以外は上記第2の較正用物体30の変形例3の1)と同じである。
なお、空洞の円柱121は空洞の直方体にすることもできる。この場合は硬性内視鏡10の先端付近の光軸上の座標は、上記第1の較正用物体20の変形例2の1)と同じ方法により取得することができる。
上記実施形態では第1の較正用物体20は直方体に球体を取り付けた構造であったが、定点または定点およびベクトルが定義できる形状を有し、ターゲット23を形成することができれば、色々な構造のものが考えられる。
以上、本発明の実施形態の一例を説明したが、本発明はこれに限定されるものではなく、特許請求の範囲に記載された技術的思想の範疇において各種の変更が可能であることは言うまでもない。
Claims (12)
- 3次元形状測定装置と、
長軸部と、長軸部先端からの光軸と、前記3次元形状測定装置により位置および姿勢が測定可能な第1標識部とを有する被測定物体と、
前記被測定物体を固定する固定手段と、
前記被測定物体の長軸部先端からの光軸が当たるターゲットと、前記3次元形状測定装置により位置および立体形状が測定可能な第2標識部とを有し、前記ターゲットの中心座標と前記第2標識部との3次元相対関係が予めわかっているか、または前記3次元形状測定装置により前記ターゲットの中心座標と前記第2標識部との3次元相対関係が測定可能な第1の較正用物体と、
前記3次元形状測定装置により位置および立体形状が測定可能な第3標識部を有し、前記固定手段により前記被測定物体が固定された状態において、前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標と前記第3標識部との3次元相対関係が予めわかっているか、または前記3次元形状測定装置により前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸上に含まれる座標と前記第3標識部との3次元相対関係が測定可能な第2の較正用物体と、
を有する3次元形状測定システムにおける、前記被測定物体の長軸部先端からの光軸と前記第1標識部との3次元相対関係を測定する3次元相対関係測定方法であって、
前記被測定物体を前記固定手段に固定し、前記第1の較正用物体のターゲットの中心に前記被測定物体の長軸部先端からの光軸が当たるようにするステップと、
前記3次元形状測定装置により、前記被測定物体、前記第1の較正用物体および前記第2の較正用物体の立体形状データ群を取得するステップと、
前記立体形状データ群から、前記第1標識部の位置および姿勢と前記第2標識部の位置および立体形状と前記第3標識部の位置および立体形状とを算出するステップと、
前記第2標識部の位置および立体形状と、前記第3標識部の位置および立体形状と、前記ターゲットの中心座標と前記第2標識部との3次元相対関係と、前記長軸部先端の座標または前記長軸部先端からの光軸上に含まれる座標と前記第3標識部との3次元相対関係とに基づいて、前記長軸部先端からの光軸の位置を算出するステップと、
前記第1標識部の位置および姿勢と、前記長軸部先端からの光軸の位置とに基づいて、前記被定物体の長軸部先端からの光軸と前記第1標識部との3次元相対関係を算出するステップと、
を有する、3次元相対関係測定方法。 - 前記第1の較正用物体の第2標識部は、予め立体形状がわかっており、
前記ターゲットの中心座標と前記第2標識部との3次元相対関係は、前記第2標識部の立体形状に関連するベクトルまたは定点位置と、前記ターゲットの中心座標との関係により定義される、請求項1記載の3次元相対関係測定方法。 - 前記第1の較正用物体のターゲットの中心座標と前記第2標識部の定点位置とは等しいか、または前記第1の較正用物体のターゲットの中心座標は前記第2標識部における2つ以上の定点を用いて定義される点であることを特徴とする請求項2記載の3次元相対関係測定方法。
- 前記第2の較正用物体の第3標識部は、予め立体形状がわかっており、
前記長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標と前記第3標識部との3次元相対関係は、前記第3標識部の立体形状に関連するベクトルまたは定点位置と、前記長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標との関係により定義される、請求項1乃至3いずれか記載の3次元相対関係測定方法。 - 前記長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標と前記第3標識部の定点位置とは等しいか、あるいは前記長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標は前記第3標識部における2つ以上の定点を用いて定義される点であることを特徴とする請求項4記載の3次元相対関係測定方法。
- 前記長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標は、前記第3標識部における2つの定点を結んだ直線または前記第3標識部における定点とベクトルから定まる直線上において、定点から前記第3標識部の立体形状に関連する長さにより定まる点であることを特徴とする請求項5記載の3次元相対関係測定方法。
- 前記長軸部先端からの光軸に含まれる座標は、光を透過する平板に形成されたターゲットの中心座標として定義されることを特徴とする請求項1乃至6いずれか記載の3次元相対関係測定方法。
- 前記第1の較正用物体および前記第2の較正用物体は、球体、多面体、円柱、円錐のいずれかの形状を一部または全体に有している、請求項1乃至7いずれか記載の3次元相対関係測定方法。
- 前記被測定物体は体内挿入器具であって、前記長軸部は体内挿入部であり、前記長軸部先端からの光軸は長軸部先端からの撮像方向であり、前記第1標識部は体内に挿入されない部分に設けられている、請求項1乃至8いずれか記載の3次元相対関係測定方法。
- 3次元形状測定装置と、
長軸部と、長軸部先端からの光軸と、前記3次元形状測定装置により位置および姿勢が測定可能な第1標識部とを有する被測定物体と、
前記被測定物体を固定する固定手段と、
前記被測定物体の長軸部先端からの光軸が当たるターゲットと、前記3次元形状測定装置により位置および立体形状が測定可能な第2標識部とを有し、前記ターゲットの中心座標と前記第2標識部との3次元相対関係が予めわかっているか、または前記3次元形状測定装置により前記ターゲットの中心座標と前記第2標識部との3次元相対関係が測定可能な第1の較正用物体と、
前記3次元形状測定装置により位置および立体形状が測定可能な第3標識部を有し、前記固定手段により前記被測定物体が固定された状態において、前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標と前記第3標識部との3次元相対関係が予めわかっているか、または前記3次元形状測定装置により前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸上に含まれる座標と前記第3標識部との3次元相対関係が測定可能な第2の較正用物体と、
を有する3次元形状測定システムにおける、前記被測定物体の長軸部先端からの光軸と前記第1標識部との3次元相対関係を測定する3次元相対関係測定プログラムであって、
前記被測定物体を前記固定手段に固定し、前記第1の較正用物体のターゲットの中心に前記被測定物体の長軸部先端からの光軸を一致させた状態で、前記3次元形状測定装置により、前記被測定物体、前記第1の較正用物体および前記第2の較正用物体の立体形状データ群を取得するステップと、
前記立体形状データ群から、前記第1標識部の位置および姿勢と前記第2標識部の位置および立体形状と前記第3標識部の位置および立体形状とを算出するステップと、
前記第2標識部の位置および立体形状と、前記第3標識部の位置および立体形状と、前記ターゲットの中心座標と前記第2標識部との3次元相対関係と、前記長軸部先端の座標または前記長軸部先端からの光軸上に含まれる座標と前記第3標識部との3次元相対関係とに基づいて、前記長軸部先端からの光軸の位置を算出するステップと、
前記第1標識部の位置および姿勢と、前記長軸部先端からの光軸の位置とに基づいて、前記被定物体の長軸部先端からの光軸と前記第1標識部との3次元相対関係を算出するステップと、
を有する、3次元相対関係測定プログラム。 - 3次元形状測定装置と、
長軸部と、長軸部先端からの光軸と、前記3次元形状測定装置により位置および姿勢が測定可能な第1標識部とを有する被測定物体と、
前記被測定物体を固定する固定手段と、
前記被測定物体の長軸部先端からの光軸が当たるターゲットと、前記3次元形状測定装置により位置および立体形状が測定可能な第2標識部とを有し、前記ターゲットの中心座標と前記第2標識部との3次元相対関係が予めわかっているか、または前記3次元形状測定装置により前記ターゲットの中心座標と前記第2標識部との3次元相対関係が測定可能な第1の較正用物体と、
前記3次元形状測定装置により位置および立体形状が測定可能な第3標識部を有し、前記固定手段により前記被測定物体が固定された状態において、前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標と前記第3標識部との3次元相対関係が予めわかっているか、または前記3次元形状測定装置により前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸上に含まれる座標と前記第3標識部との3次元相対関係が測定可能な第2の較正用物体と、を有し、
前記被測定物体を前記固定手段に固定し、前記第1の較正用物体のターゲットの中心に前記被測定物体の長軸部先端からの光軸を一致させた状態で、前記3次元形状測定装置により、前記被測定物体、前記第1の較正用物体および前記第2の較正用物体の立体形状データ群を取得する手段と、
前記立体形状データ群から、前記第1標識部の位置および姿勢と前記第2標識部の位置および立体形状と前記第3標識部の位置および立体形状とを算出する手段と、
前記第2標識部の位置および立体形状と、前記第3標識部の位置および立体形状と、前記ターゲットの中心座標と前記第2標識部との3次元相対関係と、前記長軸部先端の座標または前記長軸部先端からの光軸上に含まれる座標と前記第3標識部との3次元相対関係とに基づいて、前記長軸部先端からの光軸の位置を算出する手段と、
前記第1標識部の位置および姿勢と、前記長軸部先端からの光軸の位置とに基づいて、前記被定物体の長軸部先端からの光軸と前記第1標識部との3次元相対関係を算出する手段と、
を有する、3次元相対関係測定システム。 - 長軸部と、長軸部先端からの光軸と、3次元形状測定装置により位置および姿勢が測定可能な第1標識部とを有する被測定物体の、前記光軸と前記第1標識部との3次元相対関係を測定するための3次元相対関係測定装置であって、
前記被測定物体を固定可能な固定手段と、
前記被測定物体の長軸部先端からの光軸が当たるターゲットと、3次元形状測定装置により位置および立体形状が測定可能な第2標識部とを有し、前記ターゲットの中心座標と前記第2標識部との3次元相対関係が予めわかっているか、または3次元形状測定装置により前記ターゲットの中心座標と前記第2標識部との3次元相対関係が測定可能な第1の較正用物体と、
3次元形状測定装置により位置および立体形状が測定可能な第3標識部を有し、前記固定手段により前記被測定物体が固定された状態において、前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸に含まれる座標と前記第3標識部との3次元相対関係が予めわかっているか、または3次元形状測定装置により前記被測定物体の長軸部先端の座標または前記長軸部先端からの光軸上に含まれる座標と前記第3標識部との3次元相対関係が測定可能な第2の較正用物体と、
前記第1の較正用物体を前記被測定物体と相対的に移動させる移動手段と、を有し、
前記移動手段により前記第1の較正用物体を移動させて前記ターゲットに前記長軸部先端からの光軸が当たるようにした後、3次元形状測定装置で前記第1乃至第3標識部の3次元相対位置を測定することにより、前記被測定物体の前記光軸と前記第1標識部との3次元相対関係を測定可能な3次元相対関係測定装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009029216A JP5560424B2 (ja) | 2009-02-11 | 2009-02-11 | 長軸部先端からの光軸を有する物体における光軸位置と該物体の位置姿勢を定義する手段との3次元相対関係測定方法、装置およびシステム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009029216A JP5560424B2 (ja) | 2009-02-11 | 2009-02-11 | 長軸部先端からの光軸を有する物体における光軸位置と該物体の位置姿勢を定義する手段との3次元相対関係測定方法、装置およびシステム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2010183995A true JP2010183995A (ja) | 2010-08-26 |
JP5560424B2 JP5560424B2 (ja) | 2014-07-30 |
Family
ID=42765011
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009029216A Expired - Fee Related JP5560424B2 (ja) | 2009-02-11 | 2009-02-11 | 長軸部先端からの光軸を有する物体における光軸位置と該物体の位置姿勢を定義する手段との3次元相対関係測定方法、装置およびシステム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5560424B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112085847A (zh) * | 2020-08-21 | 2020-12-15 | 深圳市瑞立视多媒体科技有限公司 | 一种刚体标记点优化方法、设备及光学动作捕捉系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11164839A (ja) * | 1997-09-26 | 1999-06-22 | Picker Internatl Inc | 顕微鏡較正装置及び方法 |
JP2000511287A (ja) * | 1997-03-20 | 2000-08-29 | ハートフォード ホスピタル | 内視鏡の性能特性を評価する方法および装置 |
JP2008104877A (ja) * | 2006-10-23 | 2008-05-08 | Pentax Corp | スコープナビゲーション装置、方法、およびコンピュータ可読な記憶媒体 |
JP2010046276A (ja) * | 2008-08-21 | 2010-03-04 | Fujifilm Corp | 内視鏡画像のズレ量測定装置及び方法、並びに電子内視鏡及び内視鏡用画像処理装置 |
WO2010101117A1 (ja) * | 2009-03-01 | 2010-09-10 | 国立大学法人浜松医科大学 | 手術支援システム |
-
2009
- 2009-02-11 JP JP2009029216A patent/JP5560424B2/ja not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000511287A (ja) * | 1997-03-20 | 2000-08-29 | ハートフォード ホスピタル | 内視鏡の性能特性を評価する方法および装置 |
JPH11164839A (ja) * | 1997-09-26 | 1999-06-22 | Picker Internatl Inc | 顕微鏡較正装置及び方法 |
JP2008104877A (ja) * | 2006-10-23 | 2008-05-08 | Pentax Corp | スコープナビゲーション装置、方法、およびコンピュータ可読な記憶媒体 |
JP2010046276A (ja) * | 2008-08-21 | 2010-03-04 | Fujifilm Corp | 内視鏡画像のズレ量測定装置及び方法、並びに電子内視鏡及び内視鏡用画像処理装置 |
WO2010101117A1 (ja) * | 2009-03-01 | 2010-09-10 | 国立大学法人浜松医科大学 | 手術支援システム |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112085847A (zh) * | 2020-08-21 | 2020-12-15 | 深圳市瑞立视多媒体科技有限公司 | 一种刚体标记点优化方法、设备及光学动作捕捉系统 |
CN112085847B (zh) * | 2020-08-21 | 2024-05-17 | 深圳市瑞立视多媒体科技有限公司 | 一种刚体标记点优化方法、设备及光学动作捕捉系统 |
Also Published As
Publication number | Publication date |
---|---|
JP5560424B2 (ja) | 2014-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5137033B2 (ja) | 手術支援情報表示装置、手術支援情報表示方法及び手術支援情報表示プログラム | |
EP2953569B1 (en) | Tracking apparatus for tracking an object with respect to a body | |
JP6129344B2 (ja) | 内視鏡、殊に低侵襲手術のための内視鏡 | |
EP2837328B1 (en) | Unmapped region visualization | |
JP5561458B2 (ja) | 手術支援システム | |
JP5283157B2 (ja) | 手術支援情報表示装置、手術支援情報表示方法及び手術支援情報表示プログラム | |
US20200129240A1 (en) | Systems and methods for intraoperative planning and placement of implants | |
JP5569711B2 (ja) | 手術支援システム | |
CN112184653B (zh) | 一种基于双目内窥镜的病灶三维尺寸测量及显示方法 | |
CN116327079A (zh) | 内窥镜测量系统和工具 | |
WO2017170488A1 (ja) | 光軸位置測定システム、光軸位置測定方法、光軸位置測定プログラム、光軸位置測定装置 | |
KR100996826B1 (ko) | 수술용 항법 장치용 도구의 교정 방법 | |
JP5213201B2 (ja) | 体内挿入器具の種類を識別可能な手術支援システム | |
JP5560424B2 (ja) | 長軸部先端からの光軸を有する物体における光軸位置と該物体の位置姿勢を定義する手段との3次元相対関係測定方法、装置およびシステム | |
WO2009107703A1 (ja) | 体内挿入器具の種類を識別可能な手術支援システム | |
JP2017185212A (ja) | 光軸位置測定システム、光軸位置測定方法、光軸位置測定プログラム、光軸位置測定装置 | |
JP5200582B2 (ja) | 長軸部を有する物体の長軸部の先端座標と該物体の位置姿勢を定義する手段との3次元相対関係測定方法およびシステム | |
JP2016036594A (ja) | 医用装置及び超音波診断装置 | |
US20230368418A1 (en) | Accuracy check and automatic calibration of tracked instruments | |
US20230248441A1 (en) | Extended-reality visualization of endovascular navigation | |
US20100198110A1 (en) | Method and device for determining a characteristic property of an anatomical structure | |
CN204394569U (zh) | 一种带三维电子罗盘的超声心动图探头 | |
CN118251188A (zh) | 具有改进的器械跟踪的外科手术导航系统和导航方法 | |
TW202110404A (zh) | 超音波影像系統 | |
JP2005111115A (ja) | 立体マーカ、これを用いた被験体の向き位置測定方法及び装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20120114 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120114 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20131105 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20131212 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20140430 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140510 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140510 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5560424 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |