JP4468677B2 - 超音波画像生成方法及び超音波画像生成プログラム - Google Patents
超音波画像生成方法及び超音波画像生成プログラム Download PDFInfo
- Publication number
- JP4468677B2 JP4468677B2 JP2003387281A JP2003387281A JP4468677B2 JP 4468677 B2 JP4468677 B2 JP 4468677B2 JP 2003387281 A JP2003387281 A JP 2003387281A JP 2003387281 A JP2003387281 A JP 2003387281A JP 4468677 B2 JP4468677 B2 JP 4468677B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- dimensional tomographic
- ultrasonic
- reference position
- tomographic image
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
- A61B8/14—Echo-tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/12—Diagnosis using ultrasonic, sonic or infrasonic waves in body cavities or body tracts, e.g. by using catheters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
- A61B8/5276—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts due to motion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/66—Analysis of geometric attributes of image moments or centre of gravity
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30028—Colon; Small intestine
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30061—Lung
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Description
上記従来の超音波診断装置は、被検体に対して超音波を送受信して得られた超音波エコーデータに基づき、ラジアル画像(超音波振動子の走査軸に対して垂直断面)やリニア画像(超音波振動子の走査軸に対して水平断面)等の2次元断層像から構築される3次元(立体)の超音波画像を生成する超音波画像生成方法を備えている。
通常、超音波診断装置は、超音波振動子をラジアル走査(超音波振動子の走査軸に対して垂直走査)すると共に、リニア走査(超音波振動子の走査軸に対して水平走査)することにより、ヘリカルスキャンを行って上記ラジアル画像及び上記リニア画像等の2次元断層像を得ている。ここで、上記ヘリカルスキャンは、例えば、2〜4秒程度の時間を必要としている。
しかしながら、上記特開2000−316864号公報に記載の超音波画像生成方法は、拍動除去に関してある程度の効果が得られているが、拍動除去がまだ十分ではない。このため、超音波画像生成方法としては、より効果の高い方法が求められている。
また、本発明による超音波画像生成方法は、被検体に対して超音波を送受信して得られた超音波エコーデータに基づき、超音波画像を生成する超音波画像生成方法において、連続する複数の2次元断層像に対して各画像上の基準位置を求める基準位置設定ステップと、前記基準位置設定ステップで求めたそれぞれの2次元断層像の前記基準位置の不規則性を補正して規則性のある連続的な2次元断層像を得る補正ステップと、前記補正ステップで補正した前記規則性のある連続的な2次元断層像に基づいて超音波画像を生成する超音波画像生成ステップと、を具備し、前記補正ステップは、予め指定された所定の2次元断層像の基準位置と、前記所定の2次元断層像と異なる他の2次元断層像の基準位置とを一致させるように当該他の2次元断層像を移動させるステップと、前記所定の2次元断層像の基準位置と、前記他の2次元断層像の基準位置とを一致させた後に、前記所定の2次元断層像を挟む前後の2次元断層像を用いて体腔表面の位置を平滑化するステップと、前記平滑化を実施した後に、平滑前の体腔表面の位置と平滑化した体腔表面の位置との差分を算出し、この算出した差分に基づいて該当する2次元断層像を伸縮することで、連続的な超音波画像を得るステップとを具備することを特徴とし、
また、本発明による超音波画像生成プログラムは、被検体に対して超音波を送受信して得られた超音波エコーデータに基づき、超音波画像を生成する超音波画像処理装置におけるコンピュータに、前記超音波エコーデータに基づいて連続的に得られる複数の2次元断層像に対して、2次元断層像上の輝度変化に基づく表面抽出点による多角形を求める手順と、前記多角形の重心を前記2次元断層像上の基準位置とする基準位置設定手順と、前記基準位置設定手順で求めたそれぞれの2次元断層像の前記基準位置の不規則性を補正して規則性のある連続的な2次元断層像を得る補正手順と、前記補正手順で補正した前記規則性のある連続的な2次元断層像に基づいて超音波画像を生成する超音波画像生成手順と、を実行させることを特徴とし、
また、本発明による超音波画像生成プログラムは、被検体に対して超音波を送受信して得られた超音波エコーデータに基づき、超音波画像を生成する超音波画像処理装置におけるコンピュータに、連続する複数の2次元断層像に対して各画像上の基準位置を求める基準位置設定手順と、前記基準位置設定手順で求めたそれぞれの2次元断層像の前記基準位置の不規則性を補正して規則性のある連続的な2次元断層像を得る補正手順と、前記補正手順で補正した前記規則性のある連続的な2次元断層像に基づいて超音波画像を生成する超音波画像生成手順と、を実行させ、前記補正手順は、予め指定された所定の2次元断層像の基準位置と、前記所定の2次元断層像と異なる他の2次元断層像の基準位置とを一致させるように当該他の2次元断層像を移動させる手順と、前記所定の2次元断層像の基準位置と、前記他の2次元断層像の基準位置とを一致させた後に、前記所定の2次元断層像を挟む前後の2次元断層像を用いて体腔表面の位置を平滑化する手順と、前記平滑化を実施した後に、平滑前の体腔表面の位置と平滑化した体腔表面の位置との差分を算出し、この算出した差分に基づいて該当する2次元断層像を伸縮することで、連続的な超音波画像を得る手順とをコンピュータに実行させることを特徴とする。
超音波プローブ2は、体腔内に挿入可能な細長で可撓性を有する挿入部11と、この挿入部11の基端側が着脱自在に連設される駆動部12とから構成される。挿入部11は、この先端部11aに超音波を送受信する超音波振動子2aが内蔵されている。
装置本体3は、超音波振動子2aに対して超音波信号を送受信し、3次元領域の超音波エコーデータを得る超音波観測部31と、体腔内(例えば、胃20内)での超音波プローブ2の挿入部先端部11aの位置データを得る位置検出部32と、この位置検出部32及び超音波観測部31のデータをもとに超音波画像データを得るもので、位置検出部32からの位置データを2次元断層像に関連付けてラジアルスキャン面の位置を特定し、後述する超音波画像生成方法及び超音波画像生成プログラムに基づき超音波画像(3次元画像)を生成するCPUを有する画像処理部33とを備えて構成されている。
そこで、本実施例では、以下に記載する超音波画像生成方法及び超音波画像生成プログラムにより拍動の影響を除去する。
先ず、図1で説明したように超音波プローブ2は、駆動部12の第1モータと第2モータとを同期させて同時に回動駆動させることで、超音波振動子2aが所定のピッチ単位でヘリカル走査する。
同時に、位置検出部32は、体腔内での超音波プローブ2の挿入部先端部11aの位置データを取得し、この取得した位置データを画像処理部33へ順次出力する。
そして、画像処理部33は、連続した複数のラジアル画像Grに対して各画像上の基準位置を求める基準位置設定ステップとして以下のステップS2及びS3の処理を行う。
ここで、画像処理部33は、例えば、図5に示す画像40において、画像中心から探索線40aを出して輝度変化の大きいところを表面座標として抽出し、表面抽出点41を決定する。
ここで、画像処理部33は、例えば、図6に示すように表面抽出点41を結んだ多角形42を生成し、重心43を求める。
本実施例では、物理的重心を用いて多角形の重心43を求めている。尚、重心43としては、幾何的重心を用いても良い。
また、画像処理部33は、多角形42を生成した後、この多角形42がほぼ内接するような円(不図示)に近似し、この円の中心点を求めて重心とするようにしても良い。
ここで、連続した複数のラジアル画像Grは、拍動の影響や、プローブ自体の機械的なぶれで、例えば、図7に示すように画像にずれが生じる。
そこで、本実施例では、画像処理部33は、複数のラジアル画像Grをそれぞれ移動させてこれら画像の位置を調整し、S3で求めた各画像上の重心を一致させる(重心を揃える)ような処理を行う(ステップS4)。
このとき、画像処理部33は、ユーザが予め指定したラジアル画像Gr51の画像位置を固定したまま、ラジアル画像Gr52の画像位置を上下左右方向に移動させても良く、また、その逆でも良い。このことにより、画像処理部33は、機械的な大きな「ぶれ」を軽減できる。
そこで、次に、画像処理部33は、前後の画像を用いてこれら画像に挟まれる所定のラジアル画像Grの表面抽出点を平滑化し、体腔表面の位置を平滑化する処理を行う(ステップS5)。
枚数=1/ピッチ×4…(1)
尚、ここで、ピッチとは、各2次元断層像(ラジアル画像Gr)間の距離である。
例えば、平滑化するのに参照するデータの枚数は、0.25ミリピッチの場合、16枚であり、1ミリピッチの場合4枚となる。
参照される画像番号=対象とする画像番号+(n−参照画像枚数/2)…(2)
但し、nは、0から始まり、参照画像枚数−1で終わる数値である。
そして、上述のように画像処理部33は、前後の画像を用いてこれら画像に挟まれる所定のラジアル画像Grの表面抽出点41を平滑化し、体腔表面の位置を平滑化する。
このため、平滑前の体腔表面61の表面位置(抽出された体腔表面の位置)と、平滑化した体腔表面62の表面位置とが一致するように、図10で示した画像60を変形する必要がある。
このことにより、画像処理部33は、平滑化後の体腔表面と合わせることが可能なラジアル画像Grを生成することができる。
そして、画像処理部33は、図11に示すように体腔表面72が滑らかな超音波画像(3次元画像)Gsを生成することが可能となる。
しかしながら、本来、ラジアル画像Grは、超音波エコーデータの生データとして例えば、図12に示す画像80のように数百本の音線データ81で形成され、これら音線データ81のエコー輝度により表される。
そこで、画像処理部33は、図示しないが音線データの形のまま、上記平滑化及び伸縮処理を極座標系で処理するようにしても良い。
また、第1実施例の超音波画像生成方法及び超音波画像生成プログラムは、位置検出部32により超音波プローブ2の挿入部先端部11aの位置データを得るように構成しているが、本発明はこれに限定されず、手動で擬似的に3次元領域の超音波エコーデータを取得し、超音波画像として3D表示(体腔表面を表示する斜視表面構築、表面は表示しない斜視構築)する際にも、本発明の超音波画像生成方法及び超音波画像生成プログラムを適用しても画質向上に有効である。
先ず、表面抽出処理について説明する。
表面抽出処理は、超音波伝達媒体や体液等の水と生体組織との識別を行う処理である。ここで、上記ラジアル画像Grは、例えば、数百本の音線と呼ばれる線データで構成され、これら音線のエコー輝度により表されている。図14に示すように各音線をプローブ中心(超音波振動子)から周縁方向に探索し、水から生体組織に変化する部分を探索する。この探索の精度を上げるため、次の処理を用いている。
次に、上記ヒストグラムのピークを判定する。
ピークの判定は、グラフ上各点の微分値(傾き)と絶対値(高さ)とから、判断する。通常、ヒストグラム上のピークは次の2パターンとなる。
2つのピークは、水部分と組織部分である。このため、閾値は、これら2つのピークの中間値とする。
(2)1つのピークが出現する。
コントラストの低い画像の場合、ピークが一つ(組織分)のみ、現れる場合が多い。この場合、ピークと最大輝度との中間値を閾値とする。
上記ノイズを除去する最も一般的な方法は、フレーム相関であるが、極力画像自体の平滑化を避けるために、本実施例では、対象物の厚みを測定して判断するようになっている。
そこで、ある一定の厚みを設定し、対象物がその設定値以上の場合、生体組織と判定する。設定値が小さいと、誤抽出が増え、大きすぎると生体組織の層構造をノイズと判定してしまう。
そこで、図17に示すように生体組織の層構造の存在を想定して判定する。
本実施例では、輝度平均を求める対象物の厚みを2mm程度とし、この2mm程度の厚みの輝度平均が閾値以上の場合、生体組織と判定するようになっている。
本実施例では、対象となる抽出点について、プローブ中心(超音波振動子)からの距離を、その前後の音線の距離と比較し、誤抽出の補正処理を行うようになっている。
音線Nは、誤抽出により、実際の体腔表面よりも手前の体腔表面を判定している。
Lを音線の長さ,Dを次の音線Nとの差分とした場合、
音線Nの予想抽出点Lxは、
Lx=(L1+L2+Lt+L3+L4)/5+(D12+D2t+Dt3+D34)/4
但し、L1=音線N−2の長さ,L2=音線N−1の長さ,
L3=音線N+1の長さ,L4=音線N+2の長さ,
D12=音線N−2と音線N−1との差分,D2t=音線N−1と音線Nとの差分,
Dt3=音線Nと音線N+1との差分,D34=音線N+1と音線N+2との差分
ここで、上記D12とD34の平均が、プローブ中心(超音波振動子)から体腔表面までの距離の増減傾向となる。
誤抽出点は、体腔表面位置より手前もしくは、後ろにある。そこで、予想抽出点Lxと実抽出点Ltとを比較し、実抽出点Ltが上記算出した予想抽出点Lxより3mm以上離れている場合、実抽出点Ltを上記算出した予想抽出点Lxに置き換えることにより、誤抽出を補正するようにしている。尚、上記誤抽出の補正処理は、図示しないが、音線Nの前後6点の音線N−3,N−2,N−1,N+1,N+2,N+3に基づき、音線Nの予想抽出点Lxを算出するようにしても良い。
本実施例では、抽出された体腔表面を極座標(音線番号と距離)から直交座標に変換し、この直交座標上において時計回りに3、6、9、12時の位置に配置される音線の表面距離(プローブ中心(超音波振動子)から抽出点までの距離)から、体腔中心を求める。具体的には、体腔中心(X,Y)は、
X=(|3時方向の表面抽出点のX座標|+|9時方向の表面抽出点のX座標|)/2
Y=(|12時方向の表面抽出点のY座標|+|6時方向の表面抽出点のY座標|)/2
となる。
次に、上記第1実施例で説明したのと同様に、前後の画像を用いてこれら画像に挟まれる所定のラジアル画像Grの表面抽出点を(直行座標値に変換済み)Z軸方向(長手方向)に平滑化し、体腔表面の位置を平滑化する処理を行う。
これにより、3Dモデル表示(サーフェイス表示)の場合、平滑化された座標値を利用することで、より自然で滑らかな体腔表面の構築が可能となる。
先ず、平滑化した体腔表面の表面位置(平滑化後の表面座標)と平滑化前の体腔表面の表面位置(平滑化前の表面座標)との距離差分を算出する。そして、平滑化前の体腔表面の表面位置(平滑化前の表面座標)と、平滑化した体腔表面の表面位置(平滑化後の表面座標)との距離Dを算出する。
(1)前提条件
a1 = Pn.y / Pn.x
b1 = 0
a2 = ( P2.y −P1.y ) / ( P2.x −P1.x )
b2 = ( P1.y −( P2.y −P1.y ) / ( P2.x −P1.x ) × P1.x )
(2)交点Qの座標
Q.x = ( b2 −b1 ) / ( a1 −a2 )
Q.y = a1 × ( b2 −b1 ) / ( a1 −a2 ) + b1
但し、R=Pnとする。
そして、上記交点Qと上記交点Rとからこれら交点QPの差分Dが算出できる。
ここで、ラジアル画像上の任意の1点P(x,y)を決定する音線データ上の点は、次の式で求まる。
(1)前提条件
T=2π/512×音線番号(N=0〜512)
L=原点Oから体腔表面までの距離,D=差分
(2)Pの座標
P.x=cosT×(L+D)
P.y=sinT×(L+D)
これにより、L(原点Oから体腔表面までの距離)に、上記算出した平滑化した体腔表面の表面位置(平滑化後の表面座標)と、平滑化前の体腔表面の表面位置(平滑化前の表面座標)との距離の差分Dを加味することで、平滑化した体腔表面の表面位置(平滑化後の表面座標)が平滑化前の体腔表面の表面位置(平滑化前の表面座標)に合わせ込む(一致させる)ことができる。即ち、所定のラジアル画像Grを伸縮する処理が可能となる。
そして、距離補正処理されたラジアル画像Grを用いることにより、拍動を軽減したリニア画像が構築されるようになっている。
先ず、図13で説明したように超音波プローブ2Bは、駆動部12の第1モータと第2モータとを同期させて同時に回動駆動させることで、超音波振動子2aが所定のピッチ単位でヘリカル走査する。
先ず、画像処理部33Bは、連続した複数のラジアル画像Grに対して各画像の表面座標を抽出し、表面抽出点を決定する(ステップS12)。
先ず、対象とする音線Nの前後の2本ずつ(N−2,N−1,N+1,N+2)を指定する(ステップS21)。
これにより、画像処理部33Bは、連続した複数のラジアル画像Grに対して各画像の表面座標を抽出できる。
(付記項1)
被検体に対して超音波を送受信して得られた超音波エコーデータに基づき、超音波画像を生成する超音波画像生成方法において、
連続する複数の2次元断層像に対して各画像上の基準位置を求める基準位置設定ステップと、
前記基準位置設定ステップで求めたそれぞれの2次元断層像の前記基準位置の不規則性を補正して規則性のある連続的な2次元断層像を得る補正ステップと、
前記補正ステップで補正した前記規則性のある連続的な2次元断層像に基づいて超音波画像を生成する超音波画像生成ステップと、
を有することを特徴とする超音波画像生成方法。
前記補正ステップは、予め指定された所定の2次元断層像の基準位置と、前記所定の2次元断層像と異なる他の2次元断層像の基準位置とを一致させるように当該他の2次元断層像を移動させることを特徴とする付記項1に記載の超音波画像生成方法。
前記補正ステップは、更に、前記所定の2次元断層像の基準位置と、前記他の2次元断層像の基準位置とを一致させた後に、前記所定の2次元断層像を挟む前後の2次元断層像を用いて体腔表面の位置を平滑化することを特徴とする付記項2に記載の超音波画像生成方法。
前記基準位置は、抽出された体腔表面を多角形とみなしたときの重心であることを特徴とする付記項1に記載の超音波画像生成方法。
前記補正ステップは、更に、前記平滑化を実施した後に、平滑前の体腔表面の位置と平滑化した体腔表面の位置との差分を算出し、この算出した差分に基づいて該当する2次元断層像を伸縮することで、連続的な超音波画像を得ることを特徴とする付記項3に記載の超音波画像生成方法。
前記基準位置は、抽出された体腔表面を円とみなしたときの重心であることを特徴とする付記項1に記載の超音波画像生成方法。
前記補正ステップは、前記平滑化を極座標系で実施することを特徴とする付記項3に記載の超音波画像生成方法。
前記補正ステップは、前記平滑化を直交座標系で実施することを特徴とする付記項3に記載の超音波画像生成方法。
前記補正ステップは、前記伸縮を極座標系で実施することを特徴とする付記項5に記載の超音波画像生成方法。
前記補正ステップは、前記伸縮を直交座標系で実施することを特徴とする付記項5に記載の超音波画像生成方法。
連続する断層像のそれぞれについて該断層像上の位置における共通する概念の基準位置を求める基準位置設定ステップと、
前記基準位置設定ステップで求められたそれぞれの断層像の基準位置の不規則性を補正する補正ステップと、
前記それぞれの断層像を、前記補正ステップで補正されたそれぞれの断層像の基準位置に対応させて移動し、この移動したそれぞれの断層像により連続的な超音波画像を得る超音波画像生成ステップと、
を有することを特徴とする超音波画像生成方法。
前記基準位置の移動は、事前に指定されたある特定の断層像の基準位置に一致するように移動することを特徴とする付記項11に記載の超音波画像生成方法。
前記それぞれの断層像を基準位置に対応させて移動させた後、それぞれの断層像における体腔表面の位置を前後の画像の体腔表面位置を用いて平滑化したことを特徴とする付記項11に記載の超音波画像生成方法。
前記共通する概念の基準位置は、抽出された体腔表面を多角形とみなしたときの重心であることを特徴とする付記項11に記載の超音波画像生成方法。
前記共通する概念の基準位置は、抽出された体腔表面を円と見なしたときの重心であることを特徴とする付記項11に記載の超音波画像生成方法。
前記基準位置の移動は、事前に指定されたある特定の断層像の基準位置に一致するように移動することを特徴とする付記項11に記載の超音波画像生成方法。
前記それぞれの断層像を基準位置に対応させて移動させた後、それぞれの断層像における体腔表面の位置を前後の画像の体腔表面位置を用いて平滑化したことを特徴とする付記項11に記載の超音波画像生成方法。
前記平滑化は、極座標系で実施することを特徴とする付記項14に記載の超音波画像生成方法。
(付記項19)
前記平滑化は、直交座標系で実施することを特徴とする付記項14に記載の超音波画像生成方法。
前記平滑化した後に、平滑化された表面抽出位置と元の断層像の表面抽出位置との差分を求め、この差分を用いて断層像を伸縮することで、連続的な超音波画像を得ることを特徴とする付記項14に記載の超音波画像生成方法。
前記伸縮は、極座標系で実施することを特徴とする付記項20に記載の超音波画像生成方法。
(付記項22)
前記伸縮は、直交座標系で実施することを特徴とする付記項20に記載の超音波画像生成方法。
前記基準位置は、抽出された体腔表面を極座標から直交座標に変換し、この直交座標上において時計回りに3、6、9、12時の位置に配置される抽出点に基づいて算出した体腔中心であることを特徴とする付記項1又は11に記載の超音波画像生成方法。
被検体に対して超音波を送受信して得られた超音波エコーデータに基づき、超音波画像を生成する超音波画像処理装置におけるコンピュータに、
連続する複数の2次元断層像に対して各画像上の基準位置を求める基準位置設定手順と、
前記基準位置設定手順で求めたそれぞれの2次元断層像の前記基準位置の不規則性を補正して規則性のある連続的な2次元断層像を得る補正手順と、
前記補正手順で補正した前記規則性のある連続的な2次元断層像に基づいて超音波画像を生成する超音波画像生成手順と、
を実行させることを特徴とする超音波画像生成プログラム。
前記補正手順は、予め指定された所定の2次元断層像の基準位置と、前記所定の2次元断層像と異なる他の2次元断層像の基準位置とを一致させるように当該他の2次元断層像を移動させることを特徴とする付記項24に記載の超音波画像生成プログラム。
前記補正手順は、更に、前記所定の2次元断層像の基準位置と、前記他の2次元断層像の基準位置とを一致させた後に、前記所定の2次元断層像を挟む前後の2次元断層像を用いて体腔表面の位置を平滑化することを特徴とする付記項25に記載の超音波画像生成プログラム。
前記補正手順は、更に、前記平滑化を実施した後に、平滑前の体腔表面の位置と平滑化した体腔表面の位置との差分を算出し、この算出した差分に基づいて該当する2次元断層像を伸縮することで、連続的な超音波画像を得ることを特徴とする付記項26に記載の超音波画像生成プログラム。
前記基準位置は、抽出された体腔表面を多角形とみなしたときの重心であることを特徴とする付記項24に記載の超音波画像生成プログラム。
前記基準位置は、抽出された体腔表面を極座標から直交座標に変換し、この直交座標上において時計回りに3、6、9、12時の位置に配置される抽出点に基づいて算出した体腔中心であることを特徴とする付記項24に記載の超音波画像生成プログラム。
2 超音波プローブ
2a 超音波振動子
3 装置本体(超音波画像処理装置)
11 挿入部
12 駆動部
31 超音波観測部
32 位置検出部
33 画像処理部
代理人 弁理士 伊藤 進
Claims (12)
- 被検体に対して超音波を送受信して得られた超音波エコーデータに基づき、超音波画像を生成する超音波画像生成方法において、
前記超音波エコーデータに基づいて連続的に得られる複数の2次元断層像に対して、2次元断層像上の輝度変化に基づく表面抽出点による多角形を求めるステップと、
前記多角形の重心を前記2次元断層像上の基準位置とする基準位置設定ステップと、
前記基準位置設定ステップで求めたそれぞれの2次元断層像の前記基準位置の不規則性を補正して規則性のある連続的な2次元断層像を得る補正ステップと、
前記補正ステップで補正した前記規則性のある連続的な2次元断層像に基づいて超音波画像を生成する超音波画像生成ステップと、
を有することを特徴とする超音波画像生成方法。 - 前記補正ステップは、予め指定された所定の2次元断層像の基準位置と、前記所定の2次元断層像と異なる他の2次元断層像の基準位置とを一致させるように当該他の2次元断層像を移動させることを特徴とする請求項1に記載の超音波画像生成方法。
- 前記補正ステップは、更に、前記所定の2次元断層像の基準位置と、前記他の2次元断層像の基準位置とを一致させた後に、前記所定の2次元断層像を挟む前後の2次元断層像を用いて体腔表面の位置を平滑化することを特徴とする請求項2に記載の超音波画像生成方法。
- 被検体に対して超音波を送受信して得られた超音波エコーデータに基づき、超音波画像を生成する超音波画像生成方法において、
連続する複数の2次元断層像に対して各画像上の基準位置を求める基準位置設定ステップと、
前記基準位置設定ステップで求めたそれぞれの2次元断層像の前記基準位置の不規則性を補正して規則性のある連続的な2次元断層像を得る補正ステップと、
前記補正ステップで補正した前記規則性のある連続的な2次元断層像に基づいて超音波画像を生成する超音波画像生成ステップと、
を具備し、
前記補正ステップは、予め指定された所定の2次元断層像の基準位置と、前記所定の2次元断層像と異なる他の2次元断層像の基準位置とを一致させるように当該他の2次元断層像を移動させるステップと、
前記所定の2次元断層像の基準位置と、前記他の2次元断層像の基準位置とを一致させた後に、前記所定の2次元断層像を挟む前後の2次元断層像を用いて体腔表面の位置を平滑化するステップと、
前記平滑化を実施した後に、平滑前の体腔表面の位置と平滑化した体腔表面の位置との差分を算出し、この算出した差分に基づいて該当する2次元断層像を伸縮することで、連続的な超音波画像を得るステップとを具備することを特徴とする超音波画像生成方法。 - 前記多角形は、前記2次元断層像の中心から放射状に延ばした複数の探索線上において輝度変化が最大の点を前記表面抽出点とし、前記表面抽出点を結ぶことによって得られることを特徴とする請求項1に記載の超音波画像生成方法。
- 前記基準位置は、抽出された体腔表面を極座標から直交座標に変換し、この直交座標上において時計回りに3、6、9、12時の位置に配置される抽出点に基づいて算出した体腔中心であることを特徴とする請求項1に記載の超音波画像生成方法。
- 被検体に対して超音波を送受信して得られた超音波エコーデータに基づき、超音波画像を生成する超音波画像処理装置におけるコンピュータに、
前記超音波エコーデータに基づいて連続的に得られる複数の2次元断層像に対して、2次元断層像上の輝度変化に基づく表面抽出点による多角形を求める手順と、
前記多角形の重心を前記2次元断層像上の基準位置とする基準位置設定手順と、
前記基準位置設定手順で求めたそれぞれの2次元断層像の前記基準位置の不規則性を補正して規則性のある連続的な2次元断層像を得る補正手順と、
前記補正手順で補正した前記規則性のある連続的な2次元断層像に基づいて超音波画像を生成する超音波画像生成手順と、
を実行させることを特徴とする超音波画像生成プログラム。 - 前記補正手順は、予め指定された所定の2次元断層像の基準位置と、前記所定の2次元断層像と異なる他の2次元断層像の基準位置とを一致させるように当該他の2次元断層像を移動させることを特徴とする請求項7に記載の超音波画像生成プログラム。
- 前記補正手順は、更に、前記所定の2次元断層像の基準位置と、前記他の2次元断層像の基準位置とを一致させた後に、前記所定の2次元断層像を挟む前後の2次元断層像を用いて体腔表面の位置を平滑化することを特徴とする請求項8に記載の超音波画像生成プログラム。
- 被検体に対して超音波を送受信して得られた超音波エコーデータに基づき、超音波画像を生成する超音波画像処理装置におけるコンピュータに、
連続する複数の2次元断層像に対して各画像上の基準位置を求める基準位置設定手順と、
前記基準位置設定手順で求めたそれぞれの2次元断層像の前記基準位置の不規則性を補正して規則性のある連続的な2次元断層像を得る補正手順と、
前記補正手順で補正した前記規則性のある連続的な2次元断層像に基づいて超音波画像を生成する超音波画像生成手順と、
を実行させ、
前記補正手順は、予め指定された所定の2次元断層像の基準位置と、前記所定の2次元断層像と異なる他の2次元断層像の基準位置とを一致させるように当該他の2次元断層像を移動させる手順と、
前記所定の2次元断層像の基準位置と、前記他の2次元断層像の基準位置とを一致させた後に、前記所定の2次元断層像を挟む前後の2次元断層像を用いて体腔表面の位置を平滑化する手順と、
前記平滑化を実施した後に、平滑前の体腔表面の位置と平滑化した体腔表面の位置との差分を算出し、この算出した差分に基づいて該当する2次元断層像を伸縮することで、連続的な超音波画像を得る手順とをコンピュータに実行させることを特徴とする超音波画像生成プログラム。 - 前記多角形は、前記2次元断層像の中心から放射状に延ばした複数の探索線上において輝度変化が最大の点を前記表面抽出点とし、前記表面抽出点を結ぶことによって得られることを特徴とする請求項7に記載の超音波画像生成プログラム。
- 前記基準位置は、抽出された体腔表面を極座標から直交座標に変換し、この直交座標上において時計回りに3、6、9、12時の位置に配置される抽出点に基づいて算出した体腔中心であることを特徴とする請求項7に記載の超音波画像生成プログラム。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003387281A JP4468677B2 (ja) | 2003-05-19 | 2003-11-17 | 超音波画像生成方法及び超音波画像生成プログラム |
EP04733985.8A EP1625829B8 (en) | 2003-05-19 | 2004-05-19 | Ultrasonic image generating method |
PCT/JP2004/007123 WO2004100795A1 (ja) | 2003-05-19 | 2004-05-19 | 超音波画像生成方法 |
US11/280,751 US8162835B2 (en) | 2003-05-19 | 2005-11-16 | Ultrasonic image generating method |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003140823 | 2003-05-19 | ||
JP2003387281A JP4468677B2 (ja) | 2003-05-19 | 2003-11-17 | 超音波画像生成方法及び超音波画像生成プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2005000638A JP2005000638A (ja) | 2005-01-06 |
JP4468677B2 true JP4468677B2 (ja) | 2010-05-26 |
Family
ID=33455513
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003387281A Expired - Fee Related JP4468677B2 (ja) | 2003-05-19 | 2003-11-17 | 超音波画像生成方法及び超音波画像生成プログラム |
Country Status (4)
Country | Link |
---|---|
US (1) | US8162835B2 (ja) |
EP (1) | EP1625829B8 (ja) |
JP (1) | JP4468677B2 (ja) |
WO (1) | WO2004100795A1 (ja) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070299345A1 (en) * | 2004-10-27 | 2007-12-27 | Hideo Adachi | Capacitive Ultrasonic Transducer and Endo Cavity Ultrasonic Diagnosis System Using the Same |
DE102005019369B3 (de) * | 2005-04-26 | 2006-11-16 | Siemens Ag | Vorrichtung zur Aufnahme von Querschnittsbildern |
US8502821B2 (en) * | 2008-02-04 | 2013-08-06 | C Speed, Llc | System for three-dimensional rendering of electrical test and measurement signals |
EP2248122A1 (en) * | 2008-02-25 | 2010-11-10 | Inventive Medical Limited | Medical training method and apparatus |
EP2107389B1 (en) * | 2008-04-02 | 2014-05-14 | Samsung Medison Co., Ltd. | Formation of an elastic image in an ultrasound system |
US8376951B2 (en) * | 2008-07-15 | 2013-02-19 | Hitachi Medical Corporation | Ultrasonic diagnostic apparatus and method for displaying probe operation guide |
KR101182999B1 (ko) | 2009-11-25 | 2012-09-18 | 삼성메디슨 주식회사 | 초음파 영상 처리를 수행하는 초음파 시스템 및 방법 |
US9443317B2 (en) * | 2011-09-09 | 2016-09-13 | Calgary Scientific Inc. | Image display of a centerline of tubular structure |
US9443303B2 (en) * | 2011-09-09 | 2016-09-13 | Calgary Scientific Inc. | Image display of a centerline of tubular structure |
US9311720B2 (en) | 2013-11-26 | 2016-04-12 | Saudi Arabian Oil Company | Automated saw cut correction for 3D core digital modeling from computerized tomography scanner (CTS) images |
WO2016140116A1 (ja) * | 2015-03-02 | 2016-09-09 | テルモ株式会社 | 画像診断装置および画像構築方法 |
WO2018133098A1 (zh) * | 2017-01-23 | 2018-07-26 | 上海联影医疗科技有限公司 | 血管壁应力应变状态获取方法及系统 |
CN110710990B (zh) * | 2019-09-29 | 2021-07-02 | 华中科技大学 | 一种螺旋超声断层成像方法及系统 |
WO2022202310A1 (ja) * | 2021-03-26 | 2022-09-29 | テルモ株式会社 | プログラム、画像処理方法及び画像処理装置 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3263171B2 (ja) * | 1993-04-01 | 2002-03-04 | 株式会社東芝 | 超音波診断装置 |
US5622174A (en) | 1992-10-02 | 1997-04-22 | Kabushiki Kaisha Toshiba | Ultrasonic diagnosis apparatus and image displaying system |
JP3253409B2 (ja) * | 1993-03-31 | 2002-02-04 | 株式会社東芝 | 超音波ドプラ診断装置 |
JP2785636B2 (ja) | 1993-02-25 | 1998-08-13 | 株式会社エス.エス.ビー | 生体組織多次元可視装置 |
JP3413379B2 (ja) | 1993-08-05 | 2003-06-03 | オリンパス光学工業株式会社 | 超音波画像診断装置 |
JP3316268B2 (ja) | 1993-08-30 | 2002-08-19 | アロカ株式会社 | 管腔内超音波診断装置 |
JP3450937B2 (ja) * | 1995-06-08 | 2003-09-29 | アロカ株式会社 | 超音波画像処理装置 |
JPH10192A (ja) | 1996-04-15 | 1998-01-06 | Olympus Optical Co Ltd | 超音波画像診断装置 |
US5830145A (en) * | 1996-09-20 | 1998-11-03 | Cardiovascular Imaging Systems, Inc. | Enhanced accuracy of three-dimensional intraluminal ultrasound (ILUS) image reconstruction |
JPH10248844A (ja) * | 1997-03-12 | 1998-09-22 | Olympus Optical Co Ltd | 超音波画像診断装置 |
US6106466A (en) * | 1997-04-24 | 2000-08-22 | University Of Washington | Automated delineation of heart contours from images using reconstruction-based modeling |
US6095976A (en) * | 1997-06-19 | 2000-08-01 | Medinol Ltd. | Method for enhancing an image derived from reflected ultrasound signals produced by an ultrasound transmitter and detector inserted in a bodily lumen |
US5885218A (en) * | 1997-11-07 | 1999-03-23 | Scimed Life Systems, Inc. | Method and apparatus for spatial filtering in an intravascular ultrasound imaging system |
US6210328B1 (en) | 1998-10-01 | 2001-04-03 | Atl Ultrasound | Ultrasonic diagnostic imaging system with variable spatial compounding |
US6117081A (en) | 1998-10-01 | 2000-09-12 | Atl Ultrasound, Inc. | Method for correcting blurring of spatially compounded ultrasonic diagnostic images |
JP2000316864A (ja) * | 1999-05-11 | 2000-11-21 | Olympus Optical Co Ltd | 超音波診断装置 |
JP2000316865A (ja) | 1999-05-12 | 2000-11-21 | Olympus Optical Co Ltd | 超音波画像診断装置 |
DE10054814A1 (de) * | 2000-11-04 | 2002-05-08 | Philips Corp Intellectual Pty | Verfahren und Vorrichtung zur Registrierung von Bildern |
-
2003
- 2003-11-17 JP JP2003387281A patent/JP4468677B2/ja not_active Expired - Fee Related
-
2004
- 2004-05-19 EP EP04733985.8A patent/EP1625829B8/en not_active Expired - Fee Related
- 2004-05-19 WO PCT/JP2004/007123 patent/WO2004100795A1/ja active Application Filing
-
2005
- 2005-11-16 US US11/280,751 patent/US8162835B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
WO2004100795A1 (ja) | 2004-11-25 |
US8162835B2 (en) | 2012-04-24 |
JP2005000638A (ja) | 2005-01-06 |
EP1625829A1 (en) | 2006-02-15 |
EP1625829B1 (en) | 2019-02-27 |
EP1625829A4 (en) | 2010-08-04 |
US20060079772A1 (en) | 2006-04-13 |
EP1625829B8 (en) | 2019-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8162835B2 (en) | Ultrasonic image generating method | |
JP5294340B2 (ja) | 超音波診断装置 | |
JP7313359B2 (ja) | 超音波診断装置および超音波診断装置の制御方法 | |
JP2006510413A (ja) | 動脈壁の動きを決定する超音波ドップラシステム | |
US20150164482A1 (en) | Ultrasound diagnostic apparatus, ultrasound image recording method, and non-transitory computer readable recording medium | |
JP6419976B2 (ja) | 超音波診断装置および超音波診断装置の制御方法 | |
JP5166154B2 (ja) | 超音波診断装置 | |
WO2020008746A1 (ja) | 音響波診断装置および音響波診断装置の制御方法 | |
CN111970973A (zh) | 超声波诊断装置及超声波诊断装置的控制方法 | |
WO2016024422A1 (ja) | 超音波画像処理装置 | |
JP6853372B2 (ja) | 超音波診断装置および超音波診断装置の制御方法 | |
JP7313545B2 (ja) | 超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサ | |
JP4768315B2 (ja) | 超音波信号処理装置及び超音波信号処理方法 | |
CN112367921A (zh) | 声波诊断装置及声波诊断装置的控制方法 | |
US11534143B2 (en) | Acoustic wave diagnostic apparatus and control method of acoustic wave diagnostic apparatus | |
JP3283456B2 (ja) | 超音波画像診断装置及び超音波画像処理方法 | |
JP7434529B2 (ja) | 超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサ | |
JP4530834B2 (ja) | 超音波画像処理方法及び超音波画像処理装置、並びに、超音波画像処理プログラム | |
JP4445744B2 (ja) | 超音波画像生成装置の作動方法及び超音波画像生成プログラム | |
JP3905419B2 (ja) | 超音波画像処理装置 | |
EP4356843A1 (en) | Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus | |
US20240130706A1 (en) | Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus | |
US20230233186A1 (en) | Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus | |
JP4225830B2 (ja) | 超音波診断装置 | |
JP2024048667A (ja) | 超音波診断装置および超音波診断装置の制御方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20061005 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090804 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20091005 |
|
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: 20100202 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20100225 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 4468677 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130305 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140305 Year of fee payment: 4 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
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 |