JP2015037491A - 超音波診断装置、超音波画像解析方法、およびプログラム - Google Patents
超音波診断装置、超音波画像解析方法、およびプログラム Download PDFInfo
- Publication number
- JP2015037491A JP2015037491A JP2013169527A JP2013169527A JP2015037491A JP 2015037491 A JP2015037491 A JP 2015037491A JP 2013169527 A JP2013169527 A JP 2013169527A JP 2013169527 A JP2013169527 A JP 2013169527A JP 2015037491 A JP2015037491 A JP 2015037491A
- Authority
- JP
- Japan
- Prior art keywords
- joint
- bone surface
- ultrasonic
- cavity region
- ultrasonic diagnostic
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 40
- 210000000988 bone and bone Anatomy 0.000 claims abstract description 112
- 230000009266 disease activity Effects 0.000 claims abstract description 17
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 238000005259 measurement Methods 0.000 claims abstract description 8
- 201000010099 disease Diseases 0.000 claims abstract 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims abstract 2
- 210000000281 joint capsule Anatomy 0.000 claims description 34
- 238000002604 ultrasonography Methods 0.000 claims description 21
- 210000001519 tissue Anatomy 0.000 claims description 18
- 210000002435 tendon Anatomy 0.000 claims description 12
- 230000001575 pathological effect Effects 0.000 claims description 8
- 239000000523 sample Substances 0.000 claims description 7
- 238000003703 image analysis method Methods 0.000 claims description 3
- 238000010801 machine learning Methods 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 2
- JTJMJGYZQZDUJJ-UHFFFAOYSA-N phencyclidine Chemical compound C1CCCCN1C1(C=2C=CC=CC=2)CCCCC1 JTJMJGYZQZDUJJ-UHFFFAOYSA-N 0.000 claims 1
- 238000003745 diagnosis Methods 0.000 abstract description 4
- 230000017531 blood circulation Effects 0.000 description 11
- 238000012545 processing Methods 0.000 description 9
- 210000001258 synovial membrane Anatomy 0.000 description 9
- 210000000845 cartilage Anatomy 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 7
- 238000001514 detection method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 206010039073 rheumatoid arthritis Diseases 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000011002 quantification Methods 0.000 description 3
- 206010061218 Inflammation Diseases 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000003708 edge detection Methods 0.000 description 2
- 230000004054 inflammatory process Effects 0.000 description 2
- 238000011158 quantitative evaluation Methods 0.000 description 2
- 238000007637 random forest analysis Methods 0.000 description 2
- 210000003491 skin Anatomy 0.000 description 2
- 230000008719 thickening Effects 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 206010051728 Bone erosion Diseases 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 230000002491 angiogenic effect Effects 0.000 description 1
- 206010003246 arthritis Diseases 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 210000001145 finger joint Anatomy 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 208000011580 syndromic disease Diseases 0.000 description 1
- 210000001179 synovial fluid Anatomy 0.000 description 1
Images
Abstract
【課題】超音波診断装置を用いたリウマチ診断において、より高精度に骨表面と関節腔領域を検出可能な超音波診断装置を提供する。【解決手段】測定対象から得られた超音波画像から骨表面を特定する骨表面特定部と、前記骨表面に存在する関節位置を特定する関節位置特定部と、前記関節位置より浅い位置にある関節腔領域を特定する関節腔領域特定部と、前記関節腔領域の内部に着目し、前記測定対象の疾患活動性を評価する病態解析部と、を備える。【選択図】図１
Description
本発明は、超音波の送受信を行う超音波プローブを被検者の体表面、特に四肢関節に当接させて超音波診断を行う超音波診断装置、超音波画像解析方法及びそのプログラムに関するものである。
近年関節リウマチをはじめとする関節炎の疾患活動性の評価に超音波診断装置を用いることが一般的になりつつある。疾患活動性の評価にはＢモード画像とパワードプラ画像が主に用いられており、Ｂモード画像では滑膜肥厚や滑液貯留、骨びらんを、パワードプラ画像では滑膜の炎症を観察することができる。
加えて、超音波画像を用いてこれらの疾患活動性をグレード別に判定する方法が提案されている。パワードプラ画像を用いて炎症の程度をグレード化する場合、観測される血流信号が肥厚した滑膜領域のうちどれだけの割有を占めているかでグレードを判定する。その際検査者の主観で判断しているため、検査者間でグレードのばらつきが生じる。
上記課題を解決するために、疾患活動性を定量化する方法が提案されている。例えば非特許文献１においては、肥厚した滑膜を含む関節腔を超音波画像上にてフリーハンドでトレースし、このトレースした領域に占める血流信号の占有率を定量評価値として算出することが提案されている。
しかしながら、このような方法においては、関節腔を正確にトレースすることを検査者に要求しなくてはならない。診断フローにフリーハンドトレースのステップを追加することは好ましくなく、またトレース結果が検査者に依存するため、血流信号の占有率が検査者間でばらつくことが想定される。
上記課題を解決するために、例えば特許文献１では骨表面を特定してから関節腔領域を切り出し、関節腔領域を解析することにより疾患活動性を客観的に定量化する方法が提案されている。
小池隆夫、「超音波検査法を用いた関節リウマチの新しい診療」、Ｐ４０−４３、メディカルレビュー社、２０１０年３月１０日
疾患活動性をより客観的に定量評価するためには、関節腔領域の検出精度を向上させることが重要である。同時に、関節腔領域を検出するために必要な骨表面を特定する精度の向上も必要である。
本発明の目的は、より高精度に骨表面と関節腔領域を検出可能な超音波診断装置を提供することである。
上記課題を解決するために、本発明の一態様に係る超音波診断装置は、前記超音波画像から骨表面を特定する骨表面特定部と、前記骨表面に存在する関節位置を特定する関節位置特定部と、前記関節位置より浅い位置にある関節腔領域を特定する関節腔領域特定部と、前記関節腔領域の内部に着目し、前記測定対象の疾患活動性を評価する病態解析部と、を備える。
骨表面や関節包が他の組織と比較して明確な輝度差を持たない等様々なパターンの関節超音波画像に対して関節腔を検出する場合、本発明によれば画像のパターンに依存しにくい高精度な検出が可能となる。
以下、本発明の実施の形態について、図面を参照しながら説明する。
（実施の形態１）
図１は、本発明の超音波診断装置の構成を示すブロック図である。図１において、本発明の超音波診断装置は、超音波送受信部１００２、超音波画像生成部１００３、記憶部１００４、外部入力取得部１００５、関節腔領域特定部１００８、病態解析部１００９、画面作成部１０１０、から構成される。ここで、超音波プローブ１００１、外部入力部１０１１、表示部１０１２は超音波診断装置外の構成としているが、必要に応じて超音波診断装置内の構成としても構わない。超音波プローブ１００１は、超音波送受信部１００２で生成された送信波を放出する。放出された超音波は生体内の組織間の音響インピーダンスの差が異なる部分において反射される。この時、音響インピーダンスの差が顕著であるほど、反射する超音波のエネルギーも大きくなる。
図１は、本発明の超音波診断装置の構成を示すブロック図である。図１において、本発明の超音波診断装置は、超音波送受信部１００２、超音波画像生成部１００３、記憶部１００４、外部入力取得部１００５、関節腔領域特定部１００８、病態解析部１００９、画面作成部１０１０、から構成される。ここで、超音波プローブ１００１、外部入力部１０１１、表示部１０１２は超音波診断装置外の構成としているが、必要に応じて超音波診断装置内の構成としても構わない。超音波プローブ１００１は、超音波送受信部１００２で生成された送信波を放出する。放出された超音波は生体内の組織間の音響インピーダンスの差が異なる部分において反射される。この時、音響インピーダンスの差が顕著であるほど、反射する超音波のエネルギーも大きくなる。
反射された超音波は、超音波プローブ１００１で受信される。超音波プローブ１００１で受信された超音波は超音波受信信号として超音波送受信部１００２に入力される。超音波受信信号は超音波送受信部１００２でビームフォーミング、検波、対数圧縮などの処理を施された後に、音響線信号として超音波画像生成部１００３に入力される。
超音波画像生成部１００３は、超音波送受信部１００２から入力された音響線信号をまとめて、超音波画像を生成する。超音波画像生成部１００３で生成された超音波画像は記憶部１００４で一時的に記憶される。
外部入力部１０１１は、検査者が検査者名、患者名、超音波診断装置の設定情報を入力するための構成である。
外部入力取得部１００５は、検査者が外部入力部１０１１を通して入力した設定情報と、記憶部１００４に記憶されている超音波画像を対応付ける。
骨表面特定部１００６は、記憶部１００４に記憶されている超音波画像内の骨表面を特定する。骨表面が特定された超音波画像は、記憶部１００４に記憶される。
関節位置特定部１００７は、記憶部１００４に記憶されている骨表面が特定された超音波画像内の関節位置を特定する。関節位置が特定された超音波画像は、記憶部１００４に記憶される。
関節腔領域特定部１００８は、記憶部１００４に記憶されている超音波画像内の関節腔領域を特定する。関節腔領域が特定された超音波画像は、記憶部１００４に記憶される。
病態解析部１００９は、記憶部１００４に記憶されている関節腔領域が特定された超音波画像の関節腔領域を解析し、病態の疾患活動性を定量的に評価する。
表示部１０１２は、超音波画像を検査者に提示する。その際には記憶部１００４から読み出した画像に対して、画面作成部１０１０によって検査者名、患者名、時間情報、超音波診断装置の設定情報、疾患活動性の評価結果などが重畳された超音波画像を提示する。
本実施の形態における主な測定対象は関節である。関節の一例として手指関節の概略を図２に示す。図２に示すように関節は、骨２００１、骨２００２、軟骨２００３、軟骨２００４、滑膜２００５、関節包２００６からなる部位である。骨２００１と骨２００２の先端部分には軟骨２００３と２００４が付帯しており、軟骨２００３と軟骨２００４を包むように滑膜２００５が存在する。さらに滑膜２００５を囲むように関節包２００６が骨２００１と骨２００２に付着している。関節リウマチが進行すると滑膜２００５の肥厚が確認される。また、多くの症例では、滑膜２００５内部での血管新生の増殖も確認される。
関節を測定対象として超音波画像が撮像されると、図３のような画像が取得される。超音波画像には、骨表面３００１、骨表面３００２、皮膚３００３、腱３００４、関節包３００５が描画される。骨表面３００１、骨表面３００２や皮膚３００３は比較的固い組織であるため、超音波画像上でも高輝度で描画される。超音波の大部分は骨表面で反射されるため、骨の内部は描画されず、骨表面のみが高輝度で描画される。腱３００４や関節包３００５は骨表面３００１、骨表面３００２、皮膚３００３と比較すると低輝度で描画される。また滑膜や軟骨部分はほぼ輝度値を持たない。従って、関節超音波画像で比較的高い輝度で描画される組織は、皮膚、腱、関節包、骨表面である。
超音波診断装置の関節腔領域特定部１００８では、図３に代表される関節超音波画像から関節包３００５と骨表面３００１、骨表面３００２に囲まれた関節腔領域３００６が抽出される。関節腔領域特定部１００８における処理を図４を用いて説明する。
［ステップＳ４０１］
ステップＳ４０１で行われる超音波画像から骨表面を特定する処理の詳細を、図５（Ａ）、図５（Ｂ）を用いて説明する。図５（Ａ）の５００１と５００２はステップＳ４０１で最終的に特定される骨表面である。図５（Ｂ）は超音波画像内のエッジから骨表面を特定する方法を説明する図であり、超音波画像内のエッジが皮膚に属するエッジ５００３、腱に属するエッジ５００４、骨に属するエッジ５００５のいずれかに分類されることを示している。
ステップＳ４０１で行われる超音波画像から骨表面を特定する処理の詳細を、図５（Ａ）、図５（Ｂ）を用いて説明する。図５（Ａ）の５００１と５００２はステップＳ４０１で最終的に特定される骨表面である。図５（Ｂ）は超音波画像内のエッジから骨表面を特定する方法を説明する図であり、超音波画像内のエッジが皮膚に属するエッジ５００３、腱に属するエッジ５００４、骨に属するエッジ５００５のいずれかに分類されることを示している。
骨表面を特定するためには、超音波画像からエッジを抽出し、抽出したエッジを各組織に分類する方法がある。超音波画像からエッジを検出する方法としては、事前に設定した骨表面の輝度値の代表値、あるいは検査者が外部入力部１０１１を介して入力する値を閾値として二値化処理を行う方法、Ｓｏｂｅｌオペレータ、ラプラシアンフィルタ等を用いて微分画像から強勢なエッジを取り出す方法、連続したエッジを抽出するＣａｎｎｙ法などが挙げられるが、画像中のエッジを抽出する方法であれば何でもよい。
この時、骨表面は超音波画像に対して水平方向に伸展するエッジが強勢であり、その事実からエッジ検出の確度を高めることが可能である。例えば、Ｓｏｂｅｌオペレータを用いて一次微分画像を作成する場合、超音波画像に対して水平方向に強度をもつエッジが強調されるようにオペレータの微分係数を式［数１］のように設定する。そうすることで、画像水平方向に伸展する骨表面のエッジが取り出されやすくなる。
上記方法を用いて水平方向のエッジを優勢に取りだしたとしても、図５（Ｂ）のように関節超音波画像には骨表面の他に皮膚、腱、関節包ら水平方向に伸展するエッジを持つ組織が撮像されている。従って、骨表面のみを特定するためには、超音波画像内のエッジを各組織、すなわち皮膚５００３、腱５００４、骨表面５００５に分類する必要がある。腱と関節包を区別して分類することはステップＳ４０１の段階では困難であるため、同じセグメントに分類しても構わない。エッジを分類するためには、例えばＫ−ｍｅａｎｓやｍｅａｎ−ｓｈｉｆｔのように、画像上で距離が近いエッジを同じセグメントに分類することが可能な手法であれば何でもよい。
上記方法によりエッジが各組織に分類されれば、骨表面以外の組織が描画された領域を超音波画像から除外することが可能となる。骨表面以外の組織が除外された超音波画像に対して、例えばガンマ補正や輝度ヒストグラム調整などの画像処理技術を用いて骨表面がより強調されるように画質を調整する。そして、画質調整した超音波画像に対して再び骨表面を特定することで、誤り無く骨表面を特定することができ検出精度が向上する。
ステップＳ４０１の別の方法は、例えばグラフカットや領域拡張法のように、骨を代表する画素から近距離に位置し、また近い画素値をもつ領域を抽出する方法であれば何でもよい。この時、骨表面を代表する画素は自動に決定してもよいし、外部入力装置１００９を介して検査者が入力してもよい。同様の方法を用いて骨表面以外の皮膚・腱などの各組織を特定しておいてもよい。
あるいは、予め保存しておいた代表的な骨表面の超音波画像をテンプレートとしてマッチング（剛体マッチングや非剛体マッチング）を行う、予め学習しておいた骨表面の特徴量を用いてＡｄａＢｏｏｓｔ、サポートベクタマシン、ニューラルネット、ランダムフォレストなどの機械学習方法により骨表面を検出する、などの方法もある。同様の方法を用いて骨表面以外の皮膚・腱などの各組織を特定しておいてもよい。
［ステップＳ４０２］
ステップＳ４０２ではステップＳ４０１で抽出した骨表面上にある関節位置を特定する。この処理の詳細を図６（Ａ）、図６（Ｂ）を用いて説明する。図６（Ａ）はステップＳ４０１で特定された骨表面６００１と骨表面６００２が描画されており、ステップＳ４０１で特定された皮膚５００３、及び腱５００４が除外されている。６００３、６００４、６００５はステップＳ４０２で特定される関節位置を示している。
ステップＳ４０２ではステップＳ４０１で抽出した骨表面上にある関節位置を特定する。この処理の詳細を図６（Ａ）、図６（Ｂ）を用いて説明する。図６（Ａ）はステップＳ４０１で特定された骨表面６００１と骨表面６００２が描画されており、ステップＳ４０１で特定された皮膚５００３、及び腱５００４が除外されている。６００３、６００４、６００５はステップＳ４０２で特定される関節位置を示している。
関節位置６００３と関節位置６００４は、骨表面６００１と骨表面６００２が連続していない断裂部を検出することで決定される。なお、予め保存しておいた各組織の代表的な超音波画像をテンプレートとしてマッチング（剛体マッチングや非剛体マッチング）を行う、予め学習しておいた各組織の特徴量を用いてＡｄａＢｏｏｓｔ、サポートベクタマシン、ニューラルネット、ランダムフォレストなどの機械学習方法により関節位置を特定する、などの方法もある。いずれの方法においても、図６（Ａ）のように二つの骨表面毎に代表画像あるいは特徴量を保存しておき、二つの関節位置を特定してもよいし、図６（Ｂ）のように関節の全体の代表画像あるいは特徴量を保存しておき、二つの骨表面を包括する関節位置６００５を特定しても構わない。
［ステップＳ４０３］
ステップＳ４０３ではステップＳ４０２で特定した関節位置の上方に位置する関節腔領域を抽出する。この処理の詳細を図７、図８（Ａ）、（Ｂ）、（Ｃ）を用いて説明する。図７の７００２と７００５は、関節腔領域７００１と骨表面７００６および骨表面７００７よりも上方にあるエッジである。その中でもエッジ７００５は、関節包のエッジである。ステップＳ４０１で特定された骨表面が７００６と７００７であり、関節包のエッジ７００５と骨表面７００６、７００７が接する地点がそれぞれ７００３と７００４である。
ステップＳ４０３ではステップＳ４０２で特定した関節位置の上方に位置する関節腔領域を抽出する。この処理の詳細を図７、図８（Ａ）、（Ｂ）、（Ｃ）を用いて説明する。図７の７００２と７００５は、関節腔領域７００１と骨表面７００６および骨表面７００７よりも上方にあるエッジである。その中でもエッジ７００５は、関節包のエッジである。ステップＳ４０１で特定された骨表面が７００６と７００７であり、関節包のエッジ７００５と骨表面７００６、７００７が接する地点がそれぞれ７００３と７００４である。
解剖学上、関節包２００６は骨２００３と２００４に接しているため（図２を参照）、超音波画像上での関節包のエッジも骨表面に接している。超音波画像上では関節腔領域は低輝度、関節包と骨表面は比較的高輝度に描画されるため、輝度差が顕著な境界を検出することで関節腔が検出可能である。従って、二値化手法、Ｓｏｂｅｌオペレータによる微分画像からのエッジ検出、Ｃａｎｎｙ法等を用いて関節位置より上方に位置するエッジ７００２を複数検出し、そのエッジの中からエッジの両端が骨表面７００６と骨表面７００７上の地点７００３と地点７００４に接しているエッジ７００５を選択することで超音波画像上の関節包を検出することができる。骨表面７００６、骨表面７００７はステップＳ４０１ですでに抽出しているため、この骨表面７００６と骨表面７００７と関節包のエッジ７００５で囲まれた領域７００１を関節腔領域として特定することが可能である。あるいは、骨表面７００６と骨表面７００７にエッジが接する地点７００３と地点７００４を検出し、この二つの地点を結ぶエッジ７００５を関節包として検出する方法でもよい。
関節包のエッジを検出するにあたっては、超音波画像から骨表面を除外し、この画像に対してガンマ補正や輝度ヒストグラム調整などの画像処理技術を用いて関節包のエッジがより強調するように調整された超音波画像に対して上記の方法を用いると、関節包の検出精度が向上する。
他の方法としては、動的輪郭モデルや領域拡張法等を用いて低輝度に描画される関節腔領域８００５を直接検出する方法もある。この場合図８に示されるように、初期探索点８００１を設定する必要があるが、これはステップＳ４０２で決定された関節位置から自動的に導かれる。例えば、関節位置８００２と関節位置８００３の中点から上方に所定の距離離れた位置に初期探索点８００１を置くことで実現される。なお、外部入力１００９を介して検査者が初期探索点を入力しても構わない。
決定された初期探索点８００１から反復的に探索領域８００４を拡大していき、最終的には輝度差が顕著な領域８００５が関節腔領域として特定される。なお、この方法で検出された関節腔領域のうち、関節包のエッジを検出した部分のみを抽出し、骨表面に関してはステップＳ４０１の結果を流用してもよい。その場合は、超音波画像から骨領域を除外し、この画像に対してガンマ補正や輝度ヒストグラ調整などの画像処理技術を用いて関節包のエッジがより強調されるように調整された超音波画像に対して上記の方法を用いると、関節包の検出精度が向上する。特定された関節包領域と超音波画像は画面作成部１０１０にて重畳され、重畳結果は表示部１０１２で表示することが可能である。
［ステップＳ４０４］
ステップＳ４０４ではステップＳ４０３にて特定された関節腔領域内を解析することで、病態の疾患活動性を定量評価する。評価方法を図９（Ａ）、図９（Ｂ）の超音波画像を用いて説明する。図９（Ａ）には、ステップＳ４０３で特定された関節腔領域９００１の他に、説明のために関節腔領域の高さ９００２、関節腔領域の幅９００３、ｘ軸方向９００４、ｙ軸方向９００５が描かれている。また、図９（Ｂ）は、説明のためにステップＳ４０２で検出された関節位置から関節包の両端へ向かう矢印９００６、関節腔領域の角度９００７が描かれている。
ステップＳ４０４ではステップＳ４０３にて特定された関節腔領域内を解析することで、病態の疾患活動性を定量評価する。評価方法を図９（Ａ）、図９（Ｂ）の超音波画像を用いて説明する。図９（Ａ）には、ステップＳ４０３で特定された関節腔領域９００１の他に、説明のために関節腔領域の高さ９００２、関節腔領域の幅９００３、ｘ軸方向９００４、ｙ軸方向９００５が描かれている。また、図９（Ｂ）は、説明のためにステップＳ４０２で検出された関節位置から関節包の両端へ向かう矢印９００６、関節腔領域の角度９００７が描かれている。
関節腔領域内の疾患活動性の定量評価例としては、関節腔領域９００１の面積、幅、高さ、角度が挙げられる。面積は、関節腔領域内のピクセル数をカウントすることで算出される。幅は、関節腔領域９００１のｘ軸方向９００４の座標の最大値と最小値を算出、あるいはステップＳ４０３で特定した骨表面と関節包の接地点の座標差を関節腔領域の幅９００３とすることで自動定量化が可能である。高さは、関節腔領域９００１のｙ軸方向９００５の座標の最大値と最小値を算出し、それらの値の差を関節腔領域の高さ９００２とすることで自動定量化が可能である。角度は、ステップＳ４０２で特定した関節位置から関節腔領域の両端へ向かう矢印９００６が成す角９００７の角度として定量化される。関節腔領域の両端はｘ軸方向の座標の最大値と最小値、あるいはステップＳ４０３で特定した骨表面と関節包の接地点座標でよい。
定量化の別の方法を図１０を用いて説明する。図１０の超音波画像にはステップＳ４０３で特定した関節腔領域１０００１と血流信号１０００２と血流信号１０００３が描かれている。疾患活動性の定量化は関節腔領域１０００１において観測される血流信号１０００２を計測することで可能である。関節腔領域外で観測される血流信号１０００３は計測しない。定量化の指標としては、関節腔領域１０００１において観測される血流信号１０００２の総面積、関節腔領域１０００１の面積に対する血流信号１０００２の総面積の割合等が挙げられる。
本発明にかかる超音波診断装置は、関節リウマチ診断における疾患活動性の評価に利用可能性がある。
１００１ 超音波プローブ
１００２ 超音波送受信部
１００３ 超音波画像生成部
１００４ 記憶部
１００５ 外部入力取得部
１００６ 骨表面特定部
１００７ 関節位置特定部
１００８ 関節腔領域特定部
１００９ 病態解析部
１０１０ 画面作成部
１０１１ 外部入力部
１０１２ 表示部
２００１ 骨
２００２ 骨
２００３ 軟骨
２００４ 軟骨
２００５ 滑膜
２００６ 関節包
３００１ 骨表面
３００２ 骨表面
３００３ 皮膚
３００４ 腱
３００５ 関節包
３００６ 関節腔領域
５００１ 骨表面
５００２ 骨表面
５００３ 皮膚に属するエッジ
５００４ 腱に属するエッジ
５００５ 骨に属するエッジ
６００１ 骨表面
６００２ 骨表面
６００３ 関節位置
６００４ 関節位置
６００５ 関節位置
７００１ 関節腔領域
７００２ エッジ
７００３ 地点
７００４ 地点
７００５ 関節包
７００６ 骨表面
７００７ 骨表面
８００１ 初期探索点
８００２ 関節位置
８００３ 関節位置
８００４ 探索領域
８００５ 関節腔領域
９００１ 関節腔領域
９００２ 関節腔の高さ
９００３ 関節腔領域の幅
９００４ ｘ軸方向
９００５ ｙ軸方向
９００６ 矢印
９００７ 関節腔領域の角度
１０００１ 関節腔領域
１０００２ 血流信号
１０００３ 血流信号
１００２ 超音波送受信部
１００３ 超音波画像生成部
１００４ 記憶部
１００５ 外部入力取得部
１００６ 骨表面特定部
１００７ 関節位置特定部
１００８ 関節腔領域特定部
１００９ 病態解析部
１０１０ 画面作成部
１０１１ 外部入力部
１０１２ 表示部
２００１ 骨
２００２ 骨
２００３ 軟骨
２００４ 軟骨
２００５ 滑膜
２００６ 関節包
３００１ 骨表面
３００２ 骨表面
３００３ 皮膚
３００４ 腱
３００５ 関節包
３００６ 関節腔領域
５００１ 骨表面
５００２ 骨表面
５００３ 皮膚に属するエッジ
５００４ 腱に属するエッジ
５００５ 骨に属するエッジ
６００１ 骨表面
６００２ 骨表面
６００３ 関節位置
６００４ 関節位置
６００５ 関節位置
７００１ 関節腔領域
７００２ エッジ
７００３ 地点
７００４ 地点
７００５ 関節包
７００６ 骨表面
７００７ 骨表面
８００１ 初期探索点
８００２ 関節位置
８００３ 関節位置
８００４ 探索領域
８００５ 関節腔領域
９００１ 関節腔領域
９００２ 関節腔の高さ
９００３ 関節腔領域の幅
９００４ ｘ軸方向
９００５ ｙ軸方向
９００６ 矢印
９００７ 関節腔領域の角度
１０００１ 関節腔領域
１０００２ 血流信号
１０００３ 血流信号
Claims (16)
- 超音波プローブを測定対象に当て、前記測定対象の超音波画像を得る超音波診断装置であって、
前記超音波画像から骨表面を特定する骨表面特定部と、
前記骨表面に存在する関節位置を特定する関節位置特定部と、
前記関節位置より浅い位置にある関節腔領域を特定する関節腔領域特定部と、
前記関節腔領域の内部に着目し、前記測定対象の疾患活動性を評価する病態解析部と、
を備えることを特徴とする超音波診断装置。 - 前記骨表面特定部は、前記超音波画像に対して水平方向に伸展するエッジを抽出し、抽出された前記エッジを所定の組織に対応したセグメントにセグメンテーションすることを特徴とする請求項１に記載の超音波診断装置。
- 前記所定の組織は骨表面、皮膚、腱、である請求項２に記載の超音波診断装置。
- 前記骨表面特定部は、前記超音波画像から、前記骨表面以外の前記所定の組織に対応するエッジを除外することを特徴とする請求項３に記載の超音波診断装置。
- 前記骨表面特定部は、前記骨表面以外の前記所定の組織に対応するエッジが除外された前記超音波画像の輝度を調整し、前記骨表面の輝度を強調することを特徴とする請求項４に記載の超音波診断装置。
- 前記関節位置特定部は、前記骨表面の断裂部分を検出することにより前記関節位置を特定することを特徴とする請求項１〜５のいずれかに記載の超音波診断装置。
- 前記関節位置特定部は、予め用意された超音波画像を用いたテンプレートマッチングにより前記関節位置を特定することを特徴とする請求項１〜５のいずれかに記載の超音波診断装置。
- 前記関節位置特定部は、予め設定された特徴量を用いた機械学習により前記関節位置を特定することを特徴とする請求項１〜５のいずれかに記載の超音波診断装置。
- 前記関節腔領域特定部は、前記関節位置より浅い位置にあり、かつ前記骨表面に接する連続したエッジを検出することを特徴とする請求項１〜８のいずれかに記載の超音波診断装置。
- 前記関節腔領域特定部は、前記超音波画像から前記骨表面を除去し、前記骨表面が除去された前記超音波画像に対して、関節包のエッジの輝度を強調することを特徴とする請求項１〜８のいずれかに記載の超音波診断装置。
- 前記関節腔領域特定部は、前記関節位置より浅い位置に初期輪郭を設定し、動的輪郭探索手法を用いて関節包を検出することを特徴とする請求項１〜８のいずれかに記載の超音波診断装置。
- 前記関節腔領域特定部は、前記関節位置より浅い位置に初期探索点を設定し、領域拡張法を用いて関節包を検出することを特徴とする請求項１〜８のいずれかに記載の超音波診断装置。
- 前記病態解析部は、前記関節包領域内の疾患活動性を定量化することを特徴とした請求項１〜１２のいずれかに記載の超音波診断装置。
- 前記病態解析部は、前記関節腔領域の高さ、幅、面積、のいずれかまたはその組み合わせで疾患活動性を定量化することを特徴とする請求項１〜１３のいずれかに記載の超音波診断装置。
- 超音波画像を解析する方法であって、
前記超音波画像から骨表面を特定する骨表面特定ステップと、
前記骨表面に存在する関節位置を特定する関節位置特定ステップと、
前記関節位置より浅い位置にある関節腔領域を特定する関節腔領域特定ステップと、
前記関節腔領域の高さ、幅、面積のいずれかを算出する関節腔領域解析ステップと、
を備えることを特徴とする超音波画像解析方法。 - 請求項１５に記載の超音波画像解析方法をコンピュータに実行させるためのプログラム。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013169527A JP6191328B2 (ja) | 2013-08-19 | 2013-08-19 | 超音波診断装置、超音波画像解析方法、およびプログラム |
| CN201410406234.6A CN104414684B (zh) | 2013-08-19 | 2014-08-18 | 超声波诊断装置、超声波诊断装置的图像处理方法 |
| US14/462,978 US9414811B2 (en) | 2013-08-19 | 2014-08-19 | Ultrasound diagnostic device, image processing method for ultrasound diagnostic device, and non-transitory computer-readable recording medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013169527A JP6191328B2 (ja) | 2013-08-19 | 2013-08-19 | 超音波診断装置、超音波画像解析方法、およびプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2015037491A true JP2015037491A (ja) | 2015-02-26 |
| JP6191328B2 JP6191328B2 (ja) | 2017-09-06 |
Family
ID=52631039
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2013169527A Active JP6191328B2 (ja) | 2013-08-19 | 2013-08-19 | 超音波診断装置、超音波画像解析方法、およびプログラム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6191328B2 (ja) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017169793A (ja) * | 2016-03-23 | 2017-09-28 | コニカミノルタ株式会社 | 超音波診断装置及び超音波画像処理方法 |
| JP2019520954A (ja) * | 2016-04-25 | 2019-07-25 | 中慧医学成像有限公司 | 脊柱の弯曲角度を測定する方法及びデバイス |
| JP2020151270A (ja) * | 2019-03-20 | 2020-09-24 | 学校法人慶應義塾 | 関節状態値取得装置、関節状態学習装置、関節位置特定装置、関節位置学習装置、関節状態値取得方法、関節状態学習方法、関節位置特定方法、関節位置学習方法及びプログラム |
| CN112790790A (zh) * | 2020-12-23 | 2021-05-14 | 居天医疗科技(深圳)有限公司 | 基于ai识别的运动关节损伤检测设备 |
| US11934944B2 (en) | 2018-10-04 | 2024-03-19 | International Business Machines Corporation | Neural networks using intra-loop data augmentation during network training |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002024075A1 (en) * | 2000-09-25 | 2002-03-28 | Insightec-Image Guided Treatment Ltd. | Detect surface points of an ultrasound image |
| US20040193048A1 (en) * | 2001-07-24 | 2004-09-30 | Liat Tsoref | Joint analysis using ultrasound |
| JP2008006188A (ja) * | 2006-06-30 | 2008-01-17 | Fujifilm Corp | 医用画像表示処理装置、及び、医用画像表示処理プログラム |
| US20090018445A1 (en) * | 2007-06-20 | 2009-01-15 | Perception Raisonnement Action En Medecine | Ultrasonic bone motion tracking system |
| WO2012046433A1 (ja) * | 2010-10-08 | 2012-04-12 | パナソニック株式会社 | 超音波診断装置、及び超音波診断方法 |
| JP2013022463A (ja) * | 2011-07-19 | 2013-02-04 | Toshiba Corp | 運動対象輪郭抽出装置、左心室画像分離装置、運動対象輪郭抽出方法及び左心室画像分離方法 |
| JP2013051998A (ja) * | 2011-08-31 | 2013-03-21 | Toshiba Corp | 超音波診断装置及び超音波診断装置の制御プログラム |
| JP2013056156A (ja) * | 2011-09-06 | 2013-03-28 | General Electric Co <Ge> | 病状を超音波式で自動式に検出、定量化及び追跡するための方法及びシステム |
| WO2013105197A1 (ja) * | 2012-01-10 | 2013-07-18 | パナソニック株式会社 | 超音波診断装置、および、血管検出方法 |
-
2013
- 2013-08-19 JP JP2013169527A patent/JP6191328B2/ja active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002024075A1 (en) * | 2000-09-25 | 2002-03-28 | Insightec-Image Guided Treatment Ltd. | Detect surface points of an ultrasound image |
| US20040193048A1 (en) * | 2001-07-24 | 2004-09-30 | Liat Tsoref | Joint analysis using ultrasound |
| JP2008006188A (ja) * | 2006-06-30 | 2008-01-17 | Fujifilm Corp | 医用画像表示処理装置、及び、医用画像表示処理プログラム |
| US20090018445A1 (en) * | 2007-06-20 | 2009-01-15 | Perception Raisonnement Action En Medecine | Ultrasonic bone motion tracking system |
| WO2012046433A1 (ja) * | 2010-10-08 | 2012-04-12 | パナソニック株式会社 | 超音波診断装置、及び超音波診断方法 |
| JP2013022463A (ja) * | 2011-07-19 | 2013-02-04 | Toshiba Corp | 運動対象輪郭抽出装置、左心室画像分離装置、運動対象輪郭抽出方法及び左心室画像分離方法 |
| JP2013051998A (ja) * | 2011-08-31 | 2013-03-21 | Toshiba Corp | 超音波診断装置及び超音波診断装置の制御プログラム |
| JP2013056156A (ja) * | 2011-09-06 | 2013-03-28 | General Electric Co <Ge> | 病状を超音波式で自動式に検出、定量化及び追跡するための方法及びシステム |
| WO2013105197A1 (ja) * | 2012-01-10 | 2013-07-18 | パナソニック株式会社 | 超音波診断装置、および、血管検出方法 |
Non-Patent Citations (2)
| Title |
|---|
| 深江淳: "関節エコー法の有用性", 分子リウマチ治療, vol. 1, no. 3, JPN6017014418, 2008, pages 126 - 129, ISSN: 0003544873 * |
| 西森美佐子: "ＪＳＳ中部 第１５回地方会「超音波検査による関節リウマチ評価」", 超音波検査技術, vol. 36, no. 6, JPN6017014417, 2011, pages 630 - 631, ISSN: 0003544872 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017169793A (ja) * | 2016-03-23 | 2017-09-28 | コニカミノルタ株式会社 | 超音波診断装置及び超音波画像処理方法 |
| JP2019520954A (ja) * | 2016-04-25 | 2019-07-25 | 中慧医学成像有限公司 | 脊柱の弯曲角度を測定する方法及びデバイス |
| US11934944B2 (en) | 2018-10-04 | 2024-03-19 | International Business Machines Corporation | Neural networks using intra-loop data augmentation during network training |
| JP2020151270A (ja) * | 2019-03-20 | 2020-09-24 | 学校法人慶應義塾 | 関節状態値取得装置、関節状態学習装置、関節位置特定装置、関節位置学習装置、関節状態値取得方法、関節状態学習方法、関節位置特定方法、関節位置学習方法及びプログラム |
| CN112790790A (zh) * | 2020-12-23 | 2021-05-14 | 居天医疗科技(深圳)有限公司 | 基于ai识别的运动关节损伤检测设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP6191328B2 (ja) | 2017-09-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6020470B2 (ja) | 超音波診断装置、および、血管特定方法 | |
| US8411927B2 (en) | Marker detection in X-ray images | |
| JP5975803B2 (ja) | 病状を超音波式で自動式に検出、定量化及び追跡するための方法及びシステム | |
| US10217213B2 (en) | Automatic focused assessment with sonography for trauma exams | |
| JP5751815B2 (ja) | 画像処理装置、撮影システム、画像処理方法及びプログラム | |
| JP6547612B2 (ja) | 画像処理装置、画像処理方法、および、画像処理装置を備える超音波診断装置 | |
| JP6191328B2 (ja) | 超音波診断装置、超音波画像解析方法、およびプログラム | |
| JP6839328B2 (ja) | 胎児超音波画像処理 | |
| US20160042525A1 (en) | Apparatus and method for visualization of region of interest | |
| TW201347737A (zh) | 乳房超音波影像掃描及診斷輔助系統 | |
| JP2000126182A (ja) | 腫瘍診断方法 | |
| CN102247172A (zh) | 自动化胎儿胎龄评估的系统和方法 | |
| US9014452B2 (en) | Orientation-aware average intensity histogram to indicate object boundary depth in ultrasound images | |
| CN111354006A (zh) | 超声图像中目标组织的描迹方法及装置 | |
| CN112568933B (zh) | 超声成像方法、设备和存储介质 | |
| CN113712594A (zh) | 医用图像处理装置以及医用摄像装置 | |
| US9414811B2 (en) | Ultrasound diagnostic device, image processing method for ultrasound diagnostic device, and non-transitory computer-readable recording medium | |
| Rahmatullah et al. | Anatomical object detection in fetal ultrasound: computer-expert agreements | |
| JP2015039466A (ja) | 超音波診断装置、画像処理方法、およびプログラム | |
| US20220249060A1 (en) | Method for processing 3d image data and 3d ultrasonic imaging method and system | |
| JP6061984B2 (ja) | 画像処理装置、画像処理方法及びプログラム | |
| JP6303905B2 (ja) | 超音波診断装置、超音波診断装置の画像処理方法、および、画像処理プログラム | |
| CN111383323A (zh) | 一种超声成像方法和系统以及超声图像处理方法和系统 | |
| WO2023032954A1 (ja) | 情報処理方法、プログラム及び画像診断装置 | |
| Doerfler et al. | Blood vessel detection in navigated ultrasound: An assistance system for liver resections |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160720 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20170329 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20170425 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20170622 |
|
| 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: 20170711 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170724 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6191328 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |