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 PDFInfo
- 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
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.
超音波Ｂモードエコー像は、医療分野において広く一般的に使われている方法であり、従来このような像を得るための装置が多数提案されている（特許文献１を参照）。簡単に説明すると、超音波Ｂモードエコー像とは、物体に入射した超音波が反射して返って来る時の反射信号列を画像化したものである。超音波が散乱せずに真っ直ぐ進んだと仮定する場合、電気信号と同じく、進んだ先の抵抗値（固有音響インピーダンス）の違いによって反射が生じることになる。それゆえ、固有音響インピーダンスの分布が分かれば、どのような反射信号列が返って来るかを推測することが可能となる。即ち、音響物性分布が既知であれば、どのようなＢモード画像が観察されるかを推測することができる。また、その逆もしかりである。 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.
しかしながら、生体組織などの不均一で厚み（深さ）があるターゲットについては、入射して返って来る反射波形は、ターゲットに入射した超音波が散乱、吸収を経て様々な進み方をした結果と、多重反射をした結果とを反映したものとなる。このため、反射波形を固有音響インピーダンス等の音響物性に変換することは困難であると考えられており、この方法は従来検討されてこなかった。さらに、超音波Ｂモードエコー像は、生体組織内部での超音波の多重反射に起因したスペックルノイズによって乱れた画像となりやすいため、内部構造を高い精度で表示するのには不向きであるという問題もあった。ゆえに、従来装置では音響フィルタを入れるなどの対策が必要となり、構成が複雑化する等の問題があった。 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.
また、超音波Ｂモードエコー像を表示する通常の超音波診断装置では、皮膚等の生体組織内部の層情報は一応得られるものの、得られる画像は固有音響インピーダンスの異なる層同士の界面からの反射像である。よって、このような反射像では生体組織の内部構造、具体的には生体組織内部の固有音響インピーダンスの違いを把握するのには不十分であった。つまり、従来技術で得られる反射像は、層の界面の位置がどこにあるか等については感覚的に理解しやすいものである反面、界面と界面とで囲まれた中間の領域の固有音響インピーダンスがどのようになっているかが感覚的に理解しにくいものであった。ゆえに、超音波を利用して得た情報に基づいて、層構造を有する非常に薄い測定対象物の超音波断層像を、感覚的に層構造が理解しやすい態様にて構築することが望まれていた。 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
このような事情に鑑みて本願発明者らは、改良された超音波画像構築装置をすでに提案している（例えば、特許文献２を参照）。この装置は、既知の音響物性を有する基体、基材を介して超音波の送受信を行う超音波振動子、演算手段、画像構築手段等を含んで構成されている。この装置では、既知の音響物性を有する基体に、測定対象物及び既知の音響物性を有する参照物質を接して配置する。そしてこの状態で超音波を送信し、基体を介して測定対象物及び参照物質に超音波を入射させ、測定対象物及び参照物質からの超音波波形のインパルス応答を受信する。次いで、参照物質に入射した超音波波形のインパルス応答情報及び測定対象物に入射した超音波波形のインパルス応答情報から、規格化されたインパルス応答情報を得る。この規格化されたインパルス応答情報に基づき、奥行方向の音響物性分布（具体的には固有音響インピーダンス分布）を多重反射の影響を考慮して推定する演算を行う。そして、得られた奥行方向の音響物性分布に基づいて、音響物性像の画像データを構築し、所望とする超音波断層像を得るようにしている。 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.
しかしながら、上記装置では、奥行方向に沿って深くなる方向、つまり測定対象物の前面側から背面側に向かう方向に固有音響インピーダンス分布を順次推定していく手法を採用している。このため、深度が大きくなるほど固有音響インピーダンスの推定値の誤差が蓄積、拡大するという欠点があり、より精度の高い画像を構築するうえでのマイナス要因となっていた。また、このような場合、超音波断層像における奥行方向に多くの筋が入ってしまうため、微細な内部の層構造が感覚的に理解しやすいとは言い難く、いまだ改良の余地があった。 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.
上記の課題を解決するために、本願発明者らは測定対象物内の固有音響インピーダンスの推定値を適宜補正することを思い付き、鋭意検討を行った。その結果、１）測定対象物と同じく基体に接して存在する第１参照物質とは別の位置（具体的には、測定対象物において基体と接していない反対の側（奥行方向の奥側））の位置）に存在する第２参照物質の固有音響インピーダンス値を補正演算にて利用したり、２）そのような位置に第２参照物質が存在していなくても、所定の条件を満たす場合には測定対象物内の所定部位の固有音響インピーダンス値（あるいはその推定値の平均）を仮想参照部位の固有音響インピーダンス値と定義して補正演算にて利用したりすれば、深度が大きくなるほど蓄積、拡大しやすい固有音響インピーダンスの推定値の誤差を減じることができ、奥行方向の固有音響インピーダンス分布のより正確な推定が可能となることを新たに知見した。そして、本願発明者らはこれらの知見に基づいてさらに鋭意研究を進めることにより、下記に列挙する解決手段［１］～［１５］を想到するに至ったのである。 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] 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] 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] 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
［４］前記推定ステップでは、前記測定対象物内において、異なる固有音響インピーダンスを持つ無損失の微小伝送路が奥行方向に連なって伝送路の集合体をなしていると仮定して、手前側の前記微小伝送路の固有音響インピーダンスの推定結果に基づきその奥側に隣接する前記微小伝送路の固有音響インピーダンスを推定する演算を行うことを特徴とする手段１乃至３のいずれか１項に記載の超音波画像構築方法。
[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
［５］前記測定対象物は、培養細胞であり、前記第１参照物質及び前記第２参照物質は、前記培養細胞を培養するための培養液であり、前記基体は、前記培養液を収容する培養容器の一部であることを特徴とする請求項１に記載の超音波画像構築方法。
[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
［６］前記測定対象物は、組織表面に対して略平行な方向に延びる血管を含む生体の軟組織としての皮膚であり、前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値は、前記血管中に存在する血液の固有音響インピーダンス値であることを特徴とする請求項２に記載の超音波画像構築方法。
[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
［７］前記測定対象物が生体の軟組織である場合の前記測定対象物内における特定深度の固有音響インピーダンスの推定値の平均は、前記測定対象物が皮膚である場合の当該皮膚内における特定深度の固有音響インピーダンスの推定値の平均であり、前記特定深度は、前記音響物性像の画像データが構築される領域よりも深い位置に設定されることを特徴とする請求項２に記載の超音波画像構築方法。
[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
［８］基体と、前記基体に測定対象物及び第１参照物質が接して存在し、かつ前記測定対象物における前記基体とは反対の側に第２参照物質が接して存在した状態で超音波を送信し、前記基体を介して超音波を入射させたときの前記測定対象物、前記第１参照物質及び前記第２参照物質からの超音波波形のインパルス応答を受信する超音波振動子と、前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内及び前記第２参照物質内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する第１演算手段と、前記第１演算手段により得た前記第２参照物質の固有音響インピーダンスの推定値を前記第２参照物質の固有音響インピーダンスの実際値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する第２演算手段と、前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築手段とを備えたことを特徴とする超音波画像構築装置。 [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] 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] 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] 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
［１２］プロセッサに、基体に測定対象物及び第１参照物質が接して存在し、かつ前記測定対象物における前記基体とは反対の側に第２参照物質が接して存在した状態で超音波振動子に超音波を送信させ、前記基体を介して超音波を入射させたときの前記測定対象物、前記第１参照物質及び前記第２参照物質からの超音波波形のインパルス応答を前記超音波振動子に受信させる送受信ステップと、前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内及び前記第２参照物質内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する推定ステップと、前記推定ステップにて得た前記第２参照物質の固有音響インピーダンスの推定値を前記第２参照物質の固有音響インピーダンスの実際値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する補正ステップと、 前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築ステップとを実行させるための超音波画像構築プログラム。 [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] 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] 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] 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
上記手段１～１５に記載の発明によると、微細な内部構造を有する非常に薄い測定対象物の超音波断層像を、感覚的にその微細な内部構造が理解しやすい態様にて比較的簡単にかつ高い精度で構築することができる。また、超音波を利用して得た情報に基づいて、培養細胞の内部構造や、皮膚等に代表される生体の軟組織の内部構造を簡便にかつ正確に観察、評価等できる画像を構築することができる。 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.
［第１の実施形態］
以下、本発明の超音波画像構築方法及び装置を具体化した第１の実施形態を図１～図１０に基づき詳細に説明する。
[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.
図１は、本実施形態の超音波画像構築装置１を示す概略構成図である。図１に示されるように、本実施形態の超音波画像構築装置１は、超音波を用いて皮膚８を観察、診断等するための装置であって、パルス励起型超音波顕微鏡２と、パーソナルコンピュータ（ＰＣ）３とを備えている。
FIG. 1 is a schematic configuration diagram showing an ultrasonic
パルス励起型超音波顕微鏡２は、ステージ４を有する顕微鏡本体５と、ステージ４の下方に設置された超音波プローブ６とを備える。パルス励起型超音波顕微鏡２の超音波プローブ６は、ＰＣ３と電気的に接続されている。
A pulse excitation
本実施形態のステージ４は、ユーザの手動操作により、水平方向（即ちＸ方向及びＹ方向）に移動できるように構成されている。このステージ４には、測定対象物を接触させて配置するための樹脂プレート９が固定されている。ここでの測定対象物は、組織表面に対して略平行な方向に延びる血管を含む生体の軟組織（具体的には、肌組織；皮膚８）である。本実施形態では、ヒトの皮膚８を樹脂プレート９に直接押し付けることにより測定等を行っている。また、既知の音響物性を有する基体としての樹脂プレート９は、超音波を透過させることができる平板状部材であって、測定対象物である皮膚８よりも硬い材料からなる。このような形状及び硬さの部材を基体として用いた場合、測定対象物である皮膚８を確実に密着配置することが可能となり、奥行方向の固有音響インピーダンス分布を正確に推定可能となる結果、ひいては画像構築の精度が向上する。なお、本実施形態では厚さ１．４ｍｍのポリスチレン板が用いられている。勿論、ポリスチレン以外の樹脂からなる板材などを用いることも許容される。
The
この樹脂プレート９において皮膚８が接触配置される側である上面には、第１参照物質としてのリファレンス部材１０があらかじめ設置されている。リファレンス部材１０は、樹脂プレート９とは異なる既知の音響物性を有している。本実施形態では例えばアクリル樹脂（アクリル接着剤）を付着させることによりリファレンス部材１０としているが、勿論これに限定されるわけではない。リファレンス部材１０に対して密着させることが可能なものであれば、樹脂材以外のもの（例えばガラス材、金属材、セラミック材など）をリファレンス部材１０としてもよい。あるいは、このようなリファレンス部材１０を設置する代わりに、例えば樹脂プレート９の上面に接するように水等を存在させておき、これを第１参照物質として用いてもよい。なお、樹脂プレート９にリファレンス部材１０をあらかじめ設置しておくことにより、装置が置かれる環境の変化等に依存せず、リファレンス部材１０に入射した超音波波形のインパルス応答情報を正確にかつ安定的に取得することができる。
A
超音波プローブ６は、水などの超音波伝達媒体Ｗを貯留可能な貯留部１１をその先端部に有するプローブ本体１２と、プローブ本体１２の略中心部に配置される超音波トランスデューサ１３（超音波振動子）と、プローブ本体１２を前記ステージ４の面方向に沿って二次元的に走査するためのＸ－Ｙステージ１４とを備える。プローブ本体１２の貯留部１１は上部が開口しており、その貯留部１１の開口側を上向きにした状態で超音波プローブ６がステージ４の下方に設置されている。
The
超音波トランスデューサ１３は、例えば酸化亜鉛の薄膜圧電素子１６とサファイアロッドの音響レンズ１７とによって構成される。この超音波トランスデューサ１３は、パルス励起されることで樹脂プレート９の下面側から皮膚８及びリファレンス部材１０に対して超音波を照射する。超音波トランスデューサ１３が照射する超音波は、貯留部１１の超音波伝達媒体Ｗを介して円錐状に収束されて樹脂プレート９の上面（皮膚８の表面付近）で焦点を結ぶようになっている。なお本実施形態では、超音波トランスデューサ１３として、口径１．２ｍｍ、焦点距離１．５ｍｍ、中心周波数８０ＭＨｚ、帯域幅５０～１０５ＭＨｚ（－６ｄＢ）の仕様のものを用いている。
The
図２は、本実施形態の超音波画像構築装置１の電気的な構成を示すブロック図である。
FIG. 2 is a block diagram showing the electrical configuration of the ultrasonic
図２に示されるように、超音波プローブ６は、超音波トランスデューサ１３、Ｘ－Ｙステージ１４、パルス発生回路２１、受信回路２２、送受波分離回路２３、検波回路２４、Ａ／Ｄ変換回路２５、エンコーダ２６、コントローラ２７を備える。
As shown in FIG. 2, the
走査手段としてのＸ－Ｙステージ１４は、超音波の照射点を二次元的に走査させるためのＸステージ１４Ｘ及びＹステージ１４Ｙを備えるとともに、それぞれのステージ１４Ｘ，１４Ｙを駆動するモータ２８Ｘ，２８Ｙを備えている。これらのモータ２８Ｘ，２８Ｙとしては、ステッピングモータやリニアモータが使用される。
The
各モータ２８Ｘ，２８Ｙにはコントローラ２７が接続されており、該コントローラ２７の駆動信号に応答してモータ２８Ｘ，２８Ｙが駆動される。これらモータ２８Ｘ，２８Ｙの駆動により、Ｘステージ１４Ｘを連続走査（連続送り）するとともに、Ｙステージ１４Ｙを間欠送りとなるよう制御することで、Ｘ－Ｙステージ１４の高速走査が可能となっている。
A
また、本実施形態においては、Ｘステージ１４Ｘに対応してエンコーダ２６が設けられ、エンコーダ２６によりＸステージ１４Ｘの走査位置が検出される。具体的には、走査範囲を３００×３００個の測定点（ピクセル）に分割した場合、１回のＸ方向（水平方向）の走査が３００分割される。そして、各測定点の位置がエンコーダ２６によって検出されＰＣ３に取り込まれる。ＰＣ３はそのエンコーダ２６の出力に同期して駆動制御信号を生成し、その駆動制御信号をコントローラ２７に供給する。コントローラ２７は、この駆動制御信号に基づいてモータ２８Ｘを駆動する。また、コントローラ２７は、エンコーダ２６の出力信号に基づきＸ方向の１ラインの走査が終了した時点でモータ２８Ｙを駆動して、Ｙステージ１４ＹをＹ方向に１ピクセル分移動させる。
Further, in this embodiment, an
さらに、コントローラ２７は、駆動制御信号に同期してトリガ信号を生成してパルス発生回路２１に供給する。これにより、パルス発生回路２１において、そのトリガ信号に同期したタイミングで励起パルスが生成される。その励起パルスが送受波分離回路２３を介して超音波トランスデューサ１３に供給される結果、超音波トランスデューサ１３から超音波が照射される。
Furthermore, the
図３は、Ｘ－Ｙステージ１４の移動に伴う超音波の走査範囲Ｒ１の一例を示している。この例では、皮膚８を接触配置させる領域を包囲するようにリファレンス部材１０が設けられている。そして、ヒトの皮膚８を当該領域に押し付けた状態で、リファレンス部材１０がある位置から走査が開始される。そして、矢印で示すように、皮膚８の表面に沿ってＸ方向及びＹ方向に二次元的に走査が順次行われる。
FIG. 3 shows an example of an ultrasonic scanning range R1 accompanying movement of the
超音波トランスデューサ１３の薄膜圧電素子１６は、送受波兼用の超音波振動子であり、皮膚８で反射した超音波（反射波）を電気信号に変換する。そして、その反射波の信号は、送受波分離回路２３を介して受信回路２２に供給される。受信回路２２は、信号増幅回路を含んで構成されていて、反射波の信号を増幅して検波回路２４に出力する。
The thin-
検波回路２４は、皮膚８からの反射波信号を検出するための回路であり、図示しないゲート回路を含む。本実施形態の検波回路２４は、超音波トランスデューサ１３で受信した反射波信号のなかから、皮膚８からの反射波信号やリファレンス部材１０からの反射波信12月号を抽出する。そして、検波回路２４で抽出された反射波信号は、Ａ／Ｄ変換回路２５に供給されてＡ／Ｄ変換された後、ＰＣ３に転送される。
The
ＰＣ３は、ＣＰＵ３１（中央処理装置）、Ｉ／Ｆ回路３２、メモリ３３、記憶装置３４、入力装置３５、及び表示装置３６を備え、それらはバス３７を介して相互に接続されている。
The
ＣＰＵ３１は、メモリ３３を利用して制御プログラムを実行し、システム全体を統括的に制御する。制御プログラムとしては、Ｘ－Ｙステージ１４による二次元走査を制御するためのプログラム、超音波Ｂモードエコー像の元となる反射信号列のデータを固有音響インピーダンス像へ変換するためのプログラム、固有音響インピーダンス像を表示するためのプログラムなどを含む。なお、ＣＰＵ３１とは別に例えばＤＳＰ（Ｄｉｇｉｔａｌ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｏｒ：デジタル信号プロセッサ）を設けて、そこでＣＰＵ３１が行っている信号処理の一部を行わせてもよい。
The
Ｉ／Ｆ回路３２は、超音波プローブ６との間で信号の授受を行うためのインターフェース（具体的には、ＵＳＢインターフェース）である。Ｉ／Ｆ回路３２は、超音波プローブ６に制御信号（コントローラ２７への駆動制御信号）を出力したり、超音波プローブ６からの転送データ（Ａ／Ｄ変換回路２５から転送されるデータなど）を入力したりする役割を果たすものである。なお、超音波プローブ６との間で信号の授受を行う場合には、上記のような物理的なインターフェースに限定されることはなく、無線インターフェースを用いてもよい。
The I/
表示装置３６は、例えば、液晶、プラズマ、有機ＥＬ（ｅｌｅｃｔｒｏｌｕｍｉｎｅｓｃｅｎｃｅ）等のモニタディスプレイである。表示装置３６は、カラー表示、モノクロ表示を問わずに使用できるが、カラー表示であることが望ましい。この表示装置３６は、皮膚８の表層の固有音響インピーダンス像を表示したり、各種設定の入力画面を表示したりするために用いられる。
The
入力装置３５は、タッチパネル、マウス、キーボード、ポインティングデバイス等の入力ユーザインタフェースであって、ユーザからの要求や指示、パラメータの入力に用いられる。
The
記憶装置３４は、磁気ディスク装置や光ディスク装置などのハードディスクドライブであり、各種の制御プログラム及び各種のデータを記憶している。メモリ３３は、ＲＡＭ（ランダムアクセスメモリ）やＲＯＭ（リードオンリーメモリ）を含み、超音波測定のためにあらかじめ取得されたリファレンス部材１０の反射波形とその固有音響インピーダンスとを保存する。ＣＰＵ３１は、入力装置３５による指示に従い、プログラムやデータを記憶装置３４からメモリ３３へ転送し、それを逐次実行する。なお、ＣＰＵ３１が実行するプログラムとしては、メモリカード、フレキシブルディスク、光ディスクなどの記憶媒体に記憶されたプログラムや、通信媒体を介してダウンロードしたプログラムでもよく、その実行時には記憶装置３４にインストールして利用する。
The
次に、本実施形態の超音波画像構築装置１において、超音波Ｂモードエコー像の元となる反射信号列から固有音響インピーダンス像を構築する手法について説明する。
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
この超音波画像構築装置１では、第１参照物質であるリファレンス部材１０に入射した超音波波形のインパルス応答情報、及び測定対象物である皮膚８に入射した超音波波形のインパルス応答情報から、規格化されたインパルス応答情報を得て、その規格化されたインパルス応答情報に基づいて奥行方向の音響物性分布を多重反射の影響を考慮して推定するようになっている。また、このような推定を行うために、本実施形態では、測定対象物内において、異なる固有音響インピーダンスを持つ無損失の微小伝送路５１が奥行方向に連なって伝送路の集合体をなしていると仮定して、手前側の微小伝送路５１の固有音響インピーダンスの推定結果に基づきその奥側に隣接する微小伝送路５１の固有音響インピーダンスを推定する処理を順次繰り返すことにより、伝送路の奥行方向の音響物性分布（ここでは固有音響インピーダンス分布）を推定する演算を行うようになっている。このような演算は、メモリ３３内に格納された所定のアルゴリズムに基づいて実行される。
In this ultrasonic
このアルゴリズムは、超音波Ｂモードエコー像の元となる反射信号列を用いることによって、奥行方向の固有音響インピーダンスの分布の推定を行うアルゴリズムである。このアルゴリズムは、時間領域反射測定法（ＴＤＲ法：Ｔｉｍｅ Ｄｏｍａｉｎ Ｒｅｆｌｅｃｔｏｍｅｔｒｙ法）の原理を参考とするものであって、皮膚組織内部での多重反射を考慮した時間－周波数領域における解析を通じて、超音波Ｂモードエコー像の元となる反射信号列を奥行方向の固有音響インピーダンス像に変換するアルゴリズムである。以下、これについて具体的に説明する。 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.
図４（ａ）は、実際に測定を行ったときにおける測定対象物からの反射波形の取得について説明する図であり、図４（ｂ）は、測定対象物を微小伝送路５１に見立てたときにおける反射波形の取得について説明する図である。図５（ａ）は、実際に測定を行ったときにおける第１参照物質からの反射波形の取得について説明する図であり、図５（ｂ）は、第１参照物質を微小伝送路５１に見立てたときにおける反射波形の取得について説明する図である。 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
まず、図４（ａ）に示すように、超音波トランスデューサ１３を作動させ、基体としての樹脂プレート９を介して、測定対象物にとって十分な焦点深度を持った超音波の収束ビームを送信する。そして、測定対象物である皮膚８に超音波の収束ビームを入射させ、そこからの反射波形を取得する。その時のインパルス応答Γ０（ω）は、フーリエ変換を用いて、入射波Ｓ０と皮膚８からの反射波Ｓｔｇｔ（ω）とから次式１のように表される。
この場合、固有音響インピーダンスが既知かつ均一であって、測定対象物に比べて十分な厚さを有するリファレンス部材１０からの反射波形も取得する必要がある。リファレンス部材１０からの反射波Ｓｒｅｆ（ω）は、リファレンス部材１０の固有音響インピーダンスＺｒｅｆと樹脂プレート９の固有音響インピーダンスＺ０とを用いて、次式２のように表される。
また、測定対象物である皮膚８からのインパルス応答Γ０（ω）は次式３で表される。ただし、このインパルス応答Γ０（ω）には皮膚８の組織の奥の複数の界面から発生する反射が含まれているため、インパルス応答Γ０（ω）は周波数特性を持ったものとなる。なお、ここまでの式により、第１参照物質に入射した超音波波形のインパルス応答情報及び測定対象物に入射した超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報が求められる。
ここで、図６は、樹脂プレート９に接する微小伝送路５１から順に、各微小伝送路５１の特性インピーダンスＺ１、Ｚ２…Ｚｎを推定していく様子を概念的に示した図である。この図に示されるように、各微小伝送路５１の特性インピーダンスＺ１、Ｚ２…Ｚｎを、樹脂プレート９に接している微小伝送路５１から奥行方向に向かって順に推定していく。 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 .
図７は、多重反射の影響について概念的に示した図である。次式４は、インパルス応答Γ０（ω）を逆フーリエ変換したｇ０（ｔ）を表したものであるが、その第１項は多重反射の影響を受けない（図７参照）。従って、この第１項の値から樹脂プレート９に接する微小伝送路５１の特性インピーダンスＺ１を、次式５のように推定することができる。
周波数領域の皮膚８の固有音響インピーダンスＺｘ０は、Γ０を用いると次式６のように表される。
また、Ｚｘ０はさらに奥からのインパルス応答Γ１を用いて次式７のようにも表される。
ここで、次式８、９にて表すように、γは伝播定数、αは減衰定数、βは位相定数、ｆは周波数であるが、本実施形態のアルゴリズムではα＝０、及び皮膚８の全微小伝送路５１の音速をｃ＝１６００（ｍ／ｓ）と仮定している。
また、各微小伝送路５１の距離Δｌは次式１０のように表され、ここでも音速をｃ＝１６００（ｍ／ｓ）と仮定している。Δｔは皮膚８からの反射波形のサンプリング間隔の１ポイントに相当する（本実施形態ではΔｔ＝２（ｎｓ））。
そして、上記の式をもとに、さらに奥にある微小伝送路５１からのインパルス応答Γ１を求めることができる（次式１１）。即ち、Ｚｘ０及びＺ１の値をもとに、Ｚ１の終点におけるΓ１の値を推定することができる。
次式１２は、インパルス応答Γ１（ω）を逆フーリエ変換したｇ１（ｔ）を表したものであって、その第１項は多重反射の影響を受けない。従って、この第１項の値から当該微小伝送路５１に隣接するさらに奥側の微小伝送路５１の特性インピーダンスＺ２を推定することができ、同様にＺｘ１、Γ２（ω）も推定することができる（次式１３，１４，１５）。
この工程を繰り返すことによって、各微小伝送路５１の特性インピーダンス（固有音響インピーダンス）Ｚ１、Ｚ２…Ｚｎを推定することができる。
By repeating this process, the characteristic impedances (specific acoustic impedances) Z 1 , Z 2 . . . Zn of each
本実施形態のアルゴリズムでは、上記のような推定ステップを行った後、さらに以下に示す所定の補正ステップを実施する。即ち、この補正ステップでは、まず仮想参照部位の固有音響インピーダンス値を定義する。ここでの測定対象物は皮膚８（肌組織）であって（図８参照）、皮膚８を構成する表皮６１及び真皮６２よりも深い位置にある皮下組織６３の比較的浅い部位には、肌組織表面に対して略平行な方向に延びる血管６４が存在している。血管６４内に存在する血液６５の固有音響インピーダンス値は既知（１．６１ＭＮｓ／ｍ３）であるため、本実施形態ではこの値を上記「仮想参照部位の固有音響インピーダンス値」として用いている。 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.
この補正演算を図９のグラフを用いて詳細に説明する。同グラフにおいて縦軸は固有音響インピーダンスを示し、横軸は時間を示している。また、同グラフには、２つの曲線Ｃ１、Ｃ２が描かれている。実線で描かれた曲線Ｃ１は補正前の固有音響インピーダンス曲線Ｃ１を示すものであり、破線で描かれた曲線Ｃ２は補正後の固有音響インピーダンス曲線Ｃ２を示すものである。なお、これら曲線Ｃ１、Ｃ２は、左半分の領域（即ち樹脂プレート９内部に相当する領域）において差がなく重なっているが、右半分の領域（即ち皮膚８内部に相当する領域）において相違している。本実施形態では、曲線Ｃ１の右側領域よりも曲線Ｃ２の右側領域のほうが全体的に下方に位置している。 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.
ここで、基体である樹脂プレート９はポリスチレン製であるため、その固有音響インピーダンス値Ｚ（ｘ）は既知かつ一定（Ｚ（０）＝２．５０ＭＮｓ／ｍ３）である。固有音響インピーダンス曲線Ｃ１における左側の領域が、超音波トランスデューサ１３に近い位置にある樹脂プレート９の固有音響インピーダンス値Ｚ（ｘ）を示す部分である。一方、樹脂プレート９に接して配置された皮膚８は、ポリスチレン製の樹脂プレート９よりも軟らかい。このため、その固有音響インピーダンス値Ｚ（ｘ）は、固有音響インピーダンス曲線Ｃ１における点Ａを境にして急激に低くなる。固有音響インピーダンス曲線Ｃ１における右側の領域が皮膚８の固有音響インピーダンス値Ｚ（ｘ）の補正前の推定値を示す部分であり、２．５０ＭＮｓ／ｍ３よりもかなり低い値となっている。当該曲線Ｃ１において点Ｂ（即ち皮膚８内における特定深度に相当する部分）よりも右側はフラットであり、この部分が特定深度の血管内に存在する血液の固有音響インピーダンス値の補正前の推定値を示す部分となっている。ただし、当該推定値は誤差を含んでおり、実際の値とはいくぶん異なるものとして表されている（この例では１．８０ＭＮｓ／ｍ３）。そこで、特定深度の血管６４内に存在する血液６５の固有音響インピーダンス値の補正前の推定値を、仮想参照部位の固有音響インピーダンス値（１．６１ＭＮｓ／ｍ３）に置き換えるべく、１．８０ＭＮｓ／ｍ３から補正値である０．１９ＭＮｓ／ｍ３を減算する。その結果、当該曲線Ｃ１のフラットな部分が下方に移動し、これが補正された固有音響インピーダンス曲線Ｃ２のうちの一部をなすものとなる。なお、当該曲線Ｃ１における点Ｂは、曲線Ｃ２における点Ｃに移動した状態となる。
Here, since the
次に、固有音響インピーダンス曲線Ｃ１における点Ａと点Ｂとの間の領域（即ち皮膚表面から当該血管までの部分に相当する領域）を補正する。その際、先の補正値である０．１９ＭＮｓ／ｍ３を最大とし、測定対象物の前面側（即ち皮膚表面）に向かうほど小さくなる補正値を用いる。そしてこの補正値を減算することにより、当該血管から皮膚表面までの固有音響インピーダンスの推定値を補正する。換言すると、皮膚表面からの深さに応じた補正値を設定し、推定値からこれを減算する。例えば、特定深度の半分の深さ位置における推定値を補正する場合には、０．１９ＭＮｓ／ｍ３の１／２の値（即ち０．０９５ＭＮｓ／ｍ３）をその位置における補正値として設定し、これをその位置における推定値から減算する。また、特定深度の１／４の深さ位置における推定値を補正する場合には、０．１９ＭＮｓ／ｍ３の１／４の値（即ち０．０４７５ＭＮｓ／ｍ３）をその位置における補正値として設定し、これをその位置における推定値から減算する。ちなみに、図９のグラフに描かれた３本の下向き矢印は、それぞれの深度における補正値の大きさの違いを視覚的に示している。 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.
そして、このような補正ステップの演算を通じて、測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定し、最終的にはＢモードエコー像の元となる反射信号列を固有音響インピーダンス像に変換する。 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.
次に、本実施形態の超音波画像構築装置１において固有音響インピーダンス画像を構築するために、プロセッサであるＣＰＵ３１が実行する演算処理について、図１０のフローチャートを用いて説明する。
Arithmetic processing executed by the
まず、測定対象物であるヒトの皮膚８（例えば、比較的浅い位置に太い血管（けい静脈、けい動脈）が存在している首の皮膚８など）を樹脂プレート９の上面に押し付けるようにして接触配置させる。この状態で、まず超音波プローブ６に初期動作を行わせる。即ち、ＣＰＵ３１からの指示に基づいてコントローラ２７を作動させることにより、モータ２８Ｘ，２８Ｙを駆動し、リファレンス部材１０がある位置にて超音波照射が行われるようにＸ－Ｙステージ１４を移動させる。
First, the human skin 8 (for example,
またこのとき、ＣＰＵ３１からの指示に基づいて励起パルスがトランスデューサ１３に供給されると、図５（ａ）に示すように、リファレンス部材１０に超音波Ｓｏが照射され、その反射波Ｓｒｅｆ（ω）が受信回路２２を経て検波回路２４で検出される。そして、反射波取得手段としてのＣＰＵ３１は、Ａ／Ｄ変換回路２５で変換されたデジタルデータをＩ／Ｆ回路３２を介して取得し、そのデータをリファレンス部材１０からの超音波波形のインパルス応答のデータとしてメモリ３３に記憶する（ステップＳ１００）。
Also at this time, when an excitation pulse is supplied to the
その後、ＣＰＵ３１からの指示に基づいてコントローラ２７によりモータ２８Ｘ，２８Ｙが駆動され、Ｘ－Ｙステージ１４による二次元走査が開始される。ＣＰＵ３１は、エンコーダ２６の出力に基づいて測定点の座標データを取得する（ステップＳ１１０）。
After that, the
そして、図４（ａ）に示すように、ＣＰＵ３１からの指示に基づいて励起パルスがトランスデューサ１３に供給されることにより、皮膚８に超音波Ｓｏが照射され、その反射波Ｓｔｇｔ（ω）が受信回路２２を経て検波回路２４で検出される。反射波取得手段としてのＣＰＵ３１は、Ａ／Ｄ変換回路２５で変換されたデジタルデータをＩ／Ｆ回路３２を介して取得し、そのデータを皮膚８からの超音波波形のインパルス応答のデータとして座標データに関連付けてメモリ３３に記憶する（ステップＳ１２０）。
Then, as shown in FIG. 4A, an excitation pulse is supplied to the
次いで、第１演算手段としてのＣＰＵ３１は、規格化されたインパルス応答信号のデータを用いて、上記のアルゴリズムのうち、ＴＤＲ法の原理を参考とした推定ステップの演算を実行する。そしてＣＰＵ３１は、その演算により皮膚８における測定点での奥行方向の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定し、その推定結果を座標データに関連付けてメモリ３３に記憶する（ステップＳ１３０）。
Next, the
次いで、第２演算手段としてのＣＰＵ３１は、上記のアルゴリズムのうちの補正ステップの演算を実施する。即ち、仮想参照部位の固有音響インピーダンス値を定義した後、皮膚８内における特定深度の固有音響インピーダンスの推定値を前記仮想参照部位の固有音響インピーダンス値に置き換える演算を行い、この演算を通じて皮膚８内の固有音響インピーダンスの推定値を補正する。そしてＣＰＵ３１は奥行方向の固有音響インピーダンス分布を推定し、その推定結果を座標データに関連付けてメモリ３３に記憶する（ステップＳ１３２）。
Next, the
その後、画像構築手段としてのＣＰＵ３１は、奥行方向の固有音響インピーダンス分布の推定結果に基づいて、固有音響インピーダンス像（断層像）を構築するための画像処理を行う（ステップＳ１４０）。詳しくは、ＣＰＵ３１は、固有音響インピーダンス分布の推定結果に基づいてカラー変調処理を行い、固有音響インピーダンスの大きさに応じて色分けして表示した画像データを構築し、該画像データをメモリ３３に記憶する。
After that, the
次いで、ＣＰＵ３１は、全ての測定点での処理が終了して、全ての測定点で画像データが取得されたか否かを判断する（ステップＳ１５０）。ここで、全データが取得されていない場合には（ステップＳ１５０：ＮＯ）、ＣＰＵ３１は、ステップＳ１１０に戻って、ステップＳ１１０～Ｓ１４０の処理を繰り返して実行する。全データが取得された場合には（ステップＳ１５０：ＹＥＳ）、ＣＰＵ３１は、次ステップＳ１６０に移行する。
Next, the
そして、ＣＰＵ３１は、該データを表示装置３６に転送し、あらかじめ定めた直線上における固有音響インピーダンス像（断層像）を表示させた後（ステップＳ１６０）、図１０の処理を終了するようになっている。このような一連の処理により、皮膚８での固有音響インピーダンスの大きさに応じて色分けされた固有音響インピーダンス像（断層像）が表示される。
Then, the
従って、本実施の形態によれば以下の効果を得ることができる。 Therefore, according to this embodiment, the following effects can be obtained.
（１）本実施形態の超音波画像構築装置１では、リファレンス部材１０に入射した超音波波形のインパルス応答情報、及び皮膚８に入射した超音波波形のインパルス応答情報から、規格化されたインパルス応答情報、つまり装置に依存しないインパルス応答情報を得るようにしている。そして、この規格化されたインパルス応答情報に基づいて、皮膚８内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって、多重反射の影響を考慮しつつ順次推定する演算（推定ステップ）を行うようにしている。本実施形態ではこれに加えて、皮膚８内における特定深度（ここでは、横行する血管が存在する皮下組織浅部）の固有音響インピーダンスの推定値を、仮想参照部位の固有音響インピーダンス値に置き換える演算を通じて、皮膚８内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する演算（補正ステップ）も行うようにしている。以上の結果、深度が大きくなるほど蓄積、拡大しやすい固有音響インピーダンスの推定値の誤差を減じることができ、奥行方向の固有音響インピーダンス分布のより正確な推定が可能となる。また、補正ステップを行わない従来技術では超音波断層像における奥行方向に多くの筋が入ってしまうため、微細な層構造が感覚的に理解しにくかったのに対し、本実施形態ではこのような筋が少なくなり層方向の色ムラが減少する。従って、微細な層構造を有する非常に薄い皮膚８の超音波断層像を、その層構造が感覚的に理解しやすい固有音響インピーダンス像として比較的簡単にかつ高い精度で構築することができる。なお、この装置１により得られる固有音響インピーダンス像は、測定対象物を切断することなく（即ち非侵襲で）、層ごとの力学特性の断面分布（深さ分布）情報を、推定した固有音響インピーダンスの絶対値ごとに色分けして画像化したものである。それゆえ、この像は感覚的に層構造が理解しやすいものなっている。
(1) In the ultrasonic
ここで、一般に通常の超音波診断装置で得られる超音波Ｂモードエコー像では、皮膚８等の生体組織内部の層情報は一応得られるものの、得られる画像は固有音響インピーダンスがある程度以上の差を有する層同士の界面からの反射像である。即ち、固有音響インピーダンスの差がある程度小さくなると、組織学的には界面が存在しても検出されず、その構造を反映した画像を形成することが極めて困難となっていた。つまり、一般的な反射像では生体組織の微細な内部構造、微細な層構造が反映された生体組織内部の反射像（固有音響インピーダンスの違い）を把握するのには不十分であった。これに対して、この超音波画像構築装置１によると、従来の超音波Ｂモードによって全く検出することができなかった力学特性分布に基づく皮膚８の層構造を、十分な解像度を持った鮮明な断層像として捉えることができるようになった。また、このような鮮明な断層像は、他の非侵襲可視化装置（光干渉断層撮影装置（ＯＣＴ）や、ｉｎ ｖｉｖｏ共焦点顕微鏡など）では取得不能であったため、この超音波画像構築装置１が具現化されたことの意義は大きい。以上のように、本実施形態の超音波画像構築装置１によると、皮膚８の状態（皮膚８の層ごとの力学特性に関する状態）を簡便にかつ非侵襲的に評価することができる。
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
（２）本実施形態の超音波画像構築装置１では、皮膚組織内部での多重反射を考慮した時間－周波数領域における解析を通じて、超音波Ｂモードエコー像の元となる反射信号列を奥行方向の固有音響インピーダンス像に変換する手法を採用している。このため、比較的簡単に所定の演算を行うことができる。また本実施形態では、上述したアルゴリズムを用いて、手前側の微小伝送路５１の固有音響インピーダンスの推定結果に基づきその奥側に隣接する微小伝送路５１の固有音響インピーダンスを推定する処理を順次繰り返すことにより、伝送路の奥行方向の固有音響インピーダンス分布を推定する演算を行うようにしている。従って、この演算によれば、多重反射の影響を最小限にすることができ、奥行方向の固有音響インピーダンス分布を精度よく推定することができる。
(2) In the ultrasonic
［第２の実施形態］
以下、本発明の超音波画像構築方法及び装置を具体化した第２の実施形態を図１１に基づき詳細に説明する。なお、上記第１の実施形態と共通している構成については、同じ部材番号を付す代わりに詳細な説明を省略する。
[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.
本実施形態の超音波画像構築装置１では、所定の推定ステップ後に補正ステップを行ってから画像構築ステップを行っており、この点については上記第１の実施形態と共通する（上記図９、図１０等を参照）。ただし、第１の実施形態の補正ステップでは、肌組織表面に対して略平行な方向に延びる血管６４内に存在する血液６５の固有音響インピーダンス値（１．６１ＭＮｓ／ｍ３）を「仮想参照部位の固有音響インピーダンス値」として用いたのに対し、ここでは測定対象物が生体の軟組織であることに着目してその特定深度の固有音響インピーダンスの推定値の平均を「仮想参照部位の固有音響インピーダンス値」と定義して用いている。本実施形態では、肌組織表面から３５０μｍほど深い位置における組織の固有音響インピーダンスの推定値の平均（例えば１．６５ＭＮｓ／ｍ３）を求めて、その値を用いている。これは、生体の軟組織中には骨や軟骨などの硬い部分が存在せず硬さのばらつきが少ないので、その固有音響インピーダンスは例えば１．５０ＭＮｓ／ｍ３～１．７０ＭＮｓ／ｍ３の範囲から大きく外れた値とはならず、仮想参照部位として選択するのに適している、との考え方に立脚したものである。なお、脂肪を除く軟組織であれば、固有音響インピーダンスは、ほぼ１．６０ＭＮｓ／ｍ３～１．７０ＭＮｓ／ｍ３の範囲内となる。
In the ultrasonic
この場合、皮膚８における特定深度は、固有音響インピーダンス像の画像データが構築される領域よりも深い位置に設定されることが好ましい。本実施形態では、例えば肌組織表面から３００μｍほど深い位置までを画像構築領域として設定している関係上、組織の固有音響インピーダンスの推定値の平均を求めるための特定深度はそれよりも５０μｍ深い位置としている。このようにすることで、特定深度の固有音響インピーダンス像を敢えて表示しないようにし、かえって不自然なものとならないようにしている。
In this case, the specific depth in the
ここで、図１１（ａ）は補正ステップを実施することなく画像を構築して得た固有音響インピーダンス像であり、（ｂ）は補正ステップを実施したうえで画像を構築して得た固有音響インピーダンス像である。前者では、奥行方向に多くの筋が入ってしまい、層構造がわかりにくい。これに対して後者では、筋がいくぶん解消された状態で表示されるため、層方向に色が均一化される傾向が強くなり、層構造がわかりやすくなる。 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.
従って、以上説明した本実施形態の超音波画像構築装置１についても、第１の実施形態の超音波画像構築装置１と同様の作用効果を奏することができる。つまり、微細な層構造を有する非常に薄い皮膚８の超音波断層像を、その層構造が感覚的に理解しやすい態様にて比較的簡単にかつ高い精度で構築することができる。また、超音波を利用して得た情報に基づいて、皮膚８の微細な層構造を簡便にかつ正確に観察、評価等できる画像を構築することができる。特に本実施形態によれば、測定部位の内層に第１の実施形態のときのような適当な仮想参照部位（例えばけい静脈、けい動脈などの比較的太い血管）が存在していない場合であっても、固有音響インピーダンスの推定値を補正することができる。よって、首の部位に限らずいろいろな部位（例えば頬など）において観察、評価を行うことができる。
Therefore, the ultrasonic
［第３の実施形態］
以下、本発明の超音波画像構築方法及び装置を具体化した第３の実施形態を図１２～図１４に基づき詳細に説明する。なお、上記第１の実施形態と共通している構成については、同じ部材番号を付す代わりに詳細な説明を省略する。
[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.
図１２は、本実施形態の超音波画像構築装置１０１を示す概略構成図である。図１２に示されるように、この超音波画像構築装置１０１は、超音波を用いて培養細胞８Ａを観察するための装置であって、パルス励起型超音波顕微鏡２と、パーソナルコンピュータ（ＰＣ）３とを備えている。
FIG. 12 is a schematic configuration diagram showing the ultrasonic
ステージ４には、測定対象物である培養細胞８Ａを培養するための培養液Ｍを収容する培養容器１０２が設置されている。ここでは、培養細胞８Ａとして、接着性細胞の一種であるヒトのグリア細胞８Ａを用いている。グリア細胞８Ａの厚さは非常に薄く、数μｍ程度である。この培養容器１０２の底部中央部には、上面側に培養細胞８Ａが培養時に接着して支持される樹脂プレート９が固定されている。既知の音響物性を有する基体としての樹脂プレート９は、超音波を透過させることができる平板状部材であって、測定対象物である培養細胞８Ａよりも硬い材料からなる。
A
なお、培養容器１０２内において培養細胞８Ａは、培養液Ｍに完全に浸漬された状態で存在している。従って、培養液Ｍは、樹脂プレート９において培養細胞８Ａが接触配置される側である上面に接して存在するばかりでなく、培養細胞８Ａにおける樹脂プレート９とは反対の側にも接して存在していると把握できる。培養液Ｍは樹脂プレート９とは異なる既知の音響物性を有しており、本実施形態ではこれを第１参照物質及び第２参照部材として位置付けて利用している。以降の説明の便宜のため、樹脂プレート９と接して存在する培養液（即ち第１参照物質である培養液）にＭ１を付し、培養細胞８Ａにおける樹脂プレート９とは反対の側に接して存在する培養液（即ち第２参照物質である培養液）にＭ２を付して区別する。
It should be noted that the
次に、本実施形態の超音波画像構築装置１０１において、超音波Ｂモードエコー像の元となる反射信号列から固有音響インピーダンス像を構築する手法について説明する。
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
この超音波画像構築装置１０１では、基体である樹脂プレート９に測定対象物である培養細胞８Ａ及び第１参照物質である培養液Ｍ１が接して存在し、かつ培養細胞８Ａにおける樹脂プレート９とは反対の側に第２参照物質である培養液Ｍ２が接して存在した状態で超音波を送信する。そして、樹脂プレート９を介して超音波を入射させたときの培養細胞８Ａ、第１参照物質である培養液Ｍ１及び第２参照物質である培養液Ｍ２からの超音波波形のインパルス応答をそれぞれ受信する（送受信ステップ）。
In this ultrasonic
次いで、培養液Ｍ１、Ｍ２に入射した超音波波形のインパルス応答情報、及び培養細胞８Ａに入射した超音波波形のインパルス応答情報から、規格化されたインパルス応答情報を得て、その規格化されたインパルス応答情報に基づいて培養細胞８Ａ内及び第２参照部材である培養液Ｍ２内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する（推定ステップ）。このステップでは、奥行方向の固有音響インピーダンス分布を多重反射の影響を考慮して推定する演算を行うが、これについては第１の実施形態において詳しく説明したため割愛する。なお、この演算は、メモリ３３内に格納された所定のアルゴリズムに基づいて実行される。
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
本実施形態の超音波画像構築装置１０１では、上記の推定ステップ後に補正ステップを行ってから画像構築ステップを行っており、この点については上記第１の実施形態と共通する（上記図９、図１０等を参照）。ただし、本実施形態の補正ステップでは第１の実施形態のときとは異なる以下の補正ステップを実施する。即ちここでは、推定ステップにて得た第２参照物質（培養液Ｍ２）の固有音響インピーダンスについて、その「推定値」を「実際値」に置き換える補正演算を行う。培養液Ｍ２の前記「実際値」が例えば１．５５ＭＮｓ／ｍ３であり、誤差を含んでいる可能性が高い前記「推定値」が例えば１．８０ＭＮｓ／ｍ３であるとすると、その値から補正値である０．２５ＭＮｓ／ｍ３を減算することにより、前記「実際値」を補正する。その結果、先の実施形態で示した図９のグラフに示す補正前の固有音響インピーダンス曲線Ｃ１のフラットな部分が下方に移動し、これが補正された固有音響インピーダンス曲線Ｃ２のうちの一部をなすものとなる。なお、当該曲線Ｃ１における点Ｂは、曲線Ｃ２における点Ｃに移動した状態となる。
In the ultrasonic
次に、固有音響インピーダンス曲線Ｃ１における点Ａと点Ｂとの間の領域を補正する。その際、先の補正値である０．２５ＭＮｓ／ｍ３を最大とし、測定対象物の前面側（即ち培養細胞８Ａにおいて樹脂プレート９と接している面側）に向かうほど小さくなる補正値を用いる。そして前記「推定値」からこの補正値を減算することにより、培養液Ｍ２と培養細胞８Ａとの界面から、培養細胞８Ａと樹脂プレート９との界面までの前記「推定値」を補正する。
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
そして、このような補正ステップの演算を通じて、測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定し、最終的にはＢモードエコー像の元となる反射信号列を固有音響インピーダンス像に変換する。 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.
従って、以上説明した本実施形態の超音波画像構築装置１０１についても、第１の実施形態の超音波画像構築装置１と同様の作用効果を奏することができる。つまり、微細な内部構造を有する非常に薄い培養細胞８Ａの超音波断層像を、その微細な内部構造が感覚的に理解しやすい態様にて比較的簡単にかつ高い精度で構築することができる。また、超音波を利用して得た情報に基づいて、培養細胞８Ａの内部構造を簡便にかつ正確に観察、評価等できる画像を構築することができる。
Therefore, the ultrasonic
図１３はグリア細胞８Ａを概略的に示した図である。グリア細胞８Ａには、核１１２が存在するほか、細胞骨格１１３（アクチン繊維）が重合して形成される仮足１１４等が存在している。なお、仮足１１４は移動するグリア細胞８Ａの進行方向側に見られる。図１４はこの超音波画像構築装置１０１を用いて観察されたグリア細胞８Ａの超音波断層像である。これによると、グリア細胞８Ａの内部には固有音響インピーダンスの高い箇所と低い箇所とが存在していることがわかった。図１３の概略図を参照すると、固有音響インピーダンスの高い箇所、即ち比較的硬い箇所は、細胞骨格１１３を多く含む仮足１１４であると推察された。また、固有音響インピーダンスの低い箇所、即ち比較的柔らかい箇所は、核であると推察された。このような本実施形態の超音波画像構築装置１０１によれば、培養細胞８Ａの内部構造が観察できるばかりでなく、培養細胞８Ａの動き等についても観察できるものと考えられた。
FIG. 13 is a diagram schematically showing
なお、本発明の各実施の形態は以下のように変更してもよい。 Each embodiment of the present invention may be modified as follows.
・上記第１、第２実施形態の超音波画像構築装置１では、リファレンス部材１０からの反射波を参照波形として用いて演算処理を行ったが、これに限定されるものではない。この参照波形としては、樹脂プレート９の上面において皮膚８が接触していない箇所からの反射波であればよく、例えば、樹脂プレート９の上面において皮膚８もリファレンス部材１０も接触していない箇所（具体的には、リファレンス部材１０よりも外側に位置する樹脂プレート９の表面）での反射波を用いてもよい。言い換えると、樹脂プレート９と空気層との界面からの反射波を参照波形として利用してもよい。
In the ultrasonic
・上記第１、第２実施形態の超音波画像構築装置１では、下方から超音波を照射する倒立型の超音波顕微鏡２を用いて超音波の照射を行ったが、上方から超音波を照射する正立型の超音波顕微鏡を用いてもよい。
In the ultrasonic
・上記第１、第２実施形態では、基本的に疾患を有していない比較的健康な皮膚８を対象として、その状態を評価する目的で超音波画像構築装置１を用いたが、これに限定されない。例えば、皮膚がんなどの疾患に伴う皮膚の異常を早期に検出する目的で超音波画像構築装置１を用いてもよい。
- In the above-described first and second embodiments, the ultrasonic
・上記第１、第２実施形態の超音波画像構築装置１では、測定対象物がヒトの首の皮膚８であったが、首以外の他の部位（例えば頬など）における皮膚８であっても勿論よい。また、測定対象物は皮膚８でなくてもよく、例えば内臓、筋肉、脳、歯、爪、骨の表層部などであっても勿論よい。また、上記第３実施形態の超音波画像構築装置１０１では、測定対象物がヒトのグリア細胞であったが、付着性細胞であれば生物種は問わずどのような培養細胞であってもよい。さらにいうと、測定対象物は必ずしも生体組織や生物でなくてもよく、非生物（例えば塗膜など）であってもよい。換言すると、本発明の超音波画像構築装置１、１０１は医療分野、美容分野、化粧品分野のみに限定されず、例えば工業分野などの分野においても使用されることができる。
In the ultrasonic
・上記実施形態の超音波画像構築装置１、１０１は、いずれも測定対象物に対して超音波トランスデューサ１３を二次元方向に相対的に走査させる走査手段を備えていたが、これに代えて超音波トランスデューサ１３を一次元方向にのみ相対的に走査させる走査手段を備えたものとしてもよい。また、走査手段は必須の構成ではないため省略しても勿論よく、この場合には装置を小型化、簡略化、低コスト化することが可能となる。
The ultrasonic
・上記実施形態の超音波画像構築装置１、１０１では、奥行方向の固有音響インピーダンス分布の推定結果に基づいて固有音響インピーダンス像を構築したが、これに限定されない。例えば、奥行方向の音速分布を推定し、その結果に基づいて音速像を構築してもよい。
In the ultrasonic
・上記実施形態の超音波画像構築装置１、１０１では、超音波Ｂモードエコー画像の元となる反射信号列から固有音響インピーダンス像を構築してそれを表示装置３６に表示させるように構成したが、固有音響インピーダンス像ばかりでなく超音波Ｂモードエコー像も表示できるようにしても勿論よい。また、超音波Ｂモードエコー像を表示する汎用の超音波診断装置に上記実施形態のアルゴリズムを組み込むことで、超音波画像構築装置１、１０１として動作させるようにしてもよい。
In the ultrasonic
・上記実施形態では、測定対象物内での多重反射を考慮した時間－周波数領域における解析を通じて、超音波Ｂモードエコー像の元となる反射信号列を奥行方向の固有音響インピーダンス像に変換する手法を採用したが、これに限定されない。反射波形から固有音響インピーダンスの分布を推定するための方法としては、例えば、伝送路内部での多重反射を含む全ての反射経路を想定して時間軸上で随時応答を解析していく手法を採用してもよい。この手法であっても、上記実施形態の手法のときと同様の固有音響インピーダンス分布の推定結果を得ることができる。 ・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.
・上記実施形態では、超音波Ｂモードエコー像の元となる反射信号列を奥行方向の固有音響インピーダンス像に変換するにあたり、測定対象物内での多重反射を考慮した解析手法を採用したが、これに限定されない。例えば、生体の軟組織のように多重反射の影響が小さいと考えられる場合、測定対象物内での多重反射を敢えて考慮しない解析手法を採用してもよい。 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.
１、１０１…超音波画像構築装置
８…測定対象物としての生体の軟組織（皮膚）
８Ａ…測定対象物としての培養細胞（グリア細胞）
９…基体としての樹脂プレート
１０…第１参照物質としてのリファレンス部材
１３…超音波振動子としての超音波トランスデューサ
３１…第１演算手段、第２演算手段、画像構築手段、プロセッサとしてのＣＰＵ
５１…微小伝送路
６４…血管
６５…血液
Ｍ１…第１参照物質としての培養液
Ｍ２…第２参照物質としての培養液
Γ０、Γ１…インパルス応答
8A... Cultured cells (glial cells) as measurement objects
9 Resin plate as
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)
前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内及び前記第２参照物質内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する推定ステップと、
前記推定ステップにて得た前記第２参照物質の固有音響インピーダンスの推定値を前記第２参照物質の固有音響インピーダンスの実際値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する補正ステップと、
前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築ステップと
を有することを特徴とする超音波画像構築方法。 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.
前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する推定ステップと、
前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値、または、前記測定対象物が生体の軟組織である場合の前記測定対象物内における特定深度の固有音響インピーダンスの推定値の平均を仮想参照部位の固有音響インピーダンス値と定義し、前記推定ステップにて得た前記測定対象物内における前記特定深度の固有音響インピーダンスの推定値を前記仮想参照部位の固有音響インピーダンス値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する補正ステップと、
前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築ステップと
を有することを特徴とする超音波画像構築方法。 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.
前記第１参照物質及び前記第２参照物質は、前記培養細胞を培養するための培養液であり、
前記基体は、前記培養液を収容する培養容器の一部である
ことを特徴とする請求項１に記載の超音波画像構築方法。 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.
前記基体に測定対象物及び第１参照物質が接して存在し、かつ前記測定対象物における前記基体とは反対の側に第２参照物質が接して存在した状態で超音波を送信し、前記基体を介して超音波を入射させたときの前記測定対象物、前記第１参照物質及び前記第２参照物質からの超音波波形のインパルス応答を受信する超音波振動子と、
前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内及び前記第２参照物質内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する第１演算手段と、
前記第１演算手段により得た前記第２参照物質の固有音響インピーダンスの推定値を前記第２参照物質の固有音響インピーダンスの実際値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する第２演算手段と、
前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築手段と
を備えたことを特徴とする超音波画像構築装置。 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.
前記基体に測定対象物及び第１参照物質が接して存在した状態で超音波を送信し、前記基体を介して超音波を入射させたときの前記測定対象物及び前記第１参照物質からの超音波波形のインパルス応答を受信する超音波振動子と、
前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する第１演算手段と、
前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値、または、前記測定対象物が生体の軟組織である場合の前記測定対象物内における特定深度の固有音響インピーダンスの推定値の平均を仮想参照部位の固有音響インピーダンス値と定義し、前記第１演算手段により得た前記測定対象物内における前記特定深度の固有音響インピーダンスの推定値を前記仮想参照部位の固有音響インピーダンス値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する第２演算手段と、
前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築手段と
を備えたことを特徴とする超音波画像構築装置。 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.
基体に測定対象物及び第１参照物質が接して存在し、かつ前記測定対象物における前記基体とは反対の側に第２参照物質が接して存在した状態で超音波振動子に超音波を送信させ、前記基体を介して超音波を入射させたときの前記測定対象物、前記第１参照物質及び前記第２参照物質からの超音波波形のインパルス応答を前記超音波振動子に受信させる送受信ステップと、
前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内及び前記第２参照物質内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する推定ステップと、
前記推定ステップにて得た前記第２参照物質の固有音響インピーダンスの推定値を前記第２参照物質の固有音響インピーダンスの実際値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する補正ステップと、
前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築ステップと
を実行させるための超音波画像構築プログラム。 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.
基体に測定対象物及び第１参照物質が接して存在した状態で超音波振動子に超音波を送信させ、前記基体を介して超音波を入射させたときの前記測定対象物及び前記第１参照物質からの超音波波形のインパルス応答を前記超音波振動子に受信させる送受信ステップと、
前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する推定ステップと、
前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値、または、前記測定対象物が生体の軟組織である場合の前記測定対象物内における特定深度の固有音響インピーダンスの推定値の平均を仮想参照部位の固有音響インピーダンス値と定義し、前記推定ステップにて得た前記測定対象物内における前記特定深度の固有音響インピーダンスの推定値を前記仮想参照部位の固有音響インピーダンス値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する補正ステップと、
前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築ステップと
を実行させるための超音波画像構築プログラム。 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.
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)
| 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)
| 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 |
-
2019
- 2019-05-21 JP JP2019095229A patent/JP7213491B2/en active Active
Patent Citations (4)
| 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 | |
| CN103648400B (en) | Diagnostic ultrasound equipment and method | |
| JP6053512B2 (en) | Subject information acquisition apparatus and control method thereof | |
| JP6361001B1 (en) | Ultrasound image construction method, ultrasound image construction device, ultrasound image construction program, skin evaluation method for cosmetic purposes | |
| 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 | |
| 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 | |
| JP2022134005A (en) | Ultrasonic image generating method and apparatus, and signal processing method | |
| JP2017039031A (en) | Subject information acquisition device and control method thereof | |
| JP2024078031A (en) | Ultrasound imaging method and apparatus, and signal processing method | |
| JP2018199002A (en) | Information processing device and control method thereof | |
| 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 |