JP4832211B2 - Ultrasonic diagnostic apparatus and image display apparatus - Google Patents

Ultrasonic diagnostic apparatus and image display apparatus Download PDF

Info

Publication number
JP4832211B2
JP4832211B2 JP2006221165A JP2006221165A JP4832211B2 JP 4832211 B2 JP4832211 B2 JP 4832211B2 JP 2006221165 A JP2006221165 A JP 2006221165A JP 2006221165 A JP2006221165 A JP 2006221165A JP 4832211 B2 JP4832211 B2 JP 4832211B2
Authority
JP
Japan
Prior art keywords
data
image
time
series data
mode
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.)
Active
Application number
JP2006221165A
Other languages
Japanese (ja)
Other versions
JP2008043495A (en
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.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
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 GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Priority to JP2006221165A priority Critical patent/JP4832211B2/en
Publication of JP2008043495A publication Critical patent/JP2008043495A/en
Application granted granted Critical
Publication of JP4832211B2 publication Critical patent/JP4832211B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Description

本発明は、超音波診断装置及び画像表示装置に関する。   The present invention relates to an ultrasonic diagnostic apparatus and an image display apparatus.

いわゆるDモード画像を表示する超音波診断装置が知られている(例えば特許文献1)。このような超音波診断装置では、プローブからの受信信号に対して位相検波等の処理を施してドップラ情報を含む時系列データを取得し、その時系列データに対して一定のデータ点数(ウインドウサイズ:Window sizeといわれることがある。)を含む時間範囲毎にフーリエ変換を行いドプラスペクトルを得ることを、時間範囲をずらしつつ複数回行い、複数の時点毎(厳密には複数の時間範囲毎)のドプラスペクトルを得る。そして、複数のドプラスペクトルに基づいてDモード画像の画像データを生成する。
特開平8−229035号公報
An ultrasonic diagnostic apparatus that displays a so-called D-mode image is known (for example, Patent Document 1). In such an ultrasonic diagnostic apparatus, processing such as phase detection is performed on a received signal from a probe to acquire time series data including Doppler information, and a certain number of data points (window size: In order to obtain a Doppler spectrum by performing Fourier transform for each time range including “window size”), the Doppler spectrum is obtained several times while shifting the time range. Obtain a Doppler spectrum. Then, image data of the D mode image is generated based on the plurality of Doppler spectra.
JP-A-8-229035

図7は、従来のDモード画像を説明する図である。図7のDモード画像において、横軸は時間であり、縦軸は血流等の速度であり、各画素の輝度は、横軸の時間及び縦軸の速度によって特定されるドプラスペクトルのパワーが大きくなるほど高くなっている。なお、図7では、輝度が低い領域をハッチングして示している。   FIG. 7 is a diagram for explaining a conventional D-mode image. In the D-mode image of FIG. 7, the horizontal axis is time, the vertical axis is the speed of blood flow, and the luminance of each pixel is determined by the time of the horizontal axis and the power of the Doppler spectrum specified by the speed of the vertical axis. The higher it is, the higher it is. In FIG. 7, a region with low luminance is hatched.

各時点(厳密には各時間範囲)におけるドプラスペクトルを算出する際に、時系列データのデータ点数を多くすると、図7(a)に示すように、周囲に比較して輝度が低い(パワーが低い)複数の斑点501が生じる。この場合、ユーザは超音波診断装置の感度が低いと感じてしまう。   When calculating the Doppler spectrum at each time point (strictly, each time range), if the number of data points of the time-series data is increased, the luminance is lower than the surroundings (the power is lower) as shown in FIG. Multiple spots 501 are generated (low). In this case, the user feels that the sensitivity of the ultrasonic diagnostic apparatus is low.

一方、時系列データのデータ点数を少なくすると、図7(b)に示すように、周囲に比較して輝度が高い又は低い線像、すなわち、縦筋502が生じ、輝度が高い領域と低い領域との境界領域503にはギザギザが現われる。この場合、ユーザは輝度が高い領域の境界を把握しにくい。   On the other hand, when the number of data points of the time-series data is reduced, as shown in FIG. 7B, a line image having a higher or lower luminance than the surroundings, that is, a vertical streak 502 is generated, and a region having a high luminance and a region having a low luminance. Jagged edges appear in the boundary region 503. In this case, it is difficult for the user to grasp the boundary of the region with high luminance.

すなわち、時系列データのデータ点数はDモード画像の画質に影響する。そして、一般的なノイズ除去の技術思想に従えば、データ点数を多くすれば画質が向上しそうに考えられるが、実際には、時系列データのデータ点数を少なくした場合だけでなく、データ点数を多くした場合にも画質の低下が生じており、時系列データのデータ点数を調整するだけでは、高画質のDモード画像を表示することが容易ではなかった。   That is, the number of data points of time series data affects the image quality of the D-mode image. According to the general noise reduction technical idea, it seems that the image quality is likely to improve if the number of data points is increased. However, in actuality, not only when the number of data points of the time-series data is reduced, but the number of data points is also increased. In many cases, the image quality is deteriorated, and it is not easy to display a high-quality D-mode image only by adjusting the number of data of the time series data.

なお、同じ装置であっても、実験的に人工的に作成した、十分に強度があり、且つ、一定の周波数を有する連続したSin波等を受信した場合のDモード画像の場合、縦筋や斑点が目立たないことを考慮すると、縦筋や斑点が生じるのは、生体内からの反射音を拾っているために、不連続な周波数成分を含んだ受信信号を高速フーリエ変換していることや、常に十分な信号強度が得られているわけではない事が原因と推測される。   Note that even in the same device, in the case of a D-mode image that is produced artificially experimentally and has a sufficiently strong and continuous sine wave having a constant frequency, Considering that the spots are not conspicuous, vertical stripes and spots are generated because the reflected sound from the living body is picked up, so that the received signal containing discontinuous frequency components is fast Fourier transformed. The reason is that sufficient signal strength is not always obtained.

本発明の目的は、Dモード画像の画質を向上できる超音波診断装置及び画像表示装置を提供することにある。   An object of the present invention is to provide an ultrasonic diagnostic apparatus and an image display apparatus that can improve the image quality of a D-mode image.

本発明の超音波診断装置は、超音波を被検体に照射し、反射された超音波を電気信号に変換して出力するプローブと、前記プローブからの信号に基づいてドプラ情報を含む時系列データを生成する時系列データ生成手段と、前記時系列データのうち所定の時間範囲の時系列データにフーリエ変換を行いドプラスペクトルデータを生成することを、前記所定の時間範囲を一定時間ずつずらして複数回行い、複数時の前記ドプラスペクトルデータを生成するスペクトルデータ生成手段と、前記複数時のドプラスペクトルデータに基づいてDモード画像の画像データを生成する画像生成手段と、前記画像データに基づいてDモード画像を表示する表示手段と、を備え、前記スペクトルデータ生成手段は、前記所定の時間範囲の時系列データのうち所定のデータ点数の時系列データにフーリエ変換を行いドプラスペクトルデータを生成することを、前記データ点数を異ならせて複数回行い、複数種類の前記ドプラスペクトルデータを生成し、前記画像生成手段は、前記複数種類のドプラスペクトルデータに基づく複数種類のDモード画像が合成された合成Dモード画像の画像データを生成する。   The ultrasonic diagnostic apparatus of the present invention irradiates a subject with ultrasonic waves, converts the reflected ultrasonic waves into electrical signals, and outputs the time series data including Doppler information based on the signals from the probes. A plurality of time-series data generating means for generating Doppler spectrum data by performing Fourier transform on the time-series data in a predetermined time range of the time-series data to generate Doppler spectrum data. A spectrum data generating means for generating the Doppler spectrum data for a plurality of times, an image generating means for generating image data of a D-mode image based on the Doppler spectrum data for the plurality of times, and a D based on the image data. Display means for displaying a mode image, wherein the spectrum data generating means is a part of the time-series data in the predetermined time range. Generating the Doppler spectrum data by performing Fourier transform on the time series data of the number of data points, a plurality of times with different data points, generating a plurality of types of the Doppler spectrum data, and the image generating means, Image data of a combined D-mode image obtained by combining a plurality of types of D-mode images based on a plurality of types of Doppler spectrum data is generated.

好適には、前記合成Dモード画像は、前記複数種類のDモード画像の互いに対応する画素間において輝度が平均された画像である。   Preferably, the composite D-mode image is an image in which luminance is averaged between pixels corresponding to each other of the plurality of types of D-mode images.

好適には、前記合成Dモード画像は、前記複数種類のDモード画像の互いに対応する画素間において最大の輝度が選択された画像である。   Preferably, the composite D-mode image is an image in which the maximum luminance is selected between pixels corresponding to each other of the plurality of types of D-mode images.

好適には、前記合成Dモード画像は、前記複数種類のDモード画像のいずれか一のDモード画像の所定の領域において他の画像が合成された画像データである。   Preferably, the composite D-mode image is image data obtained by combining another image in a predetermined region of any one of the plurality of types of D-mode images.

好適には、前記所定の領域は、前記一のDモード画像において所定の大きさの斑点が生じる領域であり、前記他の画像は、前記一のDモード画像よりも少ないデータ点数の前記時系列データに基づくDモード画像である。   Preferably, the predetermined area is an area where a spot having a predetermined size occurs in the one D-mode image, and the other image is the time series having a smaller number of data points than the one D-mode image. It is a D-mode image based on data.

好適には、前記所定の領域は、前記一のDモード画像においてドプラ偏移周波数又はドプラ偏移周波数に相関する物理量を示す軸に沿う線が生じる領域であり、前記他の画像は、前記一のDモード画像よりも多いデータ点数の前記時系列データに基づくDモード画像である。   Preferably, the predetermined region is a region where a line along an axis indicating a physical quantity correlated with a Doppler shift frequency or a Doppler shift frequency is generated in the one D-mode image, and the other image is the one mode. The D-mode image is based on the time-series data having a larger number of data points than the D-mode image.

好適には、前記所定の時間範囲において前記データ点数が互いに異なる時系列データは、一のデータ点数の時系列データのうち一部の時間範囲の時系列データ全てが他のデータ点数の時系列データとなるように、前記所定の時間範囲の時系列データから選択されている。   Preferably, the time-series data having different data points in the predetermined time range is the time-series data in which all of the time-series data in a part of the time range among the time-series data of one data point are other data points. The time series data in the predetermined time range is selected.

好適には、前記スペクトルデータ生成手段は、前記データ点数が互いに異なる時系列データにそれぞれ、時系列データの時間中心側において値が大きくなる窓関数を掛けてからフーリエ変換を行い、前記複数種類のドプラスペクトルデータを生成し、前記所定の時間範囲において前記データ点数が互いに異なる時系列データは、各時系列データの時間中心が互いに同一になるように、前記所定の時間範囲の時系列データから選択されている。   Preferably, the spectrum data generating means performs a Fourier transform after multiplying the time series data having different data points by a window function having a larger value on the time center side of the time series data, and performing the Fourier transform, Time-series data for which Doppler spectrum data is generated and whose data points are different from each other in the predetermined time range is selected from the time-series data in the predetermined time range so that the time centers of the respective time-series data are the same. Has been.

本発明の画像表示装置は、ドプラ情報を含む時系列データを生成する時系列データ生成手段と、前記時系列データのうち所定の時間範囲の時系列データにフーリエ変換を行いドプラスペクトルデータを生成することを、前記所定の時間範囲を一定時間ずつずらして複数回行い、複数時の前記ドプラスペクトルデータを生成するスペクトルデータ生成手段と、前記複数時のドプラスペクトルデータに基づいてDモード画像の画像データを生成する画像生成手段と、前記画像データに基づいてDモード画像を表示する表示手段と、を備え、前記スペクトルデータ生成手段は、前記所定の時間範囲の時系列データのうち所定のデータ点数の時系列データにフーリエ変換を行いドプラスペクトルデータを生成することを、前記データ点数を異ならせて複数回行い、複数種類の前記ドプラスペクトルデータを生成し、前記画像生成手段は、前記複数種類のドプラスペクトルデータに基づく複数種類のDモード画像が合成された合成Dモード画像の画像データを生成する。   An image display device of the present invention generates time-series data including time-series data including Doppler information, and generates Doppler spectrum data by performing Fourier transform on the time-series data in a predetermined time range of the time-series data. A plurality of times by shifting the predetermined time range by a predetermined time and generating the Doppler spectrum data at a plurality of times; and image data of the D-mode image based on the Doppler spectrum data at the plurality of times. And a display means for displaying a D-mode image based on the image data, wherein the spectrum data generating means has a predetermined number of data points in the time series data in the predetermined time range. The time series data is subjected to Fourier transform to generate Doppler spectrum data. Several times are performed to generate a plurality of types of Doppler spectrum data, and the image generation unit generates image data of a combined D-mode image obtained by combining a plurality of types of D-mode images based on the plurality of types of Doppler spectrum data. .

好適には、前記合成Dモード画像は、前記複数種類のDモード画像の互いに対応する画素間において輝度が平均された画像である。   Preferably, the composite D-mode image is an image in which luminance is averaged between pixels corresponding to each other of the plurality of types of D-mode images.

好適には、前記合成Dモード画像は、前記複数種類のDモード画像の互いに対応する画素間において最大の輝度が選択された画像である。   Preferably, the composite D-mode image is an image in which the maximum luminance is selected between pixels corresponding to each other of the plurality of types of D-mode images.

好適には、前記合成Dモード画像は、前記複数種類のDモード画像のいずれか一のDモード画像の所定の領域において他の画像が合成された画像データである。   Preferably, the composite D-mode image is image data obtained by combining another image in a predetermined region of any one of the plurality of types of D-mode images.

好適には、前記所定の領域は、前記一のDモード画像において所定の大きさの斑点が生じる領域であり、前記他の画像は、前記一のDモード画像よりも少ないデータ点数の前記時系列データに基づくDモード画像である。   Preferably, the predetermined area is an area where a spot having a predetermined size occurs in the one D-mode image, and the other image is the time series having a smaller number of data points than the one D-mode image. It is a D-mode image based on data.

好適には、前記所定の領域は、前記一のDモード画像においてドプラ偏移周波数又はドプラ偏移周波数に相関する物理量を示す軸に沿う線が生じる領域であり、前記他の画像は、前記一のDモード画像よりも多いデータ点数の前記時系列データに基づくDモード画像である。   Preferably, the predetermined region is a region where a line along an axis indicating a physical quantity correlated with a Doppler shift frequency or a Doppler shift frequency is generated in the one D-mode image, and the other image is the one mode. The D-mode image is based on the time-series data having a larger number of data points than the D-mode image.

好適には、前記所定の時間範囲において前記データ点数が互いに異なる時系列データは、一のデータ点数の時系列データのうち一部の時間範囲の時系列データ全てが他のデータ点数の時系列データとなるように、前記所定の時間範囲の時系列データから選択されている。   Preferably, the time-series data having different data points in the predetermined time range is the time-series data in which all of the time-series data in a part of the time range among the time-series data of one data point are other data points. The time series data in the predetermined time range is selected.

好適には、前記スペクトルデータ生成手段は、前記データ点数が互いに異なる時系列データにそれぞれ、時系列データの時間中心側において値が大きくなる窓関数を掛けてからフーリエ変換を行い、前記複数種類のドプラスペクトルデータを生成し、前記所定の時間範囲において前記データ点数が互いに異なる時系列データは、各時系列データの時間中心が互いに同一になるように、前記所定の時間範囲の時系列データから選択されている。   Preferably, the spectrum data generating means performs a Fourier transform after multiplying the time series data having different data points by a window function having a larger value on the time center side of the time series data, and performing the Fourier transform, Time-series data for which Doppler spectrum data is generated and whose data points are different from each other in the predetermined time range is selected from the time-series data in the predetermined time range so that the time centers of the respective time-series data are the same. Has been.

本発明によれば、Dモード画像の画質を向上できる。   According to the present invention, the image quality of a D-mode image can be improved.

図1は、本発明の実施形態に係る超音波診断装置1の全体構成の要部を示すブロック図である。超音波診断装置1は、例えば、いわゆるパルスドプラ法でDモード画像を得る超音波診断装置として構成されている。すなわち、超音波診断装置1は、超音波の送信波をある長さのパルス波として送り、その反射波が返ってきた後に、次のパルス波を送るというように、間欠的に送受波を行う。   FIG. 1 is a block diagram showing the main part of the overall configuration of an ultrasonic diagnostic apparatus 1 according to an embodiment of the present invention. The ultrasonic diagnostic apparatus 1 is configured as, for example, an ultrasonic diagnostic apparatus that obtains a D-mode image by a so-called pulse Doppler method. That is, the ultrasonic diagnostic apparatus 1 transmits and receives waves intermittently, such as transmitting an ultrasonic transmission wave as a pulse wave of a certain length, and sending the next pulse wave after the reflected wave returns. .

超音波診断装置1は、被検体に超音波を照射し、反射された超音波を受信するプローブ11と、プローブ11との間で電気信号を送受信する送受信部12と、送受信部12からの信号に基づいて画像データの生成に必要なデータを生成するデータ処理部13と、データ処理部13により生成されたデータに基づいて画像データを生成する画像処理部14と、画像処理部14からの信号に基づいて超音波画像を表示する表示部15と、ユーザの入力操作を受け付ける操作部16と、操作部16からの信号等に基づいて各部12〜15の制御を実行する制御部17とを備えている。   The ultrasonic diagnostic apparatus 1 irradiates a subject with ultrasonic waves, receives a reflected ultrasonic wave, a transmission / reception unit 12 that transmits / receives an electrical signal between the probes 11, and a signal from the transmission / reception unit 12. A data processing unit 13 that generates data necessary for generating image data based on the data, an image processing unit 14 that generates image data based on the data generated by the data processing unit 13, and a signal from the image processing unit 14 A display unit 15 that displays an ultrasonic image based on the control unit 16, an operation unit 16 that receives a user input operation, and a control unit 17 that executes control of the units 12 to 15 based on a signal from the operation unit 16. ing.

プローブ11は、超音波と電気信号との相互変換を行なうトランスデューサを含んで構成され、入力された電気信号を超音波に変換して被検体に送信するとともに、受信した超音波(エコー信号)を電気信号に変換して出力する。   The probe 11 includes a transducer that performs mutual conversion between an ultrasonic wave and an electric signal, converts the input electric signal into an ultrasonic wave, transmits the ultrasonic wave to the subject, and receives the received ultrasonic wave (echo signal). Convert to electrical signal and output.

送受信部12は、制御部17からの制御信号に基づいて所定の繰り返し周波数で電気信号をプローブ11に出力する。また、プローブ11からの電気信号をデータ処理部13に出力する。   The transmission / reception unit 12 outputs an electrical signal to the probe 11 at a predetermined repetition frequency based on the control signal from the control unit 17. In addition, an electrical signal from the probe 11 is output to the data processing unit 13.

データ処理部13は、送受信部12からの信号に基づいてドプラ情報を含む時系列データを生成する時系列データ生成部19と、時系列データ生成部19の生成した時系列データに対してフーリエ変換を行ってドプラスペクトルデータを生成するスペクトルデータ生成部20とを備えている。データ処理部13は、例えばコンピュータにより構成されている。   The data processing unit 13 includes a time series data generation unit 19 that generates time series data including Doppler information based on a signal from the transmission / reception unit 12, and a Fourier transform on the time series data generated by the time series data generation unit 19. And a spectrum data generation unit 20 for generating Doppler spectrum data. The data processing unit 13 is configured by a computer, for example.

時系列データ生成部19は、特に図示しないが、例えば、入力された信号から所定の深度(図2のサンプルボリューム105参照)の信号を抽出する信号抽出部、入力された信号を直交検波してドプラ信号を出力する位相検波部、入力されたアナログ信号をA/D変換するAD変換部を備えている。時系列データ生成部19は、送受信部12からの信号に上述の各部による処理を施すことにより、離散的なデジタルデータである時系列データを生成する。時系列データは、ドプラ偏移周波数で振動するドプラ信号の強度の情報を一定の時間間隔毎(例えば、パルスドプラ法では繰り返し周波数に相当する周期毎)に保持している。なお、時系列データ生成部19の上述の各部の処理は、適宜な順番で行ってよい。   Although not particularly shown, the time series data generation unit 19 is, for example, a signal extraction unit that extracts a signal of a predetermined depth (see the sample volume 105 in FIG. 2) from the input signal, and orthogonally detects the input signal. A phase detector for outputting a Doppler signal and an AD converter for A / D conversion of the input analog signal are provided. The time-series data generation unit 19 generates time-series data that is discrete digital data by performing processing by the above-described units on the signal from the transmission / reception unit 12. The time-series data holds information on the intensity of the Doppler signal that vibrates at the Doppler shift frequency at regular time intervals (for example, every period corresponding to the repetition frequency in the pulse Doppler method). In addition, you may perform the process of the above-mentioned each part of the time series data generation part 19 in an appropriate order.

スペクトルデータ生成部20は、時系列データ生成部19により生成された時系列データに対してフーリエ変換を行い、ドプラ偏移周波数毎のパワーの値を保持するドプラスペクトルデータを生成する。フーリエ変換は、例えば高速フーリエ変換(FFT)により行われる。   The spectrum data generation unit 20 performs Fourier transform on the time series data generated by the time series data generation unit 19 to generate Doppler spectrum data that holds a power value for each Doppler shift frequency. The Fourier transform is performed by, for example, fast Fourier transform (FFT).

なお、データ処理部13は、時系列データ生成部19やスペクトルデータ生成部20の他にも、データ(信号)からノイズを除去するフィルタ部や、他のモード(Mモード、Bモード、カラーモード)の画像を生成するためのデータを生成する各種のデータ生成部等を備えている。   In addition to the time series data generation unit 19 and the spectrum data generation unit 20, the data processing unit 13 includes a filter unit that removes noise from data (signal) and other modes (M mode, B mode, color mode). ) Includes various data generation units that generate data for generating an image.

画像処理部14は、データ処理部13により生成されたデータに基づいて画像データを生成する。例えば、画像処理部14は、スペクトルデータ生成部20により生成されたドプラスペクトルデータに基づいてDモード画像を生成する。画像処理部14は、Dモード画像の画像データの他にも、他のモードの画像の画像データを生成可能である。画像処理部14は、生成した画像データに基づく映像信号を表示部15に出力する。画像処理部14は、例えばコンピュータにより構成されている。   The image processing unit 14 generates image data based on the data generated by the data processing unit 13. For example, the image processing unit 14 generates a D-mode image based on the Doppler spectrum data generated by the spectrum data generation unit 20. The image processing unit 14 can generate image data of an image in another mode in addition to the image data of the D mode image. The image processing unit 14 outputs a video signal based on the generated image data to the display unit 15. The image processing unit 14 is configured by a computer, for example.

表示部15は、例えばCRTディスプレイや液晶ディスプレイにより構成されており、画像処理部14から出力された映像信号に従って画像を表示する。操作部16は、例えばキーボードやポインティングディバイスにより構成されており、入力操作に応じた信号を制御部17に出力する。制御部17は、例えばコンピュータにより構成され、操作部16からの操作信号に基づいて、送受信部12、データ処理部13、画像処理部14、表示部15に制御信号を出力する。   The display unit 15 is configured by, for example, a CRT display or a liquid crystal display, and displays an image according to the video signal output from the image processing unit 14. The operation unit 16 is configured by, for example, a keyboard or a pointing device, and outputs a signal corresponding to the input operation to the control unit 17. The control unit 17 is configured by a computer, for example, and outputs control signals to the transmission / reception unit 12, the data processing unit 13, the image processing unit 14, and the display unit 15 based on an operation signal from the operation unit 16.

図2は、表示部15に表示される画像の一例を示している。画像101は、プローブ11により超音波の送受信が行われている間にリアルタイムに、又は、プローブ11による超音波の送受信が終了した後の任意の時期に、表示部15の画面に表示される。   FIG. 2 shows an example of an image displayed on the display unit 15. The image 101 is displayed on the screen of the display unit 15 in real time during transmission / reception of ultrasonic waves by the probe 11 or at an arbitrary time after the transmission / reception of ultrasonic waves by the probe 11 is completed.

画像101は、Bモード画像102と、Dモード画像103とを含んでいる。Bモード画像102は、輝度でエコー信号の強度を表した断層像の画像である。Bモード画像102では、Dモード画像103の表示対象部位を示すサンプルボリューム105が表示されている。サンプルボリューム105の位置及び大きさは、ユーザの操作部16への入力操作等により適宜に設定される。   The image 101 includes a B mode image 102 and a D mode image 103. The B-mode image 102 is a tomographic image that represents the intensity of the echo signal in terms of luminance. In the B mode image 102, a sample volume 105 indicating the display target portion of the D mode image 103 is displayed. The position and size of the sample volume 105 are appropriately set by an input operation to the operation unit 16 by the user.

Dモード画像103は、第1軸(例えば横軸)を時間とし、第1軸に直交する第2軸(例えば縦軸)をドプラ偏移周波数又はドプラ偏移周波数に相関する物理量(例えば血液の流速。以下、両者を総称して「ドプラ偏移周波数等」ということがある。)とし、輝度でドプラスペクトルのパワー又はパワーに相関する物理量(以下、両者を総称して「パワー等」ということがある。)を表した画像である。例えば、輝度は、ドプラスペクトルのパワー等が大きくなるほど高くなる。図2では、輝度が低い領域をハッチングして示している。   In the D-mode image 103, the first axis (for example, the horizontal axis) is time, and the second axis (for example, the vertical axis) orthogonal to the first axis is a Doppler shift frequency or a physical quantity (for example, blood) correlated with the Doppler shift frequency. Hereinafter, both are collectively referred to as “Doppler shift frequency, etc.”), and the physical quantity correlating with the power or power of the Doppler spectrum in luminance (hereinafter, both are collectively referred to as “power, etc.”). It is an image showing. For example, the luminance increases as the power of the Doppler spectrum increases. In FIG. 2, a region with low luminance is hatched.

Dモード画像103の大きさ、縦軸や横軸のレンジ、輝度とパワー等との対応付けは、操作部16への入力操作等により適宜に設定可能である。また、Dモード画像103とBモード画像102とを横方向に配列するなど、Dモード画像103の配置も適宜に設定可能である。   Correspondence between the size of the D-mode image 103, the ranges of the vertical and horizontal axes, the luminance and the power can be appropriately set by an input operation to the operation unit 16 or the like. Also, the arrangement of the D mode image 103 can be set as appropriate, for example, the D mode image 103 and the B mode image 102 are arranged in the horizontal direction.

設定の一例を述べる。600(縦)×800(横)ピクセルのCRTディスプレイの場合、Dモード画像103の大きさは、300(縦)×600(横)〜600(縦)×300(横)ピクセルに設定されてよい。Dモード画像103は、200(縦)×600(横)ピクセル、300(縦)×600(横)ピクセル、400(縦)×600(横)ピクセルなどの横長表示に設定されてもよいし、600(縦)×300(横)ピクセルなどの縦長表示に設定されてもよい。なお、ディスプレイの大きさ(inch)が異なる場合や液晶ディスプレイの場合も、上記の600(縦)×800(横)ピクセルのCRTディスプレイとの基本のドット数の相違を加味して若しくは加味せずに適宜に設定されてよい。Dモード画像103の横軸のレンジは1.0〜16.0(s)に、縦軸のレンジは10〜100(cm/s)に設定されてよい。輝度は、音圧を電気信号に変換して得られた電圧(パワー)に対して0〜63階調(6bit)が割り当てられてよい。なお、超音波診断装置1では、音圧を電気信号に変換して得られた数ボルトの電圧をプリアンプ増幅、デジタル変換、コンプレッション、フィルタリング等を行い、最終的に、0〜63の輝度値を得ている。すなわち、所定範囲の電圧に対して所定範囲の輝度値を割り当てることにより、間接的に所定範囲のパワーに対して所定範囲の輝度値を割り当てている。   An example of setting will be described. In the case of a 600 (vertical) × 800 (horizontal) pixel CRT display, the size of the D-mode image 103 may be set to 300 (vertical) × 600 (horizontal) to 600 (vertical) × 300 (horizontal) pixels. . The D-mode image 103 may be set to landscape display such as 200 (vertical) × 600 (horizontal) pixels, 300 (vertical) × 600 (horizontal) pixels, 400 (vertical) × 600 (horizontal) pixels, It may be set to a vertically long display such as 600 (vertical) × 300 (horizontal) pixels. In addition, in the case where the display size (inch) is different or in the case of a liquid crystal display, the difference in the basic number of dots from the 600 (vertical) × 800 (horizontal) pixel CRT display is taken into account or not. May be set appropriately. The range of the horizontal axis of the D mode image 103 may be set to 1.0 to 16.0 (s), and the range of the vertical axis may be set to 10 to 100 (cm / s). As the luminance, 0 to 63 gradations (6 bits) may be assigned to a voltage (power) obtained by converting sound pressure into an electrical signal. The ultrasonic diagnostic apparatus 1 performs preamplifier amplification, digital conversion, compression, filtering, and the like on a voltage of several volts obtained by converting sound pressure into an electrical signal, and finally sets a luminance value of 0 to 63. It has gained. That is, by assigning a predetermined range of luminance values to a predetermined range of voltage, the predetermined range of luminance values is indirectly assigned to the predetermined range of power.

図3は、データ処理部13の動作を説明する図である。なお、図3は、データ処理部13の動作を概念的に説明するためのものであり、図3において例示するデータ点数は作図上の便宜的なものに過ぎない。   FIG. 3 is a diagram for explaining the operation of the data processing unit 13. Note that FIG. 3 is for conceptually explaining the operation of the data processing unit 13, and the number of data points illustrated in FIG. 3 is merely for convenience in drawing.

まず、データ処理部13における従来と同様の動作を説明する。図3に示すように、時系列データ生成部19により、複数(図3では15個を例示)の時点データDCからなる時系列データD0が生成される。複数の時点データDCのサンプリング周期(時間間隔)T1は、繰り返し周波数の逆数であり、一定である。繰り返し周波数は、サンプリングの位置によるが、例えば100〜10kHzである。   First, an operation similar to the conventional one in the data processing unit 13 will be described. As shown in FIG. 3, the time-series data generating unit 19 generates time-series data D0 including a plurality (15 illustrated in FIG. 3) of time data DC. The sampling period (time interval) T1 of the plurality of time point data DC is the reciprocal of the repetition frequency and is constant. The repetition frequency depends on the sampling position, but is, for example, 100 to 10 kHz.

スペクトルデータ生成部20は、時系列データD0のうちデータ点数(時点データDCの数)がN1個(図3では6個を例示)の時系列データDA1に対してフーリエ変換を行い、時系列データDA1に基づくドプラスペクトルデータを生成する。同様に、時系列データDA1とは時間的な位置が異なる、データ点数がN1個の時系列データDA2、DA3(以下、単に「時系列データDA」といい、DA1〜DA3を区別しないことがある。)に順次フーリエ変換を行い、複数のドプラスペクトルデータを生成する。   The spectrum data generation unit 20 performs Fourier transform on the time series data DA1 of which the number of data points (the number of time point data DC) in the time series data D0 is N1 (six examples are illustrated in FIG. 3). Doppler spectrum data based on DA1 is generated. Similarly, time-series data DA2, DA3 (hereinafter simply referred to as “time-series data DA”), which is different in time position from time-series data DA1 and has N1 data points, may not be distinguished. .) Are sequentially subjected to Fourier transform to generate a plurality of Doppler spectrum data.

すなわち、スペクトルデータ生成部20は、所定の時間範囲T2の時系列データにフーリエ変換を行いドプラスペクトルデータを生成することを、所定の時間範囲T2を一定時間ずつずらして複数回行い、複数時のドプラスペクトルデータを生成する。なお、所定の時間範囲T2は、サンプリング周期T1×データ点数である。   That is, the spectrum data generation unit 20 performs Fourier transform on the time-series data in the predetermined time range T2 to generate Doppler spectrum data by shifting the predetermined time range T2 by a predetermined time and a plurality of times. Generate Doppler spectral data. The predetermined time range T2 is sampling period T1 × number of data points.

以上の従来と同様の動作に加え、スペクトルデータ生成部20は、時系列データDAのデータ点数N1とは異なるデータ点数N2(図3では4個を例示)の時系列データDB1〜DB3(以下、単に「時系列データDB」といい、DB1〜DB3を区別しないことがある。)に順次フーリエ変換を行い、複数のドプラスペクトルデータを生成する。   In addition to the above-described conventional operation, the spectrum data generation unit 20 uses time-series data DB1 to DB3 (hereinafter, four data points N2 are different from the data point N1 of the time-series data DA). Simply referred to as “time-series data DB”, DB1 to DB3 may not be distinguished.) Fourier transform is sequentially performed to generate a plurality of Doppler spectrum data.

すなわち、スペクトルデータ生成部20は、所定の時間範囲T2の時系列データのうち所定のデータ点数の時系列データにフーリエ変換を行いドプラスペクトルデータを生成することを、データ点数を異ならせて複数回(実施形態ではデータ点数N1の場合とN2の場合との2回)行い、複数種類(実施形態では2種類)のドプラスペクトルデータを生成する。   That is, the spectrum data generation unit 20 performs Fourier transform on the time series data of a predetermined number of data points out of the time series data of the predetermined time range T2 to generate Doppler spectrum data a plurality of times with different data points. (In the embodiment, the number of data points is N1 and N2 is twice), and multiple types (two types in the embodiment) of Doppler spectrum data are generated.

時系列データDBは、時系列データDAのうち一部の時間範囲の時点データDC全てにより構成されており、時系列データDAのサンプリング周期と時系列データDBのサンプリング周期とは互いに同一(T1)である。また、時系列データDBは、時系列データDAの中央の一部の時間範囲に設定されており、時系列データDAの時間中心Tacと時系列データDBの時間中心Tbcとは互いに一致している。   The time series data DB is composed of all the time point data DC in a part of the time range in the time series data DA, and the sampling period of the time series data DA and the sampling period of the time series data DB are the same (T1). It is. Further, the time series data DB is set to a partial time range in the center of the time series data DA, and the time center Tac of the time series data DA and the time center Tbc of the time series data DB coincide with each other. .

なお、例えば、N1=60〜90であり、N2=30〜36である。   For example, N1 = 60 to 90 and N2 = 30 to 36.

図4は、画像処理部14の動作を説明する概念図である。上述のように、スペクトルデータ生成部20により2種類のドプラスペクトルデータが生成される。従って、画像処理部14は、2種類のドプラスペクトルデータに基づいて、データ点数N1の場合のDモード画像103A(図4(a))と、データ点数N2の場合のDモード画像103B(図4(b))との2種類のDモード画像の画像データを生成可能である。ここで、上述したように、データ点数が比較的多いDモード画像103Aでは、斑点501が生じる。また、データ点数が比較的少ないDモード画像103Bでは、縦筋502が生じる。   FIG. 4 is a conceptual diagram illustrating the operation of the image processing unit 14. As described above, two types of Doppler spectrum data are generated by the spectrum data generation unit 20. Therefore, the image processing unit 14 based on the two types of Doppler spectrum data, the D-mode image 103A (FIG. 4A) in the case of N1 data points and the D-mode image 103B in the case of N2 data points (FIG. 4). (B)) and two types of D-mode image data can be generated. Here, as described above, the spot 501 is generated in the D-mode image 103A having a relatively large number of data points. In the D-mode image 103B having a relatively small number of data points, vertical stripes 502 are generated.

そこで、画像処理部14は、Dモード画像103AとDモード画像103Bとが合成された合成Dモード画像103C(図4(c))の画像データを生成する。これにより、斑点501及び縦筋502が軽減された高画質のDモード画像が得られる。   Therefore, the image processing unit 14 generates image data of a combined D-mode image 103C (FIG. 4C) obtained by combining the D-mode image 103A and the D-mode image 103B. Thereby, a high-quality D-mode image in which the spots 501 and the vertical stripes 502 are reduced is obtained.

合成は、種々の方法により行ってよい。また、合成方法は、操作部16に対する操作等により種々の合成方法から選択可能であってもよい。例えば、Dモード画像103A及びDモード画像103Bの互いに対応する画素間において輝度が平均された画像を合成Dモード画像103Cとしてよい。互いに対応する画素は、時間(横軸)及びドプラ偏移周波数等(縦軸)が互いに同一の画素である。   The synthesis may be performed by various methods. Further, the synthesis method may be selectable from various synthesis methods by operating the operation unit 16 or the like. For example, an image in which luminance is averaged between pixels corresponding to each other in the D-mode image 103A and the D-mode image 103B may be used as the combined D-mode image 103C. Pixels corresponding to each other are pixels having the same time (horizontal axis), Doppler shift frequency and the like (vertical axis).

すなわち、図4においてハッチングして示すように、Dモード画像103Aにおける時点t1及びドプラ偏移周波数等f1の画素PAと、Dモード画像103Bにおける時点t1及びドプラ偏移周波数等f1の画素PBとが互いに対応し、これらの画素PA及びPBの輝度の平均が、合成Dモード画像103Cにおける時点t1及びドプラ偏移周波数等f1の画素PCの輝度になるように合成してよい。   That is, as shown by hatching in FIG. 4, the pixel PA at the time point t1 and the Doppler shift frequency etc. f1 in the D-mode image 103A, and the pixel PB at the time point t1 and the Doppler shift frequency etc. f1 in the D-mode image 103B. Corresponding to each other, the average of the luminances of these pixels PA and PB may be combined so as to be the luminance of the pixel PC at f1 such as the time point t1 and the Doppler shift frequency in the combined D-mode image 103C.

また、例えば、Dモード画像103Aと、Dモード画像103Bとの互いに対応する画素間において大きいほうの輝度(合成されるDモード画像が3種以上の場合には最大の輝度)が選択された画像を合成Dモード画像103Cとしてよい。   In addition, for example, an image in which the higher luminance (the maximum luminance when there are three or more D-mode images to be combined) is selected between the corresponding pixels of the D-mode image 103A and the D-mode image 103B. May be the composite D-mode image 103C.

また、合成は、Dモード画像全体において行ってもよいし、Dモード画像の一部の領域に対してのみ行ってもよい。例えば、Dモード画像103Aにおいて斑点501が生じる領域505においてのみDモード画像103Bを合成してもよいし、Dモード画像103Bにおいて縦筋502が生じる領域506においてのみDモード画像103Aを合成してもよい。なお、この場合の合成も、上述のように画素毎に輝度の平均、又は、高い輝度の選択を行ってよい。   Further, the synthesis may be performed on the entire D-mode image, or may be performed only on a partial region of the D-mode image. For example, the D mode image 103B may be synthesized only in the region 505 where the spot 501 occurs in the D mode image 103A, or the D mode image 103A may be synthesized only in the region 506 where the vertical streak 502 occurs in the D mode image 103B. Good. In the synthesis in this case, the average luminance or the high luminance may be selected for each pixel as described above.

上述のような各種の合成を行うと、例えば、データ点数の少ないDモード画像103Bにおいて信号強度が弱い部分が縦方向に連続している(縦筋がある)場合、その部分は、データ点数の多いDモード画像103Aとの平均や置換により輝度が高くなり、縦筋が目立たなくなる。斑点についても同様である。そして、両画像が互いの欠点を補うことにより、高画質のDモード画像103Cが得られる。   When various combinations as described above are performed, for example, in a D-mode image 103B having a small number of data points, when a portion having a low signal intensity is continuous in the vertical direction (there is a vertical line), the portion has the number of data points. Luminance increases due to averaging and replacement with many D-mode images 103A, and vertical stripes become inconspicuous. The same applies to spots. Then, the two images compensate for each other's drawbacks, whereby a high-quality D-mode image 103C is obtained.

図5は、超音波診断装置1のDモード画像を生成する手順の一例を説明するフローチャートである。具体的には、Dモード画像103をリアルタイム表示する場合のスペクトルデータ生成部20の動作(ステップS1〜S5)及び画像処理部14の動作(ステップS6〜S10)の一例を示している。   FIG. 5 is a flowchart for explaining an example of a procedure for generating a D-mode image of the ultrasonic diagnostic apparatus 1. Specifically, an example of the operation of the spectrum data generation unit 20 (steps S1 to S5) and the operation of the image processing unit 14 (steps S6 to S10) when the D-mode image 103 is displayed in real time is shown.

超音波診断装置1は、時系列データ生成部19によりデータ点数N1の時系列データDAが新たに生成されか否か判定し(ステップS1)、データ点数N1の時系列データDAが新たに生成されたと判定すると、データ点数N1の時系列データDAに窓関数を掛け(ステップS2)、データ点数N1の時系列データDAをフーリエ変換し、ドプラスペクトルデータを生成する(ステップS3)。   The ultrasonic diagnostic apparatus 1 determines whether or not the time series data generation unit 19 newly generates the time series data DA having the data point N1 (step S1), and the time series data DA having the data point N1 is newly generated. If it is determined that the time series data DA of the data point number N1 is multiplied by a window function (step S2), the time series data DA of the data point number N1 is Fourier transformed to generate Doppler spectrum data (step S3).

窓関数は、フーリエ変換が有限区間において切り出された時系列データに対して行われることによって生じる誤差を低減するものである。窓関数には、方形窓、ハニング窓、ハミング窓、ガウス窓など、公知の適宜なものが選択されてよい。一般に、方形窓以外の窓関数は、時点データDCの時間軸に沿った順番を変数とする関数となっており、有限区間の中央(時間中心)付近において、関数の値(時系列データに乗じられる重み)が最大となる。   The window function reduces an error caused by performing Fourier transform on time series data cut out in a finite interval. As the window function, a known appropriate one such as a rectangular window, a Hanning window, a Hamming window, or a Gaussian window may be selected. In general, a window function other than a rectangular window is a function whose variable is the order along the time axis of the time point data DC, and the function value (time series data is multiplied) near the center (time center) of the finite interval. Weight) is maximized.

次に、超音波診断装置1は、ステップS2、S3と同様に、データ点数N2の時系列データDBに窓関数を掛け(ステップS4)、データ点数N2の時系列データDBをフーリエ変換し、ドプラスペクトルデータを生成する(ステップS5)。なお、ステップS4で使用される窓関数は、ステップS2で使用される窓関数と同一のものである。ただし、異なる窓関数であってもよい。   Next, as in steps S2 and S3, the ultrasound diagnostic apparatus 1 multiplies the time series data DB with the number of data points N2 by a window function (step S4), Fourier transforms the time series data DB with the number of data points N2, and performs Doppler. Spectral data is generated (step S5). Note that the window function used in step S4 is the same as the window function used in step S2. However, different window functions may be used.

その後、超音波診断装置1は、画像データを生成する時期か否かを判定する(ステップS6)。例えば、時系列データDAの取得が所定回数行われたか否か、所定のフレームレートで画面を更新すべき時期が到来したか否かを判定する。画像データ生成時期が到来していないと判定した場合は、ステップS1に戻る。画像データ生成時期が到来したと判定した場合は、Dモード画像103AとDモード画像103Bとを合成すべき領域(例えば領域505、506)を検出する(ステップS7)。なお、画像全体において合成を行う場合にはステップS7は省略される。   Thereafter, the ultrasonic diagnostic apparatus 1 determines whether it is time to generate image data (step S6). For example, it is determined whether or not the acquisition of the time series data DA has been performed a predetermined number of times and whether or not it is time to update the screen at a predetermined frame rate. If it is determined that the image data generation time has not arrived, the process returns to step S1. If it is determined that the image data generation time has arrived, areas (for example, areas 505 and 506) where the D-mode image 103A and the D-mode image 103B are to be combined are detected (step S7). Note that step S7 is omitted when the entire image is synthesized.

領域505、506の検出は、例えば、ドプラスペクトルデータの各ドプラ偏移周波数等のパワー等と、そのドプラスペクトルデータに対して時間的に前後する1〜複数のドプラスペクトルデータの対応するドプラ偏移周波数等のパワー等とを比較し、急激な変化があった場合に斑点501や縦筋502が存在すると判定することにより行ったり、生成したドプラスペクトルデータと、予め用意された、所定の大きさの斑点501や縦筋502が生じた場合の基準ドプラスペクトルとパターンマッチングを行うことにより行う。   The detection of the regions 505 and 506 is performed by, for example, the power of each Doppler shift frequency of the Doppler spectrum data and the corresponding Doppler shifts of one to a plurality of Doppler spectrum data temporally before and after the Doppler spectrum data. Compared with power such as frequency and the like, and when there is a sudden change, it is determined that there are spots 501 or vertical stripes 502, or the generated Doppler spectrum data and a predetermined size prepared in advance This is performed by performing pattern matching with the reference Doppler spectrum when the spot 501 or the vertical streak 502 occurs.

ステップS8では、超音波診断装置は、データ点数N1の時系列データをフーリエ変換したスペクトラムデータと、データ点数N2の時系列データをフーリエ変換したスペクトラムデータとを合成する。すなわち、互いに同一の時間及びドプラ周波数等の2種類のパワー等を平均し、又は、2種類のパワー等のうち一方を選択し、合成ドプラスペクトルデータを得る。検出された合成領域505又は506についてのみ合成を行う場合には、2つのドプラスペクトルデータのうち一方のデータにおいて、合成領域に対応するデータについてのみ他方のデータを合成し、合成ドプラスペクトルデータを得る。   In step S8, the ultrasonic diagnostic apparatus synthesizes the spectrum data obtained by Fourier transforming the time series data having the data points N1 and the spectrum data obtained by performing the Fourier transform on the time series data having the data points N2. That is, two types of power such as the same time and Doppler frequency are averaged, or one of the two types of power is selected to obtain combined Doppler spectrum data. In the case of performing synthesis only for the detected synthesis region 505 or 506, in one data of the two Doppler spectrum data, the other data is synthesized only for the data corresponding to the synthesis region to obtain synthesized Doppler spectrum data. .

ステップS10では、超音波診断装置は、合成ドプラスペクトルデータに基づいて合成Dモード画像103Cの画像データを生成し(ステップS9)、生成した画像データに基づく映像信号を表示部15に出力して表示部15の画面を更新する(ステップS10)。   In step S10, the ultrasound diagnostic apparatus generates image data of the combined D-mode image 103C based on the combined Doppler spectrum data (step S9), and outputs and displays a video signal based on the generated image data on the display unit 15. The screen of the unit 15 is updated (step S10).

以上の実施形態によれば、データ点数N1の時系列データDAと、データ点数N2の時系列データDBとにフーリエ変換を行い2種類のドプラスペクトルデータを生成し、その2種類のドプラスペクトルデータに基づく2種類のDモード画像103A、103Bが合成された合成Dモード画像103Cの画像データを生成することから、データ点数が多いときに生じる斑点501やデータ点数が少ないときに生じる縦筋502のうち少なくとも一方が軽減され、画質が向上する。   According to the above embodiment, the Fourier transform is performed on the time series data DA having the data points N1 and the time series data DB having the data points N2, and two types of Doppler spectrum data are generated, and the two types of Doppler spectrum data are converted into the two types of Doppler spectrum data. Since the image data of the composite D-mode image 103C obtained by combining the two types of D-mode images 103A and 103B based on it is generated, out of the spots 501 generated when the number of data points is large and the vertical stripe 502 generated when the number of data points is small At least one of them is reduced and the image quality is improved.

例えば、合成Dモード画像103Cが、Dモード画像103A、103Bの互いに対応する画素間において輝度が平均された画像である場合には、斑点501及び縦筋502の双方が低減される。   For example, when the synthesized D-mode image 103C is an image in which the luminance is averaged between pixels corresponding to each other in the D-mode images 103A and 103B, both the spots 501 and the vertical stripes 502 are reduced.

例えば、合成Dモード画像103Cが、Dモード画像103A、103Bの互いに対応する画素間において最大の輝度が選択された画像であれば、斑点501の埋まりがよくなり、感度がよい印象をユーザに与えることができる。   For example, if the composite D-mode image 103C is an image in which the maximum luminance is selected between the corresponding pixels of the D-mode images 103A and 103B, the spots 501 are well filled, giving the user a good impression of sensitivity. be able to.

例えば、合成Dモード画像103Cが、Dモード画像103A、103Bのいずれか一のDモード画像の所定の領域において他の画像が合成された画像データであれば、画質の低下を招いている部分のみを合成して効果的に画質を向上させることができる。   For example, if the synthesized D-mode image 103C is image data obtained by synthesizing another image in a predetermined region of any one of the D-mode images 103A and 103B, only the portion that causes a reduction in image quality. Can effectively improve the image quality.

時系列データDBは、時系列データDAのうち一部の時間範囲の時系列データ全てであり、時系列データDAの時系列データDAのサンプリング周期と時系列データDBのサンプリング周期とは互いに同一である。従って、生成された2種類のドプラスペクトルのサンプリング周波数及び最大周波数は互いに同一であり、データ点数の少ないドプラスペクトルデータの形は、データ点数の多いドプラスペクトルデータから低周波数側のデータを削除した形になる。このため、2種類のドプラスペクトルにおいてドプラ偏移周波数等が互いに同一のデータ点(パワー等)を特定することが容易になり、合成のアルゴリズムが簡潔になる。   The time series data DB is all the time series data in a part of the time range of the time series data DA, and the sampling period of the time series data DA of the time series data DA and the sampling period of the time series data DB are the same. is there. Accordingly, the sampling frequency and the maximum frequency of the two types of generated Doppler spectra are the same, and the Doppler spectrum data having a small number of data points is obtained by deleting the low frequency data from the Doppler spectrum data having a large number of data points. become. Therefore, it becomes easy to specify data points (power, etc.) having the same Doppler shift frequency in the two types of Doppler spectra, and the synthesis algorithm is simplified.

フーリエ変換の前には時系列データDA、DBに窓関数が掛けられ、時系列データDBは、時系列データDAの時間中心と時系列データDBの時間中心とが一致するように選択されている。一方、上述のように、窓関数は一般に中央付近の値が大きい。従って、時系列データDAに掛け合わされる窓関数と、時系列データDBに掛け合わされる窓関数とは互いにピークとなる時間が略一致する。その結果、窓関数の影響による2つのドプラスペクトル間の誤差は縮小され、合成によって新たなノイズが生じることが抑制される。   Before the Fourier transform, the time series data DA and DB are multiplied by a window function, and the time series data DB is selected so that the time center of the time series data DA and the time center of the time series data DB coincide. . On the other hand, as described above, the window function generally has a large value near the center. Accordingly, the window function multiplied by the time series data DA and the window function multiplied by the time series data DB substantially coincide with each other at the peak time. As a result, the error between the two Doppler spectra due to the influence of the window function is reduced, and generation of new noise due to synthesis is suppressed.

なお、以上の実施形態において、画像処理部14は本発明の画像生成手段の一例であり、データ処理部13、画像処理部14及び表示部15の組合せは本発明の画像表示装置の一例である。   In the above embodiment, the image processing unit 14 is an example of the image generation unit of the present invention, and the combination of the data processing unit 13, the image processing unit 14, and the display unit 15 is an example of the image display device of the present invention. .

本発明は、以上の実施形態に限定されず、種々の態様で実施してよい。   The present invention is not limited to the above embodiment, and may be implemented in various aspects.

超音波診断装置は、Dモード画像に資する時系列データを得ることができるものであればよく、パルスドプラ法により時系列データを得るものに限定されない。例えば、連続ドプラ法により時系列データを得るものであってもよい。   The ultrasonic diagnostic apparatus may be any apparatus that can obtain time-series data contributing to the D-mode image, and is not limited to an apparatus that obtains time-series data by the pulse Doppler method. For example, time series data may be obtained by a continuous Doppler method.

複数種類のドプラスペクトルデータに基づく複数種類のDモード画像が合成された合成Dモード画像の画像データの生成は、結果として合成Dモード画像の画像データが生成されればよく、実施形態に示したように、複数種類のドプラスペクトルデータを合成してから画像データを生成することにより行われるものに限定されない。例えば、複数種類のドプラスペクトルデータそれぞれに基づいてDモード画像の画像データを複数生成し、複数のDモード画像の画像データを合成することにより、合成Dモード画像の画像データを生成してもよい。また、このような場合、画像生成部は、複数種類のDモード画像の互いに同一の座標の画素を、互いに対応する画素として検出してもよい。この場合であっても、各画素の時間及びドプラ偏移周波数等が対応していることに変わりはない。Dモード画像の拡大縮小により画素数とデータ数とが一致していない場合も同様である。   The generation of the image data of the combined D-mode image obtained by combining the plurality of types of D-mode images based on the plurality of types of Doppler spectrum data is sufficient if the image data of the combined D-mode image is generated as a result. Thus, it is not limited to what is performed by combining a plurality of types of Doppler spectrum data and then generating image data. For example, a plurality of D-mode image data may be generated based on each of a plurality of types of Doppler spectrum data, and the combined D-mode image data may be generated by combining the plurality of D-mode image data. . In such a case, the image generation unit may detect pixels having the same coordinates in a plurality of types of D-mode images as corresponding pixels. Even in this case, the time of each pixel, the Doppler shift frequency, and the like still correspond. The same applies when the number of pixels does not match the number of data due to the enlargement / reduction of the D-mode image.

合成は、画素毎に複数のDモード画像の輝度(パワー等)を平均、選択するものに限定されない。例えば、一の画像の所定の領域を他の画像の対応領域に置き換えてもよい。複数のDモード画像の互いに対応する時間毎又はドプラ偏移周波数等毎に、輝度を比較して選択するなどしてもよい。3種類以上のDモード画像の互いに対応する画素の輝度を、重み係数を乗じて加算してもよい。合成が所定の領域においてのみ行われる場合、所定の領域は、実施形態のように超音波診断装置によって検出されるものに限定されず、ユーザによって設定されるものであってもよい。   The composition is not limited to the one that averages and selects the luminance (power, etc.) of a plurality of D-mode images for each pixel. For example, a predetermined area of one image may be replaced with a corresponding area of another image. You may compare and select a brightness | luminance for every time or a Doppler shift frequency etc. to which a some D mode image mutually respond | corresponds. The luminance of pixels corresponding to each other in three or more types of D-mode images may be multiplied by a weighting factor and added. When combining is performed only in a predetermined region, the predetermined region is not limited to that detected by the ultrasonic diagnostic apparatus as in the embodiment, and may be set by the user.

互いに異なるデータ点数の時系列データは、所定の時間範囲を対象とするフーリエ変換が所定の時間範囲をずらしつつ順次行われるように設定されていればよく、図3に示したものに限定されない。   The time-series data having different data points is not limited to that shown in FIG. 3 as long as Fourier transform for a predetermined time range is sequentially performed while shifting the predetermined time range.

例えば、図6(a)に示すように、所定の時間範囲T21(データ点数の多い時系列データDA11〜DA14)が互いに一部重複していてもよい。さらに、図示はしないが、データ点数の少ない時系列データDB11〜DB14が互いに一部重複していてもよい。逆に、図6(b)に示すように、所定の時間範囲T21(データ点数の多い時系列データDA21〜DA23)が互いに離間していてもよい。なお、所定の時間範囲と、その所定の時間範囲内においてデータ点数が最も多い時系列データの時間範囲とは一致している必要はない。   For example, as shown in FIG. 6A, a predetermined time range T21 (time series data DA11 to DA14 having a large number of data points) may partially overlap each other. Furthermore, although not shown, the time series data DB11 to DB14 with a small number of data points may partially overlap each other. Conversely, as shown in FIG. 6B, the predetermined time range T21 (time series data DA21 to DA23 having a large number of data points) may be separated from each other. Note that the predetermined time range and the time range of the time-series data having the largest number of data points in the predetermined time range do not have to coincide with each other.

また、例えば、図6(c)に示すように、データ点数の多い時系列データDA1〜DA3の時間中心Tacと、データ点数の少ない時系列データDB31〜DB33の時間中心Tbc′とが一致していなくてもよい。   Further, for example, as shown in FIG. 6C, the time center Tac of the time series data DA1 to DA3 having a large number of data points coincides with the time center Tbc ′ of the time series data DB31 to DB33 having a small number of data points. It does not have to be.

また、例えば、図6(d)に示すように、データ点数の多い時系列データDA1〜DA3のサンプリング周期T1と、データ点数の少ない時系列データDB41〜DB43のサンプリング周期T3(ハッチングして示すもののみが時系列データDB41〜DB43を構成する)とは一致していなくてもよい。   Further, for example, as shown in FIG. 6D, the sampling period T1 of the time series data DA1 to DA3 having a large number of data points and the sampling period T3 of the time series data DB41 to DB43 having a small number of data points (shown by hatching). Only the time-series data DB41 to DB43) may not coincide with each other.

本発明の実施形態に係る超音波診断装置の構成の要部を示すブロック図。The block diagram which shows the principal part of the structure of the ultrasound diagnosing device which concerns on embodiment of this invention. 図2の超音波診断装置の表示する画像の一例を示す図。The figure which shows an example of the image which the ultrasonic diagnostic apparatus of FIG. 2 displays. 図1の超音波診断装置のデータ処理部の動作を説明する概念図。The conceptual diagram explaining operation | movement of the data processing part of the ultrasonic diagnosing device of FIG. 図1の超音波診断装置の画像処理部の動作を説明する概念図。The conceptual diagram explaining operation | movement of the image process part of the ultrasonic diagnosing device of FIG. 図1の超音波診断装置のDモード画像を生成する手順を説明するフローチャート。The flowchart explaining the procedure which produces | generates D mode image of the ultrasound diagnosing device of FIG. 図1の超音波診断装置のデータ処理方法の変形例を示す図。The figure which shows the modification of the data processing method of the ultrasonic diagnosing device of FIG. 従来のDモード画像を説明する図。The figure explaining the conventional D mode image.

符号の説明Explanation of symbols

1…超音波診断装置、11…プローブ、19…時系列データ生成部、20…スペクトルデータ生成部、14…画像処理部、15…表示部、DA,DB…時系列データ、103C…合成Dモード画像。 DESCRIPTION OF SYMBOLS 1 ... Ultrasound diagnostic apparatus, 11 ... Probe, 19 ... Time series data generation part, 20 ... Spectral data generation part, 14 ... Image processing part, 15 ... Display part, DA, DB ... Time series data, 103C ... Composite D mode image.

Claims (16)

超音波を被検体に照射し、反射された超音波を電気信号に変換して出力するプローブと、
前記プローブからの信号に基づいてドプラ情報を含む時系列データを生成する時系列データ生成手段と、
前記時系列データのうち所定の時間範囲の時系列データにフーリエ変換を行いドプラスペクトルデータを生成することを、前記所定の時間範囲を一定時間ずつずらして複数回行い、複数時の前記ドプラスペクトルデータを生成するスペクトルデータ生成手段と、
前記複数時のドプラスペクトルデータに基づいてDモード画像の画像データを生成する画像生成手段と、
前記画像データに基づいてDモード画像を表示する表示手段と、
を備え、
前記スペクトルデータ生成手段は、前記所定の時間範囲の時系列データのうち所定のデータ点数の時系列データにフーリエ変換を行いドプラスペクトルデータを生成することを、前記データ点数を異ならせて複数回行い、複数種類の前記ドプラスペクトルデータを生成し、
前記画像生成手段は、前記複数種類のドプラスペクトルデータに基づく複数種類のDモード画像が合成された合成Dモード画像の画像データを生成する
超音波診断装置。
A probe that irradiates a subject with ultrasonic waves, converts the reflected ultrasonic waves into electrical signals, and outputs a probe;
Time-series data generating means for generating time-series data including Doppler information based on a signal from the probe;
Generating Doppler spectrum data by performing Fourier transform on the time series data in a predetermined time range among the time series data is performed a plurality of times while shifting the predetermined time range by a predetermined time, and the Doppler spectrum data at a plurality of times Spectral data generating means for generating
Image generating means for generating image data of a D-mode image based on the multiple-time Doppler spectrum data;
Display means for displaying a D-mode image based on the image data;
With
The spectrum data generation means performs Fourier transform on the time series data of a predetermined number of data points among the time series data in the predetermined time range to generate Doppler spectrum data a plurality of times with different data points. Generating a plurality of types of the Doppler spectrum data,
The ultrasonic diagnostic apparatus, wherein the image generation unit generates image data of a combined D-mode image obtained by combining a plurality of types of D-mode images based on the plurality of types of Doppler spectrum data.
前記合成Dモード画像は、前記複数種類のDモード画像の互いに対応する画素間において輝度が平均された画像である
請求項1に記載の超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 1, wherein the composite D-mode image is an image in which luminance is averaged between pixels corresponding to each other of the plurality of types of D-mode images.
前記合成Dモード画像は、前記複数種類のDモード画像の互いに対応する画素間において最大の輝度が選択された画像である
請求項1に記載の超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 1, wherein the synthesized D-mode image is an image in which maximum luminance is selected between pixels corresponding to each other of the plurality of types of D-mode images.
前記合成Dモード画像は、前記複数種類のDモード画像のいずれか一のDモード画像の所定の領域において他の画像が合成された画像データである
請求項1〜3のいずれか1項に記載の超音波診断装置。
The composite D-mode image is image data obtained by combining another image in a predetermined region of any one of the plurality of types of D-mode images. Ultrasound diagnostic equipment.
前記所定の領域は、前記一のDモード画像において所定の大きさの斑点が生じる領域であり、
前記他の画像は、前記一のDモード画像よりも少ないデータ点数の前記時系列データに基づくDモード画像である
請求項4に記載の超音波診断装置。
The predetermined region is a region where a spot of a predetermined size occurs in the one D-mode image,
The ultrasonic diagnostic apparatus according to claim 4, wherein the other image is a D-mode image based on the time-series data having a smaller number of data points than the one D-mode image.
前記所定の領域は、前記一のDモード画像においてドプラ偏移周波数又はドプラ偏移周波数に相関する物理量を示す軸に沿う線が生じる領域であり、
前記他の画像は、前記一のDモード画像よりも多いデータ点数の前記時系列データに基づくDモード画像である
請求項4に記載の超音波診断装置。
The predetermined region is a region where a line along an axis indicating a Doppler shift frequency or a physical quantity correlated with the Doppler shift frequency is generated in the one D-mode image,
The ultrasonic diagnostic apparatus according to claim 4, wherein the other image is a D-mode image based on the time-series data having more data points than the one D-mode image.
前記所定の時間範囲において前記データ点数が互いに異なる時系列データは、一のデータ点数の時系列データのうち一部の時間範囲の時系列データ全てが他のデータ点数の時系列データとなるように、前記所定の時間範囲の時系列データから選択されている
請求項1〜6のいずれか1項に記載の超音波診断装置。
The time series data having different data points in the predetermined time range is such that all of the time series data in a part of the time range among the time series data of one data point becomes time series data of other data points. The ultrasonic diagnostic apparatus according to claim 1, wherein the ultrasonic diagnostic apparatus is selected from time-series data in the predetermined time range.
前記スペクトルデータ生成手段は、前記データ点数が互いに異なる時系列データにそれぞれ、時系列データの時間中心側において値が大きくなる窓関数を掛けてからフーリエ変換を行い、前記複数種類のドプラスペクトルデータを生成し、
前記所定の時間範囲において前記データ点数が互いに異なる時系列データは、各時系列データの時間中心が互いに同一になるように、前記所定の時間範囲の時系列データから選択されている
請求項1〜7のいずれか1項に記載の超音波診断装置。
The spectrum data generation means performs a Fourier transform after multiplying the time series data having different data points by a window function having a larger value on the time center side of the time series data, and converts the plurality of types of Doppler spectrum data into Generate
The time-series data having different data points in the predetermined time range is selected from the time-series data in the predetermined time range so that the time centers of the respective time-series data are the same. The ultrasonic diagnostic apparatus according to any one of 7.
ドプラ情報を含む時系列データを生成する時系列データ生成手段と、
前記時系列データのうち所定の時間範囲の時系列データにフーリエ変換を行いドプラスペクトルデータを生成することを、前記所定の時間範囲を一定時間ずつずらして複数回行い、複数時の前記ドプラスペクトルデータを生成するスペクトルデータ生成手段と、
前記複数時のドプラスペクトルデータに基づいてDモード画像の画像データを生成する画像生成手段と、
前記画像データに基づいてDモード画像を表示する表示手段と、
を備え、
前記スペクトルデータ生成手段は、前記所定の時間範囲の時系列データのうち所定のデータ点数の時系列データにフーリエ変換を行いドプラスペクトルデータを生成することを、前記データ点数を異ならせて複数回行い、複数種類の前記ドプラスペクトルデータを生成し、
前記画像生成手段は、前記複数種類のドプラスペクトルデータに基づく複数種類のDモード画像が合成された合成Dモード画像の画像データを生成する
画像表示装置。
Time series data generating means for generating time series data including Doppler information;
Generating Doppler spectrum data by performing Fourier transform on the time series data in a predetermined time range among the time series data is performed a plurality of times while shifting the predetermined time range by a predetermined time, and the Doppler spectrum data at a plurality of times Spectral data generating means for generating
Image generating means for generating image data of a D-mode image based on the multiple-time Doppler spectrum data;
Display means for displaying a D-mode image based on the image data;
With
The spectrum data generation means performs Fourier transform on the time series data of a predetermined number of data points among the time series data in the predetermined time range to generate Doppler spectrum data a plurality of times with different data points. Generating a plurality of types of the Doppler spectrum data,
The image generation unit generates image data of a combined D-mode image obtained by combining a plurality of types of D-mode images based on the plurality of types of Doppler spectrum data.
前記合成Dモード画像は、前記複数種類のDモード画像の互いに対応する画素間において輝度が平均された画像である
請求項9に記載の画像表示装置。
The image display device according to claim 9, wherein the composite D-mode image is an image in which luminance is averaged between pixels corresponding to each other of the plurality of types of D-mode images.
前記合成Dモード画像は、前記複数種類のDモード画像の互いに対応する画素間において最大の輝度が選択された画像である
請求項9に記載の画像表示装置。
The image display device according to claim 9, wherein the composite D-mode image is an image in which maximum luminance is selected between pixels corresponding to each other of the plurality of types of D-mode images.
前記合成Dモード画像は、前記複数種類のDモード画像のいずれか一のDモード画像の所定の領域において他の画像が合成された画像データである
請求項9〜11のいずれか1項に記載の画像表示装置。
The composite D-mode image is image data in which another image is synthesized in a predetermined region of any one of the plurality of types of D-mode images. Image display device.
前記所定の領域は、前記一のDモード画像において所定の大きさの斑点が生じる領域であり、
前記他の画像は、前記一のDモード画像よりも少ないデータ点数の前記時系列データに基づくDモード画像である
請求項12に記載の画像表示装置。
The predetermined region is a region where a spot of a predetermined size occurs in the one D-mode image,
The image display device according to claim 12, wherein the other image is a D-mode image based on the time-series data having a smaller number of data points than the one D-mode image.
前記所定の領域は、前記一のDモード画像においてドプラ偏移周波数又はドプラ偏移周波数に相関する物理量を示す軸に沿う線が生じる領域であり、
前記他の画像は、前記一のDモード画像よりも多いデータ点数の前記時系列データに基づくDモード画像である
請求項12に記載の画像表示装置。
The predetermined region is a region where a line along an axis indicating a Doppler shift frequency or a physical quantity correlated with the Doppler shift frequency is generated in the one D-mode image,
The image display device according to claim 12, wherein the other image is a D-mode image based on the time-series data having a larger number of data points than the one D-mode image.
前記所定の時間範囲において前記データ点数が互いに異なる時系列データは、一のデータ点数の時系列データのうち一部の時間範囲の時系列データ全てが他のデータ点数の時系列データとなるように、前記所定の時間範囲の時系列データから選択されている
請求項9〜14のいずれか1項に記載の画像表示装置。
The time series data having different data points in the predetermined time range is such that all of the time series data in a part of the time range among the time series data of one data point becomes time series data of other data points. The image display device according to claim 9, wherein the image display device is selected from time-series data in the predetermined time range.
前記スペクトルデータ生成手段は、前記データ点数が互いに異なる時系列データにそれぞれ、時系列データの時間中心側において値が大きくなる窓関数を掛けてからフーリエ変換を行い、前記複数種類のドプラスペクトルデータを生成し、
前記所定の時間範囲において前記データ点数が互いに異なる時系列データは、各時系列データの時間中心が互いに同一になるように、前記所定の時間範囲の時系列データから選択されている
請求項9〜15のいずれか1項に記載の画像表示装置。
The spectrum data generation means performs a Fourier transform after multiplying the time series data having different data points by a window function having a larger value on the time center side of the time series data, and converts the plurality of types of Doppler spectrum data into Generate
The time-series data having different data points in the predetermined time range is selected from the time-series data in the predetermined time range so that the time centers of the time-series data are the same. The image display device according to any one of 15.
JP2006221165A 2006-08-14 2006-08-14 Ultrasonic diagnostic apparatus and image display apparatus Active JP4832211B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2006221165A JP4832211B2 (en) 2006-08-14 2006-08-14 Ultrasonic diagnostic apparatus and image display apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006221165A JP4832211B2 (en) 2006-08-14 2006-08-14 Ultrasonic diagnostic apparatus and image display apparatus

Publications (2)

Publication Number Publication Date
JP2008043495A JP2008043495A (en) 2008-02-28
JP4832211B2 true JP4832211B2 (en) 2011-12-07

Family

ID=39177818

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006221165A Active JP4832211B2 (en) 2006-08-14 2006-08-14 Ultrasonic diagnostic apparatus and image display apparatus

Country Status (1)

Country Link
JP (1) JP4832211B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107920802B (en) 2015-08-27 2022-02-25 皇家飞利浦有限公司 Spectral doppler processing with adaptive sampling window size

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61249455A (en) * 1985-04-26 1986-11-06 横河メディカルシステム株式会社 Ultrasonic doppler apparatus
JPH02195948A (en) * 1989-01-26 1990-08-02 Toshiba Corp Ultrasonic doppler diagnostic device
JPH034843A (en) * 1989-06-01 1991-01-10 Toshiba Corp Ultrasonic diagnostic device
JP2916219B2 (en) * 1990-06-19 1999-07-05 富士通株式会社 Ultrasound diagnostic equipment
JP3088586B2 (en) * 1993-05-17 2000-09-18 ジーイー横河メディカルシステム株式会社 Ultrasound diagnostic equipment
JP3446148B2 (en) * 1994-04-11 2003-09-16 株式会社日立メディコ Ultrasound Doppler diagnostic equipment
JP2004351021A (en) * 2003-05-30 2004-12-16 Ge Medical Systems Global Technology Co Llc Spectrum generating method, and ultrasonic diagnostic instrument

Also Published As

Publication number Publication date
JP2008043495A (en) 2008-02-28

Similar Documents

Publication Publication Date Title
JP4528529B2 (en) Ultrasonic diagnostic apparatus and ultrasonic image data processing method
US9814447B2 (en) Ultrasonic diagnostic apparatus
CN103300886B (en) Diagnostic ultrasound equipment and control method thereof
US10524767B2 (en) Ultrasound diagnostic apparatus and ultrasound imaging method
JP2006116149A (en) Ultrasonic diagnostic apparatus and control program thereof
JP2020110363A (en) Information processing device, information processing method, and program
JP2010119630A (en) Ultrasonograph
JPH05220138A (en) Ultrasonic diagnostic device
JP2005168585A (en) Ultrasonic doppler diagnostic apparatus
JP2019526370A (en) Phase deviation correction in ultrasonic shear wave elastography, related apparatus, system and method
JP2009268734A (en) Ultrasound observation apparatus
CN106691502B (en) Ultrasound system and method for generating elastic images
US20210298721A1 (en) Ultrasound diagnosis apparatus
JP4832211B2 (en) Ultrasonic diagnostic apparatus and image display apparatus
JPH10286256A (en) Ultrasonic diagnostic device and filter
JP2012010875A5 (en) Ultrasonic diagnostic apparatus and method for operating the same
JP5455567B2 (en) Ultrasonic diagnostic equipment
JP2005087730A (en) Simultaneous generation of spatial composite image and non-composite image
JP2008188288A (en) Ultrasonic diagnostic equipment and ultrasonic image display device
US11896435B2 (en) Ultrasound diagnostic apparatus and examination method
JP3943653B2 (en) Ultrasonic diagnostic equipment
JP4008251B2 (en) Doppler gain control device and ultrasonic imaging device
JP2008142130A (en) Ultrasonic diagnostic apparatus and its control processing program
JP6334883B2 (en) Ultrasonic diagnostic apparatus and display control program
JP2011182887A (en) Ultrasonograph

Legal Events

Date Code Title Description
A625 Written request for application examination (by other person)

Free format text: JAPANESE INTERMEDIATE CODE: A625

Effective date: 20090508

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20110727

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

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

R150 Certificate of patent or registration of utility model

Ref document number: 4832211

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20140930

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20140930

Year of fee payment: 3

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

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

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

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250