JP2013183982A - Ultrasonic diagnostic apparatus and elastic image generation method - Google Patents

Ultrasonic diagnostic apparatus and elastic image generation method Download PDF

Info

Publication number
JP2013183982A
JP2013183982A JP2012053110A JP2012053110A JP2013183982A JP 2013183982 A JP2013183982 A JP 2013183982A JP 2012053110 A JP2012053110 A JP 2012053110A JP 2012053110 A JP2012053110 A JP 2012053110A JP 2013183982 A JP2013183982 A JP 2013183982A
Authority
JP
Japan
Prior art keywords
image
measurement
elasticity
current time
compression condition
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2012053110A
Other languages
Japanese (ja)
Inventor
Takeshi Matsumura
剛 松村
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.)
Hitachi Ltd
Original Assignee
Hitachi Aloka Medical Ltd
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 Hitachi Aloka Medical Ltd filed Critical Hitachi Aloka Medical Ltd
Priority to JP2012053110A priority Critical patent/JP2013183982A/en
Publication of JP2013183982A publication Critical patent/JP2013183982A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide an elastic image capable of displaying a tissue elasticity response acquired on the basis of a pressure condition of the same size.SOLUTION: An ultrasonic diagnostic apparatus includes: an elasticity information calculation part 13 for obtaining elasticity information of an organism tissue in multiple measurement points within an image area and a measurement pressure condition of the multiple measurement points, on the basis of a reflection echo signal to be acquired by transmitting an ultrasonic wave from a probe; an image display unit 10 for displaying elastic images generated on the basis of elastic information; and a pressure condition evaluation part 18 for selecting elastic information for the generation of elastic images to be displayed in the image display unit in accordance with a measurement pressure condition at a current time or before the current time, and allowing the image display unit to display elastic images at the multiple measurement points generated on the basis of the selected elastic information at a current time and before the current time.

Description

本発明は、超音波診断装置及び弾性画像生成方法に係り、特に、被検体の生体組織に圧迫を加えて計測される組織弾性を画像化して診断に供する超音波診断装置及び弾性画像生成方法に関する。   The present invention relates to an ultrasonic diagnostic apparatus and an elastic image generation method, and more particularly, to an ultrasonic diagnostic apparatus and an elastic image generation method for imaging a tissue elasticity measured by applying pressure to a living tissue of a subject and using it for diagnosis. .

超音波診断装置は、被検体の表面に超音波探触子を当て、その探触子から被検体に超音波を送信し、被検体内部からの超音波の反射波を受信し、その受信信号である反射エコー信号に基づいて被検体の各部の生体情報を断層像などの画像により表示して診断に供するものである。   The ultrasonic diagnostic apparatus applies an ultrasonic probe to the surface of the subject, transmits ultrasonic waves from the probe to the subject, receives ultrasonic reflected waves from the inside of the subject, and receives the received signal. Based on the reflected echo signal, the biological information of each part of the subject is displayed as an image such as a tomographic image for diagnosis.

一方、特許文献1,2等に記載されているように、被検体に圧迫力を加えて内部に生じた組織の歪みや弾性率等の弾性に関する情報(以下、弾性情報という。)を求めて弾性画像を生成し、その弾性画像を表示して検者の病変部等の診断に供することが行われている。例えば、検者は、圧迫力を加減しながら計測される弾性画像を観察することにより、正常組織、ガン細胞、腫瘍などの病変を鑑別する。   On the other hand, as described in Patent Documents 1 and 2, etc., information relating to elasticity (hereinafter referred to as elasticity information) such as tissue strain and elastic modulus generated inside by applying a compression force to the subject is obtained. An elastic image is generated, and the elastic image is displayed to be used for diagnosis of a lesioned part of an examiner. For example, the examiner distinguishes lesions such as normal tissues, cancer cells, and tumors by observing elastic images measured while adjusting the compression force.

ところで、生体組織は、圧迫に対して非線形の弾性応答を示し、圧縮量などの圧迫の程度(以下、圧迫条件という。)に応じて硬さが変化することが知られている。例えば、非特許文献1によれば、図8に示すような乳腺組織の乳腺部位と脂肪部位を圧縮した際の応力(kPa)と歪み(%)の関係の実測値が報告されている。同図によれば、脂肪部位の応力−歪みは直線的な関係であるが、乳腺部位の応力−歪みの関係は非線形である。したがって、図8の応力−歪み曲線の傾きで与えられる弾性率(ヤング率)は、図9に模式的に示すように、脂肪部位では一定値となるが、乳腺部位では加えられた歪みに応じて変化することが分かる。   By the way, it is known that a living tissue exhibits a non-linear elastic response to compression, and the hardness changes according to the degree of compression such as the amount of compression (hereinafter referred to as compression condition). For example, Non-Patent Document 1 reports an actual measurement value of the relationship between stress (kPa) and strain (%) when a mammary gland region and a fat region of a mammary gland tissue are compressed as shown in FIG. According to the figure, the stress-strain in the fat region is a linear relationship, but the stress-strain relationship in the mammary gland region is non-linear. Therefore, the elastic modulus (Young's modulus) given by the slope of the stress-strain curve in FIG. 8 is a constant value in the fat region, as shown schematically in FIG. 9, but in accordance with the applied strain in the mammary region. It can be seen that changes.

一方、非特許文献2には、歪みと弾性率の関係のデータを取得した後、歪み−弾性率曲線を関数で近似し、最小二乗法などにより最もよく近似する曲線からその非線形性を表す非線形パラメータを抽出し、その組織の非線形性を表す情報として評価する試みが提案されている。   On the other hand, in Non-Patent Document 2, after acquiring data on the relationship between strain and elastic modulus, a strain-elastic modulus curve is approximated by a function, and the nonlinearity is expressed from a curve that is best approximated by a least square method or the like. Attempts have been made to extract parameters and evaluate them as information representing the nonlinearity of the tissue.

特開平5−317313号公報JP-A-5-317313 特開2000−60853号公報JP 2000-60853 A

krouskop T, et al : Elastic Moduli of Breast and Prostate TissuesUnder Compression. Ultrasonic Imaging 20 : 260-274, 1998.krouskop T, et al: Elastic Moduli of Breast and Prostate TissuesUnder Compression. Ultrasonic Imaging 20: 260-274, 1998. 「超音波による組織の非線形弾性特性の画像化」新田他、電子情報通信学会論文誌 2001/12 Vol. J84-A No.1“Imaging Nonlinear Elastic Properties of Tissues by Ultrasound” Nitta et al., IEICE Transactions 2001/12 Vol. J84-A No.1

このように、組織の硬さの情報を用いて組織の良悪性を鑑別するにあたって、例えば、弾性率がどのような圧迫条件の下で計測されたかが非常に重要な情報であり、圧縮量などの圧迫条件が大きく異なる条件で計測された診断対象については、良悪性を誤って診断してしまうことも考えられる。   Thus, in distinguishing benign and malignant tissue using tissue hardness information, for example, under what compression conditions the elastic modulus was measured is very important information, such as compression amount It may be considered that benign and malignant are erroneously diagnosed with respect to a diagnosis target measured under conditions in which the compression conditions are greatly different.

しかしながら、特許文献1,2等に記載された従来の技術では、圧迫の程度に応じて硬さが変化する生体組織の非線形性を考慮していない。
つまり、生体組織の弾性には非線形性があるため、同じ大きさの圧迫条件下で得られた弾性情報をもってして比較しない限り、正しい評価を行うことができないと考えられる。
However, the conventional techniques described in Patent Documents 1 and 2 do not consider the non-linearity of the biological tissue whose hardness changes according to the degree of compression.
In other words, since the elasticity of the living tissue has non-linearity, it is considered that a correct evaluation cannot be performed unless the comparison is performed using the elasticity information obtained under the compression condition of the same size.

本発明は、超音波診断装置において、同一の大きさの圧迫条件に基づいて得られた組織弾性の応答を表示することを可能とする弾性画像を提供することを課題とする。   It is an object of the present invention to provide an elastic image that allows a response of tissue elasticity obtained based on a compression condition of the same size to be displayed in an ultrasonic diagnostic apparatus.

上記の課題を解決するため、本発明の超音波診断装置は、探触子から超音波を送信して得られる反射エコー信号に基づいて、画像領域内の複数の計測点における生体組織の弾性情報と前記複数の計測点の計測圧迫条件を求める弾性情報演算部と、前記弾性情報に基づいて生成された弾性画像を表示する画像表示器と、前記画像表示器に表示させる前記弾性画像を生成するための前記弾性情報を、現時刻もしくは現時刻前の前記計測圧迫条件に応じて選択し、選択された現時刻及び現時刻前の前記弾性情報に基づいて生成された前記複数の計測点における前記弾性画像を前記画像表示器に表示させる圧迫条件評価部とを備えることを特徴とする。   In order to solve the above-described problem, the ultrasonic diagnostic apparatus of the present invention is based on a reflected echo signal obtained by transmitting an ultrasonic wave from a probe, and information on the elasticity of biological tissue at a plurality of measurement points in an image region. And an elasticity information calculation unit for obtaining measurement compression conditions for the plurality of measurement points, an image display for displaying an elasticity image generated based on the elasticity information, and generating the elasticity image to be displayed on the image display The elasticity information for the current measurement time is selected according to the measurement compression condition at the current time or before the current time, and at the plurality of measurement points generated based on the selected current time and the elasticity information before the current time. And a compression condition evaluation unit that displays an elastic image on the image display.

本発明によれば、同一の大きさの圧迫条件に基づいて得られた組織弾性の応答を表示することを可能とする弾性画像を提供することができる。   ADVANTAGE OF THE INVENTION According to this invention, the elasticity image which makes it possible to display the response | compatibility of the tissue elasticity obtained based on the compression conditions of the same magnitude | size can be provided.

本発明の弾性画像生成方法を実施する一実施形態の超音波診断装置のブロック構成図である。It is a block block diagram of the ultrasonic diagnostic apparatus of one Embodiment which implements the elasticity image generation method of this invention. 圧迫を加えて弾性画像を生成する生体組織を含む部位のBモード画像を示す図である。It is a figure which shows the B mode image of the site | part containing the biological tissue which applies compression and produces | generates an elastic image. 実施例1において生成される弾性画像の一例であり、ROI周辺のみが代表計測圧迫条件を満たして白黒輝度で階調化されている状態の画像を示す図である。It is an example of the elasticity image produced | generated in Example 1, and is a figure which shows the image of the state which only the periphery of ROI satisfy | fills the representative measurement compression conditions, and is gradation-ized by the monochrome luminance. 実施例1において生成される弾性画像の一例であり、現時刻の応力が代表計測圧迫条件を満たす領域のみならず、代表計測圧迫条件を満たしていない領域についても、ROI内に負荷する応力と同じ圧迫条件下で得られたヤング率の値がマッピングされた状態の画像を示す図である。It is an example of the elasticity image produced | generated in Example 1, and it is the same as the stress loaded in ROI not only in the area | region where the stress of the present time satisfy | fills representative measurement compression conditions, but the area | region which does not satisfy representative measurement compression conditions It is a figure which shows the image of the state in which the value of the Young's modulus obtained on compression conditions was mapped. 実施例1において生成される弾性画像を生成する際に、現時刻の応力は代表計測圧迫条件を満たさないが、現時刻前の応力が代表計測圧迫条件を満たした領域について、ROI内に負荷する応力と同じ圧迫条件下で得られたヤング率の値がマッピングされた状態のイメージを示す図である。When the elastic image generated in the first embodiment is generated, the stress at the current time does not satisfy the representative measurement compression condition, but the region where the stress before the current time satisfies the representative measurement compression condition is loaded in the ROI. It is a figure which shows the image of the state in which the value of the Young's modulus obtained on the same compression conditions as stress was mapped. (a)から(c)は、画像に対するROIの位置を固定した場合における、ROIに属する組織の圧迫操作による移動の状態、(d)から(f)は、ROIに属する組織の行き先を圧迫の程度に応じて追従するように、ROIを移動させた状態をそれぞれ説明するための図である。(a) to (c) are the movement states of the tissue belonging to the ROI when the position of the ROI is fixed with respect to the image, and (d) to (f) are the destinations of the tissue belonging to the ROI. It is a figure for demonstrating the state which moved ROI so that it may track according to a grade. 実施例2において生成される生体組織の応力と弾性率の関係を表示させるとともに、基準圧迫条件を設定するためのスライドバーを表示させたグラフ画像を示す図である。It is a figure which shows the graph image which displayed the relationship between the stress and elastic modulus of the biological tissue produced | generated in Example 2, and displayed the slide bar for setting a reference | standard compression condition. 乳腺組織の乳腺部位と脂肪部位を圧迫した際の歪みと応力の関係を示す図である。It is a figure which shows the relationship between the distortion at the time of compressing the mammary gland part and fat part of a mammary gland tissue. 乳腺組織の乳腺部位と脂肪部位を圧迫した際の歪みと弾性率(ヤング率)の関係を模式的に示す図である。It is a figure which shows typically the relationship between the distortion at the time of compressing the mammary gland part of a mammary gland tissue, and a fat part, and an elasticity modulus (Young's modulus).

以下、本発明の超音波診断装置及び弾性画像生成方法について、添付図面を参照して説明する。図1には、本発明の弾性画像生成方法を実施する一実施形態に係る超音波診断装置のブロック構成図を示す。図1に示すように、被検体1に当接して用いられる超音波の探触子2は、被検体1との間で超音波を送信及び受信する複数の振動子を有して形成されている。探触子2は、送信回路3から供給される超音波パルスにより駆動される。送受信制御回路4は、探触子2の複数の振動子を駆動する超音波パルスの送信タイミングを制御して、被検体1内に設定される焦点に向けて超音波ビームを形成するようになっている。また、送受信制御回路4は、探触子2の振動子の配列方向に電子的に超音波ビームを走査するようになっている。   Hereinafter, an ultrasonic diagnostic apparatus and an elastic image generation method according to the present invention will be described with reference to the accompanying drawings. FIG. 1 is a block diagram of an ultrasonic diagnostic apparatus according to an embodiment that implements the elastic image generation method of the present invention. As shown in FIG. 1, an ultrasonic probe 2 used in contact with a subject 1 is formed with a plurality of transducers that transmit and receive ultrasonic waves to and from the subject 1. Yes. The probe 2 is driven by ultrasonic pulses supplied from the transmission circuit 3. The transmission / reception control circuit 4 controls the transmission timing of ultrasonic pulses for driving the plurality of transducers of the probe 2 to form an ultrasonic beam toward the focal point set in the subject 1. ing. The transmission / reception control circuit 4 is configured to electronically scan the ultrasonic beam in the direction in which the transducers of the probe 2 are arranged.

一方、探触子2は、被検体1内から発生する反射エコー信号を受信して受信回路5に出力する。受信回路5は、送受信制御回路4から入力されるタイミング信号に従って、反射エコー信号を取り込んで増幅などの受信処理を行う。受信回路5により処理された反射エコー信号は、整相加算回路6において複数の振動子により受信された反射エコー信号の位相を合わせて加算することにより増幅される。整相加算回路6において整相加算された反射エコー信号は、信号処理部7に入力され、ゲイン補正、ログ圧縮、検波、輪郭強調、フィルタ処理等の信号処理がなされる。   On the other hand, the probe 2 receives a reflected echo signal generated from the subject 1 and outputs it to the receiving circuit 5. In accordance with the timing signal input from the transmission / reception control circuit 4, the reception circuit 5 takes in the reflected echo signal and performs reception processing such as amplification. The reflected echo signal processed by the receiving circuit 5 is amplified by adding together the phases of the reflected echo signals received by the plurality of transducers in the phasing and adding circuit 6. The reflected echo signal phased and added in the phase adjusting and adding circuit 6 is input to the signal processing unit 7 and subjected to signal processing such as gain correction, log compression, detection, contour enhancement, and filter processing.

信号処理部7により処理された反射エコー信号は、白黒スキャンコンバータ8に導かれて超音波ビームの走査面に対応した2次元の断層像データ(ディジタルデータ)に変換される。これらの信号処理部7と白黒スキャンコンバータ8によって断層像(Bモード像)の画像再構成手段が構成される。白黒スキャンコンバータ8から出力される断層像データは、切替加算部9を介して画像表示器10に供給され、図2に示すようなBモード像が表示されるようになっている。   The reflected echo signal processed by the signal processing unit 7 is guided to the monochrome scan converter 8 and converted into two-dimensional tomographic image data (digital data) corresponding to the scanning surface of the ultrasonic beam. These signal processing unit 7 and monochrome scan converter 8 constitute a tomographic image (B-mode image) image reconstruction means. The tomographic image data output from the black-and-white scan converter 8 is supplied to the image display 10 via the switching addition unit 9, and a B-mode image as shown in FIG. 2 is displayed.

一方、整相加算回路6から出力される反射エコー信号は、フレームデータ取得部11に導かれる。フレームデータ取得部11は、被検体1の生体組織に圧迫を加えて探触子2から超音波を送信して得られる反射エコー信号のフレームデータを取得する。つまり、フレームデータ取得部11は、超音波ビームの走査面(断層面)に対応する反射エコー信号群を、フレームデータとして複数フレーム分を取得してメモリなどに格納する。変位計測部12は、フレームデータ取得部11に格納されている取得時刻が異なる複数対のフレームデータを順次取り込み、取り込んだ一対のフレームデータに基づいて断層面における複数の計測点の変位ベクトルを求め、変位フレームデータとして弾性情報演算部13に出力するようになっている。   On the other hand, the reflected echo signal output from the phasing addition circuit 6 is guided to the frame data acquisition unit 11. The frame data acquisition unit 11 acquires frame data of a reflected echo signal obtained by applying pressure to the living tissue of the subject 1 and transmitting ultrasonic waves from the probe 2. That is, the frame data acquisition unit 11 acquires a plurality of frames of reflected echo signal groups corresponding to the scanning plane (tomographic plane) of the ultrasonic beam as frame data, and stores it in a memory or the like. The displacement measuring unit 12 sequentially captures a plurality of pairs of frame data with different acquisition times stored in the frame data acquiring unit 11, and obtains displacement vectors of a plurality of measurement points on the tomographic plane based on the captured pair of frame data. The displacement information is output to the elasticity information calculation unit 13 as displacement frame data.

弾性情報演算部13は、探触子2から超音波を送信して得られる反射エコー信号に基づいて、画像領域内の複数の計測点における生体組織の弾性情報と複数の計測点の計測圧迫条件(計測点に加えられる圧迫の程度を示す条件)を求める。より具体的には、弾性情報演算部13は、フレームデータ取得部11により取得されたフレームデータに基づいて、計測圧迫条件を求めている。本実施形態において、弾性情報演算部13は、変位フレームデータに基づいて各計測点の生体組織の歪み変化を求める歪み演算部と、歪みフレームデータに基づいて各計測点の生体組織の弾性率を求める弾性率演算部と、歪み演算部で求めた歪み変化を積算して計測圧迫条件の一例である歪みを求める積算歪み演算部と、圧力計測部17で求めた圧力に基づいて計測圧迫条件の一例である被検体内部の応力(分布)を求める応力演算部を有して構成されている。弾性情報演算部13にて求められた弾性率のフレームデータは弾性情報処理部14に出力され、歪み及び応力(分布)は圧迫条件評価部18に出力されるようになっている。
ここで、応力(分布)を求める方法は任意であり、圧力センサに限らず、音響カプラなどを用いて表面圧力を求め、有限要素法などを利用して応力を推定する方法など、多くの公知例がある。
Based on the reflected echo signal obtained by transmitting an ultrasonic wave from the probe 2, the elasticity information calculation unit 13 is based on the elasticity information of the biological tissue at a plurality of measurement points in the image region and the measurement compression conditions for the plurality of measurement points. (Condition indicating the degree of compression applied to the measurement point). More specifically, the elasticity information calculation unit 13 obtains the measurement compression condition based on the frame data acquired by the frame data acquisition unit 11. In the present embodiment, the elasticity information calculation unit 13 calculates a strain calculation unit for obtaining a strain change of the biological tissue at each measurement point based on the displacement frame data, and calculates an elastic modulus of the biological tissue at each measurement point based on the strain frame data. An elastic modulus calculation unit to be obtained, an integrated strain calculation unit to obtain a strain that is an example of a measurement compression condition by integrating the strain change obtained by the strain calculation unit, and a measurement compression condition based on the pressure obtained by the pressure measurement unit 17 As an example, the apparatus includes a stress calculation unit that obtains stress (distribution) inside the subject. The elastic modulus frame data obtained by the elastic information calculation unit 13 is output to the elastic information processing unit 14, and the strain and stress (distribution) are output to the compression condition evaluation unit 18.
Here, the method for obtaining the stress (distribution) is arbitrary, and is not limited to the pressure sensor, and many known methods such as a method for obtaining the surface pressure using an acoustic coupler or the like and estimating the stress using a finite element method or the like. There is an example.

弾性情報処理部14は、弾性情報演算部13から入力される各弾性情報のフレームデータに対して、座標平面内におけるスムージング処理、コントラスト最適化処理、フレーム間における時間軸方向のスムージング処理などの様々な画像処理を施して、カラースキャンコンバータ15に送出するようになっている。   The elasticity information processing unit 14 performs various processes such as smoothing processing in the coordinate plane, contrast optimization processing, and smoothing processing in the time axis direction between frames on the frame data of each elasticity information input from the elasticity information calculation unit 13. The image is processed and sent to the color scan converter 15.

カラースキャンコンバータ15は、弾性情報処理部14により処理された弾性率のフレームデータを取り込み、設定された弾性率のカラーマップに従って、フレームデータの画素ごとに色調コードを付与してカラー弾性画像を生成するようになっている。   The color scan converter 15 takes in the elastic modulus frame data processed by the elastic information processing unit 14 and generates a color elastic image by assigning a tone code to each pixel of the frame data in accordance with the color map of the set elastic modulus. It is supposed to be.

カラースキャンコンバータ15により生成されたカラー弾性画像は、図3に示すように、切替加算部9を介して画像表示器10に表示されるようになっている。また、切替加算部9は、白黒スキャンコンバータ8から出力される白黒の断層像と、カラースキャンコンバータ15から出力されるカラー弾性画像とを入力し、両画像を切り替えていずれか一方を表示させる機能と、両画像の一方を半透明にして加算合成して画像表示器10に重ねて表示させる機能と、両画像を並べて表示させる機能を有して形成されている。また、切替加算部9から出力される画像データは、装置制御インターフェイス部19の制御に従ってシネメモリ部20に格納されるようになっている。シネメモリ部20に格納された画像データは、装置制御インターフェイス部19の制御に従って画像表示器10に表示されるようになっている。このように、画像表示器10は、弾性情報に基づいて生成された弾性画像を被者が視認可能に表示する。   The color elasticity image generated by the color scan converter 15 is displayed on the image display 10 via the switching addition unit 9 as shown in FIG. The switching adder 9 also receives a black and white tomographic image output from the black and white scan converter 8 and a color elastic image output from the color scan converter 15 and switches both images to display one of them. In addition, one of the two images is translucent, added and combined, displayed on the image display 10 and displayed, and the two images are displayed side by side. Further, the image data output from the switching addition unit 9 is stored in the cine memory unit 20 under the control of the device control interface unit 19. The image data stored in the cine memory unit 20 is displayed on the image display 10 under the control of the device control interface unit 19. As described above, the image display 10 displays the elasticity image generated based on the elasticity information so that the subject can visually recognize the elasticity image.

圧迫条件評価部18は、画像表示器10に表示させる弾性画像を生成するための弾性情報を、現時刻もしくは現時刻前(現時刻より過去に遡った任意の時刻、以下、過去という。)の計測圧迫条件に応じて選択し、選択された現時刻及び過去の弾性情報に基づいて生成された複数の計測点における弾性画像を画像表示器10に表示させる。より具体的には、圧迫条件評価部18は、選択された現時刻及び過去の弾性情報に基づいて弾性情報処理部14に指令を送るとともに、カラースキャンコンバータ15で生成された複数の計測点におけるカラー弾性画像を画像表示器10に出力表示させるようになっている。   The compression condition evaluation unit 18 generates elasticity information for generating an elasticity image to be displayed on the image display 10 at the current time or before the current time (any time that goes back in the past from the current time, hereinafter referred to as the past). The image display 10 displays an elasticity image at a plurality of measurement points selected based on the measurement compression condition and generated based on the selected current time and past elasticity information. More specifically, the compression condition evaluation unit 18 sends a command to the elastic information processing unit 14 based on the selected current time and past elastic information, and at a plurality of measurement points generated by the color scan converter 15. A color elastic image is output and displayed on the image display 10.

このように構成される本実施形態の基本的な動作について説明する。まず、探触子2により被検体1に圧迫を加えて被検体1に超音波ビームを走査し、走査面からの反射エコー信号を連続的に受信する。そして、整相加算回路6から出力される反射エコー信号に基づいて、信号処理部7及び白黒スキャンコンバータ8により断層像が再構成され、切替加算部9を介して画像表示器10に表示される。   The basic operation of this embodiment configured as described above will be described. First, the probe 2 applies pressure to the subject 1 to scan the subject 1 with an ultrasonic beam, and continuously receives reflected echo signals from the scanning surface. Based on the reflected echo signal output from the phasing addition circuit 6, the tomographic image is reconstructed by the signal processing unit 7 and the black and white scan converter 8 and displayed on the image display 10 via the switching addition unit 9. .

一方、フレームデータ取得部11は、反射エコー信号を取り込んでフレームレートに同期させてフレームデータを繰り返し取得し、内蔵されたフレームメモリ内に時系列順に保存する。そして、取得時刻が異なる一対のフレームデータを単位として、連続的に複数対のフレームデータを選択して変位計測部12に出力する。変位計測部12は、選択された一対のフレームデータを1次元もしくは2次元相関処理し、走査面における複数の計測点の変位を計測して変位フレームデータを生成する。この変位ベクトルの検出法としては、例えば、特開平5−317313号公報等に記載されているブロックマッチング法、又はグラジェント法が知られている。ブロックマッチング法は、画像を例えばN×N画素からなるブロックに分け、現フレーム中の着目しているブロックに最も近似しているブロックを前フレームから探索し、これに基づいて計測点の変位を求める。また、一対のRF信号フレームデータの同一領域における自己相関を計算して変位を算出することができる。   On the other hand, the frame data acquisition unit 11 takes in the reflected echo signal, repeatedly acquires the frame data in synchronization with the frame rate, and stores it in the built-in frame memory in chronological order. Then, a plurality of pairs of frame data are successively selected and output to the displacement measuring unit 12 with a pair of frame data having different acquisition times as a unit. The displacement measuring unit 12 performs one-dimensional or two-dimensional correlation processing on the selected pair of frame data, measures the displacement of a plurality of measurement points on the scanning plane, and generates displacement frame data. As a method for detecting the displacement vector, for example, a block matching method or a gradient method described in JP-A-5-317313 is known. In the block matching method, an image is divided into blocks of N × N pixels, for example, a block closest to the target block in the current frame is searched from the previous frame, and based on this, the displacement of the measurement point is determined. Ask. Further, the displacement can be calculated by calculating the autocorrelation in the same region of the pair of RF signal frame data.

弾性情報演算部13は、変位フレームデータを取り込んで、各計測点の歪み変化を求め、求めた歪み変化に基づいて弾性情報である弾性率を演算し、弾性率フレームデータを弾性情報処理部14に出力する。歪み変化の演算は、公知のように変位を空間微分することによって計算される。また、求めた歪み変化に基づいて各計測点の弾性率を演算する。弾性率を求める場合は、圧力計測部17により計測された圧力の計測値を取り込み、これに基づいて各計測点における応力を演算する。圧力計測部17は、探触子2の超音波送受信面と被検体1との間に設けられた圧力センサ16により検出された圧力に基づいて、被検体1内部の計測点における応力を演算する。つまり、弾性情報演算部13の弾性率演算部は、各計測点における応力と、弾性情報演算部13で求めた歪みフレームデータから走査面上の各計測点の弾性率E(例えば、ヤング率)を演算し、弾性情報処理部14に出力する。   The elasticity information calculation unit 13 takes in the displacement frame data, obtains a strain change at each measurement point, calculates an elastic modulus that is elasticity information based on the obtained strain change, and obtains the elasticity frame data as an elastic information processing unit 14. Output to. The calculation of the strain change is calculated by spatially differentiating the displacement as is well known. Further, the elastic modulus at each measurement point is calculated based on the obtained strain change. When calculating | requiring an elasticity modulus, the measured value of the pressure measured by the pressure measurement part 17 is taken in, and the stress in each measurement point is calculated based on this. The pressure measurement unit 17 calculates the stress at the measurement point inside the subject 1 based on the pressure detected by the pressure sensor 16 provided between the ultrasonic transmission / reception surface of the probe 2 and the subject 1. . That is, the elastic modulus calculation unit of the elastic information calculation unit 13 calculates the elastic modulus E (for example, Young's modulus) of each measurement point on the scanning surface from the stress at each measurement point and the strain frame data obtained by the elastic information calculation unit 13. Is calculated and output to the elastic information processing unit 14.

弾性情報処理部14は、入力される弾性率をスムージング処理などの処理を施してカラースキャンコンバータ15に出力する。カラースキャンコンバータ15は、弾性情報に基づいてカラー弾性画像を生成する。カラー画像は、例えば、256階調化による色調のグラデーションで画素単位ごとにフレームデータの弾性率に応じて色付けされる。なお、カラースキャンコンバータ15に代えて、白黒スキャンコンバータを用いることができる。この場合は、弾性率が大きい領域は輝度を明るく、逆に小さい領域は輝度を暗くするなどにより、良性又は悪性を鑑別できるようにすることができる。   The elastic information processing unit 14 performs processing such as smoothing processing on the input elastic modulus and outputs the processed elastic modulus to the color scan converter 15. The color scan converter 15 generates a color elasticity image based on the elasticity information. For example, the color image is colored according to the elastic modulus of the frame data for each pixel unit with a gradation of color tone by 256 gradations. In place of the color scan converter 15, a black and white scan converter can be used. In this case, it is possible to discriminate benign or malignant, for example, by increasing the luminance in a region having a large elastic modulus and conversely decreasing the luminance in a small region.

以下に、本発明の超音波診断装置において弾性画像を生成する具体的な実施例について説明する。なお、各実施例は、本発明の特徴部である弾性情報演算部13、弾性情報処理部14、圧迫条件評価部18、装置制御インターフェイス部19及びシネメモリ部20等により実施される。   A specific embodiment for generating an elastic image in the ultrasonic diagnostic apparatus of the present invention will be described below. Each embodiment is implemented by the elastic information calculation unit 13, the elastic information processing unit 14, the compression condition evaluation unit 18, the device control interface unit 19, the cine memory unit 20, and the like, which are characteristic parts of the present invention.

本発明の弾性画像生成方法の一実施例として、図2のBモード画像に示す生体組織を含む部位に圧迫を加えて弾性画像を生成する場合について説明する。実施例1では、弾性情報演算部13において弾性フレームデータとして求められた弾性率を、生体組織の各部位における計測点の弾性情報として用いるものとする。なお、弾性情報としては、弾性率の他、粘弾性率、弾性率の非線形性に係る非線形パラメータ、ヒステリシスなども想定可能である。   As an embodiment of the elastic image generation method of the present invention, a case will be described in which an elastic image is generated by applying pressure to a part including the living tissue shown in the B-mode image of FIG. In the first embodiment, the elastic modulus obtained as the elastic frame data in the elastic information calculation unit 13 is used as the elastic information of the measurement point in each part of the living tissue. As elastic information, in addition to the elastic modulus, a viscoelastic modulus, a nonlinear parameter related to the nonlinearity of the elastic modulus, hysteresis, and the like can be assumed.

まず、圧迫条件評価部18には、装置制御インターフェイス部19から検者によって、関心領域(以下、ROIという。)が入力設定される。その際、ROIは、フレームデータ取得部11により取得されたフレームデータに対応する画像領域内に設定されるが、その対象画像は、図2に示すBモード画像でも図3に示す弾性画像でもよい。例えば、図2には、被検体内の白黒輝度の断層画像に相当する画像が表示された状態を示しており、該断層画像に示された生体組織22が属するようにROIが設定されている。また、図3には、図2の走査面に対応するカラー弾性画像に相当する画像が表示された状態を示しており、該弾性画像に示された生体組織32(図2に示す生体組織22に相当)が属するようにROIが設定されている。
次いで、圧迫条件評価部18は、弾性情報演算部13から各計測点の弾性率データ(弾性率の数値データ)を弾性情報として取り込むとともに、各計測点の応力や歪みなどの数値データを計測圧迫条件として取り込む。なお、実施例1では、圧迫の程度を評価するための計測圧迫条件として応力を用いた場合を示す。
そして、圧迫条件評価部18は、画像表示器10に表示させる弾性画像を生成するための弾性情報を、計測点における現時刻の弾性情報から選択する。併せて、圧迫条件評価部18は、画像表示器10に表示させる弾性画像を生成するための弾性情報を、弾性情報が選択された計測点における現時刻の計測圧迫条件に応じて、弾性情報が選択されなかった計測点における現時刻もしくは過去の弾性情報から該計測点ごとに選択する。例えば、弾性画像を生成するにあたって、ROIに含まれる計測点からは現時刻の弾性情報が計測点ごとに複数選択され、ROIに含まれない計測点からは現時刻もしくは過去の弾性情報が計測点ごとに複数選択される。
First, a region of interest (hereinafter referred to as ROI) is input and set to the compression condition evaluation unit 18 by the examiner from the device control interface unit 19. At that time, the ROI is set in an image region corresponding to the frame data acquired by the frame data acquisition unit 11, and the target image may be a B-mode image shown in FIG. 2 or an elastic image shown in FIG. . For example, FIG. 2 shows a state in which an image corresponding to a tomographic image of black and white luminance in the subject is displayed, and the ROI is set so that the living tissue 22 shown in the tomographic image belongs. . 3 shows a state in which an image corresponding to the color elastic image corresponding to the scanning plane in FIG. 2 is displayed, and the biological tissue 32 (the biological tissue 22 shown in FIG. 2) shown in the elastic image. The ROI is set so that the
Next, the compression condition evaluation unit 18 takes in elastic modulus data (numerical data of elastic modulus) at each measurement point from the elastic information calculation unit 13 as elastic information and measures and compresses numerical data such as stress and strain at each measurement point. Capture as a condition. In Example 1, a case where stress is used as a measurement compression condition for evaluating the degree of compression is shown.
And the compression condition evaluation part 18 selects the elasticity information for producing | generating the elasticity image displayed on the image display 10 from the elasticity information of the present time in a measurement point. In addition, the compression condition evaluation unit 18 generates elasticity information for generating an elasticity image to be displayed on the image display 10 according to the measurement compression condition at the current time at the measurement point where the elasticity information is selected. A selection is made for each measurement point from the current time or past elasticity information at the measurement point that was not selected. For example, when generating an elasticity image, a plurality of elasticity information at the current time is selected for each measurement point from the measurement points included in the ROI, and the current time or past elasticity information is measured from the measurement points not included in the ROI. Multiple selections are made for each.

実施例1の一態様として、ROIに負荷されている応力を求め、その応力下で得た弾性情報に基づいて、ROI内の画像、及び該ROI以外の領域の画像を生成する場合について説明する。
この場合、圧迫条件評価部18は、弾性情報演算部13により求められた計測圧迫条件(一例として、応力)から該計測圧迫条件を代表する代表計測圧迫条件を求め、計測点における現時刻もしくは過去の計測圧迫条件が代表計測圧迫条件を満たすか否かを判定し、代表計測圧迫条件を満たす計測点ごとに、現時刻もしくは過去の弾性情報(弾性率データ)を弾性画像を生成するための弾性情報として選択する。なお、代表計測圧迫条件は、弾性情報演算部13で求められた計測圧迫条件を代表する応力の条件であって、弾性情報演算部13から圧迫条件評価部18へ取り込まれた応力が、生体組織の弾性情報(弾性情報演算部13で求められた弾性率など)を画像化させるべき状態となっているか否かを判定する際の閾値となるものである。
As one aspect of the first embodiment, a case where the stress applied to the ROI is obtained and an image in the ROI and an image of an area other than the ROI are generated based on the elasticity information obtained under the stress will be described. .
In this case, the compression condition evaluation unit 18 obtains a representative measurement compression condition representative of the measurement compression condition from the measurement compression condition (for example, stress) obtained by the elasticity information calculation unit 13, and obtains the current time or the past at the measurement point. To determine whether or not the measurement compression condition satisfies the representative measurement compression condition, and to generate an elasticity image of the current time or past elasticity information (elastic modulus data) for each measurement point that satisfies the representative measurement compression condition Select as information. The representative measurement compression condition is a stress condition representative of the measurement compression condition obtained by the elasticity information calculation unit 13, and the stress taken into the compression condition evaluation unit 18 from the elasticity information calculation unit 13 is a biological tissue. This is a threshold for determining whether or not the elasticity information (such as the elastic modulus obtained by the elasticity information calculation unit 13) should be imaged.

実施例1において、圧迫条件評価部18は、すべての計測点のそれぞれにつき、計測された弾性率(一例として、ヤング率)とその計測時の応力区間(上限値と下限値で示される応力の範囲)の対応関係を示す情報(以下、ヤング率−応力の対応情報という。)を保存するようになっている。
例えば、応力区間は、0から1kPaまでを10分割した応力区間1〜10として設定し、計測点Pi,jのそれぞれの応力区間で得られたヤング率の値を以下のようにEi,j,mで表す。なお、Ei,j,mにおいて、i,jは計測点Pi,jの位置(例えば、図3に示す画素30の座標位置)、mは応力区間をそれぞれ示すパラメータである。
Ei,j,1:応力0kPa以上、0.1kPa未満の条件で得られたヤング率
Ei,j,2:応力0.1kPa以上、0.2kPa未満の条件で得られたヤング率
・・・
Ei,j,9:応力0.8kPa以上、0.9kPa未満の条件で得られたヤング率
Ei,j,10:応力0.9kPa以上、1.0kPa未満の条件で得られたヤング率
In Example 1, the compression condition evaluation unit 18 measures the measured elastic modulus (for example, Young's modulus) and the stress interval at the time of measurement (the stress indicated by the upper limit value and the lower limit value) for each of all measurement points. Information indicating the correspondence relationship of the range) (hereinafter referred to as Young's modulus-stress correspondence information) is stored.
For example, the stress interval is set as stress intervals 1 to 10 obtained by dividing 10 from 0 to 1 kPa, and the values of Young's modulus obtained in the respective stress intervals at the measurement point Pi, j are expressed as Ei, j, Represented by m. Note that in Ei, j, m, i, j is the position of the measurement point Pi, j (for example, the coordinate position of the pixel 30 shown in FIG. 3), and m is a parameter indicating the stress interval.
Ei, j, 1: Young's modulus obtained under conditions of stress 0 kPa or more and less than 0.1 kPa Ei, j, 2: Young's modulus obtained under conditions of stress 0.1 kPa or more and less than 0.2 kPa
Ei, j, 9: Young's modulus obtained under conditions of stress 0.8 kPa or more and less than 0.9 kPa Ei, j, 10: Young's modulus obtained under conditions of stress 0.9 kPa or more and less than 1.0 kPa

実施例1においては、代表計測圧迫条件をROI内の計測点における現時刻の応力の平均値に基づいて求める。したがって、圧迫条件評価部18は、現時刻におけるROIの領域に負荷される圧迫の大きさとして、応力を求める。例えば、現時刻におけるROI内の計測点に負荷する応力の平均値を含む応力区間を代表計測圧迫条件とし、この平均値がσroi(0)=0.86kPaと求められたとする。なお、求められた代表計測圧迫条件に対応するROI内の応力σroi(0)、及びその時のヤング率の値は、図2に示すBモード画像上や図3に示す弾性画像上に数値で表示させれば、これらを検者が容易に把握可能となる。   In the first embodiment, the representative measurement compression condition is obtained based on the average value of the stress at the current time at the measurement point in the ROI. Therefore, the compression condition evaluation unit 18 obtains stress as the size of the compression applied to the ROI region at the current time. For example, it is assumed that a stress interval including an average value of stress applied to a measurement point in the ROI at the current time is a representative measurement compression condition, and this average value is obtained as σroi (0) = 0.86 kPa. The stress σ roi (0) in the ROI corresponding to the obtained representative measurement compression condition and the value of the Young's modulus at that time are displayed numerically on the B-mode image shown in FIG. 2 or the elastic image shown in FIG. By doing so, the examiner can easily grasp these.

このように代表計測圧迫条件が求められると、次に、圧迫条件評価部18は、その応力σroi(0)の大きさに一致する応力下で得られた(換言すれば、応力σroi(0)を含む応力区間に対応する)ヤング率の値を、ヤング率−応力の対応情報を参照して求める。
例えば、応力0.86kPaは応力区間9に相当するので、計測点Pi,jのそれぞれにおけるヤング率として、Ei,j,9(0)の値を参照することで、代表計測圧迫条件下での各計測点Pi,jのヤング率の値を求めることができる。そして、各計測点Pi,jに対し、求めた値に応じて白黒輝度もしくは色相に階調を与えてヤング率のマッピングを行う。
When the representative measurement compression condition is thus obtained, the compression condition evaluation unit 18 is obtained under a stress that matches the magnitude of the stress σroi (0) (in other words, the stress σroi (0) A value of Young's modulus (corresponding to a stress section including) is obtained with reference to the Young's modulus-stress correspondence information.
For example, since the stress of 0.86 kPa corresponds to the stress section 9, referring to the value of Ei, j, 9 (0) as the Young's modulus at each measurement point Pi, j, The value of the Young's modulus at each measurement point Pi, j can be obtained. Then, for each measurement point Pi, j, the Young's modulus is mapped by giving gradation to the monochrome luminance or hue according to the obtained value.

ただし、ある計測点において現時刻における応力が代表計測圧迫条件に対応する値となっている(以下、この状態を代表計測圧迫条件を満たすという。)としても、別の計測点における現時刻の応力が代表計測圧迫条件に対応する値となっていない(同、代表計測圧迫条件を満たしていないという。)場合も起こり得る。例えば、計測点Pi,jにおいて、現時刻における応力が0.8kPa以上、0.9kPa未満ではない(応力が応力区間9にはない)場合、Ei,j,9(0)にはヤング率の値が得られていない。つまり、参照すべきヤング率−応力の対応情報が存在しない。この場合には、かかる計測点Pi,jのヤング率の画像には白黒輝度及び色相のいずれも与えられないようになっている。   However, even if the stress at the current time at a certain measurement point has a value corresponding to the representative measurement compression condition (hereinafter, this state is referred to as the representative measurement compression condition), the stress at the current time at another measurement point. May not be a value corresponding to the representative measurement compression condition (it is said that the representative measurement compression condition is not satisfied). For example, at the measurement point Pi, j, if the stress at the current time is not less than 0.8 kPa and less than 0.9 kPa (the stress is not in the stress section 9), Ei, j, 9 (0) has a Young's modulus. The value is not obtained. That is, there is no correspondence information of Young's modulus-stress to be referred to. In this case, neither the monochrome luminance nor the hue is given to the image of the Young's modulus at the measurement point Pi, j.

図3には、このようなヤング率のマッピング表示の一例が示されており、被検体内の計測点Pi,jにおけるヤング率の値が該計測点Pi,jに対応する画素30ごとに色調コードを付与したカラー画像として表示されている。この場合、生体組織32(図2に示す生体組織22に相当)が属するようにROIが設定され、代表計測圧迫条件を現時刻において満たす計測点Pi,jに対応する画素30がヤング率の値に基づいてカラー表示されている。すなわち、図3に示す弾性画像は、代表計測圧迫条件を現時刻において満たす計測点のヤング率が弾性画像を生成するための弾性情報として選択され、選択されたヤング率の値に基づいてかかる計測点に対応する画素30がカラー表示された画像となっている。その際、代表計測圧迫条件を現時刻において満たさない計測点は、そのヤング率が弾性画像を生成するための弾性情報として選択されない。   FIG. 3 shows an example of mapping display of such Young's modulus. The value of Young's modulus at the measurement point Pi, j in the subject is a color tone for each pixel 30 corresponding to the measurement point Pi, j. It is displayed as a color image with a code attached. In this case, the ROI is set so that the living tissue 32 (corresponding to the living tissue 22 shown in FIG. 2) belongs, and the pixel 30 corresponding to the measurement point Pi, j that satisfies the representative measurement compression condition at the current time is the value of Young's modulus. Color display is based on. That is, in the elasticity image shown in FIG. 3, the Young's modulus of the measurement point that satisfies the representative measurement compression condition at the current time is selected as elasticity information for generating the elasticity image, and the measurement based on the value of the selected Young's modulus. The pixel 30 corresponding to the point is an image displayed in color. At that time, a measurement point that does not satisfy the representative measurement compression condition at the current time is not selected as elasticity information for generating an elasticity image.

一例として、図3では、ヤング率の値がより大きな該計測点Pi,jほど、対応する画素30が濃色で表示されている。なお、図3においては、便宜上、ヤング率の大きさに応じて画素30ごとに白黒の濃淡を付与しているが、実際の画面上には画素30ごとに色調コードを付与したカラー画像が表示される。
したがって、図3の表示画像によれば、生体組織32に対応する領域(同図において、最も濃色で示された領域)は、代表計測圧迫条件下で生体組織32の周辺組織(同図において、生体組織32よりも淡色で示された領域)34,36よりも現時刻において高いヤング率となっていることを把握することができる。一方、かかる生体組織32の周辺組織34,36は、ROIとの重複領域を含め、現時刻において生体組織32よりも低いヤング率となっており、生体組織32から離れるに従って周辺組織34、周辺組織36の順で徐々にヤング率が低下していくことが把握できる。また、生体組織32及びその周辺組織34,36以外の組織(図3において、白色で示された領域)38は、代表計測圧迫条件を現時刻においていずれも満たしておらず、かかる組織内の計測点におけるヤング率の値が計測されていない状態であることが把握できる。
このように、ROIに負荷する現時刻の応力σroi(0)と同じ応力下で得られたヤング率の値に基づいた弾性画像を生成し、リアルタイムで表示させることができる。
As an example, in FIG. 3, the corresponding pixel 30 is displayed in a darker color at the measurement point Pi, j having a larger Young's modulus value. In FIG. 3, for the sake of convenience, black and white shading is given to each pixel 30 according to the Young's modulus, but a color image with a tone code assigned to each pixel 30 is displayed on the actual screen. Is done.
Therefore, according to the display image of FIG. 3, the region corresponding to the living tissue 32 (the region shown in the darkest color in the drawing) is the surrounding tissue (shown in the drawing) of the living tissue 32 under the representative measurement compression condition. It is possible to grasp that the Young's modulus is higher at the current time than the regions 34 and 36, which are lighter than the living tissue 32. On the other hand, the surrounding tissues 34 and 36 of the living tissue 32 have a Young's modulus lower than that of the living tissue 32 at the current time including the overlapping area with the ROI, and the surrounding tissues 34 and surrounding tissues as the distance from the living tissue 32 increases. It can be seen that the Young's modulus gradually decreases in the order of 36. Further, tissues other than the living tissue 32 and the surrounding tissues 34 and 36 (regions shown in white in FIG. 3) do not satisfy the representative measurement compression conditions at the current time, and the measurement within the tissue is not performed. It can be understood that the value of the Young's modulus at the point is not measured.
In this way, an elastic image based on the value of Young's modulus obtained under the same stress as the current time stress σroi (0) loaded on the ROI can be generated and displayed in real time.

ここで、上述した態様(以下、態様1-1という。)では、現時刻において代表計測圧迫条件を満たす計測点におけるヤング率の値に基づいた弾性画像を生成しているため、現時刻の応力が代表計測圧迫条件を満たしていない計測点については、白黒輝度もしくは色相で階調化した表示がなされない場合がある。そこで、現時刻だけでなく、過去における計測点の弾性率(一例として、ヤング率)の情報も保存しておき、この保存した情報に基づいて弾性画像を生成可能とする実施例1の態様(以下、態様1-2という。)について説明する。   Here, in the above-described aspect (hereinafter referred to as aspect 1-1), since an elastic image is generated based on the Young's modulus value at the measurement point that satisfies the representative measurement compression condition at the current time, the stress at the current time is generated. However, the measurement points that do not satisfy the representative measurement compression condition may not be displayed with gradation in black and white luminance or hue. Therefore, not only the current time but also information on the elastic modulus (as an example, Young's modulus) of the measurement point in the past is stored, and an elastic image can be generated based on the stored information (Example 1) Hereinafter, Mode 1-2 will be described.

この場合、現時刻の計測結果に加え、過去における複数回の計測結果を保存できるバッファなどを定義する。一例として、計測点Pi,jにおける応力区間ごとの現時刻及び過去10回までの計測結果を以下のようにEi,j,m(n)で表す。なお、Ei,j,m(n)において、i,jは計測点Pi,jの位置(例えば、図3に示す画素30の座標位置)、mは応力区間、nは過去のヤング率の計測結果が現時刻よりも何回前の計測によって得られたかをそれぞれ示すパラメータである。
Ei,j,1(0):計測点Pi,jにおいて、応力区間1で現時刻の計測で得られたヤング率の計測結果
Ei,j,1(1):計測点Pi,jにおいて、応力区間1で現時刻よりも1回前の計測で得られたヤング率の計測結果
Ei,j,1(2):計測点Pi,jにおいて、応力区間1で現時刻よりも2回前の計測で得られたヤング率の計測結果
・・・
Ei,j,1(10):計測点Pi,jにおいて、応力区間1で現時刻よりも10回前の計測で得られたヤング率の計測結果
・・・
そして、すべての計測点Pi,jにおいて、割り当てられた複数の応力区間のそれぞれの区間で得られた過去のヤング率の計測結果をバッファに保存し、参照可能とする。例えば、応力区間を上述した態様1-1と同様の10区間とした場合、現時刻及び過去10回の計測における11時点の計測点Pi,jのヤング率の値が、応力区間ごとに合計110パターン保存可能となる。
In this case, a buffer or the like that can store a plurality of past measurement results in addition to the measurement result at the current time is defined. As an example, the current time for each stress section at the measurement point Pi, j and the measurement results up to the past 10 times are represented by Ei, j, m (n) as follows. In Ei, j, m (n), i, j is the position of the measurement point Pi, j (for example, the coordinate position of the pixel 30 shown in FIG. 3), m is the stress interval, and n is the past Young's modulus measurement. Each parameter indicates how many times before the current time the measurement was obtained.
Ei, j, 1 (0): Measurement result of Young's modulus obtained by measurement at the current time in stress section 1 at measurement point Pi, j Ei, j, 1 (1): Stress at measurement point Pi, j Measurement result of Young's modulus obtained in measurement 1 time before current time in section 1 Ei, j, 1 (2): Measurement at measurement point Pi, j twice before current time in stress section 1 Measurement result of Young's modulus obtained in ...
Ei, j, 1 (10): Measurement result of Young's modulus obtained by measurement 10 times before the current time in stress section 1 at measurement point Pi, j
Then, at all measurement points Pi, j, the past Young's modulus measurement results obtained in each of the plurality of assigned stress sections are stored in a buffer and can be referred to. For example, when the stress interval is 10 intervals similar to the above-described aspect 1-1, the Young's modulus values of the measurement points Pi, j at 11 points in the current time and the past 10 measurements are 110 in total for each stress interval. Pattern can be saved.

上述したように態様1-1では、計測点Pi,jにおける現時刻の応力が、現時刻のROI内の応力σroi(0)の大きさに相当する応力区間に該当しない(つまり、代表計測圧迫条件を満たさない)場合、計測点Pi,jにおけるヤング率の値は得られない。これに対し、本態様1-2では、過去の計測で現時刻の代表計測圧迫条件を満たす応力値が計測され、その応力区間において得られたヤング率の値が保存されていれば、そのヤング率の値に基づいた画像を表示させることが可能である。   As described above, in aspect 1-1, the stress at the measurement point Pi, j at the current time does not correspond to the stress interval corresponding to the magnitude of the stress σroi (0) in the ROI at the current time (that is, representative measurement compression) If the condition is not satisfied, the Young's modulus value at the measurement point Pi, j cannot be obtained. On the other hand, in this aspect 1-2, if the stress value satisfying the representative measurement compression condition at the present time is measured in the past measurement and the Young's modulus value obtained in the stress interval is stored, the Young An image based on the rate value can be displayed.

例えば、現時刻においてROI内に負荷する応力の平均値が、σroi(0)=0.86kPaであったとする。計測点Pi,jにおけるヤング率として、応力区間9において過去の計測で得られた値が保存されているか否かを検索し、保存されていれば、その最新の結果(直近のヤング率)であるEi,j,9(1)の値を参照し、その値に応じた白黒輝度もしくは色相で階調化したヤング率のマッピングが行われる。なお、このように現時刻、もしくは直近のヤング率に対してマッピングを行う場合、応力区間が上述した態様1-1と同様の10区間であれば、現時刻及び過去1回の計測における2時点の計測点Pi,jのヤング率の値を、応力区間ごとに合計20パターン保存可能なバッファを定義しておけばよい。すなわち、かかるバッファ定義は、ヤング率のマッピング態様に応じて適宜変更することが可能である。   For example, it is assumed that the average value of stress loaded in the ROI at the current time is σroi (0) = 0.86 kPa. As the Young's modulus at the measurement point Pi, j, whether or not the value obtained in the past measurement is stored in the stress section 9 is searched, and if it is stored, the latest result (the latest Young's modulus) is used. A value of Ei, j, 9 (1) is referred to, and mapping of the Young's modulus is performed with gradation in black and white luminance or hue according to the value. In addition, when mapping is performed for the current time or the latest Young's modulus in this way, if the stress interval is 10 intervals similar to the above-described aspect 1-1, the current time and 2 points in the past measurement are measured. A buffer capable of storing a total of 20 patterns of Young's modulus values of the measurement points Pi, j for each stress section may be defined. That is, such a buffer definition can be changed as appropriate according to the mapping of Young's modulus.

このように本態様1-2では、現時刻の応力が代表計測圧迫条件を満たす領域のみならず、代表計測圧迫条件を満たしていない領域についても、ROI内に負荷される応力と同じ圧迫条件下(つまり、代表計測圧迫条件)で得られたヤング率の値をマッピングした一枚の画像を生成することができる。したがって、ROIを含め、表示画像(観察画像)の全領域について、ヤング率の値に基づいた弾性画像を生成して表示させることも可能となる。換言すれば、この場合、同じ圧迫条件下(代表計測圧迫条件)で得られた各計測点Pi,jにおける現時刻もしくは過去のヤング率の計測結果を全ての計測点についてマッピングし、一枚の画像に合成して画像化がなされている。   As described above, in this aspect 1-2, not only the region where the stress at the current time satisfies the representative measurement compression condition but also the region where the representative measurement compression condition is not satisfied, the same compression condition as the stress loaded in the ROI It is possible to generate a single image in which Young's modulus values obtained under (representative measurement compression conditions) are mapped. Therefore, it is possible to generate and display an elastic image based on the Young's modulus value for the entire region of the display image (observation image) including the ROI. In other words, in this case, the measurement result of the current time or the past Young's modulus at each measurement point Pi, j obtained under the same compression condition (representative measurement compression condition) is mapped for all measurement points. The image is synthesized with the image.

図4には、このように生成された画像の一表示例が示されており、被検体内の計測点Pi,jにおけるヤング率の値が該計測点Pi,jに対応する画素30ごとに色調コードを付与したカラー画像として表示されている。この場合、生体組織32(図2に示す生体組織22に相当)が属するようにROIが設定され、代表計測圧迫条件を現時刻において満たす計測点もしくは過去において満たした計測点に対応する画素30がヤング率の値に基づいてカラー表示されている。図4では、一例として、ヤング率の値がより大きな該計測点Pi,jほど、対応する画素30が濃色で表示されている。なお、図4においては、便宜上、ヤング率の大きさに応じて画素30ごとに白黒の濃淡を付与しているが、実際の画面上には画素30ごとに色調コードを付与したカラー画像が表示される。   FIG. 4 shows a display example of the image generated in this way, and the Young's modulus value at the measurement point Pi, j in the subject is determined for each pixel 30 corresponding to the measurement point Pi, j. It is displayed as a color image with a color code added. In this case, the ROI is set so that the living tissue 32 (corresponding to the living tissue 22 shown in FIG. 2) belongs, and the pixel 30 corresponding to the measurement point satisfying the representative measurement compression condition at the current time or the measurement point satisfied in the past. Color display is based on the Young's modulus value. In FIG. 4, as an example, the corresponding pixel 30 is displayed in a darker color at the measurement point Pi, j having a larger Young's modulus value. In FIG. 4, for convenience, black and white shading is given to each pixel 30 according to the Young's modulus, but a color image to which a tone code is assigned to each pixel 30 is displayed on the actual screen. Is done.

ここで、図4に示す弾性画像は、上述した態様1-1で生成された図3に示す弾性画像、すなわち、代表計測圧迫条件を現時刻において満たす計測点に対応する画素30がヤング率の値に基づいてカラー表示された画像と、図5に示す弾性画像のイメージ、すなわち、代表計測圧迫条件を過去において満たした計測点に対応する画素30がヤング率の値に基づいてカラー表示された画像が合成されて画像化された状態に相当する。したがって、かかる画像化にあたっては、図4に示す弾性画像データをカラースキャンコンバータ15へ送出させ、切替加算部9を介して画像表示器10に表示させればよいが、次のような合成処理によっても画像化が可能である。この場合、圧迫条件評価部18は、弾性情報処理部14に指令を送って、図3に示す弾性画像の画像データと図5に示す弾性画像イメージの画像データをカラースキャンコンバータ15へ送出させる。カラースキャンコンバータ15では、それぞれの画像データに対して画素30ごとに色調コードを付与してカラー画像が生成される。そして、生成されたこれらのカラー画像を切替加算部9で加算合成し、画像表示器10に重ねて表示させればよい。   Here, the elasticity image shown in FIG. 4 is the elasticity image shown in FIG. 3 generated in the above-described aspect 1-1, that is, the pixel 30 corresponding to the measurement point that satisfies the representative measurement compression condition at the current time has a Young's modulus. The image displayed in color based on the value and the image of the elastic image shown in FIG. 5, that is, the pixel 30 corresponding to the measurement point satisfying the representative measurement compression condition in the past was displayed in color based on the Young's modulus value. This corresponds to a state in which images are synthesized and imaged. Therefore, in such imaging, the elastic image data shown in FIG. 4 may be sent to the color scan converter 15 and displayed on the image display device 10 via the switching addition unit 9. Can also be imaged. In this case, the compression condition evaluation unit 18 sends a command to the elastic information processing unit 14 to send the image data of the elastic image shown in FIG. 3 and the image data of the elastic image image shown in FIG. 5 to the color scan converter 15. In the color scan converter 15, a color image is generated by assigning a tone code for each pixel 30 to each image data. Then, these generated color images may be added and combined by the switching addition unit 9 and displayed on the image display 10 in an overlapping manner.

なお、上述したように、図3に示す弾性画像は、代表計測圧迫条件を現時刻において満たす計測点のヤング率が弾性画像を生成するための弾性情報として選択され、選択されたヤング率の値に基づいてかかる計測点に対応する画素30がカラー表示された画像となっている。
これに対し、図5には、代表計測圧迫条件を現時刻においては満たさないが、過去において満たした計測点におけるその時のヤング率が弾性画像を生成するための弾性情報として選択され、選択されたヤング率の値に基づいてかかる計測点に対応する画素30に色調コードが付与された弾性画像のイメージが示されている。ただし、かかる弾性画像自体が画像表示器10に表示される訳ではない。この場合、図5において白色で示された領域は、生体組織32及びその周辺組織34,36(図3,4)に相当し、かかる領域は、現時刻において代表計測圧迫条件を満たしているため、計測点に対応する画素30には色調コードが付与されていない状態となっている。一方、図5において2つの淡色で示された領域は、周辺組織36の一部と、生体組織32及びその周辺組織34,36以外の組織38(図3,4)に相当し、かかる領域は、現時刻において代表計測圧迫条件を満たさないが、過去において満たしているため、計測点に対応する画素30には色調コードが付与された状態となっている。なお、図5には、周辺組織36の一部と、生体組織32及びその周辺組織34,36以外の組織38(図3,4)に相当する領域の全域が過去において代表計測圧迫条件を満たした状態を一例として示しているが、かかる領域の一部が過去においても代表計測圧迫条件を満たしていない場合、その計測点に対応する画素30には色調コードは付与されない。
As described above, in the elasticity image shown in FIG. 3, the Young's modulus of the measurement point that satisfies the representative measurement compression condition at the current time is selected as the elasticity information for generating the elasticity image, and the value of the selected Young's modulus is selected. The pixel 30 corresponding to the measurement point is an image in which color display is performed.
On the other hand, in FIG. 5, the representative measurement compression condition is not satisfied at the current time, but the Young's modulus at the measurement point satisfied in the past is selected and selected as elasticity information for generating an elasticity image. An image of an elastic image in which a color tone code is given to the pixel 30 corresponding to the measurement point based on the Young's modulus value is shown. However, the elastic image itself is not displayed on the image display 10. In this case, the area shown in white in FIG. 5 corresponds to the living tissue 32 and the surrounding tissues 34 and 36 (FIGS. 3 and 4), and this area satisfies the representative measurement compression condition at the current time. The color code is not assigned to the pixel 30 corresponding to the measurement point. On the other hand, two light-colored regions in FIG. 5 correspond to a part of the surrounding tissue 36 and the tissue 38 (FIGS. 3 and 4) other than the living tissue 32 and the surrounding tissues 34 and 36. The representative measurement compression condition is not satisfied at the current time, but since it has been satisfied in the past, the color code is assigned to the pixel 30 corresponding to the measurement point. In FIG. 5, a part of the surrounding tissue 36 and the entire region corresponding to the tissue 38 other than the living tissue 32 and the surrounding tissues 34 and 36 (FIGS. 3 and 4) satisfy the representative measurement compression condition in the past. The state is shown as an example. However, if a part of the region does not satisfy the representative measurement compression condition in the past, no color code is assigned to the pixel 30 corresponding to the measurement point.

このようにして生成された図4に示す表示画像によれば、生体組織32に対応する領域(同図において、最も濃色で示された領域)は、代表計測圧迫条件下で生体組織32の周辺組織(同図において、生体組織32よりも淡色で示された領域)34,36,38よりも高いヤング率となることを把握することができる。一方、かかる生体組織32の周辺組織34,36,38は、ROIとの重複領域を含め、生体組織32よりも低いヤング率となっており、生体組織32から離れるに従って周辺組織34、周辺組織36、周辺組織38の順で徐々にヤング率が低下していくことが把握できる。このように、図4においては、表示画像(観察画像)の全領域(生体組織32及び各周辺組織34,36,38)について、ヤング率の値に基づいた弾性画像が生成され、表示された状態となっている。   According to the display image shown in FIG. 4 generated in this way, the region corresponding to the living tissue 32 (the region shown in the darkest color in the drawing) is the region of the living tissue 32 under the representative measurement compression condition. It can be understood that the Young's modulus is higher than that of the surrounding tissues (regions shown in a lighter color than the living tissue 32 in the figure) 34, 36, 38. On the other hand, the surrounding tissues 34, 36, and 38 of the living tissue 32 have a lower Young's modulus than the living tissue 32 including the overlapping region with the ROI, and the surrounding tissues 34 and 36 are separated from the living tissue 32. It can be understood that the Young's modulus gradually decreases in the order of the surrounding tissue 38. As described above, in FIG. 4, the elasticity image based on the Young's modulus value is generated and displayed for the entire region of the display image (observation image) (the living tissue 32 and the surrounding tissues 34, 36, and 38). It is in a state.

なお、本態様1-2では、上述したように、各計測点Pi,jにおいて過去の計測で得られた、それぞれの応力区間における直近の弾性情報(一例として、ヤング率)の値に基づいて画像を生成する方法を示したが、圧迫条件評価部18において、弾性情報(一例として、ヤング率)が選択された計測点Pi,jにおける過去の弾性情報の計測結果を平滑処理して、もしくは現時刻及び過去の弾性情報の計測結果を平滑処理して、該平滑処理後の弾性情報に基づいて弾性画像を画像表示器10に表示させてもよい。
例えば、保存されている各応力区間において過去に計測された複数時点のヤング率の平均値を算出し、その平均値に基づいて画像化を行う場合を想定し、この場合に現時刻においてROI内に負荷する応力(代表計測圧迫条件に対応する応力)が、σroi(0)=0.86kPaであったとする。計測点Pi,jにおける応力区間9で得られたヤング率の過去の計測結果の平均値Ei,j,9(mean)は、
Ei,j,9(mean)=(Ei,j,9(1)+Ei,j,9(2)+Ei,j,9(3)+・・・+Ei,j,9(10))/10
により算出されるため、この値に応じた白黒輝度もしくは色相で階調化したヤング率のマッピングを行えばよい。なお、バッファ定義は、どの程度まで過去の計測結果を参照するかに応じて適宜変更すればよい。
In addition, in this aspect 1-2, as mentioned above, based on the value of the latest elasticity information (for example, Young's modulus) in each stress section obtained by past measurement at each measurement point Pi, j. Although the method for generating the image has been shown, the compression condition evaluation unit 18 smoothes the measurement result of the past elasticity information at the measurement point Pi, j where the elasticity information (for example, Young's modulus) is selected, or The measurement result of the current time and past elasticity information may be smoothed, and the elasticity image may be displayed on the image display 10 based on the elasticity information after the smoothing.
For example, it is assumed that an average value of Young's moduli at a plurality of points measured in the past in each stored stress interval is calculated, and imaging is performed based on the average value. Is assumed to be σroi (0) = 0.86 kPa (stress corresponding to the representative measurement compression condition). The average value Ei, j, 9 (mean) of past measurement results of Young's modulus obtained in the stress section 9 at the measurement point Pi, j is
Ei, j, 9 (mean) = (Ei, j, 9 (1) + Ei, j, 9 (2) + Ei, j, 9 (3) + ... + Ei, j, 9 (10)) / 10
Therefore, the mapping of the Young's modulus gradationd with the black and white luminance or hue according to this value may be performed. It should be noted that the buffer definition may be changed as appropriate depending on how far past measurement results are referred to.

さらに、該当する応力区間において現時刻におけるヤング率の値が得られていれば、現時刻のヤング率の値と保存されている過去のヤング率の計測結果との平均値を求めて、その平均値に応じた画像化を行っても構わない。
また、このような画像化のための処理としては、上記のようなFIR(Finite-duration Impulse Response)型の平滑処理に限らず、IIR(Infinite-duration Impulse Response)型の平滑処理に基づいた処理でもよく、過去のヤング率の計測結果の中央値(メジアン)を求める処理であっても構わない。
このように過去の弾性情報との平滑化処理を行うことで、より安定した高精度なヤング率の画像を生成することが可能となる。
Furthermore, if the value of Young's modulus at the current time is obtained in the corresponding stress section, the average value of the Young's modulus value at the current time and the stored measurement result of the past Young's modulus is obtained, and the average Imaging according to the value may be performed.
In addition, the imaging process is not limited to the above-described FIR (Finite-duration Impulse Response) type smoothing process, but is based on an IIR (Infinite-duration Impulse Response) type smoothing process. Alternatively, it may be a process for obtaining a median (median) of past Young's modulus measurement results.
Thus, by performing the smoothing process with the past elasticity information, it becomes possible to generate a more stable and highly accurate image of Young's modulus.

上述した実施例1の各態様では、検者が関心領域(ROI)を設定し、そのROI内に負荷する応力(代表計測圧迫条件に対応する応力)に基づいて弾性画像を生成する方法を示したが、圧迫条件評価部18において、画像表示器10に表示させる弾性画像を生成するための弾性情報(一例として、ヤング率)を、画像領域内のすべての計測点における現時刻の計測圧迫条件(一例として、応力)に応じて、計測点のすべてにおける現時刻もしくは過去の弾性情報から計測点ごとに選択することも可能である。すなわち、画像表示器10に表示されている弾性画像の全領域(つまり、観察している画像の全領域)をROIとし、全領域における応力(例えば、全領域の計測点に負荷する応力の平均値)に基づいて弾性画像を生成しても構わない。この場合、特にROIを設定する作業は不要である。例えば、腫瘍などの疑われる部位がないかどうか広範囲の領域をスクリーニングする検査など、特にROIの位置が固定されないような場合においては、検査効率の向上を図ることができる。   In each aspect of the first embodiment described above, a method is shown in which the examiner sets a region of interest (ROI) and generates an elasticity image based on the stress loaded in the ROI (stress corresponding to the representative measurement compression condition). However, the compression condition evaluation unit 18 uses the elastic information (as an example, Young's modulus) for generating an elastic image to be displayed on the image display 10 as the measurement compression condition at the current time at all measurement points in the image area. It is also possible to select for each measurement point from the current time or past elasticity information at all the measurement points according to (as an example, stress). That is, the entire area of the elastic image displayed on the image display 10 (that is, the entire area of the image being observed) is ROI, and the stress in the entire area (for example, the average stress applied to the measurement points in the entire area) The elastic image may be generated based on the (value). In this case, it is not particularly necessary to set the ROI. For example, particularly in the case where the position of the ROI is not fixed, such as a test for screening a wide area for a suspected site such as a tumor, the test efficiency can be improved.

また、上述した実施例1の各態様においては、ROIに属する生体組織が圧迫操作により移動することについて考慮していない。そこで、圧迫操作により移動する生体組織に追従してROIを移動させることで、弾性情報の評価精度を向上させることを可能とする実施例1の別態様(以下、態様1-3という。)について説明する。
例えば、画像に対するROIの位置を固定すると、圧迫操作の時間経過に応じてROIに属する生体組織が変化してしまうことがあり、ヤング率などの計測値の精度が悪くなる。すなわち、図6(a)から(c)に示すように、生体組織62への圧迫操作を時刻t0から時刻t2まで行った場合、時刻t0の時点では生体組織62に属するようにROIが設定されている(同図(a)に示すROI(t0)の状態)。その後、圧迫操作の時間経過とともに、生体組織62の移動によってROIが生体組織64の方向へずれ始め、時刻t1の時点では生体組織62のみならず生体組織64もROIに属するようになる(同図に示すROI(t1)の状態)。そして、時刻t2の時点に至っては、生体組織64に加えて生体組織66もROIに属するようになってしまう(同図(c)に示すROI(t2)の状態)。
Moreover, in each aspect of Example 1 mentioned above, it does not consider about the biological tissue which belongs to ROI moving by compression operation. Therefore, another aspect of the first embodiment (hereinafter referred to as aspect 1-3) that can improve the evaluation accuracy of the elasticity information by moving the ROI following the living tissue that is moved by the compression operation. explain.
For example, if the position of the ROI with respect to the image is fixed, the living tissue belonging to the ROI may change with the passage of time of the compression operation, and the accuracy of measurement values such as Young's modulus will deteriorate. That is, as shown in FIGS. 6A to 6C, when the compression operation on the living tissue 62 is performed from time t0 to time t2, the ROI is set so as to belong to the living tissue 62 at the time t0. (The state of ROI (t0) shown in FIG. 4A). Thereafter, with the passage of time of the compression operation, the ROI starts to shift toward the living tissue 64 due to the movement of the living tissue 62, and not only the living tissue 62 but also the living tissue 64 belongs to the ROI at the time t1 (FIG. State of ROI (t1) shown in FIG. At time t2, in addition to the living tissue 64, the living tissue 66 also belongs to the ROI (state of ROI (t2) shown in FIG. 10C).

このため、本態様1-3では、圧迫の程度に応じて、ROIに属する生体組織の移動先に追従してROIを移動させることにより、任意の圧迫条件の下において、同一生体組織の情報のみの抽出を可能としている。すなわち、図6(d)から(f)に示すように、生体組織62への圧迫操作を時刻t0から時刻t2まで行った場合であっても、生体組織62の時間経過による移動に追従してROIをROI(t0)、ROI(t1)、ROI(t2)の位置へそれぞれ移動させていくことで、ROIには常に生体組織62のみが属する状態を保っている(図6(d)から(f)に示す状態)。
この結果、弾性情報(一例として、ヤング率)の評価精度を向上させることが可能となる。なお、生体組織の動きに追従させてROIを移動させる処理は、上述した生体組織の変位ベクトルを求めるブロックマッチング法又はグラジェント法を適用すればよい。また、表示画像(観察画像)の全領域をROIとした場合であっても、生体組織の移動に応じて追従させることができる。
その際、同一組織点における過去のヤング率の計測結果の平均値を算出するなど、平滑処理を行ってもよい。
また、新規更新処理の場合には、現在保存されているヤング率などの弾性情報をトラッキングに応じて移動させていく処理が必要になる。
For this reason, in this aspect 1-3, according to the degree of compression, by moving the ROI following the movement destination of the living tissue belonging to the ROI, only information on the same living tissue is obtained under an arbitrary compression condition. Extraction is possible. That is, as shown in FIGS. 6D to 6F, even when the compression operation to the living tissue 62 is performed from time t0 to time t2, the movement of the living tissue 62 over time is followed. By moving the ROI to the positions of ROI (t0), ROI (t1), and ROI (t2), only the living tissue 62 always remains in the ROI (from FIG. 6 (d) to ( State shown in f).
As a result, it is possible to improve the evaluation accuracy of elasticity information (for example, Young's modulus). In addition, what is necessary is just to apply the block matching method or the gradient method which calculates | requires the displacement vector of the biological tissue mentioned above to the process which follows the motion of a biological tissue and moves ROI. Further, even when the entire region of the display image (observation image) is an ROI, it can be followed according to the movement of the living tissue.
At that time, smoothing processing may be performed such as calculating an average value of past Young's modulus measurement results at the same tissue point.
Further, in the case of the new update process, it is necessary to perform a process of moving currently stored elasticity information such as Young's modulus according to tracking.

なお、実施例1では、一例として、応力区間の範囲を0.1kPaの幅で10区間に設定した場合について説明したが、割り当てる応力区間の範囲をより狭くし、区間数をより多くして設定することで(例えば、0.01kPa幅で100区間など)、ROI外の領域のヤング率をより高精度に求めることができる。
また、実施例1では、計測圧迫条件、及び該計測圧迫条件から求められる代表計測圧迫条件として、計測点及びROI内の応力を用いているが、例えば、計測点及びROI内の歪み、変位計測部12で計測された変位量(変位ゼロからの積算を含む)、及び圧力計測部17で計測された圧迫の圧力に基づいて導出されたROI内に伝わった応力など、圧迫の状態を反映した各種の情報を用いることも可能である。
In addition, in Example 1, although the case where the range of the stress interval was set to 10 intervals with a width of 0.1 kPa was described as an example, the range of the stress interval to be allocated is made narrower and the number of intervals is set larger. By doing so (for example, 0.01 kPa width and 100 sections, etc.), the Young's modulus of the region outside the ROI can be obtained with higher accuracy.
In the first embodiment, the stress in the measurement point and the ROI is used as the measurement compression condition and the representative measurement compression condition obtained from the measurement compression condition. Reflects the state of compression, such as the amount of displacement measured by the unit 12 (including integration from zero displacement) and the stress transmitted in the ROI derived based on the pressure of the compression measured by the pressure measuring unit 17 Various types of information can also be used.

上述したように、実施例1では、生体組織の弾性情報(一例として、ヤング率)を画像表示させるか否かを、実際に生体組織を圧迫した状態で、圧迫条件評価部18において求めた代表計測圧迫条件によって判定して選択したが、かかる画像化可否の判定基準はこれに限定されない。例えば、このような画像化可否を検者が設定した基準値によって判定することも可能である。以下、画像化可否をこのように検者が設定した基準値によって判定する場合を、実施例2として説明する。なお、実施例2に係る超音波診断装置の基本的な構成は、上述した実施例1(図1)と同様であり、以下では実施例2の特有構成についての説明に止める。   As described above, in the first embodiment, whether or not the elasticity information of the living tissue (for example, Young's modulus) is displayed as an image is determined by the compression condition evaluation unit 18 in a state where the living tissue is actually pressed. Although it is determined and selected according to the measurement compression condition, the criterion for determining whether or not imaging is possible is not limited to this. For example, it is possible to determine whether or not such imaging is possible based on a reference value set by the examiner. Hereinafter, the case where the image availability is determined based on the reference value set by the examiner will be described as a second embodiment. The basic configuration of the ultrasonic diagnostic apparatus according to the second embodiment is the same as that of the above-described first embodiment (FIG. 1), and only the specific configuration of the second embodiment will be described below.

実施例2において、圧迫条件評価部18は、計測点における現時刻もしくは過去の計測圧迫条件(一例として、応力)が予め設定された基準圧迫条件を満たすか否か判定し、基準圧迫条件を満たす計測点ごとに、現時刻もしくは過去の弾性情報(一例として、ヤング率)を弾性画像を生成するための弾性情報として選択する。
したがって、実施例2においては、まず、圧迫条件評価部18に対し、装置制御インターフェイス部19から検者によって基準圧迫条件が入力設定される。かかる基準圧迫条件は、検者が実際の診断に用いるための生体組織のヤング率を画像化させるべく、応力や歪みなどの計測圧迫条件に対して予め設定される閾値となるものである。設定された基準圧迫条件は、図2に示すBモード画像上、もしくは、図3に示す弾性画像上に表示されるようになっている。
そして、圧迫条件評価部18は、弾性情報演算部13から各計測点の弾性率データ(弾性率の数値データ)と、各計測点の応力の数値データを取り込む。なお、実施例2では、圧迫の程度を評価するための計測圧迫条件として応力を用いた場合を示す。
In the second embodiment, the compression condition evaluation unit 18 determines whether the current time at the measurement point or the past measurement compression condition (for example, stress) satisfies a preset reference compression condition, and satisfies the reference compression condition. For each measurement point, the current time or past elasticity information (for example, Young's modulus) is selected as elasticity information for generating an elasticity image.
Therefore, in the second embodiment, first, the reference compression condition is input and set by the examiner from the device control interface unit 19 to the compression condition evaluation unit 18. The reference compression condition is a threshold set in advance for measurement compression conditions such as stress and strain in order to image the Young's modulus of the living tissue for use by the examiner in actual diagnosis. The set reference compression condition is displayed on the B-mode image shown in FIG. 2 or the elastic image shown in FIG.
Then, the compression condition evaluation unit 18 takes in the elastic modulus data (numerical data of elastic modulus) at each measurement point and the numerical data of stress at each measurement point from the elastic information calculation unit 13. Example 2 shows a case where stress is used as a measurement compression condition for evaluating the degree of compression.

ここでは、実施例2として、基準となる応力範囲を設定し、この範囲内で計測された弾性情報(一例として、ヤング率)に基づいて、表示(観察)している全領域の画像を生成する場合について説明する。
この場合、圧迫条件評価部18は、装置制御インターフェイス部19から検者によって入力された基準圧迫条件に基づいて、画像化可否の基準となる応力範囲を設定している。すなわち、入力された基準応力σsにある程度の幅を持たせた応力範囲、例えば、かかる基準応力σsを跨いで±10%の応力範囲(0.9×σs〜1.1×σs)を基準圧迫条件として設定する。その際、基準応力σsに対する幅の割合は特に限定されない。なお、基準圧迫条件として、このような応力範囲を直接入力してもよい。あるいは、応力範囲ではなく、入力した基準応力σsをそのまま基準圧迫条件として設定することも可能である。
Here, as Example 2, a reference stress range is set, and an image of the entire region being displayed (observed) is generated based on elasticity information (for example, Young's modulus) measured within this range. The case where it does is demonstrated.
In this case, the compression condition evaluation unit 18 sets a stress range that is a reference for whether or not imaging is possible based on the reference compression condition input by the examiner from the device control interface unit 19. That is, a stress range in which the input reference stress σs is given a certain width, for example, a ± 10% stress range (0.9 × σs to 1.1 × σs) across the reference stress σs is used as the reference compression. Set as a condition. At that time, the ratio of the width to the reference stress σs is not particularly limited. Note that such a stress range may be directly input as the reference compression condition. Alternatively, it is possible to set the input reference stress σs as the reference compression condition as it is, not the stress range.

このように基準圧迫条件が設定されると、次に、圧迫条件評価部18は、計測点において現時刻もしくは過去における応力が基準圧迫条件に対応する値となっているか否かを判定し、かかる応力が上記応力範囲(0.9×σs〜1.1×σs)に入っているとき、そのときに計測されたヤング率の値を白黒輝度もしくは色相で階調化し、弾性画像としてマッピングする。そして、圧迫条件評価部18は、弾性情報処理部14に指令を送って、カラースキャンコンバータ15から出力されるカラー弾性画像を画像表示器10に出力表示させる(例えば、図3に示すような弾性画像が表示された状態に相当)。なお、以下の説明では、計測点において現時刻もしくは過去における応力が基準圧迫条件に対応する値となっている状態を基準圧迫条件を満たすといい、かかる応力が基準圧迫条件に対応する値となっていない状態を基準圧迫条件を満たさないという。   When the reference compression condition is set in this way, the compression condition evaluation unit 18 next determines whether or not the current time or past stress at the measurement point has a value corresponding to the reference compression condition. When the stress is in the above stress range (0.9 × σs to 1.1 × σs), the Young's modulus value measured at that time is gray-scaled by monochrome luminance or hue and mapped as an elastic image. Then, the compression condition evaluation unit 18 sends a command to the elasticity information processing unit 14 to output and display the color elasticity image output from the color scan converter 15 on the image display 10 (for example, elasticity as shown in FIG. 3). Equivalent to an image being displayed). In the following explanation, it is said that the state where the stress at the measurement point at the current time or the past has a value corresponding to the reference compression condition is referred to as the reference compression condition, and the stress corresponds to the reference compression condition. It is said that it does not meet the standard compression condition.

これにより、被検体の生体組織に対する圧迫の強さを変化させることで、表示(観察)している全領域を関心領域として、基準圧迫条件を満たす計測点のヤング率に基づいた弾性画像を生成し、リアルタイムで表示させることができる。その際、圧迫操作を繰り返して行った場合に、ある計測点における応力が基準圧迫条件を再度満たしたときは、そのときに計測されたヤング率の値で弾性画像を更新させることも可能である。   As a result, by changing the strength of the compression on the living tissue of the subject, an elastic image based on the Young's modulus of the measurement point that satisfies the reference compression condition is generated with the entire region being displayed (observed) as the region of interest. And can be displayed in real time. At that time, when the compression operation is repeated, if the stress at a certain measurement point satisfies the reference compression condition again, it is possible to update the elastic image with the value of the Young's modulus measured at that time. .

なお、基準応力σsにおける計測点のヤング率を推定するために、マイナス側10%の応力範囲(0.9×σs〜σs)におけるヤング率の計測結果と、プラス側10%の応力範囲(σs〜1.1×σs)におけるヤング率の計測結果から、その間にある基準応力σsにおけるヤング率を、例えば線形補間によって求めることも可能である。この方法によれば、基準圧迫条件におけるヤング率の計測結果をより高精度に求めることができる。
その際、基準応力σsに最も近い過去のヤング率の計測結果2点を、マイナス側10%の応力範囲(0.9×σs〜σs)とプラス側10%の応力範囲(σs〜1.1×σs)からそれぞれ1点ずつ抽出し、これらから補間したヤング率の値を計測結果として保持してもよい。さらに、基準応力σsにより近い計測結果が得られた際に随時、かかる補間結果を更新するようになっていても構わない。
In order to estimate the Young's modulus of the measurement point at the reference stress σs, the measurement result of the Young's modulus in the negative 10% stress range (0.9 × σs to σs) and the positive 10% stress range (σs It is also possible to obtain the Young's modulus at the reference stress σs between them by linear interpolation, for example, from the measurement result of Young's modulus at ˜1.1 × σs). According to this method, the measurement result of the Young's modulus under the reference compression condition can be obtained with higher accuracy.
At that time, the past two Young's modulus measurement results closest to the reference stress σs are obtained by dividing the negative 10% stress range (0.9 × σs to σs) and the positive 10% stress range (σs to 1.1). Xσs) may be extracted one by one, and the Young's modulus value interpolated therefrom may be held as a measurement result. Furthermore, when a measurement result closer to the reference stress σs is obtained, the interpolation result may be updated as needed.

また、表示(観察)している全領域ではなく、特定の関心領域(ROI)を設定し、該ROI内の基準圧迫条件を満たす計測点におけるヤング率(一例として、その平均値)に基づいた弾性画像を生成して表示させることも可能である。その際には、かかるヤング率の平均値を図2に示すBモード画像上や図3に示す弾性画像上に表示させてもよい。この表示は、ROI内の計測点が基準圧迫条件を再度満たした際にヤング率の平均値が更新されるまで継続させ、その後は更新された平均値を表示させればよい。   In addition, a specific region of interest (ROI) is set instead of the entire region being displayed (observed), and based on Young's modulus (as an example, an average value) at a measurement point that satisfies the reference compression condition in the ROI It is also possible to generate and display an elastic image. At that time, the average value of the Young's modulus may be displayed on the B-mode image shown in FIG. 2 or the elastic image shown in FIG. This display may be continued until the average value of the Young's modulus is updated when the measurement point in the ROI satisfies the reference compression condition again, and then the updated average value may be displayed.

ここで、実施例2においても、上述した実施例1の態様1-1と同様に、ある計測点において現時刻における応力が基準圧迫条件を満たすとしても、別の計測点における現時刻の応力が基準圧迫条件を満たしていない場合もあり、この場合、基準圧迫条件を満たさない計測点には白黒輝度及び色相のいずれも与えず、ヤング率の画像を表示しないようになっている(図3に示すような表示状態に相当)。
すなわち、実施例2においても、現時刻においてのみ基準圧迫条件を満たす計測点におけるヤング率の値に基づいた弾性画像を生成した場合には、現時刻の応力が基準圧迫条件を満たしていない計測点については画像表示されないこととなる。
したがって、上述した実施例1の態様1-2と同様に、現時刻だけでなく、計測点において過去に基準圧迫条件を満たしたときの弾性率(一例として、ヤング率)の情報を保存しておくことで、該計測点において基準圧迫条件を満たした直近のヤング率の値、あるいは過去に計測された複数時点のヤング率の平均値などに基づいて弾性画像を生成することが可能となる。これにより、表示画像(観察画像)の全領域について、ヤング率の値に基づいた弾性画像を生成して表示させることも可能となる(例えば、図4に示すような弾性画像が表示された状態に相当)。なお、実施例2においても、上述した実施例1の態様1-2と同様に、現時刻の計測結果に加えて、過去の複数回の計測結果を保存できるバッファなどを定義しておけばよい。
Here, also in Example 2, as in the case 1-1 of Example 1 described above, even if the stress at the current time satisfies the reference compression condition at a certain measurement point, the stress at the current time at another measurement point is In some cases, the reference compression condition is not satisfied, and in this case, the measurement point that does not satisfy the reference compression condition is not provided with either monochrome luminance or hue, and the Young's modulus image is not displayed (see FIG. 3). Equivalent to the display state shown).
That is, also in Example 2, when an elastic image is generated based on the Young's modulus value at a measurement point that satisfies the reference compression condition only at the current time, the measurement point at which the stress at the current time does not satisfy the reference compression condition No image is displayed for.
Therefore, similarly to the above-described embodiment 1-2 of Example 1, not only the current time but also information on the elastic modulus (for example, Young's modulus) when the reference compression condition has been satisfied in the past at the measurement point is stored. Thus, it is possible to generate an elastic image based on the latest Young's modulus value that satisfies the reference compression condition at the measurement point, or the average value of Young's modulus measured at a plurality of points in the past. Thereby, it is also possible to generate and display an elastic image based on the Young's modulus value for the entire region of the display image (observation image) (for example, a state in which an elastic image as shown in FIG. 4 is displayed). Equivalent). In the second embodiment, a buffer that can store a plurality of past measurement results in addition to the measurement result at the current time may be defined in the same manner as in the first and second aspects of the first embodiment. .

上述した実施例2において、例えば、生成した弾性画像を表示させるとともに、計測圧迫条件と弾性情報の関係を示すグラフ画像を作成し、画像表示器10に表示させ、基準圧迫条件を表示させたグラフ画像上で設定してもよい。一例として、図7に示すような生体組織の応力と弾性率の関係を表示させ、この表示に基づいて基準圧迫条件の設定を行うことが可能である。これにより、生成した弾性画像により所定の応力下における生体組織の弾性をリアルタイムで確認しつつ、基準圧迫条件を調整しながら設定することができる。   In the second embodiment described above, for example, the generated elasticity image is displayed, and a graph image showing the relationship between the measurement compression condition and the elasticity information is created and displayed on the image display 10 to display the reference compression condition. It may be set on the image. As an example, it is possible to display the relationship between the stress and elastic modulus of the living tissue as shown in FIG. 7, and to set the reference compression condition based on this display. Accordingly, it is possible to set the reference compression condition while adjusting the elasticity of the living tissue under a predetermined stress in real time from the generated elasticity image.

この場合、圧迫条件評価部18は、各計測点の応力の数値データと各計測点の弾性率データ(弾性率の数値データ)に基づいて、これら応力と弾性率の関係をプロットしてグラフ画像を作成し、画像表示器10に表示させる。図7は、作成したグラフ画像の一表示例であり、この場合、非線形性の大きい組織と小さい組織について、応力と弾性率の関係をそれぞれ表示させている。
図7に示すように、横軸である応力軸上には、その軸方向へスライド可能なスライドバー70が表示されている。かかるスライドバー70は、基準圧迫条件を示しており、このスライドバー70を応力軸に沿ってスライドさせることで、基準圧迫条件をスライドバー70の応力軸上の所望の値に設定することができるようになっている。そして、スライドバー70の位置によって設定した基準圧迫条件を満たす計測点における弾性情報(一例として、ヤング率)の計測結果がリアルタイムにマッピングされ、弾性画像として表示される。
したがって、設定した基準圧迫条件に応じてリアルタイムに表示される弾性画像を非線形性の大きい組織と小さい組織についてそれぞれ確認しつつ、基準圧迫条件を調整することができ、検者が実際の診断に用いるために最適な生体組織のヤング率の状態を容易かつ確実に画像化させることが可能となる。
In this case, the compression condition evaluation unit 18 plots the relationship between the stress and the elastic modulus based on the numerical data of the stress at each measurement point and the elastic modulus data (the elastic modulus numerical data) at each measurement point. Is generated and displayed on the image display 10. FIG. 7 is a display example of the created graph image. In this case, the relationship between the stress and the elastic modulus is displayed for the tissue with a large nonlinearity and the tissue with a small nonlinearity.
As shown in FIG. 7, a slide bar 70 that can slide in the axial direction is displayed on the stress axis that is the horizontal axis. The slide bar 70 indicates a reference compression condition, and the reference compression condition can be set to a desired value on the stress axis of the slide bar 70 by sliding the slide bar 70 along the stress axis. It is like that. Then, the measurement result of the elasticity information (for example, Young's modulus) at the measurement point that satisfies the reference compression condition set by the position of the slide bar 70 is mapped in real time and displayed as an elasticity image.
Therefore, it is possible to adjust the reference compression condition while checking the elasticity image displayed in real time according to the set reference compression condition for each of the tissue having a large nonlinearity and the tissue having a small nonlinearity, and the examiner uses it for the actual diagnosis. Therefore, it is possible to easily and reliably image the optimal Young's modulus state of the living tissue.

なお、スライドバーは、設定される基準圧迫条件に応じて、上述した基準応力σsにある程度の幅を持たせた応力範囲(例えば、かかる基準応力σsを跨ぐ±10%の応力範囲)を示すように応力軸上の所定範囲をカバーして表示させてもよいし、基準応力σsそのものを示すように応力軸と交差して表示させてもよい。また、基準応力σsとその応力範囲をそれぞれ個別に設定することが可能となるように、複数のスライドバーを表示させてもよい。一例として、図7には、基準応力σsにある程度の幅を持たせた応力範囲を表すスライドバー70を示す。かかるスライドバー70においては、同図中に破線で表示する中心線が基準応力σsの値を示し、該中心線を挟んで応力軸方向の両側に実線で表示する境界線が応力範囲の上限値と下限値を示している。
いずれの場合であっても、スライドバー70の操作は、検者が装置制御インターフェイス部19などを介して手動で行えばよいが、圧迫条件評価部18などにより自動的に行っても構わない。
The slide bar shows a stress range in which the above-described reference stress σs is given a certain width according to the set reference compression condition (for example, a stress range of ± 10% across the reference stress σs). A predetermined range on the stress axis may be covered and displayed, or may be displayed crossing the stress axis so as to indicate the reference stress σs itself. Also, a plurality of slide bars may be displayed so that the reference stress σs and the stress range can be individually set. As an example, FIG. 7 shows a slide bar 70 representing a stress range in which a certain width is given to the reference stress σs. In such a slide bar 70, the center line indicated by a broken line in the figure indicates the value of the reference stress σs, and the boundary line indicated by a solid line on both sides in the stress axis direction across the center line indicates the upper limit value of the stress range. And the lower limit.
In either case, the operation of the slide bar 70 may be performed manually by the examiner via the device control interface unit 19 or the like, but may be automatically performed by the compression condition evaluation unit 18 or the like.

また、スライドバー70を応力軸に沿ってスライドさせた際、上記のような弾性画像に加えてもしくは代えて、ROI内や表示(観察)している全領域のヤング率の値の分布をヒストグラムとして表示させてもよい。この場合、かかるヒストグラムの統計的な特徴量(平均値や分散値など)も併せて解析し、ヒストグラムと同時に表示させても構わない。   Further, when the slide bar 70 is slid along the stress axis, in addition to or in place of the elastic image as described above, the distribution of Young's modulus values in the ROI and the entire region being displayed (observed) is a histogram. May be displayed. In this case, the statistical feature amount (average value, variance value, etc.) of the histogram may also be analyzed and displayed simultaneously with the histogram.

さらに、上述した実施例2において、圧迫条件評価部18に対して装置制御インターフェイス部19から検者によってROIが入力設定されると、該ROI内のヤング率の値が得られ、このROI内の応力とヤング率の関係がプロットされるとともに、最小二乗法などにより最もよく近似する曲線が表示され、その非線形性を表す非線形パラメータなどの特徴量が解析されるようになっていてもよい。   Furthermore, in Example 2 described above, when the ROI is input and set by the examiner from the device control interface unit 19 to the compression condition evaluation unit 18, the value of the Young's modulus in the ROI is obtained. The relationship between stress and Young's modulus may be plotted, and a curve that best approximates may be displayed by a least square method or the like, and a feature quantity such as a nonlinear parameter representing the nonlinearity may be analyzed.

超音波診断装置において表示画像をフリーズさせたタイミングで、弾性画像の生成を切換制御することも可能であり、以下、このように画像生成の切換をフリーズによって行う場合を、実施例3として説明する。なお、実施例3に係る超音波診断装置の基本的な構成は、上述した実施例1(図1)及び実施例2と同様であり、以下では実施例3の特有構成についての説明に止める。   It is also possible to switch the generation of the elastic image at the timing when the display image is frozen in the ultrasonic diagnostic apparatus. Hereinafter, a case where the image generation is switched by the freeze will be described as a third embodiment. . The basic configuration of the ultrasonic diagnostic apparatus according to the third embodiment is the same as that of the first embodiment (FIG. 1) and the second embodiment described above, and only the specific configuration of the third embodiment will be described below.

実施例3において、圧迫条件評価部18は、画像表示器10に表示させた弾性画像がフリーズされた際、画像表示器10に表示させる弾性画像を生成するための弾性情報を新たに選択し、新たに選択された弾性情報に基づいて生成された弾性画像を画像表示器10に表示させる。例えば、実施例3では、上述した実施例1(態様1-2)及び実施例2において、現時刻における代表計測圧迫条件、もしくは基準圧迫条件を満たした計測点からなる領域のヤング率の値に基づいた画像をリアルタイムで表示させ続ける。そして、超音波診断装置にフリーズ指令を与えて表示画像をフリーズさせた際に、計測点において過去に代表計測圧迫条件、もしくは基準圧迫条件を満たした時のヤング率の保存情報で、現時刻には表示されなかった計測点からなる領域に対してヤング率の値に基づいて弾性画像を生成するようになっている。   In Example 3, the compression condition evaluation unit 18 newly selects elasticity information for generating an elasticity image to be displayed on the image display 10 when the elasticity image displayed on the image display 10 is frozen. The elastic image generated based on the newly selected elasticity information is displayed on the image display 10. For example, in Example 3, the Young's modulus value of the region consisting of measurement points satisfying the representative measurement compression condition or the reference compression condition at the current time in Example 1 (Aspects 1-2) and Example 2 described above is used. Continue to display the based image in real time. Then, when a freeze command is given to the ultrasonic diagnostic apparatus to freeze the display image, the Young's modulus storage information when the representative measurement compression condition or the reference compression condition has been met at the measurement point in the past is stored at the current time. Is configured to generate an elastic image based on the Young's modulus value for an area consisting of measurement points that are not displayed.

また、実施例3において、圧迫条件評価部18は、画像表示器10に表示させた弾性画像がフリーズされた際、該フリーズまでに生成された弾性画像を画像表示器10に連続的に読み出させることで、弾性画像を動画として画像表示器10に表示させることも可能となっている。この場合、超音波診断装置にフリーズ指令を与えた際に、フリーズまでに取得された時間的に連続した所定フレーム数の弾性画像がシネメモリ部20に格納される。そして、シネメモリ部20に格納された弾性画像を連続的に読み出すことで、画像表示器10に動画として表示することができるようになっている。   In the third embodiment, when the elastic image displayed on the image display 10 is frozen, the compression condition evaluation unit 18 continuously reads the elastic image generated up to the freeze to the image display 10. By doing so, the elastic image can be displayed on the image display 10 as a moving image. In this case, when a freeze command is given to the ultrasonic diagnostic apparatus, the cine-memory unit 20 stores a predetermined number of temporally continuous elasticity images acquired until the freeze. The elastic image stored in the cine memory unit 20 is continuously read out so that it can be displayed on the image display 10 as a moving image.

なお、このようなフリーズ指令を与えた際に加えて、もしくはこれに代えて、表示可能な全画像領域に対して一定の面積以上の弾性画像が生成された時、自動的にその弾性画像(静止画)がシネメモリ部20に格納されるようになっていてもよい。   In addition to or in place of giving such a freeze command, when an elastic image having a certain area or more is generated for all displayable image areas, the elastic image ( Still image) may be stored in the cine memory unit 20.

上述した実施例1、実施例2及び実施例3では、現時刻でリアルタイムに得られた弾性情報(一例として、ヤング率)に加え、過去の弾性情報の計測結果も表示される弾性画像に反映される。したがって、保存されている過去の弾性情報と現時刻の弾性情報が、同一の計測断面(超音波走査面)から得られた弾性情報であること、換言すれば、かかる現時刻及び過去の弾性情報が同一の計測点から得られたことを担保する必要がある。このため、本実施例では、圧迫条件評価部18において、計測点における現時刻の弾性情報と過去の弾性情報が同一の計測点から得られたことを、例えば、一対のフレームデータの類似性を評価することにより、もしくは探触子2の変位をモニタリングすることにより担保している。そこで、かかる計測断面の同一性を担保する方法として、連続する2フレーム間の画像の類似性を評価する方法と、探触子2の変位をモニタリングする方法を採用した場合について、以下にそれぞれ説明する。なお、実施例4に係る超音波診断装置の基本的な構成は、上述した実施例1(図1)、実施例2及び実施例3と同様である。   In the first embodiment, the second embodiment, and the third embodiment described above, in addition to the elasticity information obtained in real time at the current time (for example, Young's modulus), the measurement results of past elasticity information are also reflected in the displayed elasticity image. Is done. Therefore, the stored past elasticity information and the elasticity information at the current time are the elasticity information obtained from the same measurement section (ultrasonic scanning plane), in other words, the current time and the past elasticity information. It is necessary to ensure that is obtained from the same measurement point. For this reason, in the present embodiment, the compression condition evaluation unit 18 determines that the current elasticity information and the past elasticity information at the measurement point are obtained from the same measurement point, for example, the similarity between a pair of frame data. This is ensured by evaluating or monitoring the displacement of the probe 2. Therefore, as a method for ensuring the identity of the measurement cross section, a method for evaluating the similarity of images between two consecutive frames and a method for monitoring the displacement of the probe 2 are described below. To do. The basic configuration of the ultrasonic diagnostic apparatus according to the fourth embodiment is the same as that of the first embodiment (FIG. 1), the second embodiment, and the third embodiment.

相関係数を用いて計測断面が変わっていないことをモニタリングする画像認識技術における方法の一つとして、例えば、2つの画像の類似性を相関係数を用いて評価する手法がある。この方法によれば、相関係数が大きいほど相関が高く、相関係数が小さいほど相関が低いと評価することができる。
したがって、現時刻と1フレーム過去の2フレームのBモード画像(つまり、計測断面)の相関係数を演算し、その相関係数が設定したある閾値よりも大きければ両画像に相関があり、計測断面の同一性が担保されているものと判定することができる。これに対し、かかる相関係数が前記閾値よりも小さければ前記両画像に相関がなく、計測断面は変化し、その同一性が担保されていないものと判定することができる。これをフレームが更新される度に繰り返すことで、現時刻と過去の複数フレームのBモード画像の相関を判定し、これら計測断面の同一性を継続的に担保することが可能となる。なお、このような処理は、例えば、圧迫条件評価部18で行えばよいが、その他の処理手段やそのための別途の処理手段を設けて行っても構わない。
As one of the methods in the image recognition technique for monitoring that the measurement cross section does not change using the correlation coefficient, for example, there is a method of evaluating the similarity between two images using the correlation coefficient. According to this method, it can be evaluated that the correlation is higher as the correlation coefficient is larger, and the correlation is lower as the correlation coefficient is smaller.
Therefore, the correlation coefficient of the B-mode image (that is, the measurement cross section) of the two frames of the current time and one frame in the past is calculated, and if the correlation coefficient is greater than a certain threshold value, both images are correlated and measurement is performed. It can be determined that the identity of the cross section is secured. On the other hand, if the correlation coefficient is smaller than the threshold value, it can be determined that there is no correlation between the two images, the measurement cross section changes, and the identity is not ensured. By repeating this every time a frame is updated, it is possible to determine the correlation between the current time and the B-mode images of a plurality of past frames, and to continuously ensure the identity of these measurement sections. Such processing may be performed by, for example, the compression condition evaluation unit 18, but may be performed by providing other processing means or a separate processing means therefor.

また、探触子2に装着したセンサから得られる変位情報を用いて計測断面の同一性を担保する方法も考えられる。
最近の超音波診断技術の一つとして、超音波診断装置の探触子2にセンサを装着し、現時刻における探触子2の位置、変位、角度、加速度などの情報をリアルタイムに検出することが可能となっている。
したがって、本技術を応用すれば、現時刻の探触子2の位置と角度の変化をモニタリングすることにより、計測断面が保持されているのか、変化したのかを判定することが可能である。
この場合、センサとしては磁気センサを想定することができるが、赤外線センサなど、位置、変位、角度、加速度などの情報を取得可能であれば、その種類は特に問わない。
また、ビデオカメラにより撮影された画像から動作を認識するモーション認識技術に基づいて、計測断面の同一性を担保する方法であっても構わない。
A method of ensuring the identity of the measurement cross section using displacement information obtained from a sensor attached to the probe 2 is also conceivable.
As one of the latest ultrasonic diagnostic technologies, a sensor is attached to the probe 2 of the ultrasonic diagnostic apparatus, and information such as the position, displacement, angle, and acceleration of the probe 2 at the current time is detected in real time. Is possible.
Therefore, if this technique is applied, it is possible to determine whether the measurement cross section is held or has changed by monitoring changes in the position and angle of the probe 2 at the current time.
In this case, a magnetic sensor can be assumed as the sensor, but the type of the sensor is not particularly limited as long as information such as the position, displacement, angle, and acceleration can be acquired.
Moreover, based on the motion recognition technique which recognizes operation | movement from the image image | photographed with the video camera, the method of ensuring the identity of a measurement cross section may be sufficient.

なお、これらの方法によって計測断面が変化して同一性が担保されていないものと判定された場合、過去の計測結果データも含め、すべての計測点における弾性情報(一例として、ヤング率)の計測結果をクリア(初期化)すればよい。これにより、計測断面の同一性の担保を確実に図ることができる。   In addition, when it is determined that the measurement cross section changes by these methods and the identity is not secured, measurement of elasticity information (as an example, Young's modulus) at all measurement points including past measurement result data The result should be cleared (initialized). As a result, it is possible to reliably ensure the identity of the measurement cross sections.

このように、本発明の実施例1〜4によれば、各計測点に加えられる現時刻もしくは現時刻前の計測圧迫条件(代表計測圧迫条件もしくは基準圧迫条件)に応じて計測点ごとに選択された現時刻及び現時刻前の弾性情報(一例として、ヤング率)に基づいて生成された弾性画像を用いているから、同一の大きさの圧迫条件に基づいて得られたヤング率などの弾性情報、換言すれば、組織弾性の応答を画像としてリアルタイムに表示させ、観察することができる。これにより、例えば、生体組織に負荷する圧迫の強さを変えた場合における硬さによる応答の違いなどを、画像により観察することができ、客観的で確定的な組織鑑別が可能となる。また、診断に適した弾性画像の選定作業が不要となり、検者の負担が軽減される。   As described above, according to the first to fourth embodiments of the present invention, the measurement time is selected for each measurement point according to the current time applied to each measurement point or the measurement compression condition (representative measurement compression condition or reference compression condition) before the current time. Elasticity images such as Young's modulus obtained based on compression conditions of the same size are used because the elasticity images generated based on the elasticity information (for example, Young's modulus) before and after the current time are used. Information, in other words, tissue elasticity response can be displayed and observed in real time as an image. Thereby, for example, a difference in response due to hardness when the strength of compression applied to the living tissue is changed can be observed with an image, and objective and definitive tissue differentiation can be performed. Further, it is not necessary to select an elastic image suitable for diagnosis, and the burden on the examiner is reduced.

なお、上述した実施例1〜4においては、対象組織を圧迫することにより、弾性情報を取得するエラストグラフィの手法を一例として説明したが、本発明の弾性画像生成方法は、例えば、加振映像法(Y.Yamakoshi et al : Ultrasonic Imaging of Internal Vibration of Soft Tissue under Forced Vibration, IEEE Trans UFFC 1990; 37; 45-53.)などの手法によって得られた弾性情報に対しても適用することができる。
あるいは、ダイナミックエラストグラフィでも、探触子2をある程度押し込むことで生体内部にはそれに応じた応力分布が生成されるため、本発明の弾性画像生成方法を適用することで同様の効果が得られる。ダイナミックエラストグラフィの技術については、“QUANTITATIVE ASSESSMENT OF BREAST LESION VISCOELASTICITY:INITIAL CLINICAL RESULTS USING SUPERSONIC SHEAR IMAGING” MICKAEL TANTER et al : Ultrasound in Med. & Biol., Vol.34, No.x,pp.xxx,2008 World Federation for Ultrasound in Medicine & Biologyや“ACOUSTIC RADIATION FORCE IMPULSE ELASTOGRAPHY FOR THE EVALUATION OF FOCAL SOLID HEPATIC LESIONS:PRELIMINARY FINDINGS” SEUNG HYUN CHO et al : Ultrasound in Med. & Biol., Vol.36, No.2,pp.202-208,2010 World Federation for Ultrasound in Medicine & Biologyなどを参照できる。
また、例えば、血流像を生成する際のカラードプラも圧迫に依存する。したがって、この場合に本発明の弾性画像生成方法を適用し、基準の圧迫条件を満たしたときのみにカラードプラの色が乗るようにすることも想定可能である。
加えて、応力条件の違いを考慮しないリアルタイムのヤング率の値に基づいた画像と、本発明の弾性画像生成方法による応力条件(代表計測圧迫条件や基準圧迫条件)を満たした計測点におけるヤング率の値に基づいた画像を2画面で表示するようになっていてもよい。
In the first to fourth embodiments described above, the elastography technique for acquiring elasticity information by compressing the target tissue has been described as an example. However, the elastic image generation method of the present invention is, for example, an excitation image. It can also be applied to elasticity information obtained by methods such as Y. Yamakoshi et al: Ultrasonic Imaging of Internal Vibration of Soft Tissue under Forced Vibration, IEEE Trans UFFC 1990; 37; 45-53. .
Alternatively, even in dynamic elastography, since the stress distribution corresponding to the inside of the living body is generated by pushing the probe 2 to some extent, the same effect can be obtained by applying the elastic image generating method of the present invention. Regarding the technology of dynamic elastography, “QUANTITATIVE ASSESSMENT OF BREAST LESION VISCOELASTICITY: INITIAL CLINICAL RESULTS USING SUPERSONIC SHEAR IMAGING” MICKAEL TANTER et al: Ultrasound in Med. & Biol., Vol.34, No.x, pp.xxx, 2008 World Federation for Ultrasound in Medicine & Biology and “ACOUSTIC RADIATION FORCE IMPULSE ELASTOGRAPHY FOR THE EVALUATION OF FOCAL SOLID HEPATIC LESIONS: PRELIMINARY FINDINGS” SEUNG HYUN CHO et al: Ultrasound in Med. & Biol., Vol.36, No.2, pp .202-208,2010 You can refer to World Federation for Ultrasound in Medicine & Biology.
Further, for example, color Doppler when generating a blood flow image also depends on compression. Therefore, in this case, it is also possible to apply the elastic image generation method of the present invention so that the color Doppler color is applied only when the standard compression condition is satisfied.
In addition, images based on real-time Young's modulus values that do not take into account differences in stress conditions, and Young's modulus at measurement points that satisfy the stress conditions (representative measurement compression conditions and reference compression conditions) according to the elastic image generation method of the present invention An image based on the value of may be displayed on two screens.

1 被検体
2 探触子
3 送信回路
4 送受信制御回路
5 受信回路
6 整相加算回路
7 信号処理部
8 白黒スキャンコンバータ
9 切替加算部
10 画像表示器
11 フレームデータ取得部
12 変位計測部
13 弾性情報演算部
14 弾性情報処理部
15 カラースキャンコンバータ
16 圧力センサ
17 圧力計測部
18 圧迫条件評価部
19 装置制御インターフェイス部
20 シネメモリ部
DESCRIPTION OF SYMBOLS 1 Subject 2 Probe 3 Transmission circuit 4 Transmission / reception control circuit 5 Reception circuit 6 Phased addition circuit 7 Signal processing part 8 Black and white scan converter 9 Switching addition part 10 Image display 11 Frame data acquisition part 12 Displacement measurement part 13 Elastic information Calculation unit 14 Elastic information processing unit 15 Color scan converter 16 Pressure sensor 17 Pressure measurement unit 18 Compression condition evaluation unit 19 Device control interface unit 20 Cine memory unit

Claims (14)

探触子から超音波を送信して得られる反射エコー信号に基づいて、画像領域内の複数の計測点における生体組織の弾性情報と前記複数の計測点の計測圧迫条件を求める弾性情報演算部と、
前記弾性情報に基づいて生成された弾性画像を表示する画像表示器と、
前記画像表示器に表示させる前記弾性画像を生成するための前記弾性情報を、現時刻もしくは現時刻前の前記計測圧迫条件に応じて選択し、選択された現時刻及び現時刻前の前記弾性情報に基づいて生成された前記複数の計測点における前記弾性画像を前記画像表示器に表示させる圧迫条件評価部とを備えることを特徴とする超音波診断装置。
Based on a reflected echo signal obtained by transmitting ultrasonic waves from the probe, an elasticity information calculation unit that obtains elasticity information of a living tissue at a plurality of measurement points in an image region and measurement compression conditions at the plurality of measurement points; ,
An image display for displaying an elasticity image generated based on the elasticity information;
The elasticity information for generating the elasticity image to be displayed on the image display is selected according to the measurement compression condition at the current time or before the current time, and the selected current time and the elasticity information before the current time An ultrasonic diagnostic apparatus comprising: a compression condition evaluation unit that causes the image display to display the elasticity image at the plurality of measurement points generated based on the measurement point.
前記圧迫条件評価部は、前記画像表示器に表示させる前記弾性画像を生成するための前記弾性情報を、前記現時刻の弾性情報から選択するとともに、前記弾性情報が選択された前記計測点における現時刻の前記計測圧迫条件に応じて、前記弾性情報が選択されなかった前記計測点における現時刻もしくは現時刻前の前記弾性情報から該計測点ごとに選択することを特徴とする請求項1に記載の超音波診断装置。   The compression condition evaluation unit selects the elasticity information for generating the elasticity image to be displayed on the image display from the elasticity information at the current time, and at the current measurement point where the elasticity information is selected. 2. The measurement point is selected for each measurement point from the current time at the measurement point where the elasticity information is not selected or the elasticity information before the current time according to the measurement compression condition of time. Ultrasound diagnostic equipment. 前記弾性情報が選択された前記計測点は、前記画像領域内に設定した関心領域内に含まれることを特徴とする請求項2に記載の超音波診断装置。   The ultrasonic diagnostic apparatus according to claim 2, wherein the measurement point for which the elasticity information is selected is included in a region of interest set in the image region. 前記圧迫条件評価部は、前記画像表示器に表示させる前記弾性画像を生成するための前記弾性情報を、前記画像領域内のすべての前記計測点における現時刻の前記計測圧迫条件に応じて、前記計測点のすべてにおける現時刻もしくは現時刻前の前記弾性情報から前記計測点ごとに選択することを特徴とする請求項1に記載の超音波診断装置。   The compression condition evaluation unit, the elasticity information for generating the elasticity image to be displayed on the image display, according to the measurement compression conditions at the current time at all the measurement points in the image area, 2. The ultrasonic diagnostic apparatus according to claim 1, wherein selection is made for each measurement point from the elasticity information at the current time or before the current time at all measurement points. 前記圧迫条件評価部は、前記弾性情報演算部により求められた前記計測圧迫条件から該計測圧迫条件を代表する代表計測圧迫条件を求め、前記計測点における現時刻もしくは現時刻前の前記計測圧迫条件が前記代表計測圧迫条件を満たすか否かを判定し、前記代表計測圧迫条件を満たす前記計測点ごとに、現時刻もしくは現時刻前の前記弾性情報を前記弾性画像を生成するための前記弾性情報として選択することを特徴とする請求項1に記載の超音波診断装置。   The compression condition evaluation unit obtains a representative measurement compression condition representing the measurement compression condition from the measurement compression condition obtained by the elasticity information calculation unit, and the measurement compression condition at the current measurement time or before the current time The elasticity information for generating the elasticity image of the elasticity information at the current time or before the current time for each of the measurement points that satisfy the representative measurement compression condition. The ultrasonic diagnostic apparatus according to claim 1, wherein: 前記代表計測圧迫条件は、前記画像領域内に設定した関心領域内の前記計測点もしくは前記画像領域内のすべての前記計測点における現時刻の前記計測圧迫条件の平均値に基づいて求められることを特徴とする請求項5に記載の超音波診断装置。   The representative measurement compression condition is obtained based on an average value of the measurement compression conditions at the current time at the measurement points in the region of interest set in the image region or all the measurement points in the image region. The ultrasonic diagnostic apparatus according to claim 5, wherein the apparatus is an ultrasonic diagnostic apparatus. 前記圧迫条件評価部は、前記計測点における現時刻もしくは現時刻前の前記計測圧迫条件が予め設定された基準圧迫条件を満たすか否か判定し、前記基準圧迫条件を満たす前記計測点ごとに、現時刻もしくは現時刻前の前記弾性情報を前記弾性画像を生成するための前記弾性情報として選択することを特徴とする請求項1に記載の超音波診断装置。   The compression condition evaluation unit determines whether or not the measurement compression condition at a current time or before the current time at the measurement point satisfies a preset reference compression condition, and for each measurement point that satisfies the reference compression condition, 2. The ultrasonic diagnostic apparatus according to claim 1, wherein the elasticity information for generating the elasticity image is selected as the elasticity information at a current time or before the current time. 前記圧迫条件評価部は、前記計測圧迫条件と前記弾性情報の関係を示すグラフ画像を作成し、前記画像表示器に表示させており、
前記基準圧迫条件は、前記グラフ画像上で設定されることを特徴とする請求項7に記載の超音波診断装置。
The compression condition evaluation unit creates a graph image indicating the relationship between the measurement compression condition and the elasticity information, and displays the graph image on the image display device.
The ultrasonic diagnostic apparatus according to claim 7, wherein the reference compression condition is set on the graph image.
前記圧迫条件評価部は、前記弾性情報が選択された前記計測点における現時刻前の前記弾性情報を平滑処理して、もしくは現時刻及び現時刻前の前記弾性情報を平滑処理して、該平滑処理後の前記弾性情報に基づいて前記弾性画像を前記画像表示器に表示させることを特徴とする請求項1から8のいずれかに記載の超音波診断装置。   The compression condition evaluation unit smoothes the elasticity information before the current time at the measurement point where the elasticity information is selected, or smoothes the elasticity information before the current time and the current time, The ultrasonic diagnostic apparatus according to claim 1, wherein the elasticity image is displayed on the image display unit based on the elasticity information after processing. 前記圧迫条件評価部は、前記画像表示器に表示させた前記弾性画像がフリーズされた際、前記画像表示器に表示させる前記弾性画像を生成するための前記弾性情報を新たに選択し、新たに選択された前記弾性情報に基づいて生成された前記弾性画像を前記画像表示器に表示させることを特徴とする請求項1から9のいずれかに記載の超音波診断装置。   The compression condition evaluation unit newly selects the elasticity information for generating the elasticity image to be displayed on the image display when the elasticity image displayed on the image display is frozen, and newly The ultrasonic diagnostic apparatus according to claim 1, wherein the elastic image generated based on the selected elastic information is displayed on the image display. 前記圧迫条件評価部は、前記画像表示器に表示させた前記弾性画像がフリーズされた際、該フリーズまでに生成された前記弾性画像を前記画像表示器に連続的に読み出させることで、前記弾性画像を動画として前記画像表示器に表示させることを特徴とする請求項1から9のいずれかに記載の超音波診断装置。   When the elastic image displayed on the image display is frozen, the compression condition evaluation unit causes the image display to continuously read the elastic image generated up to the freeze. The ultrasonic diagnostic apparatus according to claim 1, wherein an elastic image is displayed on the image display as a moving image. 前記圧迫条件評価部は、前記計測点における現時刻の前記弾性情報と現時刻前の前記弾性情報が同一の前記計測点から得られたことを、一対の前記フレームデータの類似性を評価することにより、もしくは前記探触子の変位をモニタリングすることにより担保することを特徴とする請求項1から9のいずれかに記載の超音波診断装置。   The compression condition evaluation unit evaluates the similarity between a pair of the frame data that the elasticity information at the current time at the measurement point and the elasticity information before the current time are obtained from the same measurement point. The ultrasonic diagnostic apparatus according to claim 1, wherein the ultrasonic diagnostic apparatus is secured by monitoring the displacement of the probe. 前記圧迫条件評価部は、前記画像表示器に表示させる前記弾性画像を生成するための前記弾性情報に白黒輝度又は色相に階調を与えて、前記弾性画像を前記画像表示器に表示させることを特徴とする請求項1から9のいずれかに記載の超音波診断装置。   The compression condition evaluation unit gives the elasticity information to generate black and white luminance or hue to the elasticity information for generating the elasticity image to be displayed on the image display, and causes the image display to display the elasticity image. The ultrasonic diagnostic apparatus according to claim 1, wherein the ultrasonic diagnostic apparatus is characterized in that 生体組織に超音波を送信して得られる反射エコー信号に基づいて、画像領域内の複数の計測点における前記生体組織の弾性情報と前記複数の計測点の計測圧迫条件を求め、
画像表示器に表示させる弾性画像を生成するための前記弾性情報を、現時刻もしくは現時刻前の前記計測圧迫条件に応じて選択し、
選択された現時刻及び現時刻前の前記弾性情報に基づいて前記複数の計測点における前記弾性画像を生成することを特徴とする弾性画像生成方法。
Based on the reflected echo signal obtained by transmitting ultrasonic waves to the biological tissue, the elasticity information of the biological tissue at a plurality of measurement points in the image region and the measurement compression conditions of the plurality of measurement points are obtained,
The elasticity information for generating the elasticity image to be displayed on the image display is selected according to the measurement compression condition at the current time or before the current time,
A method of generating an elasticity image, comprising: generating the elasticity images at the plurality of measurement points based on the selected current time and the elasticity information before the current time.
JP2012053110A 2012-03-09 2012-03-09 Ultrasonic diagnostic apparatus and elastic image generation method Pending JP2013183982A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2012053110A JP2013183982A (en) 2012-03-09 2012-03-09 Ultrasonic diagnostic apparatus and elastic image generation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2012053110A JP2013183982A (en) 2012-03-09 2012-03-09 Ultrasonic diagnostic apparatus and elastic image generation method

Publications (1)

Publication Number Publication Date
JP2013183982A true JP2013183982A (en) 2013-09-19

Family

ID=49385955

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012053110A Pending JP2013183982A (en) 2012-03-09 2012-03-09 Ultrasonic diagnostic apparatus and elastic image generation method

Country Status (1)

Country Link
JP (1) JP2013183982A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2023021036A (en) * 2021-07-29 2023-02-09 スーパー ソニック イマジン Ultrasonic method and system for estimating nonlinear shear wave elasticity of medium
US11995807B2 (en) 2019-05-17 2024-05-28 Supersonic Imagine Method for ultrasound determination of a corrected image of a medium, and device for implementing this method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007282932A (en) * 2006-04-19 2007-11-01 Hitachi Medical Corp Method of generating elastic image and ultrasonograph
JP2008073417A (en) * 2006-09-25 2008-04-03 Hitachi Medical Corp Ultrasonic diagnostic device

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007282932A (en) * 2006-04-19 2007-11-01 Hitachi Medical Corp Method of generating elastic image and ultrasonograph
JP2008073417A (en) * 2006-09-25 2008-04-03 Hitachi Medical Corp Ultrasonic diagnostic device

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11995807B2 (en) 2019-05-17 2024-05-28 Supersonic Imagine Method for ultrasound determination of a corrected image of a medium, and device for implementing this method
JP2023021036A (en) * 2021-07-29 2023-02-09 スーパー ソニック イマジン Ultrasonic method and system for estimating nonlinear shear wave elasticity of medium
JP7297995B2 (en) 2021-07-29 2023-06-26 スーパー ソニック イマジン Ultrasound method and system for estimating nonlinear shear wave elasticity in media

Similar Documents

Publication Publication Date Title
JP4966578B2 (en) Elastic image generation method and ultrasonic diagnostic apparatus
JP4919972B2 (en) Elastic image display method and elastic image display device
JP5160227B2 (en) Ultrasonic diagnostic apparatus and ultrasonic image display method
EP3213108B1 (en) Imaging methods and apparatuses for performing shear wave elastography imaging
JP5303147B2 (en) Ultrasonic diagnostic device for generating elastic images
JP4455003B2 (en) Ultrasonic diagnostic equipment
JP5304986B2 (en) Ultrasonic diagnostic equipment
US8538103B2 (en) Medical image processing device, medical image processing method, medical image diagnostic apparatus, operation method of medical image diagnostic apparatus, and medical image display method
US20090292205A1 (en) Ultrasonic diagnostic apparatus
JP4685633B2 (en) Ultrasonic diagnostic equipment
WO2010024168A1 (en) Ultrasonic diagnosing device
JPWO2007046272A6 (en) Ultrasonic diagnostic device for generating elastic images
KR20120044267A (en) Ultrasound diagnostic apparatus and method for tracing movement of tissue
JP5113322B2 (en) Ultrasonic diagnostic equipment
US20150294457A1 (en) Ultrasound diagnostic apparatus
JP5415669B2 (en) Ultrasonic diagnostic equipment
JP5473527B2 (en) Ultrasonic diagnostic equipment
JP5490979B2 (en) Ultrasonic diagnostic equipment
JP2016112285A (en) Ultrasonic diagnostic device
JP5890358B2 (en) Ultrasonic image pickup apparatus and ultrasonic image display method
JP5623609B2 (en) Ultrasonic diagnostic equipment
JP2013183982A (en) Ultrasonic diagnostic apparatus and elastic image generation method
CN109069117B (en) Ultrasonic diagnostic apparatus
JP5128149B2 (en) Ultrasonic diagnostic equipment
JP5638641B2 (en) Ultrasonic diagnostic equipment

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150306

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20150727

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20150727

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20160113

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160209

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20160512

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20160804