JP2006212445A - Ultrasonographic system - Google Patents

Ultrasonographic system Download PDF

Info

Publication number
JP2006212445A
JP2006212445A JP2006093397A JP2006093397A JP2006212445A JP 2006212445 A JP2006212445 A JP 2006212445A JP 2006093397 A JP2006093397 A JP 2006093397A JP 2006093397 A JP2006093397 A JP 2006093397A JP 2006212445 A JP2006212445 A JP 2006212445A
Authority
JP
Japan
Prior art keywords
image
interest
region
speckle
ultrasonic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2006093397A
Other languages
Japanese (ja)
Other versions
JP4347860B2 (en
Inventor
Noriyoshi Matsushita
典義 松下
Masaru Murashita
賢 村下
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
Aloka Co 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 Aloka Co Ltd filed Critical Aloka Co Ltd
Priority to JP2006093397A priority Critical patent/JP4347860B2/en
Publication of JP2006212445A publication Critical patent/JP2006212445A/en
Application granted granted Critical
Publication of JP4347860B2 publication Critical patent/JP4347860B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To improve a diagnosis of a tissue by quantifying fineness of a speckle (speckle pattern) contained in an ultrasonic image. <P>SOLUTION: A region of interest is set on the ultrasonic image, and binarization is performed to the image in the region of interest while a threshold value is being scanned. The number of islands of high luminance is calculated for the binarized image of each threshold value, and a speckle evaluation graph is formed as its histogram. The speckle evaluation graph reflects tissue properties, and the tissue can be diagnosed by a numerical analysis of the form of the graph. Prior to the binarizing process, it is preferable to perform a process for highlighting the speckle (for processing to remove a base component). It may simultaneously display the image before its process and the image after its process. <P>COPYRIGHT: (C)2006,JPO&NCIPI

Description

本発明は超音波診断装置に関し、特に超音波画像処理に関する。   The present invention relates to an ultrasonic diagnostic apparatus, and more particularly to ultrasonic image processing.

例えば、肝臓の超音波診断を行う場合には、腹部に超音波探触子が当接され、その状態で超音波の送受波がなされる。それにより得られる受信信号に基づいてBモード画像(二次元断層画像)が形成される。医者はそのようなBモード画像を観察することにより腫瘍の有無などを診断する。   For example, when performing ultrasonic diagnosis of the liver, an ultrasonic probe is brought into contact with the abdomen, and ultrasonic waves are transmitted and received in that state. A B-mode image (two-dimensional tomographic image) is formed based on the received signal obtained thereby. The doctor diagnoses the presence or absence of a tumor by observing such a B-mode image.

Bモード画像には、その全体にわたってスペックルというノイズが現れる。Bモード画像上において、ぼやけてやや明るい部分(粒子あるいは島状領域)がスペックルに相当する。   Speckle noise appears throughout the B-mode image. On the B-mode image, a blurry and slightly bright portion (particle or island-like region) corresponds to speckle.

このスペックルは、散乱体からの反射波(エコー)の波の相互干渉などによって引き起こされるものである。Bモード画像の画質を向上させるためには、スペックルをより低減することが望まれるが、Bモード画像には、大なり小なりスペックルが現れる。   This speckle is caused by mutual interference of reflected waves (echoes) from a scatterer. In order to improve the image quality of the B-mode image, it is desired to further reduce speckle, but speckle appears to a greater or lesser extent in the B-mode image.

特表平5−501825号公報JP-T-5-501825

ところで、従来から、Bモード画像上におけるスペックルパターン(多数のスペックルの全体的な模様)と組織の性状との間には、相関関係があることが知られている。つまり、スペックルパターンは常に同じではなく、組織の性質に依存することが知られている。しかしながら、スペックル(あるいはそのパターン)の評価分析は、従来において、医者の目視観察によるものであり、今まで、スペックルの客観的な評価あるいはスペックルの定量化を行える装置は未だ実現されていない(但し、上記特許文献1にはスペックルの分析に関して記載されている)。   Conventionally, it is known that there is a correlation between speckle patterns (overall patterns of speckles) on the B-mode image and tissue properties. That is, it is known that speckle patterns are not always the same, but depend on the nature of the organization. However, speckle (or its pattern) evaluation analysis has conventionally been based on visual observation by a doctor, and until now, an apparatus capable of objectively evaluating speckle or quantifying speckle has not been realized. None (however, the above Patent Document 1 describes the speckle analysis).

本発明は、上記従来の課題に鑑みなされたものであり、その目的は、超音波画像上に登場するスペックル(ノイズ)を積極的に疾病診断に利用することにある。本発明の他の目的は、スペックルの定量化技術を提供することにある。本発明の他の目的は、スペックルパターンの客観的分析を実現することにある。   The present invention has been made in view of the above-described conventional problems, and an object thereof is to positively use speckle (noise) appearing on an ultrasonic image for disease diagnosis. Another object of the present invention is to provide a speckle quantification technique. Another object of the present invention is to realize an objective analysis of speckle patterns.

(1)望ましくは、超音波の送受波により得られたエコーデータに基づいて超音波画像を形成する画像形成手段と、閾値レベルを変化させながら前記超音波画像を二値化処理することにより、複数の二値化画像を生成する二値化処理手段と、前記各二値化画像ごとにラベリング処理を施して高輝度又は低輝度をもった独立領域の個数を求める領域計数手段と、前記各閾値レベルごとの前記独立領域の個数を表すスペックル評価グラフを作成するグラフ作成手段と、が設けられる。 (1) Preferably, by binarizing the ultrasonic image while changing the threshold level, an image forming unit that forms an ultrasonic image based on echo data obtained by transmission and reception of ultrasonic waves, Binarization processing means for generating a plurality of binarized images, area counting means for obtaining a number of independent areas having high luminance or low luminance by applying a labeling process to each of the binarized images, Graph creating means for creating a speckle evaluation graph representing the number of independent regions for each threshold level is provided.

上記構成によれば、生体への超音波の送受波によってエコーデータが取得される。そのエコーデータに基づいて超音波画像(望ましくは二次元断層画像としてのBモード画像であるが、Mモード画像、三次元画像、ドプラ画像などであってもよい)が形成される。その超音波画像に対しては、閾値レベルを段階的にあるいは連続的に変化させつつ二値化処理がなされる。通常は、閾値のスキャンが実行されるが、複数の閾値を同時利用して並列的に二値化処理を行うことも可能である。この二値化処理後の各画像に対して白又は黒の島(独立領域)をカウントするラベリング処理が適用される。このラベリング処理に当たっては公知の各種の技術を利用できる。各画像に対してラベリング処理を行うと、各画像ごとに島の個数が求められることになり、閾値ごとの島の個数を表したヒストグラムとして、スペックル評価グラフが作成される。このスペックル評価グラフは、組織性状を反映したものであり、例えば、組織が正常か否か、腫瘍が悪性か良性か、などの診断に利用できるものである。スペックル評価グラフについては多様な分析方法を適用できる。   According to the said structure, echo data are acquired by the transmission / reception wave of the ultrasonic wave to a biological body. An ultrasonic image (desirably a B-mode image as a two-dimensional tomographic image, but may be an M-mode image, a three-dimensional image, a Doppler image, etc.) is formed based on the echo data. The ultrasound image is binarized while changing the threshold level stepwise or continuously. Normally, a threshold scan is executed, but it is also possible to perform binarization processing in parallel using a plurality of thresholds simultaneously. A labeling process for counting white or black islands (independent areas) is applied to each image after the binarization process. Various known techniques can be used for the labeling process. When the labeling process is performed on each image, the number of islands is obtained for each image, and a speckle evaluation graph is created as a histogram representing the number of islands for each threshold. This speckle evaluation graph reflects tissue properties, and can be used for diagnosis, for example, whether the tissue is normal or whether the tumor is malignant or benign. Various analysis methods can be applied to the speckle evaluation graph.

もちろん、そのスペックル評価グラフをそのまま画像表示してもよいし、そのスペックル評価グラフの分析結果から自動的に組織の性状を診断するようにしてもよい。その場合には、1又は複数のグラフ特徴量と診断内容とを対応付けたデータベースなどを利用すればよい。   Of course, the speckle evaluation graph may be displayed as an image, or the tissue characteristics may be automatically diagnosed from the analysis result of the speckle evaluation graph. In that case, a database or the like in which one or a plurality of graph feature amounts and diagnosis contents are associated with each other may be used.

望ましくは、前記超音波画像に対して関心領域を設定する関心領域設定手段を含み、前記二値化処理手段は前記関心領域内の部分画像に対して前記二値化処理を施して前記複数の二値化画像を生成する。この構成によれば、組織中の特定部分(特定組織部分)について性状診断を行える。例えば、肝臓内の腫瘍の内部に関心領域を設定することもできるし、その腫瘍を取り囲むように関心領域を設定することもできるし、更に、第1関心領域を腫瘍内に設定し、第2関心領域を腫瘍外に設定し、各関心領域ごとに作成されるスペックル評価グラフを目視的に又は自動的に比較するようにしてもよい。   Desirably, it includes a region-of-interest setting unit that sets a region of interest for the ultrasonic image, and the binarization processing unit performs the binarization processing on the partial image in the region of interest to perform the plurality of operations. A binarized image is generated. According to this configuration, a property diagnosis can be performed on a specific portion (specific tissue portion) in the tissue. For example, the region of interest can be set inside the tumor in the liver, the region of interest can be set to surround the tumor, the first region of interest is set in the tumor, and the second The region of interest may be set outside the tumor, and speckle evaluation graphs created for each region of interest may be compared visually or automatically.

望ましくは、前記超音波画像と共に前記スペックル評価グラフを表示する表示手段を含む。更に、各二値化画像を表示するようにしてもよい。   Desirably, a display means for displaying the speckle evaluation graph together with the ultrasonic image is included. Further, each binarized image may be displayed.

望ましくは、前記スペックル評価グラフを解析し、そのグラフの特徴量を表す特徴量データを出力するグラフ解析手段を含む。望ましくは、前記特徴量データには分散値、半値幅、ピークレベル、平均閾値レベル、などの1又は複数の情報が含まれる。特に複数の情報を入手して総合評価するのが望ましい。   Preferably, it includes a graph analysis means for analyzing the speckle evaluation graph and outputting feature amount data representing the feature amount of the graph. Preferably, the feature amount data includes one or more pieces of information such as a variance value, a half width, a peak level, and an average threshold level. In particular, it is desirable to obtain a plurality of information for comprehensive evaluation.

(2)望ましくは、超音波の送受波により得られたエコーデータに基づいて超音波画像を形成する画像形成手段と、組織の性状を診断するために、前記超音波画像を画像処理して当該超音波画像に含まれるスペックルの現れ方を評価するスペックル評価手段と、前記スペックル評価手段の評価結果を表示する表示手段と、が設けられる。 (2) Desirably, an image forming means for forming an ultrasonic image based on echo data obtained by transmitting and receiving ultrasonic waves, and image processing of the ultrasonic image in order to diagnose tissue properties Speckle evaluation means for evaluating the appearance of speckle included in the ultrasonic image and display means for displaying the evaluation result of the speckle evaluation means are provided.

上記構成によれば、組織の性状に応じてスペックルの現れ方(あるいはスペックルパターン)が異なる現象を利用して、組織の診断を行える。特に、自動的にスペックルの現れ方を分析評価できるので、その定量化から客観的な評価診断を行える。   According to the above configuration, diagnosis of a tissue can be performed using a phenomenon in which the appearance of speckles (or speckle patterns) varies depending on the properties of the tissue. In particular, since the speckle appearance can be automatically analyzed and evaluated, objective evaluation and diagnosis can be performed from the quantification.

望ましくは、前記スペックル評価手段はスペックルパターンの細かさを評価する。その細かさの評価に当たっては上記の閾値シフトに基づく二値化処理を適用するのが望ましいが、各種の画像処理技術を適用して、スペックルの細かさを評価することができる。望ましくは、前記組織は肝臓であるが、もちろん人体あるいは動物の他の臓器(実質臓器が望ましいが、血流などの液体であってもよい)に対して同様の手法を適用できる。   Preferably, the speckle evaluation means evaluates the fineness of the speckle pattern. In the evaluation of the fineness, it is desirable to apply the binarization processing based on the threshold shift described above, but the speckle fineness can be evaluated by applying various image processing techniques. Desirably, the tissue is the liver, but of course, the same technique can be applied to other organs of the human body or animals (which is preferably a real organ, but may be a fluid such as a bloodstream).

望ましくは、前記二値化処理に先立って、前記超音波画像に対してスペックルを強調する前処理を実行する前処理手段を含む。望ましくは、前記前処理手段は、スペックルに重畳しているベース成分を除去するフィルタである。スペックル以外の二次元的な濃淡差を解消、軽減すれば、スペックルの定量評価をより的確に行える。つまり、評価精度を高められる。   Preferably, pre-processing means for executing pre-processing for enhancing speckles on the ultrasonic image prior to the binarization processing is included. Preferably, the preprocessing means is a filter that removes a base component superimposed on speckle. Eliminating and reducing the two-dimensional shade difference other than speckle enables more accurate quantitative speckle evaluation. That is, the evaluation accuracy can be increased.

望ましくは、前記特徴量データにはグラフの一部の面積が含まれる。望ましくは、前記グラフに対して設定された基準値を基準とし、それを越える部分の面積が演算される。   Preferably, the feature amount data includes a partial area of the graph. Preferably, the reference value set for the graph is used as a reference, and the area exceeding the reference value is calculated.

(3)本発明は、超音波画像に対して、スペックルに重畳したベース成分を除去する前処理を実行するフィルタと、前記前処理後の超音波画像に対して、スペックルパターンを解析する解析手段と、を含む超音波診断装置に関する。 (3) The present invention analyzes a speckle pattern with respect to an ultrasonic image, a filter that executes preprocessing for removing a base component superimposed on speckles on the ultrasonic image, and the ultrasonic image after the preprocessing. And an ultrasonic diagnostic apparatus including the analysis means.

望ましくは、前記フィルタは、前記超音波画像上に設定されたウインド内において複数の画素値の平均値を演算する手段と、前記ウインド内の注目画素の画素値から前記平均値を減算して差分値を求める手段と、前記差分値にオフセット値を加算し、それを前記注目画像の新しい画素値とする手段と、を含む。   Preferably, the filter is configured to calculate an average value of a plurality of pixel values in a window set on the ultrasonic image, and to subtract the average value from a pixel value of a target pixel in the window. Means for obtaining a value, and means for adding an offset value to the difference value to make it a new pixel value of the image of interest.

(4)本発明は、超音波の送受波により得られたエコーデータに基づいて超音波画像を形成する画像形成手段と、前記超音波画像を元画像とし、それに対してスペックル評価を行うための前処理を実行する前処理手段と、前記前処理後の画像に対して、組織の性状を診断するために、スペックルの現れ方を評価するスペックル評価手段と、前記元画像と前記前処理後の画像とを表示する表示手段と、を含む超音波診断装置に関する。 (4) The present invention provides an image forming means for forming an ultrasonic image based on echo data obtained by transmission / reception of ultrasonic waves, and using the ultrasonic image as an original image for speckle evaluation. Pre-processing means for executing the pre-processing, speckle evaluation means for evaluating the appearance of speckles for diagnosing tissue properties with respect to the pre-processed image, the original image and the pre-processing The present invention relates to an ultrasonic diagnostic apparatus including display means for displaying a processed image.

上記構成によれば、前処理の前後の画像を併せて観察することができる。具体的に説明すると、例えば、上記のようなフィルタを用いて前処理を行った場合、スペックルの適正評価の点では有利となるが、元画像に現れていた構造などが逆に不明瞭となってしまう可能性がある。つまり、例えば画像の明るさが変動した結果、組織構造(例えば血管)の視認性が失われ、画像観察上の印象も変わる可能性がある。これに対し前処理後の画像と共に前処理前の元画像を表示すれば、前処理後の画像の表示によってスペックルの評価対象を特定できると共に、元画像の表示によって本来の組織構造をそのまま表示して診断に役立てることができる。   According to the above configuration, the images before and after the preprocessing can be observed together. More specifically, for example, when preprocessing is performed using the above-described filter, it is advantageous in terms of appropriate speckle evaluation, but the structure that appeared in the original image is unclear. There is a possibility of becoming. That is, for example, as a result of fluctuations in the brightness of the image, the visibility of the tissue structure (for example, blood vessels) is lost, and the impression on image observation may also change. On the other hand, if the pre-processed original image is displayed together with the pre-processed image, the speckle evaluation target can be specified by displaying the pre-processed image, and the original tissue structure can be displayed as it is by displaying the original image. Can be used for diagnosis.

望ましくは、前記表示手段には、更に前記スペックルの評価結果が表示される。このように、前処理前の元画像、前処理後の画像、スペックルの評価結果、二値化画像(例えば、代表となる1又は複数の二値化画像)、などを必要に応じて同時表示するのが望ましい。   Preferably, the speckle evaluation result is further displayed on the display means. In this way, an original image before preprocessing, an image after preprocessing, a speckle evaluation result, a binarized image (for example, one or a plurality of binarized images that are representative), and the like can be simultaneously used as necessary. It is desirable to display.

望ましくは、前記前処理手段は、前記超音波画像に対してスペックルに重畳したベース成分を除去するフィルタを含む。望ましくは、前記表示手段には、前記元画像と前記前処理後の画像とを識別するための情報が表示される。そのような情報は、文字、記号、色などを用いて表明できる。   Preferably, the preprocessing means includes a filter that removes a base component superimposed on speckles on the ultrasonic image. Preferably, information for identifying the original image and the pre-processed image is displayed on the display means. Such information can be expressed using letters, symbols, colors, and the like.

また、望ましくは、前記元画像に対して関心領域を設定する関心領域設定手段を含み、前記前処理手段は前記関心領域内の部分画像に対して前処理を実行し、前記表示手段には前記元画像と前記前処理後の関心領域内の部分画像とが表示される。望ましくは、前記関心領域設定手段により複数の関心領域を設定可能であり、前記表示手段には前記各関心領域を互いに識別するための情報が表示される。そのような情報は、文字、記号、色などを用いて表明できる。特に、各関心領域の枠に対して、それぞれ異なる色を付するようにしてもよい。   Preferably, the image processing apparatus includes a region of interest setting unit that sets a region of interest for the original image, the preprocessing unit performs preprocessing on the partial image in the region of interest, and the display unit includes An original image and a partial image in the region of interest after the preprocessing are displayed. Preferably, a plurality of regions of interest can be set by the region-of-interest setting unit, and information for identifying each region of interest is displayed on the display unit. Such information can be expressed using letters, symbols, colors, and the like. In particular, different frames may be assigned to the frames of each region of interest.

本発明は、超音波の送受波により得られたエコーデータに基づいて超音波画像を形成する画像形成手段と、前記超音波画像を元画像とし、それに対して関心領域を設定する関心領域設定手段と、前記元画像に対して前記関心領域を表す第1マーカーを合成してなる第1合成画像を作成する第1合成画像作成手段と、前記関心領域内の部分画像に対してスペックル評価を行うための前処理を実行する前処理手段と、前記前処理後の関心領域内の部分画像に対して前記第1マーカーに対応付けられた第2マーカーを合成してなる第2合成画像を作成する第2合成画像作成手段と、前記前処理後の関心領域内の部分画像に対して、組織の性状を診断するために、スペックルの現れ方を評価するスペックル評価手段と、前記第1合成画像と前記第2合成画像とを表示する表示手段と、を含む超音波診断装置に関する。ここで、対応関係にある第1マーカーと第2マーカーの各色を同じにして、その対応関係を視覚的に特定できるようにしてもよい。   The present invention relates to an image forming unit that forms an ultrasonic image based on echo data obtained by transmitting and receiving ultrasonic waves, and a region of interest setting unit that sets a region of interest for the ultrasonic image as an original image. First synthesized image creating means for creating a first synthesized image formed by synthesizing the first marker representing the region of interest with the original image, and speckle evaluation for the partial image in the region of interest. Pre-processing means for performing pre-processing for performing, and generating a second composite image formed by combining the second marker associated with the first marker with the partial image in the region of interest after the pre-processing A second synthesized image creating means, a speckle evaluating means for evaluating the appearance of speckles for diagnosing a tissue property with respect to the partial image in the region of interest after the preprocessing, and the first The composite image and the second combination Display means for displaying an image, an ultrasonic diagnostic apparatus including a. Here, the colors of the first marker and the second marker in the correspondence relationship may be the same so that the correspondence relationship can be visually identified.

以上説明したように、本発明によれば、スペックルあるいはスペックルパターンを診断に利用することが可能となる。   As described above, according to the present invention, speckles or speckle patterns can be used for diagnosis.

以下、本発明の好適な実施形態を図面に基づいて説明する。   DESCRIPTION OF EXEMPLARY EMBODIMENTS Hereinafter, preferred embodiments of the invention will be described with reference to the drawings.

図1〜図4には本実施形態の原理が概念的に示されている。まず、それらの図を用いてスペックル(スペックルパターン)の細かさの定量化方法について説明する。   1 to 4 conceptually show the principle of this embodiment. First, a method for quantifying the fineness of speckle (speckle pattern) will be described with reference to these drawings.

図1には、セクタ状の超音波画像100が示されている。この超音波画像100はいわゆるBモード画像であり、超音波ビームを電子セクタ走査することにより形成されるものである。もちろん、本発明は他の超音波画像に対しても適用することができる。   FIG. 1 shows a sector-shaped ultrasonic image 100. This ultrasonic image 100 is a so-called B-mode image, which is formed by scanning an ultrasonic beam with an electronic sector. Of course, the present invention can also be applied to other ultrasonic images.

超音波画像100に対して、ユーザーによって所望の形状及び所望のサイズをもった関心領域(ROI)104が設定される。例えば、超音波画像100として肝臓の断層面が表示され、その内部に腫瘍と思われる部分102が存在しているならば、その領域102内に関心領域104が設定される。もちろん、これは一例であって、複数の関心領域を設定し、それぞれについてスペックルパターンの評価を独立に行い、その評価結果を相互に対比してもよい。   A region of interest (ROI) 104 having a desired shape and a desired size is set by the user for the ultrasonic image 100. For example, if a tomographic plane of the liver is displayed as the ultrasound image 100 and a portion 102 that seems to be a tumor exists in the ultrasound image 100, a region of interest 104 is set in the region 102. Of course, this is merely an example, and a plurality of regions of interest may be set, speckle patterns may be evaluated independently for each, and the evaluation results may be compared with each other.

図2には、関心領域内において抽出された画像とその画像に対して閾値を用いた二値化処理を適用した後の画像の例が示されている。   FIG. 2 shows an example of an image extracted in the region of interest and an image after binarization processing using a threshold is applied to the image.

(A)には関心領域内の画像が示されており、その内部には上述したスペックルが含まれている。本実施形態においては、この(A)に示す画像に対して、閾値レベルを0から最大値までに段階的にスキャンさせながら、各閾値を利用して画像の二値化処理が実行される。図2においては、閾値が10に設定された場合の二値化画像が(B)で示されており、閾値が20に設定された場合の二値化画像が(C)で示され、閾値に30が設定された場合の二値化画像が(D)で示されている。例えば、画像の輝度が0から63までの範囲内にある場合には、閾値は0から63までスキャンされ、図2に示すような各閾値ごとの二値化画像が生成される。もちろん、スペックルの現れ方によっては、その閾値のスキャン幅を一定範囲内に制限し、あるいはユーザーによってそのスキャン範囲を任意に選択できるようにしてもよい。   In (A), an image in the region of interest is shown, and the above-mentioned speckle is included therein. In the present embodiment, image binarization processing is executed using each threshold value while scanning the image shown in (A) stepwise from 0 to the maximum threshold level. In FIG. 2, the binarized image when the threshold is set to 10 is shown by (B), and the binarized image when the threshold is set to 20 is shown by (C). (D) shows a binarized image when 30 is set in. For example, when the luminance of the image is in the range from 0 to 63, the threshold is scanned from 0 to 63, and a binarized image for each threshold as shown in FIG. 2 is generated. Of course, depending on how the speckle appears, the scan width of the threshold value may be limited within a certain range, or the scan range may be arbitrarily selected by the user.

各閾値ごとに二値化画像が生成されると、本実施形態においては、高輝度、すなわち画素値1(あるいは画素値0)を有する独立領域(島)の個数が計算される。これは公知のラベリング処理を適用することにより実行され、すなわち各二値化画像ごとに高輝度の領域がいくつあるかが分析される。   When a binarized image is generated for each threshold, in the present embodiment, the number of independent areas (islands) having high luminance, that is, pixel value 1 (or pixel value 0) is calculated. This is performed by applying a known labeling process, i.e., how many high-brightness regions are present for each binarized image.

そして、図3に示すように、横軸に閾値をとり、縦軸に島の個数をとったヒストグラムが作成される。ここで、各横軸は閾値に対応し、それらの各閾値ごとに上記二値化画像が生成され、その二値化処理画像ごとに島の数をカウントした結果がラベル数としてグラフにプロットされる。すると、符号200で示されるようにスペックル評価グラフがヒストグラムとして生成されることになる。   Then, as shown in FIG. 3, a histogram is created with the horizontal axis representing the threshold value and the vertical axis representing the number of islands. Here, each horizontal axis corresponds to a threshold value, and the above binarized image is generated for each threshold value, and the result of counting the number of islands for each binarized image is plotted on the graph as the number of labels. The Then, as indicated by reference numeral 200, a speckle evaluation graph is generated as a histogram.

上述したように、組織の性状に応じてスペックルパターンの細かさあるいはその現れ方は相違しており、したがって、スペックル評価グラフ200の形状は組織の性状を表すものとなる。したがって、そのグラフ自体を画像表示すれば、組織性状の診断に役立てることができ、さらに、客観的な定量化のためには、スペックル評価グラフ200を何らかの特徴量により数値化することが望まれる。   As described above, the fineness of the speckle pattern or how it appears depends on the nature of the tissue. Therefore, the shape of the speckle evaluation graph 200 represents the nature of the tissue. Therefore, if the graph itself is displayed as an image, it can be used for diagnosis of tissue properties. Further, for objective quantification, it is desirable to quantify the speckle evaluation graph 200 by some feature amount. .

この場合においては、例えば分散値、平均閾値、半値幅、最大ラベル数などの各種の特徴量データが演算され、望ましくは、それらの複数の特徴量データを組み合わせて性状が評価される。   In this case, for example, various feature amount data such as a variance value, an average threshold value, a full width at half maximum, and the maximum number of labels are calculated. Preferably, the property is evaluated by combining the plurality of feature amount data.

図4にはある線上のスペックルパターン(右側)とスペックル評価グラフ(左側)との関係が示されており、(A)に示されるように、スペックルにおける山と谷の輝度差がバラバラで、それらの差が大きいような場合には、符号200Aで示されるようにスペックル評価グラフは大きな分散値をもったなだらかな山状の分布となる。なお、(A)及び(B)の右側に示されるグラフの横軸は座標であり、縦軸は画素の輝度に相当している。   FIG. 4 shows the relationship between the speckle pattern (on the right side) and the speckle evaluation graph (on the left side) on a certain line. As shown in (A), the luminance difference between the peaks and valleys in the speckles varies. When the difference between them is large, the speckle evaluation graph has a gentle mountain distribution with a large dispersion value as indicated by reference numeral 200A. Note that the horizontal axis of the graphs shown on the right side of (A) and (B) is the coordinate, and the vertical axis corresponds to the luminance of the pixel.

一方、(B)で示されるように、スペックルにおける山と谷の輝度値がある程度そろっていれば、すなわち山と谷の差が小さければ、符号200Bで示されるように分散値の小さい急峻なあるいは先鋭な山状の分布となる。   On the other hand, as shown in (B), if the luminance values of the peaks and valleys in the speckle are uniform to some extent, that is, if the difference between the peaks and valleys is small, a steep and small dispersion value is obtained as indicated by reference numeral 200B. Or it becomes a sharp mountain-shaped distribution.

いずれにしても、組織の性状によってスペックルパターンが異なり、そのスペックルパターンが異なるとスペックル評価グラフの形態が異なるため、その形態評価をもってスペックルの定量化を実現することが可能となる。すなわち組織性状の評価を客観的に行える。   In any case, the speckle pattern differs depending on the nature of the tissue, and if the speckle pattern is different, the form of the speckle evaluation graph is different. Thus, speckle quantification can be realized by the form evaluation. That is, it is possible to objectively evaluate the tissue properties.

図5には、そのような機能をもった超音波診断装置の実施形態がブロック図として示されている。   FIG. 5 is a block diagram showing an embodiment of an ultrasonic diagnostic apparatus having such a function.

探触子10は生体内に挿入され、あるいは生体表面上に当接して用いられる超音波探触子(プローブ)である。この探触子10は本実施形態において複数の振動素子からなるアレイ振動子を有しており、そのアレイ振動子によって超音波ビームが形成される。ここで、その超音波ビームの走査方式としては、電子セクタ走査や電子リニア走査などをあげることができ、また、超音波ビームを二次元走査し、三次元データ取込領域を形成するような場合にも本発明を適用できる。この場合においては、スペックルの評価は三次元空間内において行われることになる。   The probe 10 is an ultrasonic probe (probe) that is inserted into a living body or used in contact with the surface of the living body. In this embodiment, the probe 10 has an array transducer including a plurality of transducer elements, and an ultrasonic beam is formed by the array transducer. Here, examples of the scanning method of the ultrasonic beam include electronic sector scanning and electronic linear scanning. In addition, when the ultrasonic beam is two-dimensionally scanned to form a three-dimensional data capturing area. The present invention can also be applied to. In this case, speckle evaluation is performed in a three-dimensional space.

送信器12は送信ビームフォーマーとして機能し、探触子10に対して送信信号を供給する。受信器14は受信ビームフォーマーとして機能し、探触子10から出力される受信信号に対して整相加算処理を適用し、これによって整相加算後の受信信号を出力する。   The transmitter 12 functions as a transmission beam former and supplies a transmission signal to the probe 10. The receiver 14 functions as a reception beam former, applies a phasing addition process to the reception signal output from the probe 10, and thereby outputs a reception signal after phasing addition.

DSC(デジタルスキャンコンバータ)18は、いくつかの画像処理機能を有しているものであるが、そのフレームメモリ20上には本実施形態においてBモード画像を構成する画像データが格納される。すなわち、受信器14から出力される受信信号(エコーデータ)に対して例えば対数圧縮処理などを行い、さらに座標変換などを適用して二次元断層画像が形成される。   The DSC (digital scan converter) 18 has several image processing functions, but the frame memory 20 stores image data constituting a B-mode image in this embodiment. That is, for example, a logarithmic compression process or the like is performed on the reception signal (echo data) output from the receiver 14, and further a coordinate transformation is applied to form a two-dimensional tomographic image.

本実施形態において、領域設定器22は、図1に示した関心領域104の位置及びサイズをユーザー指定するための入力手段である。   In the present embodiment, the region setting unit 22 is an input unit for designating the position and size of the region of interest 104 shown in FIG.

フレームメモリ20からはBモード画像全体の画像データが出力され、具体的には、その画像データが画像合成部32に出力される。また、フレームメモリ20からは領域設定器22によって設定された関心領域内の画像データが切り出されてそれが二値化処理部24へ出力される。二値化処理部24は、閾値スキャン回路26によってスキャン設定される閾値を基準として関心領域内の各画素値を二値化処理し、すなわち二値化画像を形成する。この場合において、閾値スキャン回路26は、0から最大値まで閾値をスキャンしており、各閾値ごとに二値化画像が生成されることになる。その様子が図2に示されている。   Image data of the entire B mode image is output from the frame memory 20, and specifically, the image data is output to the image composition unit 32. Further, the image data in the region of interest set by the region setting unit 22 is cut out from the frame memory 20 and output to the binarization processing unit 24. The binarization processing unit 24 binarizes each pixel value in the region of interest with reference to the threshold set by the threshold scan circuit 26, that is, forms a binarized image. In this case, the threshold scan circuit 26 scans the threshold from 0 to the maximum value, and a binarized image is generated for each threshold. This is shown in FIG.

ラベリング処理部28は、上述したように、各二値化画像ごとに白又は黒の独立領域(島)の個数を演算するラベリング処理を実行する回路である。そしてグラフ作成部30は図3に示したようなスペックル評価グラフをヒストグラムとして作成する回路である。   As described above, the labeling processing unit 28 is a circuit that executes a labeling process for calculating the number of white or black independent regions (islands) for each binarized image. The graph creation unit 30 is a circuit that creates a speckle evaluation graph as shown in FIG. 3 as a histogram.

画像合成部32は、表示画像を合成する手段であり、その表示画像内には二次元断層画像とスペックル評価グラフとが含まれる。もちろん他の画像を併せて合成するようにしてもよい。   The image synthesis unit 32 is a means for synthesizing a display image, and the display image includes a two-dimensional tomographic image and a speckle evaluation graph. Of course, other images may be combined together.

符号36〜符号48で示される各回路はスペックル評価グラフを定量化するためのものである。まず、サンプル数検出部36は、図3に示したスペックル評価グラフ200における各閾値ごとのラベル数を積算し、すなわちグラフの面積を演算する回路である。二乗平均値演算部40は、上記の面積の二乗を演算する。平均閾値演算部42は、以下に示す式に基づいて平均しきい値を演算する回路であり、分散値演算部44は以下に示す式に基づいて分散値を演算する回路である。   Each circuit indicated by reference numerals 36 to 48 is for quantifying the speckle evaluation graph. First, the sample number detection unit 36 is a circuit that integrates the number of labels for each threshold in the speckle evaluation graph 200 shown in FIG. 3, that is, calculates the area of the graph. The mean square value calculation unit 40 calculates the square of the area. The average threshold value calculation unit 42 is a circuit that calculates an average threshold value based on the following formula, and the variance value calculation unit 44 is a circuit that calculates a variance value based on the following formula.

Figure 2006212445
Figure 2006212445

また、半値幅演算部46は、図3に示すスペックル評価グラフ200の半値幅を演算する回路であり、最大ラベル数検出部48はスペックル評価グラフ200のピークを検出する回路である。それらのグラフ評価結果はそれぞれ表示部34に出力されており、必要に応じて数値が表示される。   The half-value width calculation unit 46 is a circuit that calculates the half-value width of the speckle evaluation graph 200 shown in FIG. 3, and the maximum label number detection unit 48 is a circuit that detects the peak of the speckle evaluation graph 200. The graph evaluation results are output to the display unit 34, and numerical values are displayed as necessary.

もちろん、1又は複数のグラフ特徴量に基づいて組織自体の評価を自動的に行うようにしてもよく、その場合にはデータベースなどに診断情報を格納し、各特徴量の大きさあるいは値の組み合わせから一定の疾病情報が画像表示されるようにしてもよい。   Of course, the organization itself may be automatically evaluated based on one or a plurality of graph feature amounts. In that case, diagnostic information is stored in a database or the like, and the size or value combination of each feature amount is stored. Thus, certain disease information may be displayed as an image.

ちなみに、上記実施形態においては、二値化画像に対するラベリング処理によって結果として図3に示したスペックル評価グラフ200を生成し、これに基づいて組織性状の評価を行うようにしたが、スペックルあるいはスペックルパターンの評価にあたっては、それが組織性状に相関付けられていることを前提として他の画像処理を適用するようにしてもよい。例えば空間周波数の分析や所定のフィルタ演算などを適用するようにしてもよい。いずれにしても、スペックルパターンを画像処理によって分析し、それを診断に利用することにより従来の超音波画像上では定量化し得なかった情報を取得して疾病診断に役立てることができる。   Incidentally, in the above-described embodiment, the speckle evaluation graph 200 shown in FIG. 3 is generated as a result of the labeling process for the binarized image, and the tissue property is evaluated based on the speckle evaluation graph 200. When the speckle pattern is evaluated, other image processing may be applied on the assumption that the speckle pattern is correlated with the tissue property. For example, spatial frequency analysis or predetermined filter calculation may be applied. In any case, by analyzing the speckle pattern by image processing and using it for diagnosis, it is possible to acquire information that could not be quantified on a conventional ultrasonic image and use it for disease diagnosis.

図6には、他の実施形態に係る超音波診断装置の全体構成がブロック図として示されている。なお、図5に示す構成と同様の構成には同一符号を付しその説明を省略する。   FIG. 6 is a block diagram showing the overall configuration of an ultrasonic diagnostic apparatus according to another embodiment. In addition, the same code | symbol is attached | subjected to the structure similar to the structure shown in FIG. 5, and the description is abbreviate | omitted.

図6に示す構成においては、DSC18と二値化処理部24との間に空間フィルタ50が設けられている。また、グラフ作成部30の後段に部分面積演算部52が設けられている。以下、それらについて詳述する。   In the configuration shown in FIG. 6, a spatial filter 50 is provided between the DSC 18 and the binarization processing unit 24. In addition, a partial area calculation unit 52 is provided after the graph creation unit 30. These will be described in detail below.

空間フィルタ50はスペックルを強調するための前処理を遂行する回路である。具体的には、超音波画像の全体あるいはその内部に設定された関心領域に対して図7に示すようなウインド54を二次元的にスキャン(ラスター走査)し、各スキャン位置において注目画素56について新しい画素値を求めるものである。図7において符号53は超音波画像を示しており、符号56は注目画素を示している。この注目画素56はウインド54の中央に設定されている。ウインド54のサイズは図7に示す例において11×11画素である。もちろん、ウインド54のサイズとしては他のものを採用できる。   The spatial filter 50 is a circuit that performs preprocessing for enhancing speckles. Specifically, the window 54 as shown in FIG. 7 is scanned two-dimensionally (raster scanning) with respect to the entire ultrasound image or the region of interest set in the ultrasound image, and the pixel of interest 56 at each scan position. A new pixel value is obtained. In FIG. 7, reference numeral 53 indicates an ultrasonic image, and reference numeral 56 indicates a target pixel. The target pixel 56 is set at the center of the window 54. The size of the window 54 is 11 × 11 pixels in the example shown in FIG. Of course, other sizes of the window 54 can be used.

図8には、空間フィルタ50における処理内容がフローチャートとして示されている。あるウインド位置において、まずS101では平均値が演算される。具体的にはウインド内に存在する複数の画素値について平均値が演算される。S102では、注目画素の画素値から上記で求められた平均値が減算され、これにより差分値が求められる。S103では、その差分値に対してオフセット値(例えば32)が加算され、その結果が注目画素についての更新画素値とされる。そして、S104では、次のウインド位置があるか否かが判断され、次のウインド位置がある場合にはウインドの位置をシフトさせて上記同様の処理が繰り返し実行される。   FIG. 8 shows a processing content in the spatial filter 50 as a flowchart. At a certain window position, an average value is first calculated in S101. Specifically, an average value is calculated for a plurality of pixel values existing in the window. In S102, the average value obtained above is subtracted from the pixel value of the target pixel, thereby obtaining a difference value. In S103, an offset value (for example, 32) is added to the difference value, and the result is used as an updated pixel value for the target pixel. In S104, it is determined whether or not there is a next window position. If there is a next window position, the window position is shifted and the same processing as described above is repeatedly executed.

その結果、図9に示されるように、超音波画像に対して前処理が実行されることになる。具体的には、図9(A)にはフィルタリング処理前の様子が示されており、図9(B)にはフィルタリング処理後の様子が示されている。なお、画像は二次元的に存在しているが、図9においては説明のためX方向についての輝度値の分布すなわち輝度値グラフが示されている。   As a result, as shown in FIG. 9, preprocessing is performed on the ultrasound image. Specifically, FIG. 9A shows a state before the filtering process, and FIG. 9B shows a state after the filtering process. Although the image exists two-dimensionally, FIG. 9 shows a luminance value distribution in the X direction, that is, a luminance value graph, for explanation.

(A)に示されるように、例えばX方向に沿ってベース成分57が単調減少している場合においてそのベース成分57上にスペックルパターンが繰り返し存在している場合、そのままその画像に対して処理を行ってしまうと、ラベリング処理において的確なスペックルの解析を行うことが困難となる。これに対し、(B)に示すように、ベース成分57を除去し、そしてオフセット値58を加算する処理を遂行すれば、ベース成分57による寄与分を排除してスペックルパターンによる変動成分を抽出することが可能となる。すなわちスペックル強調結果を得ることが可能となる。   As shown in (A), for example, when the base component 57 monotonously decreases along the X direction, if a speckle pattern is repeatedly present on the base component 57, the image is processed as it is. If this is done, it will be difficult to perform accurate speckle analysis in the labeling process. On the other hand, as shown in (B), if the process of removing the base component 57 and adding the offset value 58 is performed, the contribution due to the base component 57 is eliminated and the fluctuation component due to the speckle pattern is extracted. It becomes possible to do. That is, a speckle enhancement result can be obtained.

以上のような前処理を経た画像データが図6において二値化処理部24に供給され、上記同様の二値化処理が遂行されることになる。   The image data that has undergone the preprocessing as described above is supplied to the binarization processing unit 24 in FIG. 6, and the same binarization processing as described above is performed.

図10は図6に示した空間フィルタ50の具体的な回路構成例を示すものである。符号60は11個からなるラインメモリ62を示している。それらのラインメモリ62は互いに直列的に接続されており、それらの全体の入力及び出力と各ライン62間からデータが一画素分ずつ出力される。それらの出力は加算器64にて加算され、その加算値はラッチ65においていったん記憶された後に出力される。その出力値は演算器70のA入力端子に供給される他、タイミング調整用のディレイライン66に入力される。ここで、ディレイライン66は11個のラッチ67から構成されるものであり、11ステップ分だけ加算値の遅延が行われる。そして、そのディレイライン66の出力が演算器70のB入力端子に入力される。演算器70はA入力端子に入力された値からB入力端子に入力された値を減算する回路であり、その減算結果は加算器71に入力され、その加算器71においてはその後段に設けられたラッチ72からの出力が加算され、そして、その加算結果がラッチ72に入力されている。   FIG. 10 shows a specific circuit configuration example of the spatial filter 50 shown in FIG. Reference numeral 60 denotes an eleven line memory 62. These line memories 62 are connected in series with each other, and data is output pixel by pixel from the entire input and output and between each line 62. These outputs are added by the adder 64, and the added value is temporarily stored in the latch 65 and then output. The output value is supplied to the A input terminal of the arithmetic unit 70 and also input to the delay line 66 for timing adjustment. Here, the delay line 66 is composed of 11 latches 67, and the added value is delayed by 11 steps. Then, the output of the delay line 66 is input to the B input terminal of the computing unit 70. The arithmetic unit 70 is a circuit that subtracts the value input to the B input terminal from the value input to the A input terminal. The subtraction result is input to the adder 71, and the adder 71 is provided in the subsequent stage. The outputs from the latch 72 are added, and the addition result is input to the latch 72.

したがって、この構成により、11×11画素のウインドについての全画素値が加算されることになり、それは各ウインド位置において求められることになる。そして、平均化回路73において、その加算値が全画素数によって割られ、その結果平均値が求められる。その平均値は演算器74のB入力端子に入力される。一方、注目画素の画素値は中央のラインメモリから出力され、その画素値はディレイライン68に入力される。このディレイライン68はタイミング調整用の回路であり、6個のラッチ69から構成されている。ディレイライン68の出力は演算器74のA入力端子に供給されている。   Therefore, with this configuration, all pixel values for a window of 11 × 11 pixels are added, which is obtained at each window position. Then, in the averaging circuit 73, the added value is divided by the total number of pixels, and as a result, an average value is obtained. The average value is input to the B input terminal of the computing unit 74. On the other hand, the pixel value of the target pixel is output from the central line memory, and the pixel value is input to the delay line 68. This delay line 68 is a circuit for timing adjustment, and is composed of six latches 69. The output of the delay line 68 is supplied to the A input terminal of the arithmetic unit 74.

演算器74では、A入力端子に入力された値からB入力端子に入力された値を減算し、その減算結果を加算器75へ出力している。すなわち注目画素から平均値が減算され、それにより得られた差分値が加算器75においてオフセット値と加算されている。そしてその加算値はラッチ76を介して出力されている。   The computing unit 74 subtracts the value input to the B input terminal from the value input to the A input terminal, and outputs the subtraction result to the adder 75. That is, the average value is subtracted from the target pixel, and the difference value obtained thereby is added to the offset value in the adder 75. The added value is output via a latch 76.

もちろん、図10に示した構成は一例であって空間フィルタとしては他の回路構成を採用しうる。またハードウエア構成によらずにソフトウエアによって空間フィルタリングを実現するようにしてもよい。これは他の回路構成についても同様である。   Of course, the configuration shown in FIG. 10 is an example, and other circuit configurations may be adopted as the spatial filter. Further, spatial filtering may be realized by software regardless of the hardware configuration. The same applies to other circuit configurations.

上述したように、図6に示す構成においては部分面積演算部52が設けられている。これについて図11を用いて説明する。   As described above, the partial area calculation unit 52 is provided in the configuration shown in FIG. This will be described with reference to FIG.

図11には、スペックルの解析により得られたスペックル評価グラフ200が示されている。このスペックル評価グラフ200に対してユーザー設定によりあるいは自動的に求められた基準値80を基準として部分面積82が演算される。具体的にはその基準値80を超える部分の面積が部分面積82として演算されている。この部分面積82は、基準値80を超える各閾値に対応するラベル数の総和を意味しており、スペックル評価グラフ200の1つの特徴量を示すものである。したがって、そのような部分面積82は表示部34に出力され必要に応じて数値表示され、あるいはグラフ上において色づけなどによって表現される。もちろん、上述したように複数の特徴量を組み合わせてスペックルパターンを評価する場合においても上記の部分面積を用いるようにしてもよい。   FIG. 11 shows a speckle evaluation graph 200 obtained by speckle analysis. A partial area 82 is calculated based on a reference value 80 obtained by user setting or automatically for the speckle evaluation graph 200. Specifically, the area of the portion exceeding the reference value 80 is calculated as the partial area 82. This partial area 82 means the total number of labels corresponding to each threshold value exceeding the reference value 80, and indicates one feature amount of the speckle evaluation graph 200. Accordingly, such a partial area 82 is output to the display unit 34 and numerically displayed as necessary, or expressed by coloring or the like on the graph. Of course, as described above, the partial area may be used even when a speckle pattern is evaluated by combining a plurality of feature amounts.

上記の実施形態においては、画像上に存在するスペックルとは関係のない二次元的な濃淡差を積極的に除去してスペックルを強調し、その結果スペックルパターンを高精度に解析することが可能となる。また、その解析においては部分面積を用いることができ、スペックルパターンの他面的な分析を実現することが可能となる。なお、上述した空間フィルタ50は上記のスペックルパターンの解析以外の超音波画像の画像処理においても用いることが可能である。   In the above-described embodiment, the speckle is emphasized by positively removing the two-dimensional shading difference unrelated to the speckle existing on the image, and as a result, the speckle pattern is analyzed with high accuracy. Is possible. In addition, a partial area can be used in the analysis, and it is possible to realize another side analysis of the speckle pattern. The spatial filter 50 described above can also be used in image processing of an ultrasonic image other than the analysis of the speckle pattern.

図12には、更に他の実施形態に係る超音波診断装置の構成がブロック図として示されている。   FIG. 12 is a block diagram showing the configuration of an ultrasonic diagnostic apparatus according to still another embodiment.

図12において、図5及び図6に示した構成と同様の構成には同一符号を付し、その説明を省略する。図6においては、DSC18から1つのラインが引き出され、そのラインは超音波画像全体のデータと関心領域内の部分画像のデータとを表していたが、図12においては、技術的理解を助けるために、それらの2つのデータを2つのラインによって表した。すなわち、DSC18から出力された元画像としての超音波画像全体のデータ300は、画像合成部306へ出力されており、また、そのデータ300は、必要に応じて、空間フィルタ50に出力される。DSC18から出力された関心領域内の部分画像のデータ302は、前処理を実行する空間フィルタ50へ出力されている。   In FIG. 12, the same components as those shown in FIGS. 5 and 6 are denoted by the same reference numerals, and the description thereof is omitted. In FIG. 6, one line is drawn from the DSC 18, and the line represents the data of the entire ultrasound image and the data of the partial image in the region of interest, but in FIG. 12, to help technical understanding. The two data are represented by two lines. That is, the entire ultrasonic image data 300 as the original image output from the DSC 18 is output to the image synthesis unit 306, and the data 300 is output to the spatial filter 50 as necessary. The partial image data 302 in the region of interest output from the DSC 18 is output to the spatial filter 50 that executes preprocessing.

一方、領域設定器22によって設定された1又は複数の関心領域についての座標データは、DSC18に出力されると共に、画像合成部306へ出力されている。また、空間フィルタ50による前処理後の画像データ304が二値化処理部24へ出力され、同時に、画像合成部306にも出力されている。   On the other hand, the coordinate data for one or a plurality of regions of interest set by the region setting unit 22 is output to the DSC 18 and also to the image composition unit 306. Further, the image data 304 after the preprocessing by the spatial filter 50 is output to the binarization processing unit 24 and simultaneously output to the image composition unit 306.

以上のように、画像合成部306には、グラフ作成部30が作成したグラフのデータの他に、元画像のデータ300及び前処理後の関心領域内の部分画像のデータ304が入力されている。その画像合成部は306は、イメージ合成機能を有しており、特に、超音波画像とマーカーなどのグラフィック画像とを合成する機能を有している。画像合成部306が有する具体的な画像処理機能について図13及び図14に示される表示例を用いて説明する。   As described above, in addition to the graph data created by the graph creation unit 30, the original image data 300 and the pre-processed partial image data 304 in the region of interest are input to the image composition unit 306. . The image synthesizing unit 306 has an image synthesizing function, and particularly has a function of synthesizing an ultrasonic image and a graphic image such as a marker. A specific image processing function of the image composition unit 306 will be described with reference to display examples shown in FIGS.

図13には、表示部34における表示例が示されている。表示画面308上には、空間フィルタ50による前処理を行う前の超音波画像つまり元画像310が表示され、それと並んで、空間フィルタ50による前処理を行った後の前処理後画像312が表示される。このような表示がなされる場合には、空間フィルタ50には、超音波画像全体のデータ300が入力される。つまり、スペックルパターンの評価対象は画像全体となる。   FIG. 13 shows a display example on the display unit 34. On the display screen 308, the ultrasonic image before the pre-processing by the spatial filter 50, that is, the original image 310 is displayed, and the pre-processed image 312 after the pre-processing by the spatial filter 50 is displayed along with the ultrasonic image. Is done. When such a display is performed, the entire ultrasonic image data 300 is input to the spatial filter 50. In other words, the speckle pattern evaluation target is the entire image.

元画像310と前処理後画像312の関係を視覚的に認識可能とするために、表示画面308上には両画像310,312間に矢印マーク314が表示される。この矢印マーク314の向きを観察することによって前処理前後の各画像310,312を確認することができる。もちろん、そのような画像識別情報は、矢印マーク314に限られず、文字や他の記号であってもよい。表示画面308上には、必要に応じて、スペックルパターンの評価結果として例えばスペックル評価グラフが表示され、また、そのグラフの解析結果が表示される。   In order to visually recognize the relationship between the original image 310 and the preprocessed image 312, an arrow mark 314 is displayed between the images 310 and 312 on the display screen 308. By observing the direction of the arrow mark 314, the images 310 and 312 before and after the preprocessing can be confirmed. Of course, such image identification information is not limited to the arrow mark 314 but may be a character or another symbol. On the display screen 308, for example, a speckle evaluation graph is displayed as the speckle pattern evaluation result, and the analysis result of the graph is displayed as necessary.

この表示例によれば、特に、元画像が併せて表示されるので組織構造を明確に認識でき、元画像との対比において前処理画像を観察することができる。その上で、組織の総合的な診断を行うことができる。   According to this display example, in particular, since the original image is displayed together, the tissue structure can be clearly recognized, and the preprocessed image can be observed in comparison with the original image. In addition, a comprehensive diagnosis of the tissue can be performed.

以上のように、画像合成部306は、複数の画像やグラフなどを含む表示画像を生成している。表示画面上に、必要に応じて、1又は複数の二値化画像を併せて表示するようにしてもよい。   As described above, the image composition unit 306 generates a display image including a plurality of images and graphs. If necessary, one or a plurality of binarized images may be displayed on the display screen.

図14に示す他の表示例においては、表示画面320に元画像表示領域400及び前処理後画像表示領域402が設定されている。元画像表示領域400には、元画像322が表示される。その元画像上においては、ユーザー操作によって、1又は複数の関心領域を設定可能である。図14に示す例では、2つの関心領域324,326が設定されている。それらを視覚的に特定するために、本実施形態においては、2つの関心領域324,326を特定する矩形のマーカー(第1マーカー)324A,326Aが合成表示される。これらは各関心領域の枠に相当する。各マーカー324A,326Aの色は互いに異なり、マーカー324Aの色は例えば赤であり、マーカー326Aの色は例えば青である。もちろん、着色以外の手法によって各関心領域が区別されるようにしてもよい。例えば、線種の変更、輝度の変更、文字や記号の挿入、などの手法を用いることができる。   In another display example shown in FIG. 14, an original image display area 400 and a pre-processed image display area 402 are set on the display screen 320. An original image 322 is displayed in the original image display area 400. On the original image, one or a plurality of regions of interest can be set by a user operation. In the example shown in FIG. 14, two regions of interest 324 and 326 are set. In order to visually identify them, in the present embodiment, rectangular markers (first markers) 324A and 326A that identify the two regions of interest 324 and 326 are combined and displayed. These correspond to the frame of each region of interest. The colors of the markers 324A and 326A are different from each other, the color of the marker 324A is, for example, red, and the color of the marker 326A is, for example, blue. Of course, each region of interest may be distinguished by a technique other than coloring. For example, methods such as line type change, luminance change, and insertion of characters and symbols can be used.

前処理後画像表示領域402には、図示の例において、2つの関心領域に対応した前処理後の2つの部分画像328,330が表示される。ここで、各部分画像328,330ごとに矩形のマーカー(第2マーカー)328A,330Aが合成表示され、前者の色は例えば赤であり、後者の色は例えば青である。つまり、関心領域と部分画像との対応関係が各マーカーに同じ着色を施すことによって一目瞭然となる。この場合に、上記同様に、着色以外の手法を利用して各関心領域と前処理後の各部分画像との対応関係を表示するようにしてもよい。   In the pre-processed image display area 402, two partial images 328 and 330 after the pre-process corresponding to the two regions of interest are displayed in the illustrated example. Here, rectangular markers (second markers) 328A and 330A are combined and displayed for each of the partial images 328 and 330, and the former color is, for example, red, and the latter color is, for example, blue. That is, the correspondence between the region of interest and the partial image becomes clear at a glance by applying the same coloring to each marker. In this case, as described above, the correspondence between each region of interest and each pre-processed partial image may be displayed using a method other than coloring.

画像合成部306は、以上のような画像とマーカーを合成する機能と、各画像を含む表示画面を構成する機能とを有している。もちろん、表示画面320上には、必要に応じて、スペックル評価グラフ、グラフの解析結果、二値化画像などが表示される。   The image composition unit 306 has a function of composing the image and the marker as described above, and a function of configuring a display screen including each image. Of course, a speckle evaluation graph, a graph analysis result, a binarized image, and the like are displayed on the display screen 320 as necessary.

以上のように、図14に示す表示例によれば、元画像上において関心領域の位置及びサイズを確認でき、しかも元画像上において組織の構造を明確に把握でき、それとの対比において、スペックル評価対象となっている前処理後の部分画像を確認できる。そして、スペックル評価グラフ及びその解析結果を併せて考慮することにより、総合的な組織診断が可能となる。   As described above, according to the display example shown in FIG. 14, the position and size of the region of interest can be confirmed on the original image, and the structure of the tissue can be clearly understood on the original image. The pre-processed partial image that is the object of evaluation can be confirmed. Then, by considering the speckle evaluation graph and the analysis result together, a comprehensive tissue diagnosis can be performed.

なお、ユーザーが各関心領域を設定する場合には、各関心領域ごとに、元画像上に関心領域を示す一定サイズの初期マーカーを自動的に表示させ、それをマニュアルで移動し、あるいはそのサイズをマニュアルで変更することによって、所望の位置に所望のサイズをもった関心領域を設定させるようにすればよい。その場合において、各関心領域ごとにマーカーの色をユーザー指定させるようにしてもよいし、自動的に各マーカーの色が決定されるようにしてもよい。なお、表示部34は、実際には表示処理部及び表示器によって構成されるが、表示器は、単一のディスプレイによって構成されてもよいし、複数のディスプレイ(例えば主ディスプレイ及び補助ディスプレイ)によって構成されてもよい。   In addition, when the user sets each region of interest, for each region of interest, an initial marker of a certain size indicating the region of interest is automatically displayed on the original image, and it is moved manually or its size Is changed manually to set a region of interest having a desired size at a desired position. In that case, the color of the marker may be designated by the user for each region of interest, or the color of each marker may be automatically determined. The display unit 34 is actually configured by a display processing unit and a display unit, but the display unit may be configured by a single display or a plurality of displays (for example, a main display and an auxiliary display). It may be configured.

超音波画像と関心領域との関係を示す図である。It is a figure which shows the relationship between an ultrasonic image and a region of interest. 関心領域内の画像とそれを二値化処理した後の二値化画像を示す図である。It is a figure which shows the image in a region of interest, and the binarized image after binarizing it. スペックル評価グラフの一例を示す図である。It is a figure which shows an example of a speckle evaluation graph. スペックル評価グラフとスペックルパターンとの関係を示す図である。It is a figure which shows the relationship between a speckle evaluation graph and a speckle pattern. 本実施形態に係る超音波診断装置の全体構成を示すブロック図である。1 is a block diagram showing an overall configuration of an ultrasonic diagnostic apparatus according to the present embodiment. 他の実施形態に係る超音波診断装置の全体構成を示すブロック図である。It is a block diagram which shows the whole structure of the ultrasound diagnosing device which concerns on other embodiment. 空間フィルタが有するウインドを説明するための図である。It is a figure for demonstrating the window which a spatial filter has. 空間フィルタにおける処理内容を説明するためのフローチャートである。It is a flowchart for demonstrating the processing content in a spatial filter. ベース成分の除去とオフセット加算を説明するための図である。It is a figure for demonstrating the removal of a base component, and offset addition. 空間フィルタの具体的な構成例を説明するための図である。It is a figure for demonstrating the specific structural example of a spatial filter. グラフに対する面積演算を説明するための図である。It is a figure for demonstrating the area calculation with respect to a graph. 更に他の実施形態に係る超音波診断装置の全体構成を示すブロック図である。It is a block diagram which shows the whole structure of the ultrasonic diagnosing device which concerns on other embodiment. 画像の表示例を示す図である。It is a figure which shows the example of a display of an image. 画像の他の表示例を示す図である。It is a figure which shows the other example of a display of an image.

符号の説明Explanation of symbols

10 探触子、12 送信器、14 受信器、18 デジタルスキャンコンバータ(DSC)、20 フレームメモリ、22 領域設定器、24 二値化処理部、26 閾値スキャン回路、28 ラベリング処理部、30 グラフ作成部、32 画像合成部、34 表示部、36 サンプル数検出部、40 二乗平均値演算部、42 平均閾値演算部、44 分散値演算部、46 半値幅演算部、48 最大ラベル数検出部、50 空間フィルタ、52 部分面積演算部。   10 probes, 12 transmitters, 14 receivers, 18 digital scan converters (DSC), 20 frame memories, 22 area setting units, 24 binarization processing units, 26 threshold scanning circuits, 28 labeling processing units, 30 graph creation Unit, 32 image composition unit, 34 display unit, 36 sample number detection unit, 40 mean square value calculation unit, 42 average threshold value calculation unit, 44 variance value calculation unit, 46 half value width calculation unit, 48 maximum label number detection unit, 50 Spatial filter, 52 Partial area calculator.

Claims (9)

超音波画像に対して、スペックルに重畳したベース成分を除去する前処理を実行するフィルタと、
前記前処理後の超音波画像に対して、スペックルパターンを解析する解析手段と、
を含むことを特徴とする超音波診断装置。
A filter that performs a pre-processing for removing a base component superimposed on speckle on an ultrasonic image;
Analyzing means for analyzing a speckle pattern for the ultrasonic image after the preprocessing;
An ultrasonic diagnostic apparatus comprising:
請求項1記載の装置において、
前記フィルタは、
前記超音波画像上に設定されたウインド内において複数の画素値の平均値を演算する手段と、
前記ウインド内の注目画素の画素値から前記平均値を減算して差分値を求める手段と、
前記差分値にオフセット値を加算し、それを前記注目画像の新しい画素値とする手段と、
を含むことを特徴とする超音波診断装置。
The apparatus of claim 1.
The filter is
Means for calculating an average value of a plurality of pixel values in a window set on the ultrasonic image;
Means for subtracting the average value from the pixel value of the pixel of interest in the window to obtain a difference value;
Means for adding an offset value to the difference value and setting it as a new pixel value of the image of interest;
An ultrasonic diagnostic apparatus comprising:
超音波の送受波により得られたエコーデータに基づいて超音波画像を形成する画像形成手段と、
前記超音波画像を元画像とし、それに対してスペックル評価を行うための前処理を実行する前処理手段と、
前記前処理後の画像に対して、組織の性状を診断するために、スペックルの現れ方を評価するスペックル評価手段と、
前記元画像と前記前処理後の画像とを表示する表示手段と、
を含むことを特徴とする超音波診断装置。
Image forming means for forming an ultrasonic image based on echo data obtained by transmission and reception of ultrasonic waves;
Preprocessing means for performing preprocessing for performing speckle evaluation on the ultrasonic image as an original image;
Speckle evaluation means for evaluating the appearance of speckles for diagnosing tissue properties with respect to the image after the preprocessing;
Display means for displaying the original image and the pre-processed image;
An ultrasonic diagnostic apparatus comprising:
請求項3記載の装置において、
前記表示手段には、更に前記スペックルの評価結果が表示されることを特徴とする超音波診断装置。
The apparatus of claim 3.
The ultrasonic diagnostic apparatus, wherein the display means further displays the speckle evaluation result.
請求項3記載の装置において、
前記前処理手段は、前記超音波画像に対してスペックルに重畳したベース成分を除去するフィルタを含むことを特徴とする超音波診断装置。
The apparatus of claim 3.
The ultrasonic diagnostic apparatus, wherein the preprocessing means includes a filter that removes a base component superimposed on speckles on the ultrasonic image.
請求項3記載の装置において、
前記表示手段には、前記元画像と前記前処理後の画像とを識別するための情報が表示されることを特徴とする超音波診断装置。
The apparatus of claim 3.
Information for identifying the original image and the pre-processed image is displayed on the display means.
請求項3記載の装置において、
前記元画像に対して関心領域を設定する関心領域設定手段を含み、
前記前処理手段は前記関心領域内の部分画像に対して前処理を実行し、
前記表示手段には前記元画像と前記前処理後の関心領域内の部分画像とが表示されることを特徴とする超音波診断装置。
The apparatus of claim 3.
A region of interest setting means for setting a region of interest for the original image;
The preprocessing means performs preprocessing on a partial image in the region of interest;
The ultrasonic diagnostic apparatus according to claim 1, wherein the display means displays the original image and a partial image in the region of interest after the preprocessing.
請求項7記載の装置において、
前記関心領域設定手段により複数の関心領域を設定可能であり、
前記表示手段には前記各関心領域を互いに識別するための情報が表示されることを特徴とする超音波診断装置。
The apparatus of claim 7.
A plurality of regions of interest can be set by the region of interest setting means,
The ultrasonic diagnostic apparatus according to claim 1, wherein information for identifying each region of interest is displayed on the display unit.
超音波の送受波により得られたエコーデータに基づいて超音波画像を形成する画像形成手段と、
前記超音波画像を元画像とし、それに対して関心領域を設定する関心領域設定手段と、
前記元画像に対して前記関心領域を表す第1マーカーを合成してなる第1合成画像を作成する第1合成画像作成手段と、
前記関心領域内の部分画像に対してスペックル評価を行うための前処理を実行する前処理手段と、
前記前処理後の関心領域内の部分画像に対して前記第1マーカーに対応付けられた第2マーカーを合成してなる第2合成画像を作成する第2合成画像作成手段と、
前記前処理後の関心領域内の部分画像に対して、組織の性状を診断するために、スペックルの現れ方を評価するスペックル評価手段と、
前記第1合成画像と前記第2合成画像とを表示する表示手段と、
を含むことを特徴とする超音波診断装置。
Image forming means for forming an ultrasonic image based on echo data obtained by transmission and reception of ultrasonic waves;
Region of interest setting means for setting the region of interest for the ultrasound image as an original image;
First synthesized image creating means for creating a first synthesized image obtained by synthesizing a first marker representing the region of interest with the original image;
Preprocessing means for performing preprocessing for performing speckle evaluation on the partial image in the region of interest;
Second synthesized image creating means for creating a second synthesized image obtained by synthesizing the second marker associated with the first marker with the partial image in the region of interest after the preprocessing;
Speckle evaluation means for evaluating the appearance of speckles for diagnosing the tissue properties for the partial image in the region of interest after the preprocessing;
Display means for displaying the first composite image and the second composite image;
An ultrasonic diagnostic apparatus comprising:
JP2006093397A 2001-08-06 2006-03-30 Ultrasonic diagnostic equipment Expired - Fee Related JP4347860B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2006093397A JP4347860B2 (en) 2001-08-06 2006-03-30 Ultrasonic diagnostic equipment

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2001238124 2001-08-06
JP2002142384 2002-05-17
JP2006093397A JP4347860B2 (en) 2001-08-06 2006-03-30 Ultrasonic diagnostic equipment

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2002227928A Division JP3802462B2 (en) 2001-08-06 2002-08-05 Ultrasonic diagnostic equipment

Publications (2)

Publication Number Publication Date
JP2006212445A true JP2006212445A (en) 2006-08-17
JP4347860B2 JP4347860B2 (en) 2009-10-21

Family

ID=36976097

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006093397A Expired - Fee Related JP4347860B2 (en) 2001-08-06 2006-03-30 Ultrasonic diagnostic equipment

Country Status (1)

Country Link
JP (1) JP4347860B2 (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008206779A (en) * 2007-02-27 2008-09-11 Yamaguchi Univ Ultrasonic diagnostic apparatus and program for determining tissue property
JP2009285051A (en) * 2008-05-28 2009-12-10 Ge Medical Systems Global Technology Co Llc Image diagnostic apparatus, ultrasonic diagnostic apparatus and program
JP2010115372A (en) * 2008-11-13 2010-05-27 Toshiba Corp Three-dimensional ultrasonograph
WO2013183432A1 (en) * 2012-06-07 2013-12-12 日立アロカメディカル株式会社 Method for setting regions of interest and ultrasonic diagnostic device
JP2014036885A (en) * 2013-10-24 2014-02-27 Canon Inc Diagnostic imaging apparatus and diagnostic imaging method
US10828010B2 (en) 2009-05-01 2020-11-10 Canon Kabushiki Kaisha Image diagnosis apparatus and method for dynamically focusing tracked ultrasound probe with multimodal imaging system

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008206779A (en) * 2007-02-27 2008-09-11 Yamaguchi Univ Ultrasonic diagnostic apparatus and program for determining tissue property
JP2009285051A (en) * 2008-05-28 2009-12-10 Ge Medical Systems Global Technology Co Llc Image diagnostic apparatus, ultrasonic diagnostic apparatus and program
JP2010115372A (en) * 2008-11-13 2010-05-27 Toshiba Corp Three-dimensional ultrasonograph
US10828010B2 (en) 2009-05-01 2020-11-10 Canon Kabushiki Kaisha Image diagnosis apparatus and method for dynamically focusing tracked ultrasound probe with multimodal imaging system
WO2013183432A1 (en) * 2012-06-07 2013-12-12 日立アロカメディカル株式会社 Method for setting regions of interest and ultrasonic diagnostic device
CN104334086A (en) * 2012-06-07 2015-02-04 日立阿洛卡医疗株式会社 Method for setting regions of interest and ultrasonic diagnostic device
JPWO2013183432A1 (en) * 2012-06-07 2016-01-28 日立アロカメディカル株式会社 Region of interest setting method and ultrasonic diagnostic apparatus
JP2014036885A (en) * 2013-10-24 2014-02-27 Canon Inc Diagnostic imaging apparatus and diagnostic imaging method

Also Published As

Publication number Publication date
JP4347860B2 (en) 2009-10-21

Similar Documents

Publication Publication Date Title
EP1980210B1 (en) Elastic image display method and elastic image display
EP2016905B1 (en) Ultrasound diagnostic apparatus
US20110125016A1 (en) Fetal rendering in medical diagnostic ultrasound
JP4966578B2 (en) Elastic image generation method and ultrasonic diagnostic apparatus
JP6106190B2 (en) Visualization method of blood and blood likelihood in blood vessel image
JP4921826B2 (en) Ultrasonic diagnostic apparatus and control method thereof
US8834371B2 (en) Ultrasound diagnostic apparatus and ultrasound image processing program
US9039620B2 (en) Ultrasound diagnostic apparatus
JP2005288153A (en) Three-dimensional ultrasonic picture forming apparatus and method
JP2020531074A (en) Ultrasound system with deep learning network for image artifact identification and removal
JP4347860B2 (en) Ultrasonic diagnostic equipment
JP2005193017A (en) Method and system for classifying diseased part of mamma
JP3802462B2 (en) Ultrasonic diagnostic equipment
JP2020018694A (en) Ultrasonic diagnostic device and ultrasonic image processing method
CN111265246A (en) Ultrasonic color imaging processing method and device
JP2008073417A (en) Ultrasonic diagnostic device
EP4076208B1 (en) Systems and methods for assessing a placenta
EP3179269A1 (en) Analysis methods of ultrasound echo signals based on statistics of scatterer distributions
Chang et al. Microcalcification detection in 3-D breast ultrasound
US10105123B2 (en) Analysis methods of ultrasound echo signals based on statistics of scatterer distributions
JP2008212522A (en) Ultrasonic diagnostic apparatus
KR100885435B1 (en) Image processing system and method
JP4708839B6 (en) Biological information analyzer
JP4708839B2 (en) Biological information analyzer
CN115778435A (en) Ultrasonic imaging method and ultrasonic imaging system for fetal face

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090428

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090623

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20090714

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20090716

R150 Certificate of patent (=grant) or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20120724

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20130724

Year of fee payment: 4

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees