JP6598287B2 - 骨間距離測定装置、骨間距離測定方法、コンピュータを骨間距離測定装置として機能させるためのプログラム及び該プログラムを記憶した記録媒体 - Google Patents
骨間距離測定装置、骨間距離測定方法、コンピュータを骨間距離測定装置として機能させるためのプログラム及び該プログラムを記憶した記録媒体 Download PDFInfo
- Publication number
- JP6598287B2 JP6598287B2 JP2015022533A JP2015022533A JP6598287B2 JP 6598287 B2 JP6598287 B2 JP 6598287B2 JP 2015022533 A JP2015022533 A JP 2015022533A JP 2015022533 A JP2015022533 A JP 2015022533A JP 6598287 B2 JP6598287 B2 JP 6598287B2
- Authority
- JP
- Japan
- Prior art keywords
- edge
- joint
- interosseous
- distance
- ray image
- 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
- 238000000034 method Methods 0.000 title claims description 53
- 210000000988 bone and bone Anatomy 0.000 claims description 82
- 238000004364 calculation method Methods 0.000 claims description 21
- 230000006870 function Effects 0.000 claims description 17
- 230000035939 shock Effects 0.000 claims description 12
- 238000000691 measurement method Methods 0.000 claims description 7
- 238000000926 separation method Methods 0.000 claims description 7
- 238000010586 diagram Methods 0.000 description 13
- 238000005259 measurement Methods 0.000 description 12
- 210000003811 finger Anatomy 0.000 description 11
- 206010039073 rheumatoid arthritis Diseases 0.000 description 11
- 210000004932 little finger Anatomy 0.000 description 7
- 238000003384 imaging method Methods 0.000 description 6
- 206010023203 Joint destruction Diseases 0.000 description 4
- 210000001503 joint Anatomy 0.000 description 4
- 238000001134 F-test Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 230000002708 enhancing effect Effects 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 210000004247 hand Anatomy 0.000 description 3
- 206010061218 Inflammation Diseases 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000004054 inflammatory process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 208000023275 Autoimmune disease Diseases 0.000 description 1
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 210000000845 cartilage Anatomy 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 210000003414 extremity Anatomy 0.000 description 1
- 238000012905 input function Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 210000004934 left little finger Anatomy 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000008961 swelling Effects 0.000 description 1
- 210000003813 thumb Anatomy 0.000 description 1
Landscapes
- Length-Measuring Devices Using Wave Or Particle Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
関節のX線画像を入力するX線画像入力手段、
入力された関節のX線画像を表示する表示手段、
表示された関節のX線画像のエッジ点を選択する位置選択手段、
位置選択手段により選択されたエッジの座標情報に基づき、関節のエッジ曲線を演算する第1演算手段、及び、
第1演算手段により得られたエッジ曲線に基づき、骨間の距離を演算する第2演算手段、
を含む骨間距離測定装置。
(2)前記第2演算手段が、
第1演算手段が演算した2つのエッジ曲線の頂点間距離を少なくとも演算する、
上記(1)に記載の骨間距離測定装置。
(3)前記第2演算手段が、
第1演算手段が演算した2つのエッジ曲線間の面積に基づき骨間距離を演算する、
上記(1)に記載の骨間距離測定装置。
(4)関節のX線画像のエッジ点を選択する際の補助線を表示手段に表示するための補助線出力手段、
を含む上記(1)〜(3)の何れか一に記載の骨間距離測定装置。
(5)入力された関節のX線画像を拡大する拡大手段、
を含む上記(1)〜(4)の何れか一に記載の骨間距離測定装置。
(6)拡大手段が、
入力された関節のX線画像を骨格成分とテクスチャ成分に分離するTatal Variation正則化分離手段、及び
骨格成分のエッジ成分を鮮鋭化するshockフィルター、
を含む上記(5)に記載の骨間距離測定装置。
(7)関節のX線画像から骨間の距離を算出する骨間距離測定方法であって、
関節のX線画像を入力するX線画像入力ステップ、
入力された関節のX線画像を表示する表示ステップ、
表示された関節のX線画像のエッジ点を選択する位置選択ステップ、
位置選択ステップにより選択されたエッジ点の座標情報に基づき、関節のエッジ曲線を演算するステップ、及び、
得られたエッジ曲線に基づき、骨間の距離を演算するステップ、
を含む骨間距離測定方法。
(8)コンピュータを、上記(1)〜(6)の何れか一に記載の骨間距離測定装置として機能させるプログラム。
(9)上記(8)に記載されたプログラムが記録されたコンピュータ読み取り可能な記録媒体。
また、エッジ点を選択する際に補助線を表示することで、エッジを選択する範囲が同じになるので、測定のバラツキをより低下することができる。
更に、入力された関節のX線画像を骨格成分とテクスチャ成分に分離するTatal Variation正則化分離手段、及び骨格成分のエッジ成分を鮮鋭化するshockフィルターを用いてX線画像を拡大処理することで、解像度を極端に低下することなくX線画像を拡大できるので、患者の関節の骨間の距離の微妙な変化を正確に測定できる。
そして、本発明の骨間距離測定装置は、新規に撮影した関節のX線画像のみならず、過去のX線画像も解析できる。したがって、患者の撮影時期の異なるX線画像があれば、従来の方法では関節リウマチの発症を発見できない早期の患者であっても、X線画像を再解析することで、早期の関節リウマチ患者を発見できる。
・患者の関節付近を撮影したX線画像(以下、単に「X線画像」と記載することがある。)を入力するためのX線画像入力手段10、
・(1)入力されたX線画像、(2)表示されたX線画像に基づいて位置選択手段12により選択したエッジ点、及び(3)必要に応じて補助線、を表示する表示手段11、
・位置選択手段12により選択されたエッジ点の座標情報に基づき、関節のエッジ曲線を演算する第1演算手段13、
・第1演算手段により得られたエッジ曲線に基づき、骨間の距離を演算する第2演算手段14、
・X線画像入力手段10、表示手段11、位置選択手段12、第1演算手段13及び第2演算手段14を制御する制御手段15、
・制御手段15からアクセスされるプログラムメモリ16、
を少なくとも含んでいる。また、必要に応じて、後述する座標軸設定手段、補助線出力手段、X線画像の拡大手段を含んでいてもよい。
(2)選択したエッジ点21の各画素の列成分をx座標、行成分をy座標とする。
(3)選択した画素のうち、x座標の最小点と最大点の中点をx軸の原点とする。
(4)同様に、y軸の原点を求める。
(5)選択した各画素の座標を原点をもとに再度計算する。
(6)新たな座標値を用いて最小二乗法によってフィッティング関数(エッジ曲線22)を算出する。
(7)下の骨23についても、上記(1)〜(6)と同様の手順でエッジ点21を選択し、エッジ曲線24を算出する。
(2)選択点のうち2点の組み合わせ(図4中の点線)を考え、それぞれの距離を算出する。
(3)算出した距離の中でもっとも長い距離の2点の中点を原点32と定義する。
(4)2点を結ぶ直線をx軸33、原点32を通りx軸に直交する直線をy軸34と定義する。
(5)新たな座標軸をもとにフィッティングを実行してエッジ曲線35を算出する。
(6)下の骨36についても、(1)〜(5)と同様の手順でエッジ点31を選択し、エッジ曲線37を算出する。
(1)上記のカーブフィッティングにより得られた上の骨30のエッジ曲線pの頂点Aにおける法線lを算出する。
(2)上記のカーブフィッティングにより得られた下の骨36のエッジ曲線qの頂点Bにおける法線mを算出する。
(3)法線lとエッジ曲線qとの交点をC、法線mとエッジ曲線pとの交点をDと定義する。
(4)線分AC、BD、ABの距離をそれぞれ算出する。
(5)3つの線分の距離の平均を骨間距離として出力する。
左手中指及び小指のX線画像を装置1に入力し、先ず、TV正則化分離手段により骨格成分とテクスチャ成分に分離し、骨格成分のみをShockフィルターでエッジを強調することで、入力したX線画像を4×4倍に拡大した。図10(1)は拡大処理後の左手中指のX線画像、図10(2)は拡大処理後の左手小指のX線画像を表す。
次に、拡大処理をしていないX線画像と拡大処理したX線画像を用いた場合の違いを確認した。
図13(1)は拡大処理を行う前のX線画像、図13(2)は4×4倍に拡大処理を行った後のX線画像に対し、実施例1と同様の手順でエッジ曲線及び補助線として法線を表示した画像である。図13(1)及び(2)に示すように、大きさは全く異なっている。
次に、左右の手の親指を除く各指について、4×4倍に拡大処理を行ったX線画像と、拡大処理を行う前の通常X線画像を用いて、実施例1と同様の手順で骨間距離の測定を行った。
また、統計学的には2つの群の標準偏差の比であるF−testの値が、0.05以下であれば有意差ありとされている。拡大処理を行った場合と拡大処理を行わなかった2つの群でのF−testの差は8本の指の内、6本で0.05以下であった。したがって、拡大処理を行うことでより正確に骨間距離を測定できることが明らかとなった。
実施例2の骨間距離測定方法に代え、2つのエッジ曲線と補助線である2つの法線で囲まれている面積を、上の骨のエッジ曲線の両端部を結んだ長さで割ることで骨間距離を算出した以外は、実施例2と同様の手順で左右の指の骨間距離の測定を行った。表4に測定結果を示す。
Claims (7)
- 関節のX線画像を入力するX線画像入力手段、
入力された関節のX線画像を表示する表示手段、
表示された関節のX線画像のエッジ点を選択する位置選択手段、
位置選択手段により選択されたエッジの座標情報に基づき、関節のエッジ曲線を演算する第1演算手段、及び、
第1演算手段により得られたエッジ曲線に基づき、骨間の距離を演算する第2演算手段、
を含む、関節のX線画像から骨間の距離を算出する骨間距離測定装置であって、
前記第1演算手段は以下に記載する(1)〜(7)のアルゴリズムを用い、
(1)関節を構成する2つの骨のうち上の骨の着目するエッジ点を、前記位置選択手段を用いて複数点選択する、
(2)選択したエッジ点の各画素の列成分をx座標、行成分をy座標とする、
(3)選択した画素のうち、x座標の最小点と最大点の中点をx軸の原点とする、
(4)同様に、y軸の原点を求める、
(5)選択した各画素の座標を原点をもとに再度計算する、
(6)新たな座標値を用いて最小二乗法によってフィッティング関数(エッジ曲線)を算出する、
(7)関節を構成する2つの骨のうち下の骨についても、上記(1)〜(6)と同様の手順でエッジ点を選択し、エッジ曲線を算出する、
また、関節を形成する骨のエッジが傾いている場合には、前記(1)〜(7)のアルゴリズムの代わりに、以下に記載する(1)〜(6)のアルゴリズムを用い、
(1)関節を構成する2つの骨のうち上の骨のエッジ点を複数点選択する、
(2)選択点のうち2点の組み合わせを考え、それぞれの距離を算出する、
(3)算出した距離の中でもっとも長い距離の2点の中点を原点と定義する、
(4)前記2点を結ぶ直線をx軸、前記原点を通りx軸に直交する直線をy軸と定義する、
(5)新たな座標軸をもとにフィッティングを実行してエッジ曲線を算出する、
(6)関節を構成する2つの骨のうち下の骨についても、(1)〜(5)と同様の手順でエッジ点を選択し、エッジ曲線を算出し、
前記骨間距離は、
2つのエッジ曲線、及び2つのエッジ曲線の両端部を結んだ線で囲まれている範囲の面積を算出し、前記2つのエッジ曲線のうち一方のエッジ曲線の両端部を結んだ長さで面積を割った値を骨間距離とする、骨間距離測定装置。 - 関節のX線画像のエッジ点を選択する際の補助線を表示手段に表示するための補助線出力手段、
を含み、
前記補助手段は、エッジ曲線の両端部を通り、且つy軸と平行な直線を補助線として表示する、又は、エッジ曲線の両端部を通り、且つエッジ曲線の端部における法線となる直線を補助線として表示する手段である、請求項1に記載の骨間距離測定装置。 - 入力された関節のX線画像を拡大する拡大手段、
を含む請求項1又は2に記載の骨間距離測定装置。 - 拡大手段が、
入力された関節のX線画像を骨格成分とテクスチャ成分に分離するTatal Variation正則化分離手段、及び
骨格成分のエッジ成分を鮮鋭化するshockフィルター、
を含む請求項3に記載の骨間距離測定装置。 - 関節のX線画像を入力するX線画像入力ステップ、
入力された関節のX線画像を表示する表示ステップ、
表示された関節のX線画像のエッジ点を選択する位置選択ステップ、
位置選択ステップにより選択されたエッジ点の座標情報に基づき、関節のエッジ曲線を演算するステップ、及び、
得られたエッジ曲線に基づき、骨間の距離を演算するステップ、
を含む、関節のX線画像から骨間の距離を算出する骨間距離測定方法であって、
前記骨間の距離を演算するステップは以下に記載する(1)〜(7)のステップであり、
(1)関節を構成する2つの骨のうち上の骨の着目するエッジ点を、位置選択手段を用いて複数点選択する、
(2)選択したエッジ点の各画素の列成分をx座標、行成分をy座標とする、
(3)選択した画素のうち、x座標の最小点と最大点の中点をx軸の原点とする、
(4)同様に、y軸の原点を求める、
(5)選択した各画素の座標を原点をもとに再度計算する、
(6)新たな座標値を用いて最小二乗法によってフィッティング関数(エッジ曲線)を算出する、
(7)関節を構成する2つの骨のうち下の骨についても、上記(1)〜(6)と同様の手順でエッジ点を選択し、エッジ曲線を算出する、
また、関節を形成する骨のエッジが傾いている場合には、前記(1)〜(7)のステップの代わりに、以下に記載する(1)〜(6)ステップを用い、
(1)関節を構成する2つの骨のうち上の骨のエッジ点を複数点選択する、
(2)選択点のうち2点の組み合わせを考え、それぞれの距離を算出する、
(3)算出した距離の中でもっとも長い距離の2点の中点を原点と定義する、
(4)2点を結ぶ直線をx軸、原点を通りx軸に直交する直線をy軸と定義する、
(5)新たな座標軸をもとにフィッティングを実行してエッジ曲線を算出する、
(6)関節を構成する2つの骨のうち下の骨についても、(1)〜(5)と同様の手順でエッジ点を選択し、エッジ曲線を算出し、
前記骨間距離は、
2つのエッジ曲線、及び2つのエッジ曲線の両端部を結んだ線で囲まれている範囲の面積を算出し、前記2つのエッジ曲線のうち一方のエッジ曲線の両端部を結んだ長さで面積を割った値を骨間距離とする、骨間距離測定方法。 - コンピュータを、請求項1〜5の何れか一項に記載の骨間距離測定装置として機能させるプログラム。
- 請求項6に記載されたプログラムが記録されたコンピュータ読み取り可能な記録媒体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2015022533A JP6598287B2 (ja) | 2015-02-06 | 2015-02-06 | 骨間距離測定装置、骨間距離測定方法、コンピュータを骨間距離測定装置として機能させるためのプログラム及び該プログラムを記憶した記録媒体 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2015022533A JP6598287B2 (ja) | 2015-02-06 | 2015-02-06 | 骨間距離測定装置、骨間距離測定方法、コンピュータを骨間距離測定装置として機能させるためのプログラム及び該プログラムを記憶した記録媒体 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2016144535A JP2016144535A (ja) | 2016-08-12 |
JP6598287B2 true JP6598287B2 (ja) | 2019-10-30 |
Family
ID=56685098
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2015022533A Active JP6598287B2 (ja) | 2015-02-06 | 2015-02-06 | 骨間距離測定装置、骨間距離測定方法、コンピュータを骨間距離測定装置として機能させるためのプログラム及び該プログラムを記憶した記録媒体 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6598287B2 (ja) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6858047B2 (ja) * | 2017-03-27 | 2021-04-14 | 国立大学法人 東京大学 | 診断画像処理装置、評価支援方法及びプログラム |
JP7242041B2 (ja) * | 2019-04-05 | 2023-03-20 | 国立大学法人北海道大学 | 間隙変化検出装置、間隙変化検出方法及び間隙変化検出用プログラム |
KR20230023899A (ko) * | 2021-08-11 | 2023-02-20 | 주식회사 크레스콤 | 관절 상태를 정량화하기 위한 의료 영상 분석 방법, 의료 영상 분석 장치, 및 의료 영상 분석 시스템 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62183748A (ja) * | 1986-02-07 | 1987-08-12 | 帝人株式会社 | 骨のx線写真像の評価方法 |
US7668356B2 (en) * | 2005-08-03 | 2010-02-23 | Siemens Medical Solutions Usa, Inc. | Automatic determination of joint space width from hand radiographs |
JP2010253243A (ja) * | 2008-12-04 | 2010-11-11 | Fujifilm Corp | 関節裂隙幅計測システム及び関節裂隙幅計測方法及びプログラム |
JP2013152583A (ja) * | 2012-01-25 | 2013-08-08 | Nagoya Institute Of Technology | 画像処理装置 |
-
2015
- 2015-02-06 JP JP2015022533A patent/JP6598287B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
JP2016144535A (ja) | 2016-08-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8050734B2 (en) | Method and system for performing patient specific analysis of disease relevant changes of a disease in an anatomical structure | |
Choo et al. | Quantitative analysis of emphysema and airway measurements according to iterative reconstruction algorithms: comparison of filtered back projection, adaptive statistical iterative reconstruction and model-based iterative reconstruction | |
Safari et al. | A semi-automatic algorithm for estimating cobb angle | |
US20150262359A1 (en) | Image processing apparatus and computer-readable recording medium | |
US10140711B2 (en) | Method of analyzing vertebral edges of vertebrae in images acquired from an upright MRI system | |
JP6598287B2 (ja) | 骨間距離測定装置、骨間距離測定方法、コンピュータを骨間距離測定装置として機能させるためのプログラム及び該プログラムを記憶した記録媒体 | |
Kim et al. | Impact of radiation dose and iterative reconstruction on pulmonary nodule measurements at chest CT: a phantom study | |
Kim et al. | Effect of reconstruction parameters on the quantitative analysis of chest computed tomography | |
Anitha et al. | Automatic extraction of vertebral endplates from scoliotic radiographs using customized filter | |
Papaliodis et al. | Computer assisted Cobb angle measurements: a novel algorithm | |
US8644574B2 (en) | Measurement of thin-layered structures in X-ray computer tomography | |
Al-Bashir et al. | Computer-based Cobb angle measurement using deflection points in adolescence idiopathic scoliosis from radiographic images | |
JP2010259452A (ja) | 画像診断支援装置及び画像診断支援方法 | |
EP3849420B1 (en) | Calculation device for determining ventilation defects | |
Thalengala et al. | Computerized image understanding system for reliable estimation of spinal curvature in idiopathic scoliosis | |
Anitha et al. | Identification of apical vertebra for grading of idiopathic scoliosis using image processing | |
Lin et al. | Application of two-dimensional fractional-order convolution and bounding box pixel analysis for rapid screening of pleural effusion | |
Dovganich et al. | Automatic quality control in lung X-ray imaging with deep learning | |
Mahdavi et al. | A computer-based analysis for identification and quantification of small airway disease in lung computed tomography images: a comprehensive review for radiologists | |
JP2018192047A (ja) | 医用画像診断支援システム及び装置、並びに、プログラム | |
JP2006247293A (ja) | 画像処理方法および画像処理装置ならびにプログラム | |
JP5386160B2 (ja) | 人体構造物計測支援装置及び方法並びにプログラム | |
Levakhina et al. | A dual-axis tilt acquisition geometry for digital musculoskeletal tomosynthesis | |
JP2019053022A (ja) | 画像処理プログラム、画像処理装置、および画像処理方法 | |
JP2019052930A (ja) | 画像処理プログラム、画像処理装置、および画像処理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A80 | Written request to apply exceptions to lack of novelty of invention |
Free format text: JAPANESE INTERMEDIATE CODE: A80 Effective date: 20150227 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20180202 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20180206 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20180208 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20180202 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20181130 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20181218 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20190115 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190412 |
|
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: 20190917 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20190927 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6598287 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |