JP2010131201A - 人体構造物の解剖学的特徴点検出装置及び方法並びにプログラム - Google Patents
人体構造物の解剖学的特徴点検出装置及び方法並びにプログラム Download PDFInfo
- Publication number
- JP2010131201A JP2010131201A JP2008310122A JP2008310122A JP2010131201A JP 2010131201 A JP2010131201 A JP 2010131201A JP 2008310122 A JP2008310122 A JP 2008310122A JP 2008310122 A JP2008310122 A JP 2008310122A JP 2010131201 A JP2010131201 A JP 2010131201A
- Authority
- JP
- Japan
- Prior art keywords
- human body
- body structure
- feature point
- anatomical feature
- edge
- 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.)
- Abandoned
Links
- 238000000034 method Methods 0.000 title claims description 31
- 238000001514 detection method Methods 0.000 claims abstract description 77
- 238000012545 processing Methods 0.000 claims abstract description 40
- 230000005484 gravity Effects 0.000 claims abstract description 12
- 238000003708 edge detection Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 5
- 210000003484 anatomy Anatomy 0.000 claims description 3
- 238000003745 diagnosis Methods 0.000 abstract description 9
- 238000005259 measurement Methods 0.000 description 17
- 210000000988 bone and bone Anatomy 0.000 description 10
- 238000011156 evaluation Methods 0.000 description 8
- 210000004394 hip joint Anatomy 0.000 description 8
- 238000003384 imaging method Methods 0.000 description 6
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 210000000588 acetabulum Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 210000000629 knee joint Anatomy 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000399 orthopedic effect Effects 0.000 description 1
- 201000008482 osteoarthritis Diseases 0.000 description 1
- 210000004197 pelvis Anatomy 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 210000000689 upper leg Anatomy 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
【課題】医療画像における大腿骨骨頭中心等の人体構造物の解剖学的特徴点の位置の再現性を向上させ、より精度の高い診断を可能とする。
【解決手段】人体構造物を含む処理対象画像から人体構造物の輪郭形状を表すエッジを検出するエッジ検出部と、前記人体構造物の輪郭形状を表すエッジを閉曲線又は閉領域で近似する近似部と、前記閉曲線又は閉領域の重心を前記人体構造物の解剖学的特徴点として検出する特徴点検出部と、前記検出結果である、前記人体構造物の輪郭形状を表すエッジおよび/または、前記近似部により求められた前記閉曲線又は閉領域および/または、前記解剖学的特徴点を表示する表示部と、を備えたことを特徴とする人体構造物の解剖学的特徴点検出装置を提供することにより前記課題を解決する。
【選択図】図1
【解決手段】人体構造物を含む処理対象画像から人体構造物の輪郭形状を表すエッジを検出するエッジ検出部と、前記人体構造物の輪郭形状を表すエッジを閉曲線又は閉領域で近似する近似部と、前記閉曲線又は閉領域の重心を前記人体構造物の解剖学的特徴点として検出する特徴点検出部と、前記検出結果である、前記人体構造物の輪郭形状を表すエッジおよび/または、前記近似部により求められた前記閉曲線又は閉領域および/または、前記解剖学的特徴点を表示する表示部と、を備えたことを特徴とする人体構造物の解剖学的特徴点検出装置を提供することにより前記課題を解決する。
【選択図】図1
Description
本発明は、人体構造物の解剖学的特徴点検出装置及び方法並びにプログラムに係り、特に、大腿骨骨頭中心を検出する整形外科における計測支援技術に関する。
従来、股関節あるいは膝関節等における変形性関節症などの病気の診断は、主としてX線画像によって行われている。このとき医師の診察・判断をより正確なものとするために、例えば大腿骨骨頭中心等の解剖学的特徴点の位置などについての客観的診断情報を提供することが必要とされている。
これに対して従来、例えば骨盤等の診断対象を撮像した画像を表示部に表示し、これに対して例えば体の歪みを計測するための基準線及びその平行線などのオブジェクトを描画するための情報を入力部から入力して、表示部に表示された画像上にこのオブジェクトを描画し、描画されたオブジェクト間の幾何学的な関係を示す値を算出して、医師の診断に供するようにしたものが知られている(例えば、特許文献1等参照)。
ここでは、入力部から入力される入力情報は、ユーザ(医師等)の手入力による画像上の点であり、このユーザによって入力された点の位置情報をそのまま用いてオブジェクトが描画されている。
特開2007−151742号公報
しかしながら、例えば股関節領域の診察において用いられる評価値の一つにCE角(Center Edge Angle)があるが、これは大腿骨骨頭中心、臼蓋外端、及び左右涙痕下端等の位置情報に基づいて算出され、このCE角を算出する位置情報のうち大腿骨骨頭中心は、ある構造物の中心という定義であるため、画像上でその位置を示すための目印が存在せず、上記従来技術のようにユーザが上記位置情報を手入力する方法では、入力するユーザによって、また計測するタイミング等によってその都度大腿骨骨頭中心等の位置がばらついてしまい、正確な診断ができなくなってしまうという問題がある。
本発明は、このような事情に鑑みてなされたもので、医療画像における大腿骨骨頭中心等の人体構造物の解剖学的特徴点の位置の再現性を向上させ、より精度の高い診断を可能とする人体構造物の解剖学的特徴点検出装置及び方法並びにプログラムを提供することを目的とする。
前記目的を達成するために、請求項1に記載の発明は、人体構造物を含む処理対象画像から人体構造物の輪郭形状を表すエッジを検出するエッジ検出部と、前記人体構造物の輪郭形状を表すエッジを閉曲線又は閉領域で近似する近似部と、前記閉曲線又は閉領域の重心を前記人体構造物の解剖学的特徴点として検出する特徴点検出部と、前記検出結果である、前記人体構造物の輪郭形状を表すエッジおよび/または、前記近似部により求められた前記閉曲線又は閉領域および/または、前記解剖学的特徴点を表示する表示部と、を備えたことを特徴とする人体構造物の解剖学的特徴点検出装置を提供する。
これにより、人体構造物の解剖学的特徴点をユーザが直接入力した場合のバラツキを抑制することができ、その再現性を向上させることができる。
また、請求項2に示すように、前記表示部は、前記検出結果を前記処理対象画像に重ねて表示することを特徴とする。
このように、人体構造物の解剖学的特徴点の検出結果を処理対象画像に重ねて表示することにより、その検出結果の成否を確認することができる。
また、請求項3に示すように、請求項1または2に記載の人体構造物の解剖学的特徴点検出装置であって、前記表示部は前記処理対象画像を表示し、前記表示された処理対象画像中の人体構造物内に所定の点の入力を受け付ける入力受付手段を備え、前記所定の点に基づいて、前記人体構造物の輪郭形状を表すエッジを検出することを特徴とする。
このように、人体構造物内の任意の一点を指定するだけで、後は自動的に人体構造物の輪郭形状を表すエッジを検出することができる。
また、請求項4に示すように、請求項3に記載の人体構造物の解剖学的特徴点検出装置において、前記入力受付手段により入力を受け付ける前記所定の点が、前記解剖学的特徴点に略同一の点であり、前記エッジ検出部は、前記表示された処理対象画像中の人体構造物内に入力された略同一の点を中心とした放射状の直線上で、該直線上の他の点よりも画素濃度が大きく変化する点を、前記人体構造物の輪郭形状を表すエッジを構成する点であるとして検出することを特徴とする。
また、請求項5に示すように、請求項3に記載の人体構造物の解剖学的特徴点検出装置において、前記入力受付手段により入力を受け付ける前記所定の点が、前記人体構造物の輪郭形状を表すエッジを代表する1つ以上の点であり、該エッジを代表する1つ以上の点を用いて前記人体構造物の輪郭形状を表すエッジを検出することを特徴とする。
また、請求項6に示すように、前記エッジ検出部は、前記処理対象画像中から前記人体構造物の輪郭形状を表すエッジを自動で検出することを特徴とする。
このように、例えば公知の画像処理におけるエッジ検出技術を用いることにより、完全に自動で人体構造物の輪郭形状を表すエッジを検出することができる。
また、請求項7に示すように、前記近似部は、前記入力受付手段から入力受け付けされた1つ以上の点に基づいて、前記エッジを前記閉曲線または閉領域で近似することを特徴とする。
このように、入力手段からユーザがエッジ上の点を複数入力することにより半自動でエッジを検出することもできる。
また、請求項8に示すように、請求項1〜7のいずれかに記載の人体構造物の解剖学的特徴点検出装置であって、さらに、前記閉曲線または閉領域の再計算及び前記解剖学的特徴点の再検出を指示する手段を備え、前記再計算及び再検出を指示する手段からの指示により、前記近似部は前記閉曲線または閉領域の再計算を行うとともに、前記特徴点検出部は前記再計算された前記閉曲線または閉領域に基づいて前記解剖学的特徴点を再検出することを特徴とする。
このように、必要に応じて検出結果の修正を行うことにより、検出精度をより向上させることができる。
また、請求項9に示すように、請求項8に記載の人体構造物の解剖学的特徴点検出装置であって、さらに所定の点の入力を受け付ける入力受付手段を有する場合に、前記近似部は、前記閉曲線または閉領域の再計算を、前記入力受付手段から入力を受け付けて追加または修正された前記エッジ上の点を用いて行うことを特徴とする。
このように、検出結果の修正を行う際、ユーザが手入力することにより、検出精度をより高めることが可能となる。
また、請求項10に示すように、前記閉曲線は楕円であることを特徴とする。
このように、エッジを近似する閉曲線として楕円を用いることにより、その式による表示を利用することにより、その後の解剖学的特徴点を算出するのが容易となる。
また、請求項11に示すように、前記解剖学的特徴点は、大腿骨骨頭中心点であることを特徴とする。
これにより、特に股関節における疾病に本発明を好適に適用することができる。
また、同様に前記目的を達成するために、請求項12に記載の発明は、人体構造物を含む処理対象画像から人体構造物の輪郭形状を表すエッジを検出する工程と、前記人体構造物の輪郭形状を表すエッジを閉曲線又は閉領域で近似する工程と、前記閉曲線又は閉領域の重心を前記人体構造物の解剖学的特徴点として検出する工程と、前記検出結果である、前記人体構造物の輪郭形状を表すエッジおよび/または、前記近似部により求められた前記閉曲線又は閉領域および/または、前記解剖学的特徴点を表示する工程と、を備えたことを特徴とする人体構造物の解剖学的特徴点検出方法を提供する。
これにより、人体構造物の解剖学的特徴点をユーザが直接入力した場合のバラツキを抑制することができ、その再現性を向上させることができる。
また、請求項13に示すように、前記表示する工程は、前記検出結果を前記処理対象画像に重ねて表示することを特徴とする。
このように、人体構造物の解剖学的特徴点の検出結果を処理対象画像に重ねて表示することにより、その検出結果の成否を確認することができる。
また、請求項14に示すように、請求項12または13に記載の人体構造物の解剖学的特徴点検出方法であって、前記表示する工程は前記処理対象画像を表示し、さらに、前記表示された処理対象画像中の人体構造物内に所定の点を入力する工程を備え、前記エッジを検出する工程は、前記表示された処理対象画像中の人体構造物内に入力された所定の点を中心とした放射状の直線上で、該直線上の他の点よりも画素濃度が大きく変化する点を、前記人体構造物の輪郭形状を表すエッジを構成する点であるとして検出することを特徴とする。
これにより、人体構造物内の任意の一点を指定するだけで、後は自動的に人体構造物の輪郭形状を表すエッジを検出することができる。
また、請求項15に示すように、前記エッジを検出する工程は、前記処理対象画像中から前記人体構造物の輪郭形状を表すエッジを自動で検出することを特徴とする。
これにより、例えば公知の画像処理におけるエッジ検出技術を用いることにより、完全に自動で人体構造物の輪郭形状を表すエッジを検出することができる。
また、請求項16に示すように、前記近似する工程は、前記入力する工程において入力された点に基づき、前記検出されたエッジ上の点に基づいて、前記エッジを前記閉曲線または閉領域で近似することを特徴とする。
これにより、入力手段からユーザがエッジ上の点を複数入力することにより半自動でエッジを検出することもできる。
また、請求項17に示すように、請求項14〜16のいずれかに記載の人体構造物の解剖学的特徴点検出方法であって、さらに、前記閉曲線または閉領域の再計算を指示する工程を有し、前記近似する工程は、前記再計算の指示に基づいて、前記閉曲線または閉領域の再計算を行うとともに、前記特徴点を検出する工程は、前記再計算された前記閉曲線または閉領域に基づいて前記解剖学的特徴点を再検出することを特徴とする。
これにより、必要に応じて検出結果の修正を行うことにより、検出精度をより向上させることができる。
また、請求項18に示すように、前記近似する工程は、前記閉曲線または閉領域の再計算を、前記入力する工程において追加された前記エッジ上の点を用いて行うことを特徴とする。
このように、検出結果の修正を行う際、ユーザが手入力することにより、検出精度をより高めることが可能となる。
また、請求項19に示すように、前記閉曲線は楕円であることを特徴とする。
このように、エッジを近似する閉曲線として楕円を用いることにより、その式による表示を利用することにより、その後の解剖学的特徴点を算出するのが容易となる。
また、請求項20に示すように、前記解剖学的特徴点は、大腿骨骨頭中心であることを特徴とする。
これにより、特に股関節における疾病に本発明を好適に適用することができる。
また、同様に前記目的を達成するために、請求項21に記載の発明は、人体構造物を含む処理対象画像から人体構造物の輪郭形状を表すエッジを検出するステップと、前記人体構造物の輪郭形状を表すエッジを閉曲線又は閉領域で近似するステップと、前記閉曲線又は閉領域の重心を前記人体構造物の解剖学的特徴点として検出するステップと、前記検出結果である、前記人体構造物の輪郭形状を表すエッジおよび/または、前記近似部により求められた前記閉曲線又は閉領域および/または、前記解剖学的特徴点を表示するステップとからなる人体構造物の解剖学的特徴点検出方法をコンピュータに実行させることを特徴とするプログラムを提供する。
以上説明したように、本発明によれば、人体構造物の解剖学的特徴点をユーザが直接入力した場合のバラツキを抑制することができ、その再現性を向上させることができる。
以下、添付図面を参照して、本発明に係る人体構造物の解剖学的特徴点検出装置及び方法並びにプログラムについて詳細に説明する。
図1は、本発明に係る人体構造物の解剖学的特徴点検出装置の一実施形態の概略構成を示すブロック図である。
図1に示すように、本実施形態の人体構造物の解剖学的特徴点検出装置1(以下、単に特徴点検出装置1と言う。)は主に、画像入力部10、画像処理部12、近似部14、特徴点検出部16、計測値算出部18、情報入力部20及び表示部22から構成される。
画像入力部10は、患者の診察部位を撮像したX線画像データを入力するものである。画像入力部10へのX線画像データの入力方法は特に限定されるものではない。
例えば、X線画像データは、CR(コンピューテッド・ラジオグラフィック)システムを用いて骨部(例えば大腿骨など)の画像を含む放射線画像(CR画像)として得るようにしてもよい。すなわち、撮影装置において、X線管から発せられたX線を、被写体としての患者に向けて照射し、被写体を透過したX線を、蓄積性蛍光体を有するイメージングプレートに入射させて透過放射線画像をエネルギー情報として蓄積記録し、この透過放射線画像が蓄積記録されたイメージングプレートから透過放射線画像を読み取って画像入力部10に入力するようにしてもよい。
また画像入力部10は、FPD(Flat-Panel Detector)が搭載されたデジタルX線撮影装置からX線画像データを入力するようにしてもよいし、ネットワークを介して画像サーバからX線画像データを入力してもよいし、あるいはCD−ROM等の記録メディアからX線画像データを入力するようにしてもよい。
画像処理部12は、画像入力部10から入力されたX線画像データを受け取り、これに対して所定の画像処理を施すとともに、詳しくは後述するが、例えば大腿骨骨頭等の人体構造物のエッジ検出処理などを行うものである。
近似部14は、上で画像処理部12において検出された例えば大腿骨骨頭等の人体構造物のエッジをある曲線(閉曲線あるいは閉曲線の一部)で近似するものである。詳しくは後述するが、例えば大腿骨骨頭の場合には骨頭のエッジを近似する閉曲線としては楕円が最も好適に用いられる。
特徴点検出部16は、人体構造物の解剖学的特徴点を検出するものである。上で近似部14で人体構造物のエッジをある閉曲線で近似したとき、例えば、その閉曲線で囲まれる領域の重心を求めることにより解剖学的特徴点を検出する。詳しくは後述するが、大腿骨骨頭の場合には、解剖学的特徴点として骨頭中心を検出するが、例えば骨頭を楕円で近似したときには、その楕円の中心(重心)を骨頭中心として検出する。
計測値算出部18は、一つに絞られた計測候補点を用いて、人体構造物の所定の計測値を算出するものである。所定の計測値としては、例えば股関節の場合で言うと、CE角(Center Edge Angle)、Sharp角、AHI(Acetabular Head Index)などがある。
情報入力部20は、ユーザが、表示部22に表示された画像を見て、画像処理部12に対して、人体構造物のエッジ上の点と思われる点を入力したり、人体構造物の解剖学的特徴点であると大まかに思われる点を入力したり、また近似部14に対して、表示部22に表示された近似結果が思わしくない場合に、再度近似のやり直しを指示したり、その他システムに対する様々な指示を入力するものである。情報入力部20は、例えば、キーボードやマウス等の通常の入力のための操作部材で構成されるがこれに限定されるものではない。
表示部22は、LCD等の表示画面を有し、画像入力部10から入力され画像処理部12で所定の処理が施されたX線画像等の医用画像を表示するものである。表示部22の表示画面には、さらに、近似部14で検出した人体構造物のエッジを近似する閉曲線や、特徴点検出部16で検出した人体構造物の解剖学的特徴点などを上記医用画像に重ねて表示する。また、表示部22は、計測値算出部18で算出した計測値を表示する。
医師等のユーザは、この表示部22の表示画面を見て診断を行うことができる。また、図示を省略したが、本特徴点検出装置1は、表示画面に表示された画像等をプリントアウトするプリンタや、各種データを保存する記憶装置あるいはこれらのデータを各種記録媒体に記録する記録手段を備えていることが好ましい。
以下、本実施形態の作用について詳細に説明する。
図2は、本実施形態に係る人体構造物の解剖学的特徴点検出方法の処理の流れを示すフローチャートである。今例えば股関節のX線画像から大腿骨骨頭中心を検出する場合を例にとって説明することとする。
まず、図2のステップS100において、画像入力部10から計測対象画像を入力する。前述したように、この画像データの入力方法は様々な方法があり、特に限定されるものではない。入力された計測対象画像(X線画像データ)は、画像処理部12に送られ、所定の画像処理が施されて表示部22の表示画面に表示される。
図3に表示画面に表示された計測対象画像の例を示す。図3に示すように、今入力された計測対象画像は大腿骨骨頭が写った股関節正面画像である。
次に、ステップS110において、画像処理部12で大腿骨骨頭のエッジを検出する。
そこで、例えば、図4(a)に示すように、大腿骨骨頭が写った局所領域を切り出した画像を表示部22に表示する。このような局所領域画像を作成する方法はいろいろ考えられるが、例えば、大まかな骨頭中心位置を情報入力部20から入力し、その近傍を切り取る方法などがある。
また、図4(a)のような骨頭の局所領域を切り出した画像から骨頭のエッジを検出する方法もいろいろあり、特に限定されるものではない。一つの方法として、例えば図4(b)に示すように、骨頭内の大まかな中心点C等の任意の点を中心として、放射状に半直線Dを引き、その半直線D上における濃度変化の大きい点Eを探し、その点Eをエッジとしてもよい。このとき、骨頭内の大まかな中心点あるいは任意の点はユーザが情報入力部20から入力するようにする。
またあるいは、図4(b)の画像に対してSobelフィルタやLaplacianフィルタなどの一般的なエッジ抽出フィルタを適用してエッジを抽出する方法でもよい。
次に図2のステップS120において、上で求めたエッジを閉曲線で近似する。特にここではエッジを近似する閉曲線として楕円を用いるが、他の曲線でもよい。
図5にエッジを楕円で近似する方法を示す。
まず図5(a)に示すように、上で検出されたエッジEd上に複数個の点E1(x1,y1)、E2(x2,y2)、E3(x3,y3)、・・・、En(xn,yn)をユーザが情報入力部20から指定する。
なお、これらの点は、このようにエッジEd上で指定するのではなく、エッジEdを検出することなく、ユーザが表示部22に表示された画像を見て、情報入力部20からエッジ候補点として入力してもよい。
このエッジEdを近似する楕円(あるいはエッジ候補点がその上に乗るような楕円)の式を次のように与える。
{(x−xc)cosθ+(y−yc)sinθ}2/a2
+{(x−xc)sinθ−(y−yc)cosθ}2/b2 −1 = 0
ここで、図5(c)に示すように、(xc,yc)は楕円の中心の座標、aは長軸、bは短軸、θは中心の回りの長軸短軸の回転角を表している。なお、図5(c)では、説明のために、分かり易くするため楕円をかなり変形して表示しているが、実際の近似楕円はこれとは異なりかなり円に近い形となる。
+{(x−xc)sinθ−(y−yc)cosθ}2/b2 −1 = 0
ここで、図5(c)に示すように、(xc,yc)は楕円の中心の座標、aは長軸、bは短軸、θは中心の回りの長軸短軸の回転角を表している。なお、図5(c)では、説明のために、分かり易くするため楕円をかなり変形して表示しているが、実際の近似楕円はこれとは異なりかなり円に近い形となる。
上記式の左辺をf(xc,yc,θ,a,b,x,y)とおく。
そして、各エッジ上の点Ei(xi,yi)(i=1,・・・,n)について、次式を計算する。
g(xi,yi)=√[|f(xc,yc,θ,a,b,xi,yi)|]
ここで、記号|*|は*の絶対値であり、√[*]は、[*]の平方根を表す。
ここで、記号|*|は*の絶対値であり、√[*]は、[*]の平方根を表す。
ここでThrを所定の閾値とし、g(xi,yi)<Thrとなるエッジ上の点(もしくはエッジ候補点)Ei(xi,yi)の個数Nを算出する。
点Ei(xi,yi)が楕円上にあればf(xc,yc,θ,a,b,x,y)=0であるから、g(xi,yi)=0である。一般にはg(xi,yi)の値は点Ei(xi,yi)の楕円の周からのずれを表しており、これが所定の閾値Thrより小となる点の個数Nが大きい程、楕円f(xc,yc,θ,a,b,x,y)=0はエッジEdをよく近似することになる。
そこでg(xi,yi)<Thrとなる点Ei(xi,yi)の個数Nや、個数Nを楕円の面積で除した値や、g(xi,yi)(i=1〜N)それぞれにこの近似楕円からの距離の逆数を乗じた数の合計等を近似の評価値として算出する。
そして、xc,yc,θ,a,b等の各パラメータを変えながら、上記評価値が最大となるようなパラメータを探すことによって、図5(b)に示すように、エッジEdを近似する楕円Epを検出する。
ここで閾値Thrとしては、特に限定されるものではないが、例えば、楕円の長軸あるいは短軸の5%〜10%というようにしてもよい。
次に図2のステップS130において、骨頭の中心点を算出する。上で図5(b)に示すように、骨頭のエッジEdを近似する楕円Epが検出されているので、骨頭の中心点Cは、この近似楕円Epの中心(重心といっても同じ)として算出される。
次にステップS140において、これらの検出結果を表示部22の表示画面の画像上に重ねて表示する。
次に、ステップS150において、ユーザは表示部22に表示された結果を見て、その近似楕円でよいか判断する。結果がよければ骨頭中心点検出処理を終了する。もし上で検出された近似楕円Epから大きくはずれたエッジ上の点Ei(xi,yi)が多数あるような場合には、次のステップS160において近似結果の修正を行う。
すなわち、表示部22の表示画面で本来通るべき骨頭領域の輪郭上の点を情報入力部20から指定する。例えば、マウス等でクリックしてもよい。そして、その点の近くを通る程上記評価値が高くなるように再度パラメータxc,yc,θ,a,bを変えながら、近似楕円の評価基準を変更して再計算する。新たな評価値としては、例えば、上で算出した評価値に、いま指定した点と近似楕円との距離を除した値を用いること等が考えられる。評価値が最大となる楕円をあらたな近似楕円とする。
そして新たな近似楕円が検出されたら、ステップS170において、その近似楕円の中心として骨頭中心を検出し、ステップS180において結果を再度表示する。
なお、このように楕円でエッジを近似する場合には、低解像度の画像から徐々に解像度を上げて、近似精度を上げていくようにしてもよい。
また上で説明した例では、骨頭内の大まかな中心点と思われる点あるいは任意の点をユーザが指定し、その後自動的に、その点から放射状に半直線を引いて楕円によって近似したが、楕円以外の閉曲線でエッジを近似してもよい。この場合には閉曲線が囲む閉領域の重心として骨頭中心を検出する。またあるいは、エッジ検出から楕円近似してその中心を算出して骨頭中心を算出するまで完全自動で行ってもよい。またこのときエッジの近似を楕円以外の閉曲線で行い、その閉曲線で囲まれた閉領域の重心を自動で算出するようにしてもよい。
あるいは、このように自動的に閉曲線で近似することなく、ユーザがエッジ上の複数の候補点を手入力してエッジの近似形状を構成した後、自動で画像上のその点の座標を求めそれらの座標の重心として骨頭中心を検出するいわば半自動の方法であってもよい。
また、大腿骨骨頭のエッジEdが、図6(a)に破線30で示すように、骨頭の大腿骨側がエッジ候補点がない等の理由によって存在しない場合には、図6(b)に示すように、エッジEdの候補点がない部分を直線Lで置き換えて、骨頭のエッジEdを楕円Epの一部と直線Lで構成するようにする。
このように骨頭領域を楕円Epの一部と直線Lで囲まれた領域とすると、この領域の中心(重心)は、上で求めたような近似楕円Epの中心Cからずれていると考えられる。この場合には、例えば近似楕円Epの中心Cよりも短軸上に少しずれた点C’が骨頭の中心となると考えられる。この点は、例えば、近似楕円Epの一部及び直線Lで囲まれた領域の重心として算出される。あるいは、近似楕円Epの一部に対する直線Lの長さの割合(楕円Epの途切れ具合)から、楕円の軸上で中心をずらす量を決めるようにして算出してもよい。
このように本実施形態では、人体構造物の一例として大腿骨骨頭の形状(エッジ)を閉曲線(楕円)あるいはその一部または閉領域で近似し、その閉曲線で囲まれた領域(閉領域)の重心をその人体構造物の中心として検出するようにしている。また、検出結果を入力画像にオーバーレイ表示するようにしている。
以上のようにして大腿骨骨頭中心点Cが検出されたら、計測値算出部18は、これを用いてCE角を算出する。すなわち、CE角は、図3に示すように、左右の涙痕下端P1、P2を結ぶ直線に直交する直線と、臼蓋外端点Qと大腿骨骨頭中心点Cとを結ぶ直線とのなす角αとして算出される。
従って、本実施形態によれば、エッジを近似してそれを基に中心を検出するようにしているため、従来のようにユーザが人体構造物(大腿骨骨頭)の中心を直接入力した場合と比較して、ユーザ毎の手入力による中心点位置のバラツキを抑制することができる。
またその結果、大腿骨骨頭中心点位置の再現性が向上するため、CE角計測の再現性も向上し、より精度の高い診断が可能となる。
さらに、人体構造の近似形状とその中心点を入力画像にオーバーレイ表示することにより、中心点検出の成否を確認できるため、必要に応じて中心点検出の修正を行うことにより、中心点検出精度をより向上させることが可能となる。
また、本実施形態の人体構造物の解剖学的特徴点検出装置は、図示を省略した制御部に付属したメモリに格納されたプログラムによって制御される。すなわち、制御部によってメモリから人体構造物の解剖学的特徴点検出方法を実行するプログラムが読み出され、このプログラムに従って、例えば股関節のX線画像から大腿骨骨頭のエッジが検出され、骨頭の中心が検出される。
なお、本プログラムは、このように制御部に付属のメモリに格納されるものに限定されず、該プログラムを例えばPCカードやCD−ROMなど、本特徴点検出装置に着脱可能に構成されるメモリ媒体(リムーバブル媒体)に記録しておき、リムーバブル媒体に対応するインターフェイスを介して本装置に読み込むように構成してもよい。
以上、本発明の人体構造物の解剖学的特徴点検出装置及び方法並びにプログラムについて詳細に説明したが、本発明は、以上の例には限定されず、本発明の要旨を逸脱しない範囲において、各種の改良や変形を行ってもよいのはもちろんである。
1…人体構造物の解剖学的特徴点検出装置、10…画像入力部、12…画像処理部、14…近似部、16…特徴点検出部、18…計測値算出部、20…情報入力部、22…表示部
Claims (21)
- 人体構造物を含む処理対象画像から人体構造物の輪郭形状を表すエッジを検出するエッジ検出部と、
前記人体構造物の輪郭形状を表すエッジを閉曲線又は閉領域で近似する近似部と、
前記閉曲線又は閉領域の重心を前記人体構造物の解剖学的特徴点として検出する特徴点検出部と、
前記検出結果である、前記人体構造物の輪郭形状を表すエッジおよび/または、前記近似部により求められた前記閉曲線又は閉領域および/または、前記解剖学的特徴点を表示する表示部と、
を備えたことを特徴とする人体構造物の解剖学的特徴点検出装置。 - 前記表示部は、前記検出結果を前記処理対象画像に重ねて表示することを特徴とする請求項1に記載の人体構造物の解剖学的特徴点検出装置。
- 請求項1または2に記載の人体構造物の解剖学的特徴点検出装置であって、
前記表示部は前記処理対象画像を表示し、
前記表示された処理対象画像中の人体構造物内に所定の点の入力を受け付ける入力受付手段を備え、
前記所定の点に基づいて、前記人体構造物の輪郭形状を表すエッジを検出することを特徴とする人体構造物の解剖学的特徴点検出装置。 - 請求項3に記載の人体構造物の解剖学的特徴点検出装置において、
前記入力受付手段により入力を受け付ける前記所定の点が、前記解剖学的特徴点に略同一の点であり、
前記エッジ検出部は、前記表示された処理対象画像中の人体構造物内に入力された略同一の点を中心とした放射状の直線上で、該直線上の他の点よりも画素濃度が大きく変化する点を、前記人体構造物の輪郭形状を表すエッジを構成する点であるとして検出することを特徴とする人体構造物の解剖学的特徴点検出装置。 - 請求項3に記載の人体構造物の解剖学的特徴点検出装置において、
前記入力受付手段により入力を受け付ける前記所定の点が、前記人体構造物の輪郭形状を表すエッジを代表する1つ以上の点であり、
該エッジを代表する1つ以上の点を用いて前記人体構造物の輪郭形状を表すエッジを検出することを特徴とする人体構造物の解剖学的特徴点検出装置。 - 前記エッジ検出部は、前記処理対象画像中から前記人体構造物の輪郭形状を表すエッジを自動で検出することを特徴とする請求項1または2に記載の人体構造物の解剖学的特徴点検出装置。
- 前記近似部は、前記入力受付手段から入力受け付けされた1つ以上の点に基づいて、前記エッジを前記閉曲線または閉領域で近似することを特徴とする請求項3に記載の人体構造物の解剖学的特徴点検出装置。
- 請求項1〜7のいずれかに記載の人体構造物の解剖学的特徴点検出装置であって、さらに、前記閉曲線または閉領域の再計算及び前記解剖学的特徴点の再検出を指示する手段を備え、前記再計算及び再検出を指示する手段からの指示により、前記近似部は前記閉曲線または閉領域の再計算を行うとともに、前記特徴点検出部は前記再計算された前記閉曲線または閉領域に基づいて前記解剖学的特徴点を再検出することを特徴とする人体構造物の解剖学的特徴点検出装置。
- 請求項8に記載の人体構造物の解剖学的特徴点検出装置であって、さらに所定の点の入力を受け付ける入力受付手段を有する場合に、前記近似部は、前記閉曲線または閉領域の再計算を、前記入力受付手段から入力を受け付けて追加または修正された前記エッジ上の点を用いて行うことを特徴とする人体構造物の解剖学的特徴点検出装置。
- 前記閉曲線は楕円であることを特徴とする請求項1〜9のいずれかに記載の人体構造物の解剖学的特徴点検出装置。
- 前記解剖学的特徴点は、大腿骨骨頭中心点であることを特徴とする請求項1〜10のいずれかに記載の人体構造物の解剖学的特徴点検出装置。
- 人体構造物を含む処理対象画像から人体構造物の輪郭形状を表すエッジを検出する工程と、
前記人体構造物の輪郭形状を表すエッジを閉曲線又は閉領域で近似する工程と、
前記閉曲線又は閉領域の重心を前記人体構造物の解剖学的特徴点として検出する工程と、
前記検出結果である、前記人体構造物の輪郭形状を表すエッジおよび/または、前記近似部により求められた前記閉曲線又は閉領域および/または、前記解剖学的特徴点を表示する工程と、
を備えたことを特徴とする人体構造物の解剖学的特徴点検出方法。 - 前記表示する工程は、前記検出結果を前記処理対象画像に重ねて表示することを特徴とする請求項12に記載の人体構造物の解剖学的特徴点検出方法。
- 請求項12または13に記載の人体構造物の解剖学的特徴点検出方法であって、前記表示する工程は前記処理対象画像を表示し、さらに、前記表示された処理対象画像中の人体構造物内に所定の点を入力する工程を備え、前記エッジを検出する工程は、前記表示された処理対象画像中の人体構造物内に入力された所定の点を中心とした放射状の直線上で、該直線上の他の点よりも画素濃度が大きく変化する点を、前記人体構造物の輪郭形状を表すエッジを構成する点であるとして検出することを特徴とする人体構造物の解剖学的特徴点検出方法。
- 前記エッジを検出する工程は、前記処理対象画像中から前記人体構造物の輪郭形状を表すエッジを自動で検出することを特徴とする請求項12または13に記載の人体構造物の解剖学的特徴点検出方法。
- 前記近似する工程は、前記入力する工程において入力された点に基づき、前記検出されたエッジ上の点に基づいて、前記エッジを前記閉曲線または閉領域で近似することを特徴とする請求項14に記載の人体構造物の解剖学的特徴点検出方法。
- 請求項14〜16のいずれかに記載の人体構造物の解剖学的特徴点検出方法であって、さらに、前記閉曲線または閉領域の再計算を指示する工程を有し、前記近似する工程は、前記再計算の指示に基づいて、前記閉曲線または閉領域の再計算を行うとともに、前記特徴点を検出する工程は、前記再計算された前記閉曲線または閉領域に基づいて前記解剖学的特徴点を再検出することを特徴とする人体構造物の解剖学的特徴点検出方法。
- 前記近似する工程は、前記閉曲線または閉領域の再計算を、前記入力する工程において追加された前記エッジ上の点を用いて行うことを特徴とする請求項17に記載の人体構造物の解剖学的特徴点検出方法。
- 前記閉曲線は楕円であることを特徴とする請求項12〜18のいずれかに記載の人体構造物の解剖学的特徴点検出方法。
- 前記解剖学的特徴点は、大腿骨骨頭中心であることを特徴とする請求項12〜19のいずれかに記載の人体構造物の解剖学的特徴点検出方法。
- 人体構造物を含む処理対象画像から人体構造物の輪郭形状を表すエッジを検出するステップと、
前記人体構造物の輪郭形状を表すエッジを閉曲線又は閉領域で近似するステップと、
前記閉曲線又は閉領域の重心を前記人体構造物の解剖学的特徴点として検出するステップと、
前記検出結果である、前記人体構造物の輪郭形状を表すエッジおよび/または、前記近似部により求められた前記閉曲線又は閉領域および/または、前記解剖学的特徴点を表示するステップとからなる人体構造物の解剖学的特徴点検出方法をコンピュータに実行させることを特徴とするプログラム。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008310122A JP2010131201A (ja) | 2008-12-04 | 2008-12-04 | 人体構造物の解剖学的特徴点検出装置及び方法並びにプログラム |
US12/591,875 US8696603B2 (en) | 2008-12-04 | 2009-12-03 | System for measuring space width of joint, method for measuring space width of joint and recording medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008310122A JP2010131201A (ja) | 2008-12-04 | 2008-12-04 | 人体構造物の解剖学的特徴点検出装置及び方法並びにプログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2010131201A true JP2010131201A (ja) | 2010-06-17 |
Family
ID=42343154
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008310122A Abandoned JP2010131201A (ja) | 2008-12-04 | 2008-12-04 | 人体構造物の解剖学的特徴点検出装置及び方法並びにプログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2010131201A (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015171456A (ja) * | 2014-03-12 | 2015-10-01 | 株式会社東芝 | 医用画像表示装置、医用画像表示システムおよび医用画像表示プログラム |
JP2021115120A (ja) * | 2020-01-23 | 2021-08-10 | コニカミノルタ株式会社 | 画像診断支援装置及びプログラム |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06261894A (ja) * | 1993-01-18 | 1994-09-20 | Hitachi Ltd | 骨塩定量方法 |
JP2002541950A (ja) * | 1999-04-20 | 2002-12-10 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 関連する幾何学的物体を集合的に構成する方法及び装置 |
JP2003144454A (ja) * | 2001-11-16 | 2003-05-20 | Yoshio Koga | 関節手術支援情報算出方法、関節手術支援情報算出プログラム、及び関節手術支援情報算出システム |
JP2007111123A (ja) * | 2005-10-18 | 2007-05-10 | Hitachi Medical Corp | 医用画像表示装置 |
JP2008173236A (ja) * | 2007-01-17 | 2008-07-31 | Toshiba Corp | 医用画像表示装置 |
-
2008
- 2008-12-04 JP JP2008310122A patent/JP2010131201A/ja not_active Abandoned
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06261894A (ja) * | 1993-01-18 | 1994-09-20 | Hitachi Ltd | 骨塩定量方法 |
JP2002541950A (ja) * | 1999-04-20 | 2002-12-10 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 関連する幾何学的物体を集合的に構成する方法及び装置 |
JP2003144454A (ja) * | 2001-11-16 | 2003-05-20 | Yoshio Koga | 関節手術支援情報算出方法、関節手術支援情報算出プログラム、及び関節手術支援情報算出システム |
JP2007111123A (ja) * | 2005-10-18 | 2007-05-10 | Hitachi Medical Corp | 医用画像表示装置 |
JP2008173236A (ja) * | 2007-01-17 | 2008-07-31 | Toshiba Corp | 医用画像表示装置 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015171456A (ja) * | 2014-03-12 | 2015-10-01 | 株式会社東芝 | 医用画像表示装置、医用画像表示システムおよび医用画像表示プログラム |
JP2021115120A (ja) * | 2020-01-23 | 2021-08-10 | コニカミノルタ株式会社 | 画像診断支援装置及びプログラム |
JP7434923B2 (ja) | 2020-01-23 | 2024-02-21 | コニカミノルタ株式会社 | 画像診断支援装置及びプログラム |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8696603B2 (en) | System for measuring space width of joint, method for measuring space width of joint and recording medium | |
JP6527209B2 (ja) | 画像表示の生成方法 | |
EP3025304B1 (en) | Registration system for registering an imaging device with a tracking device | |
CN106232010B (zh) | 用于检测气管的***和方法 | |
JP5643304B2 (ja) | 胸部トモシンセシスイメージングにおけるコンピュータ支援肺結節検出システムおよび方法並びに肺画像セグメント化システムおよび方法 | |
US11490986B2 (en) | System and method for improved electronic assisted medical procedures | |
JP4265876B2 (ja) | 画像を幾何学的に計測するための計測処理装置 | |
EP2564788B1 (en) | Information processing apparatus, ultrasonic imaging apparatus, and information processing method | |
JP2010253243A (ja) | 関節裂隙幅計測システム及び関節裂隙幅計測方法及びプログラム | |
US10617474B2 (en) | System and method for optimizing an implant position in an anatomical joint | |
US20110282199A1 (en) | Ultrasound system for measuring image using figure template and method for operating ultrasound system | |
JP2007130462A (ja) | 左右対称の器官系の三次元画像の評価方法、医用画像撮影システムおよびコンピュータプログラム | |
JP7337556B2 (ja) | 医用画像処理装置、x線診断装置及び医用画像処理方法 | |
JP2012115635A (ja) | 画像処理方法、画像処理装置、画像システム及びプログラムコード | |
JP2018068400A (ja) | 動態画像処理装置 | |
US20200121279A1 (en) | Method and device for measuring spinal column curvature | |
JP2007105352A (ja) | 差分画像表示装置、差分画像表示方法およびそのプログラム | |
JP5576631B2 (ja) | 放射線撮影装置、放射線撮影方法、及びプログラム | |
JP2009112531A (ja) | 画像表示装置 | |
US10083503B2 (en) | Image area specification device and method, and X-ray image processing device and method | |
JP2010131201A (ja) | 人体構造物の解剖学的特徴点検出装置及び方法並びにプログラム | |
JP2005270635A (ja) | 画像処理方法及び処理装置 | |
US10803645B2 (en) | Visualization of anatomical cavities | |
JP5386160B2 (ja) | 人体構造物計測支援装置及び方法並びにプログラム | |
EP4128145B1 (en) | Combining angiographic information with fluoroscopic images |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20110818 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130227 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130228 |
|
A762 | Written abandonment of application |
Free format text: JAPANESE INTERMEDIATE CODE: A762 Effective date: 20130411 |