JP5222067B2 - Ultrasonic data processor - Google Patents

Ultrasonic data processor Download PDF

Info

Publication number
JP5222067B2
JP5222067B2 JP2008230162A JP2008230162A JP5222067B2 JP 5222067 B2 JP5222067 B2 JP 5222067B2 JP 2008230162 A JP2008230162 A JP 2008230162A JP 2008230162 A JP2008230162 A JP 2008230162A JP 5222067 B2 JP5222067 B2 JP 5222067B2
Authority
JP
Japan
Prior art keywords
luminance distribution
ultrasonic data
processing
luminance
overlap
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.)
Expired - Fee Related
Application number
JP2008230162A
Other languages
Japanese (ja)
Other versions
JP2009273862A (en
Inventor
雅志 中村
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
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 JP2008230162A priority Critical patent/JP5222067B2/en
Publication of JP2009273862A publication Critical patent/JP2009273862A/en
Application granted granted Critical
Publication of JP5222067B2 publication Critical patent/JP5222067B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Description

本発明は超音波データ処理装置に関し、特に対象組織の輪郭のトレース技術に関する。   The present invention relates to an ultrasonic data processing apparatus, and more particularly to a technique for tracing a contour of a target tissue.

超音波診断装置等の超音波データ処理装置は、超音波データに基づいて対象組織の面積や体積を計測する機能を有する。対象組織は、例えば、左室、胎盤、腫瘍、等である。面積演算は二次元超音波データ(断層データ)に基づいて実行され、体積演算は通常、三次元超音波データ(ボリュームデータ)に基づいて実行される。対象組織の面積や体積を計測するためには、対象組織の輪郭をトレースする必要がある。面積演算に当たっては対象組織の外形が特定され、体積演算に当たっては対象組織の表面が特定される。   An ultrasonic data processing apparatus such as an ultrasonic diagnostic apparatus has a function of measuring the area and volume of a target tissue based on ultrasonic data. The target tissue is, for example, the left ventricle, placenta, tumor, and the like. The area calculation is executed based on two-dimensional ultrasonic data (tomographic data), and the volume calculation is usually executed based on three-dimensional ultrasonic data (volume data). In order to measure the area and volume of the target tissue, it is necessary to trace the contour of the target tissue. In the area calculation, the outer shape of the target tissue is specified, and in the volume calculation, the surface of the target tissue is specified.

トレース方式としては手動トレース及び自動トレースが知られている。手動トレースによれば、経験則に基づく目視判断に基づいて輝度差があまりないような構造部分についても比較的正確にトレースを行えるが、面倒であり、負担が大きいという面を指摘できる。また、対象組織表面の全体をトレースすることはかなり困難である。自動トレースによれば、ユーザー負担が生じないが、輝度差が乏しい構造部分については精度良くトレースを行えない場合が多い。なお、特許文献1、2には組織トレースあるいは組織抽出について記載されている。   Manual tracing and automatic tracing are known as trace methods. According to manual tracing, although it is possible to trace relatively accurately even for a structural part where there is not much difference in luminance based on visual judgment based on empirical rules, it can be pointed out that it is cumbersome and burdensome. It is also quite difficult to trace the entire target tissue surface. According to the automatic trace, there is no burden on the user, but in many cases, it is not possible to accurately trace a structural portion with a small luminance difference. Patent Documents 1 and 2 describe tissue tracing or tissue extraction.

特開2002−291750号公報JP 2002-291750 A 特開平7−184892号公報JP-A-7-184892

手動トレース及び自動トレースは上記のようにトレードオフの関係にあるが、その場合に、いずれの方式を選択するのかについての判断をユーザーに委ねると、直感的な判断にならざるを得ず、結果として、各方式の利点を旨く活用できない場合が生じる。そこで、客観的な判断指標が求められている。   Manual tracing and automatic tracing are in a trade-off relationship as described above. In this case, if the user decides which method to select, the result must be intuitive. As a result, there are cases where the advantages of each method cannot be used effectively. Therefore, an objective judgment index is required.

本発明の目的は、手動トレース及び自動トレースを選択する際の客観的な指標を得られるようにすることにある。本発明の他の目的は、その客観的な指標の信頼性を高めることにある。   An object of the present invention is to obtain an objective index when selecting manual tracing and automatic tracing. Another object of the present invention is to increase the reliability of the objective index.

本発明は、二次元又は三次元の生体内空間から取得された超音波データを処理する超音波データ処理装置において、前記生体内空間に対応するデータ空間において、対象組織の外側に外側参照部を設定し、前記対象領域の内側に内側参照部を設定するための設定手段と、前記外側参照部についての外側輝度分布及び前記内側参照部についての内側輝度分布の相互関係に基づいて、前記超音波データに対するトレース方式として、自動トレース処理又は手動トレース処理を判定する判定手段と、を含むことを特徴とする超音波データ処理装置に関する。   The present invention relates to an ultrasonic data processing apparatus for processing ultrasonic data acquired from a two-dimensional or three-dimensional in-vivo space, and in the data space corresponding to the in-vivo space, an outer reference unit is provided outside a target tissue. Setting and setting means for setting an inner reference portion inside the target region, and the ultrasonic wave based on a mutual relationship between an outer luminance distribution for the outer reference portion and an inner luminance distribution for the inner reference portion The present invention relates to an ultrasonic data processing apparatus including a determination unit that determines automatic trace processing or manual trace processing as a trace method for data.

上記構成によれば、対象組織の専ら外側に、二次元又は三次元の形状をもった外側参照部が設定され、同様に、対象組織の専ら内側に、二次元又は三次元の形状をもった内側参照部が設定される。外側参照部は対象組織の完全に外側に設定されるのが望ましいが、部分的なはみ出しが生じていてもよい。同様に、内側参照部は対象組織の完全に内側に設定されるのが望ましいが、部分的なはみ出しが生じていてもよい。要求するトレース方式判定精度との兼ね合いにおいて、はみ出し量の許容度が決まる。外側輝度分布及び内側輝度分布の相互関係は、輝度軸上において、それらがどの程度近い状況にあるのか、あるいは、互いにどの程度分離されている状況にあるのか、を指標するものである。2つの輝度分布が相互に十分に離れているならば、精度良く自動トレースを行える筈であるので、その場合には自動トレースが推奨される。一方、2つの輝度分布が相互に近い関係にあるならば、自動トレースが困難であると推認することができるので、その場合にはマニュアルトレースが推奨される。このように本発明においては、対象部位の内外における輝度分布の比較から、輝度差の傾向を捉え、それをトレース方式の判定で利用するものである。対象部位の周囲全体にわたってトレース方式が一律に決定されるのが望ましいが、部分的な領域ごとにトレース方式を決定することもできる。フレーム単位あるいはボリューム単位でトレース方式を判定することもできる。超音波診断(検査)単位でトレース方式を判定することもできる。あるいは、複数のスライス面の中で、どの面を代表スライス面にするのかの判定に当たって上記手法を利用するようにしてもよい。超音波データ処理装置は超音波診断装置であってもよいし、超音波データを処理する情報処理装置であってもよい。   According to the above configuration, the outer reference portion having a two-dimensional or three-dimensional shape is set exclusively outside the target tissue, and similarly, the outer reference portion has a two-dimensional or three-dimensional shape exclusively inside the target tissue. An inner reference part is set. It is desirable that the outer reference portion is set completely outside the target tissue, but a partial protrusion may occur. Similarly, it is desirable that the inner reference portion is set completely inside the target tissue, but partial protrusion may occur. The allowance for the amount of protrusion is determined in consideration of the required trace method determination accuracy. The mutual relationship between the outer luminance distribution and the inner luminance distribution indicates how close they are on the luminance axis or how separated they are from each other. If the two luminance distributions are sufficiently separated from each other, automatic tracing should be possible with high accuracy. In that case, automatic tracing is recommended. On the other hand, if the two luminance distributions are close to each other, it can be inferred that automatic tracing is difficult. In that case, manual tracing is recommended. As described above, in the present invention, the tendency of the luminance difference is grasped from the comparison of the luminance distribution inside and outside the target region, and is used in the determination of the trace method. Although it is desirable that the tracing method is uniformly determined over the entire periphery of the target region, the tracing method can be determined for each partial region. It is also possible to determine the trace method in units of frames or volumes. The trace method can also be determined in units of ultrasonic diagnosis (examination). Or you may make it utilize the said method in determination of which surface is made into a representative slice surface among several slice surfaces. The ultrasonic data processing apparatus may be an ultrasonic diagnostic apparatus or an information processing apparatus that processes ultrasonic data.

望ましくは、前記相互関係は、前記外側輝度分布及び前記内側輝度分布の重複度である。重複度の具体的な判定方法としては各種の手法をあげることができる。いずれにしても内外の輝度の違いが大きいのか小さいのかの指標が得られる方法を採用するのが望ましい。望ましくは、前記外側参照部の参照により前記外側輝度分布を作成する外側輝度分布手段と、前記内側参照部の参照により前記内側輝度分布を作成する内側輝度分布手段と、前記外側輝度分布及び前記内側輝度分布の一方又は双方を規格化処理する規格化処理手段と、を含み、前記判定手段は、前記規格化処理後の外側輝度分布及び内側輝度分布の重複度に基づいて、自動トレース処理又は手動トレース処理を判定する。外側参照部と内側参照部のそれぞれを構成するピクセル数(あるいはボクセル数)が相違する場合、この規格化処理を適用するのが望ましい。但し、重複部分の絶対量に基づいてトレース方式の選択を行うことも可能である。その場合には規格化処理は不要となる。通常、自動トレース処理では、超音波データの二値化処理が適用される。つまり、各ピクセルの値又は各ボクセルの値(つまりエコー強度に相当する輝度)が所定の閾値と比較されることになるから、トレース方式の判別に当たって輝度情報を利用するのは合理的である。   Preferably, the interrelationship is a degree of overlap between the outer luminance distribution and the inner luminance distribution. Various methods can be cited as specific methods for determining the degree of overlap. In any case, it is desirable to adopt a method that can provide an indicator of whether the difference in brightness between the inside and outside is large or small. Preferably, an outer luminance distribution unit that creates the outer luminance distribution by referring to the outer reference unit, an inner luminance distribution unit that creates the inner luminance distribution by referring to the inner reference unit, the outer luminance distribution, and the inner luminance distribution unit. Normalization processing means for normalizing one or both of the luminance distributions, and the determination means performs automatic trace processing or manual operation based on the overlap between the outer luminance distribution and the inner luminance distribution after the normalization processing. Determine trace processing. When the number of pixels (or the number of voxels) constituting each of the outer reference portion and the inner reference portion is different, it is desirable to apply this normalization processing. However, the trace method can be selected based on the absolute amount of the overlapping portion. In that case, normalization processing becomes unnecessary. Usually, in automatic trace processing, binarization processing of ultrasonic data is applied. That is, since the value of each pixel or the value of each voxel (that is, the luminance corresponding to the echo intensity) is compared with a predetermined threshold value, it is reasonable to use the luminance information in determining the trace method.

望ましくは、前記判定手段は、輝度軸上における前記規格化処理後の外側輝度分布及び内側輝度分布の重複区間を特定する手段と、前記外側輝度分布の内で前記重複区間に属する部分を演算する手段と、前記内側輝度分布の全体に対する前記重複区間の部分の割合に基づいて前記重複度を演算する手段と、前記重複度と所定の基準値との比較結果に基づいて前記自動トレース処理又は前記手動トレース処理を判定する手段と、を含む。   Preferably, the determination means calculates a portion belonging to the overlapping section in the outer luminance distribution, and means for specifying an overlapping section of the outer luminance distribution and the inner luminance distribution after the normalization processing on the luminance axis. Means for calculating the degree of overlap based on a ratio of the overlapping section portion to the entire inner luminance distribution, and the automatic trace processing based on a comparison result between the degree of overlap and a predetermined reference value, or Means for determining a manual trace process.

望ましくは、前記設定手段は、前記超音波データに基づいて前記対象組織の代表断面を表す断層画像を生成する手段と、前記断層画像上において、前記外側参照部としての二次元閉ループのユーザー指定を受け付ける手段と、前記断層画像上において、前記内側参照部としての二次元閉ループのユーザー指定を受け付ける手段と、を含む。望ましくは、前記設定手段は、前記超音波データに基づいて前記対象組織についての複数の断面を表す複数の断層画像を生成する手段と、前記複数の断層画像上において、前記外側参照部としての三次元球状面のユーザー指定を受け付ける手段と、前記複数の断層画像上において、前記内側参照部としての三次元球状面のユーザー指定を受け付ける手段と、を含む。   Preferably, the setting unit generates a tomographic image representing a representative cross section of the target tissue based on the ultrasound data, and performs user designation of a two-dimensional closed loop as the outer reference unit on the tomographic image. Means for receiving, and means for receiving a user designation of a two-dimensional closed loop as the inner reference portion on the tomographic image. Preferably, the setting means generates a plurality of tomographic images representing a plurality of cross sections of the target tissue based on the ultrasound data, and a tertiary as the outer reference unit on the plurality of tomographic images. Means for accepting user designation of the original spherical surface, and means for accepting user designation of the three-dimensional spherical surface as the inner reference portion on the plurality of tomographic images.

また、本発明は、二次元又は三次元の生体内空間から取得された超音波データを処理する超音波データ処理装置において、前記生体内空間に対応するデータ空間において、対象組織の外側に外部参照部を設定し、前記対象組織の内側に内側参照部を設定するための設定手段と、前記外側参照部についての外側輝度分布及び前記内側参照部についての内側輝度分布に基づいて、前記超音波データに対するトレース方式として、自動トレース処理又は手動トレース処理を判定する判定手段と、を含み、前記判定手段は、前記外部参照部の参照により前記外側輝度分布を作成する外側輝度分布作成手段と、前記内側参照部の参照により前記内側輝度分布を作成する内側輝度分布作成手段と、輝度軸上における前記外側輝度分布と前記内側輝度分布の有効重複範囲を用いて求められる重複度に基づいて、自動トレース処理又は手動トレース処理を判定する処理方式判定手段と、を含み、前記処理方式判定手段は、前記外側輝度分布と前記内側輝度分布の内の少なくとも一方について裾野部分を除外することにより前記有効重複範囲を決定する、ことを特徴とする。   Further, the present invention provides an ultrasonic data processing apparatus for processing ultrasonic data acquired from a two-dimensional or three-dimensional in-vivo space, and an external reference outside the target tissue in the data space corresponding to the in-vivo space. The ultrasonic data based on a setting means for setting a portion and setting an inner reference portion inside the target tissue, an outer luminance distribution for the outer reference portion, and an inner luminance distribution for the inner reference portion Determination means for determining automatic trace processing or manual trace processing as a tracing method for the external luminance distribution creating means for creating the outer luminance distribution by referring to the external reference unit; An inner luminance distribution creating means for creating the inner luminance distribution by referring to a reference unit; and the outer luminance distribution and the inner luminance distribution on the luminance axis. Processing method determination means for determining automatic trace processing or manual trace processing based on the degree of overlap obtained using the effective overlap range, and the processing method determination means includes the outer luminance distribution and the inner luminance distribution. The effective overlapping range is determined by excluding the base portion of at least one of the two.

上記構成によれば、トレース方式の判定に当たって、外側輝度分布と内側輝度分布における少なくとも一方の裾野部分を演算対象から除外して重複区間(有効重複区間)を決定できるので、輝度分布における低レベルの部分が長く伸びてそれによって判定される重複区間が過大に見積られてしまう問題に対処することができる。一般に、裾野部分のカットは、外側輝度部分における両側、内側輝度部分における両側に適用することが可能である。但し、2つの輝度分布の重なり区間に着目するならば、外側輝度分布の低輝度側及び内側輝度分布の高輝度側において除外処理が適用されるのが望ましく、特に、後者について除外処理が適用されるのが望ましい。   According to the above configuration, when determining the trace method, it is possible to determine an overlapping section (effective overlapping section) by excluding at least one base portion of the outer luminance distribution and the inner luminance distribution from the calculation target. It is possible to cope with the problem that the overlapping section determined by the extended portion is overestimated. In general, the cut of the skirt portion can be applied to both sides of the outer luminance portion and both sides of the inner luminance portion. However, if attention is paid to the overlapping section of two luminance distributions, it is desirable that the exclusion process is applied on the low luminance side of the outer luminance distribution and the high luminance side of the inner luminance distribution, and in particular, the exclusion process is applied to the latter. Is desirable.

望ましくは、前記裾野部分は、前記内側輝度分布における高輝度側の部分であって所定レベル以下の部分である。望ましくは、前記裾野部分は、前記内側輝度分布における高輝度側の部分であって前記内側輝度分布のばらつき度合いを示す統計的指標に基づいて決定される部分である。結果として、裾野部分が除外される処理であれば上記以外の方法を採用することも可能である。   Preferably, the skirt portion is a portion on the high luminance side in the inner luminance distribution and a portion below a predetermined level. Preferably, the base portion is a portion on the high luminance side in the inner luminance distribution and is determined based on a statistical index indicating a degree of variation of the inner luminance distribution. As a result, it is possible to adopt a method other than the above as long as it is a process in which the base portion is excluded.

望ましくは、前記処理方式判定手段は、前記内側輝度部分の裾野部分を除外することにより前記有効重複区間を決定する手段と、前記内側輝度分布における前記裾野部分を除外した後の主要部と、前記外側輝度分布における前記有効重複区間に属する部分と、の比率を前記重複度として演算する手段と、前記比率を一定値を比較することにより、前記自動トレース処理又は前記手動トレース処理を判定する手段と、を含む。内側輝度分布には対象組織内の様々な構造が反映され、それが広がることも多いので、また高輝度ノイズの影響により高輝度側にノイズ波形が生じ易いので、内側輝度分布に対して裾野部分を演算上除外するのが望ましい。   Preferably, the processing method determination means determines the effective overlapping section by excluding the skirt portion of the inner luminance portion, the main portion after excluding the skirt portion in the inner luminance distribution, Means for calculating a ratio of the portion belonging to the effective overlap section in the outer luminance distribution as the degree of overlap, and means for determining the automatic trace process or the manual trace process by comparing the ratio with a constant value; ,including. Various structures within the target tissue are reflected in the inner luminance distribution and often spread, and noise waveforms are likely to occur on the high luminance side due to the influence of high luminance noise. Is preferably excluded from the calculation.

以上説明したように、本発明によれば、手動トレース及び自動トレースの両者の利点を旨く発揮させることができる。特に、それらを選択する際の客観的な指標を得られる。   As described above, according to the present invention, the advantages of both manual tracing and automatic tracing can be achieved. In particular, an objective index for selecting them can be obtained.

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

図1には、本発明に係る超音波データ処理装置の好適な実施形態が示されており、図1はその全体構成を示すブロック図である。この超音波データ処理装置は、本実施形態において超音波診断装置により構成されている。ただし、その図1に示される構成の全部または一部が情報処理装置(コンピュータ)によって構成されてもよい。   FIG. 1 shows a preferred embodiment of an ultrasonic data processing apparatus according to the present invention, and FIG. 1 is a block diagram showing the overall configuration thereof. This ultrasonic data processing apparatus is constituted by an ultrasonic diagnostic apparatus in the present embodiment. However, all or part of the configuration shown in FIG. 1 may be configured by an information processing device (computer).

3Dプローブ10は、生体内における三次元空間に対して超音波の送受波を行う送受波器である。3Dプローブ10は本実施形態において2Dアレイ振動子を有している。2Dアレイ振動子は二次元配列された複数の振動素子により構成される。それによって超音波ビームが形成され、超音波ビームは電子的に二次元走査される。これによって三次元データ取り込み空間が形成される。電子走査方式としては、電子セクタ走査、電子リニア走査等が知られている。2Dアレイ振動子に代えて、1Dアレイ振動子とそれを機械的に走査する機構とを設けてもよい。   The 3D probe 10 is a transducer that transmits and receives ultrasonic waves to and from a three-dimensional space in a living body. In this embodiment, the 3D probe 10 has a 2D array transducer. The 2D array transducer is constituted by a plurality of vibration elements arranged two-dimensionally. Thereby, an ultrasonic beam is formed, and the ultrasonic beam is electronically two-dimensionally scanned. As a result, a three-dimensional data capturing space is formed. As electronic scanning methods, electronic sector scanning, electronic linear scanning, and the like are known. Instead of the 2D array transducer, a 1D array transducer and a mechanism for mechanically scanning the 1D array transducer may be provided.

送受信部12は、送信ビームフォーマ及び受信ビームフォーマとして機能する。すなわち、送信時においては、送受信部12から2Dアレイ振動子に対して複数の送信信号が供給される。これによって送信ビームが形成される。受信時において、生体内からの反射波は2Dアレイ振動子にて受波され、これによって2Dアレイ振動子から複数の受信信号が送受信部12へ出力される。送受信部12は複数の受信信号に対する整相加算処理を実行し、これによって整相加算後の受信信号(ビームデータ)を出力する。そのビームデータは3Dメモリ14に格納される。   The transmission / reception unit 12 functions as a transmission beamformer and a reception beamformer. That is, at the time of transmission, a plurality of transmission signals are supplied from the transmission / reception unit 12 to the 2D array transducer. As a result, a transmission beam is formed. At the time of reception, the reflected wave from the living body is received by the 2D array transducer, whereby a plurality of received signals are output from the 2D array transducer to the transmission / reception unit 12. The transmission / reception unit 12 performs phasing addition processing on a plurality of reception signals, and thereby outputs a reception signal (beam data) after phasing addition. The beam data is stored in the 3D memory 14.

3Dメモリ14は、生体内の三次元空間に対応したデータ空間を有し、入力される複数のビームデータは3Dメモリ14内に格納される。その格納の際に、あるいはデータの読み出しの際に、送受波座標系からデータ座標系への座標変換が実行される。ビームデータは超音波ビーム上に沿って並んだ複数のボクセルデータにより構成される。各ボクセルデータは、本実施形態において、3Dメモリ14への書き込み時点でそれに対応するアドレスにマッピングされる。3Dメモリ14への各ボクセルデータへの書き込みに先立って、必要に応じて信号処理が実行される。その信号処理には例えば二値化処理等が含まれてもよい。   The 3D memory 14 has a data space corresponding to the three-dimensional space in the living body, and a plurality of input beam data are stored in the 3D memory 14. Coordinate conversion from the transmission / reception coordinate system to the data coordinate system is executed when the data is stored or when data is read. The beam data is composed of a plurality of voxel data arranged along the ultrasonic beam. In the present embodiment, each voxel data is mapped to an address corresponding to the voxel data at the time of writing to the 3D memory 14. Prior to writing to each voxel data in the 3D memory 14, signal processing is executed as necessary. The signal processing may include, for example, binarization processing.

データ処理部16は、3Dメモリ14に格納された複数のビームデータ(ボリュームデータ)に基づいて、各種のデータ処理を実行するモジュールである。データ処理部16は、本実施形態において、画像形成部40、自動トレース部42、及び、体積/面積演算部44を備えている。画像形成部40は、ボリュームデータに基づいてボリュームレンダリング法等の処理を適用し、これによって三次元超音波画像を形成するユニットである。また、この画像形成部40は、三次元空間に設定された3つの切断面に対応する3つの断層画像をいわゆるトリプレーン画像として形成する機能も有している。さらに、この画像形成部40はユーザーにより任意に設定された切断面に対応する任意断層画像を形成する機能も有している。それらの画像データは、表示処理部18を介して表示部20へ出力される。   The data processing unit 16 is a module that executes various types of data processing based on a plurality of beam data (volume data) stored in the 3D memory 14. In this embodiment, the data processing unit 16 includes an image forming unit 40, an automatic tracing unit 42, and a volume / area calculating unit 44. The image forming unit 40 is a unit that applies a process such as a volume rendering method based on volume data, thereby forming a three-dimensional ultrasonic image. The image forming unit 40 also has a function of forming three tomographic images corresponding to three cut planes set in a three-dimensional space as so-called triplane images. Further, the image forming unit 40 has a function of forming an arbitrary tomographic image corresponding to a cutting plane arbitrarily set by the user. Those image data are output to the display unit 20 via the display processing unit 18.

自動トレース部42は、三次元空間に設定された複数のスライス面上においてトレース処理を行う場合において、自動トレースが判定されたスライス面についてオートトレース処理を実行するユニットである。そのトレース処理に先立って、ボリュームデータに対するあるいはスライスデータに対する二値化処理が実行される。オートトレース処理はいわゆるエッジ検出法などに基づいて実現することが可能である。オートトレース法としては従来から各種の手法が提案されている。一方、本実施形態においては各スライス面上においてマニュアルトレース処理を行うことも可能である。マニュアルトレース処理は、後に説明する操作パネル38上のポインティングデバイスを利用して断層画像上において対象組織の輪郭を手動でトレースするものである。その座標データはデータ処理部16に与えられ、データ処理部16はその座標データを利用して以下に説明する計測処理を実行する。三次元的に存在する対象組織について自動トレース処理、あるいはマニュアルトレース処理が一律に設定されるようにしてもよいし、各スライス面あるいは断層面ごとに自動トレース処理またはマニュアルトレース処理が切り替え適用されてもよい。さらに、それらのトレース方法の選択が超音波診断すなわち検査を単位として行われるようにしてもよい。   The automatic trace unit 42 is a unit that executes auto-trace processing for slice planes for which automatic tracing is determined when performing trace processing on a plurality of slice planes set in a three-dimensional space. Prior to the trace processing, binarization processing for volume data or slice data is executed. The auto trace process can be realized based on a so-called edge detection method or the like. Various methods have been proposed as the auto-trace method. On the other hand, in the present embodiment, manual trace processing can also be performed on each slice plane. The manual trace process manually traces the contour of the target tissue on the tomographic image using a pointing device on the operation panel 38 described later. The coordinate data is given to the data processing unit 16, and the data processing unit 16 executes measurement processing described below using the coordinate data. Automatic trace processing or manual trace processing may be set uniformly for target tissues that exist three-dimensionally, or automatic trace processing or manual trace processing is switched and applied to each slice plane or tomographic plane. Also good. Further, the selection of the tracing method may be performed in units of ultrasonic diagnosis, that is, examination.

体積/面積演算部44は、自動トレース処理の結果及び/又はマニュアルトレース処理の結果に基づいて対象組織についての面積を演算し、また体積を演算する機能を有している。例えば、対象組織が左室であれば、その左室を横切る複数のスライス面上においてループ形状のトレースラインが形成され、各スライス面上の面積が演算された上で、それらに一定の厚みを与えつつ加算することにより左室の体積が演算される。また、胎児や胎盤などを計測対象とすることも可能である。体積/面積演算部44の演算結果は表示処理部18を経由して表示部20に送られ、その演算結果は数値として、あるいはグラフとして画面上に表示される。表示部20における画面上にトレース結果が表示されるようにしてもよい。また、表示部20の画面上には後に説明する外側ROI(関心領域)及び内側ROIの設定の際に必要に応じて断層画像や各ROIを表すグラフィックが表示される。   The volume / area calculation unit 44 has a function of calculating the area of the target tissue based on the result of the automatic trace process and / or the result of the manual trace process, and calculating the volume. For example, if the target tissue is the left ventricle, loop-shaped trace lines are formed on a plurality of slice planes crossing the left ventricle, and the area on each slice plane is calculated, and then a certain thickness is given to them. The volume of the left ventricle is calculated by adding while giving. It is also possible to measure a fetus or placenta. The calculation result of the volume / area calculation unit 44 is sent to the display unit 20 via the display processing unit 18, and the calculation result is displayed on the screen as a numerical value or as a graph. The trace result may be displayed on the screen of the display unit 20. Further, on the screen of the display unit 20, a tomographic image and a graphic representing each ROI are displayed as necessary when setting an outer ROI (region of interest) and an inner ROI, which will be described later.

次に、トレース方式判定ユニット22について説明する。このトレース方式判定ユニット22は、本実施形態においてヒストグラム作成部24,26、規格化部28,30、重複計測部32及び判定部34を備えている。   Next, the trace method determination unit 22 will be described. In this embodiment, the trace method determination unit 22 includes histogram creation units 24 and 26, normalization units 28 and 30, a duplicate measurement unit 32, and a determination unit 34.

ヒストグラム作成部24は、対象組織の外側にユーザーにより外側ROIが設定された場合に、その外側ROIの表面に相当する複数のボクセルデータに基づきヒストグラムを作成するユニットである。同様に、ヒストグラム作成部26は、対象組織の内側にユーザーにより内側ROIが設定された場合に、その表面に相当する複数のボクセルデータに基づいてヒストグラムを作成するユニットである。外側ROIは専ら対象組織の外側に設定されればよく、同様に、内側ROIは専ら対象組織の内側に設定されればよい。すなわち部分的にはみ出しが発生していても、最終的なトレース方式の判定にあたって一定の精度が得られるのであれば、そのような部分的なはみ出しは問題とはならない。したがって、各ROIを自動的に設定するようにしてもよい。各ROIは二次元形状あるいは三次元形状として設定することが可能である。その場合においては、例えば楕円あるいは楕円球のような形状を利用して各ROIが設定されるようにしてもよい。その設定にあたっては、注目しているスライス面に相当する断層画像を表示部20上に表示するようにしてもよいし、いわゆるトリプレーン画像を利用して各ROIの設定を行うようにしてもよい。さらに、三次元超音波画像上において立体的にROIを設定するようにしてもよい。ただし、三次元超音波画像上に三次元のROIを直ちに設定するのは困難であるため、トリプレーン画像の内の2つの画像を使って例えば楕円球の長軸及び短軸を指定することによりROIを設定する各ROIの表示上のボクセル群ではなく、各ROIの表層部内のボクセル群に基づいてヒストグラムを作成してもよい。   The histogram creation unit 24 is a unit that creates a histogram based on a plurality of voxel data corresponding to the surface of the outer ROI when an outer ROI is set by the user outside the target tissue. Similarly, the histogram creation unit 26 is a unit that creates a histogram based on a plurality of voxel data corresponding to the surface when an inner ROI is set by the user inside the target tissue. The outer ROI only needs to be set outside the target tissue, and similarly, the inner ROI only needs to be set inside the target tissue. That is, even if partial protrusion occurs, such partial protrusion is not a problem as long as a certain accuracy can be obtained in the determination of the final trace method. Therefore, each ROI may be set automatically. Each ROI can be set as a two-dimensional shape or a three-dimensional shape. In that case, each ROI may be set using a shape such as an ellipse or an elliptic sphere. In the setting, a tomographic image corresponding to the slice plane of interest may be displayed on the display unit 20, or each ROI may be set using a so-called triplane image. . Furthermore, the ROI may be set three-dimensionally on the three-dimensional ultrasonic image. However, since it is difficult to immediately set a three-dimensional ROI on a three-dimensional ultrasound image, for example, by specifying the major axis and minor axis of an elliptical sphere using two of the triplane images. A histogram may be created based on the voxel group in the surface layer portion of each ROI instead of the voxel group on the display of each ROI for setting the ROI.

本実施形態においては、ROIの表面あるいは輪郭が利用される。すなわち三次元的な表面あるいは二次元的な輪郭上に存在する複数のボクセルデータが参照され、各ボクセルデータの値に基づいて外側ヒストグラム及び内側ヒストグラムが作成される。具体例については後に説明する。規格化ユニット27は、2つの規格化部28,30により構成される。規格化ユニット27は、作成された2つのヒストグラムのうちで少なくとも一方について規格化処理を適用することにより、2つのヒストグラムの互いの関係を計測するための前処理を実行する。本実施形態においては、それぞれのヒストグラムが規格化されており、すなわち互いのヒストグラム上の面積が同一とされている。   In this embodiment, the surface or contour of the ROI is used. That is, a plurality of voxel data existing on a three-dimensional surface or a two-dimensional contour is referred to, and an outer histogram and an inner histogram are created based on the value of each voxel data. Specific examples will be described later. The normalization unit 27 is composed of two normalization units 28 and 30. The normalization unit 27 performs preprocessing for measuring the relationship between the two histograms by applying a normalization process to at least one of the two generated histograms. In the present embodiment, the respective histograms are standardized, that is, the areas on the histograms are the same.

重複計測部32は、上記のように作成された外側ヒストグラム及び内側ヒストグラムにおいて相互間のオーバーラップ部分である重複部分を特定し、それを判定基準として演算するユニットである。2つのヒストグラムの重複部分が少ない場合には、対象組織の内外の分離が良好であることが推認でき、その一方において、2つのヒストグラムの重複部分が大きい場合には対象組織の内外が互いに近い輝度状態になっていることが推認できる。したがって、本実施形態においてはそのような物理的な状況あるいは経験則に基づいて、重複部分を指標として利用し、それを基準としてトレース方式の判定を行っている。その判定は具体的には判定部34において実行されている。重複度の演算方法及び具体的な判定方法については後に詳述する。いずれにしても、判定部34は、2つのヒストグラムの重複関係に基づいて適切なトレース方式を選択する。具体的にはマニュアルトレースが適当であるか、あるいは自動トレースが適当であるかを状況に応じて判定している。その判定結果はデータ処理部16に送られている。   The overlap measurement unit 32 is a unit that specifies an overlap portion that is an overlap portion between the outer histogram and the inner histogram created as described above, and calculates this as a determination criterion. When there is little overlap between the two histograms, it can be inferred that the inside and outside separation of the target tissue is good. On the other hand, when the two histograms have a large overlap, the inside and outside of the target tissue are close to each other. It can be inferred that it is in a state. Therefore, in this embodiment, based on such a physical situation or empirical rule, the overlapping portion is used as an index, and the trace method is determined based on the overlapping portion. Specifically, the determination is performed in the determination unit 34. A method for calculating the degree of overlap and a specific determination method will be described in detail later. In any case, the determination unit 34 selects an appropriate trace method based on the overlapping relationship between the two histograms. Specifically, it is determined depending on the situation whether manual tracing is appropriate or automatic tracing is appropriate. The determination result is sent to the data processing unit 16.

データ処理部16では、判定部の判定結果に応じて対象となるスライス面について自動トレースを実行するのか、あるいはマニュアルトレースの入力を受け付けるのかを判断している。データ処理部16は、上記判定部34の判定結果がマニュアルトレースを示す場合に、表示画面上にマニュアルトレースを促す旨を表示する。その表示を契機として、ユーザーは、表示されている対象画像上においてマニュアルトレースを実行する。その際に入力される座標データがデータ処理部16において認識される。もっとも、この段階においてユーザーの意志を尊重し、ユーザーが選択した方式に応じて実際のトレースが実行されてもよい。その場合においては、ユーザーに推奨されるべき判定方式の情報を提供するだけでよい。   The data processing unit 16 determines whether to execute automatic tracing for a target slice plane or to accept manual trace input according to the determination result of the determination unit. When the determination result of the determination unit 34 indicates a manual trace, the data processing unit 16 displays on the display screen a prompt for manual tracing. In response to the display, the user performs manual tracing on the displayed target image. The data processing unit 16 recognizes the coordinate data input at that time. However, the user's will is respected at this stage, and the actual trace may be executed according to the method selected by the user. In that case, it is only necessary to provide information on a determination method to be recommended to the user.

制御部36は、CPU及び動作プログラムによって構成され、図1に示される各構成の動作制御を行っている。制御部36には操作パネル38が接続されている。操作パネル38はキーボードやトラックボール等により構成される。操作パネル38を利用して、ユーザーは外側ROI及び内側ROIの設定を行うことができ、またマニュアルトレースを行うことができ、さらに必要に応じて判定部34に入力する閾値αの入力指定を行える。   The control unit 36 includes a CPU and an operation program, and performs operation control of each configuration illustrated in FIG. An operation panel 38 is connected to the control unit 36. The operation panel 38 is composed of a keyboard, a trackball, and the like. Using the operation panel 38, the user can set the outer ROI and the inner ROI, can perform manual tracing, and can specify the input of the threshold value α to be input to the determination unit 34 as necessary. .

次に、図2乃至図7を用いて、本実施形態に係るトレース方式の判別方法について詳述する。   Next, a tracing method determination method according to the present embodiment will be described in detail with reference to FIGS.

図2には、断層画像50が示されている。断層画像50には対象組織(正確には対象組織の断面)52が含まれている。符号52Aは対象組織52の輪郭を表している。その輪郭52Aの外側にユーザーによって外側ROI54が設定される。同様に、輪郭52Aの内側にユーザーにより内側ROI56が設定される。それらのROI54,56は二次元平面状の形状であってもよいし、三次元空間内の立体であってもよい。例えば楕円の長軸、短軸及び中心座標を指定することにより、各ROI54,56が設定されるようにしてもよい。また、さらにもう1つの長軸を指定することにより、それぞれのROIを楕円球として定義することも可能である。本実施形態においては楕円あるいは楕円球が用いられていたが、ROIの形状としてはそれ以外にも各種の形状を採用しうる。   A tomographic image 50 is shown in FIG. The tomographic image 50 includes a target tissue (more precisely, a cross section of the target tissue) 52. Reference numeral 52 </ b> A represents the contour of the target tissue 52. An outer ROI 54 is set by the user outside the contour 52A. Similarly, the inner ROI 56 is set by the user inside the contour 52A. These ROIs 54 and 56 may have a two-dimensional planar shape or a solid in a three-dimensional space. For example, the ROIs 54 and 56 may be set by designating the major axis, minor axis, and center coordinates of the ellipse. Further, by specifying another major axis, each ROI can be defined as an elliptical sphere. In the present embodiment, an ellipse or an elliptic sphere is used, but various shapes can be adopted as the shape of the ROI.

図3には、外側ROIの表面上に存在する複数のボクセルの値によって構成されるヒストグラム100が示されている。横軸は輝度すなわちボクセルの値を表しており、縦軸は規格化後の度数を表している。すなわち、ヒストグラム100は規格化後のヒストグラムである。このヒストグラム100は符号108で示される輝度から符号109で示される範囲内に分布している。一方、図4には、内側ROIの表面上に存在する複数のボクセルの値により構成されるヒストグラム102が示されている。このヒストグラム102も規格化後のヒストグラムであり、それは符号104で示される輝度から符号106で示される輝度の間に分布している。   FIG. 3 shows a histogram 100 composed of a plurality of voxel values existing on the surface of the outer ROI. The horizontal axis represents the luminance, that is, the voxel value, and the vertical axis represents the frequency after normalization. That is, the histogram 100 is a normalized histogram. The histogram 100 is distributed within the range indicated by reference numeral 109 from the luminance indicated by reference numeral 108. On the other hand, FIG. 4 shows a histogram 102 composed of a plurality of voxel values existing on the surface of the inner ROI. This histogram 102 is also a normalized histogram, which is distributed between the luminance indicated by reference numeral 104 and the luminance indicated by reference numeral 106.

図5には、2つのヒストグラム100,102の相互関係が示されている。すなわちそれらのヒストグラム100,102が同一の輝度軸上に表されている。ここでAはヒストグラム102が分布している範囲を示すものである。S1はヒストグラム100の面積を表しており、S2はヒストグラム102の面積を表している。ヒストグラム100において区間Aに属する部分の面積がΔS1で表されている。   FIG. 5 shows the mutual relationship between the two histograms 100 and 102. That is, the histograms 100 and 102 are represented on the same luminance axis. Here, A indicates a range in which the histogram 102 is distributed. S 1 represents the area of the histogram 100, and S 2 represents the area of the histogram 102. In the histogram 100, the area of the portion belonging to the section A is represented by ΔS1.

本実施形態においては、2つのヒストグラム100,102の重複関係を評価するため、重複度を表す指標εとしてΔS1/S2が用いられている。すなわち内側ヒストグラムに対応する面積S2に対する、外側ヒストグラムにおける重複部分の面積ΔS1の比εが用いられている。その比εが大きければ2つのヒストグラム100,102が大きくオーバーラップしていることを推認でき、その一方、その比εが小さければ2つのヒストグラムがほどよく分離していることを推認できる。したがって、その比εを上述した所定の閾値αと比較することにより、指標εがαよりも大きければ自動トレースが困難であることが予想されるためにマニュアルトレースが判定され、逆に、その比εがαよりも小さければ対象部位の内外において輝度差が明瞭であることを推認できるために自動トレースが判定される。上述した判定方法は一例であり、いずれにしても2つのヒストグラムの相互関係が定量化されるように判定方法を定めればよい。   In the present embodiment, ΔS1 / S2 is used as an index ε representing the degree of overlap in order to evaluate the overlapping relationship between the two histograms 100 and 102. That is, the ratio ε of the area ΔS1 of the overlapping portion in the outer histogram to the area S2 corresponding to the inner histogram is used. If the ratio ε is large, it can be inferred that the two histograms 100 and 102 are largely overlapped. On the other hand, if the ratio ε is small, it can be inferred that the two histograms are well separated. Therefore, by comparing the ratio ε with the above-described predetermined threshold value α, it is predicted that automatic tracing is difficult if the index ε is larger than α, so that manual tracing is determined. If ε is smaller than α, it can be inferred that the brightness difference is clear inside and outside the target region, and therefore automatic tracing is determined. The determination method described above is an example, and in any case, the determination method may be determined so that the mutual relationship between the two histograms is quantified.

図6及び図7には、別の例が示されている。判定方法としては同一のものが利用されている。断層画像60上には組織62が存在し、その組織62は対象組織62Aとそれに連なる組織62Bとにより構成されている。例えばこのような状況は胎盤と胎児との関係において認められる。このような場合、対象組織62Aの輪郭63の外側に外側ROI64が設定され、その内側に内側ROI66が設定される。ただし、この例では外側ROI64は本来計測の対象とはならない組織62Bまでくい込んでいる。このような場合において、2つのヒストグラムを作成すると、図7に示すようなものとなる。すなわち、外側ヒストグラム110と内側ヒストグラム112とが形成され、両者はかなりオーバーラップしている。外側ヒストグラム110においては2つのピークが認められ、低い方のピークは内側ヒストグラム112内に位置している。   Another example is shown in FIG. 6 and FIG. The same determination method is used. A tissue 62 exists on the tomographic image 60, and the tissue 62 includes a target tissue 62A and a tissue 62B connected to the target tissue 62A. For example, this situation is observed in the relationship between the placenta and the fetus. In such a case, the outer ROI 64 is set outside the contour 63 of the target tissue 62A, and the inner ROI 66 is set inside thereof. However, in this example, the outer ROI 64 penetrates to the tissue 62B that is not originally a measurement target. In such a case, when two histograms are created, the result is as shown in FIG. That is, the outer histogram 110 and the inner histogram 112 are formed, and both are substantially overlapped. Two peaks are observed in the outer histogram 110, and the lower peak is located in the inner histogram 112.

この場合において、上記同様の判定手法を適用すると、内側ヒストグラム112が存在している輝度114から輝度116の間が区間Aとして特定される。一方、外側ヒストグラム110における輝度118から輝度119の範囲において、Aに属する部分の面積がΔS10として特定される。内側ヒストグラム112の面積はS12である。このような前提のもと、比εとしてΔS10/S12が求められ、その比εが所定の閾値αと比較されることになる。この場合においては、かなりオーバーラップが生じているために、αよりもεの方が大きくなり、その結果自動トレースよりもマニュアルトレースの方がふさわしいと判定されることになる。αは予めユーザーにより設定することもできるし、画像の状況に応じて適応的にすなわち自動的に設定することも可能である。   In this case, when the same determination method as described above is applied, a section between the luminance 114 and the luminance 116 where the inner histogram 112 exists is specified as the section A. On the other hand, in the range from the luminance 118 to the luminance 119 in the outer histogram 110, the area of the portion belonging to A is specified as ΔS10. The area of the inner histogram 112 is S12. Under such a premise, ΔS10 / S12 is obtained as the ratio ε, and the ratio ε is compared with a predetermined threshold value α. In this case, since there is a considerable overlap, ε is larger than α, and as a result, it is determined that manual tracing is more suitable than automatic tracing. α can be set in advance by the user, or can be set adaptively, that is, automatically according to the state of the image.

2つのヒストグラムの関係を分析する手法としては各種の方法を採用することができる。例えば、図8に示すように2つのヒストグラム120,122が形成されている状況において、各ヒストグラム120,122のピークの輝度P1,P2と各ヒストグラム120,122のヒストグラム幅σ1,σ2を利用して2つのヒストグラムの重複関係を評価するようにしてもよい。ただし、それらのパラメーターによると必ずしもオーバーラップ部分の大小を直接的に求められない可能性があるため、2つのヒストグラム120,122における重複部分をより直接的に評価するのが望ましい。例えば、2つのヒストグラム120,122において重複している面積ΔS’を求め、それを絶対的に利用してトレース方式の判別に利用するようにしてもよい。あるいは重複部分が生じている区間ΔB’と一方のヒストグラムの幅Bとの比を使って重複関係を評価するようにしてもよい。この場合において、ΔB’を単独で利用して重複関係の大小を評価するようにしてもよい。いずれにしても、2つのヒストグラムの交わりの大きさを定量的に評価することにより、その評価結果からトレース方式を客観的に判断することが可能となる。   Various methods can be employed as a method for analyzing the relationship between the two histograms. For example, in a situation where two histograms 120 and 122 are formed as shown in FIG. 8, the peak luminances P1 and P2 of the histograms 120 and 122 and the histogram widths σ1 and σ2 of the histograms 120 and 122 are used. You may make it evaluate the duplication relation of two histograms. However, since it may not always be possible to directly determine the size of the overlap portion according to these parameters, it is desirable to evaluate the overlap portions in the two histograms 120 and 122 more directly. For example, the overlapping area ΔS ′ in the two histograms 120 and 122 may be obtained and used absolutely for the discrimination of the trace method. Alternatively, the overlapping relationship may be evaluated using the ratio between the interval ΔB ′ where the overlapping portion is generated and the width B of one histogram. In this case, ΔB ′ may be used alone to evaluate the size of the overlapping relationship. In any case, by quantitatively evaluating the size of the intersection of the two histograms, it is possible to objectively determine the trace method from the evaluation result.

図9には、図5を用いて説明した方法の変形例が示されている。図9を用いて説明する手法は区間の決定において図5を用いて説明した方法と相違する。なお、図9において、図5に示した構成と同様の構成には同一符号を付してその説明を省略する。   FIG. 9 shows a modification of the method described with reference to FIG. The method described using FIG. 9 is different from the method described using FIG. 5 in determining the section. In FIG. 9, the same components as those shown in FIG.

図9において、規格化後の外側ヒストグラム100は、輝度108から輝度109の間に存在し、同じく規格化後の内側ヒストグラム102は、輝度104から輝度106の間に存在している。内側ヒストグラム102においては、対象組織の構造や高輝度ノイズ等に起因して、その高輝度側に低レベルの裾野部分102Aが生じており、それがかなり外側ヒストグラム100内にくい込んでいる。このような場合において、図5を用いて説明した手法をそのまま適用すると、重複部分が過大に見積られてしまうという問題が生じる。そこで、裾野部分102Aの全部又は一部を除外して有効な重複区間を定めるのが望ましい。   In FIG. 9, the normalized outer histogram 100 exists between the luminance 108 and the luminance 109, and the normalized inner histogram 102 exists between the luminance 104 and the luminance 106. In the inner histogram 102, due to the structure of the target tissue, high luminance noise, and the like, a low level skirt portion 102A is generated on the high luminance side, which is considerably difficult to insert in the outer histogram 100. In such a case, if the method described with reference to FIG. 5 is applied as it is, there arises a problem that an overlapping portion is excessively estimated. Therefore, it is desirable to define an effective overlapping section by excluding all or part of the base portion 102A.

具体的には、所定値(度数)C1が固定的にあるいは適応的に定められ、それよりも低いレベルをもった部分が区間判定対象から除外される。図9に示す例では、内側ヒストグラム102における、高輝度側におけるC1以下の部分(点102Bよりも右側の部分であって輝度107と輝度106との間の部分)が低レベル裾野部分とみなされ、内側ヒストグラムの全体区間A1から、区間A3を除いた区間A2が有効な区間とされ、内側ヒストグラム102における区間A2に属する主要部について面積S2が演算される。一方、外側ヒストグラム100についても、有効重複区間は輝度108から輝度107までの区間であり、外側ヒストグラム100における同区間内の面積ΔS1が演算される。そして、面積S2と面積ΔS1との比率をもって、上記同様に妥当なトレース方針が判定される。   Specifically, the predetermined value (frequency) C1 is fixedly or adaptively determined, and a portion having a lower level is excluded from the section determination targets. In the example shown in FIG. 9, the portion below C1 on the high luminance side (the portion on the right side of the point 102B and between the luminance 107 and the luminance 106) in the inner histogram 102 is regarded as the low level base portion. The section A2 excluding the section A3 from the entire section A1 of the inner histogram is regarded as an effective section, and the area S2 is calculated for the main part belonging to the section A2 in the inner histogram 102. On the other hand, also for the outer histogram 100, the effective overlapping section is a section from the luminance 108 to the luminance 107, and the area ΔS1 in the outer histogram 100 in the same section is calculated. Then, an appropriate trace policy is determined in the same manner as described above, with the ratio between the area S2 and the area ΔS1.

上記の所定値C1は固定値として定めることが可能である。それに代えて、内側ヒストグラム102に対する統計処理の演算結果(例えば分散、標準偏差、半値幅)に基づいて所定値C1を適用的に定めるようにしてもよい。例えば、ピークレベルCから一定割合だけ下がった地点として所定値C1を定めるようにしてもよい。これが図9においてDで示されている。なお、統計処理の演算結果から閾値レベルを定めるのではなく、面積割合等を基準として除外部分を定めることも可能である。つまり、一定割合の面積をもった部分を面積演算対象となる主要部として認定し、同時に、それ以外の面積をもった部分から外側ヒストグラム100についての有効重複区間を認定するようにしてもよい。   The predetermined value C1 can be determined as a fixed value. Instead, the predetermined value C <b> 1 may be appropriately determined based on the calculation result (for example, variance, standard deviation, half width) of the statistical processing for the inner histogram 102. For example, the predetermined value C1 may be determined as a point that is lower than the peak level C by a certain percentage. This is indicated by D in FIG. Instead of determining the threshold level from the calculation result of the statistical processing, it is possible to determine the excluded portion based on the area ratio or the like. That is, a portion having a certain percentage of area may be recognized as a main portion that is an area calculation target, and at the same time, an effective overlapping section for the outer histogram 100 may be recognized from a portion having other areas.

上記の説明では、内側ヒストグラム102の高域側における裾野部分について除外処理が適用されていたが、更に、内側ヒストグラムの低域側についても同様の除外処理を適用してもよい。また、外側ヒストグラム100における高輝度側の裾野部分に対して上記同様の除外処理を適用してもよい。更に、一方又は両方のヒストグラムについて一定レベル以下の部分を全部除外する処理を適用するようにしてもよい。その場合においても、裾野部分を除外する結果が得られるので、上記で指摘した裾野部分が輝度軸上に広がっていることに起因する問題を解消することができる。   In the above description, the exclusion process is applied to the bottom part of the inner histogram 102 on the high frequency side, but the same exclusion process may be applied to the low frequency side of the inner histogram. Further, the same exclusion process as described above may be applied to the bottom portion on the high luminance side in the outer histogram 100. Furthermore, a process of excluding all portions below a certain level for one or both histograms may be applied. Even in such a case, a result of excluding the skirt portion can be obtained, so that the problem caused by the above-described skirt portion spreading on the luminance axis can be solved.

本発明に係る超音波データ処理装置の好適な実施形態を示すブロック図である。1 is a block diagram showing a preferred embodiment of an ultrasonic data processing apparatus according to the present invention. 2つのROI(関心領域)の設定方法を説明するための図である。It is a figure for demonstrating the setting method of two ROI (region of interest). 外側ヒストグラムを示す図である。It is a figure which shows an outer side histogram. 内側ヒストグラムを示す図である。It is a figure which shows an inner side histogram. 2つのヒストグラムの重複関係を示す図である。It is a figure which shows the duplication relationship of two histograms. 他の例に係る2つのROIの設定を説明するための図である。It is a figure for demonstrating the setting of two ROI which concerns on another example. 図6に示した例における2つのヒストグラムの重複関係を示す図である。It is a figure which shows the duplication relationship of two histograms in the example shown in FIG. 重複関係の他の判定方法を説明するための図である。It is a figure for demonstrating the other determination method of duplication relationship. 図5に示した方法の変形例を説明するための図である。It is a figure for demonstrating the modification of the method shown in FIG.

符号の説明Explanation of symbols

10 3Dプローブ、12 送受信部、14 3Dメモリ、16 データ処理部、22 トレース方式判定ユニット、24,26 ヒストグラム作成部、28,30 規格化部、32 重複計測部、34 判定部、40 画像形成部、42 自動トレース部、44 体積/面積演算部。   10 3D probe, 12 transmission / reception unit, 14 3D memory, 16 data processing unit, 22 trace method determination unit, 24, 26 histogram creation unit, 28, 30 normalization unit, 32 overlap measurement unit, 34 determination unit, 40 image formation unit 42 Automatic trace section, 44 Volume / area calculation section.

Claims (10)

二次元又は三次元の生体内空間から取得された超音波データを処理する超音波データ処理装置において、
前記生体内空間に対応するデータ空間において、対象組織の外側に外側参照部を設定し、前記対象領域の内側に内側参照部を設定するための設定手段と、
前記外側参照部についての外側輝度分布及び前記内側参照部についての内側輝度分布の相互関係に基づいて、前記超音波データに対するトレース方式として、自動トレース処理又は手動トレース処理を判定する判定手段と、
を含むことを特徴とする超音波データ処理装置。
In an ultrasonic data processing apparatus for processing ultrasonic data acquired from a two-dimensional or three-dimensional in-vivo space,
In a data space corresponding to the in-vivo space, setting means for setting an outer reference portion outside a target tissue and setting an inner reference portion inside the target region;
Based on the correlation between the outer luminance distribution for the outer reference portion and the inner luminance distribution for the inner reference portion, as a tracing method for the ultrasonic data, a determination unit that determines automatic tracing processing or manual tracing processing;
An ultrasonic data processing apparatus comprising:
請求項1記載の装置において、
前記相互関係は、前記外側輝度分布及び前記内側輝度分布の重複度である、ことを特徴とする超音波データ処理装置。
The apparatus of claim 1.
The ultrasonic data processing apparatus according to claim 1, wherein the correlation is a degree of overlap between the outer luminance distribution and the inner luminance distribution.
請求項2記載の装置において、
前記外側参照部の参照により前記外側輝度分布を作成する外側輝度分布手段と、
前記内側参照部の参照により前記内側輝度分布を作成する内側輝度分布手段と、
前記外側輝度分布及び前記内側輝度分布の一方又は双方を規格化処理する規格化処理手段と、
を含み、
前記判定手段は、前記規格化処理後の外側輝度分布及び内側輝度分布の重複度に基づいて、自動トレース処理又は手動トレース処理を判定する、
ことを特徴とする超音波データ処理装置。
The apparatus of claim 2.
Outer luminance distribution means for creating the outer luminance distribution by referring to the outer reference section;
An inner luminance distribution means for creating the inner luminance distribution by referring to the inner reference portion;
Normalization processing means for normalizing one or both of the outer luminance distribution and the inner luminance distribution;
Including
The determination means determines automatic trace processing or manual trace processing based on the degree of overlap between the outer luminance distribution and the inner luminance distribution after the normalization processing.
An ultrasonic data processing apparatus.
請求項3記載の装置において、
前記判定手段は、
輝度軸上における前記規格化処理後の外側輝度分布及び内側輝度分布の重複区間を特定する手段と、
前記外側輝度分布の内で前記重複区間に属する部分を演算する手段と、
前記内側輝度分布の全体に対する前記重複区間に属する部分の割合に基づいて前記重複度を演算する手段と、
前記重複度と所定の基準値との比較結果に基づいて前記自動トレース処理又は前記手動トレース処理を判定する手段と、
を含むことを特徴とする超音波データ処理装置。
The apparatus of claim 3.
The determination means includes
Means for specifying an overlapping section of the outer luminance distribution and the inner luminance distribution after the normalization processing on the luminance axis;
Means for calculating a portion belonging to the overlapping section in the outer luminance distribution;
Means for calculating the degree of overlap based on a ratio of a portion belonging to the overlap section with respect to the entire inner luminance distribution;
Means for determining the automatic trace process or the manual trace process based on a comparison result between the degree of overlap and a predetermined reference value;
An ultrasonic data processing apparatus comprising:
請求項1記載の装置において、
前記設定手段は、
前記超音波データに基づいて前記対象組織の代表断面を表す断層画像を生成する手段と、
前記断層画像上において、前記外側参照部としての二次元閉ループのユーザー指定を受け付ける手段と、
前記断層画像上において、前記内側参照部としての二次元閉ループのユーザー指定を受け付ける手段と、
を含むことを特徴とする超音波データ処理装置。
The apparatus of claim 1.
The setting means includes
Means for generating a tomographic image representing a representative cross section of the target tissue based on the ultrasound data;
On the tomographic image, means for receiving a user designation of a two-dimensional closed loop as the outer reference unit;
On the tomographic image, means for receiving a user designation of a two-dimensional closed loop as the inner reference portion;
An ultrasonic data processing apparatus comprising:
請求項1記載の装置において、
前記設定手段は、
前記超音波データに基づいて前記対象組織についての複数の断面を表す複数の断層画像を生成する手段と、
前記複数の断層画像上において、前記外側参照部としての三次元球状面のユーザー指定を受け付ける手段と、
前記複数の断層画像上において、前記内側参照部としての三次元球状面のユーザー指定を受け付ける手段と、
を含むことを特徴とする超音波データ処理装置。
The apparatus of claim 1.
The setting means includes
Means for generating a plurality of tomographic images representing a plurality of cross sections of the target tissue based on the ultrasound data;
On the plurality of tomographic images, means for accepting user designation of a three-dimensional spherical surface as the outer reference portion;
On the plurality of tomographic images, means for accepting user designation of a three-dimensional spherical surface as the inner reference portion;
An ultrasonic data processing apparatus comprising:
二次元又は三次元の生体内空間から取得された超音波データを処理する超音波データ処理装置において、
前記生体内空間に対応するデータ空間において、対象組織の外側に外部参照部を設定し、前記対象組織の内側に内側参照部を設定するための設定手段と、
前記外側参照部についての外側輝度分布及び前記内側参照部についての内側輝度分布に基づいて、前記超音波データに対するトレース方式として、自動トレース処理又は手動トレース処理を判定する判定手段と、
を含み、
前記判定手段は、
前記外部参照部の参照により前記外側輝度分布を作成する外側輝度分布作成手段と、
前記内側参照部の参照により前記内側輝度分布を作成する内側輝度分布作成手段と、
輝度軸上における前記外側輝度分布と前記内側輝度分布の有効重複範囲を用いて求められる重複度に基づいて、自動トレース処理又は手動トレース処理を判定する処理方式判定手段と、
を含み、
前記処理方式判定手段は、前記外側輝度分布と前記内側輝度分布の内の少なくとも一方について裾野部分を除外することにより前記有効重複範囲を決定する、ことを特徴とする超音波データ処理装置。
In an ultrasonic data processing apparatus for processing ultrasonic data acquired from a two-dimensional or three-dimensional in-vivo space,
In a data space corresponding to the in-vivo space, a setting means for setting an external reference portion outside the target tissue and setting an inner reference portion inside the target tissue;
Based on the outer luminance distribution for the outer reference portion and the inner luminance distribution for the inner reference portion, as a tracing method for the ultrasonic data, a determination unit that determines automatic trace processing or manual trace processing;
Including
The determination means includes
An outer luminance distribution creating means for creating the outer luminance distribution by referring to the external reference unit;
An inner luminance distribution creating means for creating the inner luminance distribution by referring to the inner reference portion;
Processing method determination means for determining automatic trace processing or manual trace processing based on the degree of overlap obtained using the effective overlap range of the outer luminance distribution and the inner luminance distribution on the luminance axis;
Including
The ultrasonic data processing apparatus according to claim 1, wherein the processing method determination unit determines the effective overlapping range by excluding a tail portion of at least one of the outer luminance distribution and the inner luminance distribution.
請求項7記載の装置において、
前記裾野部分は、前記内側輝度分布における高輝度側の部分であって所定レベル以下の部分である、ことを特徴とする超音波データ処理装置。
The apparatus of claim 7.
The ultrasonic data processing apparatus according to claim 1, wherein the base portion is a portion on a high luminance side in the inner luminance distribution and a portion below a predetermined level.
請求項7記載の装置において、
前記裾野部分は、前記内側輝度分布における高輝度側の部分であって前記内側輝度分布のばらつき度合いを示す統計的指標に基づいて決定される部分である、ことを特徴とする超音波データ処理装置。
The apparatus of claim 7.
The base portion is a portion on the high luminance side in the inner luminance distribution, and is a portion determined based on a statistical index indicating a variation degree of the inner luminance distribution. .
請求項7乃至9のいずれか1項に記載の装置において、
前記処理方式判定手段は、
前記内側輝度部分の裾野部分を除外することにより前記有効重複区間を決定する手段と、
前記内側輝度分布における前記裾野部分を除外した後の主要部と、前記外側輝度分布における前記有効重複区間に属する部分と、の比率を前記重複度として演算する手段と、
前記比率を一定値と比較することにより、前記自動トレース処理又は前記手動トレース処理を判定する手段と、
を含むことを特徴とする超音波データ処理装置。
The device according to any one of claims 7 to 9,
The processing method determination means includes
Means for determining the effective overlap section by excluding the skirt portion of the inner luminance portion;
Means for calculating, as the degree of overlap, a ratio between a main part after excluding the skirt part in the inner luminance distribution and a part belonging to the effective overlapping section in the outer luminance distribution;
Means for determining the automatic trace process or the manual trace process by comparing the ratio with a constant value;
An ultrasonic data processing apparatus comprising:
JP2008230162A 2008-04-14 2008-09-08 Ultrasonic data processor Expired - Fee Related JP5222067B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008230162A JP5222067B2 (en) 2008-04-14 2008-09-08 Ultrasonic data processor

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2008104549 2008-04-14
JP2008104549 2008-04-14
JP2008230162A JP5222067B2 (en) 2008-04-14 2008-09-08 Ultrasonic data processor

Publications (2)

Publication Number Publication Date
JP2009273862A JP2009273862A (en) 2009-11-26
JP5222067B2 true JP5222067B2 (en) 2013-06-26

Family

ID=41439814

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008230162A Expired - Fee Related JP5222067B2 (en) 2008-04-14 2008-09-08 Ultrasonic data processor

Country Status (1)

Country Link
JP (1) JP5222067B2 (en)

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01254151A (en) * 1988-04-04 1989-10-11 Fujitsu Ltd Ultrasonic diagnosing device
JPH0496735A (en) * 1990-08-10 1992-03-30 Toshiba Corp Image processing device
JPH0732773B2 (en) * 1991-03-27 1995-04-12 アロカ株式会社 Ultrasonic image display device
JP4373698B2 (en) * 2003-04-25 2009-11-25 株式会社東芝 Ultrasonic diagnostic apparatus and ultrasonic diagnostic support program
JP4725461B2 (en) * 2006-08-29 2011-07-13 ソニー株式会社 Image signal processing apparatus, imaging apparatus, and image signal processing method

Also Published As

Publication number Publication date
JP2009273862A (en) 2009-11-26

Similar Documents

Publication Publication Date Title
US20210177373A1 (en) Ultrasound system with an artificial neural network for guided liver imaging
US11793483B2 (en) Target probe placement for lung ultrasound
US8582848B2 (en) System and method for detection of acoustic shadows and automatic assessment of image usability in 3D ultrasound images
JP4899837B2 (en) Ultrasound imaging system and method
US8204288B2 (en) Apparatus and method for processing an ultrasound spectrum image
CN109788939B (en) Method and system for enhancing visualization and selection of representative ultrasound images by automatically detecting B-lines and scoring images of ultrasound scans
US20190105013A1 (en) Ultrasound system and method for detecting lung sliding
JP7336443B2 (en) Ultrasound imaging system, device, method and storage medium
JP7010948B2 (en) Fetal ultrasound imaging
CN110072466B (en) Prenatal ultrasound imaging
US20200202551A1 (en) Fetal development monitoring
JP2020511250A (en) Volume rendered ultrasound image
JP2019528881A (en) Apparatus and method for detecting interventional tools
KR20210081243A (en) Methods and systems for automatic measurement of strains and strain-ratio calculation for sonoelastography
KR102017069B1 (en) Method for measuring fetal body and device for measuring fetal body using the same
JP5222067B2 (en) Ultrasonic data processor
JP5406484B2 (en) Ultrasonic data processor
JP4616465B2 (en) Region detection method and ultrasonic image processing apparatus
US20240119705A1 (en) Systems, methods, and apparatuses for identifying inhomogeneous liver fat
JP6286926B2 (en) Ultrasonic diagnostic apparatus, image processing method, and program
JP7299100B2 (en) ULTRASOUND DIAGNOSTIC DEVICE AND ULTRASOUND IMAGE PROCESSING METHOD
US20220277834A1 (en) Methods and systems for utilizing histogram views for improved visualization of three-dimensional (3d) medical images
CN115886876A (en) Fetal posture evaluation method, ultrasonic imaging method and ultrasonic imaging system
CN113040822A (en) Method for measuring endometrial peristalsis and device for measuring endometrial peristalsis
JP2024522310A (en) SYSTEM AND METHOD FOR PROVIDING USER GUIDANCE FOR ACQUIRING HEPATIC-RENAL INDEX MEASUREMENTS - Patent application

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110809

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130221

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130308

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

Free format text: PAYMENT UNTIL: 20160315

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

LAPS Cancellation because of no payment of annual fees