JP4768495B2

JP4768495B2 - Ultrasonic diagnostic equipment

Info

Publication number: JP4768495B2
Application number: JP2006102660A
Authority: JP
Inventors: 浩二廣田
Original assignee: Hitachi Aloka Medical Ltd
Current assignee: Hitachi Ltd
Priority date: 2006-04-04
Filing date: 2006-04-04
Publication date: 2011-09-07
Anticipated expiration: 2026-04-04
Also published as: JP2007275150A

Description

本発明は超音波診断装置に関し、特に三次元空間から得られたデータの重み付け処理に関する。 The present invention relates to an ultrasonic diagnostic apparatus, and more particularly to weighting processing of data obtained from a three-dimensional space.

超音波診断装置は、超音波の送受波によって得られたデータに基づいて超音波画像を形成する装置である。近時、生体内の三次元領域に対して超音波を送受波し、それにより得られたデータに基づいて三次元超音波画像を形成することが行われている。 An ultrasonic diagnostic apparatus is an apparatus that forms an ultrasonic image based on data obtained by transmission and reception of ultrasonic waves. Recently, ultrasonic waves are transmitted / received to / from a three-dimensional region in a living body, and a three-dimensional ultrasonic image is formed based on data obtained thereby.

特許文献１には、二次元の超音波画像から所定組織の輪郭を抽出する超音波診断装置が記載されている。かかる装置では、超音波画像上に関心領域が設定され、その形状を基準として関心領域内に段階的に重み付けパターンが設定されている。その重み付けパターンに従って画像に対する重み付け処理が施される。 Patent Document 1 describes an ultrasonic diagnostic apparatus that extracts a contour of a predetermined tissue from a two-dimensional ultrasonic image. In such an apparatus, a region of interest is set on an ultrasonic image, and a weighting pattern is set stepwise in the region of interest based on the shape thereof. A weighting process is performed on the image according to the weighting pattern.

特開平１１−２６７１２７号公報JP 11-267127 A

上記特許文献１に記載された手法によれば、関心領域内のデータ群に対して画一的な処理を行うのではなく、各データの位置に応じて適応的に重み付け処理を施した上で、画像処理を行える。しかし、特許文献１には三次元データに対する重み付け処理は記載されていない。また、特許文献１の手法では、複数の重み付けパターンからなるパターンセットを利用するので、複数の重み付けパターンの形状を個別的に保存、管理する必要があるので、それを三次元に拡張適用することは難しいという面を指摘できる。 According to the method described in the above-mentioned Patent Document 1, instead of performing uniform processing on the data group in the region of interest, after performing adaptive weighting processing according to the position of each data, Image processing can be performed. However, Patent Document 1 does not describe a weighting process for three-dimensional data. In addition, since the method of Patent Document 1 uses a pattern set composed of a plurality of weighting patterns, it is necessary to individually store and manage the shapes of the plurality of weighting patterns. Can point out that it is difficult.

本発明の目的は、三次元空間から得られたデータに対して簡便かつ適切に重み付け処理を施すことにある。 An object of the present invention is to perform weighting processing simply and appropriately on data obtained from a three-dimensional space.

本発明の目的は、三次元位置に応じて重み値を連続的に変動させるようにして重み付け処理の結果を自然なものにできるようにすることにある。 An object of the present invention is to make the weighting result natural by continuously changing the weight value according to the three-dimensional position.

本発明は、生体内の三次元空間に対して超音波を送受波する送受波手段と、前記超音波の送受波によって得られるデータ集合を構成する各データに対し、各データごとに特定される重み係数を用いて、前処理としての重み付け処理を施す重み付け処理手段と、前記重み付け処理後のデータ集合に対してデータ処理を実行するデータ処理手段と、を含み、前記重み付け処理手段は、前記三次元空間において立体形状を定義する数式を利用して、前記各データごとに前記重み係数を演算する重み係数演算手段と、前記各データに対して前記重み係数を用いた重み付け処理を実行する重み付け処理実行手段と、を含むことを特徴とする超音波診断装置に関する。 The present invention is specified for each data with respect to each data constituting a data set obtained by transmitting and receiving ultrasonic waves to and from a three-dimensional space in a living body and a data set obtained by the ultrasonic wave transmission and reception. Weighting processing means for performing weighting processing as preprocessing using a weighting factor, and data processing means for executing data processing on the data set after the weighting processing, wherein the weighting processing means includes the tertiary processing A weighting factor calculating means for calculating the weighting factor for each piece of data using a mathematical formula that defines a solid shape in the original space, and a weighting step for executing a weighting process using the weighting factor for each piece of data And an executing means.

上記構成によれば、三次元空間内で取得されたデータを処理する場合には、それに先立って前処理として重み付け処理が適用される。重み付け処理では、三次元空間において立体的形状を定義する数式（つまり関数）を用いて、データの座標から重み係数が演算される。関数演算を利用するので、複雑な重み付けパターンを保有しておくことによる問題を解消できる。勿論、関数演算の結果として直ちに重み係数が求められてもよいし、関数演算の結果に基づいて重み係数の区分を特定するようにしてもよいし、関数演算の結果に線形変換等を施して重み係数を特定するようにしてもよい。 According to the above configuration, when data acquired in a three-dimensional space is processed, a weighting process is applied as a preprocessing prior to that. In the weighting process, a weighting coefficient is calculated from the coordinates of data using a mathematical expression (that is, a function) that defines a three-dimensional shape in a three-dimensional space. Since the function operation is used, the problem caused by having a complicated weighting pattern can be solved. Of course, the weighting factor may be obtained immediately as a result of the function calculation, the classification of the weighting factor may be specified based on the result of the function calculation, or linear transformation or the like may be performed on the result of the function calculation. A weighting factor may be specified.

望ましくは、前記重み係数演算手段は、前記各データについて、送受波座標系における座標からデータ処理用の直交座標系における座標を演算する手段と、前記数式に対して前記各データについての直交座標系での座標を与える関数演算を実行することにより、前記重み係数を演算する手段と、を含む。この構成によれば、データの座標が送受波座標系（例えば極座標系）で表される場合でも、関数演算を直交座標系上で行ってその演算内容を簡易化できる。もちろん、すべて送受波座標系で処理するようにしてもよい。 Preferably, the weighting factor calculating means calculates means for calculating the coordinates in the orthogonal coordinate system for data processing from the coordinates in the transmission / reception coordinate system for each of the data, and the orthogonal coordinate system for the respective data with respect to the mathematical expression. And means for calculating the weighting coefficient by executing a function calculation that gives coordinates at. According to this configuration, even when the data coordinates are expressed in a transmission / reception coordinate system (for example, a polar coordinate system), it is possible to perform the function calculation on the orthogonal coordinate system and simplify the calculation content. Of course, all may be processed in the transmission / reception coordinate system.

望ましくは、前記数式に含まれる複数のパラメータ値をユーザーによって指定するためのユーザー指定手段を含む。この場合に、互いに交差関係にある複数の断層画像を表示し、それらの断層画像上で三次元形状あるいはそれについての基本的パラメータ（中心座標、半径等）を指定させるようにしてもよい。望ましくは、前記数式は楕円体又は球体を定義する数式である。生体の臓器（あるいは組織）の抽出や弁別を前提とする場合、多くの臓器が丸みをもって塊として存在していることから、上記のような関数を利用するのが合理的である。本発明に係る前処理は特に三次元輪郭抽出に先立って実行されるのが望ましい。 Preferably, a user specifying means for specifying a plurality of parameter values included in the mathematical expression by a user is included. In this case, a plurality of tomographic images that intersect each other may be displayed, and a three-dimensional shape or basic parameters (center coordinates, radius, etc.) about the three-dimensional shape may be designated on these tomographic images. Preferably, the mathematical formula is a mathematical formula that defines an ellipsoid or a sphere. When it is assumed that a living organ (or tissue) is extracted or discriminated, it is reasonable to use the above function because many organs exist as a lump with roundness. It is desirable that the pre-processing according to the present invention be performed prior to the three-dimensional contour extraction.

以上説明したように、本発明によれば、三次元空間から得られたデータに対して簡便かつ適切に重み付け処理を施せる。あるいは、三次元位置に応じて重み値を連続的に変動させるようにして重み付け処理の結果を自然なものにできる。 As described above, according to the present invention, weighting processing can be performed simply and appropriately on data obtained from a three-dimensional space. Alternatively, the weighting result can be made natural by continuously changing the weight value according to the three-dimensional position.

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

図１には、本発明に係る超音波診断装置の好適な実施形態が示されており、図１はその全体構成を示すブロック図である。 FIG. 1 shows a preferred embodiment of an ultrasonic diagnostic apparatus according to the present invention, and FIG. 1 is a block diagram showing the overall configuration thereof.

３Ｄプローブ１０は、生体内における三次元空間に対して超音波を送受波する送受波器である。３Ｄプローブ１０は、本実施形態において２Ｄアレイ振動子を有している。この２Ｄアレイ振動子は複数の振動素子を二次元配列してなるものである。この場合、超音波ビームは二次元的に電子走査される。もちろん、１Ｄアレイ振動子を機械的に走査することにより、三次元空間を形成するようにしてもよい。 The 3D probe 10 is a transducer that transmits and receives ultrasonic waves to and from a three-dimensional space in a living body. The 3D probe 10 has a 2D array transducer in this embodiment. This 2D array vibrator is formed by two-dimensionally arranging a plurality of vibration elements. In this case, the ultrasonic beam is electronically scanned two-dimensionally. Of course, a three-dimensional space may be formed by mechanically scanning the 1D array transducer.

図１には、三次元空間１２が概念図として模式的に示されている。三次元空間１２は、三次元エコーデータ取込空間であり、超音波ビーム１４を二次元走査することにより形成される。具体的には、超音波ビーム１４における深さ方向がｒ方向であり、超音波ビーム１４はθ方向に電子走査される。これによって走査面が形成される。走査面はφ方向に順次移動され、その結果として、三次元空間１２が形成される。三次元空間１２の形状としてはこれ以外にも各種のものがあげられる。送受信部１６は、送信ビームフォーマー及び受信ビームフォーマーとして機能する。送受信部１６から３Ｄプローブ１０における複数の振動素子に対して複数の送信信号が並列的に供給され、これによって３Ｄプローブ１０において送信ビームが形成される。一方、超音波の受波時においては、複数の振動素子から出力された複数の受信信号が送受信部１６に入力され、それらの受信信号に対する整相加算処理によって整相加算後の受信信号が得られる。すなわち、これによって電子的に受信ビームが形成される。送受信部１６は、各超音波ビームごとに受信信号としてのビームデータを出力している。 In FIG. 1, a three-dimensional space 12 is schematically shown as a conceptual diagram. The three-dimensional space 12 is a three-dimensional echo data capturing space, and is formed by two-dimensional scanning with the ultrasonic beam 14. Specifically, the depth direction in the ultrasonic beam 14 is the r direction, and the ultrasonic beam 14 is electronically scanned in the θ direction. As a result, a scanning plane is formed. The scanning plane is sequentially moved in the φ direction, and as a result, a three-dimensional space 12 is formed. Various other shapes can be used as the shape of the three-dimensional space 12. The transmission / reception unit 16 functions as a transmission beam former and a reception beam former. A plurality of transmission signals are supplied in parallel to the plurality of vibration elements in the 3D probe 10 from the transmission / reception unit 16, whereby a transmission beam is formed in the 3D probe 10. On the other hand, when receiving ultrasonic waves, a plurality of reception signals output from a plurality of vibration elements are input to the transmission / reception unit 16, and a reception signal after phasing addition is obtained by phasing addition processing on the reception signals. It is done. That is, a reception beam is electronically formed by this. The transmission / reception unit 16 outputs beam data as a reception signal for each ultrasonic beam.

ビームデータ処理部１８は、検波、対数圧縮等の必要なビームデータ処理を実行するモジュールである。このビームデータ処理部１８は必要に応じて設けられる。３Ｄメモリ２０は、上記の三次元空間１２内において取得されたデータ群をいったん格納するメモリであり、３Ｄメモリ２０は上記の三次元空間１２に相当する記憶空間を有している。３Ｄメモリ２０に格納された各データは必要に応じて読み出されて、以下に説明する断層画像形成部２２、座標変換部２６及び三次元重み付け処理部３０へ出力される。 The beam data processing unit 18 is a module that executes necessary beam data processing such as detection and logarithmic compression. The beam data processing unit 18 is provided as necessary. The 3D memory 20 is a memory that temporarily stores a data group acquired in the three-dimensional space 12. The 3D memory 20 has a storage space corresponding to the three-dimensional space 12. Each data stored in the 3D memory 20 is read out as necessary, and is output to a tomographic image forming unit 22, a coordinate conversion unit 26, and a three-dimensional weighting processing unit 30 described below.

断層画像形成部２２は、いわゆるトリプレーン画像を形成するために、三次元空間内においてユーザーにより任意に設定された複数の切断面に相当する複数の断層画像を形成する。その場合においては、３Ｄメモリ２０から必要なデータが読み出されている。断層画像形成部２２は座標変換機能及び補間機能などを有している。断層画像形成部２２がいわゆる任意切断面に対応する任意断層画像を形成するものであってもよい。断層画像形成部２２から出力される断層画像の画像データは表示処理部２４へ出力される。表示処理部２４は、後に説明するように複数の画像を同一画面上に表示する画像合成機能等を有している。 The tomographic image forming unit 22 forms a plurality of tomographic images corresponding to a plurality of cut planes arbitrarily set by the user in the three-dimensional space in order to form a so-called triplane image. In that case, necessary data is read from the 3D memory 20. The tomographic image forming unit 22 has a coordinate conversion function and an interpolation function. The tomographic image forming unit 22 may form an arbitrary tomographic image corresponding to a so-called arbitrary cut surface. The image data of the tomographic image output from the tomographic image forming unit 22 is output to the display processing unit 24. The display processing unit 24 has an image composition function for displaying a plurality of images on the same screen, as will be described later.

座標変換部２６は、三次元空間１２において取得されたデータ群に基づいて三次元画像を形成する前に送受波座標系から直交座標系への座標変換を実行している。すなわち上述したようにデータの取込においては、ｒθφ座標系で各データが特定されており、その一方において、画像処理あるいは三次元画像の形成はＸＹＺ直交座標系において行われており、座標変換部２６は各データに対して送受波座標系から直交座標系への座標変換を実行している。また、その場合において必要に応じて補間処理が実行される。座標変換後のデータは合成部２８へ出力される。ちなみに、合成部２８は必要に応じて３Ｄメモリを有する。 The coordinate conversion unit 26 performs coordinate conversion from the transmission / reception coordinate system to the orthogonal coordinate system before forming a three-dimensional image based on the data group acquired in the three-dimensional space 12. In other words, as described above, in data acquisition, each data is specified in the rθφ coordinate system, and on the other hand, image processing or three-dimensional image formation is performed in the XYZ orthogonal coordinate system, and the coordinate conversion unit 26 performs coordinate conversion from the transmission / reception coordinate system to the orthogonal coordinate system for each data. In that case, an interpolation process is executed as necessary. The data after the coordinate conversion is output to the combining unit 28. Incidentally, the synthesis unit 28 includes a 3D memory as necessary.

符号２１はデータ処理部を表しており、このデータ処理部は、以下に説明する三次元重み付け処理部３０等の各構成を具備している。三次元重み付け処理部３０は、入力されたデータに対してその三次元空間内における位置に応じた重み係数（重み値）を特定し、当該データに対してその重み係数を乗算する演算を実行している。すなわち、各データごとに重み付け処理を行っている。その具体的な演算方法については後に図２を用いて説明するが、本実施形態においては楕円体を表す関数を用いてデータの座標から関数演算の実行により重み係数を瞬時に演算している。 Reference numeral 21 represents a data processing unit, and this data processing unit includes various components such as a three-dimensional weighting processing unit 30 described below. The three-dimensional weighting processing unit 30 specifies a weighting factor (weight value) corresponding to the position in the three-dimensional space for the input data, and executes an operation for multiplying the data by the weighting factor. ing. That is, a weighting process is performed for each data. The specific calculation method will be described later with reference to FIG. 2. In the present embodiment, the weight coefficient is instantaneously calculated by executing the function calculation from the data coordinates using a function representing an ellipsoid.

二値化処理部３２は、重み付け処理後のデータに対して所定の閾値を用いた二値化処理を実行する。この場合において、例えば対象となる臓器が心臓であれば、心筋と心腔とを弁別する閾値が設定され、例えば心腔のみを抽出するように二値化処理が実行される。この場合においては必要に応じていわゆる反転二値化処理が実行される。座標変換部３４は、上記の座標変換部２６と同様に、送受波座標系から直交座標系への座標変換を実行する。この場合においては必要に応じて補間処理が適用される。ちなみに、座標変換は、任意の箇所において実行することが可能であり、例えば３Ｄメモリ２０の前段において座標変換を実行するようにしてもよいし、あるいは合成部２８の後段において座標変換を実行するようにしてもよい。 The binarization processing unit 32 executes binarization processing using a predetermined threshold for the data after the weighting processing. In this case, for example, if the target organ is the heart, a threshold value for discriminating the myocardium from the heart chamber is set, and for example, the binarization processing is executed so as to extract only the heart chamber. In this case, so-called inversion binarization processing is executed as necessary. Similar to the coordinate conversion unit 26 described above, the coordinate conversion unit 34 performs coordinate conversion from the transmission / reception coordinate system to the orthogonal coordinate system. In this case, interpolation processing is applied as necessary. Incidentally, the coordinate transformation can be executed at an arbitrary place. For example, the coordinate transformation may be executed in the former stage of the 3D memory 20, or the coordinate transformation may be executed in the latter stage of the synthesizing unit 28. It may be.

境界抽出部３６は、座標変換後の二値化データに基づいて、組織の境界を抽出する処理を実行する。例えば心筋と心腔との間の面すなわち心内膜を抽出する処理を実行する。これによって心内膜を表すデータが合成部２８へ出力されることになる。 The boundary extraction unit 36 executes a process of extracting the tissue boundary based on the binarized data after the coordinate conversion. For example, a process of extracting a plane between the myocardium and the heart chamber, that is, the endocardium is executed. As a result, data representing the endocardium is output to the synthesis unit 28.

合成部２８は、三次元空間内におけるデータ群を背景としてそれに対して心内膜を表すデータ等を合成することにより心臓を表すデータ群を生成する。この場合に心内膜のみが輝度強調されたデータ群が得られる。ただし、本実施形態における重み付け処理は、上述したような境界抽出処理以外のデータ処理に対しても適用でき、すなわちデータ処理に先だって前処理として特定の三次元データを強調する場合において利用することが可能である。 The synthesizing unit 28 generates a data group representing the heart by synthesizing data representing the endocardium with the data group in the three-dimensional space as a background. In this case, a data group in which only the endocardium is brightness enhanced is obtained. However, the weighting process in the present embodiment can also be applied to data processing other than the boundary extraction process as described above, that is, it can be used when emphasizing specific three-dimensional data as preprocessing prior to data processing. Is possible.

レンダリング部３８は、合成後のデータ群に対して所定の視点を設定した上で、その視点から出る複数の視線（レイ）を設定し、各レイごとに公知のボリュームレンダリング等を実行することにより三次元超音波画像を形成する。この場合においては、必要に応じてカラー演算などを適用するようにしてもよい。本実施形態においては、心内膜が輝度強調された心臓の三次元画像がいわゆるボリュームレンダリング画像として形成されている。この場合において心内膜のみに着色処理を施すようにしてもよい。 The rendering unit 38 sets a predetermined viewpoint for the combined data group, sets a plurality of lines of sight (rays) from the viewpoint, and executes known volume rendering or the like for each ray. A three-dimensional ultrasound image is formed. In this case, color calculation or the like may be applied as necessary. In the present embodiment, a three-dimensional image of the heart in which the endocardium is brightness enhanced is formed as a so-called volume rendering image. In this case, only the endocardium may be colored.

表示処理部２４は、超音波三次元画像と上述した１又は複数の断層画像とを並べて表示画像を構成し、その画像データを表示部４４に出力する。例えば、重み付け関数をユーザーによって定義する際においては、表示部４４にトリプレーン画像が表示され、そのトリプレーン画像を用いて重み付け関数として楕円体等が定義され、そのようなユーザー入力にしたがって上述した重み付け処理が実行される。 The display processing unit 24 configures a display image by arranging the ultrasonic three-dimensional image and the one or more tomographic images described above, and outputs the image data to the display unit 44. For example, when the weighting function is defined by the user, a triplane image is displayed on the display unit 44, and an ellipsoid or the like is defined as the weighting function using the triplane image. A weighting process is executed.

制御部４０は、プログラム動作するＣＰＵなどによって構成され、制御部４０は図１に示される各構成の動作制御を行っている。特に、本実施形態において、制御部４０は三次元重み付け処理部３０に対して必要なパラメータ値を渡している。入力部４２は例えば操作パネルによって構成され、その入力部４２を用いてユーザーが重み付け関数を自在に定義することが可能である。本実施形態においては後に説明するように楕円体関数として重み付け関数が定義されており、特に楕円体における３つの軸上における各半径と、その楕円体の中心座標とがユーザーによって指定される。楕円体関数に代えて例えば球体を表す関数を定義するようにしてもよい。 The control unit 40 is configured by a CPU or the like that performs a program operation, and the control unit 40 performs operation control of each configuration illustrated in FIG. In particular, in the present embodiment, the control unit 40 passes necessary parameter values to the three-dimensional weighting processing unit 30. The input unit 42 is configured by, for example, an operation panel, and the user can freely define a weighting function using the input unit 42. In this embodiment, as will be described later, a weighting function is defined as an ellipsoid function. In particular, each radius on three axes of the ellipsoid and the center coordinates of the ellipsoid are designated by the user. For example, a function representing a sphere may be defined instead of the ellipsoid function.

図２には、図１に示した三次元重み付け処理部３０の構成例が概念的に示されている。三次元重み付け処理部３０における具体的な処理のプロセスは例えばソフトウエアの機能によって実現することが可能であり、あるいは、デジタルシグナルプロセッサにおいて行わせるようにしてもよい。重み付け器５０は、今処理対象となっている注目データｄに対して重み係数αを乗算する乗算器として構成され、その処理結果として重み付け後の注目データｄ’が得られている。 FIG. 2 conceptually shows a configuration example of the three-dimensional weighting processing unit 30 shown in FIG. A specific processing process in the three-dimensional weighting processing unit 30 can be realized by a function of software, for example, or may be performed by a digital signal processor. The weighting unit 50 is configured as a multiplier that multiplies the attention data d that is currently processed by a weighting coefficient α, and the weighted attention data d ′ is obtained as a result of the processing.

重み係数αは重み演算器５４によって生成され、その重み演算器５４の前段には座標演算器５２が設けられている。座標演算器５２は超音波ビームの走査条件にしたがって、各注目データごとにそれを特定する座標を演算するモジュールであり、本実施形態においては、表示座標系における座標データか直交座標系における座標データが生成されている。例えば、座標変換器５２には超音波ビーム上における深さｒ，各走査面における初期ビーム角度θ₀，各走査面におけるビーム間ピッチ（ビーム間角度）Δθ，各走査面におけるビーム本数ｍ，初期走査面の角度φ₀，走査面間のピッチ（角度）Δφ，走査面数ｎ等が入力され、座標演算器５２はそれらの情報にしたがって、ＸＹＺ座標系すなわち直交座標系における注目データの座標（ｘ₁，ｙ₁，ｚ₁）を演算する。 The weight coefficient α is generated by the weight calculator 54, and a coordinate calculator 52 is provided in the preceding stage of the weight calculator 54. The coordinate calculator 52 is a module that calculates coordinates for specifying each piece of attention data according to the scanning condition of the ultrasonic beam. In this embodiment, the coordinate data in the display coordinate system or the coordinate data in the orthogonal coordinate system is used. Has been generated. For example, the coordinate converter 52 includes a depth r on the ultrasonic beam, an initial beam angle θ ₀ on each scanning plane, an inter-beam pitch (inter-beam angle) Δθ on each scanning plane, the number of beams m on each scanning plane, and an initial value. The angle φ _{0 of} the scanning plane, the pitch (angle) Δφ between scanning planes, the number n of scanning planes, etc. are input, and the coordinate calculator 52 follows the information to coordinate the data of interest ( x ₁ , y ₁ , z ₁ ) are calculated.

一方、重み演算器５４には、ユーザーによって指定された楕円体関数における中心座標（ｘ₀，ｙ₀，ｚ₀）が入力されており、更に当該楕円体における各軸上の半径Ａ，Ｂ，Ｃが入力されている。それらのパラメータ値はユーザーによって指定されたものである。 On the other hand, the center point (x ₀ , y ₀ , z ₀ ) in the ellipsoid function designated by the user is input to the weight calculator 54, and the radii A, B, C is entered. Those parameter values are specified by the user.

重み演算器５４は、基本的には楕円体を表す関数（Ｘ²／Ａ²＋Ｙ²／Ｂ²＋Ｚ²／Ｃ²）を利用してその関数値として重み係数αを求めるものである。ここで、ＸはＸ軸上において中心座標と注目データの座標との差分によって定義され、これはＹ及びＺについても同様である。本実施形態では、楕円体の関数が直交座標系上において定義されており、関数の形式は周知のように非常にシンプルなものである。一般に、上記の（Ｘ²／Ａ²＋Ｙ²／Ｂ²＋Ｚ²／Ｃ²）を左辺とし、その関数の右辺を１とした場合、それはユーザーによって定義された楕円体関数の三次元表面形状を定義するものである。本実施形態においては、仮に注目データがその表面形状に一致した場合には重みとして１が求まり、その表面形状よりも内側に注目データが存在している場合には中心座標からの距離をパラメータとして当該注目データについて右辺の数値として重み係数αが求められている。すなわち、一般的には楕円体関数はその表面形状を定義し、あるいはその内外判別に利用される関数であるが、本実施形態においては楕円体を定義する関数をそのまま重み値を求めるための関数として利用する点に特徴を有している。 The weight calculator 54 basically obtains a weight coefficient α as a function value using a function (X ² / A ² + Y ² / B ² + Z ² / C ² ) representing an ellipsoid. Here, X is defined by the difference between the center coordinates and the coordinates of the data of interest on the X axis, and the same applies to Y and Z. In this embodiment, the function of the ellipsoid is defined on the orthogonal coordinate system, and the format of the function is very simple as is well known. In general, when (X ² / A ² + Y ² / B ² + Z ² / C ² ) is set as the left side and the right side of the function is set as 1, it represents the three-dimensional surface shape of the ellipsoid function defined by the user. To define. In this embodiment, if the attention data matches the surface shape, 1 is obtained as the weight, and if the attention data exists inside the surface shape, the distance from the center coordinate is used as a parameter. A weighting factor α is obtained as a numerical value on the right side of the attention data. That is, in general, the ellipsoid function is a function that defines the surface shape or is used for the inside / outside discrimination, but in this embodiment, the function that defines the ellipsoid is a function for obtaining the weight value as it is. It has the feature in using as.

もちろん、求められたαをそのまま重み係数として利用してもよいし、αの値が属する区間ごとに別途重み値を定めておいて、属する区間に対応する重み値を重み係数として利用するようにしてもよい。あるいは、求められたαに対して線形変換等の他の関数演算を適用し、これによって一意に重み値を求めるようにしてもよい。いずれにしても、三次元的に複数の重み係数パターンを保有しておくと、その保有のための記憶容量が非常に増大することになるが、本実施形態においては、三次元形状を定義する関数を基礎として各データの座標から関数演算によって合理的な重み係数を迅速かつ容易に求められるという利点がある。ちなみに、図２において、符号５５は、上記関数式における右辺としてある値αが求められる場合における表面形態を表している。 Of course, the obtained α may be used as it is as a weighting factor, or a separate weighting value is set for each section to which the value of α belongs, and the weighting value corresponding to the section to which the value belongs is used as the weighting coefficient. May be. Alternatively, other function operations such as linear transformation may be applied to the obtained α, thereby uniquely obtaining the weight value. In any case, if a plurality of weight coefficient patterns are held in three dimensions, the storage capacity for holding them will greatly increase. In this embodiment, a three-dimensional shape is defined. There is an advantage that a rational weight coefficient can be obtained quickly and easily by function calculation from the coordinates of each data based on the function. Incidentally, in FIG. 2, reference numeral 55 represents a surface form in the case where a value α is obtained as the right side in the above functional expression.

図３には、組織との関係における重み付け関数が示されている。ここでは、説明の都合上、二次元断面として重み付け関数が表されている。符号１００は例えば心筋の断面を表しており、符号１００Ｂは心外膜であり、符号１００Ａは心内膜である。符号１０２は心腔を表している。その心腔１０２には例えば血液が満たされており、それは左室に相当する。本実施形態においては、このような組織との関係において、ユーザーによって中心座標Ｏが指定され、また直交座標系における各軸上の半径が定義される。そのような定義にしたがって、各αの値をマッピングするならば、図１においてＳ０〜Ｓ６で示されるような多重の三次元形状を想定できる。それぞれの表面上においては同じ重み係数αが与えられることになる。中心座標からの距離に応じて重み係数は連続的に変化することになり、ここで、例えばＳ０で表される最も外側の形状が上記関数式における右辺において１を表す場合に対応する。同様に、Ｓ１〜Ｓ６の各形状においてはそれぞれの形状に応じた重み係数が与えられることになる。また、関数式の右辺の値がどの区間（立体的区間）に属するかに基づいてその区間に対応する重み係数を求めるようにしてもよい。 FIG. 3 shows a weighting function in relation to the tissue. Here, for convenience of explanation, the weighting function is represented as a two-dimensional section. Reference numeral 100 represents, for example, a cross section of the myocardium, reference numeral 100B represents the epicardium, and reference numeral 100A represents the endocardium. Reference numeral 102 represents a heart chamber. The heart chamber 102 is filled with blood, for example, which corresponds to the left ventricle. In this embodiment, the center coordinate O is designated by the user in such a relationship with the tissue, and the radius on each axis in the orthogonal coordinate system is defined. If the values of each α are mapped according to such a definition, a multiple three-dimensional shape as shown by S0 to S6 in FIG. 1 can be assumed. The same weighting factor α is given on each surface. The weighting coefficient changes continuously according to the distance from the center coordinate, and corresponds to the case where the outermost shape represented by S0, for example, represents 1 on the right side in the above functional expression. Similarly, a weighting coefficient corresponding to each shape is given to each shape of S1 to S6. Further, based on which section (stereoscopic section) the value of the right side of the function expression belongs to, a weighting coefficient corresponding to the section may be obtained.

図４には、図３に示したＲ−Ｒ’における重み係数の変化が示されている。三次元形状を表す関数として各種のものを使用できるが、例えば符号１０４で示されるように、最も外側の位置から中心に向かうにしたがって徐々に重み係数の値を減少させるようにしてもよいし、また符号１０６で示されるように非線形にその減少が行われるようにしてもよい。あるいは、符号１０８で示されるように、各立体的区間ごとに段階的に重み係数の値が小さくなるようにしてもよい。あるいは、ユーザーにより設定された基本となる関心領域が例えばここではＳ３である場合、その内側及び外側に離れるにしたがって徐々に重み係数の値が小さくなるようにしてもよい。この場合において、例えば符号１０４で示される関数に対して一定の場合分けを行った上で一次変換を行うことにより、上記のような符号１１０で示される特性を容易に生成することが可能である。 FIG. 4 shows changes in the weighting factor in R-R ′ shown in FIG. 3. Various functions can be used as a function representing the three-dimensional shape. For example, as indicated by reference numeral 104, the value of the weighting factor may be gradually decreased from the outermost position toward the center. Further, as indicated by reference numeral 106, the reduction may be performed non-linearly. Alternatively, as indicated by reference numeral 108, the value of the weighting factor may be decreased step by step for each three-dimensional section. Alternatively, when the basic region of interest set by the user is, for example, S3 here, the value of the weighting factor may gradually decrease as the distance from the inner region to the outer region increases. In this case, it is possible to easily generate the characteristic indicated by reference numeral 110 as described above, for example, by performing a primary transformation after performing a predetermined case division on the function indicated by reference numeral 104. .

図５に示されるように、（Ａ）で示される元のデータ列１１２に対して、符号１１４で示されるような重み係数の特性を作用させると、図５における（Ｂ）で示すように、エッジに相当する部分が符号１３０で示されるように持ち上がり、その結果として重み付け後のデータ列１１２Ａを得られる。したがって、例えば閾値１１６として心筋と心腔を弁別する値を設定しておけば、心腔と心筋を明瞭に区別して、心内膜の画像を明確に表現することが可能となる。 As shown in FIG. 5, when the characteristic of the weight coefficient as shown by reference numeral 114 is applied to the original data string 112 shown in (A), as shown in (B) in FIG. 5, A portion corresponding to the edge is lifted as indicated by reference numeral 130, and as a result, a weighted data string 112A is obtained. Therefore, for example, if a value for discriminating the myocardium from the heart chamber is set as the threshold value 116, it is possible to clearly distinguish the heart chamber from the heart muscle and express the endocardial image clearly.

上記実施形態においては、三次元形状を定義するシンプルな関数を基礎として重み係数の値をその都度発生させるようにしたので、複雑なパターンを保有しておくための記憶容量を削減でき、また迅速に重み係数の演算を行って、三次元画像のリアルタイム処理を実現できるという利点がある。特に、楕円体を表す関数あるいは球体を表す関数を重み係数を生成するための基礎となる関数として利用したので、新しい重み値発生方式を実現できるという利点がある。 In the above embodiment, the value of the weighting coefficient is generated each time on the basis of a simple function that defines a three-dimensional shape, so that the storage capacity for holding a complex pattern can be reduced, and more quickly There is an advantage that real-time processing of a three-dimensional image can be realized by calculating a weighting factor. In particular, since a function representing an ellipsoid or a function representing a sphere is used as a function serving as a basis for generating a weighting coefficient, there is an advantage that a new weight value generation method can be realized.

本発明に係る超音波診断装置の全体構成を示すブロック図である。1 is a block diagram showing an overall configuration of an ultrasonic diagnostic apparatus according to the present invention. 図１に示す三次元重み付け処理部の具体的な構成例を示す概念図である。It is a conceptual diagram which shows the specific structural example of the three-dimensional weighting process part shown in FIG. 組織との関係における重み関数の変化を示す概念図である。It is a conceptual diagram which shows the change of the weight function in the relationship with an structure | tissue. 各種の重み関数を示す図である。It is a figure which shows various weight functions. 輪郭強調処理の内容を示す図である。It is a figure which shows the content of the outline emphasis process.

Explanation of symbols

１０３Ｄプローブ、２０３Ｄメモリ、２１データ処理部、２８合成部、３０三次元重み付け処理部、３２二値化処理部、３４座標変換部、３６境界抽出部、５０重み付け器、５２座標演算器、５４重み演算器。 10 3D probe, 20 3D memory, 21 data processing unit, 28 synthesis unit, 30 three-dimensional weighting processing unit, 32 binarization processing unit, 34 coordinate conversion unit, 36 boundary extraction unit, 50 weighting unit, 52 coordinate computing unit, 54 Weight calculator.

Claims

Transmitting and receiving means for transmitting and receiving ultrasonic waves to and from a three-dimensional space in the living body;
Weighting processing means for performing weighting processing as preprocessing, using a weighting coefficient specified for each data, for each data constituting a data set obtained by transmission and reception of the ultrasonic waves,
Data processing means for executing data processing on the data set after the weighting processing;
Including
The weighting processing means includes
A weighting factor calculating means for calculating the weighting factor for each of the data, using a mathematical formula that defines a three-dimensional shape in the three-dimensional space;
Weighting processing execution means for executing weighting processing using the weighting coefficient for each data;
An ultrasonic diagnostic apparatus comprising:

The apparatus of claim 1.
The weighting factor calculating means includes
For each of the data, means for calculating coordinates in the orthogonal coordinate system for data processing from coordinates in the transmission / reception coordinate system;
Means for calculating the weighting factor by performing a function calculation that gives coordinates in an orthogonal coordinate system for each data with respect to the mathematical formula;
An ultrasonic diagnostic apparatus comprising:

The apparatus of claim 1.
An ultrasonic diagnostic apparatus comprising user specifying means for specifying a plurality of parameter values included in the mathematical expression by a user.

The apparatus of claim 1.
The ultrasonic diagnostic apparatus, wherein the mathematical expression is a mathematical expression defining an ellipsoid or a sphere.