JP2011200282A - Ultrasonic diagnostic apparatus and control program for the same - Google Patents

Ultrasonic diagnostic apparatus and control program for the same Download PDF

Info

Publication number
JP2011200282A
JP2011200282A JP2010067702A JP2010067702A JP2011200282A JP 2011200282 A JP2011200282 A JP 2011200282A JP 2010067702 A JP2010067702 A JP 2010067702A JP 2010067702 A JP2010067702 A JP 2010067702A JP 2011200282 A JP2011200282 A JP 2011200282A
Authority
JP
Japan
Prior art keywords
elasticity
data
region
contour
diagnostic apparatus
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.)
Granted
Application number
JP2010067702A
Other languages
Japanese (ja)
Other versions
JP2011200282A5 (en
JP5624345B2 (en
Inventor
Haruji Funaya
晴二 船矢
Mayumi Ito
真由美 伊藤
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 JP2010067702A priority Critical patent/JP5624345B2/en
Publication of JP2011200282A publication Critical patent/JP2011200282A/en
Publication of JP2011200282A5 publication Critical patent/JP2011200282A5/ja
Application granted granted Critical
Publication of JP5624345B2 publication Critical patent/JP5624345B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

PROBLEM TO BE SOLVED: To provide an ultrasonic diagnostic apparatus which accurately extracts the contour of a region having different elasticity within a living tissue.SOLUTION: The ultrasonic diagnostic apparatus comprises a physical quantity data-forming unit which forms physical amount data of a living tissue on the basis of two temporary different echo signals on the same sound ray acquired by transmitting and receiving ultrasonic wave to the living tissue, a display unit which displays an elastic image EG on the basis of the physical amount data formed by the physical amount data-forming unit, and a contour extractor which calculates a threshold of the elastic data on the basis of change rate of the physical amount data between adjacent pixels and extracts the contour of a region C of a tumor by comparing the threshold with the elastic data of each pixel.

Description

本発明は、生体組織の硬さ又は軟らかさを表す弾性画像を表示する超音波診断装置及びその制御プログラムに関する。   The present invention relates to an ultrasonic diagnostic apparatus that displays an elastic image representing the hardness or softness of a living tissue and a control program therefor.

通常のBモード画像と、生体組織の硬さ又は軟らかさを表す弾性画像とを合成して表示させる超音波診断装置が、例えば特許文献1などに開示されている。この種の超音波診断装置において、弾性画像は次のようにして作成される。先ず、超音波プローブを体表面に押し当てて圧迫とその弛緩を繰り返すなどして生体組織を変形させながら超音波の送受信を行い、エコー信号を取得する。そして、得られたエコー信号に基づいて、生体組織の各部について弾性に関する物理量を算出し、この物理量を色相情報に変換してカラーの弾性画像を作成する。ちなみに、生体組織の弾性に関する物理量としては、例えば生体組織の歪みなどを算出している。   For example, Patent Literature 1 discloses an ultrasonic diagnostic apparatus that synthesizes and displays a normal B-mode image and an elastic image representing the hardness or softness of a living tissue. In this type of ultrasonic diagnostic apparatus, the elasticity image is created as follows. First, an ultrasonic probe is pressed against the body surface, and compression and relaxation thereof are repeated, and ultrasonic waves are transmitted and received while deforming a living tissue, thereby obtaining an echo signal. Based on the obtained echo signal, a physical quantity related to elasticity is calculated for each part of the living tissue, and the physical quantity is converted into hue information to create a color elastic image. Incidentally, as a physical quantity related to the elasticity of the living tissue, for example, a strain of the living tissue is calculated.

また、前記特許文献1では、生体組織の各部について算出した歪みと閾値とを比較して、弾性が異なる領域を抽出し、抽出された領域の輪郭を弾性画像に重ねて表示している。   Moreover, in the said patent document 1, the distortion calculated about each part of a biological tissue and a threshold value are compared, the area | region where elasticity differs is extracted, and the outline of the extracted area | region is superimposed and displayed on an elasticity image.

特開2004−135929号公報JP 2004-135929 A

ところで、前記超音波プローブによる圧迫力の強さにより、生体組織の変形量は変化する。従って、前記特許文献1のように、生体組織の各部について算出した歪みと閾値とを比較して、弾性が異なる領域を抽出する場合、生体組織への圧迫力の強さに応じた前記閾値を設定する必要があった。   By the way, the amount of deformation of the living tissue changes depending on the strength of the compression force by the ultrasonic probe. Therefore, as in Patent Document 1, when comparing the strain calculated for each part of the biological tissue with the threshold and extracting a region having different elasticity, the threshold corresponding to the strength of the compression force on the biological tissue is set. It was necessary to set.

上述の課題を解決するためになされた第1の観点の発明は、生体組織に対する超音波の送受信によって得られた同一音線上の時間的に異なる二つのエコー信号に基づいて、生体組織の弾性データを作成する弾性データ作成部と、この弾性データ作成部によって作成された弾性データに基づく弾性画像を表示する表示部と、隣り合う画素における弾性データの変化率に基づいて弾性データの閾値を算出し、この閾値と各画素の弾性データとを比較して、生体組織における異なる弾性を有する領域の輪郭を抽出する輪郭抽出部と、を備えることを特徴とする超音波診断装置である。   The invention of the first aspect made to solve the above-mentioned problem is that the elasticity data of the living tissue is based on two temporally different echo signals on the same sound ray obtained by transmitting and receiving ultrasonic waves to the living tissue. The elasticity data creation unit for creating the elastic data, the display unit for displaying the elasticity image based on the elasticity data created by the elasticity data creation unit, and the elasticity data threshold value based on the rate of change of the elasticity data in the adjacent pixels. An ultrasonic diagnostic apparatus comprising: a contour extracting unit that compares the threshold value and elasticity data of each pixel and extracts a contour of a region having different elasticity in a living tissue.

第2の観点の発明は、第1の観点の発明において、前記輪郭抽出部は、前記弾性画像に複数のサンプリングラインを設定し、各サンプリングライン上において、隣り合う画素間における弾性データの変化率が最大の部分を検出し、この最大の部分における隣り合う画素のうちいずれか一方の画素についての弾性データに基づいて前記閾値を算出することを特徴とする超音波診断装置である。   According to a second aspect of the invention, in the first aspect of the invention, the contour extraction unit sets a plurality of sampling lines in the elastic image, and a change rate of elastic data between adjacent pixels on each sampling line. Is an ultrasonic diagnostic apparatus that detects the maximum portion and calculates the threshold value based on elasticity data for any one of adjacent pixels in the maximum portion.

第3の観点の発明は、第2の観点の発明において、前記輪郭抽出部は、前記最大の部分における隣り合う画素のうち、注目領域の弾性により近い弾性を有する画素の弾性データについて、前記各サンプリングラインにおける平均を算出して、この平均を前記閾値とすることを特徴とする超音波診断装置である。   According to a third aspect of the present invention, in the second aspect of the invention, the contour extraction unit may determine the elasticity data of pixels having elasticity closer to the elasticity of the region of interest among adjacent pixels in the maximum portion. An ultrasonic diagnostic apparatus characterized in that an average in a sampling line is calculated and the average is used as the threshold value.

第4の観点の発明は、第3の観点の発明において、前記注目領域の弾性により近い弾性を有する画素は、前記最大の部分における隣り合う画素のうち、より硬い生体組織に対応する画素であることを特徴とする超音波診断装置である。   According to a fourth aspect of the invention, in the third aspect of the invention, the pixel having elasticity closer to the elasticity of the region of interest is a pixel corresponding to a harder biological tissue among adjacent pixels in the maximum portion. This is an ultrasonic diagnostic apparatus.

第5の観点の発明は、第2〜4のいずれか一の観点の発明において、前記輪郭抽出部は、前記弾性画像に設定された所定の領域内に前記サンプリングラインを設定することを特徴とする超音波診断装置である。   According to a fifth aspect of the invention, in the invention according to any one of the second to fourth aspects, the contour extraction unit sets the sampling line within a predetermined region set in the elastic image. This is an ultrasonic diagnostic apparatus.

第6の観点の発明は、第1〜5のいずれか一の観点の発明において、前記弾性データは、前記二つのエコー信号に基づいて算出された生体組織の弾性に関する物理量データであることを特徴とする超音波診断装置である。   The invention according to a sixth aspect is the invention according to any one of the first to fifth aspects, wherein the elasticity data is physical quantity data relating to elasticity of a biological tissue calculated based on the two echo signals. This is an ultrasonic diagnostic apparatus.

第7の観点の発明は、第1〜5のいずれか一の観点の発明において、前記弾性データは、前記二つのエコー信号に基づいて算出された生体組織の弾性に関する物理量データを階調化して得られた弾性画像データであることを特徴とする超音波診断装置である。   The invention according to a seventh aspect is the invention according to any one of the first to fifth aspects, wherein the elasticity data is obtained by gradationizing physical quantity data relating to elasticity of a biological tissue calculated based on the two echo signals. An ultrasonic diagnostic apparatus characterized by the obtained elastic image data.

第8の観点の発明は、第1〜7のいずれか一の観点の発明において、前記輪郭抽出部によって抽出された輪郭で囲まれる領域を示す表示を前記弾性画像上に表示させる表示設定部を備えることを特徴とする超音波診断装置である。   The invention according to an eighth aspect is the invention according to any one of the first to seventh aspects, wherein the display setting unit is configured to display a display indicating a region surrounded by the contour extracted by the contour extraction unit on the elastic image. An ultrasonic diagnostic apparatus comprising:

第9の観点の発明は、第1〜8のいずれか一の観点の発明において、前記輪郭抽出部によって抽出された輪郭で囲まれる領域の面積又は体積を算出する算出部を備えることを特徴とする超音波診断装置である。   A ninth aspect of the invention is the invention of any one of the first to eighth aspects, further comprising a calculation unit that calculates an area or volume of a region surrounded by the contour extracted by the contour extraction unit. This is an ultrasonic diagnostic apparatus.

第10の観点の発明は、コンピュータに、生体組織に対する超音波の送受信によって得られた同一音線上の時間的に異なる二つのエコー信号に基づいて、生体組織の弾性データを作成する弾性データ作成機能と、隣り合う画素における弾性データの変化率に基づいて弾性データの閾値を算出し、該閾値と各画素の弾性データとを比較して、生体組織における異なる弾性を有する領域の輪郭を抽出する輪郭抽出機能と、を実行させることを特徴とする超音波診断装置の制御プログラムである。   The tenth aspect of the invention is an elasticity data creation function for creating elasticity data of a living tissue based on two echo signals different in time on the same sound line obtained by transmitting and receiving ultrasonic waves to the living tissue. And calculating a threshold value of elasticity data based on a rate of change of elasticity data in adjacent pixels, comparing the threshold value with the elasticity data of each pixel, and extracting a contour of a region having different elasticity in a living tissue A control program for an ultrasonic diagnostic apparatus, characterized by executing an extraction function.

上記観点の発明によれば、隣り合う画素における弾性データの変化率に基づいて設定された閾値と各画素の弾性データとを比較して、生体組織における異なる弾性を有する領域の輪郭が抽出される。このように前記閾値は、隣り合う画素における弾性データの変化率に基づいて算出されるので、生体組織の変形量が異なっていても、生体組織における異なる弾性を有する領域の輪郭を正確に抽出することができる。   According to the invention of the above aspect, the threshold value set based on the change rate of the elasticity data in the adjacent pixels is compared with the elasticity data of each pixel, and the contour of the region having different elasticity in the living tissue is extracted. . As described above, the threshold value is calculated based on the rate of change of the elasticity data in adjacent pixels, so that even if the deformation amount of the living tissue is different, the contour of the region having different elasticity in the living tissue is accurately extracted. be able to.

本発明に係る超音波診断装置の実施形態の概略構成の一例を示すブロック図である。It is a block diagram which shows an example of schematic structure of embodiment of the ultrasonic diagnosing device which concerns on this invention. 物理量データの作成の説明図である。It is explanatory drawing of creation of physical quantity data. 図1に示す超音波診断装置における表示制御部の構成を示すブロック図である。It is a block diagram which shows the structure of the display control part in the ultrasonic diagnosing device shown in FIG. 図1に示す超音波診断装置における表示部の表示の一例を示す図である。It is a figure which shows an example of the display of the display part in the ultrasonic diagnosing device shown in FIG. 弾性画像に輪郭線を表示させる際のフローチャートを示す図である。It is a figure which shows the flowchart at the time of displaying an outline on an elasticity image. 図5に示すフローチャートにおける閾値の算出についてのフローチャートを示す図である。It is a figure which shows the flowchart about calculation of the threshold value in the flowchart shown in FIG. 弾性画像に所定の領域が設定された表示部を示す図である。It is a figure which shows the display part by which the predetermined area | region was set to the elasticity image. 所定の領域内へのサンプリングラインの設定を説明する図である。It is a figure explaining the setting of the sampling line in a predetermined area | region. サンプリングライン上における隣り合う画素間で歪みの変化率が最大である部分を説明する図である。It is a figure explaining the part where the change rate of distortion is the maximum between the adjacent pixels on a sampling line. サンプリングライン上における画素を示す概念図である。It is a conceptual diagram which shows the pixel on a sampling line. 斜め方向のサンプリングラインにおける隣り合う画素を説明する図である。It is a figure explaining the adjacent pixel in the sampling line of a diagonal direction. 腫瘍の領域の輪郭の抽出を説明する図である。It is a figure explaining the extraction of the outline of the area | region of a tumor. 弾性画像に輪郭線が表示された表示部を示す図である。It is a figure which shows the display part by which the outline was displayed on the elasticity image. 互いに直交する三断面を示す概念図である。It is a conceptual diagram which shows three cross sections orthogonal to each other. 第一変形例における表示部を示す図である。It is a figure which shows the display part in a 1st modification. 二次元の超音波画像に二次元の領域が設定された表示部を示す図である。It is a figure which shows the display part by which the two-dimensional area | region was set to the two-dimensional ultrasonic image. 所定の断面について二次元の領域が設定されたことを説明する図である。It is a figure explaining that the two-dimensional area | region was set about the predetermined cross section. 三次元の所定の領域の設定を説明する図である。It is a figure explaining the setting of a three-dimensional predetermined area | region. 第二変形例の表示制御部の構成を示すブロック図である。It is a block diagram which shows the structure of the display control part of a 2nd modification.

以下、本発明の実施形態について図1〜図13に基づいて説明する。図1に示す超音波診断装置1は、超音波プローブ2、送受信部3、Bモードデータ作成部4、物理量データ作成部5、表示制御部6、表示部7、制御部8及び操作部9を備える。   Hereinafter, embodiments of the present invention will be described with reference to FIGS. 1 includes an ultrasonic probe 2, a transmission / reception unit 3, a B-mode data creation unit 4, a physical quantity data creation unit 5, a display control unit 6, a display unit 7, a control unit 8, and an operation unit 9. Prepare.

前記超音波プローブ2は、生体組織に対して超音波を送信しそのエコーを受信する。この超音波プローブ2における超音波の送受信面を体表面に当接させた状態で、例えば圧迫と弛緩を繰り返しながら超音波の送受信を行なって取得されたエコー信号に基づいて、後述のように弾性画像が作成される。   The ultrasonic probe 2 transmits an ultrasonic wave to a living tissue and receives an echo thereof. With the ultrasonic wave transmission / reception surface of the ultrasonic probe 2 in contact with the body surface, for example, based on an echo signal obtained by performing ultrasonic wave transmission / reception while repeating compression and relaxation, elasticity is obtained as described later. An image is created.

前記送受信部3は、前記超音波プローブ2を所定の走査条件で駆動させて音線毎の超音波の走査を行なう。また、送受信部3は、前記超音波プローブ2で受信したエコー信号について、整相加算処理等の信号処理を行なう。前記送受信部3で信号処理されたエコー信号は、前記Bモードデータ作成部4及び前記物理量データ作成部5に出力される。   The transmission / reception unit 3 drives the ultrasonic probe 2 under a predetermined scanning condition to perform ultrasonic scanning for each sound ray. The transmission / reception unit 3 performs signal processing such as phasing addition processing on the echo signal received by the ultrasonic probe 2. The echo signal signal-processed by the transmission / reception unit 3 is output to the B-mode data creation unit 4 and the physical quantity data creation unit 5.

ちなみに、前記送受信部3は、Bモード画像を作成するためのBモード画像用走査と、弾性画像を作成するための弾性画像用走査とを別に行なう。弾性画像用走査としては、被検体における弾性画像を作成する領域(弾性画像作成領域)において、同一音線上に少なくとも二回の走査を行なう。   Incidentally, the transmission / reception unit 3 separately performs a B-mode image scan for creating a B-mode image and an elastic image scan for creating an elastic image. As elastic image scanning, at least two scans are performed on the same sound ray in a region (elastic image creation region) in which an elastic image is created in the subject.

前記Bモードデータ作成部4は、前記送受信部3から出力されたエコー信号に対し、対数圧縮処理、包絡線検波処理等のBモード処理を行い、Bモードデータを作成する。   The B-mode data creation unit 4 performs B-mode processing such as logarithmic compression processing and envelope detection processing on the echo signal output from the transmission / reception unit 3 to create B-mode data.

前記物理量データ作成部5は、前記送受信部3から出力されたエコー信号に基づいて、生体組織における各部の弾性に関する物理量データを作成する(物理量データ作成機能)。もう少し詳しく説明すると、この物理量データ作成部5は、生体組織における各部の弾性に関する物理量として、前記超音波プローブ2による圧迫とその弛緩などによって生じた生体組織における各部の歪みStを算出することにより前記物理量データを作成する。前記物理量データ作成部5は、図2に示すように時間的に異なる二つのフレーム(i),(ii)に属する同一音線上における二つのエコー信号に基づいて歪みStを算出し、歪みStのデータからなる物理量データを作成する。前記物理量データ作成部5は、本発明における弾性データ作成部の実施の形態の一例である。また、前記物理量データは、本発明における弾性データの実施の形態の一例である。さらに、前記物理量データ作成機能は、本発明における弾性データ作成機能の実施の形態の一例である。   The physical quantity data creation unit 5 creates physical quantity data related to the elasticity of each part in the living tissue based on the echo signal output from the transmission / reception unit 3 (physical quantity data creation function). More specifically, the physical quantity data creation unit 5 calculates the strain St of each part in the living tissue caused by the compression and relaxation of the ultrasonic probe 2 as a physical quantity related to the elasticity of each part in the living tissue. Create physical quantity data. The physical quantity data creation unit 5 calculates a distortion St based on two echo signals on the same sound ray belonging to two temporally different frames (i) and (ii) as shown in FIG. Create physical quantity data consisting of data. The physical quantity data creation unit 5 is an example of an embodiment of an elasticity data creation unit in the present invention. The physical quantity data is an example of an embodiment of elasticity data in the present invention. Furthermore, the physical quantity data creation function is an example of an embodiment of the elasticity data creation function in the present invention.

具体的には、前記物理量データ作成部5は、例えば、特開2008−126079号公報に記載されているように、前記二つのエコー信号に相関ウィンドウを設定し、各相関ウィンドウ間で複素相関関数の虚数部の演算を行なって得られた値を歪みStとする。前記特開2008−126079号公報では、前記二つのエコー信号の波形のずれを歪みと推定している。   Specifically, the physical quantity data creation unit 5 sets a correlation window for the two echo signals, as described in, for example, Japanese Patent Application Laid-Open No. 2008-126079, and a complex correlation function between the correlation windows. The value obtained by performing the calculation of the imaginary part is taken as distortion St. In the Japanese Patent Application Laid-Open No. 2008-126079, the deviation of the waveform of the two echo signals is estimated as distortion.

ちなみに、前記相関ウィンドウは、弾性画像の一画素に相当する。従って、前記二つのエコー信号に設定された一対の相関ウィンドウについて得られた歪みStのデータは一画素についてのデータである。   Incidentally, the correlation window corresponds to one pixel of the elastic image. Therefore, the distortion St data obtained for the pair of correlation windows set for the two echo signals is data for one pixel.

前記表示制御部6には、前記Bモードデータ作成部4からのBモードデータ及び前記物理量データ作成部5からの物理量データが入力されるようになっている。前記表示制御部6は、図3に示すように、Bモード画像データ作成部61、弾性画像データ作成部62、合成部63、輪郭抽出部64、輪郭抽出用領域設定部65及び輪郭線設定部66を有している。   The display control unit 6 is inputted with B mode data from the B mode data creation unit 4 and physical quantity data from the physical quantity data creation unit 5. As shown in FIG. 3, the display control unit 6 includes a B-mode image data creation unit 61, an elastic image data creation unit 62, a synthesis unit 63, a contour extraction unit 64, a contour extraction region setting unit 65, and a contour line setting unit. 66.

前記Bモード画像データ作成部61及び前記弾性画像データ作成部62は、スキャンコンバータ(scan converter)を有している。そして、前記Bモード画像データ作成部61は、前記Bモードデータを、エコーの信号強度に応じた輝度情報を有するBモード画像データに変換する。また、前記弾性画像データ作成部62は、前記物理量データを歪みに応じた色相情報を有するカラー弾性画像データに変換する(弾性画像作成機能)。   The B-mode image data creation unit 61 and the elastic image data creation unit 62 have a scan converter. The B-mode image data creation unit 61 converts the B-mode data into B-mode image data having luminance information corresponding to the echo signal intensity. The elastic image data creation unit 62 converts the physical quantity data into color elastic image data having hue information corresponding to distortion (elastic image creation function).

ちなみに、前記Bモード画像データにおける輝度情報及び前記カラー弾性画像データにおける色相情報は所定の階調(例えば256階調)からなる。すなわち前記Bモード画像データは前記Bモードデータを階調化して得られたデータであり、また前記カラー弾性画像データは前記物理量データを階調化して得られたデータである。   Incidentally, the luminance information in the B-mode image data and the hue information in the color elastic image data have predetermined gradations (for example, 256 gradations). That is, the B-mode image data is data obtained by gradation of the B-mode data, and the color elastic image data is data obtained by gradation of the physical quantity data.

なお、Bモード画像データに変換される前のBモードデータ及びカラー弾性画像データに変換される前の物理量データを、ローデータ(Raw Data)と云う。   The B-mode data before being converted into B-mode image data and the physical quantity data before being converted into color elastic image data are referred to as raw data.

前記合成部63は、前記Bモード画像データ及び前記カラー弾性画像データを加算処理することによって合成し、前記表示部7に表示する二次元の超音波画像の画像データを作成する。この画像データは、図4に示すように白黒のBモード画像BGとカラーの弾性画像EGとが合成された二次元の超音波画像Gとして前記表示部7に表示される。本例では、前記弾性画像EGは、関心領域R内に半透明で(背景のBモード画像が透けた状態で)表示される。   The synthesizing unit 63 synthesizes the B-mode image data and the color elastic image data by adding them, and creates image data of a two-dimensional ultrasonic image displayed on the display unit 7. The image data is displayed on the display unit 7 as a two-dimensional ultrasonic image G in which a monochrome B-mode image BG and a color elastic image EG are combined as shown in FIG. In this example, the elastic image EG is displayed in the region of interest R in a translucent manner (with the background B-mode image transparent).

前記輪郭抽出部64は、隣り合う画素における物理量データの変化率、すなわち隣り合う画素における歪みStの変化率に基づいて歪みStの閾値StTHを設定し、この閾値StTHと各画素の物理量データ(歪みStの値)とを比較して、生体組織における異なる弾性を有する領域の輪郭を抽出する。前記輪郭抽出部64は、本発明における輪郭抽出部の実施の形態の一例である。 The contour extraction unit 64 sets a threshold value St TH for the distortion St based on the rate of change of physical quantity data in adjacent pixels, that is, the rate of change of distortion St in adjacent pixels, and this threshold value St TH and the physical quantity data of each pixel. The contour of the region having different elasticity in the living tissue is extracted by comparing (the value of strain St). The contour extracting unit 64 is an example of an embodiment of a contour extracting unit in the present invention.

前記輪郭抽出部64は、本例では、生体組織における異なる弾性を有する領域として、正常な組織よりも硬い腫瘍の領域Cの輪郭を抽出する。詳細は後述する。前記腫瘍の領域Cは、本発明における注目領域の実施の形態の一例である。   In this example, the contour extraction unit 64 extracts the contour of the region C of the tumor that is harder than normal tissue as a region having different elasticity in the living tissue. Details will be described later. The tumor region C is an example of an embodiment of a region of interest in the present invention.

前記輪郭抽出用領域設定部65は、前記関心領域R内に表示された弾性画像EGに、所定の領域r(後述の図7等参照)を設定する。前記輪郭抽出用領域設定部64は、前記所定の領域rを前記操作部9からの入力に基づいて設定する。前記所定の領域rは、前記閾値StTHを算出して前記腫瘍の領域Cの輪郭を抽出する際に設定される。詳細は後述する。 The contour extraction region setting unit 65 sets a predetermined region r (see FIG. 7 described later) in the elastic image EG displayed in the region of interest R. The contour extraction region setting unit 64 sets the predetermined region r based on an input from the operation unit 9. The predetermined region r is set when the threshold value St TH is calculated and the contour of the tumor region C is extracted. Details will be described later.

前記輪郭線設定部66は、前記輪郭抽出部64によって抽出された腫瘍の領域Cの輪郭を示す輪郭線Oを前記弾性画像EGに表示する。詳細は後述する。前記輪郭線設定部66は、本発明における表示設定部の実施の形態の一例である。また、前記輪郭線Oは、本発明において、輪郭抽出部によって抽出された輪郭で囲まれる領域を示す表示の実施の形態の一例である。   The contour setting unit 66 displays a contour O indicating the contour of the tumor region C extracted by the contour extracting unit 64 on the elastic image EG. Details will be described later. The outline setting unit 66 is an example of an embodiment of a display setting unit in the present invention. Further, the contour line O is an example of an embodiment of a display showing a region surrounded by the contour extracted by the contour extracting unit in the present invention.

前記表示部7は、LCD(Liquid Crystal Display)やCRT(Cathode Ray Tube)などで構成される。前記表示部7は、本発明における表示部の実施の形態の一例である。   The display unit 7 includes an LCD (Liquid Crystal Display), a CRT (Cathode Ray Tube), or the like. The display unit 7 is an example of an embodiment of a display unit in the present invention.

前記制御部8は、CPU(Central Processing Unit)で構成され、図示しない記憶部に記憶された制御プログラムを読み出し、前記超音波診断装置1の各部における機能を実行させる。また、前記操作部9は、操作者が指示や情報を入力するためのキーボード及びポインティングデバイス(図示省略)などを含んで構成されている。   The control unit 8 is constituted by a CPU (Central Processing Unit), reads a control program stored in a storage unit (not shown), and executes functions in each unit of the ultrasonic diagnostic apparatus 1. The operation unit 9 includes a keyboard and a pointing device (not shown) for the operator to input instructions and information.

さて、本例の超音波診断装置1の作用について説明する。先ず、前記送受信部3は、前記超音波プローブ2から被検体の生体組織へ超音波を送信させ、そのエコー信号を取得する。このとき、前記超音波プローブ2により、例えば被検体への圧迫とその弛緩を繰り返すなどして生体組織を変形させながら超音波の送受信を行う。   Now, the operation of the ultrasonic diagnostic apparatus 1 of this example will be described. First, the transmission / reception unit 3 transmits an ultrasonic wave from the ultrasonic probe 2 to a living tissue of a subject and acquires an echo signal thereof. At this time, the ultrasonic probe 2 transmits and receives ultrasonic waves while deforming the living tissue by repeatedly pressing and relaxing the subject, for example.

エコー信号が取得されると、前記Bモードデータ作成部4がBモードデータを作成し、前記物理量データ作成部5が物理量データを作成する。そして、これらBモードデータ及び物理量データに基づいて、前記Bモード画像データ作成部61及び前記弾性画像データ作成部62がBモード画像データ及びカラー弾性画像データを作成し、これらを前記合成部63が合成して得られた画像データに基づく超音波画像Gが前記表示部7に表示される。   When the echo signal is acquired, the B-mode data creation unit 4 creates B-mode data, and the physical quantity data creation unit 5 creates physical quantity data. Based on the B-mode data and physical quantity data, the B-mode image data creation unit 61 and the elastic image data creation unit 62 create B-mode image data and color elastic image data. An ultrasonic image G based on the combined image data is displayed on the display unit 7.

前記超音波画像Gが表示されると、腫瘍の領域Cの輪郭抽出を行なった後に、前記超音波画像Gにおける弾性画像EGに前記輪郭線Oを表示させる。具体的な処理について、図5のフローチャートに基づいて説明する。   When the ultrasonic image G is displayed, the contour line O is displayed on the elastic image EG in the ultrasonic image G after extracting the contour of the tumor region C. Specific processing will be described based on the flowchart of FIG.

先ず、ステップS1では、前記輪郭抽出部64は、腫瘍の領域Cの輪郭抽出を行なう際に用いる歪みの閾値StTHを算出する。前記輪郭抽出部64は、前記弾性画像EGに複数のサンプリングラインSL(後述の図8参照)を設定し、各サンプリングラインSL上において、隣り合う画素間における歪みStの変化率が最大の部分を検出する。そして、この最大の部分における隣り合う画素のうちいずれか一方の画素についての歪みStのデータに基づいて前記閾値StTHを算出する。このような閾値StTHの算出について、図6のフローチャートに基づいて以下詳しく説明する。 First, in step S1, the contour extraction unit 64 calculates a distortion threshold St TH used when performing contour extraction of a tumor region C. The contour extraction unit 64 sets a plurality of sampling lines SL (see FIG. 8 described later) in the elastic image EG, and on each sampling line SL, a portion where the rate of change of the distortion St between adjacent pixels is the maximum. To detect. Then, the threshold value St TH is calculated based on the distortion St data for one of the adjacent pixels in the maximum portion. The calculation of the threshold value St TH will be described in detail below based on the flowchart of FIG.

図6のフローチャートにおいて、先ず、ステップS11では、前記輪郭抽出部64は、図7に示すように弾性画像EGにおいて、所定の領域rを設定する。この所定の領域rは、本例では四角形になっており、操作者が弾性画像EG(或いは弾性画像EGとBモード画像BGの両方)を参照して、腫瘍と思われる部分を含むように設定される。また、前記所定の領域rは、その中心が腫瘍と思われる部分の中心付近になるように設定される。   In the flowchart of FIG. 6, first, in step S11, the contour extraction unit 64 sets a predetermined region r in the elastic image EG as shown in FIG. The predetermined region r is a quadrangle in this example, and the operator refers to the elastic image EG (or both the elastic image EG and the B-mode image BG) and is set to include a portion that seems to be a tumor. Is done. Further, the predetermined region r is set so that the center thereof is near the center of the portion considered to be a tumor.

ちなみに、図7では説明の便宜上弾性画像EGのみが示されている(図8、図9、図13においても同様)。   Incidentally, only the elastic image EG is shown in FIG. 7 for convenience of explanation (the same applies to FIGS. 8, 9, and 13).

次にステップS12では、前記輪郭抽出部64は、図8に示すように前記所定の領域r内にサンプリングラインSLを設定する。ただし、このサンプリングラインSLは、弾性画像EG上に表示されるものではなく、仮想のラインである。   In step S12, the contour extraction unit 64 sets a sampling line SL in the predetermined region r as shown in FIG. However, this sampling line SL is not displayed on the elastic image EG, but is a virtual line.

本例では、前記サンプリングラインSLとして、前記所定の領域rの中心から放射状に8本のサンプリングラインSL1,SL2,SL3,SL4,SL5,SL6,SL7,SL8が設定される。前記サンプリングラインSL1,SL3,SL5,SL7は、前記所定の領域rの中心と、この所定の領域rの各辺の中点とを結ぶラインである。また、前記サンプリングラインSL2,SL4,SL6,SL8は、前記所定の領域の中心と、この所定の領域rの頂点とを結ぶラインである。   In this example, eight sampling lines SL1, SL2, SL3, SL4, SL5, SL6, SL7, and SL8 are set as the sampling line SL radially from the center of the predetermined region r. The sampling lines SL1, SL3, SL5, and SL7 are lines connecting the center of the predetermined region r and the midpoint of each side of the predetermined region r. The sampling lines SL2, SL4, SL6, and SL8 are lines that connect the center of the predetermined area and the apex of the predetermined area r.

次に、ステップS13では、前記輪郭抽出部64は、前記各サンプリングラインSL1〜SL8上において、隣り合う画素間で歪みStの変化率が最大である部分を検出する。ここで、腫瘍の領域Cの輪郭の部分においては、隣り合う画素間の歪みStの変化率が大きくなるため、図9に示すように、前記各サンプリングラインSL1〜SL8上における腫瘍の領域Cの輪郭の部分に、隣り合う画素間で歪みStの変化率が最大である部分M1、M2,M3,M4,M5、M6,M7,M8が検出される。ちなみに、前記部分M1〜M8は、隣り合う画素で構成される。   Next, in step S13, the contour extraction unit 64 detects a portion on the sampling lines SL1 to SL8 where the rate of change of the distortion St is maximum between adjacent pixels. Here, in the contour portion of the tumor region C, the rate of change of the distortion St between the adjacent pixels becomes large. Therefore, as shown in FIG. 9, the tumor region C on each of the sampling lines SL1 to SL8. In the contour portion, portions M1, M2, M3, M4, M5, M6, M7, and M8 where the change rate of the distortion St between adjacent pixels is the maximum are detected. Incidentally, the portions M1 to M8 are composed of adjacent pixels.

具体的に、サンプリングラインSL1上において、図10に示すように画素p1,p2,p3,p4,・・・,p(m−1),pm,・・・,p(n−1),pnが並んでいる場合を挙げて説明する。前記画素p1,pnは前記サンプリングラインSLの両端に相当する画素であり、画素p1が所定の領域rの中心側の画素である。   Specifically, the pixels p1, p2, p3, p4,..., P (m−1), pm,..., P (n−1), pn on the sampling line SL1 as shown in FIG. A description will be given of the case where the symbols are lined up. The pixels p1 and pn are pixels corresponding to both ends of the sampling line SL, and the pixel p1 is a pixel on the center side of the predetermined region r.

前記輪郭抽出部64は、隣り合う画素p1,p2間、画素p2,p3間、画素p3,p4間、・・・、画素p(m−1),pm間、・・・、画素p(n−1),pn間の歪みの変化率を算出する。歪みの変化率とは、隣り合う画素のそれぞれについて算出された歪み値の差である。そして、前記輪郭抽出部64は、歪みStの変化率が最大になっている部分を検出する。ここでは、例えば、画素p(m−1),pm間の歪みの変化率が最大であり、歪みの変化率が最大である部分は画素p(m−1),pmであったとする。この場合、これら画素p(m−1),pmは、前記部分M1を構成する画素となる。   The contour extraction unit 64 is arranged between adjacent pixels p1 and p2, between pixels p2 and p3, between pixels p3 and p4,..., Between pixels p (m−1) and pm,. -1) The rate of change in strain between pn is calculated. The distortion change rate is a difference between distortion values calculated for adjacent pixels. Then, the contour extraction unit 64 detects a portion where the rate of change of the strain St is maximum. Here, for example, it is assumed that the distortion change rate between the pixels p (m−1) and pm is the maximum, and the portion where the distortion change rate is the maximum is the pixels p (m−1) and pm. In this case, these pixels p (m−1) and pm are pixels constituting the portion M1.

ちなみに、斜め方向のサンプリングラインS2,S4,S6,S8上における隣り合う画素について図11に基づいて説明する。図11には、サンプリングラインS2上における画素p1′,p2′,p3′,・・・,p(n−1)′,pn′が示されている。サンプリングラインS2においては、このサンプリングラインS2上における隣り合う画素p1′,p2′間、画素p2′,p3′間、・・・、画素p(n−1)′,pn′間の歪みの変化率を算出して、歪みの変化率が最大である部分の検出を行なう。   Incidentally, adjacent pixels on the diagonal sampling lines S2, S4, S6, and S8 will be described with reference to FIG. FIG. 11 shows pixels p1 ′, p2 ′, p3 ′,..., P (n−1) ′, pn ′ on the sampling line S2. In the sampling line S2, a change in distortion between the adjacent pixels p1 'and p2', between the pixels p2 'and p3' on the sampling line S2, ..., and between the pixels p (n-1) 'and pn'. The rate is calculated, and the portion where the distortion change rate is maximum is detected.

次に、ステップS14では、前記輪郭抽出部64は、各サンプリングラインSL1〜SL8において歪みStの変化率が最大である部分M1〜M8を構成する隣り合う画素のうち、前記所定の領域rの中心側の画素の歪みStについて、前記部分M1〜M8における平均値を算出し、この平均値を前記閾値StTHとする。例えば、前記部分M1を構成する画素p(m−1),pmのうち、所定の領域rの中心側の画素は画素p(m−1)である。従って、前記部分M1については、画素p(m−1)の歪みStの値を平均値算出の際に用いる。 Next, in step S14, the contour extraction unit 64 determines the center of the predetermined region r among the adjacent pixels constituting the portions M1 to M8 where the change rate of the distortion St is maximum in each sampling line SL1 to SL8. For the pixel side distortion St, an average value in the portions M1 to M8 is calculated, and this average value is set as the threshold value St TH . For example, among the pixels p (m−1) and pm constituting the portion M1, the pixel on the center side of the predetermined region r is the pixel p (m−1). Therefore, for the portion M1, the value of the distortion St of the pixel p (m−1) is used when calculating the average value.

ちなみに、前記部分M1〜M8を構成する隣り合う画素のうち、所定の領域rの中心側の画素は、腫瘍の領域Cの範囲内である可能性が高い画素であり、腫瘍の領域Cの歪みにより近い弾性を有する画素であり、またより硬い生体組織に対応する画素である。   Incidentally, among the adjacent pixels constituting the portions M1 to M8, the pixel on the center side of the predetermined region r is a pixel that is highly likely to be within the range of the tumor region C, and is a distortion of the tumor region C. It is a pixel that has an elasticity closer to that and a pixel that corresponds to a harder biological tissue.

ただし、前記輪郭抽出部64は、前記閾値StTHの算出において、前記部分M1〜M8を構成する隣り合う画素の歪み値Stを比較して、歪みStの値が高い(生体組織が硬い)方の画素の歪みStの平均値を算出するようにしてもよい。 However, in the calculation of the threshold value St TH , the contour extraction unit 64 compares the distortion values St of adjacent pixels constituting the portions M1 to M8, and the value of the distortion St is high (the living tissue is hard). The average value of the pixel distortion St may be calculated.

ステップS1において歪みの閾値StTHが算出されると、ステップS2では、腫瘍の領域Cとそれ以外の領域の境界を抽出して、腫瘍の領域Cの輪郭を抽出する。具体的には、前記所定の領域r内における各画素の歪みStと前記閾値StTHとを比較し、閾値StTH以上である画素と閾値StTH未満の画素とを特定する。そして、閾値StTH以上である画素と閾値StTH未満の画素とが互いに隣り合う部分において、閾値StTH未満の画素の部分を腫瘍の領域Cの輪郭とする。例えば、図12に示すように、画素p15,p18が閾値StTH以上の画素であり、一方で画素p11、p12、p13、p14、p16、p17、p19が閾値未満の画素であるとする。この場合、画素p11、p12、p13、p14、p16、p17、p19を腫瘍の領域Cの輪郭とする。 When the distortion threshold value St TH is calculated in step S1, the boundary between the tumor region C and the other region is extracted in step S2, and the contour of the tumor region C is extracted. Specifically, the comparing in a predetermined region r and strain St of each pixel and the threshold value St TH, identifies the pixel below the pixel and the threshold value St TH is equal to or greater than the threshold St TH. Then, in a portion where pixels that are equal to or greater than the threshold St TH and pixels that are less than the threshold St TH are adjacent to each other, the portion of the pixels that is less than the threshold St TH is defined as the contour of the tumor region C. For example, as shown in FIG. 12, it is assumed that the pixels p15 and p18 are pixels that are equal to or greater than the threshold St TH , while the pixels p11, p12, p13, p14, p16, p17, and p19 are pixels that are less than the threshold. In this case, the pixels p11, p12, p13, p14, p16, p17, and p19 are defined as the contour of the tumor region C.

次に、ステップS3では、図13に示すように、腫瘍の領域Cの輪郭に輪郭線Oを表示する。   Next, in step S3, as shown in FIG. 13, an outline O is displayed on the outline of the tumor region C.

以上説明した本例の超音波診断装置1によれば、隣り合う画素における歪みStの変化率に基づいて前記閾値StTHが設定される。従って、生体組織の変形量が異なっていても、前記腫瘍の領域Cを正確に抽出することができる。 According to the ultrasonic diagnostic apparatus 1 of the present example described above, the threshold value St TH is set based on the rate of change of the distortion St in adjacent pixels. Therefore, even if the amount of deformation of the living tissue is different, the region C of the tumor can be accurately extracted.

また、前記腫瘍の領域Cの輪郭部分に前記輪郭線Oが表示されるので、腫瘍の形状を正確に、なおかつ容易に把握することができる。   Further, since the contour line O is displayed at the contour portion of the tumor region C, the shape of the tumor can be accurately and easily grasped.

次に、上記実施形態の変形例について説明する。先ず第一変形例について説明する。上記実施形態では、二次元の弾性画像EGについて説明したが、三次元の弾性画像についても同様に適用することができる。具体的には、前記超音波プローブ2により三次元領域について超音波の送受信を行なって得られたボリュームデータに基づいて、図14に示すように互いに直交する断面XY、断面YZ、断面ZXの三断面について、図15に示すように二次元の超音波画像G1,G2,G3を前記表示部7に表示させるとともに、三次元弾性画像EG4を表示させる。前記超音波画像G1は、前記断面XYについての画像であり、Bモード画像BG1と弾性画像EG1とが合成された画像であり、前記超音波画像G2は、前記断面YZについての画像であり、Bモード画像BG2と弾性画像EG2とが合成された画像である。また、前記超音波画像G3は、前記断面ZXについての画像であり、Bモード画像BG3と弾性画像EG3とが合成された画像である。さらに、前記三次元弾性画像EG4は、例えば生体組織の硬さが増すほど不透明度が大きくなるようなテーブルを用いたボリュームレンダリングによって作成された弾性画像である(例えば特開2008−259605号公報)。ただし、前記三次元弾性画像EG4は、図15では簡略化して立方体で示されている。   Next, a modification of the above embodiment will be described. First, the first modification will be described. In the above-described embodiment, the two-dimensional elastic image EG has been described, but the same can be applied to a three-dimensional elastic image. Specifically, based on volume data obtained by transmitting / receiving ultrasonic waves in a three-dimensional region by the ultrasonic probe 2, as shown in FIG. 14, three cross sections XY, YZ, and ZX that are orthogonal to each other are shown. Regarding the cross section, as shown in FIG. 15, two-dimensional ultrasonic images G1, G2, and G3 are displayed on the display unit 7 and a three-dimensional elastic image EG4 is displayed. The ultrasonic image G1 is an image of the cross section XY, is an image obtained by combining a B-mode image BG1 and an elastic image EG1, and the ultrasonic image G2 is an image of the cross section YZ, and B This is a composite image of the mode image BG2 and the elastic image EG2. The ultrasonic image G3 is an image of the cross section ZX, and is an image obtained by combining the B-mode image BG3 and the elastic image EG3. Further, the three-dimensional elastic image EG4 is an elastic image created by volume rendering using a table whose opacity increases as the hardness of the living tissue increases (for example, Japanese Patent Application Laid-Open No. 2008-259605). . However, the three-dimensional elasticity image EG4 is shown as a cube in a simplified manner in FIG.

ステップS11における所定の領域rの設定にあっては、前記表示部7に表示された前記断面XY、前記断面YZ、断面ZXのうち、少なくとも二断面において二次元の領域r′を指定する。この二次元の領域r′は、腫瘍の領域Cを含むように設定される。これにより、三次元の所定の領域rが設定される。具体的に説明すると、例えば、図16に示すように、超音波画像G1及び超音波画像G2に二次元の領域r1′,r2′を設定する。これにより、図17に示すように、前記断面XYについて前記領域r1′が設定され、また前記断面YZについて前記領域r2′が設定されたことになる。このように領域r1′,r2′が設定されると、図18(A)に示すように、二次元の前記領域r1′を断面とし、z方向(前記断面YZの面方向)を奥行とする四角柱の領域r1′′と、図18(B)に示すように、二次元の領域r2′を断面とし、x方向(前記断面XYの面方向)を奥行とする四角柱の領域r2′′とが重なり合う領域が、前記三次元の所定の領域rになる(特に図示せず)。   In setting the predetermined region r in step S11, a two-dimensional region r ′ is specified in at least two of the cross sections XY, YZ, and ZX displayed on the display unit 7. The two-dimensional region r ′ is set so as to include the tumor region C. Thereby, a three-dimensional predetermined region r is set. More specifically, for example, as shown in FIG. 16, two-dimensional regions r1 ′ and r2 ′ are set in the ultrasonic image G1 and the ultrasonic image G2. As a result, as shown in FIG. 17, the region r1 ′ is set for the cross section XY, and the region r2 ′ is set for the cross section YZ. When the regions r1 ′ and r2 ′ are thus set, as shown in FIG. 18A, the two-dimensional region r1 ′ is a cross section, and the z direction (the surface direction of the cross section YZ) is the depth. As shown in FIG. 18 (B), a quadrangular prism region r2 ″ having a cross-section as a cross section and a depth in the x direction (plane direction of the cross section XY) as shown in FIG. The region where and overlap is the predetermined three-dimensional region r (not shown).

このようにして三次元の所定の領域rが設定されると、前記輪郭抽出部64は、この三次元の所定の領域r内にサンプリングラインSLを設定する。そして、上記と同様にして前記サンプリングラインSL上において隣り合う画素で歪みStの変化率が最大である部分を検出して歪みの閾値StTHを算出した後、腫瘍の領域Cの輪郭を抽出する。腫瘍の領域Cの輪郭が抽出されると、特に図示しないが前記三次元弾性画像EG4に、ある方向から見た場合における腫瘍の領域Cの輪郭線Oを表示する。 When the three-dimensional predetermined region r is set in this way, the contour extracting unit 64 sets the sampling line SL in the three-dimensional predetermined region r. Then, in the same manner as described above, after detecting a portion where the rate of change of the distortion St is maximum in adjacent pixels on the sampling line SL and calculating the distortion threshold St TH , the contour of the tumor region C is extracted. . When the contour of the tumor region C is extracted, the contour line O of the tumor region C when viewed from a certain direction is displayed on the three-dimensional elastic image EG4, although not particularly illustrated.

次に第二変形例について説明する。本例では、前記腫瘍の領域Cの面積や体積を算出してもよい。面積を算出する場合、前記表示制御部6は、図19に示すように、前記Bモード画像データ作成部61、前記弾性画像データ作成部62、前記合成部63、前記輪郭抽出部64、前記輪郭抽出用領域設定部65、前記輪郭線設定部66のほかに、面積算出部67を有する。また、体積を算出する場合には、前記表示制御部6は、特に図示しないが、前記面積算出部67の代わりに体積算出部を有する。前記面積算出部67及び前記体積算出部は、前記輪郭抽出部64で抽出された腫瘍の領域Cの輪郭で囲まれた領域の面積及び体積を算出する。算出された面積及び体積は、前記表示部7に表示するようにしてもよい。   Next, a second modification will be described. In this example, the area and volume of the tumor region C may be calculated. When calculating the area, the display control unit 6, as shown in FIG. 19, displays the B-mode image data creation unit 61, the elastic image data creation unit 62, the synthesis unit 63, the contour extraction unit 64, and the contour. In addition to the extraction region setting unit 65 and the contour line setting unit 66, an area calculating unit 67 is provided. When calculating the volume, the display control unit 6 includes a volume calculation unit instead of the area calculation unit 67, although not particularly illustrated. The area calculation unit 67 and the volume calculation unit calculate the area and volume of the region surrounded by the contour of the tumor region C extracted by the contour extraction unit 64. The calculated area and volume may be displayed on the display unit 7.

以上、本発明を前記各実施形態によって説明したが、本発明はその主旨を変更しない範囲で種々変更実施可能なことはもちろんである。例えば、前記輪郭抽出部64は、物理量データ(歪みSt)の代わりに、この物理量データを階調化して得られたカラー弾性画像データを用いて前記閾値StTHを算出してもよい。すなわち、前記輪郭抽出部64は、サンプリングラインSLにおいて隣り合う画素におけるカラー弾性画像データの変化率に基づいて上述と同様にして前記閾値StTHを算出してもよい。 As mentioned above, although this invention was demonstrated by each said embodiment, of course, this invention can be variously implemented in the range which does not change the main point. For example, the contour extraction unit 64 may calculate the threshold value St TH using color elastic image data obtained by gradationing the physical quantity data instead of the physical quantity data (distortion St). That is, the contour extraction unit 64 may calculate the threshold value St TH in the same manner as described above based on the rate of change of the color elastic image data in the adjacent pixels in the sampling line SL.

また、前記物理量データ作成部5は、生体組織の弾性に関する物理量として、歪みの代わりに生体組織の変形による変位や弾性率などを算出してもよい。   In addition, the physical quantity data creation unit 5 may calculate a displacement due to deformation of the living tissue, an elastic modulus, etc. instead of strain as a physical quantity related to the elasticity of the living tissue.

さらに、前記輪郭線Oの代わりに、腫瘍の領域Cを示す表示として、例えば腫瘍の領域Cに斜線を表示するなどしてもよい。   Further, instead of the outline O, as a display indicating the tumor region C, for example, a diagonal line may be displayed in the tumor region C.

1 超音波診断装置
5 物理量データ作成部(弾性データ作成部)
62 弾性画像データ作成部(弾性データ作成部)
64 輪郭抽出部
66 輪郭線設定部(表示設定部)
67 面積算出部(算出部)
7 表示部
O 輪郭線
SL,SL1〜SL8 サンプリングライン
M1〜M8 変化率が最大の部分
1 Ultrasonic diagnostic device 5 Physical quantity data creation unit (elastic data creation unit)
62 Elastic image data creation unit (elastic data creation unit)
64 Contour extraction unit 66 Contour line setting unit (display setting unit)
67 Area calculation part (calculation part)
7 Display part O Contour line SL, SL1-SL8 Sampling line M1-M8 The part with the highest rate of change

Claims (10)

生体組織に対する超音波の送受信によって得られた同一音線上の時間的に異なる二つのエコー信号に基づいて、生体組織の弾性データを作成する弾性データ作成部と、
該弾性データ作成部によって作成された弾性データに基づく弾性画像を表示する表示部と、
隣り合う画素における弾性データの変化率に基づいて弾性データの閾値を算出し、該閾値と各画素の弾性データとを比較して、生体組織における異なる弾性を有する領域の輪郭を抽出する輪郭抽出部と、
を備えることを特徴とする超音波診断装置。
Based on two temporally different echo signals on the same sound line obtained by transmission / reception of ultrasonic waves to / from a living tissue, an elasticity data creating unit that creates elasticity data of the living tissue;
A display unit for displaying an elasticity image based on the elasticity data created by the elasticity data creation unit;
A contour extraction unit that calculates a threshold value of elasticity data based on a rate of change of elasticity data in adjacent pixels, compares the threshold value with the elasticity data of each pixel, and extracts a contour of a region having different elasticity in a living tissue. When,
An ultrasonic diagnostic apparatus comprising:
前記輪郭抽出部は、前記弾性画像に複数のサンプリングラインを設定し、各サンプリングライン上において、隣り合う画素間における弾性データの変化率が最大の部分を検出し、該最大の部分における隣り合う画素のうちいずれか一方の画素についての弾性データに基づいて前記閾値を算出することを特徴とする請求項1に記載の超音波診断装置。   The contour extraction unit sets a plurality of sampling lines in the elastic image, detects a portion having the maximum rate of change in elasticity data between adjacent pixels on each sampling line, and detects adjacent pixels in the maximum portion. The ultrasonic diagnostic apparatus according to claim 1, wherein the threshold value is calculated based on elasticity data for any one of the pixels. 前記輪郭抽出部は、前記最大の部分における隣り合う画素のうち、注目領域の弾性により近い弾性を有する画素の弾性データについて、前記各サンプリングラインにおける平均を算出して、該平均を前記閾値とすることを特徴とする請求項2に記載の超音波診断装置。   The contour extraction unit calculates an average of the sampling lines for the elasticity data of pixels having elasticity closer to that of the region of interest among adjacent pixels in the maximum portion, and uses the average as the threshold value. The ultrasonic diagnostic apparatus according to claim 2. 前記注目領域の弾性により近い弾性を有する画素は、前記最大の部分における隣り合う画素のうち、より硬い生体組織に対応する画素であることを特徴とする請求項3に記載の超音波診断装置。   The ultrasonic diagnostic apparatus according to claim 3, wherein the pixel having elasticity closer to the elasticity of the region of interest is a pixel corresponding to a harder biological tissue among adjacent pixels in the maximum portion. 前記輪郭抽出部は、前記弾性画像に設定された所定の領域内に前記サンプリングラインを設定することを特徴とする請求項2〜4のいずれか一項に記載の超音波診断装置。   The ultrasonic diagnostic apparatus according to claim 2, wherein the contour extraction unit sets the sampling line in a predetermined region set in the elasticity image. 前記弾性データは、前記二つのエコー信号に基づいて算出された生体組織の弾性に関する物理量データであることを特徴とする請求項1〜5のいずれか一項に記載の超音波診断装置。   The ultrasonic diagnostic apparatus according to claim 1, wherein the elasticity data is physical quantity data relating to elasticity of a biological tissue calculated based on the two echo signals. 前記弾性データは、前記二つのエコー信号に基づいて算出された生体組織の弾性に関する物理量データを階調化して得られた弾性画像データであることを特徴とする請求項1〜5のいずれか一項に記載の超音波診断装置。   6. The elasticity image data according to claim 1, wherein the elasticity data is elasticity image data obtained by gradationizing physical quantity data relating to elasticity of a living tissue calculated based on the two echo signals. The ultrasonic diagnostic apparatus according to item. 前記輪郭抽出部によって抽出された輪郭で囲まれる領域を示す表示を前記弾性画像上に表示させる表示設定部を備えることを特徴とする請求項1〜7のいずれか一項に記載の超音波診断装置。   The ultrasonic diagnosis according to claim 1, further comprising: a display setting unit configured to display a display indicating an area surrounded by the contour extracted by the contour extraction unit on the elastic image. apparatus. 前記輪郭抽出部によって抽出された輪郭で囲まれる領域の面積又は体積を算出する算出部を備えることを特徴とする請求項1〜8のいずれか一項に記載の超音波診断装置。   The ultrasonic diagnostic apparatus according to claim 1, further comprising a calculation unit that calculates an area or volume of a region surrounded by the contour extracted by the contour extraction unit. コンピュータに、
生体組織に対する超音波の送受信によって得られた同一音線上の時間的に異なる二つのエコー信号に基づいて、生体組織の弾性データを作成する弾性データ作成機能と、
隣り合う画素における弾性データの変化率に基づいて弾性データの閾値を算出し、該閾値と各画素の弾性データとを比較して、生体組織における異なる弾性を有する領域の輪郭を抽出する輪郭抽出機能と、
を実行させることを特徴とする超音波診断装置の制御プログラム。
On the computer,
Based on two echo signals that differ in time on the same sound line obtained by transmitting and receiving ultrasonic waves to and from the living tissue, an elasticity data creation function that creates elasticity data of the living tissue,
A contour extraction function that calculates a threshold value of elasticity data based on a rate of change of elasticity data in adjacent pixels, compares the threshold value with the elasticity data of each pixel, and extracts a contour of a region having different elasticity in a living tissue. When,
A control program for an ultrasonic diagnostic apparatus, characterized in that
JP2010067702A 2010-03-24 2010-03-24 Ultrasonic diagnostic apparatus and control program therefor Active JP5624345B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2010067702A JP5624345B2 (en) 2010-03-24 2010-03-24 Ultrasonic diagnostic apparatus and control program therefor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010067702A JP5624345B2 (en) 2010-03-24 2010-03-24 Ultrasonic diagnostic apparatus and control program therefor

Publications (3)

Publication Number Publication Date
JP2011200282A true JP2011200282A (en) 2011-10-13
JP2011200282A5 JP2011200282A5 (en) 2012-07-19
JP5624345B2 JP5624345B2 (en) 2014-11-12

Family

ID=44877717

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010067702A Active JP5624345B2 (en) 2010-03-24 2010-03-24 Ultrasonic diagnostic apparatus and control program therefor

Country Status (1)

Country Link
JP (1) JP5624345B2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015119819A (en) * 2013-12-24 2015-07-02 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ultrasound diagnostic device
JP2016047474A (en) * 2016-01-08 2016-04-07 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ultrasonic diagnostic device
JP7299100B2 (en) 2019-08-08 2023-06-27 富士フイルムヘルスケア株式会社 ULTRASOUND DIAGNOSTIC DEVICE AND ULTRASOUND IMAGE PROCESSING METHOD

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004135929A (en) * 2002-10-18 2004-05-13 Hitachi Medical Corp Ultrasonic diagnostic apparatus
WO2005025425A1 (en) * 2003-09-12 2005-03-24 Hitachi Medical Corporation Ultrasonograph
JP2007275456A (en) * 2006-04-11 2007-10-25 Aloka Co Ltd Ultrasonic diagnosing system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004135929A (en) * 2002-10-18 2004-05-13 Hitachi Medical Corp Ultrasonic diagnostic apparatus
WO2005025425A1 (en) * 2003-09-12 2005-03-24 Hitachi Medical Corporation Ultrasonograph
JP2007275456A (en) * 2006-04-11 2007-10-25 Aloka Co Ltd Ultrasonic diagnosing system

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015119819A (en) * 2013-12-24 2015-07-02 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ultrasound diagnostic device
JP2016047474A (en) * 2016-01-08 2016-04-07 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ultrasonic diagnostic device
JP7299100B2 (en) 2019-08-08 2023-06-27 富士フイルムヘルスケア株式会社 ULTRASOUND DIAGNOSTIC DEVICE AND ULTRASOUND IMAGE PROCESSING METHOD

Also Published As

Publication number Publication date
JP5624345B2 (en) 2014-11-12

Similar Documents

Publication Publication Date Title
US20110194748A1 (en) Ultrasonic diagnostic apparatus and ultrasonic image display method
JP5038304B2 (en) Ultrasonic diagnostic equipment
JP5436533B2 (en) Ultrasonic diagnostic apparatus and elastic image display method
JP5087341B2 (en) Ultrasonic diagnostic equipment
US20110306884A1 (en) Ultrasonic diagnostic apparatus
CN105122299A (en) Ultrasound image displaying apparatus and method for displaying ultrasound image
JP6290336B2 (en) Ultrasonic diagnostic equipment
JP5074097B2 (en) Ultrasonic diagnostic equipment
KR101629541B1 (en) Ultrasonic diagnostic apparatus and control program thereof
JP2010124946A (en) Ultrasonic diagnosing apparatus, and program
JP5624345B2 (en) Ultrasonic diagnostic apparatus and control program therefor
JP2011120614A (en) Ultrasonic diagnostic apparatus and ultrasonic image display method
JP5345477B2 (en) Ultrasonic diagnostic apparatus and control program therefor
JP5710566B2 (en) Ultrasonic diagnostic apparatus and control program therefor
JP5334103B2 (en) Ultrasonic diagnostic equipment
JP2013240721A (en) Ultrasonic diagnostic apparatus
JP5606998B2 (en) Ultrasonic diagnostic apparatus and control program therefor
JP5394693B2 (en) Ultrasonic diagnostic equipment
JP5725732B2 (en) Ultrasonic diagnostic apparatus and control program therefor
JP5879230B2 (en) Ultrasonic diagnostic apparatus and control program therefor
KR101574821B1 (en) Ultrasonic diagnosis apparatus and program for controlling the same
JP5626986B2 (en) Ultrasonic diagnostic apparatus and control program therefor
JP5922521B2 (en) Ultrasonic diagnostic apparatus and control program therefor
JP2017042179A (en) Ultrasonic diagnostic device
JP2019050961A (en) Ultrasonic diagnostic apparatus

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120605

A625 Written request for application examination (by other person)

Free format text: JAPANESE INTERMEDIATE CODE: A625

Effective date: 20121016

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20131122

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20131125

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140122

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140926

R150 Certificate of patent or registration of utility model

Ref document number: 5624345

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

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