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.
超音波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.
しかしながら、上記装置では、奥行方向に沿って深くなる方向、つまり測定対象物の前面側から背面側に向かう方向に固有音響インピーダンス分布を順次推定していく手法を採用している。このため、深度が大きくなるほど固有音響インピーダンスの推定値の誤差が蓄積、拡大するという欠点があり、より精度の高い画像を構築するうえでのマイナス要因となっていた。また、このような場合、超音波断層像における奥行方向に多くの筋が入ってしまうため、微細な内部の層構造が感覚的に理解しやすいとは言い難く、いまだ改良の余地があった。 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
[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
[5]前記測定対象物は、培養細胞であり、前記第1参照物質及び前記第2参照物質は、前記培養細胞を培養するための培養液であり、前記基体は、前記培養液を収容する培養容器の一部であることを特徴とする請求項1に記載の超音波画像構築方法。
[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]前記測定対象物は、組織表面に対して略平行な方向に延びる血管を含む生体の軟組織としての皮膚であり、前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値は、前記血管中に存在する血液の固有音響インピーダンス値であることを特徴とする請求項2に記載の超音波画像構築方法。
[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]前記測定対象物が生体の軟組織である場合の前記測定対象物内における特定深度の固有音響インピーダンスの推定値の平均は、前記測定対象物が皮膚である場合の当該皮膚内における特定深度の固有音響インピーダンスの推定値の平均であり、前記特定深度は、前記音響物性像の画像データが構築される領域よりも深い位置に設定されることを特徴とする請求項2に記載の超音波画像構築方法。
[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]基体と、前記基体に測定対象物及び第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
[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
上記手段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の実施形態を図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
パルス励起型超音波顕微鏡2は、ステージ4を有する顕微鏡本体5と、ステージ4の下方に設置された超音波プローブ6とを備える。パルス励起型超音波顕微鏡2の超音波プローブ6は、PC3と電気的に接続されている。
A pulse excitation
本実施形態のステージ4は、ユーザの手動操作により、水平方向(即ちX方向及びY方向)に移動できるように構成されている。このステージ4には、測定対象物を接触させて配置するための樹脂プレート9が固定されている。ここでの測定対象物は、組織表面に対して略平行な方向に延びる血管を含む生体の軟組織(具体的には、肌組織;皮膚8)である。本実施形態では、ヒトの皮膚8を樹脂プレート9に直接押し付けることにより測定等を行っている。また、既知の音響物性を有する基体としての樹脂プレート9は、超音波を透過させることができる平板状部材であって、測定対象物である皮膚8よりも硬い材料からなる。このような形状及び硬さの部材を基体として用いた場合、測定対象物である皮膚8を確実に密着配置することが可能となり、奥行方向の固有音響インピーダンス分布を正確に推定可能となる結果、ひいては画像構築の精度が向上する。なお、本実施形態では厚さ1.4mmのポリスチレン板が用いられている。勿論、ポリスチレン以外の樹脂からなる板材などを用いることも許容される。
The
この樹脂プレート9において皮膚8が接触配置される側である上面には、第1参照物質としてのリファレンス部材10があらかじめ設置されている。リファレンス部材10は、樹脂プレート9とは異なる既知の音響物性を有している。本実施形態では例えばアクリル樹脂(アクリル接着剤)を付着させることによりリファレンス部材10としているが、勿論これに限定されるわけではない。リファレンス部材10に対して密着させることが可能なものであれば、樹脂材以外のもの(例えばガラス材、金属材、セラミック材など)をリファレンス部材10としてもよい。あるいは、このようなリファレンス部材10を設置する代わりに、例えば樹脂プレート9の上面に接するように水等を存在させておき、これを第1参照物質として用いてもよい。なお、樹脂プレート9にリファレンス部材10をあらかじめ設置しておくことにより、装置が置かれる環境の変化等に依存せず、リファレンス部材10に入射した超音波波形のインパルス応答情報を正確にかつ安定的に取得することができる。
A
超音波プローブ6は、水などの超音波伝達媒体Wを貯留可能な貯留部11をその先端部に有するプローブ本体12と、プローブ本体12の略中心部に配置される超音波トランスデューサ13(超音波振動子)と、プローブ本体12を前記ステージ4の面方向に沿って二次元的に走査するためのX-Yステージ14とを備える。プローブ本体12の貯留部11は上部が開口しており、その貯留部11の開口側を上向きにした状態で超音波プローブ6がステージ4の下方に設置されている。
The
超音波トランスデューサ13は、例えば酸化亜鉛の薄膜圧電素子16とサファイアロッドの音響レンズ17とによって構成される。この超音波トランスデューサ13は、パルス励起されることで樹脂プレート9の下面側から皮膚8及びリファレンス部材10に対して超音波を照射する。超音波トランスデューサ13が照射する超音波は、貯留部11の超音波伝達媒体Wを介して円錐状に収束されて樹脂プレート9の上面(皮膚8の表面付近)で焦点を結ぶようになっている。なお本実施形態では、超音波トランスデューサ13として、口径1.2mm、焦点距離1.5mm、中心周波数80MHz、帯域幅50~105MHz(-6dB)の仕様のものを用いている。
The
図2は、本実施形態の超音波画像構築装置1の電気的な構成を示すブロック図である。
FIG. 2 is a block diagram showing the electrical configuration of the ultrasonic
図2に示されるように、超音波プローブ6は、超音波トランスデューサ13、X-Yステージ14、パルス発生回路21、受信回路22、送受波分離回路23、検波回路24、A/D変換回路25、エンコーダ26、コントローラ27を備える。
As shown in FIG. 2, the
走査手段としてのX-Yステージ14は、超音波の照射点を二次元的に走査させるためのXステージ14X及びYステージ14Yを備えるとともに、それぞれのステージ14X,14Yを駆動するモータ28X,28Yを備えている。これらのモータ28X,28Yとしては、ステッピングモータやリニアモータが使用される。
The
各モータ28X,28Yにはコントローラ27が接続されており、該コントローラ27の駆動信号に応答してモータ28X,28Yが駆動される。これらモータ28X,28Yの駆動により、Xステージ14Xを連続走査(連続送り)するとともに、Yステージ14Yを間欠送りとなるよう制御することで、X-Yステージ14の高速走査が可能となっている。
A
また、本実施形態においては、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
さらに、コントローラ27は、駆動制御信号に同期してトリガ信号を生成してパルス発生回路21に供給する。これにより、パルス発生回路21において、そのトリガ信号に同期したタイミングで励起パルスが生成される。その励起パルスが送受波分離回路23を介して超音波トランスデューサ13に供給される結果、超音波トランスデューサ13から超音波が照射される。
Furthermore, the
図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
超音波トランスデューサ13の薄膜圧電素子16は、送受波兼用の超音波振動子であり、皮膚8で反射した超音波(反射波)を電気信号に変換する。そして、その反射波の信号は、送受波分離回路23を介して受信回路22に供給される。受信回路22は、信号増幅回路を含んで構成されていて、反射波の信号を増幅して検波回路24に出力する。
The thin-
検波回路24は、皮膚8からの反射波信号を検出するための回路であり、図示しないゲート回路を含む。本実施形態の検波回路24は、超音波トランスデューサ13で受信した反射波信号のなかから、皮膚8からの反射波信号やリファレンス部材10からの反射波信12月号を抽出する。そして、検波回路24で抽出された反射波信号は、A/D変換回路25に供給されてA/D変換された後、PC3に転送される。
The
PC3は、CPU31(中央処理装置)、I/F回路32、メモリ33、記憶装置34、入力装置35、及び表示装置36を備え、それらはバス37を介して相互に接続されている。
The
CPU31は、メモリ33を利用して制御プログラムを実行し、システム全体を統括的に制御する。制御プログラムとしては、X-Yステージ14による二次元走査を制御するためのプログラム、超音波Bモードエコー像の元となる反射信号列のデータを固有音響インピーダンス像へ変換するためのプログラム、固有音響インピーダンス像を表示するためのプログラムなどを含む。なお、CPU31とは別に例えばDSP(Digital Signal Processor:デジタル信号プロセッサ)を設けて、そこでCPU31が行っている信号処理の一部を行わせてもよい。
The
I/F回路32は、超音波プローブ6との間で信号の授受を行うためのインターフェース(具体的には、USBインターフェース)である。I/F回路32は、超音波プローブ6に制御信号(コントローラ27への駆動制御信号)を出力したり、超音波プローブ6からの転送データ(A/D変換回路25から転送されるデータなど)を入力したりする役割を果たすものである。なお、超音波プローブ6との間で信号の授受を行う場合には、上記のような物理的なインターフェースに限定されることはなく、無線インターフェースを用いてもよい。
The I/
表示装置36は、例えば、液晶、プラズマ、有機EL(electroluminescence)等のモニタディスプレイである。表示装置36は、カラー表示、モノクロ表示を問わずに使用できるが、カラー表示であることが望ましい。この表示装置36は、皮膚8の表層の固有音響インピーダンス像を表示したり、各種設定の入力画面を表示したりするために用いられる。
The
入力装置35は、タッチパネル、マウス、キーボード、ポインティングデバイス等の入力ユーザインタフェースであって、ユーザからの要求や指示、パラメータの入力に用いられる。
The
記憶装置34は、磁気ディスク装置や光ディスク装置などのハードディスクドライブであり、各種の制御プログラム及び各種のデータを記憶している。メモリ33は、RAM(ランダムアクセスメモリ)やROM(リードオンリーメモリ)を含み、超音波測定のためにあらかじめ取得されたリファレンス部材10の反射波形とその固有音響インピーダンスとを保存する。CPU31は、入力装置35による指示に従い、プログラムやデータを記憶装置34からメモリ33へ転送し、それを逐次実行する。なお、CPU31が実行するプログラムとしては、メモリカード、フレキシブルディスク、光ディスクなどの記憶媒体に記憶されたプログラムや、通信媒体を介してダウンロードしたプログラムでもよく、その実行時には記憶装置34にインストールして利用する。
The
次に、本実施形態の超音波画像構築装置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
この超音波画像構築装置1では、第1参照物質であるリファレンス部材10に入射した超音波波形のインパルス応答情報、及び測定対象物である皮膚8に入射した超音波波形のインパルス応答情報から、規格化されたインパルス応答情報を得て、その規格化されたインパルス応答情報に基づいて奥行方向の音響物性分布を多重反射の影響を考慮して推定するようになっている。また、このような推定を行うために、本実施形態では、測定対象物内において、異なる固有音響インピーダンスを持つ無損失の微小伝送路51が奥行方向に連なって伝送路の集合体をなしていると仮定して、手前側の微小伝送路51の固有音響インピーダンスの推定結果に基づきその奥側に隣接する微小伝送路51の固有音響インピーダンスを推定する処理を順次繰り返すことにより、伝送路の奥行方向の音響物性分布(ここでは固有音響インピーダンス分布)を推定する演算を行うようになっている。このような演算は、メモリ33内に格納された所定のアルゴリズムに基づいて実行される。
In this ultrasonic
このアルゴリズムは、超音波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に超音波の収束ビームを入射させ、そこからの反射波形を取得する。その時のインパルス応答Γ0(ω)は、フーリエ変換を用いて、入射波S0と皮膚8からの反射波Stgt(ω)とから次式1のように表される。
この場合、固有音響インピーダンスが既知かつ均一であって、測定対象物に比べて十分な厚さを有するリファレンス部材10からの反射波形も取得する必要がある。リファレンス部材10からの反射波Sref(ω)は、リファレンス部材10の固有音響インピーダンスZrefと樹脂プレート9の固有音響インピーダンスZ0とを用いて、次式2のように表される。
また、測定対象物である皮膚8からのインパルス応答Γ0(ω)は次式3で表される。ただし、このインパルス応答Γ0(ω)には皮膚8の組織の奥の複数の界面から発生する反射が含まれているため、インパルス応答Γ0(ω)は周波数特性を持ったものとなる。なお、ここまでの式により、第1参照物質に入射した超音波波形のインパルス応答情報及び測定対象物に入射した超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報が求められる。
ここで、図6は、樹脂プレート9に接する微小伝送路51から順に、各微小伝送路51の特性インピーダンスZ1、Z2…Znを推定していく様子を概念的に示した図である。この図に示されるように、各微小伝送路51の特性インピーダンスZ1、Z2…Znを、樹脂プレート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は、インパルス応答Γ0(ω)を逆フーリエ変換したg0(t)を表したものであるが、その第1項は多重反射の影響を受けない(図7参照)。従って、この第1項の値から樹脂プレート9に接する微小伝送路51の特性インピーダンスZ1を、次式5のように推定することができる。
周波数領域の皮膚8の固有音響インピーダンスZx0は、Γ0を用いると次式6のように表される。
また、Zx0はさらに奥からのインパルス応答Γ1を用いて次式7のようにも表される。
ここで、次式8、9にて表すように、γは伝播定数、αは減衰定数、βは位相定数、fは周波数であるが、本実施形態のアルゴリズムではα=0、及び皮膚8の全微小伝送路51の音速をc=1600(m/s)と仮定している。
また、各微小伝送路51の距離Δlは次式10のように表され、ここでも音速をc=1600(m/s)と仮定している。Δtは皮膚8からの反射波形のサンプリング間隔の1ポイントに相当する(本実施形態ではΔt=2(ns))。
そして、上記の式をもとに、さらに奥にある微小伝送路51からのインパルス応答Γ1を求めることができる(次式11)。即ち、Zx0及びZ1の値をもとに、Z1の終点におけるΓ1の値を推定することができる。
次式12は、インパルス応答Γ1(ω)を逆フーリエ変換したg1(t)を表したものであって、その第1項は多重反射の影響を受けない。従って、この第1項の値から当該微小伝送路51に隣接するさらに奥側の微小伝送路51の特性インピーダンスZ2を推定することができ、同様にZx1、Γ2(ω)も推定することができる(次式13,14,15)。
この工程を繰り返すことによって、各微小伝送路51の特性インピーダンス(固有音響インピーダンス)Z1、Z2…Znを推定することができる。
By repeating this process, the characteristic impedances (specific acoustic impedances) Z 1 , Z 2 . . . Zn of each
本実施形態のアルゴリズムでは、上記のような推定ステップを行った後、さらに以下に示す所定の補正ステップを実施する。即ち、この補正ステップでは、まず仮想参照部位の固有音響インピーダンス値を定義する。ここでの測定対象物は皮膚8(肌組織)であって(図8参照)、皮膚8を構成する表皮61及び真皮62よりも深い位置にある皮下組織63の比較的浅い部位には、肌組織表面に対して略平行な方向に延びる血管64が存在している。血管64内に存在する血液65の固有音響インピーダンス値は既知(1.61MNs/m3)であるため、本実施形態ではこの値を上記「仮想参照部位の固有音響インピーダンス値」として用いている。 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/m3)である。固有音響インピーダンス曲線C1における左側の領域が、超音波トランスデューサ13に近い位置にある樹脂プレート9の固有音響インピーダンス値Z(x)を示す部分である。一方、樹脂プレート9に接して配置された皮膚8は、ポリスチレン製の樹脂プレート9よりも軟らかい。このため、その固有音響インピーダンス値Z(x)は、固有音響インピーダンス曲線C1における点Aを境にして急激に低くなる。固有音響インピーダンス曲線C1における右側の領域が皮膚8の固有音響インピーダンス値Z(x)の補正前の推定値を示す部分であり、2.50MNs/m3よりもかなり低い値となっている。当該曲線C1において点B(即ち皮膚8内における特定深度に相当する部分)よりも右側はフラットであり、この部分が特定深度の血管内に存在する血液の固有音響インピーダンス値の補正前の推定値を示す部分となっている。ただし、当該推定値は誤差を含んでおり、実際の値とはいくぶん異なるものとして表されている(この例では1.80MNs/m3)。そこで、特定深度の血管64内に存在する血液65の固有音響インピーダンス値の補正前の推定値を、仮想参照部位の固有音響インピーダンス値(1.61MNs/m3)に置き換えるべく、1.80MNs/m3から補正値である0.19MNs/m3を減算する。その結果、当該曲線C1のフラットな部分が下方に移動し、これが補正された固有音響インピーダンス曲線C2のうちの一部をなすものとなる。なお、当該曲線C1における点Bは、曲線C2における点Cに移動した状態となる。
Here, since the
次に、固有音響インピーダンス曲線C1における点Aと点Bとの間の領域(即ち皮膚表面から当該血管までの部分に相当する領域)を補正する。その際、先の補正値である0.19MNs/m3を最大とし、測定対象物の前面側(即ち皮膚表面)に向かうほど小さくなる補正値を用いる。そしてこの補正値を減算することにより、当該血管から皮膚表面までの固有音響インピーダンスの推定値を補正する。換言すると、皮膚表面からの深さに応じた補正値を設定し、推定値からこれを減算する。例えば、特定深度の半分の深さ位置における推定値を補正する場合には、0.19MNs/m3の1/2の値(即ち0.095MNs/m3)をその位置における補正値として設定し、これをその位置における推定値から減算する。また、特定深度の1/4の深さ位置における推定値を補正する場合には、0.19MNs/m3の1/4の値(即ち0.0475MNs/m3)をその位置における補正値として設定し、これをその位置における推定値から減算する。ちなみに、図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
まず、測定対象物であるヒトの皮膚8(例えば、比較的浅い位置に太い血管(けい静脈、けい動脈)が存在している首の皮膚8など)を樹脂プレート9の上面に押し付けるようにして接触配置させる。この状態で、まず超音波プローブ6に初期動作を行わせる。即ち、CPU31からの指示に基づいてコントローラ27を作動させることにより、モータ28X,28Yを駆動し、リファレンス部材10がある位置にて超音波照射が行われるようにX-Yステージ14を移動させる。
First, the human skin 8 (for example,
またこのとき、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
その後、CPU31からの指示に基づいてコントローラ27によりモータ28X,28Yが駆動され、X-Yステージ14による二次元走査が開始される。CPU31は、エンコーダ26の出力に基づいて測定点の座標データを取得する(ステップS110)。
After that, the
そして、図4(a)に示すように、CPU31からの指示に基づいて励起パルスがトランスデューサ13に供給されることにより、皮膚8に超音波Soが照射され、その反射波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
次いで、第1演算手段としてのCPU31は、規格化されたインパルス応答信号のデータを用いて、上記のアルゴリズムのうち、TDR法の原理を参考とした推定ステップの演算を実行する。そしてCPU31は、その演算により皮膚8における測定点での奥行方向の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定し、その推定結果を座標データに関連付けてメモリ33に記憶する(ステップS130)。
Next, the
次いで、第2演算手段としてのCPU31は、上記のアルゴリズムのうちの補正ステップの演算を実施する。即ち、仮想参照部位の固有音響インピーダンス値を定義した後、皮膚8内における特定深度の固有音響インピーダンスの推定値を前記仮想参照部位の固有音響インピーダンス値に置き換える演算を行い、この演算を通じて皮膚8内の固有音響インピーダンスの推定値を補正する。そしてCPU31は奥行方向の固有音響インピーダンス分布を推定し、その推定結果を座標データに関連付けてメモリ33に記憶する(ステップS132)。
Next, the
その後、画像構築手段としてのCPU31は、奥行方向の固有音響インピーダンス分布の推定結果に基づいて、固有音響インピーダンス像(断層像)を構築するための画像処理を行う(ステップS140)。詳しくは、CPU31は、固有音響インピーダンス分布の推定結果に基づいてカラー変調処理を行い、固有音響インピーダンスの大きさに応じて色分けして表示した画像データを構築し、該画像データをメモリ33に記憶する。
After that, the
次いで、CPU31は、全ての測定点での処理が終了して、全ての測定点で画像データが取得されたか否かを判断する(ステップS150)。ここで、全データが取得されていない場合には(ステップS150:NO)、CPU31は、ステップS110に戻って、ステップS110~S140の処理を繰り返して実行する。全データが取得された場合には(ステップS150:YES)、CPU31は、次ステップS160に移行する。
Next, the
そして、CPU31は、該データを表示装置36に転送し、あらかじめ定めた直線上における固有音響インピーダンス像(断層像)を表示させた後(ステップS160)、図10の処理を終了するようになっている。このような一連の処理により、皮膚8での固有音響インピーダンスの大きさに応じて色分けされた固有音響インピーダンス像(断層像)が表示される。
Then, the
従って、本実施の形態によれば以下の効果を得ることができる。 Therefore, according to this embodiment, the following effects can be obtained.
(1)本実施形態の超音波画像構築装置1では、リファレンス部材10に入射した超音波波形のインパルス応答情報、及び皮膚8に入射した超音波波形のインパルス応答情報から、規格化されたインパルス応答情報、つまり装置に依存しないインパルス応答情報を得るようにしている。そして、この規格化されたインパルス応答情報に基づいて、皮膚8内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって、多重反射の影響を考慮しつつ順次推定する演算(推定ステップ)を行うようにしている。本実施形態ではこれに加えて、皮膚8内における特定深度(ここでは、横行する血管が存在する皮下組織浅部)の固有音響インピーダンスの推定値を、仮想参照部位の固有音響インピーダンス値に置き換える演算を通じて、皮膚8内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する演算(補正ステップ)も行うようにしている。以上の結果、深度が大きくなるほど蓄積、拡大しやすい固有音響インピーダンスの推定値の誤差を減じることができ、奥行方向の固有音響インピーダンス分布のより正確な推定が可能となる。また、補正ステップを行わない従来技術では超音波断層像における奥行方向に多くの筋が入ってしまうため、微細な層構造が感覚的に理解しにくかったのに対し、本実施形態ではこのような筋が少なくなり層方向の色ムラが減少する。従って、微細な層構造を有する非常に薄い皮膚8の超音波断層像を、その層構造が感覚的に理解しやすい固有音響インピーダンス像として比較的簡単にかつ高い精度で構築することができる。なお、この装置1により得られる固有音響インピーダンス像は、測定対象物を切断することなく(即ち非侵襲で)、層ごとの力学特性の断面分布(深さ分布)情報を、推定した固有音響インピーダンスの絶対値ごとに色分けして画像化したものである。それゆえ、この像は感覚的に層構造が理解しやすいものなっている。
(1) In the ultrasonic
ここで、一般に通常の超音波診断装置で得られる超音波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
(2)本実施形態の超音波画像構築装置1では、皮膚組織内部での多重反射を考慮した時間-周波数領域における解析を通じて、超音波Bモードエコー像の元となる反射信号列を奥行方向の固有音響インピーダンス像に変換する手法を採用している。このため、比較的簡単に所定の演算を行うことができる。また本実施形態では、上述したアルゴリズムを用いて、手前側の微小伝送路51の固有音響インピーダンスの推定結果に基づきその奥側に隣接する微小伝送路51の固有音響インピーダンスを推定する処理を順次繰り返すことにより、伝送路の奥行方向の固有音響インピーダンス分布を推定する演算を行うようにしている。従って、この演算によれば、多重反射の影響を最小限にすることができ、奥行方向の固有音響インピーダンス分布を精度よく推定することができる。
(2) In the ultrasonic
[第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/m3)を「仮想参照部位の固有音響インピーダンス値」として用いたのに対し、ここでは測定対象物が生体の軟組織であることに着目してその特定深度の固有音響インピーダンスの推定値の平均を「仮想参照部位の固有音響インピーダンス値」と定義して用いている。本実施形態では、肌組織表面から350μmほど深い位置における組織の固有音響インピーダンスの推定値の平均(例えば1.65MNs/m3)を求めて、その値を用いている。これは、生体の軟組織中には骨や軟骨などの硬い部分が存在せず硬さのばらつきが少ないので、その固有音響インピーダンスは例えば1.50MNs/m3~1.70MNs/m3の範囲から大きく外れた値とはならず、仮想参照部位として選択するのに適している、との考え方に立脚したものである。なお、脂肪を除く軟組織であれば、固有音響インピーダンスは、ほぼ1.60MNs/m3~1.70MNs/m3の範囲内となる。
In the ultrasonic
この場合、皮膚8における特定深度は、固有音響インピーダンス像の画像データが構築される領域よりも深い位置に設定されることが好ましい。本実施形態では、例えば肌組織表面から300μmほど深い位置までを画像構築領域として設定している関係上、組織の固有音響インピーダンスの推定値の平均を求めるための特定深度はそれよりも50μm深い位置としている。このようにすることで、特定深度の固有音響インピーダンス像を敢えて表示しないようにし、かえって不自然なものとならないようにしている。
In this case, the specific depth in the
ここで、図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
[第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
ステージ4には、測定対象物である培養細胞8Aを培養するための培養液Mを収容する培養容器102が設置されている。ここでは、培養細胞8Aとして、接着性細胞の一種であるヒトのグリア細胞8Aを用いている。グリア細胞8Aの厚さは非常に薄く、数μm程度である。この培養容器102の底部中央部には、上面側に培養細胞8Aが培養時に接着して支持される樹脂プレート9が固定されている。既知の音響物性を有する基体としての樹脂プレート9は、超音波を透過させることができる平板状部材であって、測定対象物である培養細胞8Aよりも硬い材料からなる。
A
なお、培養容器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
次に、本実施形態の超音波画像構築装置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
この超音波画像構築装置101では、基体である樹脂プレート9に測定対象物である培養細胞8A及び第1参照物質である培養液M1が接して存在し、かつ培養細胞8Aにおける樹脂プレート9とは反対の側に第2参照物質である培養液M2が接して存在した状態で超音波を送信する。そして、樹脂プレート9を介して超音波を入射させたときの培養細胞8A、第1参照物質である培養液M1及び第2参照物質である培養液M2からの超音波波形のインパルス応答をそれぞれ受信する(送受信ステップ)。
In this ultrasonic
次いで、培養液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
本実施形態の超音波画像構築装置101では、上記の推定ステップ後に補正ステップを行ってから画像構築ステップを行っており、この点については上記第1の実施形態と共通する(上記図9、図10等を参照)。ただし、本実施形態の補正ステップでは第1の実施形態のときとは異なる以下の補正ステップを実施する。即ちここでは、推定ステップにて得た第2参照物質(培養液M2)の固有音響インピーダンスについて、その「推定値」を「実際値」に置き換える補正演算を行う。培養液M2の前記「実際値」が例えば1.55MNs/m3であり、誤差を含んでいる可能性が高い前記「推定値」が例えば1.80MNs/m3であるとすると、その値から補正値である0.25MNs/m3を減算することにより、前記「実際値」を補正する。その結果、先の実施形態で示した図9のグラフに示す補正前の固有音響インピーダンス曲線C1のフラットな部分が下方に移動し、これが補正された固有音響インピーダンス曲線C2のうちの一部をなすものとなる。なお、当該曲線C1における点Bは、曲線C2における点Cに移動した状態となる。
In the ultrasonic
次に、固有音響インピーダンス曲線C1における点Aと点Bとの間の領域を補正する。その際、先の補正値である0.25MNs/m3を最大とし、測定対象物の前面側(即ち培養細胞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
そして、このような補正ステップの演算を通じて、測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定し、最終的には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
図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
なお、本発明の各実施の形態は以下のように変更してもよい。 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
・上記第1、第2実施形態の超音波画像構築装置1では、下方から超音波を照射する倒立型の超音波顕微鏡2を用いて超音波の照射を行ったが、上方から超音波を照射する正立型の超音波顕微鏡を用いてもよい。
In the ultrasonic
・上記第1、第2実施形態では、基本的に疾患を有していない比較的健康な皮膚8を対象として、その状態を評価する目的で超音波画像構築装置1を用いたが、これに限定されない。例えば、皮膚がんなどの疾患に伴う皮膚の異常を早期に検出する目的で超音波画像構築装置1を用いてもよい。
- In the above-described first and second embodiments, the ultrasonic
・上記第1、第2実施形態の超音波画像構築装置1では、測定対象物がヒトの首の皮膚8であったが、首以外の他の部位(例えば頬など)における皮膚8であっても勿論よい。また、測定対象物は皮膚8でなくてもよく、例えば内臓、筋肉、脳、歯、爪、骨の表層部などであっても勿論よい。また、上記第3実施形態の超音波画像構築装置101では、測定対象物がヒトのグリア細胞であったが、付着性細胞であれば生物種は問わずどのような培養細胞であってもよい。さらにいうと、測定対象物は必ずしも生体組織や生物でなくてもよく、非生物(例えば塗膜など)であってもよい。換言すると、本発明の超音波画像構築装置1、101は医療分野、美容分野、化粧品分野のみに限定されず、例えば工業分野などの分野においても使用されることができる。
In the ultrasonic
・上記実施形態の超音波画像構築装置1、101は、いずれも測定対象物に対して超音波トランスデューサ13を二次元方向に相対的に走査させる走査手段を備えていたが、これに代えて超音波トランスデューサ13を一次元方向にのみ相対的に走査させる走査手段を備えたものとしてもよい。また、走査手段は必須の構成ではないため省略しても勿論よく、この場合には装置を小型化、簡略化、低コスト化することが可能となる。
The ultrasonic
・上記実施形態の超音波画像構築装置1、101では、奥行方向の固有音響インピーダンス分布の推定結果に基づいて固有音響インピーダンス像を構築したが、これに限定されない。例えば、奥行方向の音速分布を推定し、その結果に基づいて音速像を構築してもよい。
In the ultrasonic
・上記実施形態の超音波画像構築装置1、101では、超音波Bモードエコー画像の元となる反射信号列から固有音響インピーダンス像を構築してそれを表示装置36に表示させるように構成したが、固有音響インピーダンス像ばかりでなく超音波Bモードエコー像も表示できるようにしても勿論よい。また、超音波Bモードエコー像を表示する汎用の超音波診断装置に上記実施形態のアルゴリズムを組み込むことで、超音波画像構築装置1、101として動作させるようにしてもよい。
In the ultrasonic
・上記実施形態では、測定対象物内での多重反射を考慮した時間-周波数領域における解析を通じて、超音波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参照物質としての培養液
Γ0、Γ1…インパルス応答
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)
前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内及び前記第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.
前記超音波波形のインパルス応答情報から得た規格化されたインパルス応答情報に基づき、前記測定対象物内の固有音響インピーダンスを奥行方向の手前側から奥側に向かって順次推定する推定ステップと、
前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値、または、前記測定対象物が生体の軟組織である場合の前記測定対象物内における特定深度の固有音響インピーダンスの推定値の平均を仮想参照部位の固有音響インピーダンス値と定義し、前記推定ステップにて得た前記測定対象物内における前記特定深度の固有音響インピーダンスの推定値を前記仮想参照部位の固有音響インピーダンス値に置き換える演算を通じて、前記測定対象物内の固有音響インピーダンスの推定値を補正したうえで、奥行方向の固有音響インピーダンス分布を推定する補正ステップと、
前記奥行方向の固有音響インピーダンス分布に基づいて、音響物性像の画像データを構築する画像構築ステップと
を有することを特徴とする超音波画像構築方法。 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.
前記第1参照物質及び前記第2参照物質は、前記培養細胞を培養するための培養液であり、
前記基体は、前記培養液を収容する培養容器の一部である
ことを特徴とする請求項1に記載の超音波画像構築方法。 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.
前記測定対象物内における特定深度にて既知の固有音響インピーダンス値を有する物質が層方向にわたり均一に存在している場合の当該固有音響インピーダンス値は、前記血管中に存在する血液の固有音響インピーダンス値である
ことを特徴とする請求項2に記載の超音波画像構築方法。 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:
前記特定深度は、前記音響物性像の画像データが構築される領域よりも深い位置に設定される
ことを特徴とする請求項2に記載の超音波画像構築方法。 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.
基体に測定対象物及び第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.
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 | |
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 |