JP2009045097A - 三次元画像生成装置及び三次元画像生成方法 - Google Patents
三次元画像生成装置及び三次元画像生成方法 Download PDFInfo
- Publication number
- JP2009045097A JP2009045097A JP2007211126A JP2007211126A JP2009045097A JP 2009045097 A JP2009045097 A JP 2009045097A JP 2007211126 A JP2007211126 A JP 2007211126A JP 2007211126 A JP2007211126 A JP 2007211126A JP 2009045097 A JP2009045097 A JP 2009045097A
- Authority
- JP
- Japan
- Prior art keywords
- probe
- dimensional
- image
- sensor
- magnetic field
- 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.)
- Pending
Links
Landscapes
- Measurement Of Length, Angles, Or The Like Using Electric Or Magnetic Means (AREA)
- Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
【課題】安全で測定精度の良い三次元画像を生成する三次元画像生成装置を提供する。
【解決手段】この三次元画像生成装置110は、被検査対象物6に超音波を送信する送信部と被検査対象物6から反射した超音波を受信する受信部とを備えたプローブ1と、プローブ1の三次元位置情報及び姿勢情報を生成するセンサ2と、所定のエリアに磁界ベクトルを形成するトランスミッタ4と、三次元位置情報及び姿勢情報に基づいてプローブ1から取得された反射画像の三次元形状を演算するPC(演算部)12と、演算結果を三次元画像として表示するモニタ(表示部)14と、センサ2からの座標をリアルタイムに演算する位置検出回路10と、プローブ1から得られた反射画像データを処理するイメージプロセッサ11と、を備えて構成されている。
【選択図】図2
【解決手段】この三次元画像生成装置110は、被検査対象物6に超音波を送信する送信部と被検査対象物6から反射した超音波を受信する受信部とを備えたプローブ1と、プローブ1の三次元位置情報及び姿勢情報を生成するセンサ2と、所定のエリアに磁界ベクトルを形成するトランスミッタ4と、三次元位置情報及び姿勢情報に基づいてプローブ1から取得された反射画像の三次元形状を演算するPC(演算部)12と、演算結果を三次元画像として表示するモニタ(表示部)14と、センサ2からの座標をリアルタイムに演算する位置検出回路10と、プローブ1から得られた反射画像データを処理するイメージプロセッサ11と、を備えて構成されている。
【選択図】図2
Description
本発明は、三次元画像生成装置に関し、さらに詳しくは、超音波プローブにより取得した被測定対象物の反射画像の各座標点を、プローブの位置情報と姿勢情報に基づいて、三次元座標を演算して三次元画像を生成する三次元画像生成装置とその生成方法に関するものである。
近年、メタボリックシンドロームと呼ばれる肥満と脂質代謝異常、高血圧、高血糖を伴う症候群が増加しており、我が国では予備群を含めると2000万人が該当すると言われている。我が国におけるメタボリックシンドローム診断基準は、2005年4月に策定され、腹部肥満、内臓脂肪蓄積を必須項目としている。必須項目としての腹部肥満については、BMI(体格指数)の基準がなく、ウエスト周囲径のみで判定しており、腹部CT検査で測定した内臓脂肪蓄積の基準値である100cm2に相当するウエスト周囲径が基準となっており、男性85cm以上、女性80cm以上を腹部肥満と判定している。
また、内臓脂肪の測定には、コンピューター断層撮影法(CT)が一般的に用いられているが、人体に多方向からX線を照射するために、放射線被曝の問題が懸念されている。また、CT法では、得られた断面画像の脂肪断面の長さから脂肪量を推定することから、視点が固定される上に、推定誤差が生じる可能性がある。一方、MRI法では、体内の水素原子核の核磁気共鳴を利用して組織を断面画像化するが、CT法と同様に断面画像から脂肪量を推定するために、測定精度が劣る場合がある。また、超音波パルスのエコーを検出して測定する場合は、皮下脂肪組織と筋肉組織との境界面において反射して戻ってくるエコーを検出することから、生体の表面にある皮下脂肪は測定できるが、より生体の深層にある内臓脂肪は他の方法により求めるのが一般的である。
特許文献1には、レーザ光源とカメラに三次元センサを取り付け、対象物に照射したレーザビームをカメラにより撮影して、その画像とセンサから得られた位置情報と姿勢情報に基づいて対象物の三次元画像を生成する三次元計測装置について開示されている。
また、特許文献2は、生体組織(脂肪組織と除脂肪組織)の数以上の複数の異なる周波数を有する超音波信号を、一表面から生体内に入射して他の表面において生体内を透過した信号を検出し、その検出時の強度比を求め、また、生体の一表面から他の表面までの距離を求め、超音波信号の入射時の強度比の対数値と出力時の強度比の対数値とに基づいて各生体組織の厚さを演算する超音波生体組織測定方法について開示されている。
また、非特許文献1は、三次元デジタイザにより計測断面の位置と向きを取得することで、超音波三次元画像を再構成する技術について報告されている。
特開2005−98978公報
特開2000−237189公報
[位置の同時計測による超音波2次元断層像の3次元再構成]田野島英司、眞渓歩、大城理、千原國宏、奈良先端科学技術大学院大学
また、特許文献2は、生体組織(脂肪組織と除脂肪組織)の数以上の複数の異なる周波数を有する超音波信号を、一表面から生体内に入射して他の表面において生体内を透過した信号を検出し、その検出時の強度比を求め、また、生体の一表面から他の表面までの距離を求め、超音波信号の入射時の強度比の対数値と出力時の強度比の対数値とに基づいて各生体組織の厚さを演算する超音波生体組織測定方法について開示されている。
また、非特許文献1は、三次元デジタイザにより計測断面の位置と向きを取得することで、超音波三次元画像を再構成する技術について報告されている。
しかし、内臓脂肪の測定法として腹部コンピューター断層撮影やMRI法が利用されているが、機器が大型で高価であるといった問題がある。また妊婦等での使用には注意が必要となる他、断面画像の脂肪断面の長さから脂肪量を推定するために、推定誤差が生じる可能性がある等、汎用性に問題がある。
また、特許文献1に開示されている従来技術は、三角測量の原理に基づくため、計測精度を上げるために、レーザ光源、カメラ及び計測対象の配置を広く取る必要がある。
また、特許文献2に開示されている従来技術は、生体の部位(例えば骨等)によっては、入射した超音波が反射して検出できないといった問題がある。
また、非特許文献1に開示されている従来技術は、三次元デジタイザにより検査対象物をなぞって三次元画像を生成するがアームに取り付けてあるため、アームが邪魔になることや、計測器の姿勢に制限があるため操作が煩わしいといった問題がある。
本発明は、かかる課題に鑑み、超音波プローブに三次元の位置情報と姿勢情報を検出する磁気センサを取り付け、超音波プローブで反射画像情報を取得すると共に、リアルタイムにプローブの三次元情報を取得することにより、安全で測定精度の良い三次元画像を生成する三次元画像生成装置を提供するものである。
また、特許文献1に開示されている従来技術は、三角測量の原理に基づくため、計測精度を上げるために、レーザ光源、カメラ及び計測対象の配置を広く取る必要がある。
また、特許文献2に開示されている従来技術は、生体の部位(例えば骨等)によっては、入射した超音波が反射して検出できないといった問題がある。
また、非特許文献1に開示されている従来技術は、三次元デジタイザにより検査対象物をなぞって三次元画像を生成するがアームに取り付けてあるため、アームが邪魔になることや、計測器の姿勢に制限があるため操作が煩わしいといった問題がある。
本発明は、かかる課題に鑑み、超音波プローブに三次元の位置情報と姿勢情報を検出する磁気センサを取り付け、超音波プローブで反射画像情報を取得すると共に、リアルタイムにプローブの三次元情報を取得することにより、安全で測定精度の良い三次元画像を生成する三次元画像生成装置を提供するものである。
本発明はかかる課題を解決するために、請求項1は、被検査対象物に超音波を送信する送信部及び該被検査対象物から反射した超音波を受信する受信部を備えたプローブと、該プローブの三次元位置情報及び姿勢情報を生成するセンサと、所定のエリアに磁界ベクトルを形成するトランスミッタと、前記三次元位置情報及び前記姿勢情報に基づいて前記プローブから取得された反射画像の三次元形状を演算する演算部と、該演算部による演算結果を三次元画像として表示する表示部と、を備えたことを特徴とする。
超音波プローブにより取得された画像は、超音波が被検査対象物内に進行して、固いものに当たって反射した反射音波を測定し、反射音波が返ってくるまでの時間から距離を計算して可視化したものである。そして、複数の超音波を発生することで二次元の反射画像を生成することができる。この二次元の反射画像は被検査対象物の1つの断層画像を表示しているに過ぎない。そこで本発明では、超音波プローブにセンサを取り付け、超音波プローブを移動したときの位置情報と、姿勢情報を取得し、移動することにより得られた位置情報と姿勢情報に基づいて二次元の断層画像を演算部により演算して三次元形状を得る。この結果を可視化して表示部に表示する。これにより、従来使用していた超音波プローブにセンサを取り付けるだけで容易に二次元の断層画像から三次元画像を生成することができる。
超音波プローブにより取得された画像は、超音波が被検査対象物内に進行して、固いものに当たって反射した反射音波を測定し、反射音波が返ってくるまでの時間から距離を計算して可視化したものである。そして、複数の超音波を発生することで二次元の反射画像を生成することができる。この二次元の反射画像は被検査対象物の1つの断層画像を表示しているに過ぎない。そこで本発明では、超音波プローブにセンサを取り付け、超音波プローブを移動したときの位置情報と、姿勢情報を取得し、移動することにより得られた位置情報と姿勢情報に基づいて二次元の断層画像を演算部により演算して三次元形状を得る。この結果を可視化して表示部に表示する。これにより、従来使用していた超音波プローブにセンサを取り付けるだけで容易に二次元の断層画像から三次元画像を生成することができる。
請求項2は、前記演算部は、前記プローブを移動させることで前記被検査対象物からの前記反射画像を取得すると共に、前記センサにより前記磁界ベクトルを読み取って前記トランスミッタの原点に対する当該プローブの三次元位置情報及び姿勢情報を検出し、前記反射画像の各座標点における前記三次元位置情報及び前記姿勢情報に基づいて前記各座標点の三次元座標をリアルタイムに演算することを特徴とする。
プローブは被検査対象物に接触することで反射画像を取得することができる。また、プローブを移動させることで、センサからプローブの三次元位置情報及び姿勢情報が検出される。演算部は反射画像の各座標点における三次元位置情報及び姿勢情報に基づいて、各座標点の三次元座標をリアルタイムに演算する。これにより、プローブを被検査対象物に対して移動させるだけで、反射画像の三次元形状をリアルタイムに演算することができる。
プローブは被検査対象物に接触することで反射画像を取得することができる。また、プローブを移動させることで、センサからプローブの三次元位置情報及び姿勢情報が検出される。演算部は反射画像の各座標点における三次元位置情報及び姿勢情報に基づいて、各座標点の三次元座標をリアルタイムに演算する。これにより、プローブを被検査対象物に対して移動させるだけで、反射画像の三次元形状をリアルタイムに演算することができる。
請求項3は、前記演算部は、前記センサにより受信した磁界ベクトルに基づいて、前記プローブの三次元位置情報、及び該プローブのロール角、ピッチ角、ヨー角に係る姿勢情報を演算することを特徴とする。
本発明では半球状の磁界を発生するトランスミッタを備え、その磁界の磁力線に対するベクトルを計算することにより三次元位置情報及び姿勢情報を取得するものである。従って、被検査対象物とプローブの距離はこの磁界のエリア内に限定される。これにより、磁界のエリア内であればプローブの姿勢がどのように変化しても正確な姿勢情報を得ることができる。
本発明では半球状の磁界を発生するトランスミッタを備え、その磁界の磁力線に対するベクトルを計算することにより三次元位置情報及び姿勢情報を取得するものである。従って、被検査対象物とプローブの距離はこの磁界のエリア内に限定される。これにより、磁界のエリア内であればプローブの姿勢がどのように変化しても正確な姿勢情報を得ることができる。
請求項4は、前記センサは磁気センサであり、前記プローブの適所に固定されていることを特徴とする。
1つの磁気センサではトランスミッタの位置における磁界の強さおよび方向をAとし、磁界の強さおよび方向AのX軸方向磁界の強さ、Y軸方向磁界の強さ、Z軸方向磁界の強さをそれぞれAx、Ay、Azとし、cosα、cosβ、cosγを磁界の強さおよび方向Aの方向余弦とすれば、Ax=Acosα、Ay=Acosβ、Az=Acosγであり、磁気センサのX軸方向の磁界検出コイルからはAxの出力が、Y軸方向の磁界検出コイルからはAyの出力が、Z軸方向の磁界検出コイルからはAzの出力が送出される。そして磁界の強さおよび方向Aは、A=(Ax2+Ay2+Az2)1/2で与えられる。従って、プローブの位置と姿勢を正確に取得するためには、センサとプローブが常に一体であることが必要である。これにより、単純な構成で各座標毎に独立した信号を同時に検出することができる。
1つの磁気センサではトランスミッタの位置における磁界の強さおよび方向をAとし、磁界の強さおよび方向AのX軸方向磁界の強さ、Y軸方向磁界の強さ、Z軸方向磁界の強さをそれぞれAx、Ay、Azとし、cosα、cosβ、cosγを磁界の強さおよび方向Aの方向余弦とすれば、Ax=Acosα、Ay=Acosβ、Az=Acosγであり、磁気センサのX軸方向の磁界検出コイルからはAxの出力が、Y軸方向の磁界検出コイルからはAyの出力が、Z軸方向の磁界検出コイルからはAzの出力が送出される。そして磁界の強さおよび方向Aは、A=(Ax2+Ay2+Az2)1/2で与えられる。従って、プローブの位置と姿勢を正確に取得するためには、センサとプローブが常に一体であることが必要である。これにより、単純な構成で各座標毎に独立した信号を同時に検出することができる。
請求項5は、前記トランスミッタは、前記プローブが移動する範囲に磁界が存在するように配置されることを特徴とする。
半球状の磁界を発生するトランスミッタを備え、演算部は、磁気センサにより検出した磁界の磁力線に対するベクトルを計算することにより、三次元位置情報及び姿勢情報を取得する。従って、センサ(プローブ)が常に磁界内に存在することが重要である。そこで本発明では、予め、プローブが移動する範囲を決定し、その範囲に磁界が存在するようにトランスミッタを配置する。これにより、プローブが磁界の範囲から逸脱することを事前に防止することができる。
半球状の磁界を発生するトランスミッタを備え、演算部は、磁気センサにより検出した磁界の磁力線に対するベクトルを計算することにより、三次元位置情報及び姿勢情報を取得する。従って、センサ(プローブ)が常に磁界内に存在することが重要である。そこで本発明では、予め、プローブが移動する範囲を決定し、その範囲に磁界が存在するようにトランスミッタを配置する。これにより、プローブが磁界の範囲から逸脱することを事前に防止することができる。
請求項6は、被検査対象物に超音波を送信する送信部及び該被検査対象物から反射した超音波を受信する受信部を備えたプローブと、該プローブの三次元位置情報及び姿勢情報を生成するセンサと、所定のエリアに磁界ベクトルを形成するトランスミッタと、前記プローブから取得された反射画像に基づいて該反射画像の三次元形状を演算する演算部と、を備え、前記演算部は、前記プローブを前記被検査対象物に対して移動させることで前記反射画像を取得すると共に、前記センサにより前記磁界ベクトルを読み取って前記トランスミッタの原点に対する当該プローブの三次元位置情報及び姿勢情報を検出し、前記反射画像の各座標点における前記三次元位置情報及び前記姿勢情報に基づいて前記各座標点の三次元座標をリアルタイムに演算することを特徴とする。
本発明は、請求項1及び2と同様の作用効果を奏する。
本発明は、請求項1及び2と同様の作用効果を奏する。
以上記載のごとく本発明によれば、超音波プローブにセンサを取り付け、超音波プローブを移動したときの位置情報と、姿勢情報を取得し、移動することにより得られた位置情報と姿勢情報に基づいて二次元の断層画像を演算部により演算して三次元形状を取得して、この結果を可視化して表示部に表示するので、従来使用していた超音波プローブにセンサを取り付けるだけで容易に二次元の断層画像から三次元画像を生成することができる。
また、演算部はリアルタイムに反射画像の各座標点における三次元位置情報及び姿勢情報に基づいて、各座標点の三次元座標をリアルタイムに演算するので、プローブを移動させるだけで、反射画像の三次元形状をリアルタイムに演算することができる。
また、演算部はリアルタイムに反射画像の各座標点における三次元位置情報及び姿勢情報に基づいて、各座標点の三次元座標をリアルタイムに演算するので、プローブを移動させるだけで、反射画像の三次元形状をリアルタイムに演算することができる。
また、被検査対象物とプローブの距離はこの磁界のエリア内に限定されるので、磁界のエリア内であればプローブの姿勢がどのように変化しても正確な姿勢情報を得ることができる。
また、センサは磁気センサであり、前記プローブの適所に固定されているので、単純な構成で各座標毎に独立した信号を同時に検出することができる。
また、予め、プローブが移動する範囲を決定し、その範囲に磁界が存在するようにトランスミッタを配置するので、プローブが磁界の範囲から逸脱することを事前に防止することができる。
また、センサは磁気センサであり、前記プローブの適所に固定されているので、単純な構成で各座標毎に独立した信号を同時に検出することができる。
また、予め、プローブが移動する範囲を決定し、その範囲に磁界が存在するようにトランスミッタを配置するので、プローブが磁界の範囲から逸脱することを事前に防止することができる。
以下、本発明を図に示した実施形態を用いて詳細に説明する。但し、この実施形態に記載される構成要素、種類、組み合わせ、形状、その相対配置などは特定的な記載がない限り、この発明の範囲をそれのみに限定する主旨ではなく単なる説明例に過ぎない。
図1は本発明の一実施形態に係る三次元画像生成装置を模式的に示す図である。三次元磁気センサ(以下、単にセンサと呼ぶ)2は、トランスミッタ4によって周囲に磁場をつくり、センサ2が磁場を読み取ることにより、トランスミッタ4の中心(原点)Oに対するセンサ2の三次元位置と姿勢を検出する。手に保持したプローブ1を被検査対象物6に沿って移動させることで、被検査対象物6全体の画像を生成する。この際センサ2は、プローブ1の三次元位置と姿勢をリアルタイムで検出する。
即ち、図1において、反射画像3上の点をu,vとおくと、反射画像座標(u,v)に対するセンサ2の中心Qに原点をおいたセンサ座標(xr,yr)の拡大、回転および平行移動を考慮するとu,vとxr,yrの関係は以下の式で表される。
ここで、超音波の発振平面5とセンサ座標のxr,yr軸が平面が平行になるようにセンサ2を固定し、zrを一定または0としている。
図1は本発明の一実施形態に係る三次元画像生成装置を模式的に示す図である。三次元磁気センサ(以下、単にセンサと呼ぶ)2は、トランスミッタ4によって周囲に磁場をつくり、センサ2が磁場を読み取ることにより、トランスミッタ4の中心(原点)Oに対するセンサ2の三次元位置と姿勢を検出する。手に保持したプローブ1を被検査対象物6に沿って移動させることで、被検査対象物6全体の画像を生成する。この際センサ2は、プローブ1の三次元位置と姿勢をリアルタイムで検出する。
即ち、図1において、反射画像3上の点をu,vとおくと、反射画像座標(u,v)に対するセンサ2の中心Qに原点をおいたセンサ座標(xr,yr)の拡大、回転および平行移動を考慮するとu,vとxr,yrの関係は以下の式で表される。
ここで、超音波の発振平面5とセンサ座標のxr,yr軸が平面が平行になるようにセンサ2を固定し、zrを一定または0としている。
式(1)の係数k11〜k32までを次のキャリブレーションによって求める。
即ち、式(1)の係数を求めるために、あらかじめ座標がわかっている(x,y)と(u,v)の複数の組み合わせを式(1)に代入し、係数k11〜k32を算出する。式(1)より、
sを消去して、
これを変形して、
即ち、式(1)の係数を求めるために、あらかじめ座標がわかっている(x,y)と(u,v)の複数の組み合わせを式(1)に代入し、係数k11〜k32を算出する。式(1)より、
sを消去して、
これを変形して、
エコー画面上での座標(ui,vi)をマウスでクリックすることで読みとる。また、超音波シート上のセンサ座標(xi,yi)をキーボードで入力する。これを繰り返し、4点の座標を読みとり式(4)に代入する。尚、なるべく広い範囲で点を指示してやる方が望ましい。これらのデータより、以下の行列を作成し、k11〜k32を計算する。
k11〜k32が決まると
式(4)を変形して、
行列で表すと
よって(xr,yr)は
式(8)が、u,vからxr,yrへの変換式である。
k11〜k32が決まると
式(4)を変形して、
行列で表すと
よって(xr,yr)は
式(8)が、u,vからxr,yrへの変換式である。
トランスミッタ4の中心(原点)Oに対するセンサ2の三次元位置を(xow,yow,zow)、姿勢を(φ,θ,ψ)とした場合、磁気レシーバの中心に原点を置いた磁気レシーバ座標(センサ座標)上の点(xr,yr,zr)は、以下の式で、トランスミッタを原点とした座標系(xrw,yrw,zrw)に変換される。zrは一定または0である。
ここで
C:Cos、S:Sinとする。
ここで
C:Cos、S:Sinとする。
図2は本発明の一実施形態に係る三次元画像生成装置の構成を示す図である。同じ構成要素には図1と同じ参照番号を付して説明する。この三次元画像生成装置110は、被検査対象物6に超音波を送信する送信部と被検査対象物6から反射した超音波を受信する受信部とを備えたプローブ1と、プローブ1の三次元位置情報及び姿勢情報を生成するセンサ2と、所定のエリアに磁界ベクトルを形成するトランスミッタ4と、三次元位置情報及び姿勢情報に基づいてプローブ1から取得された反射画像の三次元形状を演算するPC(演算部)12と、演算結果を三次元画像として表示するモニタ(表示部)14と、センサ2からの座標をリアルタイムに演算する位置検出回路10と、プローブ1から得られた反射画像データを処理するイメージプロセッサ11と、を備えて構成されている。尚、位置検出回路10はセンサ2の座標をリアルタイムに演算するために回路をFPGA(Field Programmable Gate Array)により構成している。またPC12には装置全体を制御するプログラムを格納するROM(Read Only Memory)13が接続されている。
即ち、プローブ1により取得された画像は、超音波が被検査対象物6内に進行して、固いものに当たって反射した反射音波を測定し、反射音が返ってくるまでの時間から距離を計算して可視化したものである。そして、複数の超音波を発生することで二次元の反射画像を生成することができる。この二次元の反射画像は被検査対象物6の1つの断層画像を表示しているに過ぎない。そこで本実施形態では、プローブ1にセンサ2を取り付け、プローブ1を移動したときの位置情報と姿勢情報を取得し、移動することにより得られた位置情報と姿勢情報に基づいて二次元の断層画像をPC12により演算して三次元形状を得る。この結果を可視化してモニタ14に表示する。これにより、従来使用していたプローブ1にセンサ2を取り付けるだけで容易に二次元の断層画像から三次元画像を生成することができる。
即ち、プローブ1により取得された画像は、超音波が被検査対象物6内に進行して、固いものに当たって反射した反射音波を測定し、反射音が返ってくるまでの時間から距離を計算して可視化したものである。そして、複数の超音波を発生することで二次元の反射画像を生成することができる。この二次元の反射画像は被検査対象物6の1つの断層画像を表示しているに過ぎない。そこで本実施形態では、プローブ1にセンサ2を取り付け、プローブ1を移動したときの位置情報と姿勢情報を取得し、移動することにより得られた位置情報と姿勢情報に基づいて二次元の断層画像をPC12により演算して三次元形状を得る。この結果を可視化してモニタ14に表示する。これにより、従来使用していたプローブ1にセンサ2を取り付けるだけで容易に二次元の断層画像から三次元画像を生成することができる。
また、PC12は、プローブ1を被検査対象物6に対して移動させることで反射画像を取得すると共に、センサ2により磁界ベクトルを読み取ってトランスミッタ4の原点Oに対する当該プローブ1の三次元位置情報及び姿勢情報を検出し、反射画像の各座標点における三次元位置情報及び姿勢情報に基づいて、各座標点の三次元座標をリアルタイムに演算する。
即ち、プローブ1は被検査対象物6に接触することで反射画像を取得することができる。また、プローブ1を被検査対象物6に沿って移動させることで、センサ2からプローブ1の三次元位置情報及び姿勢情報が検出される。PC12は反射画像の各座標点における三次元位置情報及び姿勢情報に基づいて、各座標点の三次元座標をリアルタイムに演算する。これにより、プローブ1を被検査対象物6に対して移動させるだけで、反射画像の三次元形状をリアルタイムに演算することができる。
また、PC12は、センサ2により受信した磁界ベクトルに基づいて、当該プローブ1の三次元位置情報及びロール角、ピッチ角、及びヨー角の姿勢情報を演算する。本発明では半球状の磁界を発生するトランスミッタ4を備え、その磁界の磁力線に対するベクトルを計算することにより三次元位置情報及び姿勢情報を取得するものである。従って、被検査対象物6とプローブ1の距離はこの磁界のエリア内に限定される。これにより、磁界のエリア内であればプローブ1の姿勢がどのように変化しても正確な姿勢情報を得ることができる(詳細は後述する)。
即ち、プローブ1は被検査対象物6に接触することで反射画像を取得することができる。また、プローブ1を被検査対象物6に沿って移動させることで、センサ2からプローブ1の三次元位置情報及び姿勢情報が検出される。PC12は反射画像の各座標点における三次元位置情報及び姿勢情報に基づいて、各座標点の三次元座標をリアルタイムに演算する。これにより、プローブ1を被検査対象物6に対して移動させるだけで、反射画像の三次元形状をリアルタイムに演算することができる。
また、PC12は、センサ2により受信した磁界ベクトルに基づいて、当該プローブ1の三次元位置情報及びロール角、ピッチ角、及びヨー角の姿勢情報を演算する。本発明では半球状の磁界を発生するトランスミッタ4を備え、その磁界の磁力線に対するベクトルを計算することにより三次元位置情報及び姿勢情報を取得するものである。従って、被検査対象物6とプローブ1の距離はこの磁界のエリア内に限定される。これにより、磁界のエリア内であればプローブ1の姿勢がどのように変化しても正確な姿勢情報を得ることができる(詳細は後述する)。
図3は本発明の三次元画像生成装置110により三次元画像を生成する過程を説明する模式図である。図3(A)は、例えば、四角錘(ピラミッド型)20をプローブ1により上部から読み取り、例えば、a、b、cの3箇所の反射画像について図示している。図3(B)は図3(A)のa部分の反射画像の模式図である。この図からa部分は厚み21がある大きな台形として捉えられる。図3(C)は図3(A)のb部分の反射画像の模式図である。この図からb部分は厚み21がある三角形として捉えられる。図3(D)は図3(A)のc部分の反射画像の模式図である。この図からc部分は厚み21がある小さな台形として捉えられる。図3(E)はプローブ1により上部から読み取った反射画像をプローブ1の位置情報と姿勢情報に基づいて、三次元座標を演算して三次元画像を表示したものである。尚、図3(E)の画像は、マウスにより指定することにより、回転させたり、部分的に切断して切断面を表示することもできる。
図4は本発明の三次元画像生成装置110により三次元画像を生成する過程を説明する模式図である。図4(A)は、例えば、トマト25を水中に沈めプローブ1を水面に接触させて上部から読み取り、例えば、a、b、c、d、eの5箇所の反射画像について図示している。図4(B)は図4(A)のa〜eの各部分の反射画像の模式図である。そして、図3と同様にプローブ1により上部から読み取った反射画像をプローブ1の位置情報と姿勢情報に基づいて、三次元座標を演算して三次元画像を生成することができる。図4(C)は実際にトマトを水中に入れて、プローブ1を水面に接触させて反射画像を取得し、それらの画像の三次元座標を演算して三次元画像を生成した例である。この図から、反射画像f〜kにより三次元画像が形成されている様子が分かる。
図5は本発明のセンサ44、45が位置情報と姿勢情報を生成する原理を説明する図である。図5(a)は磁界ベクトルと各センサの位置関係を表す図であり、トランスミッタ41から放射された磁界ベクトル42が半球状に形成され、その磁界ベクトル42上に例えばセンサ44、45が存在し、被検査対象物40が磁界ベクトル42の内側に含まれるものとする。図5(b)は磁界ベクトルを説明する模式図である。例えばセンサ44ではトランスミッタ41の位置における磁界の強さおよび方向をAとし、磁界の強さおよび方向AのX軸方向磁界の強さ、Y軸方向磁界の強さ、Z軸方向磁界の強さをそれぞれAx、Ay、Azとし、cosα、cosβ、cosγを磁界の強さおよび方向Aの方向余弦とすれば、Ax=Acosα、Ay=Acosβ、Az=Acosγであり、センサ44のX軸方向の磁界検出コイルからはAxの出力が、Y軸方向の磁界検出コイルからはAyの出力が、Z軸方向の磁界検出コイルからはAzの出力が送出される。そして磁界の強さおよび方向Aは、A=(Ax2+Ay2+Az2)1/2で与えられる。同じくセンサ45の磁界の強さおよび方向Bは、B=(Bx2+By2+Bz2)1/2で与えられる。
次に、本発明の三次元画像生成装置を超音波診断装置に応用した場合について言及する。超音波画像の画像分解能は、プローブ(超音波探触子)の材質、周波数等によって決まるが、そのほかプローブの型によっても若干違いが出る。通常、医療用では扇形に画像が得られるコンベックスが用いられるが、検体が比較的平であれば、プローブの幅でストレートに画像が得られるリニアモデルが利用される。そのほか、ペンシルビームのようなタイプで、検体との接触面がかなり小さくても比較的広い視野で画像を得るもの(セクター)もある。生体の周囲をトレースすることによって、例えば、胴体の腹側、背側、側面などのうち、どの部分に脂肪が蓄積しやすいのか、また内臓脂肪の動態を生体に何ら損傷を与えることなく診断できる利点がある。また、例えば、この輪切りの状態を首の側から腰の側へトレースするとこれまで、体構成割合(骨・筋肉・脂肪)としてしか捉えられなかったものが、リアルタイムの立体画像として捉えることができるようになる。現在は、リニアモデル(リアルタイムの2MHzの周波数)で図6(a)、(b)に示すような部分的な画像を得ているが、本発明では、胴体の周囲を回転させた立体画像として捉えられるようになる。ただ、超音波の伝播距離や画像の分解能の点から、適確な深度を得られる周波数や、プローブのモデル、さらにはプローブの材質の開発が必要となる。
また、これまで皮下脂肪30と内臓脂肪31の割合を測定するには腹部CTあるいはMRIを用いて臍高を基準として図7(a)のように測定している。この場合、臍高の横断面のみの脂肪量を反映するため、体型(下腹部肥満など図7(b))によっては脂肪全体の用量を正確に反映していなかった。これを本発明の三次元画像生成装置により三次元化することで、部位による脂肪蓄積の偏りが反映され、より正確な脂肪量を測定することが可能とる。また、これまで皮下脂肪については皮下の厚みを測っていたが、腹部での分布や偏りについてはあまり言及されていなかった。今後、立体的に分布を反映することができれば、体型による偏りが反映されることで、今まであまり議論されていない皮下脂肪と疾患についての研究に結びつく可能性があり、また、ダイエットをする際にも立体画像を提示することで、目に見える効果をより理解しやすく励みにもなるという効果がある。
また、測定装置について、現在のCTやMRIでは測定面積に制限があるため、高度肥満の患者は実際に測定器に入れない場合がある。本発明により測定器の測定面積が、既存のCTやMRIよりも広く設定できれば、これまで測定できなかった患者も検査が可能となる。
また、測定装置について、現在のCTやMRIでは測定面積に制限があるため、高度肥満の患者は実際に測定器に入れない場合がある。本発明により測定器の測定面積が、既存のCTやMRIよりも広く設定できれば、これまで測定できなかった患者も検査が可能となる。
1 プローブ、2 センサ、3 反射画像、4 トランスミッタ、5 発振平面、6 被検査対象物、10 位置検出回路、11 イメージプロセッサ、12 PC、13 ROM、14 モニタ、110 三次元画像生成装置
Claims (6)
- 被検査対象物に超音波を送信する送信部及び該被検査対象物から反射した超音波を受信する受信部を備えたプローブと、該プローブの三次元位置情報及び姿勢情報を生成するセンサと、所定のエリアに磁界ベクトルを形成するトランスミッタと、前記三次元位置情報及び前記姿勢情報に基づいて前記プローブから取得された反射画像の三次元形状を演算する演算部と、該演算部による演算結果を三次元画像として表示する表示部と、を備えたことを特徴とする三次元画像生成装置。
- 前記演算部は、前記プローブを移動させることで前記被検査対象物からの前記反射画像を取得すると共に、前記センサにより前記磁界ベクトルを読み取って前記トランスミッタの原点に対する当該プローブの三次元位置情報及び姿勢情報を検出し、前記反射画像の各座標点における前記三次元位置情報及び前記姿勢情報に基づいて前記各座標点の三次元座標をリアルタイムに演算することを特徴とする請求項1に記載の三次元画像生成装置。
- 前記演算部は、前記センサにより受信した磁界ベクトルに基づいて、前記プローブの三次元位置情報、及び該プローブのロール角、ピッチ角、ヨー角に係る姿勢情報を演算することを特徴とする請求項1又は2に記載の三次元画像生成装置。
- 前記センサは磁気センサであり、前記プローブの適所に固定されていることを特徴とする請求項1、2又は3に記載の三次元画像生成装置。
- 前記トランスミッタは、前記プローブが移動する範囲に磁界が存在するように配置されることを特徴とする請求項1乃至4の何れか一項に記載の三次元画像生成装置。
- 被検査対象物に超音波を送信する送信部及び該被検査対象物から反射した超音波を受信する受信部を備えたプローブと、該プローブの三次元位置情報及び姿勢情報を生成するセンサと、所定のエリアに磁界ベクトルを形成するトランスミッタと、前記プローブから取得された反射画像に基づいて該反射画像の三次元形状を演算する演算部と、を備え、
前記演算部は、前記プローブを前記被検査対象物に対して移動させることで前記反射画像を取得すると共に、前記センサにより前記磁界ベクトルを読み取って前記トランスミッタの原点に対する当該プローブの三次元位置情報及び姿勢情報を検出し、前記反射画像の各座標点における前記三次元位置情報及び前記姿勢情報に基づいて前記各座標点の三次元座標をリアルタイムに演算することを特徴とする三次元画像生成方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007211126A JP2009045097A (ja) | 2007-08-13 | 2007-08-13 | 三次元画像生成装置及び三次元画像生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007211126A JP2009045097A (ja) | 2007-08-13 | 2007-08-13 | 三次元画像生成装置及び三次元画像生成方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2009045097A true JP2009045097A (ja) | 2009-03-05 |
Family
ID=40497905
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007211126A Pending JP2009045097A (ja) | 2007-08-13 | 2007-08-13 | 三次元画像生成装置及び三次元画像生成方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2009045097A (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010259536A (ja) * | 2009-04-30 | 2010-11-18 | Canon Inc | 画像処理装置及びその制御方法 |
JP2014506990A (ja) * | 2010-10-25 | 2014-03-20 | ロッキード マーティン コーポレイション | 水中構造物の3次元モデルの構築 |
CN113611276A (zh) * | 2021-07-08 | 2021-11-05 | 北京小唱科技有限公司 | 声反馈抑制方法、装置及存储介质 |
-
2007
- 2007-08-13 JP JP2007211126A patent/JP2009045097A/ja active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010259536A (ja) * | 2009-04-30 | 2010-11-18 | Canon Inc | 画像処理装置及びその制御方法 |
JP2014506990A (ja) * | 2010-10-25 | 2014-03-20 | ロッキード マーティン コーポレイション | 水中構造物の3次元モデルの構築 |
CN113611276A (zh) * | 2021-07-08 | 2021-11-05 | 北京小唱科技有限公司 | 声反馈抑制方法、装置及存储介质 |
CN113611276B (zh) * | 2021-07-08 | 2024-06-11 | 北京小唱科技有限公司 | 声反馈抑制方法、装置及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Prager et al. | Rapid calibration for 3-D freehand ultrasound | |
Huang et al. | Fully automatic three-dimensional ultrasound imaging based on conventional B-scan | |
JP6430498B2 (ja) | 超音波剪断波エラストグラフィ測定のマッピングのためのシステムおよび方法 | |
Pagoulatos et al. | A fast calibration method for 3-D tracking of ultrasound images using a spatial localizer | |
JP6675305B2 (ja) | エラストグラフィ測定システム及びその方法 | |
Boctor et al. | Rapid calibration method for registration and 3D tracking of ultrasound images using spatial localizer | |
CN100469321C (zh) | 三维超声波检测方法 | |
US20120203107A1 (en) | Ultrasound measuring apparatus and control method thereof | |
US20160073980A1 (en) | Imaging apparatus and positioning apparatus | |
Huang et al. | 2.5-D extended field-of-view ultrasound | |
CN105120760B (zh) | 用于通过应力测量在成像期间测量超声探头上的力和扭矩的系统 | |
KR20140020486A (ko) | 초음파를 이용하여 조직의 탄성을 분석하는 방법 및 장치 | |
CN110167447A (zh) | 用于快速和自动超声探头校准的系统和方法 | |
CN104510495A (zh) | 被检体信息获取装置及其控制方法 | |
CN111772677A (zh) | 具有维度的生物组织弹性检测方法,检测系统和存储介质 | |
JP2014124319A (ja) | 超音波キャリブレーションシステム及び超音波キャリブレーション方法 | |
Prager et al. | Automatic calibration for 3-D free-hand ultrasound | |
JP2009045097A (ja) | 三次元画像生成装置及び三次元画像生成方法 | |
JP6518116B2 (ja) | 超音波診断システム | |
Abeysekera | Three dimensional ultrasound elasticity imaging | |
Dandekar et al. | A phantom with reduced complexity for spatial 3-D ultrasound calibration | |
JP2005160704A (ja) | 超音波診断装置 | |
Martin et al. | A miniature position and orientation locator for three dimensional echocardiography | |
WO2015099835A1 (en) | System and method for displaying ultrasound images | |
JP2006334273A (ja) | 超音波診断装置 |