JP7213491B2 - Ultrasound Image Construction Method, Ultrasound Image Construction Apparatus, Ultrasound Image Construction Program - Google Patents

Ultrasound Image Construction Method, Ultrasound Image Construction Apparatus, Ultrasound Image Construction Program Download PDF

Info

Publication number
JP7213491B2
JP7213491B2 JP2019095229A JP2019095229A JP7213491B2 JP 7213491 B2 JP7213491 B2 JP 7213491B2 JP 2019095229 A JP2019095229 A JP 2019095229A JP 2019095229 A JP2019095229 A JP 2019095229A JP 7213491 B2 JP7213491 B2 JP 7213491B2
Authority
JP
Japan
Prior art keywords
acoustic impedance
ultrasonic
specific
value
measurement object
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.)
Active
Application number
JP2019095229A
Other languages
Japanese (ja)
Other versions
JP2020190454A (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.)
Honda Electronics Co Ltd
Toyohashi University of Technology NUC
Original Assignee
Honda Electronics Co Ltd
Toyohashi University of Technology NUC
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 Honda Electronics Co Ltd, Toyohashi University of Technology NUC filed Critical Honda Electronics Co Ltd
Priority to JP2019095229A priority Critical patent/JP7213491B2/en
Publication of JP2020190454A publication Critical patent/JP2020190454A/en
Application granted granted Critical
Publication of JP7213491B2 publication Critical patent/JP7213491B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

本発明は、超音波を利用して得た情報に基づき、皮膚等に代表される生体の軟組織や培養細胞などの画像を構築する方法、装置及びそのためのプログラムに関するものである。 The present invention relates to a method, an apparatus, and a program for constructing an image of a living body's soft tissue, such as skin, or cultured cells, based on information obtained using ultrasound.

超音波Bモードエコー像は、医療分野において広く一般的に使われている方法であり、従来このような像を得るための装置が多数提案されている(特許文献1を参照)。簡単に説明すると、超音波Bモードエコー像とは、物体に入射した超音波が反射して返って来る時の反射信号列を画像化したものである。超音波が散乱せずに真っ直ぐ進んだと仮定する場合、電気信号と同じく、進んだ先の抵抗値(固有音響インピーダンス)の違いによって反射が生じることになる。それゆえ、固有音響インピーダンスの分布が分かれば、どのような反射信号列が返って来るかを推測することが可能となる。即ち、音響物性分布が既知であれば、どのようなBモード画像が観察されるかを推測することができる。また、その逆もしかりである。 An ultrasonic B-mode echo image is a method widely used in the medical field, and many apparatuses for obtaining such an image have conventionally been proposed (see Patent Document 1). Briefly, an ultrasonic B-mode echo image is an image of a reflected signal train when ultrasonic waves incident on an object are reflected and returned. Assuming that the ultrasonic wave travels straight without scattering, reflection occurs due to the difference in the resistance value (specific acoustic impedance) where it travels, just like an electric signal. Therefore, if the distribution of the characteristic acoustic impedance is known, it becomes possible to estimate what kind of reflected signal train is returned. That is, if the acoustic physical property distribution is known, it is possible to estimate what kind of B-mode image will be observed. And vice versa.

しかしながら、生体組織などの不均一で厚み(深さ)があるターゲットについては、入射して返って来る反射波形は、ターゲットに入射した超音波が散乱、吸収を経て様々な進み方をした結果と、多重反射をした結果とを反映したものとなる。このため、反射波形を固有音響インピーダンス等の音響物性に変換することは困難であると考えられており、この方法は従来検討されてこなかった。さらに、超音波Bモードエコー像は、生体組織内部での超音波の多重反射に起因したスペックルノイズによって乱れた画像となりやすいため、内部構造を高い精度で表示するのには不向きであるという問題もあった。ゆえに、従来装置では音響フィルタを入れるなどの対策が必要となり、構成が複雑化する等の問題があった。 However, for a non-uniform and thick (depth) target such as a biological tissue, the reflected waveform that is incident and returned is the result of the ultrasonic waves that have entered the target traveling in various ways after being scattered and absorbed. , and the result of multiple reflection. For this reason, it is considered difficult to convert the reflected waveform into an acoustic physical property such as a specific acoustic impedance, and this method has not been studied in the past. Furthermore, since the ultrasonic B-mode echo image tends to be disturbed by speckle noise caused by multiple reflections of ultrasonic waves inside the living tissue, it is not suitable for displaying the internal structure with high accuracy. There was also Therefore, in the conventional apparatus, measures such as inserting an acoustic filter are required, which causes problems such as complicating the configuration.

また、超音波Bモードエコー像を表示する通常の超音波診断装置では、皮膚等の生体組織内部の層情報は一応得られるものの、得られる画像は固有音響インピーダンスの異なる層同士の界面からの反射像である。よって、このような反射像では生体組織の内部構造、具体的には生体組織内部の固有音響インピーダンスの違いを把握するのには不十分であった。つまり、従来技術で得られる反射像は、層の界面の位置がどこにあるか等については感覚的に理解しやすいものである反面、界面と界面とで囲まれた中間の領域の固有音響インピーダンスがどのようになっているかが感覚的に理解しにくいものであった。ゆえに、超音波を利用して得た情報に基づいて、層構造を有する非常に薄い測定対象物の超音波断層像を、感覚的に層構造が理解しやすい態様にて構築することが望まれていた。 In addition, with a normal ultrasonic diagnostic apparatus that displays an ultrasonic B-mode echo image, although layer information inside a living tissue such as skin can be obtained to some extent, the obtained image is the reflection from the interface between layers with different specific acoustic impedances. is a statue. Therefore, such a reflected image is insufficient to grasp the internal structure of the living tissue, specifically, the difference in the characteristic acoustic impedance inside the living tissue. In other words, in the reflected image obtained by the conventional technique, it is intuitively easy to understand the position of the interface between the layers, but on the other hand, the characteristic acoustic impedance of the intermediate region surrounded by the interfaces is It was difficult to intuitively understand how things were going. Therefore, it is desired to construct an ultrasonic tomographic image of a very thin measurement object having a layered structure in a manner that facilitates intuitive understanding of the layered structure, based on information obtained using ultrasonic waves. was

このような事情に鑑みて本願発明者らは、改良された超音波画像構築装置をすでに提案している(例えば、特許文献2を参照)。この装置は、既知の音響物性を有する基体、基材を介して超音波の送受信を行う超音波振動子、演算手段、画像構築手段等を含んで構成されている。この装置では、既知の音響物性を有する基体に、測定対象物及び既知の音響物性を有する参照物質を接して配置する。そしてこの状態で超音波を送信し、基体を介して測定対象物及び参照物質に超音波を入射させ、測定対象物及び参照物質からの超音波波形のインパルス応答を受信する。次いで、参照物質に入射した超音波波形のインパルス応答情報及び測定対象物に入射した超音波波形のインパルス応答情報から、規格化されたインパルス応答情報を得る。この規格化されたインパルス応答情報に基づき、奥行方向の音響物性分布(具体的には固有音響インピーダンス分布)を多重反射の影響を考慮して推定する演算を行う。そして、得られた奥行方向の音響物性分布に基づいて、音響物性像の画像データを構築し、所望とする超音波断層像を得るようにしている。 In view of such circumstances, the inventors of the present application have already proposed an improved ultrasonic image construction device (see, for example, Patent Document 2). This apparatus includes a substrate having known acoustic properties, an ultrasonic transducer for transmitting and receiving ultrasonic waves through the substrate, computing means, image constructing means, and the like. In this apparatus, an object to be measured and a reference substance having known acoustic properties are placed in contact with a substrate having known acoustic properties. Then, in this state, an ultrasonic wave is transmitted, the ultrasonic wave is made incident on the object to be measured and the reference substance through the substrate, and the impulse response of the ultrasonic waveform is received from the object to be measured and the reference substance. Next, normalized impulse response information is obtained from the impulse response information of the ultrasonic waveform incident on the reference substance and the impulse response information of the ultrasonic waveform incident on the measurement object. Based on this normalized impulse response information, calculations are performed to estimate the acoustic physical property distribution in the depth direction (specifically, the characteristic acoustic impedance distribution) in consideration of the influence of multiple reflections. Then, image data of an acoustic physical property image is constructed based on the obtained acoustic physical property distribution in the depth direction, and a desired ultrasonic tomographic image is obtained.

特開2006-271765号公報JP 2006-271765 A 特許第6361001号公報Japanese Patent No. 6361001

しかしながら、上記装置では、奥行方向に沿って深くなる方向、つまり測定対象物の前面側から背面側に向かう方向に固有音響インピーダンス分布を順次推定していく手法を採用している。このため、深度が大きくなるほど固有音響インピーダンスの推定値の誤差が蓄積、拡大するという欠点があり、より精度の高い画像を構築するうえでのマイナス要因となっていた。また、このような場合、超音波断層像における奥行方向に多くの筋が入ってしまうため、微細な内部の層構造が感覚的に理解しやすいとは言い難く、いまだ改良の余地があった。 However, the above apparatus adopts a method of sequentially estimating the characteristic acoustic impedance distribution in the direction of increasing depth along the depth direction, that is, in the direction from the front side to the back side of the object to be measured. For this reason, there is a drawback that the error in the estimated value of the characteristic acoustic impedance accumulates and expands as the depth increases, which is a negative factor in constructing a more accurate image. Moreover, in such a case, many streaks appear in the depth direction in the ultrasonic tomographic image, so it is difficult to say that the fine internal layer structure is intuitively understandable, and there is still room for improvement.

なお、培養細胞なども非常に薄い測定対象物の一種であって、その内部には核や細胞骨格などの微細な構造を有しているが、このようなものを測定対象物とした場合において、従来その微細構造を感覚的に理解しやすい態様にて超音波断層像化できる方法は提案されていなかった。 Cultured cells are also a type of very thin measurement object, and have fine structures such as nuclei and cytoskeleton inside them. However, there has been no proposal for a method capable of producing an ultrasonic tomographic image of the fine structure in a manner that is intuitively understandable.

本発明は上記の課題に鑑みてなされたものであり、その主たる目的は、微細な内部構造を有する非常に薄い測定対象物の超音波断層像を、感覚的にその微細な内部構造が理解しやすい態様にて比較的簡単にかつ高い精度で構築することができる超音波画像構築方法、超音波画像構築装置、超音波画像構築プログラムを提供することにある。 SUMMARY OF THE INVENTION The present invention has been made in view of the above problems, and its main purpose is to sensuously understand the fine internal structure of an ultrasonic tomographic image of a very thin measurement object having a fine internal structure. An object of the present invention is to provide an ultrasonic image construction method, an ultrasonic image construction apparatus, and an ultrasonic image construction program capable of constructing an image relatively easily and with high accuracy in an easy manner.

また、本発明のさらなる目的は、超音波を利用して得た情報に基づいて、培養細胞の内部構造や、皮膚等に代表される生体の軟組織の内部構造を簡便にかつ正確に観察、評価等できる画像を構築することにある。 A further object of the present invention is to simply and accurately observe and evaluate the internal structure of cultured cells and the internal structure of soft tissues of a living body such as skin, based on information obtained using ultrasound. It is to construct an image that can be equivalent.

上記の課題を解決するために、本願発明者らは測定対象物内の固有音響インピーダンスの推定値を適宜補正することを思い付き、鋭意検討を行った。その結果、1)測定対象物と同じく基体に接して存在する第1参照物質とは別の位置(具体的には、測定対象物において基体と接していない反対の側(奥行方向の奥側))の位置)に存在する第2参照物質の固有音響インピーダンス値を補正演算にて利用したり、2)そのような位置に第2参照物質が存在していなくても、所定の条件を満たす場合には測定対象物内の所定部位の固有音響インピーダンス値(あるいはその推定値の平均)を仮想参照部位の固有音響インピーダンス値と定義して補正演算にて利用したりすれば、深度が大きくなるほど蓄積、拡大しやすい固有音響インピーダンスの推定値の誤差を減じることができ、奥行方向の固有音響インピーダンス分布のより正確な推定が可能となることを新たに知見した。そして、本願発明者らはこれらの知見に基づいてさらに鋭意研究を進めることにより、下記に列挙する解決手段[1]~[15]を想到するに至ったのである。 In order to solve the above problems, the inventors of the present application came up with the idea of appropriately correcting the estimated value of the specific acoustic impedance in the object to be measured, and conducted earnest studies. As a result, 1) a position different from the first reference substance existing in contact with the substrate like the measurement object (specifically, the opposite side of the measurement object that is not in contact with the substrate (deep side in the depth direction) ) is used in the correction calculation, or 2) when a predetermined condition is satisfied even if the second reference material does not exist at such a position. For example, if the specific acoustic impedance value (or the average of the estimated values) of a predetermined part in the measurement object is defined as the specific acoustic impedance value of the virtual reference part and used in the correction calculation, the deeper the depth, the more accumulated , it is possible to reduce the error in the estimated value of the characteristic acoustic impedance, which tends to expand, and to make it possible to more accurately estimate the characteristic acoustic impedance distribution in the depth direction. Based on these findings, the inventors of the present application have further pursued research, and have arrived at the solutions [1] to [15] listed below.

[1]基体に測定対象物及び第1参照物質が接して存在し、かつ前記測定対象物における前記基体とは反対の側に第2参照物質が接して存在した状態で超音波を送信し、前記基体を介して超音波を入射させたときの前記測定対象物、前記第1参照物質及び前記第2参照物質からの超音波波形のインパルス応答を受信する送受信ステップと、前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内及び前記第2参照物質内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する推定ステップと、前記推定ステップにて得た前記第2参照物質の固有音響インピーダンスの推定値を前記第2参照物質の固有音響インピーダンスの実際値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する補正ステップと、前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築ステップとを有することを特徴とする超音波画像構築方法。 [1] transmitting ultrasonic waves in a state in which an object to be measured and a first reference substance are in contact with a substrate, and a second reference substance is in contact with the side of the object to be measured opposite to the substrate; a transmitting/receiving step of receiving an impulse response of an ultrasonic waveform from the measurement object, the first reference material, and the second reference material when ultrasonic waves are incident through the base; and an impulse of the ultrasonic waveform. an estimation step of sequentially estimating the intrinsic acoustic impedance in the measurement object and in the second reference material from the near side to the far side in the depth direction based on the normalized impulse response information obtained from the response information; The estimated value of the specific acoustic impedance of the object to be measured is obtained by replacing the estimated value of the specific acoustic impedance of the second reference material obtained in the estimation step with the actual value of the specific acoustic impedance of the second reference material. After correction, it has a correction step of estimating the characteristic acoustic impedance distribution in the depth direction, and an image construction step of constructing image data of an acoustic physical property image based on the characteristic acoustic impedance distribution in the depth direction. ultrasound image construction method.

[2]基体に測定対象物及び第1参照物質が接して存在した状態で超音波を送信し、前記基体を介して超音波を入射させたときの前記測定対象物及び前記第1参照物質からの超音波波形のインパルス応答を受信する送受信ステップと、前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する推定ステップと、前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値、または、前記測定対象物が生体の軟組織である場合の前記測定対象物内における特定深度の固有音響インピーダンスの推定値の平均を仮想参照部位の固有音響インピーダンス値と定義し、前記推定ステップにて得た前記測定対象物内における前記特定深度の固有音響インピーダンスの推定値を前記仮想参照部位の固有音響インピーダンス値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する補正ステップと、前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築ステップとを有することを特徴とする超音波画像構築方法。 [2] Ultrasonic waves are transmitted in a state in which the object to be measured and the first reference substance are in contact with a substrate, and when the ultrasonic waves are incident through the substrate, from the object to be measured and the first reference substance Based on the transmission/reception step of receiving the impulse response of the ultrasonic waveform and the normalized impulse response information obtained from the impulse response information of the ultrasonic waveform, the specific acoustic impedance in the measurement object is measured in the depth direction on the front side an estimating step of sequentially estimating toward the far side from the measuring object, and the specific acoustic impedance value when a substance having a known specific acoustic impedance value is uniformly present in the layer direction at a specific depth in the measurement object, Alternatively, when the measurement object is a soft tissue of a living body, an average of estimated values of specific acoustic impedance at a specific depth in the measurement object is defined as a specific acoustic impedance value of the virtual reference site, and obtained in the estimation step. After correcting the estimated value of the specific acoustic impedance in the measurement object through a calculation that replaces the estimated value of the specific acoustic impedance at the specific depth in the measurement object with the specific acoustic impedance value of the virtual reference part, An ultrasonic image construction characterized by comprising a correction step of estimating a characteristic acoustic impedance distribution in the depth direction and an image construction step of constructing image data of an acoustic physical property image based on the characteristic acoustic impedance distribution in the depth direction. Method.

[3]前記補正ステップでは、前記第2参照物質の固有音響インピーダンスの実際値または前記仮想参照部位の固有音響インピーダンス値に置き換えるときの値を最大とし、前記測定対象物の前面側に向かうほど小さくなる補正値を加算または減算することにより、前記測定対象物内の固有音響インピーダンスの推定値を補正することを特徴とする手段1または2に記載の超音波画像構築方法。 [3] In the correcting step, the actual value of the specific acoustic impedance of the second reference material or the value when replacing with the specific acoustic impedance value of the virtual reference part is maximized, and decreases toward the front side of the measurement object. 3. The ultrasonic image construction method according to means 1 or 2, wherein the estimated value of the intrinsic acoustic impedance in the object to be measured is corrected by adding or subtracting a correction value of .

[4]前記推定ステップでは、前記測定対象物内において、異なる固有音響インピーダンスを持つ無損失の微小伝送路が奥行方向に連なって伝送路の集合体をなしていると仮定して、手前側の前記微小伝送路の固有音響インピーダンスの推定結果に基づきその奥側に隣接する前記微小伝送路の固有音響インピーダンスを推定する演算を行うことを特徴とする手段1乃至3のいずれか1項に記載の超音波画像構築方法。 [4] In the estimation step, it is assumed that lossless minute transmission lines having different characteristic acoustic impedances are connected in the depth direction to form a collection of transmission lines in the object to be measured. 4. The method according to any one of means 1 to 3, characterized in that calculation for estimating the characteristic acoustic impedance of the minute transmission line adjacent to the rear side thereof is performed based on the estimation result of the characteristic acoustic impedance of the minute transmission line. Ultrasound image construction method.

[5]前記測定対象物は、培養細胞であり、前記第1参照物質及び前記第2参照物質は、前記培養細胞を培養するための培養液であり、前記基体は、前記培養液を収容する培養容器の一部であることを特徴とする請求項に記載の超音波画像構築方法。 [5] The measurement object is a cultured cell, the first reference substance and the second reference substance are a culture medium for culturing the cultured cell, and the substrate contains the culture medium. 2. The method of claim 1 , wherein the ultrasonic image construction method is part of a culture vessel.

[6]前記測定対象物は、組織表面に対して略平行な方向に延びる血管を含む生体の軟組織としての皮膚であり、前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値は、前記血管中に存在する血液の固有音響インピーダンス値であることを特徴とする請求項に記載の超音波画像構築方法。 [6] The measurement object is skin as a soft tissue of a living body containing blood vessels extending in a direction substantially parallel to the tissue surface, and has a known specific acoustic impedance value at a specific depth within the measurement object. 3. The ultrasonic image constructing method according to claim 2 , wherein the characteristic acoustic impedance value when the substance exists uniformly in the layer direction is the characteristic acoustic impedance value of blood existing in the blood vessel. .

[7]前記測定対象物が生体の軟組織である場合の前記測定対象物内における特定深度の固有音響インピーダンスの推定値の平均は、前記測定対象物が皮膚である場合の当該皮膚内における特定深度の固有音響インピーダンスの推定値の平均であり、前記特定深度は、前記音響物性像の画像データが構築される領域よりも深い位置に設定されることを特徴とする請求項に記載の超音波画像構築方法。 [7] When the measurement object is a soft tissue of a living body, the average of the estimated values of the specific acoustic impedance at a specific depth in the measurement object is the specific depth in the skin when the measurement object is the skin 3. The ultrasonic wave according to claim 2 , wherein the specific depth is set at a position deeper than the region in which the image data of the acoustic physical property image is constructed. Image construction method.

[8]基体と、前記基体に測定対象物及び第1参照物質が接して存在し、かつ前記測定対象物における前記基体とは反対の側に第2参照物質が接して存在した状態で超音波を送信し、前記基体を介して超音波を入射させたときの前記測定対象物、前記第1参照物質及び前記第2参照物質からの超音波波形のインパルス応答を受信する超音波振動子と、前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内及び前記第2参照物質内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する第1演算手段と、前記第1演算手段により得た前記第2参照物質の固有音響インピーダンスの推定値を前記第2参照物質の固有音響インピーダンスの実際値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する第2演算手段と、前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築手段とを備えたことを特徴とする超音波画像構築装置。 [8] ultrasonic waves in a state in which a substrate, an object to be measured and a first reference substance are in contact with the substrate, and a second reference substance is in contact with the side of the object to be measured opposite to the substrate; and receives impulse responses of ultrasonic waveforms from the measurement object, the first reference material, and the second reference material when ultrasonic waves are incident through the base; Based on the normalized impulse response information obtained from the impulse response information of the ultrasonic waveform, the intrinsic acoustic impedance in the measurement object and in the second reference material is sequentially measured from the near side to the far side in the depth direction. Through a first calculation means for estimating and a calculation for replacing the estimated value of the specific acoustic impedance of the second reference material obtained by the first calculation means with the actual value of the specific acoustic impedance of the second reference material, the measurement object A second calculating means for estimating the characteristic acoustic impedance distribution in the depth direction after correcting the estimated value of the characteristic acoustic impedance in the second calculation means, and constructing image data of the acoustic physical property image based on the characteristic acoustic impedance distribution in the depth direction. An ultrasound image construction device comprising: image construction means for:

[9]基体と、前記基体に測定対象物及び第1参照物質が接して存在した状態で超音波を送信し、前記基体を介して超音波を入射させたときの前記測定対象物及び前記第1参照物質からの超音波波形のインパルス応答を受信する超音波振動子と、前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する第1演算手段と、前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値、または、前記測定対象物が生体の軟組織である場合の前記測定対象物内における特定深度の固有音響インピーダンスの推定値の平均を仮想参照部位の固有音響インピーダンス値と定義し、前記第1演算手段により得た前記測定対象物内における前記特定深度の固有音響インピーダンスの推定値を前記仮想参照部位の固有音響インピーダンス値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する第2演算手段と、 前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築手段とを備えたことを特徴とする超音波画像構築装置。 [9] Ultrasonic waves are transmitted in a state in which a substrate, an object to be measured and a first reference substance are in contact with the substrate, and the object to be measured and the first reference substance when the ultrasonic waves are incident through the substrate. 1. An ultrasonic transducer that receives an impulse response of an ultrasonic waveform from a reference substance, and a specific acoustic impedance in the measurement object based on normalized impulse response information obtained from the impulse response information of the ultrasonic waveform and a first computing means for estimating sequentially from the front side to the back side in the depth direction, and a substance having a known specific acoustic impedance value at a specific depth in the measurement object is uniformly present in the layer direction defined as the intrinsic acoustic impedance value of the virtual reference site Then, through calculation for replacing the estimated value of the specific acoustic impedance at the specific depth in the measurement object obtained by the first calculation means with the specific acoustic impedance value of the virtual reference part, the specific acoustic impedance in the measurement object a second calculation means for estimating the characteristic acoustic impedance distribution in the depth direction after correcting the estimated value of; and an image constructing means for constructing image data of an acoustic physical property image based on the characteristic acoustic impedance distribution in the depth direction. An ultrasound image construction device comprising:

[10]前記第2演算手段は、前記第2参照物質の固有音響インピーダンスの実際値または前記仮想参照部位の固有音響インピーダンス値に置き換えるときの値を最大し、前記測定対象物の前面側に向かうほど小さくなる補正値を加算または減算することにより、前記測定対象物内の固有音響インピーダンスの推定値を補正することを特徴とする手段8または9に記載の超音波画像構築装置。 [10] The second calculating means maximizes the actual value of the specific acoustic impedance of the second reference material or the value when replacing it with the specific acoustic impedance value of the virtual reference part, and directs it toward the front side of the measurement object. 10. The ultrasonic image constructing apparatus according to means 8 or 9, wherein the estimated value of the intrinsic acoustic impedance in the object to be measured is corrected by adding or subtracting a correction value that becomes smaller as the correction value decreases.

[11]前記第1演算手段は、前記測定対象物内において、異なる固有音響インピーダンスを持つ無損失の微小伝送路が奥行方向に連なって伝送路の集合体をなしていると仮定して、手前側の前記微小伝送路の固有音響インピーダンスの推定結果に基づきその奥側に隣接する前記微小伝送路の固有音響インピーダンスを推定する演算を行うことを特徴とする手段8乃至10のいずれか1項に記載の超音波画像構築装置。 [11] The first computing means assumes that lossless minute transmission lines having different characteristic acoustic impedances are connected in the depth direction to form a collection of transmission lines in the object to be measured. 11. The method according to any one of means 8 to 10, wherein calculation is performed for estimating the characteristic acoustic impedance of the minute transmission line adjacent to the back side based on the estimation result of the characteristic acoustic impedance of the minute transmission line on the side. The ultrasound image construction device described.

[12]プロセッサに、基体に測定対象物及び第1参照物質が接して存在し、かつ前記測定対象物における前記基体とは反対の側に第2参照物質が接して存在した状態で超音波振動子に超音波を送信させ、前記基体を介して超音波を入射させたときの前記測定対象物、前記第1参照物質及び前記第2参照物質からの超音波波形のインパルス応答を前記超音波振動子に受信させる送受信ステップと、前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内及び前記第2参照物質内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する推定ステップと、前記推定ステップにて得た前記第2参照物質の固有音響インピーダンスの推定値を前記第2参照物質の固有音響インピーダンスの実際値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する補正ステップと、 前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築ステップとを実行させるための超音波画像構築プログラム。 [12] The processor performs ultrasonic vibration in a state in which an object to be measured and a first reference substance are in contact with a substrate, and a second reference substance is in contact with the side of the object to be measured opposite to the substrate. ultrasonic wave is transmitted from the element, and the impulse response of the ultrasonic waveform from the measurement object, the first reference substance and the second reference substance when the ultrasonic wave is incident through the base is the ultrasonic vibration Based on the transmission/reception step of causing the child to receive and the normalized impulse response information obtained from the impulse response information of the ultrasonic waveform, the characteristic acoustic impedance in the measurement object and in the second reference material is measured in the depth direction. Through an estimation step of sequentially estimating from the side to the back side, and an operation of replacing the estimated value of the specific acoustic impedance of the second reference material obtained in the estimation step with the actual value of the specific acoustic impedance of the second reference material a correction step of estimating a specific acoustic impedance distribution in the depth direction after correcting the estimated value of the specific acoustic impedance in the measurement object; and an image of an acoustic physical property image based on the specific acoustic impedance distribution in the depth direction. An ultrasound image construction program for performing an image construction step of constructing data.

[13]プロセッサに、基体に測定対象物及び第1参照物質が接して存在した状態で超音波振動子に超音波を送信させ、前記基体を介して超音波を入射させたときの前記測定対象物及び前記第1参照物質からの超音波波形のインパルス応答を前記超音波振動子に受信させる送受信ステップと、前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する推定ステップと、前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値、または、前記測定対象物が生体の軟組織である場合の前記測定対象物内における特定深度の固有音響インピーダンスの推定値の平均を仮想参照部位の固有音響インピーダンス値と定義し、前記推定ステップにて得た前記測定対象物内における前記特定深度の固有音響インピーダンスの推定値を前記仮想参照部位の固有音響インピーダンス値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する補正ステップと、 前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築ステップとを実行させるための超音波画像構築プログラム。 [13] When the processor causes the ultrasonic transducer to transmit ultrasonic waves in a state in which the object to be measured and the first reference substance are in contact with the substrate, and the ultrasonic waves are incident through the substrate, the object to be measured a transmitting/receiving step of causing the ultrasonic transducer to receive impulse responses of ultrasonic waveforms from the object and the first reference material; and based on normalized impulse response information obtained from the impulse response information of the ultrasonic waveforms, an estimating step of sequentially estimating the intrinsic acoustic impedance within the object to be measured from the near side to the far side in the depth direction; The average of the specific acoustic impedance value when the measurement object exists uniformly or the estimated value of the specific acoustic impedance at a specific depth in the measurement object when the measurement object is a soft tissue of a living body is the virtual reference part defined as a specific acoustic impedance value, and through a calculation that replaces the estimated value of the specific acoustic impedance at the specific depth in the measurement object obtained in the estimation step with the specific acoustic impedance value of the virtual reference part, the measurement object A correcting step of estimating the characteristic acoustic impedance distribution in the depth direction after correcting the estimated value of the characteristic acoustic impedance in the image, and constructing image data of the acoustic physical property image based on the characteristic acoustic impedance distribution in the depth direction. An ultrasound image construction program for performing the construction steps.

[14]前記補正ステップでは、前記第2参照物質の固有音響インピーダンスの実際値または前記仮想参照部位の固有音響インピーダンス値に置き換えるときの値を最大とし、前記測定対象物の前面側に向かうほど小さくなる補正値を加算または減算することにより、前記測定対象物内の固有音響インピーダンスの推定値を補正することを特徴とする請求項12または13に記載の超音波画像構築プログラム。 [14] In the correcting step, the actual value of the specific acoustic impedance of the second reference material or the value when replacing with the specific acoustic impedance value of the virtual reference part is maximized, and decreases toward the front side of the measurement object. 14. The program for constructing an ultrasonic image according to claim 12 or 13, wherein the estimated value of the characteristic acoustic impedance in the object to be measured is corrected by adding or subtracting a correction value of .

[15]前記推定ステップでは、前記測定対象物内において、異なる固有音響インピーダンスを持つ無損失の微小伝送路が奥行方向に連なって伝送路の集合体をなしていると仮定して、手前側の前記微小伝送路の固有音響インピーダンスの推定結果に基づきその奥側に隣接する前記微小伝送路の固有音響インピーダンスを推定する演算を行うことを特徴とする請求項12乃至14のいずれか1項に記載の超音波画像構築プログラム。 [15] In the estimating step, it is assumed that lossless minute transmission lines having different characteristic acoustic impedances are connected in the depth direction to form a collection of transmission lines in the object to be measured. 15. The method according to any one of claims 12 to 14, wherein calculation is performed for estimating the specific acoustic impedance of said micro-transmission channel adjacent on the far side based on the result of estimating the specific acoustic impedance of said micro-transmission channel. ultrasound image construction program.

上記手段1~15に記載の発明によると、微細な内部構造を有する非常に薄い測定対象物の超音波断層像を、感覚的にその微細な内部構造が理解しやすい態様にて比較的簡単にかつ高い精度で構築することができる。また、超音波を利用して得た情報に基づいて、培養細胞の内部構造や、皮膚等に代表される生体の軟組織の内部構造を簡便にかつ正確に観察、評価等できる画像を構築することができる。 According to the inventions described in the above means 1 to 15, an ultrasonic tomographic image of a very thin measurement object having a fine internal structure can be relatively easily obtained in a manner in which the fine internal structure can be intuitively understood. And it can be constructed with high accuracy. In addition, based on the information obtained using ultrasound, constructing an image that enables simple and accurate observation and evaluation of the internal structure of cultured cells and the internal structure of soft tissues of the living body such as the skin. can be done.

本発明を具体化した第1の実施形態の超音波画像構築装置を示す概略構成図。1 is a schematic configuration diagram showing an ultrasonic image constructing apparatus of a first embodiment embodying the present invention; FIG. 第1の実施形態の超音波画像構築装置の電気的構成を示すブロック図。FIG. 2 is a block diagram showing the electrical configuration of the ultrasonic image construction apparatus of the first embodiment; FIG. X-Yステージの移動に伴う超音波の走査範囲の一例を示す概略図。FIG. 4 is a schematic diagram showing an example of an ultrasonic scanning range accompanying movement of the XY stage; (a)は実際に測定を行ったときにおける測定対象物からの反射波形の取得についての説明図、(b)は測定対象物を微小伝送路に見立てたときにおける反射波形の取得についての説明図。(a) is an explanatory diagram of acquisition of the reflected waveform from the measurement object when the measurement is actually performed, and (b) is an explanatory diagram of acquisition of the reflected waveform when the measurement object is regarded as a minute transmission path. . (a)は実際に測定を行ったときにおける参照物質からの反射波形の取得についての説明図、(b)は参照物質を微小伝送路に見立てたときにおける反射波形の取得についての説明図。(a) is an explanatory diagram of acquisition of a reflected waveform from a reference substance when actually performing measurement, and (b) is an explanatory diagram of acquisition of a reflected waveform when the reference substance is regarded as a minute transmission path. 各微小伝送路の特性インピーダンスを推定していく様子を概念的に示した図。FIG. 4 is a diagram conceptually showing how the characteristic impedance of each minute transmission path is estimated. 多重反射の影響について概念的に示した図。FIG. 2 is a diagram conceptually showing the influence of multiple reflection; ヒトの首の皮膚の層構造を示す概略断面図。Schematic cross-sectional view showing the layered structure of human neck skin. 第1実施形態における補正ステップを説明するためのグラフ。5 is a graph for explaining correction steps in the first embodiment; 第1実施形態において固有音響インピーダンス像の構築についての演算処理を説明するためのフローチャート。4 is a flowchart for explaining arithmetic processing for constructing a specific acoustic impedance image in the first embodiment; 本発明を具体化した第2の実施形態の超音波画像構築装置における補正ステップを説明するためのものであって、(a)は補正ステップを実施していないときの固有音響インピーダンス像、(b)は補正ステップを実施したときの固有音響インピーダンス像。It is for explaining the correction step in the ultrasonic image construction apparatus of the second embodiment embodying the present invention, (a) is a specific acoustic impedance image when the correction step is not performed, (b ) is the characteristic acoustic impedance image when the correction step is performed. 本発明を具体化した第3の実施形態の超音波画像構築装置を示す概略構成図。FIG. 3 is a schematic block diagram showing an ultrasonic image constructing apparatus of a third embodiment embodying the present invention; 測定対象物である培養細胞(グリア細胞)を概略的に示した図。A diagram schematically showing a cultured cell (glial cell), which is an object to be measured. 第3の実施形態の超音波画像構築装置を用いて観察された培養細胞(グリア細胞)の固有音響インピーダンス像。FIG. 10 is a specific acoustic impedance image of cultured cells (glial cells) observed using the ultrasonic image construction apparatus of the third embodiment; FIG.

[第1の実施形態]
以下、本発明の超音波画像構築方法及び装置を具体化した第1の実施形態を図1~図10に基づき詳細に説明する。
[First embodiment]
A first embodiment embodying an ultrasonic image constructing method and apparatus of the present invention will be described in detail below with reference to FIGS. 1 to 10. FIG.

図1は、本実施形態の超音波画像構築装置1を示す概略構成図である。図1に示されるように、本実施形態の超音波画像構築装置1は、超音波を用いて皮膚8を観察、診断等するための装置であって、パルス励起型超音波顕微鏡2と、パーソナルコンピュータ(PC)3とを備えている。 FIG. 1 is a schematic configuration diagram showing an ultrasonic image constructing apparatus 1 of this embodiment. As shown in FIG. 1, an ultrasonic image constructing apparatus 1 of the present embodiment is an apparatus for observing, diagnosing, etc. skin 8 using ultrasonic waves, and includes a pulse excitation ultrasonic microscope 2 and a personal A computer (PC) 3 is provided.

パルス励起型超音波顕微鏡2は、ステージ4を有する顕微鏡本体5と、ステージ4の下方に設置された超音波プローブ6とを備える。パルス励起型超音波顕微鏡2の超音波プローブ6は、PC3と電気的に接続されている。 A pulse excitation ultrasonic microscope 2 includes a microscope main body 5 having a stage 4 and an ultrasonic probe 6 installed below the stage 4 . The ultrasonic probe 6 of the pulse excitation ultrasonic microscope 2 is electrically connected to the PC3.

本実施形態のステージ4は、ユーザの手動操作により、水平方向(即ちX方向及びY方向)に移動できるように構成されている。このステージ4には、測定対象物を接触させて配置するための樹脂プレート9が固定されている。ここでの測定対象物は、組織表面に対して略平行な方向に延びる血管を含む生体の軟組織(具体的には、肌組織;皮膚8)である。本実施形態では、ヒトの皮膚8を樹脂プレート9に直接押し付けることにより測定等を行っている。また、既知の音響物性を有する基体としての樹脂プレート9は、超音波を透過させることができる平板状部材であって、測定対象物である皮膚8よりも硬い材料からなる。このような形状及び硬さの部材を基体として用いた場合、測定対象物である皮膚8を確実に密着配置することが可能となり、奥行方向の固有音響インピーダンス分布を正確に推定可能となる結果、ひいては画像構築の精度が向上する。なお、本実施形態では厚さ1.4mmのポリスチレン板が用いられている。勿論、ポリスチレン以外の樹脂からなる板材などを用いることも許容される。 The stage 4 of this embodiment is configured to be movable in the horizontal direction (that is, the X direction and the Y direction) by a user's manual operation. A resin plate 9 is fixed to the stage 4 for placing an object to be measured in contact therewith. The object to be measured here is a living body's soft tissue (specifically, skin tissue; skin 8) containing blood vessels extending in a direction substantially parallel to the tissue surface. In this embodiment, the measurement or the like is performed by directly pressing the human skin 8 against the resin plate 9 . The resin plate 9 as a substrate having known acoustic properties is a plate-like member that can transmit ultrasonic waves, and is made of a material harder than the skin 8 that is the object to be measured. When a member having such a shape and hardness is used as a substrate, it becomes possible to securely place the skin 8, which is the object to be measured, in close contact with the body, and as a result, it becomes possible to accurately estimate the characteristic acoustic impedance distribution in the depth direction. As a result, the accuracy of image construction is improved. A polystyrene plate having a thickness of 1.4 mm is used in this embodiment. Of course, it is permissible to use a plate member made of a resin other than polystyrene.

この樹脂プレート9において皮膚8が接触配置される側である上面には、第1参照物質としてのリファレンス部材10があらかじめ設置されている。リファレンス部材10は、樹脂プレート9とは異なる既知の音響物性を有している。本実施形態では例えばアクリル樹脂(アクリル接着剤)を付着させることによりリファレンス部材10としているが、勿論これに限定されるわけではない。リファレンス部材10に対して密着させることが可能なものであれば、樹脂材以外のもの(例えばガラス材、金属材、セラミック材など)をリファレンス部材10としてもよい。あるいは、このようなリファレンス部材10を設置する代わりに、例えば樹脂プレート9の上面に接するように水等を存在させておき、これを第1参照物質として用いてもよい。なお、樹脂プレート9にリファレンス部材10をあらかじめ設置しておくことにより、装置が置かれる環境の変化等に依存せず、リファレンス部材10に入射した超音波波形のインパルス応答情報を正確にかつ安定的に取得することができる。 A reference member 10 as a first reference substance is previously installed on the upper surface of the resin plate 9 on which the skin 8 is placed in contact. The reference member 10 has known acoustic properties different from those of the resin plate 9 . In this embodiment, for example, the reference member 10 is formed by adhering an acrylic resin (acrylic adhesive), but the present invention is, of course, not limited to this. The reference member 10 may be made of a material other than a resin material (for example, a glass material, a metal material, a ceramic material, etc.) as long as it can be brought into close contact with the reference member 10 . Alternatively, instead of installing such a reference member 10, for example, water or the like may be present so as to be in contact with the upper surface of the resin plate 9, and this may be used as the first reference substance. By previously setting the reference member 10 on the resin plate 9, the impulse response information of the ultrasonic waveform incident on the reference member 10 can be obtained accurately and stably without depending on changes in the environment in which the apparatus is placed. can be obtained.

超音波プローブ6は、水などの超音波伝達媒体Wを貯留可能な貯留部11をその先端部に有するプローブ本体12と、プローブ本体12の略中心部に配置される超音波トランスデューサ13(超音波振動子)と、プローブ本体12を前記ステージ4の面方向に沿って二次元的に走査するためのX-Yステージ14とを備える。プローブ本体12の貯留部11は上部が開口しており、その貯留部11の開口側を上向きにした状態で超音波プローブ6がステージ4の下方に設置されている。 The ultrasonic probe 6 includes a probe body 12 having a storage portion 11 capable of storing an ultrasonic transmission medium W such as water at its tip, and an ultrasonic transducer 13 (ultrasonic wave and an XY stage 14 for two-dimensionally scanning the probe body 12 along the surface direction of the stage 4 . The storage part 11 of the probe main body 12 is open at the top, and the ultrasonic probe 6 is installed below the stage 4 with the opening side of the storage part 11 facing upward.

超音波トランスデューサ13は、例えば酸化亜鉛の薄膜圧電素子16とサファイアロッドの音響レンズ17とによって構成される。この超音波トランスデューサ13は、パルス励起されることで樹脂プレート9の下面側から皮膚8及びリファレンス部材10に対して超音波を照射する。超音波トランスデューサ13が照射する超音波は、貯留部11の超音波伝達媒体Wを介して円錐状に収束されて樹脂プレート9の上面(皮膚8の表面付近)で焦点を結ぶようになっている。なお本実施形態では、超音波トランスデューサ13として、口径1.2mm、焦点距離1.5mm、中心周波数80MHz、帯域幅50~105MHz(-6dB)の仕様のものを用いている。 The ultrasonic transducer 13 is composed of, for example, a zinc oxide thin film piezoelectric element 16 and a sapphire rod acoustic lens 17 . The ultrasonic transducer 13 is pulse-excited to irradiate the skin 8 and the reference member 10 with ultrasonic waves from the lower surface side of the resin plate 9 . The ultrasonic waves emitted by the ultrasonic transducer 13 are converged into a conical shape via the ultrasonic transmission medium W of the reservoir 11 and focused on the upper surface of the resin plate 9 (near the surface of the skin 8). . In this embodiment, as the ultrasonic transducer 13, one with a diameter of 1.2 mm, a focal length of 1.5 mm, a center frequency of 80 MHz, and a bandwidth of 50 to 105 MHz (-6 dB) is used.

図2は、本実施形態の超音波画像構築装置1の電気的な構成を示すブロック図である。 FIG. 2 is a block diagram showing the electrical configuration of the ultrasonic image constructing apparatus 1 of this embodiment.

図2に示されるように、超音波プローブ6は、超音波トランスデューサ13、X-Yステージ14、パルス発生回路21、受信回路22、送受波分離回路23、検波回路24、A/D変換回路25、エンコーダ26、コントローラ27を備える。 As shown in FIG. 2, the ultrasonic probe 6 includes an ultrasonic transducer 13, an XY stage 14, a pulse generation circuit 21, a reception circuit 22, a transmission/reception separation circuit 23, a detection circuit 24, and an A/D conversion circuit 25. , an encoder 26 and a controller 27 .

走査手段としてのX-Yステージ14は、超音波の照射点を二次元的に走査させるためのXステージ14X及びYステージ14Yを備えるとともに、それぞれのステージ14X,14Yを駆動するモータ28X,28Yを備えている。これらのモータ28X,28Yとしては、ステッピングモータやリニアモータが使用される。 The XY stage 14 as scanning means includes an X stage 14X and a Y stage 14Y for two-dimensionally scanning the irradiation point of the ultrasonic waves, and motors 28X and 28Y for driving the respective stages 14X and 14Y. I have. A stepping motor or a linear motor is used as these motors 28X and 28Y.

各モータ28X,28Yにはコントローラ27が接続されており、該コントローラ27の駆動信号に応答してモータ28X,28Yが駆動される。これらモータ28X,28Yの駆動により、Xステージ14Xを連続走査(連続送り)するとともに、Yステージ14Yを間欠送りとなるよう制御することで、X-Yステージ14の高速走査が可能となっている。 A controller 27 is connected to each of the motors 28X and 28Y, and the motors 28X and 28Y are driven in response to drive signals from the controller 27. FIG. By driving these motors 28X and 28Y, the X stage 14X is continuously scanned (continuously fed), and the Y stage 14Y is controlled to be intermittently fed, thereby enabling high-speed scanning of the XY stage 14. .

また、本実施形態においては、Xステージ14Xに対応してエンコーダ26が設けられ、エンコーダ26によりXステージ14Xの走査位置が検出される。具体的には、走査範囲を300×300個の測定点(ピクセル)に分割した場合、1回のX方向(水平方向)の走査が300分割される。そして、各測定点の位置がエンコーダ26によって検出されPC3に取り込まれる。PC3はそのエンコーダ26の出力に同期して駆動制御信号を生成し、その駆動制御信号をコントローラ27に供給する。コントローラ27は、この駆動制御信号に基づいてモータ28Xを駆動する。また、コントローラ27は、エンコーダ26の出力信号に基づきX方向の1ラインの走査が終了した時点でモータ28Yを駆動して、Yステージ14YをY方向に1ピクセル分移動させる。 Further, in this embodiment, an encoder 26 is provided corresponding to the X stage 14X, and the scanning position of the X stage 14X is detected by the encoder 26. FIG. Specifically, when the scanning range is divided into 300×300 measurement points (pixels), one scanning in the X direction (horizontal direction) is divided into 300. Then, the position of each measuring point is detected by the encoder 26 and taken into the PC3. The PC 3 generates a drive control signal in synchronization with the output of the encoder 26 and supplies the drive control signal to the controller 27 . The controller 27 drives the motor 28X based on this drive control signal. Further, the controller 27 drives the motor 28Y based on the output signal of the encoder 26 when scanning of one line in the X direction is completed, and moves the Y stage 14Y by one pixel in the Y direction.

さらに、コントローラ27は、駆動制御信号に同期してトリガ信号を生成してパルス発生回路21に供給する。これにより、パルス発生回路21において、そのトリガ信号に同期したタイミングで励起パルスが生成される。その励起パルスが送受波分離回路23を介して超音波トランスデューサ13に供給される結果、超音波トランスデューサ13から超音波が照射される。 Furthermore, the controller 27 generates a trigger signal in synchronization with the drive control signal and supplies it to the pulse generation circuit 21 . As a result, the pulse generation circuit 21 generates an excitation pulse at a timing synchronized with the trigger signal. As a result of the excitation pulse being supplied to the ultrasonic transducer 13 via the transmission/reception wave separation circuit 23, the ultrasonic transducer 13 emits ultrasonic waves.

図3は、X-Yステージ14の移動に伴う超音波の走査範囲R1の一例を示している。この例では、皮膚8を接触配置させる領域を包囲するようにリファレンス部材10が設けられている。そして、ヒトの皮膚8を当該領域に押し付けた状態で、リファレンス部材10がある位置から走査が開始される。そして、矢印で示すように、皮膚8の表面に沿ってX方向及びY方向に二次元的に走査が順次行われる。 FIG. 3 shows an example of an ultrasonic scanning range R1 accompanying movement of the XY stage 14. As shown in FIG. In this example, a reference member 10 is provided so as to surround the area where the skin 8 is to be placed in contact. Scanning is then started from a position where the reference member 10 is located while the human skin 8 is pressed against the area. Then, as indicated by arrows, scanning is sequentially performed two-dimensionally along the surface of the skin 8 in the X and Y directions.

超音波トランスデューサ13の薄膜圧電素子16は、送受波兼用の超音波振動子であり、皮膚8で反射した超音波(反射波)を電気信号に変換する。そして、その反射波の信号は、送受波分離回路23を介して受信回路22に供給される。受信回路22は、信号増幅回路を含んで構成されていて、反射波の信号を増幅して検波回路24に出力する。 The thin-film piezoelectric element 16 of the ultrasonic transducer 13 is an ultrasonic transducer for transmitting and receiving waves, and converts ultrasonic waves (reflected waves) reflected by the skin 8 into electrical signals. A signal of the reflected wave is supplied to the receiving circuit 22 via the transmitting/receiving wave separating circuit 23 . The receiving circuit 22 includes a signal amplifying circuit, amplifies the signal of the reflected wave, and outputs the amplified signal to the detecting circuit 24 .

検波回路24は、皮膚8からの反射波信号を検出するための回路であり、図示しないゲート回路を含む。本実施形態の検波回路24は、超音波トランスデューサ13で受信した反射波信号のなかから、皮膚8からの反射波信号やリファレンス部材10からの反射波信12月号を抽出する。そして、検波回路24で抽出された反射波信号は、A/D変換回路25に供給されてA/D変換された後、PC3に転送される。 The detection circuit 24 is a circuit for detecting reflected wave signals from the skin 8 and includes a gate circuit (not shown). The detection circuit 24 of this embodiment extracts the reflected wave signal from the skin 8 and the reflected wave signal from the reference member 10 from the reflected wave signals received by the ultrasonic transducer 13 . Then, the reflected wave signal extracted by the detection circuit 24 is supplied to the A/D conversion circuit 25 to be A/D converted, and then transferred to the PC 3 .

PC3は、CPU31(中央処理装置)、I/F回路32、メモリ33、記憶装置34、入力装置35、及び表示装置36を備え、それらはバス37を介して相互に接続されている。 The PC 3 includes a CPU 31 (central processing unit), an I/F circuit 32, a memory 33, a storage device 34, an input device 35, and a display device 36, which are interconnected via a bus 37.

CPU31は、メモリ33を利用して制御プログラムを実行し、システム全体を統括的に制御する。制御プログラムとしては、X-Yステージ14による二次元走査を制御するためのプログラム、超音波Bモードエコー像の元となる反射信号列のデータを固有音響インピーダンス像へ変換するためのプログラム、固有音響インピーダンス像を表示するためのプログラムなどを含む。なお、CPU31とは別に例えばDSP(Digital Signal Processor:デジタル信号プロセッサ)を設けて、そこでCPU31が行っている信号処理の一部を行わせてもよい。 The CPU 31 executes a control program using the memory 33 to centrally control the entire system. The control programs include a program for controlling two-dimensional scanning by the XY stage 14, a program for converting the data of the reflected signal train that is the basis of the ultrasonic B-mode echo image into a specific acoustic impedance image, a specific acoustic Includes programs for displaying impedance images. In addition, for example, a DSP (Digital Signal Processor) may be provided separately from the CPU 31 to perform part of the signal processing performed by the CPU 31 there.

I/F回路32は、超音波プローブ6との間で信号の授受を行うためのインターフェース(具体的には、USBインターフェース)である。I/F回路32は、超音波プローブ6に制御信号(コントローラ27への駆動制御信号)を出力したり、超音波プローブ6からの転送データ(A/D変換回路25から転送されるデータなど)を入力したりする役割を果たすものである。なお、超音波プローブ6との間で信号の授受を行う場合には、上記のような物理的なインターフェースに限定されることはなく、無線インターフェースを用いてもよい。 The I/F circuit 32 is an interface (specifically, a USB interface) for exchanging signals with the ultrasonic probe 6 . The I/F circuit 32 outputs control signals (driving control signals to the controller 27) to the ultrasonic probe 6, and transfers data from the ultrasonic probe 6 (data transferred from the A/D conversion circuit 25, etc.). It plays the role of inputting When transmitting and receiving signals to and from the ultrasonic probe 6, a wireless interface may be used without being limited to the physical interface as described above.

表示装置36は、例えば、液晶、プラズマ、有機EL(electroluminescence)等のモニタディスプレイである。表示装置36は、カラー表示、モノクロ表示を問わずに使用できるが、カラー表示であることが望ましい。この表示装置36は、皮膚8の表層の固有音響インピーダンス像を表示したり、各種設定の入力画面を表示したりするために用いられる。 The display device 36 is, for example, a monitor display such as liquid crystal, plasma, or organic EL (electroluminescence). The display device 36 can be used regardless of color display or monochrome display, but color display is desirable. This display device 36 is used to display the characteristic acoustic impedance image of the surface layer of the skin 8 and to display input screens for various settings.

入力装置35は、タッチパネル、マウス、キーボード、ポインティングデバイス等の入力ユーザインタフェースであって、ユーザからの要求や指示、パラメータの入力に用いられる。 The input device 35 is an input user interface such as a touch panel, mouse, keyboard, pointing device, etc., and is used to input requests, instructions, and parameters from the user.

記憶装置34は、磁気ディスク装置や光ディスク装置などのハードディスクドライブであり、各種の制御プログラム及び各種のデータを記憶している。メモリ33は、RAM(ランダムアクセスメモリ)やROM(リードオンリーメモリ)を含み、超音波測定のためにあらかじめ取得されたリファレンス部材10の反射波形とその固有音響インピーダンスとを保存する。CPU31は、入力装置35による指示に従い、プログラムやデータを記憶装置34からメモリ33へ転送し、それを逐次実行する。なお、CPU31が実行するプログラムとしては、メモリカード、フレキシブルディスク、光ディスクなどの記憶媒体に記憶されたプログラムや、通信媒体を介してダウンロードしたプログラムでもよく、その実行時には記憶装置34にインストールして利用する。 The storage device 34 is a hard disk drive such as a magnetic disk device or an optical disk device, and stores various control programs and various data. The memory 33 includes a RAM (random access memory) and a ROM (read only memory), and stores the reflected waveform of the reference member 10 and its characteristic acoustic impedance acquired in advance for ultrasonic measurement. The CPU 31 transfers programs and data from the storage device 34 to the memory 33 according to instructions from the input device 35 and executes them sequentially. The program executed by the CPU 31 may be a program stored in a storage medium such as a memory card, flexible disk, or optical disk, or a program downloaded via a communication medium. do.

次に、本実施形態の超音波画像構築装置1において、超音波Bモードエコー像の元となる反射信号列から固有音響インピーダンス像を構築する手法について説明する。 Next, a method for constructing a specific acoustic impedance image from a reflected signal train that is the basis of an ultrasonic B-mode echo image in the ultrasonic image constructing apparatus 1 of this embodiment will be described.

この超音波画像構築装置1では、第1参照物質であるリファレンス部材10に入射した超音波波形のインパルス応答情報、及び測定対象物である皮膚8に入射した超音波波形のインパルス応答情報から、規格化されたインパルス応答情報を得て、その規格化されたインパルス応答情報に基づいて奥行方向の音響物性分布を多重反射の影響を考慮して推定するようになっている。また、このような推定を行うために、本実施形態では、測定対象物内において、異なる固有音響インピーダンスを持つ無損失の微小伝送路51が奥行方向に連なって伝送路の集合体をなしていると仮定して、手前側の微小伝送路51の固有音響インピーダンスの推定結果に基づきその奥側に隣接する微小伝送路51の固有音響インピーダンスを推定する処理を順次繰り返すことにより、伝送路の奥行方向の音響物性分布(ここでは固有音響インピーダンス分布)を推定する演算を行うようになっている。このような演算は、メモリ33内に格納された所定のアルゴリズムに基づいて実行される。 In this ultrasonic image constructing apparatus 1, from the impulse response information of the ultrasonic waveform incident on the reference member 10, which is the first reference material, and the impulse response information of the ultrasonic waveform incident on the skin 8, which is the measurement object, the standard After obtaining normalized impulse response information, the acoustic physical property distribution in the depth direction is estimated based on the normalized impulse response information, taking into consideration the influence of multiple reflections. In order to perform such estimation, in the present embodiment, lossless minute transmission lines 51 having different characteristic acoustic impedances are connected in the depth direction to form a collection of transmission lines in the object to be measured. Assuming that, by sequentially repeating the process of estimating the characteristic acoustic impedance of the minute transmission line 51 adjacent to the back side based on the result of estimating the characteristic acoustic impedance of the minute transmission line 51 on the front side, the depth direction of the transmission line A calculation for estimating the acoustic physical property distribution (specific acoustic impedance distribution in this case) is performed. Such operations are performed based on predetermined algorithms stored in memory 33 .

このアルゴリズムは、超音波Bモードエコー像の元となる反射信号列を用いることによって、奥行方向の固有音響インピーダンスの分布の推定を行うアルゴリズムである。このアルゴリズムは、時間領域反射測定法(TDR法:Time Domain Reflectometry法)の原理を参考とするものであって、皮膚組織内部での多重反射を考慮した時間-周波数領域における解析を通じて、超音波Bモードエコー像の元となる反射信号列を奥行方向の固有音響インピーダンス像に変換するアルゴリズムである。以下、これについて具体的に説明する。 This algorithm is an algorithm for estimating the distribution of the characteristic acoustic impedance in the depth direction by using the reflected signal train that is the basis of the ultrasonic B-mode echo image. This algorithm is based on the principle of the time domain reflectometry method (TDR method: Time Domain Reflectometry method), and considers multiple reflections inside the skin tissue. Through analysis in the time-frequency domain, ultrasound B This is an algorithm that transforms the reflected signal train, which is the basis of the mode echo image, into the characteristic acoustic impedance image in the depth direction. This will be described in detail below.

図4(a)は、実際に測定を行ったときにおける測定対象物からの反射波形の取得について説明する図であり、図4(b)は、測定対象物を微小伝送路51に見立てたときにおける反射波形の取得について説明する図である。図5(a)は、実際に測定を行ったときにおける第1参照物質からの反射波形の取得について説明する図であり、図5(b)は、第1参照物質を微小伝送路51に見立てたときにおける反射波形の取得について説明する図である。 FIG. 4(a) is a diagram for explaining acquisition of the reflected waveform from the object to be measured when the measurement is actually performed, and FIG. FIG. 10 is a diagram for explaining acquisition of a reflected waveform in . FIG. 5(a) is a diagram for explaining acquisition of a reflected waveform from the first reference substance when actually performing measurement, and FIG. FIG. 10 is a diagram for explaining acquisition of a reflected waveform when

まず、図4(a)に示すように、超音波トランスデューサ13を作動させ、基体としての樹脂プレート9を介して、測定対象物にとって十分な焦点深度を持った超音波の収束ビームを送信する。そして、測定対象物である皮膚8に超音波の収束ビームを入射させ、そこからの反射波形を取得する。その時のインパルス応答Γ(ω)は、フーリエ変換を用いて、入射波S0と皮膚8からの反射波Stgt(ω)とから次式1のように表される。

Figure 0007213491000001
First, as shown in FIG. 4A, the ultrasonic transducer 13 is operated to transmit a focused ultrasonic beam having a sufficient depth of focus for the object to be measured through the resin plate 9 as a base. Then, a convergent beam of ultrasonic waves is made incident on the skin 8, which is an object to be measured, and a reflected waveform therefrom is acquired. The impulse response Γ 0 (ω) at that time is represented by the following equation 1 from the incident wave S0 and the reflected wave S tgt (ω) from the skin 8 using Fourier transform.
Figure 0007213491000001

この場合、固有音響インピーダンスが既知かつ均一であって、測定対象物に比べて十分な厚さを有するリファレンス部材10からの反射波形も取得する必要がある。リファレンス部材10からの反射波Sref(ω)は、リファレンス部材10の固有音響インピーダンスZrefと樹脂プレート9の固有音響インピーダンスZとを用いて、次式2のように表される。

Figure 0007213491000002
In this case, it is also necessary to acquire a reflected waveform from the reference member 10 having a known and uniform specific acoustic impedance and having a sufficient thickness compared to the object to be measured. A reflected wave S ref (ω) from the reference member 10 is expressed by the following equation 2 using the specific acoustic impedance Z ref of the reference member 10 and the specific acoustic impedance Z 0 of the resin plate 9 .
Figure 0007213491000002

また、測定対象物である皮膚8からのインパルス応答Γ(ω)は次式3で表される。ただし、このインパルス応答Γ(ω)には皮膚8の組織の奥の複数の界面から発生する反射が含まれているため、インパルス応答Γ(ω)は周波数特性を持ったものとなる。なお、ここまでの式により、第1参照物質に入射した超音波波形のインパルス応答情報及び測定対象物に入射した超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報が求められる。

Figure 0007213491000003
Also, the impulse response Γ 0 (ω) from the skin 8 which is the object to be measured is expressed by the following equation (3). However, since this impulse response Γ 0 (ω) includes reflections generated from a plurality of interfaces deep within the tissue of the skin 8, the impulse response Γ 0 (ω) has frequency characteristics. By the above formula, the normalized impulse response information obtained from the impulse response information of the ultrasonic waveform incident on the first reference substance and the impulse response information of the ultrasonic waveform incident on the measurement object is obtained.
Figure 0007213491000003

ここで、図6は、樹脂プレート9に接する微小伝送路51から順に、各微小伝送路51の特性インピーダンスZ、Z…Zを推定していく様子を概念的に示した図である。この図に示されるように、各微小伝送路51の特性インピーダンスZ、Z…Zを、樹脂プレート9に接している微小伝送路51から奥行方向に向かって順に推定していく。 Here, FIG. 6 is a diagram conceptually showing how the characteristic impedances Z 1 , Z 2 . . As shown in this figure, the characteristic impedances Z 1 , Z 2 .

図7は、多重反射の影響について概念的に示した図である。次式4は、インパルス応答Γ(ω)を逆フーリエ変換したg(t)を表したものであるが、その第1項は多重反射の影響を受けない(図7参照)。従って、この第1項の値から樹脂プレート9に接する微小伝送路51の特性インピーダンスZを、次式5のように推定することができる。

Figure 0007213491000004

Figure 0007213491000005
FIG. 7 is a diagram conceptually showing the influence of multiple reflections. The following equation 4 expresses g 0 (t) obtained by inverse Fourier transforming the impulse response Γ 0 (ω), the first term of which is not affected by multiple reflections (see FIG. 7). Therefore, the characteristic impedance Z1 of the minute transmission line 51 in contact with the resin plate 9 can be estimated from the value of the first term as shown in the following equation (5).
Figure 0007213491000004

Figure 0007213491000005

周波数領域の皮膚8の固有音響インピーダンスZxは、Γを用いると次式6のように表される。

Figure 0007213491000006
The characteristic acoustic impedance Zx 0 of the skin 8 in the frequency domain is expressed by the following Equation 6 using Γ 0 .
Figure 0007213491000006

また、Zxはさらに奥からのインパルス応答Γを用いて次式7のようにも表される。

Figure 0007213491000007
Zx 0 is also expressed by the following equation 7 using the impulse response Γ 1 from the back.
Figure 0007213491000007

ここで、次式8、9にて表すように、γは伝播定数、αは減衰定数、βは位相定数、fは周波数であるが、本実施形態のアルゴリズムではα=0、及び皮膚8の全微小伝送路51の音速をc=1600(m/s)と仮定している。

Figure 0007213491000008

Figure 0007213491000009
Here, as represented by the following equations 8 and 9, γ is the propagation constant, α is the attenuation constant, β is the phase constant, and f is the frequency. It is assumed that the speed of sound in all microtransmission channels 51 is c=1600 (m/s).
Figure 0007213491000008

Figure 0007213491000009

また、各微小伝送路51の距離Δlは次式10のように表され、ここでも音速をc=1600(m/s)と仮定している。Δtは皮膚8からの反射波形のサンプリング間隔の1ポイントに相当する(本実施形態ではΔt=2(ns))。

Figure 0007213491000010
Also, the distance .DELTA.l of each minute transmission path 51 is represented by the following equation 10, again assuming that the speed of sound is c=1600 (m/s). Δt corresponds to one sampling interval point of the reflected waveform from the skin 8 (Δt=2 (ns) in this embodiment).
Figure 0007213491000010

そして、上記の式をもとに、さらに奥にある微小伝送路51からのインパルス応答Γを求めることができる(次式11)。即ち、Zx及びZの値をもとに、Zの終点におけるΓの値を推定することができる。

Figure 0007213491000011
Then, based on the above equation, the impulse response Γ1 from the micro-transmission path 51 located further inside can be obtained (equation 11 below). That is, the value of Γ 1 at the end point of Z 1 can be estimated based on the values of Zx 0 and Z 1 .
Figure 0007213491000011

次式12は、インパルス応答Γ(ω)を逆フーリエ変換したg(t)を表したものであって、その第1項は多重反射の影響を受けない。従って、この第1項の値から当該微小伝送路51に隣接するさらに奥側の微小伝送路51の特性インピーダンスZを推定することができ、同様にZx、Γ(ω)も推定することができる(次式13,14,15)。

Figure 0007213491000012

Figure 0007213491000013

Figure 0007213491000014

Figure 0007213491000015
The following equation 12 expresses g 1 (t) obtained by inverse Fourier transforming the impulse response Γ 1 (ω), the first term of which is not affected by multiple reflections. Therefore, from the value of the first term, it is possible to estimate the characteristic impedance Z2 of the microtransmission line 51 adjacent to the microtransmission line 51 on the farther side, and similarly estimate Zx1 and Γ2 (ω). (Formulas 13, 14 and 15 below).
Figure 0007213491000012

Figure 0007213491000013

Figure 0007213491000014

Figure 0007213491000015

この工程を繰り返すことによって、各微小伝送路51の特性インピーダンス(固有音響インピーダンス)Z、Z…Zを推定することができる。 By repeating this process, the characteristic impedances (specific acoustic impedances) Z 1 , Z 2 . . . Zn of each minute transmission line 51 can be estimated.

本実施形態のアルゴリズムでは、上記のような推定ステップを行った後、さらに以下に示す所定の補正ステップを実施する。即ち、この補正ステップでは、まず仮想参照部位の固有音響インピーダンス値を定義する。ここでの測定対象物は皮膚8(肌組織)であって(図8参照)、皮膚8を構成する表皮61及び真皮62よりも深い位置にある皮下組織63の比較的浅い部位には、肌組織表面に対して略平行な方向に延びる血管64が存在している。血管64内に存在する血液65の固有音響インピーダンス値は既知(1.61MNs/m)であるため、本実施形態ではこの値を上記「仮想参照部位の固有音響インピーダンス値」として用いている。 In the algorithm of this embodiment, after performing the estimation step as described above, the following predetermined correction step is further performed. That is, in this correction step, first, the intrinsic acoustic impedance value of the virtual reference site is defined. The object to be measured here is the skin 8 (skin tissue) (see FIG. 8). There are blood vessels 64 extending in a direction substantially parallel to the tissue surface. Since the specific acoustic impedance value of the blood 65 present in the blood vessel 64 is known (1.61 MNs/m 3 ), this value is used as the "specific acoustic impedance value of the virtual reference site" in this embodiment.

そして次に、推定ステップにて得た測定対象物内における特定深度(ここでは、横行する血管が存在する皮下組織浅部)の固有音響インピーダンスの推定値を、仮想参照部位の固有音響インピーダンス値に置き換える補正演算を行う。 Next, the estimated value of the specific acoustic impedance at a specific depth (here, shallow subcutaneous tissue where transverse blood vessels exist) in the measurement object obtained in the estimation step is applied to the specific acoustic impedance value of the virtual reference site. Perform correction calculation to replace.

この補正演算を図9のグラフを用いて詳細に説明する。同グラフにおいて縦軸は固有音響インピーダンスを示し、横軸は時間を示している。また、同グラフには、2つの曲線C1、C2が描かれている。実線で描かれた曲線C1は補正前の固有音響インピーダンス曲線C1を示すものであり、破線で描かれた曲線C2は補正後の固有音響インピーダンス曲線C2を示すものである。なお、これら曲線C1、C2は、左半分の領域(即ち樹脂プレート9内部に相当する領域)において差がなく重なっているが、右半分の領域(即ち皮膚8内部に相当する領域)において相違している。本実施形態では、曲線C1の右側領域よりも曲線C2の右側領域のほうが全体的に下方に位置している。 This correction calculation will be described in detail using the graph of FIG. In the graph, the vertical axis indicates the specific acoustic impedance, and the horizontal axis indicates time. Also, two curves C1 and C2 are drawn in the same graph. A curve C1 drawn with a solid line indicates the characteristic acoustic impedance curve C1 before correction, and a curve C2 drawn with a broken line indicates the characteristic acoustic impedance curve C2 after correction. Although these curves C1 and C2 overlap without difference in the left half region (that is, the region corresponding to the inside of the resin plate 9), they are different in the right half region (that is, the region corresponding to the inside of the skin 8). ing. In the present embodiment, the right side area of the curve C2 is positioned generally lower than the right side area of the curve C1.

ここで、基体である樹脂プレート9はポリスチレン製であるため、その固有音響インピーダンス値Z(x)は既知かつ一定(Z(0)=2.50MNs/m)である。固有音響インピーダンス曲線C1における左側の領域が、超音波トランスデューサ13に近い位置にある樹脂プレート9の固有音響インピーダンス値Z(x)を示す部分である。一方、樹脂プレート9に接して配置された皮膚8は、ポリスチレン製の樹脂プレート9よりも軟らかい。このため、その固有音響インピーダンス値Z(x)は、固有音響インピーダンス曲線C1における点Aを境にして急激に低くなる。固有音響インピーダンス曲線C1における右側の領域が皮膚8の固有音響インピーダンス値Z(x)の補正前の推定値を示す部分であり、2.50MNs/mよりもかなり低い値となっている。当該曲線C1において点B(即ち皮膚8内における特定深度に相当する部分)よりも右側はフラットであり、この部分が特定深度の血管内に存在する血液の固有音響インピーダンス値の補正前の推定値を示す部分となっている。ただし、当該推定値は誤差を含んでおり、実際の値とはいくぶん異なるものとして表されている(この例では1.80MNs/m)。そこで、特定深度の血管64内に存在する血液65の固有音響インピーダンス値の補正前の推定値を、仮想参照部位の固有音響インピーダンス値(1.61MNs/m)に置き換えるべく、1.80MNs/mから補正値である0.19MNs/mを減算する。その結果、当該曲線C1のフラットな部分が下方に移動し、これが補正された固有音響インピーダンス曲線C2のうちの一部をなすものとなる。なお、当該曲線C1における点Bは、曲線C2における点Cに移動した状態となる。 Here, since the base resin plate 9 is made of polystyrene, its intrinsic acoustic impedance value Z(x) is known and constant (Z(0)=2.50 MNs/m 3 ). The area on the left side of the characteristic acoustic impedance curve C 1 is a portion that indicates the characteristic acoustic impedance value Z(x) of the resin plate 9 located near the ultrasonic transducer 13 . On the other hand, the skin 8 placed in contact with the resin plate 9 is softer than the polystyrene resin plate 9 . Therefore, the characteristic acoustic impedance value Z(x) sharply decreases across point A on the characteristic acoustic impedance curve C1. The area on the right side of the characteristic acoustic impedance curve C1 is the part showing the pre-correction estimated value of the characteristic acoustic impedance value Z(x) of the skin 8, which is much lower than 2.50 MNs/m 3 . In the curve C1, the right side of the point B (that is, the portion corresponding to the specific depth in the skin 8) is flat, and this portion is the pre-correction estimated value of the characteristic acoustic impedance value of the blood present in the blood vessel at the specific depth. It is the part that shows However, the estimated value contains an error and is presented as somewhat different from the actual value (1.80 MNs/m 3 in this example). Therefore, in order to replace the pre-correction estimated value of the specific acoustic impedance value of the blood 65 present in the blood vessel 64 at a specific depth with the specific acoustic impedance value (1.61 MNs/m 3 ) of the virtual reference site, 1.80 MNs/ A correction value of 0.19 MNs / m3 is subtracted from m3. As a result, the flat portion of the curve C1 moves downwards and forms part of the corrected specific acoustic impedance curve C2. Note that the point B on the curve C1 is moved to the point C on the curve C2.

次に、固有音響インピーダンス曲線C1における点Aと点Bとの間の領域(即ち皮膚表面から当該血管までの部分に相当する領域)を補正する。その際、先の補正値である0.19MNs/mを最大とし、測定対象物の前面側(即ち皮膚表面)に向かうほど小さくなる補正値を用いる。そしてこの補正値を減算することにより、当該血管から皮膚表面までの固有音響インピーダンスの推定値を補正する。換言すると、皮膚表面からの深さに応じた補正値を設定し、推定値からこれを減算する。例えば、特定深度の半分の深さ位置における推定値を補正する場合には、0.19MNs/mの1/2の値(即ち0.095MNs/m)をその位置における補正値として設定し、これをその位置における推定値から減算する。また、特定深度の1/4の深さ位置における推定値を補正する場合には、0.19MNs/mの1/4の値(即ち0.0475MNs/m)をその位置における補正値として設定し、これをその位置における推定値から減算する。ちなみに、図9のグラフに描かれた3本の下向き矢印は、それぞれの深度における補正値の大きさの違いを視覚的に示している。 Next, the area between point A and point B on the characteristic acoustic impedance curve C1 (that is, the area corresponding to the portion from the skin surface to the blood vessel) is corrected. At that time, a correction value is used in which the previous correction value of 0.19 MNs/m 3 is the maximum, and the correction value decreases toward the front side of the measurement object (that is, the skin surface). By subtracting this correction value, the estimated value of the specific acoustic impedance from the blood vessel to the skin surface is corrected. In other words, a correction value corresponding to the depth from the skin surface is set and subtracted from the estimated value. For example, when correcting the estimated value at a depth position that is half the specific depth, set a half value of 0.19 MNs/m 3 (that is, 0.095 MNs/m 3 ) as the correction value at that position. , subtract this from the estimate at that location. Also, when correcting the estimated value at a depth position of 1/4 of the specific depth, the value of 1/4 of 0.19 MNs/m 3 (that is, 0.0475 MNs/m 3 ) is used as the correction value at that position set and subtract this from the estimate at that location. Incidentally, the three downward arrows drawn in the graph of FIG. 9 visually indicate the difference in magnitude of the correction value at each depth.

そして、このような補正ステップの演算を通じて、測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定し、最終的にはBモードエコー像の元となる反射信号列を固有音響インピーダンス像に変換する。 Then, through the calculation of such a correction step, the estimated value of the specific acoustic impedance in the object to be measured is corrected, and then the distribution of the specific acoustic impedance in the depth direction is estimated. is converted into a characteristic acoustic impedance image.

次に、本実施形態の超音波画像構築装置1において固有音響インピーダンス画像を構築するために、プロセッサであるCPU31が実行する演算処理について、図10のフローチャートを用いて説明する。 Arithmetic processing executed by the CPU 31 as a processor in order to construct a specific acoustic impedance image in the ultrasonic image constructing apparatus 1 of the present embodiment will now be described with reference to the flowchart of FIG.

まず、測定対象物であるヒトの皮膚8(例えば、比較的浅い位置に太い血管(けい静脈、けい動脈)が存在している首の皮膚8など)を樹脂プレート9の上面に押し付けるようにして接触配置させる。この状態で、まず超音波プローブ6に初期動作を行わせる。即ち、CPU31からの指示に基づいてコントローラ27を作動させることにより、モータ28X,28Yを駆動し、リファレンス部材10がある位置にて超音波照射が行われるようにX-Yステージ14を移動させる。 First, the human skin 8 (for example, neck skin 8 having relatively shallow thick blood vessels (cavitary veins, carotid arteries), etc.), which is the object to be measured, is pressed against the upper surface of the resin plate 9. Place in contact. In this state, first, the ultrasonic probe 6 is caused to perform an initial operation. That is, by operating the controller 27 based on instructions from the CPU 31, the motors 28X and 28Y are driven, and the XY stage 14 is moved so that ultrasonic waves are applied to the position where the reference member 10 is located.

またこのとき、CPU31からの指示に基づいて励起パルスがトランスデューサ13に供給されると、図5(a)に示すように、リファレンス部材10に超音波Soが照射され、その反射波Sref(ω)が受信回路22を経て検波回路24で検出される。そして、反射波取得手段としてのCPU31は、A/D変換回路25で変換されたデジタルデータをI/F回路32を介して取得し、そのデータをリファレンス部材10からの超音波波形のインパルス応答のデータとしてメモリ33に記憶する(ステップS100)。 Also at this time, when an excitation pulse is supplied to the transducer 13 based on an instruction from the CPU 31, as shown in FIG. passes through the receiving circuit 22 and is detected by the detecting circuit 24 . Then, the CPU 31 as reflected wave acquisition means acquires the digital data converted by the A/D conversion circuit 25 through the I/F circuit 32, and converts the data into an impulse response of the ultrasonic waveform from the reference member 10. It is stored in the memory 33 as data (step S100).

その後、CPU31からの指示に基づいてコントローラ27によりモータ28X,28Yが駆動され、X-Yステージ14による二次元走査が開始される。CPU31は、エンコーダ26の出力に基づいて測定点の座標データを取得する(ステップS110)。 After that, the controller 27 drives the motors 28X and 28Y based on instructions from the CPU 31, and the XY stage 14 starts two-dimensional scanning. The CPU 31 acquires coordinate data of the measurement point based on the output of the encoder 26 (step S110).

そして、図4(a)に示すように、CPU31からの指示に基づいて励起パルスがトランスデューサ13に供給されることにより、皮膚8に超音波Sが照射され、その反射波Stgt(ω)が受信回路22を経て検波回路24で検出される。反射波取得手段としてのCPU31は、A/D変換回路25で変換されたデジタルデータをI/F回路32を介して取得し、そのデータを皮膚8からの超音波波形のインパルス応答のデータとして座標データに関連付けてメモリ33に記憶する(ステップS120)。 Then, as shown in FIG. 4A, an excitation pulse is supplied to the transducer 13 based on an instruction from the CPU 31, so that the skin 8 is irradiated with an ultrasonic wave S o and its reflected wave S tgt (ω) passes through the receiving circuit 22 and is detected by the detecting circuit 24 . The CPU 31 as reflected wave acquisition means acquires the digital data converted by the A/D conversion circuit 25 via the I/F circuit 32, and coordinates the data as impulse response data of the ultrasonic waveform from the skin 8. It is stored in the memory 33 in association with the data (step S120).

次いで、第1演算手段としてのCPU31は、規格化されたインパルス応答信号のデータを用いて、上記のアルゴリズムのうち、TDR法の原理を参考とした推定ステップの演算を実行する。そしてCPU31は、その演算により皮膚8における測定点での奥行方向の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定し、その推定結果を座標データに関連付けてメモリ33に記憶する(ステップS130)。 Next, the CPU 31 as the first computing means uses the data of the normalized impulse response signal to perform the computation of the estimation step referring to the principle of the TDR method among the above algorithms. Then, the CPU 31 sequentially estimates the intrinsic acoustic impedance in the depth direction at the measurement points on the skin 8 from the near side to the far side in the depth direction by the calculation, and stores the estimation results in the memory 33 in association with the coordinate data. (Step S130).

次いで、第2演算手段としてのCPU31は、上記のアルゴリズムのうちの補正ステップの演算を実施する。即ち、仮想参照部位の固有音響インピーダンス値を定義した後、皮膚8内における特定深度の固有音響インピーダンスの推定値を前記仮想参照部位の固有音響インピーダンス値に置き換える演算を行い、この演算を通じて皮膚8内の固有音響インピーダンスの推定値を補正する。そしてCPU31は奥行方向の固有音響インピーダンス分布を推定し、その推定結果を座標データに関連付けてメモリ33に記憶する(ステップS132)。 Next, the CPU 31 as second computing means performs the computation of the correction step in the above algorithm. That is, after defining the specific acoustic impedance value of the virtual reference site, a calculation is performed to replace the estimated value of the specific acoustic impedance at a specific depth in the skin 8 with the specific acoustic impedance value of the virtual reference site. Correct the estimate of the intrinsic acoustic impedance of Then, the CPU 31 estimates the characteristic acoustic impedance distribution in the depth direction, associates the estimation result with the coordinate data, and stores it in the memory 33 (step S132).

その後、画像構築手段としてのCPU31は、奥行方向の固有音響インピーダンス分布の推定結果に基づいて、固有音響インピーダンス像(断層像)を構築するための画像処理を行う(ステップS140)。詳しくは、CPU31は、固有音響インピーダンス分布の推定結果に基づいてカラー変調処理を行い、固有音響インピーダンスの大きさに応じて色分けして表示した画像データを構築し、該画像データをメモリ33に記憶する。 After that, the CPU 31 as image constructing means performs image processing for constructing a specific acoustic impedance image (tomographic image) based on the estimation result of the specific acoustic impedance distribution in the depth direction (step S140). Specifically, the CPU 31 performs color modulation processing based on the estimated result of the characteristic acoustic impedance distribution, constructs image data displayed in different colors according to the magnitude of the characteristic acoustic impedance, and stores the image data in the memory 33. do.

次いで、CPU31は、全ての測定点での処理が終了して、全ての測定点で画像データが取得されたか否かを判断する(ステップS150)。ここで、全データが取得されていない場合には(ステップS150:NO)、CPU31は、ステップS110に戻って、ステップS110~S140の処理を繰り返して実行する。全データが取得された場合には(ステップS150:YES)、CPU31は、次ステップS160に移行する。 Next, the CPU 31 determines whether the processing at all the measurement points has been completed and the image data has been acquired at all the measurement points (step S150). Here, if all the data has not been acquired (step S150: NO), CPU 31 returns to step S110 and repeats the processes of steps S110 to S140. When all the data have been acquired (step S150: YES), the CPU 31 proceeds to the next step S160.

そして、CPU31は、該データを表示装置36に転送し、あらかじめ定めた直線上における固有音響インピーダンス像(断層像)を表示させた後(ステップS160)、図10の処理を終了するようになっている。このような一連の処理により、皮膚8での固有音響インピーダンスの大きさに応じて色分けされた固有音響インピーダンス像(断層像)が表示される。 Then, the CPU 31 transfers the data to the display device 36, displays a specific acoustic impedance image (tomographic image) on a predetermined straight line (step S160), and then terminates the processing of FIG. there is Through such a series of processes, a specific acoustic impedance image (tomographic image) is displayed that is color-coded according to the magnitude of the specific acoustic impedance of the skin 8 .

従って、本実施の形態によれば以下の効果を得ることができる。 Therefore, according to this embodiment, the following effects can be obtained.

(1)本実施形態の超音波画像構築装置1では、リファレンス部材10に入射した超音波波形のインパルス応答情報、及び皮膚8に入射した超音波波形のインパルス応答情報から、規格化されたインパルス応答情報、つまり装置に依存しないインパルス応答情報を得るようにしている。そして、この規格化されたインパルス応答情報に基づいて、皮膚8内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって、多重反射の影響を考慮しつつ順次推定する演算(推定ステップ)を行うようにしている。本実施形態ではこれに加えて、皮膚8内における特定深度(ここでは、横行する血管が存在する皮下組織浅部)の固有音響インピーダンスの推定値を、仮想参照部位の固有音響インピーダンス値に置き換える演算を通じて、皮膚8内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する演算(補正ステップ)も行うようにしている。以上の結果、深度が大きくなるほど蓄積、拡大しやすい固有音響インピーダンスの推定値の誤差を減じることができ、奥行方向の固有音響インピーダンス分布のより正確な推定が可能となる。また、補正ステップを行わない従来技術では超音波断層像における奥行方向に多くの筋が入ってしまうため、微細な層構造が感覚的に理解しにくかったのに対し、本実施形態ではこのような筋が少なくなり層方向の色ムラが減少する。従って、微細な層構造を有する非常に薄い皮膚8の超音波断層像を、その層構造が感覚的に理解しやすい固有音響インピーダンス像として比較的簡単にかつ高い精度で構築することができる。なお、この装置1により得られる固有音響インピーダンス像は、測定対象物を切断することなく(即ち非侵襲で)、層ごとの力学特性の断面分布(深さ分布)情報を、推定した固有音響インピーダンスの絶対値ごとに色分けして画像化したものである。それゆえ、この像は感覚的に層構造が理解しやすいものなっている。 (1) In the ultrasonic image construction apparatus 1 of the present embodiment, the impulse response normalized from the impulse response information of the ultrasonic waveform incident on the reference member 10 and the impulse response information of the ultrasonic waveform incident on the skin 8 Information, ie device independent impulse response information, is obtained. Then, based on the normalized impulse response information, calculation (estimation step) for sequentially estimating the specific acoustic impedance in the skin 8 from the front side to the back side in the depth direction while considering the influence of multiple reflections. I am trying to do In addition to this, in the present embodiment, an operation of replacing the estimated value of the specific acoustic impedance at a specific depth in the skin 8 (here, the superficial subcutaneous tissue where the transverse blood vessel exists) with the specific acoustic impedance value of the virtual reference site. After correcting the estimated value of the specific acoustic impedance in the skin 8 through the skin 8, calculation (correction step) for estimating the specific acoustic impedance distribution in the depth direction is also performed. As a result, it is possible to reduce the error in the estimated value of the characteristic acoustic impedance, which tends to accumulate and expand as the depth increases, and it is possible to more accurately estimate the characteristic acoustic impedance distribution in the depth direction. In addition, in the conventional technique that does not perform the correction step, many streaks appear in the depth direction in the ultrasonic tomographic image, making it difficult to intuitively understand the fine layer structure. Streaks are reduced and color unevenness in the layer direction is reduced. Therefore, an ultrasonic tomographic image of the very thin skin 8 having a fine layer structure can be constructed relatively easily and with high precision as a specific acoustic impedance image whose layer structure is intuitively understandable. The specific acoustic impedance image obtained by this apparatus 1 is obtained by estimating the cross-sectional distribution (depth distribution) information of the mechanical properties for each layer without cutting the measurement object (that is, noninvasively). The image is color-coded for each absolute value of . Therefore, the layer structure of this image is intuitively easy to understand.

ここで、一般に通常の超音波診断装置で得られる超音波Bモードエコー像では、皮膚8等の生体組織内部の層情報は一応得られるものの、得られる画像は固有音響インピーダンスがある程度以上の差を有する層同士の界面からの反射像である。即ち、固有音響インピーダンスの差がある程度小さくなると、組織学的には界面が存在しても検出されず、その構造を反映した画像を形成することが極めて困難となっていた。つまり、一般的な反射像では生体組織の微細な内部構造、微細な層構造が反映された生体組織内部の反射像(固有音響インピーダンスの違い)を把握するのには不十分であった。これに対して、この超音波画像構築装置1によると、従来の超音波Bモードによって全く検出することができなかった力学特性分布に基づく皮膚8の層構造を、十分な解像度を持った鮮明な断層像として捉えることができるようになった。また、このような鮮明な断層像は、他の非侵襲可視化装置(光干渉断層撮影装置(OCT)や、in vivo共焦点顕微鏡など)では取得不能であったため、この超音波画像構築装置1が具現化されたことの意義は大きい。以上のように、本実施形態の超音波画像構築装置1によると、皮膚8の状態(皮膚8の層ごとの力学特性に関する状態)を簡便にかつ非侵襲的に評価することができる。 Here, in an ultrasonic B-mode echo image generally obtained by a normal ultrasonic diagnostic apparatus, although layer information inside a living tissue such as the skin 8 can be obtained, the obtained image has a difference of more than a certain amount in intrinsic acoustic impedance. It is a reflected image from the interface between the layers. That is, when the difference in specific acoustic impedance becomes small to some extent, even if an interface exists, it is not detected histologically, making it extremely difficult to form an image that reflects the structure. In other words, a general reflected image is insufficient for grasping the fine internal structure of the living tissue and the reflected image (difference in specific acoustic impedance) inside the living tissue reflecting the fine layered structure. On the other hand, according to this ultrasonic image construction device 1, the layer structure of the skin 8 based on the mechanical property distribution, which could not be detected at all by the conventional B-mode ultrasonic waves, can be clearly visualized with sufficient resolution. It became possible to capture it as a tomogram. In addition, such a clear tomographic image could not be obtained with other noninvasive visualization devices (optical coherence tomography (OCT), in vivo confocal microscope, etc.), so this ultrasonic image construction device 1 The significance of being embodied is great. As described above, according to the ultrasonic image construction apparatus 1 of the present embodiment, the state of the skin 8 (the state regarding the mechanical properties of each layer of the skin 8) can be easily and noninvasively evaluated.

(2)本実施形態の超音波画像構築装置1では、皮膚組織内部での多重反射を考慮した時間-周波数領域における解析を通じて、超音波Bモードエコー像の元となる反射信号列を奥行方向の固有音響インピーダンス像に変換する手法を採用している。このため、比較的簡単に所定の演算を行うことができる。また本実施形態では、上述したアルゴリズムを用いて、手前側の微小伝送路51の固有音響インピーダンスの推定結果に基づきその奥側に隣接する微小伝送路51の固有音響インピーダンスを推定する処理を順次繰り返すことにより、伝送路の奥行方向の固有音響インピーダンス分布を推定する演算を行うようにしている。従って、この演算によれば、多重反射の影響を最小限にすることができ、奥行方向の固有音響インピーダンス分布を精度よく推定することができる。 (2) In the ultrasonic image construction apparatus 1 of the present embodiment, through analysis in the time-frequency domain considering multiple reflections inside the skin tissue, the reflected signal train that is the basis of the ultrasonic B-mode echo image is obtained in the depth direction. A method of converting to a characteristic acoustic impedance image is adopted. Therefore, the predetermined calculation can be performed relatively easily. Further, in this embodiment, using the algorithm described above, the process of estimating the specific acoustic impedance of the minute transmission line 51 adjacent to the rear side based on the result of estimating the characteristic acoustic impedance of the minute transmission line 51 on the front side is sequentially repeated. By doing so, a calculation for estimating the characteristic acoustic impedance distribution in the depth direction of the transmission path is performed. Therefore, according to this calculation, the influence of multiple reflections can be minimized, and the characteristic acoustic impedance distribution in the depth direction can be accurately estimated.

[第2の実施形態]
以下、本発明の超音波画像構築方法及び装置を具体化した第2の実施形態を図11に基づき詳細に説明する。なお、上記第1の実施形態と共通している構成については、同じ部材番号を付す代わりに詳細な説明を省略する。
[Second embodiment]
A second embodiment embodying the ultrasonic image constructing method and apparatus of the present invention will now be described in detail with reference to FIG. It should be noted that the same member numbers are assigned to the configurations that are common to the first embodiment, and detailed description thereof will be omitted.

本実施形態の超音波画像構築装置1では、所定の推定ステップ後に補正ステップを行ってから画像構築ステップを行っており、この点については上記第1の実施形態と共通する(上記図9、図10等を参照)。ただし、第1の実施形態の補正ステップでは、肌組織表面に対して略平行な方向に延びる血管64内に存在する血液65の固有音響インピーダンス値(1.61MNs/m)を「仮想参照部位の固有音響インピーダンス値」として用いたのに対し、ここでは測定対象物が生体の軟組織であることに着目してその特定深度の固有音響インピーダンスの推定値の平均を「仮想参照部位の固有音響インピーダンス値」と定義して用いている。本実施形態では、肌組織表面から350μmほど深い位置における組織の固有音響インピーダンスの推定値の平均(例えば1.65MNs/m)を求めて、その値を用いている。これは、生体の軟組織中には骨や軟骨などの硬い部分が存在せず硬さのばらつきが少ないので、その固有音響インピーダンスは例えば1.50MNs/m~1.70MNs/mの範囲から大きく外れた値とはならず、仮想参照部位として選択するのに適している、との考え方に立脚したものである。なお、脂肪を除く軟組織であれば、固有音響インピーダンスは、ほぼ1.60MNs/m~1.70MNs/mの範囲内となる。 In the ultrasonic image constructing apparatus 1 of the present embodiment, the correction step is performed after the predetermined estimation step, and then the image constructing step is performed. 10, etc.). However, in the correction step of the first embodiment, the characteristic acoustic impedance value (1.61 MNs/m 3 ) of the blood 65 present in the blood vessel 64 extending in a direction substantially parallel to the skin tissue surface is defined as the “virtual reference site However, here, focusing on the fact that the object to be measured is the soft tissue of a living body, the average of the estimated values of the specific acoustic impedance at a specific depth is taken as the "specific acoustic impedance value of the virtual reference site value”. In this embodiment, an average value (for example, 1.65 MNs/m 3 ) of the estimated values of the intrinsic acoustic impedance of the tissue at a depth of about 350 μm from the surface of the skin tissue is obtained and used. This is because the soft tissues of a living body do not contain hard parts such as bones or cartilage and have little variation in hardness . This is based on the idea that it is suitable for selection as a virtual reference site because the value does not deviate greatly. For soft tissue excluding fat, the specific acoustic impedance is approximately within the range of 1.60 MNs/m 3 to 1.70 MNs/m 3 .

この場合、皮膚8における特定深度は、固有音響インピーダンス像の画像データが構築される領域よりも深い位置に設定されることが好ましい。本実施形態では、例えば肌組織表面から300μmほど深い位置までを画像構築領域として設定している関係上、組織の固有音響インピーダンスの推定値の平均を求めるための特定深度はそれよりも50μm深い位置としている。このようにすることで、特定深度の固有音響インピーダンス像を敢えて表示しないようにし、かえって不自然なものとならないようにしている。 In this case, the specific depth in the skin 8 is preferably set to a position deeper than the region where the image data of the specific acoustic impedance image is constructed. In this embodiment, for example, since the image construction region is set to a position about 300 μm deep from the surface of the skin tissue, the specific depth for obtaining the average of the estimated values of the specific acoustic impedance of the tissue is a position 50 μm deeper than that. and By doing so, the specific acoustic impedance image at a specific depth is intentionally not displayed, so that the image does not look unnatural.

ここで、図11(a)は補正ステップを実施することなく画像を構築して得た固有音響インピーダンス像であり、(b)は補正ステップを実施したうえで画像を構築して得た固有音響インピーダンス像である。前者では、奥行方向に多くの筋が入ってしまい、層構造がわかりにくい。これに対して後者では、筋がいくぶん解消された状態で表示されるため、層方向に色が均一化される傾向が強くなり、層構造がわかりやすくなる。 Here, FIG. 11(a) is the characteristic acoustic impedance image obtained by constructing the image without performing the correction step, and (b) is the characteristic acoustic impedance image obtained by constructing the image after performing the correction step. It is an impedance image. In the former, many streaks appear in the depth direction, making it difficult to understand the layer structure. On the other hand, in the latter, since the streaks are somewhat eliminated, the color tends to be more uniform in the layer direction, making the layer structure easier to understand.

従って、以上説明した本実施形態の超音波画像構築装置1についても、第1の実施形態の超音波画像構築装置1と同様の作用効果を奏することができる。つまり、微細な層構造を有する非常に薄い皮膚8の超音波断層像を、その層構造が感覚的に理解しやすい態様にて比較的簡単にかつ高い精度で構築することができる。また、超音波を利用して得た情報に基づいて、皮膚8の微細な層構造を簡便にかつ正確に観察、評価等できる画像を構築することができる。特に本実施形態によれば、測定部位の内層に第1の実施形態のときのような適当な仮想参照部位(例えばけい静脈、けい動脈などの比較的太い血管)が存在していない場合であっても、固有音響インピーダンスの推定値を補正することができる。よって、首の部位に限らずいろいろな部位(例えば頬など)において観察、評価を行うことができる。 Therefore, the ultrasonic image construction apparatus 1 of this embodiment described above can also achieve the same effects as the ultrasonic image construction apparatus 1 of the first embodiment. That is, it is possible to construct an ultrasonic tomographic image of the very thin skin 8 having a fine layer structure relatively easily and with high accuracy in a manner in which the layer structure is intuitively understandable. Also, based on the information obtained using ultrasound, it is possible to construct an image that allows easy and accurate observation and evaluation of the fine layer structure of the skin 8 . In particular, according to this embodiment, there is no appropriate virtual reference site (for example, relatively thick blood vessels such as cervical veins and arteries) in the inner layer of the measurement site as in the first embodiment. can also correct the estimate of the intrinsic acoustic impedance. Therefore, it is possible to observe and evaluate not only the neck region but also various regions (for example, cheeks).

[第3の実施形態]
以下、本発明の超音波画像構築方法及び装置を具体化した第3の実施形態を図12~図14に基づき詳細に説明する。なお、上記第1の実施形態と共通している構成については、同じ部材番号を付す代わりに詳細な説明を省略する。
[Third embodiment]
A third embodiment embodying the ultrasonic image constructing method and apparatus of the present invention will be described in detail below with reference to FIGS. 12 to 14. FIG. It should be noted that the same member numbers are assigned to the configurations that are common to the first embodiment, and detailed description thereof will be omitted.

図12は、本実施形態の超音波画像構築装置101を示す概略構成図である。図12に示されるように、この超音波画像構築装置101は、超音波を用いて培養細胞8Aを観察するための装置であって、パルス励起型超音波顕微鏡2と、パーソナルコンピュータ(PC)3とを備えている。 FIG. 12 is a schematic configuration diagram showing the ultrasonic image constructing apparatus 101 of this embodiment. As shown in FIG. 12, this ultrasonic image construction device 101 is a device for observing cultured cells 8A using ultrasonic waves, and includes a pulse excitation ultrasonic microscope 2 and a personal computer (PC) 3. and

ステージ4には、測定対象物である培養細胞8Aを培養するための培養液Mを収容する培養容器102が設置されている。ここでは、培養細胞8Aとして、接着性細胞の一種であるヒトのグリア細胞8Aを用いている。グリア細胞8Aの厚さは非常に薄く、数μm程度である。この培養容器102の底部中央部には、上面側に培養細胞8Aが培養時に接着して支持される樹脂プレート9が固定されている。既知の音響物性を有する基体としての樹脂プレート9は、超音波を透過させることができる平板状部材であって、測定対象物である培養細胞8Aよりも硬い材料からなる。 A culture container 102 containing a culture medium M for culturing cultured cells 8A, which are objects to be measured, is installed on the stage 4 . Here, human glial cells 8A, which are a type of adherent cells, are used as the cultured cells 8A. The thickness of the glial cell 8A is very thin, approximately several μm. At the center of the bottom of the culture container 102, a resin plate 9 is fixed on the upper surface side to which the cultured cells 8A are adhered and supported during culture. The resin plate 9 as a substrate having known acoustic properties is a plate-like member that can transmit ultrasonic waves, and is made of a material harder than the cultured cells 8A that are the objects to be measured.

なお、培養容器102内において培養細胞8Aは、培養液Mに完全に浸漬された状態で存在している。従って、培養液Mは、樹脂プレート9において培養細胞8Aが接触配置される側である上面に接して存在するばかりでなく、培養細胞8Aにおける樹脂プレート9とは反対の側にも接して存在していると把握できる。培養液Mは樹脂プレート9とは異なる既知の音響物性を有しており、本実施形態ではこれを第1参照物質及び第2参照部材として位置付けて利用している。以降の説明の便宜のため、樹脂プレート9と接して存在する培養液(即ち第1参照物質である培養液)にM1を付し、培養細胞8Aにおける樹脂プレート9とは反対の側に接して存在する培養液(即ち第2参照物質である培養液)にM2を付して区別する。 It should be noted that the cultured cells 8A exist in the culture container 102 in a state of being completely immersed in the culture medium M. Therefore, the culture medium M is present not only in contact with the upper surface of the resin plate 9 on which the cultured cells 8A are arranged, but also in contact with the side of the cultured cells 8A opposite to the resin plate 9. It can be understood that The culture medium M has known acoustic properties different from those of the resin plate 9, and is used as the first reference material and the second reference member in this embodiment. For the convenience of the following description, the culture solution existing in contact with the resin plate 9 (that is, the culture solution that is the first reference substance) is denoted by M1, and the cultured cell 8A is attached to the side opposite to the resin plate 9. The existing culture medium (ie the culture medium which is the second reference substance) is marked with M2 to distinguish it.

次に、本実施形態の超音波画像構築装置101において、超音波Bモードエコー像の元となる反射信号列から固有音響インピーダンス像を構築する手法について説明する。 Next, a method for constructing a specific acoustic impedance image from a reflected signal train that is the basis of an ultrasonic B-mode echo image in the ultrasonic image constructing apparatus 101 of this embodiment will be described.

この超音波画像構築装置101では、基体である樹脂プレート9に測定対象物である培養細胞8A及び第1参照物質である培養液M1が接して存在し、かつ培養細胞8Aにおける樹脂プレート9とは反対の側に第2参照物質である培養液M2が接して存在した状態で超音波を送信する。そして、樹脂プレート9を介して超音波を入射させたときの培養細胞8A、第1参照物質である培養液M1及び第2参照物質である培養液M2からの超音波波形のインパルス応答をそれぞれ受信する(送受信ステップ)。 In this ultrasonic image construction apparatus 101, the cultured cells 8A that are the objects to be measured and the culture solution M1 that is the first reference substance are in contact with the resin plate 9 that is the base, and the resin plate 9 in the cultured cells 8A is Ultrasonic waves are transmitted while the culture medium M2, which is the second reference substance, is in contact with the opposite side. Then, impulse responses of ultrasonic waveforms from the cultured cells 8A, the culture solution M1 as the first reference substance, and the culture solution M2 as the second reference substance when ultrasonic waves are incident through the resin plate 9 are received respectively. (transmission/reception step).

次いで、培養液M1、M2に入射した超音波波形のインパルス応答情報、及び培養細胞8Aに入射した超音波波形のインパルス応答情報から、規格化されたインパルス応答情報を得て、その規格化されたインパルス応答情報に基づいて培養細胞8A内及び第2参照部材である培養液M2内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する(推定ステップ)。このステップでは、奥行方向の固有音響インピーダンス分布を多重反射の影響を考慮して推定する演算を行うが、これについては第1の実施形態において詳しく説明したため割愛する。なお、この演算は、メモリ33内に格納された所定のアルゴリズムに基づいて実行される。 Next, normalized impulse response information is obtained from the impulse response information of the ultrasonic waveforms incident on the culture solutions M1 and M2 and the impulse response information of the ultrasonic waveforms incident on the cultured cells 8A, and the normalized impulse response information is obtained. Based on the impulse response information, the characteristic acoustic impedance in the cultured cells 8A and in the culture medium M2, which is the second reference member, is sequentially estimated from the near side to the far side in the depth direction (estimating step). In this step, an operation is performed to estimate the characteristic acoustic impedance distribution in the depth direction in consideration of the influence of multiple reflections, but this has been explained in detail in the first embodiment, so its explanation is omitted. Note that this calculation is executed based on a predetermined algorithm stored in the memory 33 .

本実施形態の超音波画像構築装置101では、上記の推定ステップ後に補正ステップを行ってから画像構築ステップを行っており、この点については上記第1の実施形態と共通する(上記図9、図10等を参照)。ただし、本実施形態の補正ステップでは第1の実施形態のときとは異なる以下の補正ステップを実施する。即ちここでは、推定ステップにて得た第2参照物質(培養液M2)の固有音響インピーダンスについて、その「推定値」を「実際値」に置き換える補正演算を行う。培養液M2の前記「実際値」が例えば1.55MNs/mであり、誤差を含んでいる可能性が高い前記「推定値」が例えば1.80MNs/mであるとすると、その値から補正値である0.25MNs/mを減算することにより、前記「実際値」を補正する。その結果、先の実施形態で示した図9のグラフに示す補正前の固有音響インピーダンス曲線C1のフラットな部分が下方に移動し、これが補正された固有音響インピーダンス曲線C2のうちの一部をなすものとなる。なお、当該曲線C1における点Bは、曲線C2における点Cに移動した状態となる。 In the ultrasonic image constructing apparatus 101 of the present embodiment, the correction step is performed after the estimation step, and then the image constructing step is performed. 10, etc.). However, in the correction step of this embodiment, the following correction steps that are different from those in the first embodiment are performed. That is, here, correction calculation is performed to replace the "estimated value" with the "actual value" for the specific acoustic impedance of the second reference substance (culture medium M2) obtained in the estimation step. Assuming that the "actual value" of the culture solution M2 is, for example, 1.55 MNs/m3 , and the "estimated value" that is likely to contain an error is, for example, 1.80 MNs/m3 , from that value The "actual value" is corrected by subtracting a correction value of 0.25 MNs /m3. As a result, the flat part of the characteristic acoustic impedance curve C1 before correction shown in the graph of FIG. 9 shown in the previous embodiment moves downward, forming part of the corrected characteristic acoustic impedance curve C2. become a thing. Note that the point B on the curve C1 is moved to the point C on the curve C2.

次に、固有音響インピーダンス曲線C1における点Aと点Bとの間の領域を補正する。その際、先の補正値である0.25MNs/mを最大とし、測定対象物の前面側(即ち培養細胞8Aにおいて樹脂プレート9と接している面側)に向かうほど小さくなる補正値を用いる。そして前記「推定値」からこの補正値を減算することにより、培養液M2と培養細胞8Aとの界面から、培養細胞8Aと樹脂プレート9との界面までの前記「推定値」を補正する。 Next, the area between points A and B on the characteristic acoustic impedance curve C1 is corrected. At that time, the correction value is set to a maximum value of 0.25 MNs/m 3 , which is the previous correction value, and a correction value that decreases toward the front side of the measurement object (that is, the side of the cultured cell 8A that is in contact with the resin plate 9) is used. . By subtracting this correction value from the "estimated value", the "estimated value" from the interface between the culture fluid M2 and the cultured cells 8A to the interface between the cultured cells 8A and the resin plate 9 is corrected.

そして、このような補正ステップの演算を通じて、測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定し、最終的にはBモードエコー像の元となる反射信号列を固有音響インピーダンス像に変換する。 Then, through the calculation of such a correction step, the estimated value of the specific acoustic impedance in the object to be measured is corrected, and then the distribution of the specific acoustic impedance in the depth direction is estimated. is converted into a characteristic acoustic impedance image.

従って、以上説明した本実施形態の超音波画像構築装置101についても、第1の実施形態の超音波画像構築装置1と同様の作用効果を奏することができる。つまり、微細な内部構造を有する非常に薄い培養細胞8Aの超音波断層像を、その微細な内部構造が感覚的に理解しやすい態様にて比較的簡単にかつ高い精度で構築することができる。また、超音波を利用して得た情報に基づいて、培養細胞8Aの内部構造を簡便にかつ正確に観察、評価等できる画像を構築することができる。 Therefore, the ultrasonic image construction apparatus 101 of this embodiment described above can also achieve the same effects as the ultrasonic image construction apparatus 1 of the first embodiment. That is, it is possible to construct an ultrasonic tomographic image of the very thin cultured cells 8A having fine internal structures relatively easily and with high precision in a manner in which the fine internal structures are intuitively understandable. Also, based on the information obtained using ultrasonic waves, it is possible to construct an image that allows easy and accurate observation and evaluation of the internal structure of the cultured cells 8A.

図13はグリア細胞8Aを概略的に示した図である。グリア細胞8Aには、核112が存在するほか、細胞骨格113(アクチン繊維)が重合して形成される仮足114等が存在している。なお、仮足114は移動するグリア細胞8Aの進行方向側に見られる。図14はこの超音波画像構築装置101を用いて観察されたグリア細胞8Aの超音波断層像である。これによると、グリア細胞8Aの内部には固有音響インピーダンスの高い箇所と低い箇所とが存在していることがわかった。図13の概略図を参照すると、固有音響インピーダンスの高い箇所、即ち比較的硬い箇所は、細胞骨格113を多く含む仮足114であると推察された。また、固有音響インピーダンスの低い箇所、即ち比較的柔らかい箇所は、核であると推察された。このような本実施形態の超音波画像構築装置101によれば、培養細胞8Aの内部構造が観察できるばかりでなく、培養細胞8Aの動き等についても観察できるものと考えられた。 FIG. 13 is a diagram schematically showing glial cells 8A. The glial cell 8A has a nucleus 112 and a pseudopodia 114 formed by polymerization of a cytoskeleton 113 (actin filaments). In addition, the pseudopodia 114 can be seen on the traveling direction side of the migrating glial cells 8A. FIG. 14 is an ultrasonic tomographic image of glial cells 8A observed using this ultrasonic image constructing apparatus 101. FIG. According to this, it was found that there are high and low specific acoustic impedance points inside the glial cells 8A. Referring to the schematic diagram of FIG. 13, it was inferred that the sites with high specific acoustic impedance, ie, relatively hard sites, were pseudopodia 114 containing many cytoskeleton 113 . Also, it was presumed that the portion with low specific acoustic impedance, ie, the relatively soft portion, was the nucleus. According to such an ultrasonic image constructing apparatus 101 of this embodiment, not only the internal structure of the cultured cell 8A can be observed, but also the movement of the cultured cell 8A can be observed.

なお、本発明の各実施の形態は以下のように変更してもよい。 Each embodiment of the present invention may be modified as follows.

・上記第1、第2実施形態の超音波画像構築装置1では、リファレンス部材10からの反射波を参照波形として用いて演算処理を行ったが、これに限定されるものではない。この参照波形としては、樹脂プレート9の上面において皮膚8が接触していない箇所からの反射波であればよく、例えば、樹脂プレート9の上面において皮膚8もリファレンス部材10も接触していない箇所(具体的には、リファレンス部材10よりも外側に位置する樹脂プレート9の表面)での反射波を用いてもよい。言い換えると、樹脂プレート9と空気層との界面からの反射波を参照波形として利用してもよい。 In the ultrasonic image constructing apparatus 1 of the first and second embodiments, arithmetic processing is performed using the reflected wave from the reference member 10 as the reference waveform, but the present invention is not limited to this. The reference waveform may be a reflected wave from a portion of the upper surface of the resin plate 9 that is not in contact with the skin 8. For example, a portion ( Specifically, a reflected wave from the surface of the resin plate 9 located outside the reference member 10 may be used. In other words, the reflected wave from the interface between the resin plate 9 and the air layer may be used as the reference waveform.

・上記第1、第2実施形態の超音波画像構築装置1では、下方から超音波を照射する倒立型の超音波顕微鏡2を用いて超音波の照射を行ったが、上方から超音波を照射する正立型の超音波顕微鏡を用いてもよい。 In the ultrasonic image constructing apparatus 1 of the first and second embodiments, ultrasonic waves are emitted using the inverted ultrasonic microscope 2 that emits ultrasonic waves from below, but ultrasonic waves are emitted from above. An upright ultrasonic microscope may be used.

・上記第1、第2実施形態では、基本的に疾患を有していない比較的健康な皮膚8を対象として、その状態を評価する目的で超音波画像構築装置1を用いたが、これに限定されない。例えば、皮膚がんなどの疾患に伴う皮膚の異常を早期に検出する目的で超音波画像構築装置1を用いてもよい。 - In the above-described first and second embodiments, the ultrasonic image construction apparatus 1 is used for the purpose of evaluating the condition of relatively healthy skin 8 that is basically disease-free. Not limited. For example, the ultrasonic image construction device 1 may be used for the purpose of early detection of skin abnormalities associated with diseases such as skin cancer.

・上記第1、第2実施形態の超音波画像構築装置1では、測定対象物がヒトの首の皮膚8であったが、首以外の他の部位(例えば頬など)における皮膚8であっても勿論よい。また、測定対象物は皮膚8でなくてもよく、例えば内臓、筋肉、脳、歯、爪、骨の表層部などであっても勿論よい。また、上記第3実施形態の超音波画像構築装置101では、測定対象物がヒトのグリア細胞であったが、付着性細胞であれば生物種は問わずどのような培養細胞であってもよい。さらにいうと、測定対象物は必ずしも生体組織や生物でなくてもよく、非生物(例えば塗膜など)であってもよい。換言すると、本発明の超音波画像構築装置1、101は医療分野、美容分野、化粧品分野のみに限定されず、例えば工業分野などの分野においても使用されることができる。 In the ultrasonic image construction apparatus 1 of the first and second embodiments, the measurement object was the skin 8 of the human neck, but the skin 8 in other parts other than the neck (such as cheeks) is of course good. Moreover, the object to be measured may not be the skin 8, but may of course be internal organs, muscles, brains, teeth, nails, surface layers of bones, and the like. In addition, in the ultrasonic image construction apparatus 101 of the third embodiment, the object to be measured was human glial cells, but any cultured cells of any species can be used as long as they are adherent cells. . Furthermore, the object to be measured does not necessarily have to be a living tissue or living organism, and may be a non-living organism (for example, a coating film). In other words, the ultrasonic image construction apparatus 1, 101 of the present invention is not limited to the medical field, beauty field, and cosmetic field, but can also be used in the industrial field, for example.

・上記実施形態の超音波画像構築装置1、101は、いずれも測定対象物に対して超音波トランスデューサ13を二次元方向に相対的に走査させる走査手段を備えていたが、これに代えて超音波トランスデューサ13を一次元方向にのみ相対的に走査させる走査手段を備えたものとしてもよい。また、走査手段は必須の構成ではないため省略しても勿論よく、この場合には装置を小型化、簡略化、低コスト化することが可能となる。 The ultrasonic image constructing apparatuses 1 and 101 of the above-described embodiments both have scanning means for relatively scanning the measurement object with the ultrasonic transducer 13 in two-dimensional directions. Scanning means may be provided for relatively scanning the acoustic wave transducer 13 only in one-dimensional direction. Further, since the scanning means is not an essential component, it may of course be omitted. In this case, the size and cost of the apparatus can be reduced.

・上記実施形態の超音波画像構築装置1、101では、奥行方向の固有音響インピーダンス分布の推定結果に基づいて固有音響インピーダンス像を構築したが、これに限定されない。例えば、奥行方向の音速分布を推定し、その結果に基づいて音速像を構築してもよい。 In the ultrasonic image constructing apparatuses 1 and 101 of the above-described embodiments, the characteristic acoustic impedance image is constructed based on the estimation result of the characteristic acoustic impedance distribution in the depth direction, but the present invention is not limited to this. For example, a sound velocity distribution in the depth direction may be estimated, and a sound velocity image may be constructed based on the result.

・上記実施形態の超音波画像構築装置1、101では、超音波Bモードエコー画像の元となる反射信号列から固有音響インピーダンス像を構築してそれを表示装置36に表示させるように構成したが、固有音響インピーダンス像ばかりでなく超音波Bモードエコー像も表示できるようにしても勿論よい。また、超音波Bモードエコー像を表示する汎用の超音波診断装置に上記実施形態のアルゴリズムを組み込むことで、超音波画像構築装置1、101として動作させるようにしてもよい。 In the ultrasonic image constructing apparatuses 1 and 101 of the above-described embodiments, a specific acoustic impedance image is constructed from the reflected signal train that is the basis of the ultrasonic B-mode echo image, and the image is displayed on the display device 36. Of course, not only the characteristic acoustic impedance image but also the ultrasonic B-mode echo image may be displayed. Further, by incorporating the algorithm of the above-described embodiment into a general-purpose ultrasonic diagnostic apparatus that displays an ultrasonic B-mode echo image, the apparatus may be operated as the ultrasonic image constructing apparatus 1 or 101 .

・上記実施形態では、測定対象物内での多重反射を考慮した時間-周波数領域における解析を通じて、超音波Bモードエコー像の元となる反射信号列を奥行方向の固有音響インピーダンス像に変換する手法を採用したが、これに限定されない。反射波形から固有音響インピーダンスの分布を推定するための方法としては、例えば、伝送路内部での多重反射を含む全ての反射経路を想定して時間軸上で随時応答を解析していく手法を採用してもよい。この手法であっても、上記実施形態の手法のときと同様の固有音響インピーダンス分布の推定結果を得ることができる。 ・In the above embodiment, a method of converting a reflected signal train, which is the basis of an ultrasonic B-mode echo image, into a specific acoustic impedance image in the depth direction through analysis in the time-frequency domain that considers multiple reflections within the measurement object. was adopted, but is not limited to this. As a method for estimating the distribution of the characteristic acoustic impedance from the reflected waveform, for example, a method is adopted in which all reflection paths, including multiple reflections inside the transmission path, are assumed and the response is analyzed at any time on the time axis. You may Even with this method, it is possible to obtain the same estimation result of the characteristic acoustic impedance distribution as with the method of the above embodiment.

・上記実施形態では、超音波Bモードエコー像の元となる反射信号列を奥行方向の固有音響インピーダンス像に変換するにあたり、測定対象物内での多重反射を考慮した解析手法を採用したが、これに限定されない。例えば、生体の軟組織のように多重反射の影響が小さいと考えられる場合、測定対象物内での多重反射を敢えて考慮しない解析手法を採用してもよい。 In the above embodiment, an analysis method that considers multiple reflections within the object to be measured is used to convert the reflected signal train, which is the basis of the ultrasonic B-mode echo image, into the characteristic acoustic impedance image in the depth direction. It is not limited to this. For example, when the influence of multiple reflections is considered to be small, such as soft tissue of a living body, an analysis method that does not consider multiple reflections within the measurement object may be employed.

1、101…超音波画像構築装置
8…測定対象物としての生体の軟組織(皮膚)
8A…測定対象物としての培養細胞(グリア細胞)
9…基体としての樹脂プレート
10…第1参照物質としてのリファレンス部材
13…超音波振動子としての超音波トランスデューサ
31…第1演算手段、第2演算手段、画像構築手段、プロセッサとしてのCPU
51…微小伝送路
64…血管
65…血液
M1…第1参照物質としての培養液
M2…第2参照物質としての培養液
Γ、Γ…インパルス応答
Reference Signs List 1, 101: Ultrasonic image construction device 8: Soft tissue (skin) of a living body as an object to be measured
8A... Cultured cells (glial cells) as measurement objects
9 Resin plate as base 10 Reference member as first reference material 13 Ultrasonic transducer as ultrasonic transducer 31 First calculation means, second calculation means, image constructing means, CPU as processor
51 Microtransmission channel 64 Blood vessel 65 Blood M1 Culture medium M2 as first reference material Culture medium Γ 0 , Γ 1 as second reference material Impulse response

Claims (15)

基体に測定対象物及び第1参照物質が接して存在し、かつ前記測定対象物における前記基体とは反対の側に第2参照物質が接して存在した状態で超音波を送信し、前記基体を介して超音波を入射させたときの前記測定対象物、前記第1参照物質及び前記第2参照物質からの超音波波形のインパルス応答を受信する送受信ステップと、
前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内及び前記第2参照物質内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する推定ステップと、
前記推定ステップにて得た前記第2参照物質の固有音響インピーダンスの推定値を前記第2参照物質の固有音響インピーダンスの実際値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する補正ステップと、
前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築ステップと
を有することを特徴とする超音波画像構築方法。
An ultrasonic wave is transmitted in a state in which an object to be measured and a first reference substance are in contact with a substrate, and a second reference substance is in contact with the side of the object to be measured opposite to the substrate, and the substrate is a transmitting/receiving step of receiving impulse responses of ultrasonic waveforms from the measurement object, the first reference material, and the second reference material when ultrasonic waves are incident through the
Based on the normalized impulse response information obtained from the impulse response information of the ultrasonic waveform, the intrinsic acoustic impedance in the measurement object and in the second reference material is sequentially measured from the near side to the far side in the depth direction. an estimation step to estimate;
The estimated value of the specific acoustic impedance of the object to be measured is obtained by replacing the estimated value of the specific acoustic impedance of the second reference material obtained in the estimation step with the actual value of the specific acoustic impedance of the second reference material. a correction step of estimating the characteristic acoustic impedance distribution in the depth direction after correction;
and an image construction step of constructing image data of an acoustic physical property image based on the characteristic acoustic impedance distribution in the depth direction.
基体に測定対象物及び第1参照物質が接して存在した状態で超音波を送信し、前記基体を介して超音波を入射させたときの前記測定対象物及び前記第1参照物質からの超音波波形のインパルス応答を受信する送受信ステップと、
前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する推定ステップと、
前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値、または、前記測定対象物が生体の軟組織である場合の前記測定対象物内における特定深度の固有音響インピーダンスの推定値の平均を仮想参照部位の固有音響インピーダンス値と定義し、前記推定ステップにて得た前記測定対象物内における前記特定深度の固有音響インピーダンスの推定値を前記仮想参照部位の固有音響インピーダンス値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する補正ステップと、
前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築ステップと
を有することを特徴とする超音波画像構築方法。
Ultrasonic waves from the measurement object and the first reference substance when the ultrasonic waves are transmitted while the measurement object and the first reference substance are in contact with the substrate, and the ultrasonic waves are incident through the substrate. a transmitting/receiving step of receiving an impulse response of the waveform;
an estimation step of sequentially estimating the intrinsic acoustic impedance in the measurement object from the near side to the far side in the depth direction based on the normalized impulse response information obtained from the impulse response information of the ultrasonic waveform;
A specific acoustic impedance value when a substance having a known specific acoustic impedance value exists uniformly in the layer direction at a specific depth in the measurement object, or when the measurement object is a soft tissue of a living body The average of the estimated values of the specific acoustic impedance at a specific depth in the measurement object is defined as the specific acoustic impedance value of the virtual reference part, and the specific acoustic impedance at the specific depth in the measurement object obtained in the estimation step a correcting step of estimating a specific acoustic impedance distribution in the depth direction after correcting the estimated value of the specific acoustic impedance in the measurement object through a calculation that replaces the estimated value of the impedance with the specific acoustic impedance value of the virtual reference part; ,
and an image construction step of constructing image data of an acoustic physical property image based on the characteristic acoustic impedance distribution in the depth direction.
前記補正ステップでは、前記第2参照物質の固有音響インピーダンスの実際値または前記仮想参照部位の固有音響インピーダンス値に置き換えるときの値を最大とし、前記測定対象物の前面側に向かうほど小さくなる補正値を加算または減算することにより、前記測定対象物内の固有音響インピーダンスの推定値を補正することを特徴とする請求項1または2に記載の超音波画像構築方法。 In the correction step, the correction value is set such that the actual value of the specific acoustic impedance of the second reference material or the value when replacing with the specific acoustic impedance value of the virtual reference portion is maximized, and the correction value decreases toward the front side of the measurement object. 3. The method of constructing an ultrasound image according to claim 1, wherein the estimated value of the intrinsic acoustic impedance within the object is corrected by adding or subtracting . 前記推定ステップでは、前記測定対象物内において、異なる固有音響インピーダンスを持つ無損失の微小伝送路が奥行方向に連なって伝送路の集合体をなしていると仮定して、手前側の前記微小伝送路の固有音響インピーダンスの推定結果に基づきその奥側に隣接する前記微小伝送路の固有音響インピーダンスを推定する演算を行うことを特徴とする請求項1乃至3のいずれか1項に記載の超音波画像構築方法。 In the estimating step, it is assumed that lossless minute transmission lines having different characteristic acoustic impedances are connected in the depth direction to form a collection of transmission lines in the object to be measured. 4. The ultrasonic wave according to any one of claims 1 to 3, wherein an operation for estimating the specific acoustic impedance of said micro-transmission channel adjacent to the back side thereof is performed based on the result of estimating the specific acoustic impedance of the channel. Image construction method. 前記測定対象物は、培養細胞であり、
前記第1参照物質及び前記第2参照物質は、前記培養細胞を培養するための培養液であり、
前記基体は、前記培養液を収容する培養容器の一部である
ことを特徴とする請求項に記載の超音波画像構築方法。
The object to be measured is a cultured cell,
The first reference substance and the second reference substance are culture solutions for culturing the cultured cells,
2. The method of constructing an ultrasonic image according to claim 1 , wherein the substrate is a part of a culture vessel containing the culture medium.
前記測定対象物は、組織表面に対して略平行な方向に延びる血管を含む生体の軟組織としての皮膚であり、
前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値は、前記血管中に存在する血液の固有音響インピーダンス値である
ことを特徴とする請求項に記載の超音波画像構築方法。
The object to be measured is the skin as a soft tissue of a living body containing blood vessels extending in a direction substantially parallel to the tissue surface,
When a substance having a known specific acoustic impedance value exists uniformly in the layer direction at a specific depth in the measurement object, the specific acoustic impedance value is the specific acoustic impedance value of blood present in the blood vessel. 3. The method of constructing an ultrasound image according to claim 2 , wherein:
前記測定対象物が生体の軟組織である場合の前記測定対象物内における特定深度の固有音響インピーダンスの推定値の平均は、前記測定対象物が皮膚である場合の当該皮膚内における特定深度の固有音響インピーダンスの推定値の平均であり、
前記特定深度は、前記音響物性像の画像データが構築される領域よりも深い位置に設定される
ことを特徴とする請求項に記載の超音波画像構築方法。
When the measurement object is the soft tissue of a living body, the average of the estimated values of the specific acoustic impedance at the specific depth in the measurement object is the specific acoustic impedance at the specific depth in the skin when the measurement object is the skin. is the average of the impedance estimates,
3. The ultrasonic image constructing method according to claim 2 , wherein the specific depth is set to a position deeper than a region in which the image data of the acoustic physical property image is constructed.
基体と、
前記基体に測定対象物及び第1参照物質が接して存在し、かつ前記測定対象物における前記基体とは反対の側に第2参照物質が接して存在した状態で超音波を送信し、前記基体を介して超音波を入射させたときの前記測定対象物、前記第1参照物質及び前記第2参照物質からの超音波波形のインパルス応答を受信する超音波振動子と、
前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内及び前記第2参照物質内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する第1演算手段と、
前記第1演算手段により得た前記第2参照物質の固有音響インピーダンスの推定値を前記第2参照物質の固有音響インピーダンスの実際値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する第2演算手段と、
前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築手段と
を備えたことを特徴とする超音波画像構築装置。
a substrate;
An ultrasonic wave is transmitted in a state in which an object to be measured and a first reference substance are in contact with the base, and a second reference substance is in contact with the side of the object to be measured opposite to the base, and the base is an ultrasonic transducer that receives impulse responses of ultrasonic waveforms from the measurement object, the first reference material, and the second reference material when ultrasonic waves are incident via
Based on the normalized impulse response information obtained from the impulse response information of the ultrasonic waveform, the intrinsic acoustic impedance in the measurement object and in the second reference material is sequentially measured from the near side to the far side in the depth direction. a first computing means for estimating;
Estimated value of specific acoustic impedance in the object to be measured through calculation of replacing the estimated value of specific acoustic impedance of the second reference material obtained by the first calculating means with the actual value of specific acoustic impedance of the second reference material. a second computing means for estimating the characteristic acoustic impedance distribution in the depth direction after correcting the
An ultrasonic image constructing apparatus, comprising: image constructing means for constructing image data of an acoustic physical property image based on the characteristic acoustic impedance distribution in the depth direction.
基体と、
前記基体に測定対象物及び第1参照物質が接して存在した状態で超音波を送信し、前記基体を介して超音波を入射させたときの前記測定対象物及び前記第1参照物質からの超音波波形のインパルス応答を受信する超音波振動子と、
前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する第1演算手段と、
前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値、または、前記測定対象物が生体の軟組織である場合の前記測定対象物内における特定深度の固有音響インピーダンスの推定値の平均を仮想参照部位の固有音響インピーダンス値と定義し、前記第1演算手段により得た前記測定対象物内における前記特定深度の固有音響インピーダンスの推定値を前記仮想参照部位の固有音響インピーダンス値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する第2演算手段と、
前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築手段と
を備えたことを特徴とする超音波画像構築装置。
a substrate;
Ultrasonic waves are transmitted while the object to be measured and the first reference substance are in contact with the substrate, and the ultrasonic waves from the object to be measured and the first reference substance when the ultrasonic waves are incident through the substrate an ultrasonic transducer that receives an impulse response of a sound wave waveform;
a first computing means for sequentially estimating the intrinsic acoustic impedance in the object to be measured from the near side to the far side in the depth direction based on the normalized impulse response information obtained from the impulse response information of the ultrasonic waveform; ,
A specific acoustic impedance value when a substance having a known specific acoustic impedance value exists uniformly in the layer direction at a specific depth in the measurement object, or when the measurement object is a soft tissue of a living body defined as the intrinsic acoustic impedance value of the virtual reference part, and the intrinsic acoustic impedance at the specific depth in the measurement object obtained by the first computing means A second method for estimating a specific acoustic impedance distribution in the depth direction after correcting the estimated value of the specific acoustic impedance in the measurement object through a calculation for replacing the estimated value of the acoustic impedance with the specific acoustic impedance value of the virtual reference part. computing means;
An ultrasonic image constructing apparatus, comprising: image constructing means for constructing image data of an acoustic physical property image based on the characteristic acoustic impedance distribution in the depth direction.
前記第2演算手段は、前記第2参照物質の固有音響インピーダンスの実際値または前記仮想参照部位の固有音響インピーダンス値に置き換えるときの値を最大し、前記測定対象物の前面側に向かうほど小さくなる補正値を加算または減算することにより、前記測定対象物内の固有音響インピーダンスの推定値を補正することを特徴とする請求項8または9に記載の超音波画像構築装置。 The second computing means maximizes the actual value of the specific acoustic impedance of the second reference material or the value when replacing with the specific acoustic impedance value of the virtual reference part, and decreases toward the front side of the measurement object. 10. The ultrasonic image construction apparatus according to claim 8 or 9, wherein the estimated value of the intrinsic acoustic impedance within the measurement object is corrected by adding or subtracting a correction value. 前記第1演算手段は、前記測定対象物内において、異なる固有音響インピーダンスを持つ無損失の微小伝送路が奥行方向に連なって伝送路の集合体をなしていると仮定して、手前側の前記微小伝送路の固有音響インピーダンスの推定結果に基づきその奥側に隣接する前記微小伝送路の固有音響インピーダンスを推定する演算を行うことを特徴とする請求項8乃至10のいずれか1項に記載の超音波画像構築装置。 The first computing means assumes that in the object to be measured, lossless minute transmission lines having different characteristic acoustic impedances are connected in the depth direction to form a collection of transmission lines, and the front side of the 11. The method according to any one of claims 8 to 10, wherein calculation is performed for estimating the specific acoustic impedance of said micro-transmission channel adjacent on the far side based on the result of estimating the specific acoustic impedance of the micro-transmission channel. Ultrasound imaging equipment. プロセッサに、
基体に測定対象物及び第1参照物質が接して存在し、かつ前記測定対象物における前記基体とは反対の側に第2参照物質が接して存在した状態で超音波振動子に超音波を送信させ、前記基体を介して超音波を入射させたときの前記測定対象物、前記第1参照物質及び前記第2参照物質からの超音波波形のインパルス応答を前記超音波振動子に受信させる送受信ステップと、
前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内及び前記第2参照物質内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する推定ステップと、
前記推定ステップにて得た前記第2参照物質の固有音響インピーダンスの推定値を前記第2参照物質の固有音響インピーダンスの実際値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する補正ステップと、
前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築ステップと
を実行させるための超音波画像構築プログラム。
to the processor,
An ultrasonic wave is transmitted to an ultrasonic transducer in a state in which an object to be measured and a first reference substance are in contact with a substrate, and a second reference substance is in contact with the side of the object to be measured opposite to the substrate. and a transmitting/receiving step of causing the ultrasonic transducer to receive an impulse response of an ultrasonic waveform from the object to be measured, the first reference substance, and the second reference substance when ultrasonic waves are incident through the substrate. When,
Based on the normalized impulse response information obtained from the impulse response information of the ultrasonic waveform, the intrinsic acoustic impedance in the measurement object and in the second reference material is sequentially measured from the near side to the far side in the depth direction. an estimation step to estimate;
The estimated value of the specific acoustic impedance of the object to be measured is obtained by replacing the estimated value of the specific acoustic impedance of the second reference material obtained in the estimation step with the actual value of the specific acoustic impedance of the second reference material. a correction step of estimating the characteristic acoustic impedance distribution in the depth direction after correction;
An ultrasonic image construction program for executing an image construction step of constructing image data of an acoustic physical property image based on the characteristic acoustic impedance distribution in the depth direction.
プロセッサに、
基体に測定対象物及び第1参照物質が接して存在した状態で超音波振動子に超音波を送信させ、前記基体を介して超音波を入射させたときの前記測定対象物及び前記第1参照物質からの超音波波形のインパルス応答を前記超音波振動子に受信させる送受信ステップと、
前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する推定ステップと、
前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値、または、前記測定対象物が生体の軟組織である場合の前記測定対象物内における特定深度の固有音響インピーダンスの推定値の平均を仮想参照部位の固有音響インピーダンス値と定義し、前記推定ステップにて得た前記測定対象物内における前記特定深度の固有音響インピーダンスの推定値を前記仮想参照部位の固有音響インピーダンス値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する補正ステップと、
前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築ステップと
を実行させるための超音波画像構築プログラム。
to the processor,
The object to be measured and the first reference when the object to be measured and the first reference substance exist in contact with the substrate and the ultrasonic wave is transmitted through the ultrasonic transducer and the ultrasonic wave is incident through the substrate. a transmitting/receiving step of causing the ultrasonic transducer to receive an impulse response of an ultrasonic waveform from a substance;
an estimation step of sequentially estimating the intrinsic acoustic impedance in the measurement object from the near side to the far side in the depth direction based on the normalized impulse response information obtained from the impulse response information of the ultrasonic waveform;
A specific acoustic impedance value when a substance having a known specific acoustic impedance value exists uniformly in the layer direction at a specific depth in the measurement object, or when the measurement object is a soft tissue of a living body The average of the estimated values of the specific acoustic impedance at a specific depth in the measurement object is defined as the specific acoustic impedance value of the virtual reference part, and the specific acoustic impedance at the specific depth in the measurement object obtained in the estimation step a correcting step of estimating a specific acoustic impedance distribution in the depth direction after correcting the estimated value of the specific acoustic impedance in the measurement object through a calculation that replaces the estimated value of the impedance with the specific acoustic impedance value of the virtual reference part; ,
An ultrasonic image construction program for executing an image construction step of constructing image data of an acoustic physical property image based on the characteristic acoustic impedance distribution in the depth direction.
前記補正ステップでは、前記第2参照物質の固有音響インピーダンスの実際値または前記仮想参照部位の固有音響インピーダンス値に置き換えるときの値を最大とし、前記測定対象物の前面側に向かうほど小さくなる補正値を加算または減算することにより、前記測定対象物内の固有音響インピーダンスの推定値を補正することを特徴とする請求項12または13に記載の超音波画像構築プログラム。 In the correction step, the correction value is set such that the actual value of the specific acoustic impedance of the second reference material or the value when replacing with the specific acoustic impedance value of the virtual reference portion is maximized, and the correction value decreases toward the front side of the measurement object. 14. The ultrasonic image construction program according to claim 12 or 13, wherein the estimated value of the intrinsic acoustic impedance in the object to be measured is corrected by adding or subtracting . 前記推定ステップでは、前記測定対象物内において、異なる固有音響インピーダンスを持つ無損失の微小伝送路が奥行方向に連なって伝送路の集合体をなしていると仮定して、手前側の前記微小伝送路の固有音響インピーダンスの推定結果に基づきその奥側に隣接する前記微小伝送路の固有音響インピーダンスを推定する演算を行うことを特徴とする請求項12乃至14のいずれか1項に記載の超音波画像構築プログラム。 In the estimating step, it is assumed that lossless minute transmission lines having different characteristic acoustic impedances are connected in the depth direction to form a collection of transmission lines in the object to be measured. 15. The ultrasonic wave according to any one of claims 12 to 14, wherein an operation is performed for estimating the specific acoustic impedance of said micro-transmission channel adjacent on the far side based on the result of estimating the specific acoustic impedance of the channel. Image construction program.
JP2019095229A 2019-05-21 2019-05-21 Ultrasound Image Construction Method, Ultrasound Image Construction Apparatus, Ultrasound Image Construction Program Active JP7213491B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2019095229A JP7213491B2 (en) 2019-05-21 2019-05-21 Ultrasound Image Construction Method, Ultrasound Image Construction Apparatus, Ultrasound Image Construction Program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019095229A JP7213491B2 (en) 2019-05-21 2019-05-21 Ultrasound Image Construction Method, Ultrasound Image Construction Apparatus, Ultrasound Image Construction Program

Publications (2)

Publication Number Publication Date
JP2020190454A JP2020190454A (en) 2020-11-26
JP7213491B2 true JP7213491B2 (en) 2023-01-27

Family

ID=73453629

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019095229A Active JP7213491B2 (en) 2019-05-21 2019-05-21 Ultrasound Image Construction Method, Ultrasound Image Construction Apparatus, Ultrasound Image Construction Program

Country Status (1)

Country Link
JP (1) JP7213491B2 (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000005180A (en) 1998-06-25 2000-01-11 Olympus Optical Co Ltd Acoustic impedance measuring device
WO2008069250A1 (en) 2006-12-08 2008-06-12 Nihon University Device for measuring physical property of cell
JP2011172611A (en) 2010-02-23 2011-09-08 Canon Inc Ultrasonic imaging device and delay control method
WO2018207276A1 (en) 2017-05-10 2018-11-15 本多電子株式会社 Ultrasound image construction method, ultrasound image construction apparatus, ultrasound image construction program, and method for evaluating skin

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6083645A (en) * 1983-10-17 1985-05-11 テルモ株式会社 Ultrasonic tomographic measuring method and apparatus
JPH02215449A (en) * 1989-02-16 1990-08-28 Fujitsu Ltd Ultrasonic diagnostic device

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000005180A (en) 1998-06-25 2000-01-11 Olympus Optical Co Ltd Acoustic impedance measuring device
WO2008069250A1 (en) 2006-12-08 2008-06-12 Nihon University Device for measuring physical property of cell
JP2011172611A (en) 2010-02-23 2011-09-08 Canon Inc Ultrasonic imaging device and delay control method
WO2018207276A1 (en) 2017-05-10 2018-11-15 本多電子株式会社 Ultrasound image construction method, ultrasound image construction apparatus, ultrasound image construction program, and method for evaluating skin

Also Published As

Publication number Publication date
JP2020190454A (en) 2020-11-26

Similar Documents

Publication Publication Date Title
EP2926739B1 (en) Acquisition control for elasticity ultrasound imaging
KR102170630B1 (en) Measuring acoustic absorption or attenuation of ultrasound
JP6169707B2 (en) Ultrasonic diagnostic apparatus and elasticity evaluation method
CN102458260B (en) Measure the method and apparatus of cardiac systolic function
Lawrence Physics and instrumentation of ultrasound
JP3723663B2 (en) Ultrasonic diagnostic equipment
JP6361001B1 (en) Ultrasound image construction method, ultrasound image construction device, ultrasound image construction program, skin evaluation method for cosmetic purposes
CN103648400B (en) Diagnostic ultrasound equipment and method
JP6053512B2 (en) Subject information acquisition apparatus and control method thereof
US20160192906A1 (en) Method for producing elastic image and ultrasonic diagnostic apparatus
JP2017104526A (en) Analysis device
KR20140086626A (en) Method for measuring the displacement of shear wave and mechanical parameters in tissue by using shear wave and the system comprising the same
JP5588924B2 (en) Ultrasonic diagnostic equipment
JP2017006213A (en) Ultrasonic diagnostic apparatus and control program
JP6257892B2 (en) Ultrasonic diagnostic apparatus and control program
JP7213491B2 (en) Ultrasound Image Construction Method, Ultrasound Image Construction Apparatus, Ultrasound Image Construction Program
JP7539641B2 (en) Ultrasound imaging method and apparatus, and signal processing method
JP2005324072A (en) Ultrasonic diagnostic apparatus
JP6081744B2 (en) Ultrasonic diagnostic equipment
JP7295534B2 (en) Ultrasound image construction method, device and program, and signal processing method
JP2014068755A (en) Ultrasonic inspection apparatus, signal processing method of ultrasonic inspection apparatus, and program
JP2018199002A (en) Information processing device and control method thereof
JP2017039031A (en) Subject information acquisition device and control method thereof
JP2024078031A (en) Ultrasonic image generation method and device, and signal processing method
JP2020065799A (en) Ultrasound diagnostic apparatus

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20220215

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220215

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220316

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20221129

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20221130

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20221208

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230106

R150 Certificate of patent or registration of utility model

Ref document number: 7213491

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150