Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP4127827B2 - 超音波診断装置 - Google Patents
[go: Go Back, main page]

JP4127827B2 - 超音波診断装置 - Google Patents

超音波診断装置 Download PDF

Info

Publication number
JP4127827B2
JP4127827B2 JP2004194888A JP2004194888A JP4127827B2 JP 4127827 B2 JP4127827 B2 JP 4127827B2 JP 2004194888 A JP2004194888 A JP 2004194888A JP 2004194888 A JP2004194888 A JP 2004194888A JP 4127827 B2 JP4127827 B2 JP 4127827B2
Authority
JP
Japan
Prior art keywords
speed
speed range
interest
range
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
Application number
JP2004194888A
Other languages
English (en)
Other versions
JP2006014891A (ja
Inventor
浩仲 宮木
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Olympus Corp
Original Assignee
Olympus Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Olympus Corp filed Critical Olympus Corp
Priority to JP2004194888A priority Critical patent/JP4127827B2/ja
Priority to EP05753303A priority patent/EP1769746A4/en
Priority to PCT/JP2005/011763 priority patent/WO2006003866A1/ja
Priority to CNB200580021841XA priority patent/CN100463658C/zh
Publication of JP2006014891A publication Critical patent/JP2006014891A/ja
Priority to US11/648,285 priority patent/US20070167770A1/en
Application granted granted Critical
Publication of JP4127827B2 publication Critical patent/JP4127827B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8979Combined Doppler and pulse-echo imaging systems
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52053Display arrangements
    • G01S7/52057Cathode ray tube displays
    • G01S7/52071Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Clinical applications
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/488Diagnostic techniques involving Doppler signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/02Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
    • G01S15/50Systems of measurement, based on relative movement of the target
    • G01S15/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S15/582Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse-modulated waves and based upon the Doppler effect resulting from movement of targets
    • G01S15/584Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse-modulated waves and based upon the Doppler effect resulting from movement of targets with measures taken for suppressing velocity ambiguities, i.e. anti-aliasing

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Acoustics & Sound (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • General Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Hematology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

この発明は、生体内に超音波を複数回照射するとともにそのエコーを順次受波して同一音線方向毎に超音波スキャンを繰り返し行い、この超音波スキャンによって得られた複数の超音波データをもとに、この生体内にて動く運動体の速度を示すカラー画像である速度画像を生成出力する超音波診断装置に関するものである。
従来から、生体内に超音波を照射するとともにそのエコーを受波する超音波スキャンを行ってこの生体内の超音波断層像とこの生体内にて動く運動体の速度を示す速度画像とを生成出力する超音波診断装置が、生体内の病変部等の関心領域の断層像または血液等の運動体の速度をリアルタイムに観察可能な医療用診断装置として普及されている。たとえば、超音波診断装置は、生体内の運動体に対する超音波スキャンを行って得られた超音波データを用いてドプラ法に基づく演算処理等を行えば、この運動体の速度を求めることができる。また、超音波診断装置は、操作者によって検査目的とする所望の関心速度レンジが予め設定され、この関心速度レンジ内の速度毎に輝度または色相等が割り当てられたカラースケールデータを用いることによって、この運動体の速度を示す速度画像を生成出力することができる。
しかし、超音波診断装置によって検出された運動体の速度がこの設定された関心速度レンジ外である場合、超音波診断装置には、この運動体の速度を誤った値および方向を示す色相または輝度の速度画像にて表示出力する現象(エリアジング)が発生する。このエリアジングは、サンプリングの定理によって決まる現象であり、この関心速度レンジを狭くするほどその発生頻度が高くなる。このエリアジングが発生した場合、超音波診断装置が運動体の実際の速度と異なる値および方向を示す速度画像を表示出力するので、操作者は、検査目的とする運動体の速度たとえば血流の流速値および血流方向を正確に認識できない。このエリアジングの発生を抑制可能な超音波診断装置として、関心領域内の運動体の速度たとえば関心血流の流速値が関心速度レンジ(すなわち流速レンジ)の上限値に達することによって飽和している画素数をもとにこの流速レンジがこの関心血流に対して適正か否かを判定し、その判定結果に応じて流速レンジを拡大する処理を自動的に行う超音波診断装置がある(特許文献1参照)。
特開平11−146879号公報
しかしながら、上述した特許文献1に記載された超音波診断装置では、予め設定された関心速度レンジとしての流速レンジを拡大することによって上述したエリアジングの発生を抑制しているので、検査目的外の運動体の速度に対応する周波数成分を除去するフィルタ処理のカットオフ周波数がこの流速レンジの拡大に応じて増加する場合が多い。この場合、このフィルタ処理によって検査目的外の運動体の速度を除去する速度レンジ(除去速度レンジ)がこの流速レンジの拡大すなわちカットオフ周波数の増加に応じて拡大し、これによって、本来の検査目的としての流速レンジ内の低速度成分を検出することが困難になるという問題点があった。
この発明は、上記事情に鑑みてなされたものであって、検査目的とする関心速度レンジ内の低速度成分の検出能力を損なうことなく、エリアジングの発生を抑制できる超音波診断装置を提供することを目的とする。
上述した課題を解決し、目的を達成するために、本発明にかかる超音波診断装置は、被検体内に対して超音波の送受信を複数回繰り返し行って得た複数の超音波データをもとに、前記被検体内の超音波断層像を生成出力するとともに、前記被検体内にて動く運動体の速度を所定速度レンジ内の速度として算出し、該算出した速度とカラースケールデータとをもとに前記運動体の速度を示す速度画像を生成出力する超音波診断装置において、前記運動体の関心速度レンジを示す情報を指示入力する入力手段と、前記入力手段から入力された情報をもとに、前記関心速度レンジ以上に広い速度レンジであって該関心速度レンジを含む可変の検出可能速度レンジを前記所定速度レンジとして設定する速度レンジ設定制御手段と、前記検出可能速度レンジ内の速度毎に前記カラースケールデータを割り当て、該割り当てたカラースケールデータと前記算出した速度とをもとに前記速度画像を生成する画像処理制御手段と、を備えたことを特徴とする。
また、本発明にかかる超音波診断装置は、上記発明において、前記速度レンジ設定制御手段は、前記関心速度レンジに応じた零近傍の速度レンジ部分を除去する前記検出可能速度レンジを前記所定速度レンジとして設定することを特徴とする。
また、本発明にかかる超音波診断装置は、被検体内に対して超音波の送受信を複数回繰り返し行って得た複数の超音波データをもとに、前記被検体内の超音波断層像を生成出力するとともに、前記被検体内にて動く運動体の速度を所定速度レンジ内の速度として算出し、該算出した速度とカラースケールデータとをもとに前記運動体の速度を示す速度画像を生成出力する超音波診断装置において、前記運動体の関心速度レンジを示す情報を指示入力する入力手段と、前記入力手段から入力された情報をもとに、前記関心速度レンジ以上に広い速度レンジであって該関心速度レンジを含む可変の検出可能速度レンジを前記所定速度レンジとして設定し、さらに前記関心速度レンジに応じた零近傍の速度レンジ部分を除去する前記検出可能速度レンジを設定する速度レンジ設定制御手段と、を備えたことを特徴とする。
また、本発明にかかる超音波診断装置は、上記発明において、前記検出可能速度レンジ内の速度毎に割り当てられたカラースケールデータに対応するカラースケール画像と前記速度画像とを含む各種画像を複数同時に表示出力する表示手段を備え、前記画像処理制御手段は、前記表示手段に対して前記各種画像の表示制御を行うことを特徴とする。
また、本発明にかかる超音波診断装置は、上記発明において、前記画像処理制御手段は、前記表示手段に対し、前記関心速度レンジ内の前記カラースケール画像よりも前記関心速度レンジ外であって前記検出可能速度レンジ内の前記カラースケール画像を小さくする表示制御を行うことを特徴とする。
また、本発明にかかる超音波診断装置は、上記発明において、前記画像処理制御手段は、前記検出可能速度レンジ内の速度毎に前記カラースケールデータを割り当てる場合、前記関心速度レンジ外であって前記検出可能速度レンジ内の速度変化に対する前記カラースケールデータの色相変化または輝度変化よりも、前記関心速度レンジ内の速度変化に対する前記カラースケールデータの色相変化または輝度変化を大きくすることを特徴とする。
また、本発明にかかる超音波診断装置は、上記発明において、前記速度レンジ設定制御手段は、前記入力手段から入力された情報をもとに前記関心速度レンジを設定するとともに該関心速度レンジの最大速度値に対応する基準繰り返し周波数と前記超音波の送受信の繰り返し回数に関する仮の繰り返し周波数とを算出し、該仮の繰り返し周波数をスイープする周波数スイープ処理を行うとともに該周波数スイープ処理毎の仮の繰り返し周波数を用いて前記超音波の送受信制御を順次行う送受信制御手段と、前記送受信制御手段による前記超音波の送受信制御によって得られた前記複数の超音波データをもとに前記運動体の速度を順次算出する速度演算制御手段と、を備え、前記送受信制御手段は、前記速度演算制御手段から前記運動体の速度を順次検出するとともに、該順次検出した運動体の速度をもとにエリアジング現象の発生有無を順次判定するエリアジング判定処理を行い、該エリアジング判定処理の結果をもとに前記超音波の送受信制御を行う実際の繰り返し周波数と前記検出可能速度レンジとを設定し、前記速度演算制御手段は、前記基準繰り返し周波数を少なくとも用い、前記検出可能速度レンジから前記関心速度レンジに応じた零近傍の速度レンジ部分を除去する速度レンジを設定することを特徴とする。
この発明によれば、所望の低速度範囲の速度を検出する能力を損なうことなく、被検体内の所望の運動体の検出速度に対するエリアジングの発生を抑制可能な超音波診断装置を実現できるという効果を奏する。
また、この発明によれば、このエリアジングの発生を抑制できるとともに、被検体内から検出した運動体の速度を示す速度画像を確実に画面表示できる超音波診断装置を実現することができるという効果を奏する。
以下、添付図面を参照して、この発明にかかる超音波診断装置の好適な実施の形態を詳細に説明する。なお、この実施の形態によって、この発明が限定されるものではない。
(実施の形態1)
図1は、この発明の実施の形態1である超音波診断装置の一構成例を例示するブロック図である。図1において、この超音波診断装置1は、入力部2、超音波振動子3、送受信回路4、Bモードデータ演算部5、グレー画像データ生成部6、速度データ演算部7、カラー画像データ生成部8、画像合成部9、モニタ10、記憶部11、および制御部12を有する。
入力部2は、キーボード、タッチパネル、トラックボール、マウス、またはロータリスイッチ等を単一または組み合わせて用いることによって実現され、制御部12と電気的に接続される。入力部2は、操作者による情報入力操作に応じ、超音波診断装置1の各構成部が行う動作の開始、終了、または切り替え等を指示する各種指示情報、超音波診断装置1の各構成部が行う演算処理の各種パラメータ情報、グレー画像データを生成する場合に用いられるグレースケールデータに関するグレースケール情報、あるいはカラー画像データを生成する場合に用いられるカラースケールデータに関するカラースケール情報等を制御部12に入力する。
たとえば、入力部2は、操作者による情報入力操作に応じ、制御部12の動作モードをBモード、カラードプラ画像モード、またはティシュードプラ画像モードに切替指示する動作モード指示情報を制御部12に入力する。また、入力部2は、操作者による情報入力操作に応じ、被検体内の関心領域にて動く所望の運動体の速度(関心速度)の速度レンジ(関心速度レンジ)を指示する関心速度レンジ指示情報を制御部12に入力する。なお、このBモードは、被検体内の超音波断層像すなわちBモード画像をモニタ10に表示出力する動作モードである。また、このカラードプラ画像モードは、血液等の比較的高速度で動く運動体の速度を検出するとともに速度画像として表示出力する速度表示モードであり、この速度画像としてカラードプラ画像を表示出力する動作モードである。さらに、このティシュードプラ画像モードは、生体組織等の比較的低速度で動く運動体の速度を検出するとともに速度画像として表示出力する速度表示モードであり、この速度画像としてティシュードプラ画像を表示出力する動作モードである。一方、この被検体内にて動く運動体は、たとえば血液、生体組織、または被検体内に入れられた超音波造影剤等であって、被検体内において超音波振動子3に対して動く各種物体である。
超音波振動子3は、チタン酸バリウムまたはチタン酸ジルコン酸鉛等の圧電材料が複数配列されたアレイ振動子を用いて実現され、送受信回路4と電気的に接続される。超音波振動子3は、送受信回路4から受信した電気的パルス信号を逆圧電効果によって音響的パルス信号すなわち超音波に変換出力する機能と、この変換出力した音響的パルス信号の反射信号(エコー信号)を圧電効果によって電気的パルス信号に変換するとともにこの電気的パルス信号を送受信回路4に出力する機能とを有する。この場合、超音波振動子3は、送受信回路4から順次受信した電気的パルス信号に基づき、たとえば被検体内に音響的パルス信号を順次送信するとともにこの被検体内からのエコー信号を順次受信し、このエコー信号に対応する電気的パルス信号を送受信回路4に順次送信する。すなわち、超音波振動子3は、被検体内の各音線方向について送受信回路4から電気的パルス信号を複数回繰り返し受信することによって、被検体内の音線方向毎に音響的パルス信号を同じ複数回送信するとともにこのエコー信号を同じ複数回受信する動作を繰り返す。この場合、超音波振動子3は、送受信回路4の制御のもと、被検体内の1断層面(1フレーム)毎に複数回の超音波スキャンを行うことができる。
送受信回路4は、上述した電気的パルス信号を超音波振動子3に順次送信するとともに超音波振動子3によって変換出力された電気的パルス信号を順次受信する送受信制御を行うビームフォーミング回路を用いて実現され、超音波振動子3、Bモードデータ演算部5、および速度データ演算部7とそれぞれ電気的に接続される。送受信回路4は、制御部12から受信した制御信号に基づき、上述した被検体内の各音線方向について複数回繰り返し送信する電気的パルス信号の繰り返し周波数を設定する。また、送受信回路4は、制御部12の制御のもと、上述した各音線方向について、予め設定された可変の繰り返し回数から電気的パルス信号の繰り返し回数を決定するとともに、この決定した繰り返し回数の電気的パルス信号を繰り返し送受信する。これによって、送受信回路4は、制御部12の制御のもと、被検体内の1フレーム毎に上述した複数回の超音波スキャンによる複数の超音波データを得ることができる。さらに、送受信回路4は、制御部12の動作モードがBモードである場合、制御部12の制御のもと、これら複数の超音波データをBモードデータ演算部5に送信する。一方、送受信回路4は、制御部12の動作モードが速度表示モードである場合、制御部12の制御のもと、Bモードデータ演算部5および速度データ演算部7に対し、これら複数の超音波データを被検体内の1フレーム毎に交互に送信する。
なお、送受信回路4は、特公平6−2134号公報に開示されている方法とほぼ同様に、同一の音線方向に送受信される音響的パルス信号に対応する電気的パルス信号の送受信を繰り返し行うとともに他の音線方向に送受信される音響的パルス信号に対応する電気的パルス信号の送受信を行ってもよいし、同一の音線方向に送受信される音響的パルス信号に対応する電気的パルス信号の送受信を所定回数繰り返し行い、その後、別の音線方向に送受信される音響的パルス信号に対応する電気的パルス信号の送受信を繰り返すようにしてもよい。
Bモードデータ演算部5は、送受信回路4から受信した超音波データをもとに被検体内の超音波断層像(Bモード画像)に対応するBモードデータを演算出力する周知の演算処理回路を用いて実現され、送受信回路4およびグレー画像データ生成部6とそれぞれ電気的に接続される。具体的には、Bモードデータ演算部5は、制御部12の制御のもと、被検体内の1フレーム毎に送受信回路4から順次受信する複数の超音波データを用いてバンドパスフィルタ処理、ログ圧縮処理、ゲイン調整処理、コントラスト調整処理、およびフレーム相関処理等を行い、被検体内の1フレーム毎のBモード画像に対応するBモードデータを演算出力する。この場合、Bモードデータ演算部5は、被検体内の1フレーム毎に送受信回路4から受信した複数の超音波データを用い、被検体内の1フレーム毎のBモードデータを順次演算出力してもよいし、被検体内の3次元領域に配列された複数のBモード画像に対応するBモードデータを順次演算出力してもよい。また、Bモードデータ演算部5は、得られたBモードデータをグレー画像データ生成部6に送信する。
グレー画像データ生成部6は、Bモードデータ演算部5から受信したBモードデータ、所定のグレースケールデータ、および所定のルックアップテーブルをもとにグレー画像データを生成出力する周知の演算処理回路を用いて実現され、Bモードデータ演算部5および画像合成部9とそれぞれ電気的に接続される。具体的には、グレー画像データ生成部6は、制御部12に制御のもと、このグレースケールデータおよびルックアップテーブルを用い、Bモードデータ演算部5から受信したBモードデータをグレー画像データに順次変換出力する。このグレー画像データは、このBモードデータに対応するBモード画像をグレー画像としてモニタ9に表示出力するための画像データである。グレー画像データ生成部6は、変換出力したグレー画像データを画像合成部9に順次送信する。
また、グレー画像データ生成部6は、RAM(Random Access Memory)またはROM(Read Only Memory)等を用いた記憶部(図示せず)を有し、このグレースケールデータおよびルックアップテーブルを更新可能に保持する。この場合、グレー画像データ生成部6は、操作者が入力部2を用いて所望のグレースケール情報の入力操作を行うことによって、制御部12を介してこの所望のグレースケール情報に対応する所望のグレースケールデータを更新可能に保持することができる。このグレースケールデータは、上述したBモードデータの値に応じて光の3原色(赤、緑、青)の各輝度を変えた色のデータであって、このBモードデータに対する色の変化特性を変更可能である。すなわち、グレー画像データ生成部6は、この所望のグレースケールデータを用いることによって、上述したBモードデータに対応する所望輝度のグレー画像データを生成出力できる。
速度データ演算部7は、送受信回路4およびカラー画像データ生成部8とそれぞれ電気的に接続され、制御部12の制御のもと、送受信回路4から受信した超音波データと制御部12から受信したパラメータ信号とをもとに上述した運動体の速度を演算出力するように機能する。すなわち、速度データ演算部7は、制御部12の動作モードが速度表示モードである場合、制御部12の制御のもと、被検体内の1フレーム毎に送受信回路4から順次受信する複数の超音波データと制御部12からのパラメータ信号に基づく各種パラメータとをもとに、被検体内の1フレーム毎の各空間位置について、この運動体の速度を演算出力する。この場合、速度データ演算部7は、被検体内の1フレーム毎の各空間位置について、検査目的の運動体の速度に対応する速度データ(関心速度データ)または検査目的外の運動体の速度に対応する速度データ(目的外速度データ)をカラー画像データ生成部8に送信する。なお、速度データ演算部7の詳細な構成については、後述する。
カラー画像データ生成部8は、速度データ演算部7から受信した各種速度データ、所定のカラースケールデータ、および所定のルックアップテーブルをもとに各種カラー画像データを生成出力する演算処理回路を用いて実現され、速度データ演算部7および画像合成部9とそれぞれ電気的に接続される。すなわち、カラー画像データ生成部8は、制御部12の動作モードが速度表示モードである場合、制御部12の制御のもと、被検体内の1フレーム毎の各空間位置について、このカラースケールデータおよびルックアップテーブルを用い、速度データ演算部7から受信した関心速度データをカラー画像データに順次変換出力するとともに、速度データ演算部7から受信した目的外速度データを目的外カラー画像データに順次変換出力する。このカラー画像データは、この関心速度データに対応する運動体の速度をカラー表示するカラー画像すなわち速度画像をモニタ9に表示出力するための画像データである。一方、この目的外カラー画像データは、この目的外速度データに対応する運動体の速度に対して所定色たとえば黒色が割り当てられた画像データであって、モニタ9に表示出力させない画像データである。カラー画像データ生成部8は、被検体内の1フレーム毎の各空間位置について、変換出力したカラー画像データと変換出力した目的外カラー画像データとを画像合成部9に順次送信する。この場合、カラー画像データ生成部8は、この目的外速度データに対応する運動体の速度のカラー画像をモニタ9に表示出力させない制御信号として、この目的外カラー画像データを画像構成部9に送信する。
また、カラー画像データ生成部8は、RAMまたはROM等を用いた記憶部(図示せず)を有し、このカラースケールデータおよびルックアップテーブルを更新可能に保持する。この場合、カラー画像データ生成部8は、操作者が入力部2を用いて所望のカラースケール情報の入力操作を行うことによって、制御部12を介してこの所望のカラースケール情報に対応する所望のカラースケールデータを更新可能に保持することができる。このカラースケールデータは、光の3原色(赤、緑、青)の輝度または色相等に関する所定の組合せによって構成される色のデータであって、上述したカラースケール情報に基づきこの組合せが変更可能である。カラー画像データ生成部8は、保持するカラースケールデータと制御部12からのパラメータ信号とを用い、上述したエリアジングを起こさずに被検体内の所望の運動体の速度を検出可能な速度レンジ(検出可能速度レンジ)内の速度毎にこのカラースケールデータの各輝度または各色相をそれぞれ割り当てる。この検出可能速度レンジは、上述した関心速度レンジ以上に広い速度レンジであって、この関心速度レンジを少なくとも含む速度レンジである。すなわち、カラー画像データ生成部8は、このパラメータ信号に基づく検出可能速度レンジ内の速度毎にこのカラースケールデータの各輝度または各色相を割り当てることによって、この検出可能速度レンジに含まれる関心速度レンジ内の速度毎にこのカラースケールデータの各輝度または各色相を割り当てるとともに、この検出可能速度レンジに含まれる関心速度レンジ外の速度毎にこのカラースケールデータの残りの各輝度または各色相を割り当てることができる。
画像合成部9は、制御部12の動作モードが速度表示モードである場合、制御部12の制御のもと、グレー画像データ生成部6から受信したグレー画像データとカラー画像データ生成部8から受信したカラー画像データまたは目的外カラー画像データとを上述した被検体内の1フレーム毎の各空間位置について合成し、合成画像データを得る。この場合、画像合成部9は、被検体内の1フレーム毎の各空間位置において、グレー画像データにカラー画像データを上書きするとともに目的外カラー画像データにグレー画像データを上書きする。すなわち、この合成画像データは、被検体内の1フレーム毎の各空間位置において、各カラー画像データに対応する各空間位置にカラー画像データがそれぞれ構成され、それ以外の各空間位置に各グレー画像データがそれぞれ構成される。その後、画像合成部9は、得られた合成画像データを表示画像データに変換するとともに、この表示画像データをモニタ10に送信する。モニタ10は、画像合成部9と電気的に接続され、画像合成部9から受信した表示画像データをもとに、この合成画像データに対応する超音波断層像および速度画像を表示出力する。すなわち、モニタ10は、画像合成部9から順次受信した表示画像データ毎に合成画像データに対応する超音波断層像および速度画像をリアルタイムに順次更新する。
一方、画像合成部9は、制御部12の動作モードがBモードである場合、制御部12の制御のもと、グレー画像データ生成部6から受信したグレー画像データを表示画像データに変換するとともに、この表示画像データをモニタ10に送信する。この場合、モニタ10は、画像合成部9から受信した表示画像データをもとに、このグレー画像データに対応する超音波断層像を表示出力する。すなわち、モニタ10は、画像合成部9から順次受信した表示画像データ毎にグレー画像データに対応する超音波断層像をリアルタイムに順次更新する。
また、画像合成部9は、グレー画像データ生成部6から直接または制御部12を介してグレースケールデータが入力された場合、このグレースケールデータを表示画像データに変換するとともに、この表示画像データをモニタ10に送信する。この場合、モニタ10は、画像合成部9から受信した表示画像データをもとに、このグレースケールデータに対応するグレースケールを表示出力する。これと同様に、画像合成部9は、カラー画像データ生成部8から直接または制御部12を介してカラースケールデータが入力された場合、このカラースケールデータを表示画像データに変換するとともに、この表示画像データをモニタ10に送信する。この場合、モニタ10は、画像合成部9から受信した表示画像データをもとに、このカラースケールデータに対応するカラースケールを表示出力する。これによって、モニタ10は、グレースケールと超音波断層像とを同一画面上に表示出力でき、またはグレースケールおよびカラースケールと超音波断層像および速度画像とを同一画面上に表示出力できる。
記憶部11は、EEPROMまたはフラッシュメモリ等の各種ICメモリ、ハードディスクドライブ、あるいは光磁気ディスクドライブ等のデータの書き込みおよび読み出しが可能な各種記憶装置を用いて実現される。記憶部11は、制御部12の制御のもと、制御部12から入力された合成画像データ、グレー画像データ、またはカラー画像データ等の各種画像データを記憶する。また、記憶部11は、制御部12の制御のもと、各種パラメータ情報、グレースケール情報、またはカラースケール情報等の制御部12から入力された各種情報を記憶する。さらに、記憶部11は、制御部12の制御のもと、これらの記憶した各種情報を制御部12に送信する。
制御部12は、処理プログラム等の各種データが予め記憶されたROMと、演算パラメータ等を一時的に記憶するRAMと、この処理プログラムを実行するCPUとを用いて実現される。制御部12は、入力部2、送受信回路4、Bモードデータ演算部5、グレー画像データ生成部6、速度データ演算部7、カラー画像データ生成部8、画像合成部9、および記憶部11と電気的に接続され、上述したように、これらの各構成部の動作および各種情報の入出力を制御する。
また、制御部12は、入力部2から入力された動作モード指示情報に基づき、Bモード、カラードプラ画像モード、またはティシュードプラ画像モードに動作モードを切り替える。その後、制御部12は、この動作モードに応じ、上述したように、送受信回路4、Bモードデータ演算部5、グレー画像データ生成部6、速度データ演算部7、カラー画像データ生成部8、および画像合成部9の各動作と情報入出力とを制御する。
さらに、制御部12は、入力部2から入力された関心速度レンジ指示情報に基づき、被検体内の関心領域にて動く所望の運動体の関心速度レンジ±Viを一意的に設定するとともに、この関心速度レンジ±Viを関心速度の検出可能な速度レンジとして設定した場合の繰り返し周波数である基準繰り返し周波数fiを算出する。なお、関心速度レンジ±Viは、その最小速度−Viからその最大速度Viまでの各速度を含む速度レンジとして定義する。その後、制御部12は、この基準繰り返し周波数fiを用い、送受信回路4が送受信する電気的パルス信号の繰り返し回数を制御するための実繰り返し周波数frを設定するとともに、この実繰り返し周波数frに対応する制御信号を送受信回路4に送信する。この場合、送受信回路4は、この制御部12からの制御信号に基づき、上述した電気的パルス信号の繰り返し周波数を設定する。また、送受信回路4は、予め設定された可変の繰り返し回数をもとに、送受信する電気的パルス信号の繰り返し回数を決定する。
また、制御部12は、この関心速度レンジ±Viを少なくとも用い、被検体内の関心領域にて動く所望の運動体の検出可能速度レンジ±Vrを設定する。この場合、制御部12は、検出可能速度レンジ±Vrを関心速度レンジ±Viに対して可変に設定する。なお、検出可能速度レンジ±Vrは、その最小速度−Vrからその最大速度Vrまでの各速度を含む速度レンジとして定義する。さらに、制御部12は、被検体内の関心領域にて動く目的外の運動体の速度成分を除去するためのカットオフ周波数fcを算出する。その後、制御部12は、実繰り返し周波数frおよびカットオフ周波数fcに対応する各パラメータ信号を速度データ演算部7に送信する。また、制御部12は、関心速度レンジ±Viおよび検出可能速度レンジ±Vrに対応する各パラメータ信号をカラー画像データ生成部8に送信する。
つぎに、速度データ演算部7の構成について詳細に説明する。図2は、速度データ演算部7の一構成を詳細に例示するブロック図である。図2において、速度データ演算部7は、複素信号化回路71、フィルタ72、自己相関回路73、運動情報演算部74、およびしきい値処理回路75を有する。
複素信号化回路71は、直交検波器を用いて実現され、送受信回路4から受信した超音波データに対応する電気的パルス信号を複素信号に変換するように機能する。具体的には、複素信号化回路71は、互いに位相が90度異なる正弦波信号とこの送受信回路4から受信した電気的パルス信号とを用いて乗算処理を行い、得られた電気信号をローパスフィルタに通すことによってこの複素信号を取得する。その後、複素信号化回路71は、この取得した複素信号をフィルタ72に送信する。
たとえば、制御部12の動作モードがカラードプラ画像モードまたはティシュードプラ画像モードである場合に、送受信回路4は、制御部12の制御のもと、所望の運動体を検出する音線方向毎に電気的パルス信号の送受信を上述した繰り返し回数(たとえば8回程度)ずつ繰り返し行う。この場合、複素信号化回路71は、所望の運動体を検出する音線方向毎に、この繰り返し回数と同数(たとえば8つ)の超音波データ群に対応する電気的パルス信号群を送受信回路4から受信するとともに、この電気的パルス信号群を変換した複素信号群を取得する。すなわち、複素信号化回路71は、被検体内の関心領域における2次元空間または3次元空間の各位置について、所望の運動体を検出した超音波データ群に対応する複素信号群を取得する。この場合、複素信号化回路71は、この取得した複素信号群をフィルタ72に送信する。
なお、複素信号化回路71は、RAM等のメモリ(図示せず)を有するようにし、この取得した複素信号群を保持してもよい。また、複素信号化回路71は、この取得した複素信号群を制御部12に送信し、制御部12がこの複素信号群を保持管理してもよい。
フィルタ72は、DSP(Digital Signal Processor)またはFPGA(Field Programmable Gate Array)等を用いたデジタルFIR(Finite Impulse Response)フィルタまたはデジタルIIR(Infinite Impulse Response)フィルタを用いて実現される。フィルタ72は、制御部12の制御のもと、複素信号化回路71から順次受信した複素信号群の実数信号群および虚数信号群についてそれぞれ個別にフィルタ処理を行う。
たとえば、制御部12の動作モードがカラードプラ画像モードである場合、制御部12は、上述したカットオフ周波数fcに対応するパラメータ信号をフィルタ72に送信する。フィルタ72は、制御部12からこのパラメータ信号を受信するとともに、受信したパラメータ信号に基づき、このフィルタ処理のカットオフ周波数fcを設定する。この場合、フィルタ72は、血液等に例示される比較的高速度運動の運動体の速度を検出する場合にこの複素信号群に対してフィルタ処理を行う周知のMTIフィルタとして機能し、この複素信号群から低周波成分すなわち変動の少ない成分をノイズとして除去する。このことは、この複素信号群から比較的低速度運動の運動体の速度に対応する成分を除去することに相当する。その後、フィルタ72は、このフィルタ処理によって低周波成分を除去した実数信号群および虚数信号群からなるフィルタ処理後の複素信号群を自己相関回路73に送信する。
一方、制御部12の動作モードがティシュードプラ画像モードである場合、フィルタ72は、制御部12の制御のもと、所定のフィルタ係数を設定するとともに、生体組織等に例示される比較的低速度運動の運動体の速度を検出する場合にこの複素信号群に対してフィルタ処理を行う周知のローパスフィルタとして機能する。この場合、フィルタ72は、この複素信号群から高周波成分すなわち変動の大きい成分をノイズとして除去する。このことは、この複素信号群から比較的高速度運動の運動体の速度に対応する成分を除去することに相当する。その後、フィルタ72は、このフィルタ処理によって低周波成分を除去した実数信号群および虚数信号群からなるフィルタ処理後の複素信号群を自己相関回路73に送信する。なお、フィルタ72は、制御部12の動作モードがティシュードプラ画像モードである場合、制御部12の制御のもと、フィルタとしての機能を停止させてもよい。この場合、フィルタ72は、複素信号化回路71から受信した複素信号群にフィルタ処理を行わずに、この複素信号群を自己相関回路73に送信する。
自己相関回路73は、DSPまたはFPGA等を用いて実現され、フィルタ72から受信した複素信号群をもとに、この複素信号群の複素自己相関値Rを算出する。たとえば、この複素信号群に含まれるN個(N:2以上の整数)の複素信号のうちのa番目の複素信号を示す複素数Zaは、次式(1)によって表される。
a=xa+jya (a=1〜N) ・・・(1)
この場合、自己相関回路73は、次式(2)に基づき、この複素信号群の複素自己相関値Rを算出する。
Figure 0004127827

なお、式(2)の複素数Za *は、複素数Zaと共役な複素数である。自己相関回路73は、式(2)に基づき算出した複素自己相関値Rに対応する電気信号を運動情報演算部74に出力する。
運動情報演算部74は、DSPまたはFPGA等を用いて実現され、制御部12の制御のもと、被検体内の1フレーム毎の各空間位置について所望の運動体の速度Vとエコー強度Iとをそれぞれ算出する。具体的には、運動情報演算部74は、自己相関回路73から受信した電気信号に基づく複素自己相関値Rと、制御部12から受信したパラメータ信号に基づく実繰り返し周波数frと、音速cと、送受信回路4が送受信する電気的パルス信号の中心周波数f0とを用い、次式(3)に基づいて速度Vを算出する。また、運動情報演算部74は、次式(4)に基づき、エコー強度Iを算出する。
Figure 0004127827

I=|R| ・・・(4)
なお、式(3)において、運動情報演算部74は、複素自己相関値Rに基づきその実数成分Rxとその虚数成分Ryとを取得する。また、運動情報演算部74は、実繰り返し周波数frの逆数として周期Tを取得する。すなわち、周期Tは、送受信回路4が被検体内の同一音線方向毎に電気的パルス信号を繰り返し送受信する場合の動作周期である。
その後、運動情報演算部74は、被検体内の1フレーム毎の各空間位置について、式(3)に基づき算出した速度Vに対応する電気信号と式(4)に基づき算出したエコー強度Iに対応する電気信号とをしきい値処理回路75に出力する。また、運動情報演算部74は、制御部12の制御のもと、式(3)に基づき算出した速度Vに対応する電気信号を制御部12に送信する。これによって、制御部12は、運動情報演算部74が算出した運動体の速度Vをリアルタイムにて検出できる。
なお、運動情報演算部74は、RAM等のメモリ(図示せず)を有するようにし、音速cおよび中心周波数f0等の演算パラメータを予め保持してもよい。また、運動情報演算部74は、音速cおよび中心周波数f0等の演算パラメータを制御部12からのパラメータ信号に基づき取得してもよい。
しきい値処理回路75は、DSPまたはFPGA等を用いて実現され、制御部12の制御のもと、運動情報演算部74が算出した運動体の速度Vが速度画像としてモニタ10に表示すべき運動体の速度であるか否かを判定する表示判定処理を行う。この場合、しきい値処理回路75は、運動情報演算部74から受信した各電気信号に基づき、被検体内の1フレーム毎の各空間位置の速度Vおよびエコー強度Iをそれぞれ取得し、取得した各速度Vと速度に関する所定の速度しきい値とを比較するとともに、取得した各エコー強度Iとエコー強度に関する所定の強度しきい値とを比較することによって、この表示判定処理を行う。
たとえば、制御部12の動作モードがカラードプラ画像モードである場合、しきい値処理回路75は、制御部12の制御のもと、被検体内の1フレーム毎の各空間位置の各速度Vと速度しきい値VTH1とを比較し、これら各速度Vが次式(5)を満足するか否かを判定する。
|V|>VTH1 ・・・(5)
これと同時に、しきい値処理回路75は、制御部12の制御のもと、被検体内の1フレーム毎の各空間位置の各エコー強度Iと強度しきい値ITH1,ITH2とを比較し、これら各エコー強度Iが次式(6)を満足するか否かを判定する。
TH1<I<ITH2 ・・・(6)
なお、速度しきい値VTH1は、これらの各速度Vが血液等の比較的高速度にて動く運動体の速度であるか否かを判定するためのしきい値である。強度しきい値ITH1は、取得したエコー強度Iがノイズであるか否かを判定するためのしきい値である。強度しきい値ITH2は、これらの各速度Vにて動く運動体が生体組織等の固体であるか血液等の流体であるかを判定するためのしきい値である。
この場合、しきい値処理回路75は、式(5)を満足する速度Vを検査目的の比較的高速度にて動く運動体の速度と判定する。また、しきい値処理回路75は、強度しきい値ITH1以下のエコー強度Iをノイズと判定し、強度しきい値ITH2未満のエコー強度Iに関連する速度Vを血液等の流体の速度と判定し、強度しきい値ITH2より大きいエコー強度Iに関連する速度Vを生体組織等の固体の速度と判定する。したがって、しきい値処理回路75は、カラードプラ画像モードにおいて、式(5)を満足しかつ式(6)を満足するエコー強度Iに関連する速度Vを速度画像としてモニタ10に表示すべき所望の運動体の速度と判定する。これによって、しきい値処理回路75は、被検体内の1フレーム毎の各空間位置の速度Vが速度画像としてモニタ10に表示すべき所望の運動体たとえば血液の速度であるか否かを判定できる。その後、しきい値処理回路75は、被検体内の1フレーム毎の各空間位置について、式(5)を満足しかつ式(6)を満足するエコー強度Iに関連する速度Vを上述した関心速度データとしてカラー画像データ生成部8に送信し、それ以外の速度Vを上述した目的外速度データとしてカラー画像データ生成部8に送信する。この場合、しきい値処理回路75は、式(5)を満足しない速度Vまたは式(6)を満足しないエコー強度Iに関連する速度Vを速度零に置き換え、この速度零のデータを上述した目的外速度データとしてカラー画像データ生成部8に送信してもよい。
なお、これらの速度しきい値VTH1または強度しきい値ITH1,ITH2の各最適値は実験的に取得することができる。また、しきい値処理回路75は、これらの速度しきい値VTH1または強度しきい値ITH1,ITH2を適切に設定すれば、血液以外の運動体についても同様に表示判定処理を行うことができる。さらに、しきい値処理回路75は、RAM等のメモリ(図示せず)を有するようにし、速度しきい値VTH1または強度しきい値ITH1,ITH2を予め保持してもよい。また、しきい値処理回路75は、速度しきい値VTH1または強度しきい値ITH1,ITH2を制御部12からのパラメータ信号に基づき取得してもよい。
一方、制御部12の動作モードがティシュードプラ画像モードである場合、しきい値処理回路75は、制御部12の制御のもと、被検体内の1フレーム毎の各空間位置の各速度Vと速度しきい値VTH2とを比較し、これら各速度Vが次式(7)を満足するか否かを判定する。
|V|<VTH2 ・・・(7)
これと同時に、しきい値処理回路75は、制御部12の制御のもと、被検体内の1フレーム毎の各空間位置の各エコー強度Iと強度しきい値ITH3とを比較し、これら各エコー強度Iが次式(8)を満足するか否かを判定する。
I>ITH3 ・・・(8)
なお、速度しきい値VTH2は、これらの各速度Vが生体組織等の比較的低速度にて動く運動体の速度であるか否かを判定するためのしきい値である。強度しきい値ITH3は、これらの各速度Vにて動く運動体が生体組織等の固体であるか否かを判定するためのしきい値である。
この場合、しきい値処理回路75は、式(7)を満足する速度Vを検査目的の比較的低速度にて動く運動体の速度と判定する。また、しきい値処理回路75は、強度しきい値ITH3より大きいエコー強度Iに関連する速度Vを生体組織等の固体の速度と判定する。したがって、しきい値処理回路75は、ティシュードプラ画像モードにおいて、式(7)を満足しかつ式(8)を満足するエコー強度Iに関連する速度Vを速度画像としてモニタ10に表示すべき所望の運動体の速度と判定する。これによって、しきい値処理回路75は、被検体内の1フレーム毎の各空間位置の速度Vが速度画像としてモニタ10に表示すべき所望の運動体たとえば生体組織の速度であるか否かを判定できる。その後、しきい値処理回路75は、被検体内の1フレーム毎の各空間位置について、式(7)を満足しかつ式(8)を満足するエコー強度Iに関連する速度Vを上述した関心速度データとしてカラー画像データ生成部8に送信し、それ以外の速度Vを上述した目的外速度データとしてカラー画像データ生成部8に送信する。この場合、しきい値処理回路75は、式(7)を満足しない速度Vまたは式(8)を満足しないエコー強度Iに関連する速度Vを速度零に置き換え、この速度零のデータを上述した目的外速度データとしてカラー画像データ生成部8に送信してもよい。
なお、これらの速度しきい値VTH2または強度しきい値ITH3の各最適値は実験的に取得することができる。また、しきい値処理回路75は、速度しきい値または強度しきい値を適切に設定すれば、生体組織以外の運動体についても同様に表示判定処理を行うことができる。さらに、しきい値処理回路75は、速度しきい値VTH2または強度しきい値ITH3を予め保持してもよいし、制御部12から受信したパラメータ信号に基づき取得してもよい。
つぎに、制御部12が、カラードプラ画像モードにおいて被検体内の運動体の速度を示す速度画像すなわちカラードプラ画像をモニタ10に表示出力させるまでの処理手順について詳細に説明する。図3は、被検体内の運動体の速度画像をモニタ10に表示出力させるまでに制御部が行う処理手順を例示するフローチャートである。図4は、被検体内のBモード画像とカラードプラ画像とを含むモニタ画像の一表示例を模式的に例示する模式図である。
図3において、まず、操作者は、入力部2を用い、被検体内にて動く所望の運動体についての関心のある速度レンジすなわち上述した関心速度レンジを選択する操作を行う。たとえば、入力部2には、この選択される関心速度レンジとして、超音波振動子3の種類、被検体の観察部位、または送受信する音響的パルス信号の周波数等によって決定する所望数の選択肢が設定される。操作者は、この入力部2に設定された選択肢たとえば5cm/s、10cm/s、20cm/s、および40cm/s等の関心速度レンジを示す所望の速度間隔の選択肢の中から、所望の関心速度レンジを選択する操作を行う。この場合、入力部2は、操作者が選択した関心速度レンジを設定指示する関心速度レンジ指示情報を制御部12に入力する。制御部12は、入力部2から入力された関心速度レンジ指示情報を検知し(ステップS101,Yes)、検知した関心速度レンジ指示情報に基づき、上述した関心速度レンジ±Viを設定する。これと同時に、制御部12は、次式(9)に基づき、上述した基準繰り返し周波数fiを算出する(ステップS102)。
Figure 0004127827

なお、式(9)において、最大速度Viは、上述したように、関心速度レンジ±Vi内の最大速度である。
一方、操作者が入力部2を用いて関心速度レンジを選択する操作を行わなければ、制御部12は、関心速度レンジ指示情報を検知せずに(ステップS101,No)、このステップS101の処理手順を繰り返す。
つぎに、制御部12は、ステップS102にて算出した基準繰り返し周波数fiと予め設定された可変の係数パラメータα(α:1以上の実数)とを用い、上述した実繰り返し周波数frを設定する実繰り返し周波数設定処理を行う(ステップS103)。この場合、制御部12は、次式(10)に基づく実繰り返し周波数frを取得するとともに、上述したように、この取得した実繰り返し周波数frを示すパラメータ信号を運動情報演算部74に送信する。
fr=α×fi ・・・(10)
その後、制御部12は、ステップS103にて設定した実繰り返し周波数frと上述した中心周波数f0および音速cとを用いて上述した検出可能速度レンジ±Vrの最大速度Vrを算出し、この算出した最大速度Vrとこの最大速度Vrの符号を反転した最小速度−Vrとによる検出可能速度レンジ±Vrを設定する(ステップS104)。なお、制御部12によって設定された関心速度レンジ±Viの最大速度Viが式(9)に基づき次式(11)によって表されるので、制御部12は、次式(12)に基づき、この最大速度Vrを算出できる。
Figure 0004127827
Figure 0004127827
つぎに、制御部12は、上述した関心速度レンジ±Viおよび検出可能速度レンジ±Vrを取得した場合、この関心速度レンジ±Viの最大速度Viと検出可能速度レンジ±Vrの最大速度Vrと予め設定された係数パラメータβ(β:正の小数)とを用い、次式(13)に基づき、上述したカットオフ周波数fcを算出する。これと同時に、制御部12は、この算出したカットオフ周波数fcを示すパラメータ信号をフィルタ72に送信することによって、フィルタ72にカットオフ周波数fcの設定指示を行う(ステップS105)。この場合、フィルタ72は、このカットオフ周波数fcをフィルタ処理のカットオフ周波数として設定するとともに、上述したMTIフィルタとして機能するようになる。
Figure 0004127827
なお、係数パラメータβは、速度検出の対象となる運動体に応じて可変に設定される。制御部12は、操作者による入力部2の操作に応じて係数パラメータβを可変に設定する。たとえば、係数パラメータβは、速度検出の対象となる運動体が血液等の比較的高速度にて動く運動体である場合、0.1〜0.2程度の値に設定されることが望ましい。
その後、制御部12は、上述した関心速度レンジ±Viおよび検出可能速度レンジ±Vrを示す各パラメータ信号をカラー画像データ生成部8に送信することによって、カラー画像データ生成部8に対し、上述したカラースケールデータと速度との対応付けを指示する(ステップS106)。この場合、カラー画像データ生成部8は、上述したように、保持するカラースケールデータとこれらの各パラメータ信号に基づく関心速度レンジ±Viおよび検出可能速度レンジ±Vrとを用い、関心速度レンジ±Vi内の各速度にこのカラースケールデータの各輝度または各色相をそれぞれ割り当てる。これにと同時に、カラー画像データ生成部8は、関心速度レンジ±Vi外の速度であって検出可能速度レンジ±Vr内の各速度にこのカラースケールデータの残りの各輝度または各色相をそれぞれ割り当てる。これによって、カラー画像データ生成部8は、このカラースケールデータと速度との対応付けを達成する。
上述したカットオフ周波数fcを示すパラメータ信号がフィルタ72に送信されかつ上述した関心速度レンジ±Viおよび検出可能速度レンジ±Vrを示す各パラメータ信号がカラー画像データ生成部8に送信された場合、制御部12は、ステップS105によるフィルタ72のカットオフ周波数fcの設定とステップS106によるカラー画像データ生成部8の速度と対応付けられたカラースケールデータの設定とがともに完了したか否かを確認する。この場合、制御部12は、フィルタ72がカットオフ周波数fcを設定完了したことを示す応答信号に基づき、このカットオフ周波数fcの設定完了を確認し、カラー画像データ生成部8がカラースケールデータと速度との対応付けを完了したことを示す応答信号に基づき、この速度と対応付けられたカラースケールデータの設定完了を確認する。したがって、制御部12は、このフィルタ72からの応答信号またはこのカラー画像データ生成部8からの応答信号を受信しなければ、このカットオフ周波数fcの設定完了とこの速度と対応付けられたカラースケールデータの設定完了とをともに検知せず(ステップS107,No)、このステップS107の処理手順を繰り返す。
一方、制御部12は、制御部12は、このフィルタ72からの応答信号とこのカラー画像データ生成部8からの応答信号とを受信した場合、このカットオフ周波数fcの設定完了とこの速度と対応付けられたカラースケールデータの設定完了とをともに検知し(ステップS107,Yes)、上述したステップS103によって設定した実繰り返し周波数frを示す制御信号を送受信回路4に送信することによって、送受信回路4に上述した電気的パルス信号の送受信を指示する(ステップS108)。これによって、送受信回路4は、上述したように、この実繰り返し周波数frに基づく回数の電気的パルス信号を繰り返し送受信する。
つぎに、制御部12は、画像合成部9に対し、速度データ演算部7が算出した運動体の速度Vを示すカラードプラ画像と被検体内のBモード画像とを少なくとも含むモニタ画像をモニタ10に表示するように指示する(ステップS109)。この場合、画像合成部9は、制御部12の制御のもと、上述した合成画像データを生成するとともに、得られた合成画像データを上述した表示画像データに変換し、さらに、この表示画像データをモニタ10に送信する。モニタ10は、この画像合成部9から受信した表示画像データをもとに、図4に例示するモニタ画像100を表示出力する。たとえば、モニタ10は、図4に示すように、被検体内を示すBモード画像101と、この被検体内の所望領域すなわち関心領域にて動く運動体のカラードプラ画像102と、Bモード画像101のグレースケール103と、カラードプラ画像102のカラースケール104とを含むモニタ画像100を表示出力する。この場合、操作者は、カラードプラ画像102とカラースケール104とを参照することによって、被検体内を比較的高速度にて動く所望の運動体の速度たとえば血液等の流速およびその向きを把握できる。なお、操作者は、入力部2を操作することによって、Bモード画像101上のカラードプラ画像102の表示領域を所望の領域に設定できる。
その後、操作者が入力部2を用いて終了指示情報または関心速度レンジ指示情報の入力操作を行わなければ、制御部12は、これらの指示情報を検知せず(ステップS110,No)、上述したステップS108以降の処理手順を繰り返す。なお、この終了指示情報は、この所望の運動体の速度検出を終了指示するための指示情報であり、たとえば送受信回路4に対して上述した電気的パルス信号の送受信終了を指示する指示情報である。
一方、制御部12は、入力部2から入力された指示情報を検知し(ステップS110,Yes)、さらにこの検知した指示情報が上述した関心速度レンジ指示情報である場合(ステップS111、関心速度レンジ指示情報)、上述したステップS102以降の処理手順を繰り返す。また、制御部12は、入力部2から入力された指示情報を検知し(ステップS110,Yes)、この検知した指示情報が上述した終了指示情報である場合(ステップS111、終了指示情報)、たとえば送受信回路4に対して上述した電気的パルス信号の送受信終了を指示し、この所望の運動体の速度検出に関する各種処理を終了させる。
なお、制御部12は、ティシュードプラ画像モードにおいて上述したステップS101〜ステップS111の各処理手順を行うことによって、被検体内の運動体の速度を示す速度画像すなわちティシュードプラ画像をモニタ10に表示出力させることができる。この場合、制御部12は、上述したステップS105の処理手順に代えて、フィルタ72に所定のフィルタ係数を設定させるとともにフィルタ72をローパスフィルタとして機能させる処理手順またはフィルタ72のフィルタとしての機能を停止させる処理手順を行う。これによって、モニタ10は、図4に示したモニタ画像100のカラードプラ画像102に代えてティシュードプラ画像を表示出力し、かつカラースケール104に代えてこのティシュードプラ画像のカラースケールを表示出力する。操作者は、このティシュードプラ画像とこのカラースケールとを参照することによって、被検体内を比較的低速度にて動く所望の運動体の速度たとえば生体組織等の運動速度を把握できる。
つぎに、制御部12が上述したステップS103の実繰り返し周波数設定処理を達成するまでの処理手順について詳細に説明する。図5は、ステップS103の実繰り返し周波数設定処理を達成するまでの処理手順を例示するフローチャートである。図5において、制御部12は、上述したステップS102において基準繰り返し周波数fiを算出した場合、上述した式(10)と同様に、この得られた基準繰り返し周波数fiと可変の係数パラメータαの最大値である係数パラメータαmaxとを乗算することによって、仮の実繰り返し周波数fr'を算出する(ステップS201)。
その後、制御部12は、得られた仮の実繰り返し周波数fr'を徐々に基準繰り返し周波数fiに近づけるとともに、この得られた仮の実繰り返し周波数fr'を示す制御信号を送受信回路4に送信することによって、上述した送受信回路4による電気的パルス信号の送受信を制御する周波数スイープ処理を行う(ステップS202)。この場合、制御部12は、この基準繰り返し周波数fiに乗じた係数パラメータαを最大値(すなわちα=αmax)から所定の数値間隔にてα=1に近づけるように変化させ、仮の実繰り返し周波数fr'を順次変化させる。
さらに、制御部12は、このステップS202の周波数スイープ処理を行う毎に、上述したステップS104とほぼ同様に、仮の実繰り返し周波数fr'に対応する検出可能速度レンジの最大値である仮の最大速度Vr'を算出し、この得られた仮の最大速度Vr'に基づく仮の検出可能速度レンジ±Vr'を設定する(ステップS203)。
なお、制御部12は、カラードプラ画像モードであれば、仮の検出可能速度レンジ±Vr'を設定する毎に、上述したステップS105とほぼ同様の処理手順を行って、フィルタ72にカットオフ周波数fcを仮に設定させてもよい。また、制御部12は、ティシュードプラ画像モードであれば、仮の検出可能速度レンジ±Vr'を設定した場合にフィルタ72に対し、フィルタとしての機能を停止させるように制御する。
一方、制御部12が仮の実繰り返し周波数fr'を示す制御信号を送受信回路4に送信した場合、送受信回路4は、制御部12の制御のもと、この仮の実繰り返し周波数fr'に基づく回数の電気的パルス信号を繰り返し送受信する。この場合、運動情報演算部74は、上述したように、この仮の実繰り返し周波数fr'に基づく回数の電気的パルス信号を繰り返し送受信することによって得られる超音波データ群に基づく被検体内の運動体の速度Vを算出する。制御部12は、運動情報演算部74に対してこの算出した速度Vをフィードバックするように制御し、運動情報演算部74からこの算出した速度Vを検出する(ステップS204)。
その後、制御部12は、運動情報演算部74から検出した速度Vを用い、ステップS203において設定した仮の検出可能速度レンジ±Vr'において上述したエリアジングが発生したか否かを判定する。この場合、制御部12は、この速度Vに符号の反転が発生したか否かによって、エリアジングが発生したか否かを判定する。ここで、係数パラメータαがステップS202の周波数スイープ処理において最大値(=αmax)の近傍である場合、この係数パラメータαに基づく仮の実繰り返し周波数fr'が基準繰り返し周波数fiよりも十分に大きいので、この仮の実繰り返し周波数fr'に基づく仮の検出可能速度レンジ±Vr'は、関心速度レンジ±Viに比して十分大きい速度レンジすなわち速度幅を有する。この場合、関心速度レンジ±Vi内の速度として推定される運動体の速度は、この仮の検出可能速度レンジ±Vr'の範囲内であると考えられる。したがって、制御部12は、サンプリングの定理に基づき、この場合に運動情報演算部74が算出した速度Vを正しい符号の速度として検出できる。
制御部12は、この正しい符号の速度として検出した速度Vについて、ステップS202の周波数スイープ処理毎にその符号の反転が生じたか否かを確認する。この速度Vの符号の反転が確認されていない場合、制御部12は、この仮の検出可能速度レンジ±Vr'においてエリアジング発生を検知せず(ステップS205,No)、上述したステップS202以降の処理手順を繰り返す。一方、この速度Vの符号の反転が確認された場合、制御部12は、この仮の検出可能速度レンジ±Vr'においてエリアジング発生を検知し(ステップS205,Yes)、このエリアジング発生を検知していない最後の周波数スイープ処理によって設定した係数パラメータαと基準繰り返し周波数fiとを乗算した仮の実繰り返し周波数fr'を上述した実繰り返し周波数frとして設定する(ステップS206)。
なお、制御部12は、上述したステップS201に代えて、上述した式(10)と同様に、基準繰り返し周波数fiと可変の係数パラメータαの最小値(すなわち1)とを乗算することによって仮の実繰り返し周波数fr'を算出し、さらに、上述したステップS202に代えて、この仮の実繰り返し周波数fr'を徐々に大きくするとともに、この仮の実繰り返し周波数fr'を示す制御信号を送受信回路4に送信することによって、上述した送受信回路4による電気的パルス信号の送受信を制御する周波数スイープ処理を行ってもよい。すなわち、制御部12は、この基準繰り返し周波数fiに乗じた係数パラメータαを最小値(すなわちα=1)から所定の数値間隔にて徐々に大きくなるように変化させ、仮の実繰り返し周波数fr'を順次変化させてもよい。
この場合、制御部12は、係数パラメータαが最小値近傍であれば、反転した符号である場合が多い速度Vを運動情報演算部74から検出する。したがって、制御部12は、上述したステップS205に代えて、この運動情報演算部74から検出した速度Vの符号の反転を確認した場合すなわちエリアジング発生を検知した場合に、上述した周波数スイープ処理以降の処理手順を繰り返し、この速度Vの符号の反転を確認しなくなった場合すなわちエリアジング発生を検知しなくなった場合に、実繰り返し周波数frを設定する。また、制御部12は、この実繰り返し周波数frを設定する場合、上述したステップS206に代えて、このエリアジング発生を検知しなくなった最初の周波数スイープ処理によって設定した係数パラメータαと基準繰り返し周波数fiとを乗算した仮の実繰り返し周波数fr'を上述した実繰り返し周波数frとして設定する。
ここで、制御部12は、上述したように、入力部2から入力された関心速度レンジ指示情報に基づき関心速度レンジ±Viを一意的に設定するとともに、この設定した関心速度レンジ±Viについて一意的に基準繰り返し周波数fiを取得する。また、制御部12は、上述した周波数スイープ処理を行って仮の実繰り返し周波数fr'を徐々に変化させ、エリアジング発生を検知しない最後の周波数スイープ処理による係数パラメータα言い換えればエリアジング発生を検知しなくなった最初の周波数スイープ処理による係数パラメータαを検出するとともに、この基準繰り返し周波数fiとこの係数パラメータαとを乗算することによって実繰り返し周波数frを取得している。したがって、制御部12は、この実繰り返し周波数frをもとに最大速度Vrを算出することによって、関心速度レンジ±Viを変化させずに、関心速度レンジ±Viに比して適度に大きい検出可能速度レンジ±Vrを設定することができる。この場合、検出可能速度レンジ±Vrは、関心速度レンジ±Viに対して過度に大きくなく、この関心速度レンジ±Vi内と推定される速度がオーバーレンジしない程度に十分大きい速度幅を有する。したがって、カラードプラ画像、ティシュードプラ画像、またはカラースケール等の各画像表示に支障を来たすことなく、運動体の速度を示すカラードプラ画像またはティシュードプラ画像を表示出力する場合のエリアジング発生を抑制できる。
なお、制御部12は、係数パラメータαを2から4までの範囲の実数として設定することが望ましく、これによって、上述した適度な大きさの速度幅を有する検出可能速度レンジを設定することができる。この場合、制御部12は、上述した周波数スイープ処理を行う場合、係数パラメータαを1から5程度の範囲において変化させることが望ましい。
また、カットオフ周波数fcは、式(13)に示すように、基準繰り返し周波数fiと係数パラメータα,βとを用いて表される。すなわち、制御部12は、フィルタ72にカットオフ周波数fcを設定させることによって、検出可能速度レンジ±Vrのノイズとして除去する速度範囲として、関心速度レンジ±Viのノイズとして除去する本来の速度範囲を設定することができる。これによって、制御部12は、関心速度レンジ±Viの本来検出すべき低速度範囲の速度の検出能力を損なうことなく、関心速度レンジ±Vi以上に広い速度幅を有する検出可能速度レンジ±Vrを設定できる。
つぎに、カラー画像データ生成部8が上述したカラースケールデータと速度とを対応付ける処理について詳細に説明する。図6は、検出可能速度レンジ±Vr内の各速度と対応付けられたカラースケールデータの一例を模式的に例示する模式図である。図7は、速度変化に対するカラースケールデータの輝度変化の一例を示す模式図である。図8は、速度変化に対するカラースケールデータの別の輝度変化を例示する模式図である。図9は、速度変化に対するカラースケールデータの色相変化の一例を例示する模式図である。
カラー画像データ生成部8は、上述したように、保持するカラースケールデータと制御部12からの各パラメータ信号に基づく関心速度レンジ±Viおよび検出可能速度レンジ±Vrとを用い、検出可能速度レンジ±Vr内の各速度すなわち関心速度レンジ±Vi内の各速度と関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の各速度とに対し、このカラースケールデータの各輝度または各色相をそれぞれ割り当てる。これによって、カラー画像データ生成部8は、たとえば図6に示すカラースケールデータ110を生成する。
カラースケールデータ110は、図6に示すように、運動体の正の速度すなわち検出可能速度レンジ±Vrの正の速度に対応するカラースケール要素110aと、運動体の負の速度すなわち検出可能速度レンジ±Vrの負の速度に対応するカラースケール要素110bとによって構成される。この場合、カラースケール要素110aは、0〜Viの速度範囲すなわち関心速度レンジ±Vi内の正の速度とVi〜Vrの速度範囲すなわち関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の正の速度とに対応する。また、カラースケール要素110bは、0〜−Viの速度範囲すなわち関心速度レンジ±Vi内の負の速度と−Vi〜−Vrの速度範囲すなわち関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の負の速度とに対応する。
たとえば、カラースケール要素110aは、速度0近傍すなわちノイズとして除去される速度範囲に黒色が割り当てられ、検出可能速度レンジ±Vr内の正の速度範囲に赤色から黄色に変化する過程の各色が割り当てられる。また、カラースケール要素110bは、速度0近傍すなわちノイズとして除去される速度範囲に黒色が割り当てられ、検出可能速度レンジ±Vr内の負の速度範囲に濃紫色から水色に変化する過程の各色が割り当てられる。
ここで、カラー画像データ生成部8は、検出可能速度レンジ±Vrの正の各速度に対してカラースケール要素110aの各輝度または各色相をそれぞれ割り当てる場合、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の正の速度の速度変化に対する輝度変化または色相変化よりも、関心速度レンジ±Vi内の正の速度の速度変化に対する輝度変化または色相変化を大きくする。
たとえば、カラー画像データ生成部8は、図7に示すように、速度0から速度Viまで変化する速度範囲において、カラースケール要素110aの輝度Lを0から輝度Liまでリニアに単調増加させ、速度Viから速度Vrまで変化する速度範囲において、輝度Lを輝度Liから輝度Lrまでリニアに単調増加させる。この場合、カラー画像データ生成部8は、速度Viを境界として、速度Viから速度Vrまで変化する速度範囲の速度変化に対する輝度変化よりも、速度0から速度Viまで変化する速度範囲の速度変化に対する輝度変化を大きくする。これによって、カラー画像データ生成部8は、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の正の各速度にそれぞれ割り当てる各輝度の隣接する速度間の輝度変化を小さくするとともに、関心速度レンジ±Vi内の正の各速度にそれぞれ割り当てる各輝度の隣接する速度間の輝度変化を大きくすることができる。
また、カラー画像データ生成部8は、図8に示すように、速度0から速度Viまで変化する速度範囲において、カラースケール要素110aの輝度Lを0から輝度Liまで曲線的に単調増加させ、速度Viから速度Vrまで変化する速度範囲において、輝度Lを輝度Liから輝度Lrまで曲線的に単調増加させ、さらに、速度Viを境界として、速度Viから速度Vrまで変化する速度範囲の速度変化に対する輝度変化よりも、速度0から速度Viまで変化する速度範囲の速度変化に対する輝度変化を大きくしてもよい。この場合も、カラー画像データ生成部8は、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の正の各速度にそれぞれ割り当てる各輝度の隣接する速度間の輝度変化を小さくするとともに、関心速度レンジ±Vi内の正の各速度にそれぞれ割り当てる各輝度の隣接する速度間の輝度変化を大きくすることができる。
さらに、カラー画像データ生成部8は、図9に示すように、速度0から速度Viまで変化する速度範囲において、カラースケール要素110aの色相CLを色相CL1から色相CL2までリニアに単調変化させ、速度Viから速度Vrまで変化する速度範囲において、色相CLを色相CL2から色相CL3までリニアに単調変化させ、さらに、速度Viを境界として、速度Viから速度Vrまで変化する速度範囲の速度変化に対する色相変化よりも、速度0から速度Viまで変化する速度範囲の速度変化に対する色相変化を大きくしてもよい。この場合、カラー画像データ生成部8は、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の正の各速度にそれぞれ割り当てる各色相の隣接する速度間の色相変化を小さくするとともに、関心速度レンジ±Vi内の正の各速度にそれぞれ割り当てる各色相の隣接する速度間の色相変化を大きくすることができる。
なお、カラースケール要素110bは、カラースケール要素110aと対応付けられる各速度の符号を反転したデータである。したがって、カラー画像データ生成部8は、カラースケール要素110bについてもカラースケール要素110aの場合とほぼ同様に、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の負の各速度にそれぞれ割り当てる各輝度または各色相の隣接する速度間の輝度変化または色相変化を小さくするとともに、関心速度レンジ±Vi内の負の各速度にそれぞれ割り当てる各輝度または各色相の隣接する速度間の輝度変化または色相変化を大きくすることができる。
カラー画像データ生成部8は、かかるカラースケール要素110aとカラースケール要素110bとによって構成したカラースケールデータ110を用いることによって、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の各速度に比較的緩やかな変化の各輝度または各色相を対応付けるとともに、関心速度レンジ±Vi内の各速度にダイナミックに変化する各輝度または各色相を対応付けることができる。これによって、カラー画像データ生成部8は、モニタ10に表示出力された場合に操作者が運動体の関心速度を認識し易い速度画像に対応するカラー画像データを生成できる。
なお、この発明の実施の形態1では、検出可能速度レンジ±Vrの全速度範囲すなわち関心速度レンジ±Viの速度範囲と関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の速度範囲とにおいて速度変化のスケールを同じにした場合を例示したが、この発明はこれに限定されるものではなく、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の速度範囲内の速度変化のスケールよりも、関心速度レンジ±Viの速度範囲内の速度変化のスケールを大きくしてもよい。
図10は、関心速度レンジ±Viの速度範囲内の速度変化のスケールを大きくしたカラースケールデータの一例を模式的に例示する模式図である。このカラースケールデータ120は、図10に示すように、関心速度レンジ±Vi内の正の速度に対応するカラースケール要素120aと、関心速度レンジ±Vi内の負の速度に対応するカラースケール要素120bと、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の正の速度に対応するカラースケール要素120cと、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の負の速度に対応するカラースケール要素120dとによって構成される。
カラー画像データ生成部8は、検出可能速度レンジ±Vrの各速度にカラースケールデータ120の各輝度または各色相をそれぞれ割り当てる場合、カラースケール要素120c,120dと対応付ける速度の速度変化のスケールを縮小するとともに、カラースケール要素120a,120bと対応付ける速度の速度変化のスケールを拡大する。これによって、カラー画像データ生成部8は、図10に示すように、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の速度に対応する部分に比して、関心速度レンジ±Vi内の速度に対応する部分の速度幅を十分大きくしたカラースケールデータ120を生成できる。すなわち、カラー画像データ生成部8は、関心速度レンジ±Viに比してさらに大きい検出可能速度レンジ±Vrが設定された場合であっても、関心速度レンジ±Vi内の速度に対応する部分の占める割合が他の部分よりも十分高いカラースケールデータを確実に生成できる。
この場合、モニタ10には、カラースケールデータ120に例示される関心速度レンジ±Vi内の速度に対応する部分の占める割合が他の部分よりも十分高いカラースケールデータを示すカラースケールの画像を表示出力できる。操作者は、この表示出力されたカラースケールを参照することによって、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の速度を認識できるとともに、関心速度レンジ±Vi内の速度を容易にかつ確実に認識できる。
なお、カラー画像データ生成部8は、最大速度Viを境界にしてカラースケール要素120aとカラースケール要素120cとを分割してもよい。これと同様に、カラー画像データ生成部8は、最小速度−Viを境界にしてカラースケール要素120bとカラースケール要素120dとを分割してもよい。
また、この発明の実施の形態1では、アレイ振動子を用いて実現された超音波振動子3が設けられた場合を例示したが、この発明はこれに限定されるものではなく、回転駆動系が設けられ、機械的な駆動を伴った超音波スキャンを行う超音波振動子を超音波振動子3に代えて用いてもよい。
さらに、この発明の実施の形態1では、演算パラメータ等の各種情報が各構成部にそれぞれ格納されていたが、この発明はこれに限定されるものではなく、制御部12がこれらの各種情報を一括して記憶し管理してもよい。
また、この発明の実施の形態1では、制御部12が、カットオフ周波数fcを算出するとともに、このカットオフ周波数fcを示すパラメータ信号をフィルタ72に送信していたが、この発明はこれに限定されるものではなく、制御部12が基準繰り返し周波数fiおよび最大速度Vi,Vrを示す各パラメータ信号をフィルタ72に送信し、フィルタ72が、制御部12から受信した各パラメータ信号に基づき、カットオフ周波数fcを算出してもよい。
さらに、この発明の実施の形態1では、制御部12は、フィルタ72にカットオフ周波数fcを設定指示し、その後、カラー画像データ生成部8にカラースケールデータと速度との対応付けを指示していたが、この発明はこれに限定されるものではなく、制御部12は、カラー画像データ生成部8にカラースケールデータと速度との対応付けを指示し、その後またはこれと同時に、フィルタ72にカットオフ周波数fcを設定指示してもよい。
また、この発明の実施の形態1では、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の速度変化に応じてカラースケールデータの輝度または色相を変化させていたが、この発明はこれに限定されるものではなく、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の速度変化に対してカラースケールデータの輝度または色相を一定にしてもよい。
さらに、この発明の実施の形態1では、検出可能速度レンジ±Vr内の速度変化に対して色相をリニアに変化させていたが、この発明はこれに限定されるものではなく、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の速度変化に対する色相変化に比して関心速度レンジ±Vi内の速度変化に対する色相変化が大きければ、検出可能速度レンジ±Vr内の速度変化に対して色相を曲線的に変化させてもよい。
以上に説明したように、この発明の実施の形態1では、操作者の入力操作によって指示入力された関心速度レンジ指示情報をもとに関心速度レンジを一意的に設定するとともに、この関心速度レンジ以上に広い速度レンジであってこの関心速度レンジを含む可変の検出可能速度レンジを設定し、さらに、この関心速度レンジに応じた速度零近傍の速度レンジを除くこの検出可能速度レンジ内の速度を被検体内にて動く所望の運動体の速度として算出するように構成したので、この関心速度レンジ内の所望の低速度範囲の速度を検出する能力を損なうことなく、この所望の運動体の速度を検出する場合のエリアジングの発生を抑制できる超音波診断装置を実現することができる。
また、この検出可能速度レンジ内の各速度に所望のカラースケールデータの各輝度または各色相をそれぞれ割り当て、この割り当てたカラースケールデータと所望の運動体の速度として算出した速度とをもとに、この所望の運動体の速度を示す速度画像を生成出力するように構成したので、エリアジングを発生させずに、この検出可能速度レンジ内の速度を示す速度画像を確実に画面表示できる超音波診断装置を実現することができる。
さらに、この関心速度レンジ外であってこの検出可能速度レンジ内の各速度にそれぞれ割り当てる各輝度または各色相よりも、この関心速度レンジ内の各速度にそれぞれ割り当てる各輝度または各色相の隣接する速度間の輝度変化または色相変化を大きくしているので、この関心速度レンジ内の速度を容易に認識できる速度画像を画面表示することができる。
また、関心速度レンジ外であって検出可能速度レンジ内の速度変化のスケールを縮小するとともに、関心速度レンジ内の速度変化のスケールを拡大することによって、関心速度レンジ外であって検出可能速度レンジ内の速度に対応するカラースケールデータに比して、関心速度レンジ内の速度に対応するカラースケールデータの速度幅を十分大きくすることができ、関心速度レンジ内の速度に対応するカラースケールデータの占める割合が他の部分よりも十分高いカラースケールの画像を画面表示できる。操作者は、このカラースケールの画像を参照することによって、関心速度レンジ外であって検出可能速度レンジ内の速度を認識できるとともに、関心速度レンジ内の速度を容易にかつ確実に認識できる。
(実施の形態2)
つぎに、この発明の実施の形態2について詳細に説明する。上述した実施の形態1では、検出可能速度レンジ±Vr内の全ての速度にカラースケールデータの各輝度または各色相をそれぞれ割り当て、この検出可能速度レンジ±Vr内の全ての速度に対応するカラースケールの画像を表示出力していたが、この実施の形態2では、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の各速度にカラースケールデータの各輝度および各色相を割り当てずに、関心速度レンジ±Vi内の各速度にカラースケールデータの各輝度または各色相をそれぞれ割り当て、この関心速度レンジ±Vi内の全ての速度に対応するカラースケールの画像を表示出力するように構成している。
図11は、この発明の実施の形態2である超音波診断装置の一構成例を例示するブロック図である。この超音波診断装置21は、カラー画像データ生成部8に代えてカラー画像データ生成部22が設けられる。その他の構成は実施の形態1と同じであり、同一構成部分には同一符号を付している。
図12は、関心速度レンジ±Vi内の各速度と対応付けられたカラースケールデータの一例を模式的に例示する模式図である。カラー画像データ生成部22は、上述したカラー画像データ生成部8とほぼ同様の機能および構成を有する。また、カラー画像データ生成部22は、関心速度レンジ±Viおよび検出可能速度レンジ±Vrを示す各パラメータ信号を制御部12から受信した場合、制御部12の制御のもと、保持するカラースケールデータとこのパラメータ信号に基づき、関心速度レンジ±Vi内の各速度にこのカラースケールデータの各輝度または各色相をそれぞれ割り当てる。これによって、カラー画像データ生成部22は、たとえば図12に示すカラースケールデータ130を生成する。
カラースケールデータ130は、図12に示すように、0〜Viの速度範囲すなわち関心速度レンジ±Viの正の速度に対応するカラースケール要素130aと、0〜−Viの速度範囲すなわち関心速度レンジ±Virの負の速度に対応するカラースケール要素110bとによって構成される。たとえば、カラースケール要素130aは、速度0近傍すなわちノイズとして除去される速度範囲に黒色が割り当てられ、関心速度レンジ±Vi内の正の速度範囲に赤色から黄色に変化する過程の各色が割り当てられる。また、カラースケール要素130bは、速度0近傍すなわちノイズとして除去される速度範囲に黒色が割り当てられ、関心速度レンジ±Vi内の負の速度範囲に濃紫色から水色に変化する過程の各色が割り当てられる。この場合、カラー画像データ生成部22は、関心速度レンジ±Vi内に対し、保持するカラースケールデータのほぼ全ての輝度またはほぼ全ての色相を割り当てる。
また、カラー画像データ生成部22は、上述した関心速度データを速度データ演算部7から受信した場合、この関心速度データに基づく運動体の速度を関心速度レンジ±Vi内の速度と関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の速度とに分類する。さらに、カラー画像データ生成部22は、この関心速度データに基づく運動体の速度を関心速度レンジ±Vi内の速度として取得した場合、この関心速度レンジ±Vi内の各速度と対応付けたカラースケールデータとこの関心速度データとを用い、この取得した速度に対応するカラー画像データを取得する。この場合、カラー画像データ生成部22は、この取得したカラー画像データを画像合成部9に送信する。これによって、モニタ10は、この取得した速度を示す速度画像すなわちカラードプラ画像またはティシュードプラ画像を表示出力できる。さらに、カラー画像データ生成部22がカラースケールデータ130に例示されるこの関心速度レンジ±Vi内の各速度と対応付けたカラースケールデータに対応する画像データを画像合成部9に送信することによって、モニタ10は、この取得した速度を示すカラードプラ画像またはティシュードプラ画像と同一画面上に、このカラースケールデータを示すカラースケールの画像を表示出力できる。
一方、カラー画像データ生成部22は、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の各速度に所定色たとえば黒色を割り当てる。したがって、カラー画像データ生成部22は、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の速度として、速度データ演算部7から受信した関心速度データに基づく運動体の速度を取得した場合、この関心速度データを変換し、この取得した速度に黒色を割り当てた画像データを取得する。この場合、カラー画像データ生成部22は、この運動体の速度を関心速度レンジ±Vi内の速度と関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の速度とに分類するとともに、上述した関心速度レンジ±Vi内の各速度と対応付けたカラースケールデータを用いず、この関心速度データに基づく速度に黒色を割り当てる。したがって、カラー画像データ生成部22は、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の速度を速度データ演算部7から取得した場合であっても上述したエリアジングを発生させない。
また、カラー画像データ生成部22は、この取得した画像データを上述した目的外カラー画像データとして画像合成部9に送信する。画像合成部9は、上述したように、この目的外カラー画像データに上述したグレー画像データを上書きすることによって合成画像データを取得する。この場合、モニタ10は、関心速度レンジ±Vi外であって検出可能速度レンジ±Vr内の運動体の速度を示すカラードプラ画像またはティシュードプラ画像を表示しない。
以上に説明したように、この発明の実施の形態2では、上述した実施の形態1とほぼ同様の機能および構成を有し、また、設定した関心速度レンジ内の速度に対してカラースケールデータの全ての輝度または全ての色相を割り当てるとともに、この関心速度レンジ外であって設定した検出可能速度レンジ内の各速度に所定色たとえば黒色を割り当てるようにし、さらに、算出した運動体の速度をこの関心速度レンジ内の速度とこの関心速度レンジ外であってこの検出可能速度レンジ内の速度とに分類し、この関心速度レンジ内の速度として分類した運動体の速度を示す速度画像を生成出力し、かつこの関心速度レンジ外であってこの検出可能速度レンジ内の速度として分類した運動体の速度を速度画像として表示出力しないように構成している。したがって、この算出した運動体の速度についてエリアジングの発生を抑制できるとともに、この関心速度レンジ外であってこの検出可能速度レンジ内の運動体の速度を示す速度画像を画面表示せずに、この関心速度レンジ内の運動体の速度を示す速度画像を画面表示でき、上述した実施の形態1とほぼ同様の作用効果を享受するとともに、この関心速度レンジ内である所望の運動体の速度を容易に認識できる超音波診断装置を実現することができる。
この発明の実施の形態1である超音波診断装置の一構成例を例示するブロック図である。 速度データ演算部の一構成を詳細に例示するブロック図である。 運動体の速度画像をモニタ表示するまでの処理手順を例示するフローチャートである。 被検体内のBモード画像とカラードプラ画像とを含むモニタ画像の一表示例を模式的に例示する模式図である。 実繰り返し周波数設定処理を達成するまでの処理手順を例示するフローチャートである。 検出可能速度レンジ内の各速度と対応付けられたカラースケールデータの一例を模式的に例示する模式図である。 速度変化に対するカラースケールデータの輝度変化の一例を示す模式図である。 速度変化に対するカラースケールデータの別の輝度変化を例示する模式図である。 速度変化に対するカラースケールデータの色相変化の一例を例示する模式図である。 関心速度レンジ内の速度変化のスケールを大きくしたカラースケールデータの一例を模式的に例示する模式図である。 この発明の実施の形態2である超音波診断装置の一構成例を例示するブロック図である。 関心速度レンジ内の各速度と対応付けられたカラースケールデータの一例を模式的に例示する模式図である。
符号の説明
1,21 超音波診断装置
2 入力部
3 超音波振動子
4 送受信回路
5 Bモードデータ演算部
6 グレー画像データ生成部
7 速度データ演算部
8,22 カラー画像データ生成部
9 画像合成部
10 モニタ
11 記憶部
12 制御部
71 複素信号化回路
72 フィルタ
73 自己相関回路
74 運動情報演算部
75 しきい値処理回路
100 モニタ画像
101 Bモード画像
102 カラードプラ画像
103 グレースケール
104 カラースケール
110,120,130 カラースケールデータ
110a,110b,120a〜120d,130a,130b カラースケール要素

Claims (6)

  1. 被検体内に対して超音波の送受信を複数回繰り返し行って得た複数の超音波データをもとに、前記被検体内の超音波断層像を生成出力するとともに、前記被検体内にて動く運動体の速度を所定速度レンジ内の速度として算出し、該算出した速度とカラースケールデータとをもとに前記運動体の速度を示す速度画像を生成出力する超音波診断装置において、
    前記運動体の関心速度レンジを示す情報を指示入力する入力手段と、
    前記入力手段から入力された情報によって指示される関心速度レンジに比して広い速度レンジであって該関心速度レンジを含む可変の検出可能速度レンジを前記所定速度レンジとして設定し、この設定した前記検出可能速度レンジの零近傍に、前記関心速度レンジに応じた零近傍の速度レンジ部分を除去する速度レンジを設定する速度レンジ設定制御手段と、
    を備えたことを特徴とする超音波診断装置。
  2. 前記検出可能速度レンジ内の速度毎に前記カラースケールデータを割り当て、該割り当てたカラースケールデータと前記算出した速度とをもとに前記速度画像を生成する画像処理制御手段を備えたことを特徴とする請求項1に記載の超音波診断装置。
  3. 前記検出可能速度レンジ内の速度毎に割り当てられたカラースケールデータに対応するカラースケール画像と前記速度画像とを含む各種画像を複数同時に表示出力する表示手段を備え、前記画像処理制御手段は、前記表示手段に対して前記各種画像の表示制御を行うことを特徴とする請求項2に記載の超音波診断装置。
  4. 前記画像処理制御手段は、前記表示手段に対し、前記関心速度レンジ内の前記カラースケール画像よりも前記関心速度レンジ外であって前記検出可能速度レンジ内の前記カラースケール画像を小さくする表示制御を行うことを特徴とする請求項3に記載の超音波診断装置。
  5. 前記画像処理制御手段は、前記検出可能速度レンジ内の速度毎に前記カラースケールデータを割り当てる場合、前記関心速度レンジ外であって前記検出可能速度レンジ内の速度変化に対する前記カラースケールデータの色相変化または輝度変化よりも、前記関心速度レンジ内の速度変化に対する前記カラースケールデータの色相変化または輝度変化を大きくすることを特徴とする請求項2〜4のいずれか一つに記載の超音波診断装置。
  6. 前記速度レンジ設定制御手段は、
    前記入力手段から入力された情報をもとに前記関心速度レンジを設定するとともに該関心速度レンジの最大速度値に対応する基準繰り返し周波数と前記超音波の送受信の繰り返し回数に関する仮の繰り返し周波数とを算出し、該仮の繰り返し周波数をスイープする周波数スイープ処理を行うとともに該周波数スイープ処理毎の仮の繰り返し周波数を用いて前記超音波の送受信制御を順次行う送受信制御手段と、
    前記送受信制御手段による前記超音波の送受信制御によって得られた前記複数の超音波データをもとに前記運動体の速度を順次算出する速度演算制御手段と、
    を備え、前記送受信制御手段は、前記速度演算制御手段から前記運動体の速度を順次検出するとともに、該順次検出した運動体の速度をもとにエリアジング現象の発生有無を順次判定するエリアジング判定処理を行い、該エリアジング判定処理の結果をもとに前記超音波の送受信制御を行う実際の繰り返し周波数と前記検出可能速度レンジとを設定し、前記速度演算制御手段は、前記基準繰り返し周波数を少なくとも用い、前記検出可能速度レンジから前記関心速度レンジに応じた零近傍の速度レンジ部分を除去する速度レンジを設定することを特徴とする請求項1〜5のいずれか一つに記載の超音波診断装置。
JP2004194888A 2004-06-30 2004-06-30 超音波診断装置 Expired - Fee Related JP4127827B2 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2004194888A JP4127827B2 (ja) 2004-06-30 2004-06-30 超音波診断装置
EP05753303A EP1769746A4 (en) 2004-06-30 2005-06-27 ULTRASOUND DIAGNOSTIC TOOL
PCT/JP2005/011763 WO2006003866A1 (ja) 2004-06-30 2005-06-27 超音波診断装置
CNB200580021841XA CN100463658C (zh) 2004-06-30 2005-06-27 超声波诊断装置
US11/648,285 US20070167770A1 (en) 2004-06-30 2006-12-29 Ultrasonic diagnostic apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2004194888A JP4127827B2 (ja) 2004-06-30 2004-06-30 超音波診断装置

Publications (2)

Publication Number Publication Date
JP2006014891A JP2006014891A (ja) 2006-01-19
JP4127827B2 true JP4127827B2 (ja) 2008-07-30

Family

ID=35782675

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004194888A Expired - Fee Related JP4127827B2 (ja) 2004-06-30 2004-06-30 超音波診断装置

Country Status (5)

Country Link
US (1) US20070167770A1 (ja)
EP (1) EP1769746A4 (ja)
JP (1) JP4127827B2 (ja)
CN (1) CN100463658C (ja)
WO (1) WO2006003866A1 (ja)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101404940A (zh) * 2006-03-21 2009-04-08 皇家飞利浦电子股份有限公司 用于分析心功能不全的超声心动描记设备和方法
CN100448409C (zh) * 2007-02-01 2009-01-07 上海交通大学 三维超声心动图四腔切面图像自动检测的方法
DE202008003245U1 (de) 2007-02-21 2008-06-05 Universität Duisburg-Essen Vorrichtung zur Ultraschall-Messung von Blutfluß
JP4907382B2 (ja) * 2007-02-23 2012-03-28 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波画像表示方法および超音波診断装置
US7656750B2 (en) * 2008-02-26 2010-02-02 Semiconductor Components Industries, Llc Echo detection
JP5732461B2 (ja) * 2009-08-11 2015-06-10 日立アロカメディカル株式会社 超音波画像診断方法および超音波画像診断装置
KR101117879B1 (ko) * 2010-05-27 2012-03-07 삼성메디슨 주식회사 컬러 재구성 영상을 제공하는 초음파 시스템 및 방법
KR101820497B1 (ko) * 2011-07-27 2018-01-22 삼성디스플레이 주식회사 표시 장치 및 그 구동 방법
US9924142B2 (en) * 2012-11-21 2018-03-20 Omnivision Technologies, Inc. Camera array systems including at least one bayer type camera and associated methods
KR102243022B1 (ko) * 2014-03-05 2021-04-21 삼성메디슨 주식회사 선택 정보에 기초하여 관심 영역에 포함된 혈류에 대한 정보를 출력하는 방법, 장치 및 시스템.
JP6265933B2 (ja) * 2015-03-05 2018-01-24 富士フイルム株式会社 音響波診断装置およびその制御方法
CN106137257A (zh) * 2016-07-29 2016-11-23 飞依诺科技(苏州)有限公司 一种超声波彩色成像方法及装置
JP7157649B2 (ja) * 2018-12-19 2022-10-20 富士フイルムヘルスケア株式会社 超音波撮像装置およびその制御方法
JP7343342B2 (ja) 2019-09-25 2023-09-12 キヤノンメディカルシステムズ株式会社 超音波診断装置、及び画像処理装置
JP7555719B2 (ja) * 2020-03-31 2024-09-25 キヤノンメディカルシステムズ株式会社 超音波診断装置
JP7747548B2 (ja) * 2022-02-18 2025-10-01 キヤノンメディカルシステムズ株式会社 超音波診断装置及び表示制御方法

Family Cites Families (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63317142A (ja) * 1987-06-22 1988-12-26 Toshiba Corp 超音波診断装置
JPH03188841A (ja) * 1989-09-20 1991-08-16 Toshiba Corp 超音波診断装置
US5211169A (en) * 1990-11-08 1993-05-18 Prism Imaging, Inc. Blood pool imaging and analysis technique using ultrasound
JPH06186A (ja) * 1992-06-18 1994-01-11 Toshiba Corp 超音波血流イメージング装置
US5622174A (en) * 1992-10-02 1997-04-22 Kabushiki Kaisha Toshiba Ultrasonic diagnosis apparatus and image displaying system
JP3524129B2 (ja) * 1992-12-28 2004-05-10 株式会社東芝 超音波カラードプラ診断装置
US5441052A (en) * 1992-12-28 1995-08-15 Kabushiki Kaisha Toshiba Color doppler-type ultrasonic diagnostic apparatus
US5419331A (en) * 1994-02-10 1995-05-30 The University Of Rochester System for estimating target velocity from pulse echoes in response to their correspondence with predetermined delay trajectories corresponding to different distinct velocities
JP3833282B2 (ja) 1994-06-24 2006-10-11 株式会社東芝 超音波診断装置
JPH10314170A (ja) 1997-05-15 1998-12-02 Toshiba Corp 超音波診断装置
JP4077913B2 (ja) 1997-11-19 2008-04-23 東芝医用システムエンジニアリング株式会社 超音波診断装置
US5961462A (en) 1998-05-18 1999-10-05 Atl Ultrasound Ultrasonic doppler imaging at high frame rates of display
US6547736B1 (en) * 1998-11-11 2003-04-15 Spentech, Inc. Doppler ultrasound method and apparatus for monitoring blood flow and detecting emboli
US6176828B1 (en) * 1998-12-24 2001-01-23 General Electric Company Method and apparatus for optimal data mapping of power doppler images
US6162176A (en) * 1998-12-31 2000-12-19 General Electric Company Ultrasound color flow display optimization
JP3696763B2 (ja) * 1999-11-05 2005-09-21 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波撮影装置
US6508767B2 (en) * 2000-02-16 2003-01-21 Koninklijke Philips Electronics N.V. Ultrasonic harmonic image segmentation
US20030045797A1 (en) * 2001-08-28 2003-03-06 Donald Christopher Automatic optimization of doppler display parameters
US6663566B2 (en) * 2002-02-19 2003-12-16 Ge Medical Systems Global Technology Company, Llc Method and apparatus for automatic control of spectral doppler imaging
JP4253494B2 (ja) * 2002-11-08 2009-04-15 株式会社東芝 超音波診断装置
US6733454B1 (en) * 2003-02-26 2004-05-11 Siemens Medical Solutions Usa, Inc. Automatic optimization methods and systems for doppler ultrasound imaging

Also Published As

Publication number Publication date
CN1976637A (zh) 2007-06-06
US20070167770A1 (en) 2007-07-19
JP2006014891A (ja) 2006-01-19
EP1769746A1 (en) 2007-04-04
CN100463658C (zh) 2009-02-25
WO2006003866A1 (ja) 2006-01-12
EP1769746A4 (en) 2009-08-05

Similar Documents

Publication Publication Date Title
JP4127827B2 (ja) 超音波診断装置
JP3935433B2 (ja) 空間的複合のための適応画像処理
CN102283681B (zh) 超声波诊断装置及超声波诊断装置控制方法
CN101170947B (zh) 超声波诊断装置及超声波图像显示方法
US8206301B2 (en) Ultrasonic imaging apparatus and ultrasonic image processing apparatus, method and program
CN1182357A (zh) 超声波诊断装置
JP6345895B2 (ja) 超音波診断装置
WO2007114375A1 (ja) 超音波診断装置及び超音波診断装置制御方法
JP4928801B2 (ja) 超音波診断装置
JP2006187594A (ja) 超音波診断システム及び方法
JP2005253827A (ja) 超音波撮像方法及び装置
JP2000354595A (ja) カラードップラー映像システムにおけるカラー映像表示方法及び装置
KR20090023090A (ko) 초음파 진단 장치
JP4712931B2 (ja) パワー・ドップラー・イメージの最適データ・マッピングのための方法および装置
US20100137715A1 (en) Ultrasonic imaging apparatus and control method for ultrasonic imaging apparatus
JP4590256B2 (ja) 超音波撮像装置、超音波画像処理方法、及び、超音波画像処理プログラム
JP2004016612A (ja) 超音波診断装置
US8398548B2 (en) Ultrasound diagnostic apparatus and ultrasound diagnostic method
JP2000149015A (ja) 画像のエッジ強調方法及びイメ―ジング装置
JP4530834B2 (ja) 超音波画像処理方法及び超音波画像処理装置、並びに、超音波画像処理プログラム
JP7343342B2 (ja) 超音波診断装置、及び画像処理装置
JP4651379B2 (ja) 超音波画像処理装置及び超音波画像処理方法、並びに、超音波画像処理プログラム
JP4248673B2 (ja) 超音波診断装置
JP2006280520A (ja) 超音波撮像装置、超音波画像処理方法、及び、超音波画像処理プログラム
JP2006212054A (ja) 超音波観測装置、画像処理装置、及び、画像処理プログラム

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070911

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20071026

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: 20080507

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: 20080512

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110523

Year of fee payment: 3

R151 Written notification of patent or utility model registration

Ref document number: 4127827

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: 20110523

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120523

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130523

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140523

Year of fee payment: 6

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

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees