JP2004222754A - 超音波診断装置 - Google Patents
超音波診断装置 Download PDFInfo
- Publication number
- JP2004222754A JP2004222754A JP2003010531A JP2003010531A JP2004222754A JP 2004222754 A JP2004222754 A JP 2004222754A JP 2003010531 A JP2003010531 A JP 2003010531A JP 2003010531 A JP2003010531 A JP 2003010531A JP 2004222754 A JP2004222754 A JP 2004222754A
- Authority
- JP
- Japan
- Prior art keywords
- flow velocity
- time
- value
- ultrasonic
- maximum flow
- 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
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
【解決手段】超音波プローブ11で得た受信超音波信号を、受信プリアンプ21、A/D変換器22、ビームフォーマ23を経て検波回路24に送ってドプラ効果による周波数シフトを求め、さらにこれをレンジゲート25、クラッタ除去回路26を経てDSP27に送り、そのドプラ周波数シフトの周波数スペクトルを求めるFFTを行い、これを超音波信号の受信ごとに繰り返して順次得た周波数スペクトルを時間方向にトレースすることによって流速波形を求め、この流速波形をFIRフィルタで微分して極値を検知し、その極値が前に求めた最大流速の類似範囲に入っていることを判定して、前に求めた最大流速の時点から今回最大流速が得られた時点までの時間間隔を1心拍時間として得、この1心拍時間や最大流速などからPI、RIなどを自動的に演算する。
【選択図】 図1
Description
【発明の属する技術分野】
この発明は、医学的な診断に用いられる超音波診断装置に関し、とくに被検体内の流体運動を画像化する超音波診断装置に関する。
【0002】
【従来の技術】
超音波診断装置は、超音波ビームを被検体(被診察者の身体)内に入射し、その反射波を受波することによって身体内の断層像を得たり、ドプラ現象を利用して血流などの速度を表す画像を描出したりするものであり、医学的診断用途に広く普及している。超音波ドプラの基本的な原理は、流体からの反射波が受けたドプラ効果を超音波キャリア周波数の偏移として検出することにより流体の運動情報を得るというものである。
【0003】
具体的には、通常、超音波振動子を備えたプローブ(探触子)に、所定周波数の送信信号を送ってプローブ内の振動子を励振し、超音波を被検体内部に放射し、被検体内部で反射した超音波をそのプローブの振動子で受ける。プローブで検出した反射波信号は検波してキャリア周波数からの偏移成分のみ取り出し、それを周波数分析しスペクトルを得る。この周波数スペクトルは流体の速度に対応する。この周波数スペクトルを取得する超音波スキャンを繰り返し行い、このスペクトルを時間的にスクロールしながら表示することにより、血流などの流速の時間的変化を表す波形を得ることができる。
【0004】
臨床医療では、1心拍内の波形に対して、最大流速・最小流速・平均流速を求め、PI、RIと呼ばれる値を計算し、これらの値が適性範囲に入っているかを調べて疾患の診断を行う。
PI、RIは
PI=(収縮期最高血流速度−拡張終末期血流速度)/平均血流速度
RI=(収縮期最高血流速度−拡張終末期血流速度)/収縮期最高血流速度
と定義されており、これらの値を演算するためには収縮期最高血流速度、拡張終末期血流速度、平均血流速度が必要であるが、通常、拡張終末期血流速度=最小血流速度と近似している。
【0005】
【発明が解決しようとする課題】
しかしながら、従来の超音波診断装置では、PI、RIなどを求めることが簡単にできないという問題がある。これらの値を求めるためにはいったん超音波スキャンを停止し、表示されているスペクトル画像からデータを求めて演算する必要があり、リアルタイム性が損なわれる。加えて、ユーザは最大流速がどこで、最小流速がどこかを計算ソフトに指示したり、1心拍がどれだけの幅なのかを指示したりする必要もある。どの心拍期間に対して演算を行うかも指示する必要がある。
【0006】
この発明は、上記に鑑み、超音波スキャンを行っている最中にリアルタイムでPI、RIなどを自動的に求めることを可能とする超音波診断装置を提供することを目的とする。
【0007】
【課題を解決するための手段】
上記の目的を達成するため、この発明による超音波診断装置においては、超音波プローブと、該超音波プローブを介して超音波の送信および受信を行う超音波送受信手段と、受信した超音波反射波信号のドプラシフトを求める検波手段と、上記ドプラシフトの周波数スペクトルを得るフーリエ変換手段と、該スペクトルの時間トレースデータを時間方向に微分する手段と、微分値の符号が変化したときのデータを極値として求め、この極値により最大流速および最小流速を判定する手段とが備えられることが特徴となっている。
【0008】
超音波送受信ごとに受信信号のドプラシフトの周波数スペクトルがリアルタイムで得られ、その時間トレースデータが時間方向に微分される。最大流速や最小流速は、このスペクトルの時間トレースデータにおける極値として現れるので、上記の微分値の符号が変化する極値のデータのみに着目すれば足りる。このように極値のデータにのみ着目して最大流速、最小流速を判定することとすれば、演算が簡潔になる。判定された最大流速間の時間間隔により1心拍期間を求めることができ、その1心拍期間内での平均流速、最大流速、最小流速よりPI、RIなどを求めることができるので、超音波送受を連続して行いながらリアルタイムでのPI、RIなどの自動的な算出が可能となる。
【0009】
さらに、ドプラシフトの周波数スペクトルの時間トレースデータをスムージングする手段を備えるようにすれば、ノイズ等の影響を除去してより正確な測定ができる。
【0010】
このスムージングする手段はFIRフィルタで構成することができる。
【0011】
上記のドプラシフトの周波数スペクトルの時間トレースデータを時間方向に微分する手段は、FIRフィルタで構成することができる。
【0012】
スムージングする手段と微分する手段とをそれぞれFIRフィルタで構成した場合には、一つのFIRフィルタで両手段と兼用することが可能である。
【0013】
さらに、上記の、微分値の符号が変化したときのデータを極値として求め、この極値により最大流速および最小流速を判定する手段は、今回求めた極値が前回の極値として求めた最大流速の類似範囲に入っている条件で最大流速と判定して、前回の最大流速から今回の最大流速の間の時間間隔を1心拍時間と判定する機能をも備えるようにすることができる。
【0014】
【発明の実施の形態】
つぎに、この発明の実施の形態について図面を参照しながら詳細に説明する。
図1はこの発明の実施の形態を示す。この図1において、超音波プローブ11は被検体(被検者の腹部などの身体)にあてられるもので、超音波振動子を有し、超音波ビームの身体内への送信と身体内部からの反射波の受信を行う。送信ドライバ回路13がパルス制御回路12からのパルスに応じて所定周波数の駆動信号をパルス状に発生する。このパルス状の駆動信号が超音波プローブ11に送られ、その超音波振動子がパルス的に駆動されて励振させられ、所定周波数の超音波がパルス状に発生して、被検体内に放射される。
【0015】
被検体内で反射した超音波はプローブ11に戻ってきてその振動子を振動させて受信信号を出力させる。この受信信号は受信プリアンプ21を経てA/D変換器22に送られデジタル信号に変換される。この受信超音波反射波のデジタル信号はビームフォーマ23に送られ、鋭い指向性の受信超音波ビームの電子的な形成がなされることにより、そのビーム上の信号が取り出され、検波回路24に送られて、超音波キャリアの周波数で複素検波され、IとQの信号が得られる。このI,Qの検波信号はレンジゲート25を経て所定の時間範囲内のみ通過させられることにより、所定の反射深さの信号のみ取り出される。この受信超音波ビームの指向性と反射深さの限定により、所望の血管部分などの特定領域からの反射波のキャリア周波数からのドップラシフトを表す信号が取り出されることになる。
【0016】
このレンジゲート25を経た信号はクラッタ除去回路26に送られ、流体以外の信号が低減された後、DSP(Digtal Signal Processor)27に送られ、FFT(高速フーリエ変換)処理やFIR(Finite Impulse Response)フィルタ処理などを受ける。ビデオ表示装置28およびディスプレイ装置29は、FFT処理によって得られた周波数スペクトルを時間的にスクロールしながら表示することによって、図3に示すような流速の時間的な変化を表す波形を表示する。さらに、DSP27での演算によって求められた1心拍時間、最大流速、最小流速、PI、RIなどの各値を表示する。
【0017】
DSP27での処理のアルゴリズムは図2のフローチャートのとおりであり、FFTなどの各処理は、超音波スキャンに対応する一定時間周期で行われ、その演算ルーチンはタイマー31の割り込みによって起動させられる。タイマー割り込み時間間隔はここでは2ミリ秒としている。このタイマー割り込みで、まずFFT処理32が行われる。すなわち、レンジゲート25内の信号のクラッタ除去回路26を経たものをリングバッファ等のメモリに積算しておき、そのデータをタイマー割り込みで取り出してFFT処理する。FFT処理32によって得られたスペクトルはリアルタイムトレース処理33を受ける。
【0018】
このドプラスペクトルのデータは割り込み時間間隔で順次得られ、これらが一定数リングバッファ等で構成されたメモリに蓄積されることにより、ドプラスペクトルデータの時間トレースがなされる。このメモリ内容をそのまま表示すれば図3のような速度波形が表示される。メモリをリングバッファで構成し、現時点から一定割り込み回数前までのスペクトルデータをつねにメモリに蓄え、これを表示するようにすれば、図3のような波形が左方向にスクロールして表示されることになる。図3において横軸は時間tであり、縦軸は周波数シフトすなわち速度vとなる。
【0019】
ドプラスペクトルは、各周波数シフトごとの強度分布スペクトルであるから、厳密にいえば1点の値をとるわけではなく、図3においても縦軸方向にある幅を有する。このドプラスペクトルにおいて、最も強度の高い周波数シフトは特定領域で最も支配的な速度に対応し、周波数シフトの平均値は特定領域での速度の平均値に対応するので、これらを取り出せば、図3に示すような細い曲線で示される速度波形が得られることになる。
【0020】
この一定数のトレースデータに対してスムージングフィルタをかけてノイズを除去するとともに、そのスムージング後のデータをAとし、これを別のリングバッファ等に蓄積する(図2のステップ34)。スムージングフィルタはFIRローパスフィルタで構成する。
【0021】
つぎにこのリングバッファ等に蓄積された一定数のデータAに対して微分フィルタ(FIRフィルタで構成される)をかけて微分値Dを取得する(ステップ35)。前のタイマー割り込みルーチンで演算された値Dに対して、今回の割り込みで得られた値Dの符号が変っているかどうかを調べ(ステップ36)、もし変っているなら、今回の割り込みで得られたAを極値として認識し、前回の割り込みまでに得られた値を入れておく最大流速レジスタRmaxの値の類似の範囲に今回のAが入っているかどうかを調べる(ステップ41、42、43)。ここで極値のみを処理していることにより、演算を簡潔にすることができる。最大流速や最小流速は極値をとるはずであるから、極値が得られた割り込みで最大流速・最小流速を判定すれば足りるからである。
【0022】
ここでは、この類似の範囲は前回までのRmaxの値の1.2倍〜0.8倍としており、今回のAがその範囲内であれば類似であるとしている。そのため、まずステップ41で今回のAが前回までのRmaxの値より大きいかを調べ、大きければ、ステップ42で前回までのRmaxの値の1.2倍と今回のAとを比較し、今回のAの方が大きくなければ、ステップ43で前回までのRmaxの値の0.8倍と今回のAとを比較して今回のAの方が大きければ類似性ありとする。類似性ありと判定されたとき、ステップ44で心拍時間レジスタRhbtに、ハートレートカウンタHRCの値に割り込み時間間隔をかけた値を入れる。ここで、ハートレートカウンタHRCは各割り込みで順次インクリメントしていくカウンタである。
【0023】
さらに、このRhbtの値が200ミリ秒より小さくないことをステップ45で判定したとき、1心拍が認識できたものとして、このRhbtの値を1心拍時間とする。このように、前回の割り込みまでの間に得られたRmaxの値に類似するAが今回の割り込みで得られ、かつRmaxの値が得られた割り込みから今回の割り込みまでの時間として求めた1心拍時間が200ミリ秒以上の正常値であるとの判定がなされたとき、1心拍が認識されたとしている。ここでは1心拍が認識できたことを「ロックした」と表現し、ロックした割り込みで、最大流速レジスタRmaxの値をその心拍内での最大流速、最小流速レジスタRminの値をその心拍内での最小流速とする。そして、平均流速レジスタRaveは後述のように各割り込みで得たA値を順次積算しているので、そのRaveの値を心拍時間レジスタRhbtの値で割った値をその心拍内での平均流速とする。これら最大流速、最小流速、平均流速および心拍時間の各値をビデオ表示装置28が持つプロセッサに送り(ステップ46)、PI、RI等を演算してディスプレイ装置29の画面に表示する。このロックした割り込みでは、その後、Rmaxの値をAに入れ替えるとともにRminの値もAに入れ替え、Rhbtは今回得られた値のままとし、さらにRaveの値は0として、各レジスタを初期化する(ステップ47)。この後割り込みを終了する(ステップ48)。
【0024】
ステップ45で、Rhbtの値が200ミリ秒より小さいと判定されたときは、システムの仕様として定めた最短心拍時間の200ミリ秒より短い心拍時間が検出されたこととなるので、ステップ49に進み、心拍が正常に検知できなかったと判断し(これをアンロック状態と呼ぶ)、Rmax=0、Rmin=0、Rave=0、Rhbt=6秒(仕様内の検出しうる最大心拍時間)、HRC=0と初期化し、その後割り込みを終了する(ステップ48)。
【0025】
もし、極値が得られた割り込みでAが前回までのRmaxの値の類似の範囲に入らなかった場合、大きい方に外れたときはステップ42からステップ51に進み、小さい方に外れたときはステップ43からステップ61へと進む。大きい方に外れた場合は、連続した安定な検出が検出系で行われていないと判断し(アンロック)、ステップ51に進んで、Rmax=0、Rmin=0、Rave=0、Rhbt=6秒、HRC=0と初期化し、その後割り込みを終了する(ステップ52)。
【0026】
小さい方に外れた場合は、ステップ61でその値Aが最小流速であるかどうかを判定する。そのため、このステップ61では、(Rmax−A)の絶対値が(Rmax−Rmin)の絶対値より大きいかどうかを判定する。ここでRmaxおよびRminは前回までの最大流速レジスタおよび最小流速レジスタの値である。これが大きいと判定されたときは、今回のAが前回までのRminより小さいということであるから、最小流速であると判定されたこととなり、そのため、ステップ62でRminの値をこのAに更新する。その後、HRCをカウントアップし(ステップ63)、RaveにAを積算して(ステップ64)、割り込みを終了する(ステップ48)。最小流速ではないと判定されたときは、Rminの更新は行わずに、HRCをカウントアップし(ステップ63)、RaveにAを積算して(ステップ64)、割り込みを終了する(ステップ48)。
【0027】
ステップ36で今回のDが前のDと符号が違うことを判定したが、符号が違わない場合は、ステップ71に進んでHRCをカウントアップし、つぎにステップ72でRaveにAを積算する。さらに、ステップ73で、今回の割り込みまでの積算時間つまりHRC*2ミリ秒が前回までのRhbtの値プラス2秒よりも大きいかどうかを判定する。大であれば、安定な心拍の検出がなされていないと判断し、ステップ74でアンロック状態とし、Rmax=0、Rmin=0、Rave=0、Rhbt=6秒、HRC=0と初期化し、その後割り込みを終了する(ステップ75)。小である場合は、このようなアンロックのステップ74を経ずに終了する(ステップ75)。
【0028】
なお、この例では、上記からもわかるとおり、検出可能な心拍時間の仕様を6秒〜200ミリ秒とするとともに、1心拍検知の指標となる最大流速の前回の値に対する類似範囲を1.2倍〜0.8倍としているが、これらは適宜変更可能である。また、最小流速の条件は前回までのRmaxの値と今回のA値との差の絶対値が最大になることとしているが、最小流速に関しては、本来最小流速でなく拡張終末期血流速度を求めるべきであるところ、拡張終末期血流速度は最小流速で近似できるためこれを用いているので、最大流速の一つ前の極値とすることも可能である。通常、最大流速の一つ前の極値は、ほぼ最小流速と同じであることが多いからである。そのため、これについては、システムの定義上の問題ともいえるが、このように最大流速の一つ前の極値をもって最小流速とすれば、この値をRminに反映することは非常に簡単になることは指摘できる。さらに、タイマー割り込みの時間間隔は2ミリ秒とし、心拍時間が前回の心拍時間プラス2秒を超えるとアンロック状態としているが、これらの時間設定についても変更可能である。また、数心拍の平均をとって各値を求めることにより、より確度の高い値にすることも自明であろう。その他、微分フィルタとスムージングフィルタとは一つのFIRフィルタで行うことも可能であるなど、アルゴリズムおよび具体的な構成に関してはいろいろのバリエーションがあり得ることはもちろんである。
【0029】
【発明の効果】
以上説明したように、この発明の超音波診断装置によれば、超音波スキャンを行っている最中にリアルタイムでPI、RIなどを自動的に算出してこれらを流速波形と同様にディスプレイ画面に表示することが可能となり、ユーザは画質調整以外にとくに何らの操作も要することなく、これらの値を読み取ることができるため、疾患の診断がきわめて容易になる。
【図面の簡単な説明】
【図1】この発明の実施の形態を示すブロック図。
【図2】同実施形態のDSPの動作を説明するためのフローチャート。
【図3】流速波形を示すタイムチャート。
【符号の説明】
11 超音波プローブ
12 パルス制御回路
13 送信ドライバ回路
21 受信プリアンプ
22 A/D変換器
23 ビームフォーマ
24 検波回路
25 レンジゲート
26 クラッタ除去回路
27 DSP
28 ビデオ表示装置
29 ディスプレイ装置
Claims (1)
- 超音波プローブと、該超音波プローブを介して超音波の送信および受信を行う超音波送受信手段と、受信した超音波反射波信号の周波数シフトを求める検波手段と、上記周波数シフトのスペクトルを得るフーリエ変換手段と、該スペクトルの時間トレースデータを時間方向に微分する手段と、微分値の符号が変化したときのデータを極値として求め、この極値により最大流速および最小流速を判定する手段とを備えることを特徴とする超音波診断装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003010531A JP2004222754A (ja) | 2003-01-20 | 2003-01-20 | 超音波診断装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003010531A JP2004222754A (ja) | 2003-01-20 | 2003-01-20 | 超音波診断装置 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008157440A Division JP2008212746A (ja) | 2008-06-17 | 2008-06-17 | 超音波診断装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2004222754A true JP2004222754A (ja) | 2004-08-12 |
Family
ID=32899696
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003010531A Pending JP2004222754A (ja) | 2003-01-20 | 2003-01-20 | 超音波診断装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2004222754A (ja) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006102489A (ja) * | 2004-09-07 | 2006-04-20 | Toshiba Corp | 超音波ドプラ診断装置及び診断パラメータ計測方法 |
JP2007007200A (ja) * | 2005-06-30 | 2007-01-18 | Toshiba Corp | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム |
JP2007205736A (ja) * | 2006-01-31 | 2007-08-16 | Hioki Ee Corp | 表示制御装置 |
US8303507B2 (en) | 2004-09-07 | 2012-11-06 | Kabushiki Kaisha Toshiba | Ultrasonic doppler diagnostic apparatus and measuring method of diagnostic parameter |
CN105286920A (zh) * | 2015-11-05 | 2016-02-03 | 无锡祥生医学影像有限责任公司 | 超声多普勒信息自适应干扰抑制方法 |
CN111801052A (zh) * | 2018-02-07 | 2020-10-20 | 挪威科技大学 | 超声血流监测 |
-
2003
- 2003-01-20 JP JP2003010531A patent/JP2004222754A/ja active Pending
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006102489A (ja) * | 2004-09-07 | 2006-04-20 | Toshiba Corp | 超音波ドプラ診断装置及び診断パラメータ計測方法 |
JP4744994B2 (ja) * | 2004-09-07 | 2011-08-10 | 株式会社東芝 | 超音波ドプラ診断装置及び診断パラメータ計測方法 |
US8303507B2 (en) | 2004-09-07 | 2012-11-06 | Kabushiki Kaisha Toshiba | Ultrasonic doppler diagnostic apparatus and measuring method of diagnostic parameter |
JP2007007200A (ja) * | 2005-06-30 | 2007-01-18 | Toshiba Corp | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム |
JP2007205736A (ja) * | 2006-01-31 | 2007-08-16 | Hioki Ee Corp | 表示制御装置 |
CN105286920A (zh) * | 2015-11-05 | 2016-02-03 | 无锡祥生医学影像有限责任公司 | 超声多普勒信息自适应干扰抑制方法 |
CN111801052A (zh) * | 2018-02-07 | 2020-10-20 | 挪威科技大学 | 超声血流监测 |
CN111801052B (zh) * | 2018-02-07 | 2024-06-07 | 西蒙医疗公司 | 超声血流监测 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2020044769A1 (ja) | 超音波診断装置および超音波診断装置の制御方法 | |
US6770034B2 (en) | Ultrasonic diagnostic apparatus | |
JP4829960B2 (ja) | 超音波診断装置 | |
KR102025328B1 (ko) | 평면파 합성을 이용한 초음파 벡터 도플러 영상의 생성 장치 및 방법 | |
JP4189405B2 (ja) | 超音波診断装置 | |
JPH1133024A (ja) | 超音波ドプラ診断装置 | |
JP2001178720A (ja) | 連続データ獲得を用いた超音波フロー・イメージングにおける運動の可視化のための方法および装置 | |
JP3410821B2 (ja) | 超音波診断装置 | |
US20050119573A1 (en) | Method and system for quantification of arterial stenosis | |
JP2003126090A (ja) | 生体信号計測装置及び超音波診断装置 | |
JP2009039277A (ja) | 超音波診断装置 | |
JP2004222754A (ja) | 超音波診断装置 | |
JP2008212746A (ja) | 超音波診断装置 | |
JP2001309918A (ja) | 超音波診断装置 | |
JP3578680B2 (ja) | 超音波診断装置 | |
JP3820162B2 (ja) | 循環動態測定装置 | |
JP3668687B2 (ja) | 脈波伝播速度計測装置及び超音波診断装置 | |
JP2004329609A (ja) | 超音波診断装置 | |
JP2007090003A (ja) | 超音波診断装置及びその制御方法 | |
US9737283B2 (en) | Ultrasound measurement apparatus and ultrasound measurement method | |
JP3238467B2 (ja) | 超音波ドップラ診断装置 | |
JP4357260B2 (ja) | 加速度脈波測定装置 | |
JP3151463B2 (ja) | 超音波診断装置 | |
KR100842234B1 (ko) | 도플러 스펙트럼의 베이스라인 및 스케일을 조절하는 영상처리 시스템 및 방법 | |
JP4679141B2 (ja) | 超音波診断装置および超音波診断画像の表示方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Effective date: 20050308 Free format text: JAPANESE INTERMEDIATE CODE: A621 |
|
RD02 | Notification of acceptance of power of attorney |
Effective date: 20050401 Free format text: JAPANESE INTERMEDIATE CODE: A7422 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20050401 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20050404 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070814 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080129 |
|
A521 | Written amendment |
Effective date: 20080326 Free format text: JAPANESE INTERMEDIATE CODE: A523 |
|
A131 | Notification of reasons for refusal |
Effective date: 20080422 Free format text: JAPANESE INTERMEDIATE CODE: A131 |
|
A02 | Decision of refusal |
Effective date: 20080826 Free format text: JAPANESE INTERMEDIATE CODE: A02 |