JPS5947A - 超音波診断装置 - Google Patents
超音波診断装置Info
- Publication number
- JPS5947A JPS5947A JP10791982A JP10791982A JPS5947A JP S5947 A JPS5947 A JP S5947A JP 10791982 A JP10791982 A JP 10791982A JP 10791982 A JP10791982 A JP 10791982A JP S5947 A JPS5947 A JP S5947A
- Authority
- JP
- Japan
- Prior art keywords
- points
- curvature
- ultrasonic diagnostic
- arbitrary
- arcuate portion
- 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.)
- Pending
Links
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
本発明は生体等の内部の診断を行なう際に用いられる超
音波診断装置に関し、超音波画像を形成する曲線部分の
曲率半径を簡単な操作により極めて容易な演算によって
高速に算出することのできる超音波診断装置を提供する
ものである０超音波診断装置のブロック図を第１図に示
す。
音波診断装置に関し、超音波画像を形成する曲線部分の
曲率半径を簡単な操作により極めて容易な演算によって
高速に算出することのできる超音波診断装置を提供する
ものである０超音波診断装置のブロック図を第１図に示
す。
１は超音波探触子４に印加する電気信号を送信す用して
いるため送信、受信のタイミング切換を行なう送受信切
換回路、４は電気信号と音響信号間の変換を行ない超音
波パルスを送受信する超音波探触子である。６は高周波
増幅回路２で増幅された高周波からその包絡線をとり出
す検波回路、６は検波回路５の出力を映像信号に変換す
るための映像増幅回路、７は超音波画像を表示するため
のＣＲＴなどの表示部、８は送信、受信のタイミングコ
ントロールを行なったり、音響走査からＴＶ走査への変
換用クロック信号の発生などを行なうコントロール回路
、９はキーボード、亀子キャリパー１画像フリーズなど
の操作を行なうパネル操作部である。
いるため送信、受信のタイミング切換を行なう送受信切
換回路、４は電気信号と音響信号間の変換を行ない超音
波パルスを送受信する超音波探触子である。６は高周波
増幅回路２で増幅された高周波からその包絡線をとり出
す検波回路、６は検波回路５の出力を映像信号に変換す
るための映像増幅回路、７は超音波画像を表示するため
のＣＲＴなどの表示部、８は送信、受信のタイミングコ
ントロールを行なったり、音響走査からＴＶ走査への変
換用クロック信号の発生などを行なうコントロール回路
、９はキーボード、亀子キャリパー１画像フリーズなど
の操作を行なうパネル操作部である。
パネル操作部９の一部を構成する、電子キャリバー回路
のブロック図を第２図に示す。１１．１２はそれぞれＣ
ｆｔＴなどの表示画面上の点をＸ方向。
のブロック図を第２図に示す。１１．１２はそれぞれＣ
ｆｔＴなどの表示画面上の点をＸ方向。
Ｙ方向に移動させるためのバイアス眼圧を発生させるボ
リュームであって、一般にはＸ方向、Ｙ方向の移動を同
時に行なうために１本のコントロールスティックが用い
られている。１３はＸ、Ｙ座標に対応した出力電圧を交
互に取り込むためのタイミング切換回路、１４はボリュ
ーム１１．１２によって発生されタイミング回路１３を
介して得られたＸ、Ｙ座標に対応したアナログ出力電圧
をデジタル信号に変換するためのＡ／Ｄコンバータ、１
６はデータ送受のコントロールや演算を行なうマイクロ
プロセッサ、１６はキャラクタジェネレータ１７に文字
コードおよび制御信号を送るＣＲＴコントローラ、１８
はＣＲＴなどの表示部である。１９はＡ／Ｄコンバータ
１３のＡ／Ｄ変換動作をスタート、あるいは中断させる
信号を与えるスイッチ、２０はＡ／Ｄコンバータ１３の
変換動作が中断された時にＡ／Ｄコンバータに保持され
ているデータを記憶するＲＡＭζ２１はデータ処理のた
めのデータ変換テーブルを内蔵するＲＯＭである。
リュームであって、一般にはＸ方向、Ｙ方向の移動を同
時に行なうために１本のコントロールスティックが用い
られている。１３はＸ、Ｙ座標に対応した出力電圧を交
互に取り込むためのタイミング切換回路、１４はボリュ
ーム１１．１２によって発生されタイミング回路１３を
介して得られたＸ、Ｙ座標に対応したアナログ出力電圧
をデジタル信号に変換するためのＡ／Ｄコンバータ、１
６はデータ送受のコントロールや演算を行なうマイクロ
プロセッサ、１６はキャラクタジェネレータ１７に文字
コードおよび制御信号を送るＣＲＴコントローラ、１８
はＣＲＴなどの表示部である。１９はＡ／Ｄコンバータ
１３のＡ／Ｄ変換動作をスタート、あるいは中断させる
信号を与えるスイッチ、２０はＡ／Ｄコンバータ１３の
変換動作が中断された時にＡ／Ｄコンバータに保持され
ているデータを記憶するＲＡＭζ２１はデータ処理のた
めのデータ変換テーブルを内蔵するＲＯＭである。
従来の超音波診断装置における電子キャリパ−を用いた
超音波画像の計測としては、生体組織の長さ、直径１面
積などの計測に多く用いられているだけで、曲線部分の
曲率半径の測定は行なわれていない。したがって、頚動
脈の超音波診断における頚動脈の蛇行の程度が有用な診
断情報を与えるものであるにもかかわらず、蛇行の程度
を定量的にとらえることは行なわれず、目視での観察に
よる診断がなされているに過ぎない。そこで、定量的な
計測を行なおうとする場合、曲線の曲率にはＸ−Ｖ平面
内の曲線ｙ ＝ ｆ （ｘ）に対して２ｙ ｘ２ で与えられる事が一般に知られているが、超音波診断画
像上における円弧状部の曲率を第（り式の数値計算によ
り求めることは、演算が複雑となり、超音波診断装置を
構成するマイクロプロセッサにこの演算を行なわせるこ
とは難しく、また演算時間も長くなり、素早い診断が要
求される超音波診断装置としては不都合が生じる結果と
なる。
超音波画像の計測としては、生体組織の長さ、直径１面
積などの計測に多く用いられているだけで、曲線部分の
曲率半径の測定は行なわれていない。したがって、頚動
脈の超音波診断における頚動脈の蛇行の程度が有用な診
断情報を与えるものであるにもかかわらず、蛇行の程度
を定量的にとらえることは行なわれず、目視での観察に
よる診断がなされているに過ぎない。そこで、定量的な
計測を行なおうとする場合、曲線の曲率にはＸ−Ｖ平面
内の曲線ｙ ＝ ｆ （ｘ）に対して２ｙ ｘ２ で与えられる事が一般に知られているが、超音波診断画
像上における円弧状部の曲率を第（り式の数値計算によ
り求めることは、演算が複雑となり、超音波診断装置を
構成するマイクロプロセッサにこの演算を行なわせるこ
とは難しく、また演算時間も長くなり、素早い診断が要
求される超音波診断装置としては不都合が生じる結果と
なる。
本発明は以上のような問題に鑑みてなされたもので、超
音波画像上の円弧状部における３点の座標を用いて曲率
半径を求めることによシ、簡単な演算で、高速で、しか
も容易な操作で曲率半径を得ることのできる超音波診断
装置を提供することを目的とするものである。
音波画像上の円弧状部における３点の座標を用いて曲率
半径を求めることによシ、簡単な演算で、高速で、しか
も容易な操作で曲率半径を得ることのできる超音波診断
装置を提供することを目的とするものである。
以下に本発明の一実施例を説明する。
本発明の第１の実施例は、超音波画像の円弧状部上の任
意の３点から曲率半径を求めることである０ ｘ−ｙ平面内の円の方程式はｇ、ｆ、ａを定数として ｘ２＋ ７２＋ ２ｑｘ＋２ｆｙ＋ ｃ ＝Ｏ−−−−
・−（２）で与えられるので、同一円弧上に任意の３点
（ｘ、＋ ７１ ）＋ （”２１１２）＋ （”３＋
７３）を設定すれば、これら３組のＸ、７座標の値を第
（２）式に代入することによって未知定数ｇ、ｆ、ｃを
求めるととが出来、従ってｒ＝↓ｒ２ ’＋ １２
ｃ で与えられる曲率半径を算出することは可能であ
る。第（２）式に（”１１７１ ）＋ （”２１ ｙ２
）ｙ （ｘ３＋ ３’３）を代入すると ｘ２１＋ｙ２１＋２ｑｘ１＋２１ｙ、＋ｃ＝＝Ｏ・・−
・−（３）ｘ’、＋ｙ’、＋２ｑｘ２＋２１ｙ２＋ｃ＝
ｏ −・−−−−＝−（４）Ｘ’ ＋ｙ２＋２ｑｘ
＋２ｆｙ３＋ｃ＝Ｏ・−・−・−（ｔｓ）３ ３
３ となり、第（３）、 （４）、 （５）式よりここに
・・・（８）従って
、第（６）、 （力、（８）式とｒ ＝ ｇ２＋ ｆ２
− ｃ −−−−−−−−−（１０）−上り
、３点の座標（Ｘｌ、 ３’１ ）ｙ （”２Ｔ Ｖ２
）ｔ（ｘ：３＋ ｙ ３ ）を用いて曲率半径ｒを求め
ることができる。
意の３点から曲率半径を求めることである０ ｘ−ｙ平面内の円の方程式はｇ、ｆ、ａを定数として ｘ２＋ ７２＋ ２ｑｘ＋２ｆｙ＋ ｃ ＝Ｏ−−−−
・−（２）で与えられるので、同一円弧上に任意の３点
（ｘ、＋ ７１ ）＋ （”２１１２）＋ （”３＋
７３）を設定すれば、これら３組のＸ、７座標の値を第
（２）式に代入することによって未知定数ｇ、ｆ、ｃを
求めるととが出来、従ってｒ＝↓ｒ２ ’＋ １２
ｃ で与えられる曲率半径を算出することは可能であ
る。第（２）式に（”１１７１ ）＋ （”２１ ｙ２
）ｙ （ｘ３＋ ３’３）を代入すると ｘ２１＋ｙ２１＋２ｑｘ１＋２１ｙ、＋ｃ＝＝Ｏ・・−
・−（３）ｘ’、＋ｙ’、＋２ｑｘ２＋２１ｙ２＋ｃ＝
ｏ −・−−−−＝−（４）Ｘ’ ＋ｙ２＋２ｑｘ
＋２ｆｙ３＋ｃ＝Ｏ・−・−・−（ｔｓ）３ ３
３ となり、第（３）、 （４）、 （５）式よりここに
・・・（８）従って
、第（６）、 （力、（８）式とｒ ＝ ｇ２＋ ｆ２
− ｃ −−−−−−−−−（１０）−上り
、３点の座標（Ｘｌ、 ３’１ ）ｙ （”２Ｔ Ｖ２
）ｔ（ｘ：３＋ ｙ ３ ）を用いて曲率半径ｒを求め
ることができる。
次に本発明の第２の実施例を説明する。第２の実施例は
、任意の２点とその中点とから曲率半径を求めるもので
ある。
、任意の２点とその中点とから曲率半径を求めるもので
ある。
説明を簡単にするため、第３図においてＸ −７平面上
の円弧Ｋ）の曲率半径を求める方法を説明する。ここで
、超音波診断装置における表示画面をＰとし、この平面
をＸ−３Ｆ平面としたとき、画面Ｐの左下隅を原点０（
ｏ、ｏ）としても−膜性を失わないので、第３図のよう
に座標軸を設定するものとする。
の円弧Ｋ）の曲率半径を求める方法を説明する。ここで
、超音波診断装置における表示画面をＰとし、この平面
をＸ−３Ｆ平面としたとき、画面Ｐの左下隅を原点０（
ｏ、ｏ）としても−膜性を失わないので、第３図のよう
に座標軸を設定するものとする。
第２図に示したボリューム１１．１２を用いて第３図に
おける矢印Ａ１で示すように任意の位置からマーカー３
１をまずＡ点の位置に移動させ、マーカー３１がＡ点と
一致したときスイッチ１９を用いてＡ／Ｄコンバータ１
３の変換動作を中断させ、Ａ点のＸ、７座標の値（”１
＋ｙ１） をデータ記憶用ＲＡＭ２０に格納する。次
に矢印Ａ２で示すように同様にしてマーカー３１をＡ点
からＢ点に移動させ、Ｂ点に一致したときのｘ、ｙ座標
の値（ｘ２． ｙ２） をＲＡＭ２０に格納する。
おける矢印Ａ１で示すように任意の位置からマーカー３
１をまずＡ点の位置に移動させ、マーカー３１がＡ点と
一致したときスイッチ１９を用いてＡ／Ｄコンバータ１
３の変換動作を中断させ、Ａ点のＸ、７座標の値（”１
＋ｙ１） をデータ記憶用ＲＡＭ２０に格納する。次
に矢印Ａ２で示すように同様にしてマーカー３１をＡ点
からＢ点に移動させ、Ｂ点に一致したときのｘ、ｙ座標
の値（ｘ２． ｙ２） をＲＡＭ２０に格納する。
ここでＸ、とＸ２．”１とｙ２をそれぞれマイクロプロ
セッサ１６によって比較し、ｙ１＝ｙ２のときは手順Ｃ
Ｉ）へ、ｙ、＼ｙ２のときは手順〔ｌ）へ進むようにマ
イクロプロセッサ１５で制御する。当然”１””２．Ｖ
、”Ｖ２が同時に成シ立つことはあり得ない。
セッサ１６によって比較し、ｙ１＝ｙ２のときは手順Ｃ
Ｉ）へ、ｙ、＼ｙ２のときは手順〔ｌ）へ進むようにマ
イクロプロセッサ１５で制御する。当然”１””２．Ｖ
、”Ｖ２が同時に成シ立つことはあり得ない。
〔手順１ ］ ｙ、＝ｙ２のとき
ｙ１＝ｙ２．となる場合の一例を第４図に示す。矢印Ａ
３で示す方向に円弧ＡＢに添ってマーカー３１をＢ点か
ら中点Ｍへ移動させる。中点ＭはＸ、 Ｘ座標を（ｘ３
．ｙ３） とすると ｘ３−（ｘ１＋ｘ２）／２ ・・・・・・・・・
（１１）なる関係が成立・する点であり、マーカー３１
を移動させながら各点におけるＸ座標の値と（ｘ１＋Ｘ
２ ） ／２の値とをマイクロプロセッサ１５によっ
てリアルタイムで判定を行ない、第（１１）式が成立す
る点を中点Ｍとして設定し、同時にＸ座標の値ｘ３およ
びその時のｙ座標の値ｙ３をＲＡＭ２０に格納する。Ｍ
点ではｙ座標の値ｙ３は、当然ｙ４．ｙ２と等しくない
。もしｙ３＝ｙ１＝ｙ２であればＡ、Ｂ、Ｍ点は一直線
上に存在することになり円弧ＡＢの曲率半径の無限大と
なる０また明らかにｘｌとｘ２は等しくない。（第（３
）式−第（４）式）にｙ １＝ｙ２なる条件および中点
条件第（１１）式を代入すると、 （ｘ １！ ２ ） （Ｘ ３＋ ｇ ） −〇！、＼
ｘ２ゆえ ｇ＝−！３ ・・・・・・（１２）
第（１２）式を第（５）式に代入するとｃ＝ノーｙ２−
２ｆｙ３ ・・・・・・ （１３） ３ 第（３）式に第（１２）？ （１３）式を代入して未知
定数量を求めると、 第（１２）、（１３）、（１４）式を第（１ｏ）式に代
入して曲率半径ｒを求めると、 となる。従ってＲＡＭ２０に格納した２点の座標Ａ（ｘ
１． ｙ １） およびＭ （Ｋ３． Ｋ３）を用
いて第（１５）式に従って、マイクロプロセッサ１６で
計算を行なえば曲率半径ｒを算出することができる。
３で示す方向に円弧ＡＢに添ってマーカー３１をＢ点か
ら中点Ｍへ移動させる。中点ＭはＸ、 Ｘ座標を（ｘ３
．ｙ３） とすると ｘ３−（ｘ１＋ｘ２）／２ ・・・・・・・・・
（１１）なる関係が成立・する点であり、マーカー３１
を移動させながら各点におけるＸ座標の値と（ｘ１＋Ｘ
２ ） ／２の値とをマイクロプロセッサ１５によっ
てリアルタイムで判定を行ない、第（１１）式が成立す
る点を中点Ｍとして設定し、同時にＸ座標の値ｘ３およ
びその時のｙ座標の値ｙ３をＲＡＭ２０に格納する。Ｍ
点ではｙ座標の値ｙ３は、当然ｙ４．ｙ２と等しくない
。もしｙ３＝ｙ１＝ｙ２であればＡ、Ｂ、Ｍ点は一直線
上に存在することになり円弧ＡＢの曲率半径の無限大と
なる０また明らかにｘｌとｘ２は等しくない。（第（３
）式−第（４）式）にｙ １＝ｙ２なる条件および中点
条件第（１１）式を代入すると、 （ｘ １！ ２ ） （Ｘ ３＋ ｇ ） −〇！、＼
ｘ２ゆえ ｇ＝−！３ ・・・・・・（１２）
第（１２）式を第（５）式に代入するとｃ＝ノーｙ２−
２ｆｙ３ ・・・・・・ （１３） ３ 第（３）式に第（１２）？ （１３）式を代入して未知
定数量を求めると、 第（１２）、（１３）、（１４）式を第（１ｏ）式に代
入して曲率半径ｒを求めると、 となる。従ってＲＡＭ２０に格納した２点の座標Ａ（ｘ
１． ｙ １） およびＭ （Ｋ３． Ｋ３）を用
いて第（１５）式に従って、マイクロプロセッサ１６で
計算を行なえば曲率半径ｒを算出することができる。
〔手順）ｙ１＼ｙ２のとき
このとき、まずｘにｘ２が成立するかどうかをマイクロ
プロセッサ１６によって判定し、ｘ１＝ｘ２の時は〔手
順■−１〕に、ｘ１＼ｘ２のときは〔手順１−２〕に進
むＯ 〔手順ｙ、 −１） ｘ１＝ｘ２のときＸ １” ！
２となる場合の一例を第５図に示す。この時も〔手順１
．〕と同様に矢印Ａ３で示す方向に円弧ＡＢに添ってマ
ーカー３１をＢ点から中点Ｍへ移動させる。ここで中点
ＭはＸ、７座標を（Ｋ３．Ｋ３） とすると、 ｙ３＝（ｙ１＋ｙ２）／２ ・・・・・・（１
６）の値とをマイクロプロセッサ１６によってリアルタ
イムで判定を行ない、第（１６）式が成立する点をＭ点
として設定し、同時にＸ座標の値ｙ３およびその時の１
座標の値ｘ３をＲＡＭ２０に格納する。（第（３）式−
第（４）式）にｘ、−Ｋ２および中点条件第（１６）式
を代入するとｙ１＼ｙ２ゆえ□、７）えオ、シ）えよい
よ、よ −（°″ゝ（＝−ｙ ａ ＝＝ Ｋ２−Ｋ２−２ｇｘ３ ・・・・・・
・・・（１８） ３ 第（３）式に第（１７）、（１８）式を代入して未知定
数９を求めると 第（１７）、（１Ｂ）、（１９）式を第（１０）式に代
入して曲率半径ｒを求めると、 ”となる。従ってＲＡＭ２０に格納した２点の座標Ａ（
ｘｌ；ｙｌ）およびＭ（ｘ３＋ Ｋ３） を用いて第
（２ｏ）式に従ってマイクロプロセッサ１５で計算を行
なえば曲率半径ｒを算出することができる。
プロセッサ１６によって判定し、ｘ１＝ｘ２の時は〔手
順■−１〕に、ｘ１＼ｘ２のときは〔手順１−２〕に進
むＯ 〔手順ｙ、 −１） ｘ１＝ｘ２のときＸ １” ！
２となる場合の一例を第５図に示す。この時も〔手順１
．〕と同様に矢印Ａ３で示す方向に円弧ＡＢに添ってマ
ーカー３１をＢ点から中点Ｍへ移動させる。ここで中点
ＭはＸ、７座標を（Ｋ３．Ｋ３） とすると、 ｙ３＝（ｙ１＋ｙ２）／２ ・・・・・・（１
６）の値とをマイクロプロセッサ１６によってリアルタ
イムで判定を行ない、第（１６）式が成立する点をＭ点
として設定し、同時にＸ座標の値ｙ３およびその時の１
座標の値ｘ３をＲＡＭ２０に格納する。（第（３）式−
第（４）式）にｘ、−Ｋ２および中点条件第（１６）式
を代入するとｙ１＼ｙ２ゆえ□、７）えオ、シ）えよい
よ、よ −（°″ゝ（＝−ｙ ａ ＝＝ Ｋ２−Ｋ２−２ｇｘ３ ・・・・・・
・・・（１８） ３ 第（３）式に第（１７）、（１８）式を代入して未知定
数９を求めると 第（１７）、（１Ｂ）、（１９）式を第（１０）式に代
入して曲率半径ｒを求めると、 ”となる。従ってＲＡＭ２０に格納した２点の座標Ａ（
ｘｌ；ｙｌ）およびＭ（ｘ３＋ Ｋ３） を用いて第
（２ｏ）式に従ってマイクロプロセッサ１５で計算を行
なえば曲率半径ｒを算出することができる。
〔手順ＩＩ ２ ） ｘ １＼ｘ２のときこの場合の
一例を第６図に示す。Ａ点、Ｂ点を結ぶ直線とＸ軸のな
す角をθとすると、θは次式１式％ 従って角θを求めるにはＲＡＭ２０に格納した２点の座
標Ａ （Ｋ１． ｙｌ ）ｔ Ｂ （Ｋ２． Ｋ２）お
よびＲＯＭ２１に格納された逆三角関数変換テーブルか
ら第（２１）式、に従ってマイクロプロセッサ１５で計
算することができる。次に、 Ｘ ＝＝ ＸＣＯ，θ＋ｙｌｌｉｎθ ） ・・・・・・・・・（２２）Ｙ ＝ −ｘ
ｓ ｉｎθ＋Ｙ ｃｏｓなる座標変換により第６図に
示すようにもとのＸ −７座標に対して角θだけ回転し
たＸ−Ｘ座標を設定すれば第４図に示したｙｌ−＝ｙ２
の場合と同様に扱うことができる。Ａ、Ｂ点の座標軸回
転後ノ新座４ｉ１ソレソｉ （Ｘｌ、 ｙｌ）、 （Ｋ
２． Ｋ２） とすると、ＲＡＭ２０に格納した２点
の座標（ｘｌ。
一例を第６図に示す。Ａ点、Ｂ点を結ぶ直線とＸ軸のな
す角をθとすると、θは次式１式％ 従って角θを求めるにはＲＡＭ２０に格納した２点の座
標Ａ （Ｋ１． ｙｌ ）ｔ Ｂ （Ｋ２． Ｋ２）お
よびＲＯＭ２１に格納された逆三角関数変換テーブルか
ら第（２１）式、に従ってマイクロプロセッサ１５で計
算することができる。次に、 Ｘ ＝＝ ＸＣＯ，θ＋ｙｌｌｉｎθ ） ・・・・・・・・・（２２）Ｙ ＝ −ｘ
ｓ ｉｎθ＋Ｙ ｃｏｓなる座標変換により第６図に
示すようにもとのＸ −７座標に対して角θだけ回転し
たＸ−Ｘ座標を設定すれば第４図に示したｙｌ−＝ｙ２
の場合と同様に扱うことができる。Ａ、Ｂ点の座標軸回
転後ノ新座４ｉ１ソレソｉ （Ｘｌ、 ｙｌ）、 （Ｋ
２． Ｋ２） とすると、ＲＡＭ２０に格納した２点
の座標（ｘｌ。
ｙｌ）＋ （”２．７２）およびＲＯＭ２１に格納した
三角関数変換テーブルより として新座標の値を求めることができる。また中点Ｍの
新座標を（Ｋ３．Ｋ３）とするとＭ点はＫ３−（Ｘ、＋
Ｋ２）／２ ・・・・・・・・・（２３）が成立
する点であり、マーカー３１を円弧ＡＢに添って矢印Ａ
３の方向へ移動させながら各点におけるＸ座標の値と（
Ｘ １＋ Ｘ ２ ） ／２の値とをマイクロプロセッ
サ１５によってリアルタイムで判定を行ない第（２４）
式が成立する点を中点Ｍとして設定し、同時にＸ座標の
値ｘ３およびＸ座標の値Ｙ をＲＡＭ２０に格納する。
三角関数変換テーブルより として新座標の値を求めることができる。また中点Ｍの
新座標を（Ｋ３．Ｋ３）とするとＭ点はＫ３−（Ｘ、＋
Ｋ２）／２ ・・・・・・・・・（２３）が成立
する点であり、マーカー３１を円弧ＡＢに添って矢印Ａ
３の方向へ移動させながら各点におけるＸ座標の値と（
Ｘ １＋ Ｘ ２ ） ／２の値とをマイクロプロセッ
サ１５によってリアルタイムで判定を行ない第（２４）
式が成立する点を中点Ｍとして設定し、同時にＸ座標の
値ｘ３およびＸ座標の値Ｙ をＲＡＭ２０に格納する。
この時のＭ点の新座標（Ｋ３．Ｋ３）は座標軸回転前の
座標（Ｋ３゜ｙ３）を用いて第（２３）式と同様に として求めることができる。以上のようにして新座標を
求めれば第（１６）式を用いて曲率半径Ｔはとして求め
ることができる。
座標（Ｋ３゜ｙ３）を用いて第（２３）式と同様に として求めることができる。以上のようにして新座標を
求めれば第（１６）式を用いて曲率半径Ｔはとして求め
ることができる。
以上の、丁−イクロプロセッサ１６で行なう手順を第７
図にフローチャートとして示す。
図にフローチャートとして示す。
一般にはｘ１＼Ｘ ２ 、かつｙ１＼ｙ２の場合が最も
多く、従って〔手順Ｉｔ−，２）の手順を用いることが
ほとんどである。また実際に使用する場合はマーカーを
Ａ点に設定したときにＡ点の座標（ｘｌ。
多く、従って〔手順Ｉｔ−，２）の手順を用いることが
ほとんどである。また実際に使用する場合はマーカーを
Ａ点に設定したときにＡ点の座標（ｘｌ。
ｙｌ）を（０，Ｏ）と初期設定することが多いので削算
はより簡単な形になる。本実施例ではまずｙ１＝ｙ２の
判定を行なったが、第８図に示すようにｘ、＝ｘ２の判
定から始めてもよい。
はより簡単な形になる。本実施例ではまずｙ１＝ｙ２の
判定を行なったが、第８図に示すようにｘ、＝ｘ２の判
定から始めてもよい。
以上述べたように本発明によれば、超音波画像の任意の
円弧ＡＢの曲率半径を、まず二点Ａ、 Ｂを設定し、次
に円弧Ａ″Ｂ上をマーカーを移動させるという簡単な操
作によシ、高速で円弧ＡＢの曲率半径を容易な演算で求
めることができ、頚動脈の曲率半径をも容易に定量的に
求めることができるなど診断に有用な情報を得ることが
できる。
円弧ＡＢの曲率半径を、まず二点Ａ、 Ｂを設定し、次
に円弧Ａ″Ｂ上をマーカーを移動させるという簡単な操
作によシ、高速で円弧ＡＢの曲率半径を容易な演算で求
めることができ、頚動脈の曲率半径をも容易に定量的に
求めることができるなど診断に有用な情報を得ることが
できる。
なお、本発明では、Ａ点の座標を初期に（ｏ。
０）と設定する必要はない。さらに中点Ｍを求める場合
、マーカーは中点近傍だけ円弧上に添って移動させるだ
けでよく、必ずしもＢ点から円弧上を移動させる必要は
ない。
、マーカーは中点近傍だけ円弧上に添って移動させるだ
けでよく、必ずしもＢ点から円弧上を移動させる必要は
ない。
第１図は超音波診断装置のブロック図、第２図は電子キ
ャリパ−回路のブロック図、第３図は超音波診断装置の
表示画像上に表示された円弧を説明するための図、第４
図〜第６図は本発明の一実施例を説明するための図、第
７図は同実施例の手順を示すフローチャート、第８図は
本発明の他の実施例の手順を示すフローチャートである
。 ＡＩ、Ａ２．Ａ３・・・・・・マーカーの移動方向、Ａ
Ｂ・・・・・・円弧、Ｍ・・・・・・中点、Ｐ・・・・
・・表示画面、１１．１２・・・・・・ボリューム、１
３・・・・・・タイミング切換回路、１４・°・°・・
Ａ／Ｄコン・（−タ、１６・・・；・・マイクロプロセ
ッサ、１６・・・・・・ＣＲＴコントローラ、１７・・
・・・・キャラクタジェネレータ、１８・・・・・・表
示部、１９・・・・・・スイッチ、２０・・・・・・Ｒ
ＡＭ、２１・・・・・・ＲＯＭ、３１・・・・・・マー
カー。 代理人の氏名 弁理士 中 尾 敏 男 ほか１名第３
図 第４図 第６図 第７図
ャリパ−回路のブロック図、第３図は超音波診断装置の
表示画像上に表示された円弧を説明するための図、第４
図〜第６図は本発明の一実施例を説明するための図、第
７図は同実施例の手順を示すフローチャート、第８図は
本発明の他の実施例の手順を示すフローチャートである
。 ＡＩ、Ａ２．Ａ３・・・・・・マーカーの移動方向、Ａ
Ｂ・・・・・・円弧、Ｍ・・・・・・中点、Ｐ・・・・
・・表示画面、１１．１２・・・・・・ボリューム、１
３・・・・・・タイミング切換回路、１４・°・°・・
Ａ／Ｄコン・（−タ、１６・・・；・・マイクロプロセ
ッサ、１６・・・・・・ＣＲＴコントローラ、１７・・
・・・・キャラクタジェネレータ、１８・・・・・・表
示部、１９・・・・・・スイッチ、２０・・・・・・Ｒ
ＡＭ、２１・・・・・・ＲＯＭ、３１・・・・・・マー
カー。 代理人の氏名 弁理士 中 尾 敏 男 ほか１名第３
図 第４図 第６図 第７図
Claims (1)
- 【特許請求の範囲】 （１）被検体の超音波画像を表示する表示部と、前記表
示部に表示された前記超音波画像の円弧状部」二の３点
を指示する指示手段と、指示された前記３点の前記表示
部上における座標値より前記円弧状部の曲率半径を求め
る演算手段とを備えたことを特徴とする超音波診断装置
。 （２）指示手段が、電子キャリバーであることを特徴と
する特許請求の範囲第１項記載の超音波診断装置。 （３）円弧状部上の３点が、それぞれ任意の点であるこ
とを特徴とする特許請求の範囲第１項記載の超音波診断
装置。 （４）円弧状部上の３点が、任意の２点と、前記２点に
挾まれた前記円弧状部の中点であることを特徴とする特
許請求の範囲第１項記載の超音波診断装置。 （６）演算手段が、任意の点のそれぞれの′座標値を読
み取る手順と、前記任意の２点を結ぶ直線の傾きを求め
る手順と、前記直線の傾きに応じて前記任意の２点の座
標値を変換する手順と、前記円弧状部の中点を求める手
順と、前記任意の２点および前記中点の座標値から前記
円弧状部の曲率半径を計算する手順とを有することを特
徴とする特許請求の範囲第４項記載の超音波診断装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP10791982A JPS5947A (ja) | 1982-06-23 | 1982-06-23 | 超音波診断装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP10791982A JPS5947A (ja) | 1982-06-23 | 1982-06-23 | 超音波診断装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPS5947A true JPS5947A (ja) | 1984-01-05 |
Family
ID=14471367
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP10791982A Pending JPS5947A (ja) | 1982-06-23 | 1982-06-23 | 超音波診断装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS5947A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61155755A (ja) * | 1984-12-27 | 1986-07-15 | Toshiba Corp | 超音波探傷装置 |
| US4904460A (en) * | 1988-04-13 | 1990-02-27 | Mitsubishi Kasei Corporation | Process for producing monosilane |
-
1982
- 1982-06-23 JP JP10791982A patent/JPS5947A/ja active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61155755A (ja) * | 1984-12-27 | 1986-07-15 | Toshiba Corp | 超音波探傷装置 |
| US4904460A (en) * | 1988-04-13 | 1990-02-27 | Mitsubishi Kasei Corporation | Process for producing monosilane |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPS6290144A (ja) | 超音波診断装置 | |
| JPH0428354A (ja) | 超音波診断装置 | |
| JPS5947A (ja) | 超音波診断装置 | |
| US20030199766A1 (en) | Ultrasonic image display apparatus and ultrasonic image display method | |
| JP3267739B2 (ja) | 超音波カラードプラ診断装置 | |
| JP2002282250A (ja) | 超音波画像診断装置 | |
| JP2002330966A (ja) | 超音波診断装置 | |
| JPS636014B2 (ja) | ||
| JP3532129B2 (ja) | 超音波診断装置 | |
| JPH04295346A (ja) | 超音波断層装置 | |
| JPH01230346A (ja) | 超音波診断装置 | |
| JPH09322893A (ja) | 超音波診断装置 | |
| JPH02193655A (ja) | 超音波診断装置 | |
| JPH0417548B2 (ja) | ||
| JPH0549640A (ja) | 超音波診断装置 | |
| JPS63186634A (ja) | 超音波診断装置 | |
| JP2874932B2 (ja) | 超音波診断装置 | |
| JPH0677584B2 (ja) | 超音波診断装置 | |
| JPH069614Y2 (ja) | 超音波診断装置 | |
| JPS62139631A (ja) | 超音波断層装置 | |
| JPH05103784A (ja) | 超音波診断装置 | |
| JPH0471543A (ja) | 断層診断装置 | |
| JPH066123B2 (ja) | 超音波ドプラ診断装置 | |
| CN111557687A (zh) | 超声波诊断装置、记录介质及控制台的向导显示方法 | |
| JPS632613B2 (ja) |