JP7242852B2 - 関節投影面設定装置、方法およびプログラム - Google Patents
関節投影面設定装置、方法およびプログラム Download PDFInfo
- Publication number
- JP7242852B2 JP7242852B2 JP2021524681A JP2021524681A JP7242852B2 JP 7242852 B2 JP7242852 B2 JP 7242852B2 JP 2021524681 A JP2021524681 A JP 2021524681A JP 2021524681 A JP2021524681 A JP 2021524681A JP 7242852 B2 JP7242852 B2 JP 7242852B2
- Authority
- JP
- Japan
- Prior art keywords
- joint
- projection
- image
- projection plane
- cartilage
- 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 42
- PWPJGUXAGUPAHP-UHFFFAOYSA-N lufenuron Chemical compound C1=C(Cl)C(OC(F)(F)C(C(F)(F)F)F)=CC(Cl)=C1NC(=O)NC(=O)C1=C(F)C=CC=C1F PWPJGUXAGUPAHP-UHFFFAOYSA-N 0.000 title 1
- 210000000845 cartilage Anatomy 0.000 claims description 150
- 210000002303 tibia Anatomy 0.000 claims description 49
- 238000011002 quantification Methods 0.000 claims description 43
- 238000013507 mapping Methods 0.000 claims description 31
- 230000005484 gravity Effects 0.000 claims description 18
- 230000007717 exclusion Effects 0.000 claims description 17
- 239000000284 extract Substances 0.000 claims description 12
- 238000009795 derivation Methods 0.000 claims description 11
- 235000010893 Bischofia javanica Nutrition 0.000 claims description 10
- 240000005220 Bischofia javanica Species 0.000 claims description 10
- 230000007547 defect Effects 0.000 claims description 9
- 238000003384 imaging method Methods 0.000 claims description 8
- 238000010586 diagram Methods 0.000 description 34
- 238000012545 processing Methods 0.000 description 27
- 238000003860 storage Methods 0.000 description 19
- 210000000988 bone and bone Anatomy 0.000 description 14
- 210000000629 knee joint Anatomy 0.000 description 14
- 238000002595 magnetic resonance imaging Methods 0.000 description 8
- 210000000689 upper leg Anatomy 0.000 description 8
- 238000003745 diagnosis Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 7
- 210000005065 subchondral bone plate Anatomy 0.000 description 5
- 238000000605 extraction Methods 0.000 description 4
- 210000000056 organ Anatomy 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 3
- 238000009877 rendering Methods 0.000 description 3
- 239000003086 colorant Substances 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000012447 hatching Effects 0.000 description 2
- 230000005499 meniscus Effects 0.000 description 2
- 201000008482 osteoarthritis Diseases 0.000 description 2
- 210000000426 patellar ligament Anatomy 0.000 description 2
- 238000002600 positron emission tomography Methods 0.000 description 2
- 210000002967 posterior cruciate ligament Anatomy 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 102000008186 Collagen Human genes 0.000 description 1
- 108010035532 Collagen Proteins 0.000 description 1
- 208000003947 Knee Osteoarthritis Diseases 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229920001436 collagen Polymers 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 210000002436 femur neck Anatomy 0.000 description 1
- 210000003127 knee Anatomy 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 210000004417 patella Anatomy 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Optics & Photonics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Description
本開示は、関節を含む３次元画像から関節の投影画像を生成する際の投影面を設定する関節投影面設定装置、方法およびプログラムに関する。
近年、ＣＴ（Computed Tomography）装置およびＭＲＩ（Magnetic Resonance Imaging）装置等の医療機器の進歩により、質の高い高解像度の３次元画像が画像診断に用いられるようになってきている。ここで、３次元画像は多数の２次元画像から構成され情報量が多いため、医師が所望の観察部位を見つけて診断することに時間を要する場合がある。そこで、注目する臓器を認識し、注目する臓器を含む３次元画像から、例えば最大値投影法（ＭＩＰ（Maximum Intensity Projection）法）および最小値投影法（ＭｉｎＩＰ（Minimum Intensity Projection）法）等の方法を用いて、注目する臓器を抽出してＭＩＰ表示等を行ったり、３次元画像のボリュームレンダリング（ＶＲ（Volume Rendering））表示を行ったりすることにより、臓器全体や病変の視認性を高め、診断の効率化を図ることが行われている。
一方、変形性関節症は、高齢者に多く発症する疾病である。とくに変形性膝関節症は、膝関節の痛みおよび動作範囲の減少を引き起こすため、症状が進行すると歩けなくなることもある。このような変形性関節症の診断には、膝関節の軟骨を質的に評価する必要がある。このため、３次元画像を用いて関節の軟骨を定量化する手法が種々提案されている。例えば、特許文献１においては、ＭＲＩ画像から抽出された軟骨領域の投影方向を決定し、決定した投影方向に軟骨領域を投影して投影画像を生成し、投影画像上において軟骨領域の定量値を導出する手法が提案されている。特許文献１に記載された手法を用いることにより、軟骨領域を定量化するための領域を適切に定めることができるため、軟骨について安定した診断結果を得ることができる。とくに、軟骨の厚さを定量値とすることにより、軟骨の厚さの評価を行うことが可能となる。
また、特許文献２には、大腿骨頸部のような球状の関節を識別するために、関節の３次元画像から関節骨の頭部中心および半径並びに頸部軸を判定する手法が提案されている。特許文献２に記載の手法においては、骨の頸部部分の３次元表面モデル上で３次元頸部最小曲線を判定し、３次元頸部最小曲線に対する最小二乗法適合平面を判定し、最小二乗適合平面に直交する方向を、精密な頸部軸の方向として計算し、最小二乗適合平面上の３次元頸部最小曲線の投影の中心を、精密な頸部軸の点として計算している。
一方、関節は曲面を含む。とくに、脛骨の関節は、平面であるように見えるが、内側顆および外側顆で傾きが異なり、かつ関節面が凹形状を有する。このため、特許文献１に記載された手法を用いたとしても、投影方向を適切に設定しないと、軟骨領域の定量化の誤差が大きくなる。例えば、図２１に示すように、脛骨９０の外側顆側の関節面にある軟骨９１が、脛骨９０の中心軸Ｘ０に対して傾斜する矢印Ａ方向に投影されたとする。この場合、図２２に示すように、斜線で示す軟骨９１の欠損部分９２が、投影画像９３において軟骨９１が存在するように投影されてしまう。また、特許文献２に記載された手法は、球形状の関節における軸を判定しており、関節を投影するための投影面を適切に設定していない。
本開示は上記事情に鑑みなされ、関節を含む３次元画像から関節の投影画像を生成する際の投影面を適切に設定することを目的とする。
本開示による関節投影面設定装置は、関節の３次元画像を取得する画像取得部と、
関節の関節面を近似する平面を、３次元画像を投影して投影画像を生成するための投影面に設定する面設定部とを備える。
関節の関節面を近似する平面を、３次元画像を投影して投影画像を生成するための投影面に設定する面設定部とを備える。
なお、本開示による関節投影面設定装置においては、面設定部は、関節面を構成する複数のボクセルのうち、投影面からの距離がボクセル除外距離以上となるボクセルを除外して、新たな投影面を設定してもよい。
また、本開示による関節投影面設定装置においては、面設定部は、ボクセルの除外および新たな投影面の設定を複数回繰り返して、投影面を設定してもよい。
また、本開示による関節投影面設定装置においては、面設定部は、最小二乗法により関節面を近似する平面を投影面に設定してもよい。
また、本開示による関節投影面設定装置においては、関節は脛骨の関節であってもよい。
また、本開示による関節投影面設定装置においては、面設定部は、脛骨の関節面を顆間領域設定方向から見た画像において、顆間の領域を除いた領域を関節面領域として抽出し、関節面領域に含まれる関節面を近似する平面を、投影面に設定してもよい。
また、本開示による関節投影面設定装置においては、３次元画像の関節にある軟骨を投影面に直交する方向に投影することにより投影画像を生成する投影部をさらに備えてもよい。
また、本開示による関節投影面設定装置においては、脛骨の内側顆の関節面の重心および外側顆の関節面の重心を結ぶ線が、水平方向を向くように、３次元画像の関節にある軟骨を投影面に直交する方向に投影することにより投影画像を生成する投影部をさらに備えてもよい。
また、本開示による関節投影面設定装置においては、投影面に垂直な方向から見た脛骨における赤木ラインが、垂直方向を向くように、３次元画像の関節にある軟骨を投影面に直交する方向に投影することにより投影画像を生成する投影部をさらに備えてもよい。
「赤木ライン」とは、後十字靱帯付着部と膝蓋腱付着部内側縁を結んだ線である。
また、本開示による関節投影面設定装置においては、投影画像を表示部に表示する表示制御部をさらに備えてもよい。
また、本開示による関節投影面設定装置においては、投影画像上において、軟骨の定量値を導出する定量化部をさらに備えてもよい。
また、本開示による関節投影面設定装置においては、定量化部は、投影画像上の関心領域において、定量値を導出してもよい。
また、本開示による関節投影面設定装置においては、投影部は、関節における軟骨が存在すべき領域を規定する輪郭に基づいて関心領域を設定してもよい。
「軟骨が存在すべき領域を規定する輪郭」とは、関節において実際に軟骨が存在する領域ではなく、関節において解剖学的に軟骨が存在すべき領域を規定する輪郭を意味する。例えば、脛骨においては、関節において軟骨が存在すべき領域の輪郭が、関節面に凸状部分として含まれる。このため、脛骨においては、この凸状部分を「軟骨が存在すべき領域を規定する輪郭」とすることができる。
また、本開示による関節投影面設定装置においては、定量化部は、関心領域内における軟骨の被覆率を、定量値として導出してもよい。
また、本開示による関節投影面設定装置においては、定量化部は、関心領域における軟骨の欠損面積を、定量値として導出してもよい。
また、本開示による関節投影面設定装置においては、定量化部は、関心領域内の各位置における軟骨の厚さの代表値を、定量値として導出してもよい。
また、本開示による関節投影面設定装置においては、定量化部は、関心領域内の各位置における軟骨の厚さを、定量値として導出してもよい。
また、本開示による関節投影面設定装置においては、定量化部は、関心領域における軟骨の厚さマップを生成してもよい。
また、本開示による関節投影面設定装置においては、厚さマップを表示部に表示する表示制御部をさらに備えてもよい。
また、本開示による関節投影面設定装置においては、定量化部は、関心領域における軟骨の厚さが、定量値導出対象厚さ以上となる領域においてのみ、定量値を導出してもよい。
また、本開示による関節投影面設定装置においては、定量化部は、投影画像上の関心領域を分割し、分割により得られる各領域において、定量値を導出してもよい。
また、本開示による関節投影面設定装置においては、３次元画像を取得した被検体と同一の被検体についての撮影時期が異なる他の３次元画像から導出した定量値の導出結果が存在する場合、投影部は、他の３次元画像から定量値を導出した際と同一位置に関心領域を設定してもよい。
また、本開示による関節投影面設定装置においては、３次元画像を取得した被検体と同一の被検体についての撮影時期が異なる他の３次元画像の投影画像に関心領域が設定されている場合、投影部は、３次元画像の投影画像に対して、他の３次元画像の投影画像の関心領域と同一の関心領域を設定してもよい。
また、本開示による関節投影面設定装置においては、３次元画像を取得した被検体と同一の被検体についての撮影時期が異なる他の３次元画像の投影画像に関心領域が設定されている場合、投影部は、３次元画像の投影画像に対して、他の３次元画像の投影画像の関心領域と同一の関心領域を設定してもよい。
また、本開示による関節投影面設定装置においては、定量化部は、投影画像上の軟骨の面積を定量値として導出してもよい。
また、本開示による関節投影面設定装置においては、定量化部は、投影画像上の軟骨の体積を定量値として導出してもよい。
また、本開示による関節投影面設定装置においては、３次元画像を取得した被検体と同一の被検体についての撮影時期が異なる他の３次元画像から生成した他の投影画像が存在する場合、面設定部は、他の投影画像を生成した投影面を、３次元画像を投影して投影画像を生成するための投影面に設定してもよい。
また、本開示による関節投影面設定装置においては、関節についての機能画像を投影画像にマッピングしたマッピング画像を生成するマッピング部をさらに備えてもよい。
「機能画像」とは、被検体の関節において、水およびコラーゲン等の含有される物質の量または相関性等に応じた信号値を有する画像である。具体的には、ＭＲＩにより取得されるＴ１マップ、Ｔ２マップ、Ｔ２＊マップ、Ｔ１ρマップ、およびＣＥＳＴ(Chemical Exchange Saturation Transfer)等を機能画像として用いることができる。
また、本開示による関節投影面設定装置においては、マッピング画像を表示部に表示する表示制御部をさらに備えてもよい。
本開示による関節投影面設定方法は、関節の３次元画像を取得し、
関節の関節面を近似する平面を、３次元画像を投影して投影画像を生成するための投影面に設定する。
関節の関節面を近似する平面を、３次元画像を投影して投影画像を生成するための投影面に設定する。
なお、本開示による関節投影面設定方法をコンピュータに実行させるためのプログラムとして提供してもよい。
本開示による他の関節投影面設定装置は、
コンピュータに実行させるための命令を記憶するメモリと、
記憶された命令を実行するよう構成されたプロセッサとを備え、プロセッサは、
関節の３次元画像を取得し、
関節の関節面を近似する平面を、３次元画像を投影して投影画像を生成するための投影面に設定する処理を実行する。
コンピュータに実行させるための命令を記憶するメモリと、
記憶された命令を実行するよう構成されたプロセッサとを備え、プロセッサは、
関節の３次元画像を取得し、
関節の関節面を近似する平面を、３次元画像を投影して投影画像を生成するための投影面に設定する処理を実行する。
本開示によれば、関節を含む３次元画像から投影画像を生成するための投影面を適切に設定することができる。
以下、図面を参照して本開示の実施形態について説明する。図１は、本開示の第１の実施形態による関節投影面設定装置を適用した、診断支援システムの概要を示すハードウェア構成図である。図１に示すように、診断支援システムでは、本実施形態による関節投影面設定装置１、３次元画像撮影装置２、および画像保管サーバ３が、ネットワーク４を経由して通信可能な状態で接続されている。
３次元画像撮影装置２は、被検体の診断対象となる部位を撮影することにより、その部位を表す３次元画像を生成する装置である。具体的には、３次元画像撮影装置２は、ＣＴ装置、ＭＲＩ装置、およびＰＥＴ（Positron Emission Tomography）装置等である。３次元画像撮影装置２により生成された３次元画像は画像保管サーバ３に送信され、保存される。なお、本実施形態においては、被検体である患者の診断対象部位は膝関節であり、３次元画像撮影装置２はＭＲＩ装置であり、３次元画像撮影装置２において、被検体の膝関節のＭＲＩ画像が３次元画像として生成されるとする。
画像保管サーバ３は、各種データを保存して管理するコンピュータであり、大容量外部記憶装置およびデータベース管理用ソフトウェアを備えている。画像保管サーバ３は、有線または無線のネットワーク４を介して他の装置と通信を行い、画像データ等を送受信する。具体的には３次元画像撮影装置２で生成された３次元画像等の画像データを含む各種データをネットワーク経由で取得し、大容量外部記憶装置等の記録媒体に保存して管理する。なお、画像データの格納形式およびネットワーク４経由での各装置間の通信は、ＤＩＣＯＭ（Digital Imaging and Communication in Medicine）等のプロトコルに基づいている。
関節投影面設定装置１は、１台のコンピュータにインストールした本開示の関節投影面設定プログラムを実行することにより実現される。コンピュータは、診断を行う医師が直接操作するワークステーションまたはパーソナルコンピュータでもよいし、それらとネットワークを介して接続されたサーバコンピュータでもよい。関節投影面設定プログラムは、ネットワークに接続されたサーバコンピュータの記憶装置、またはネットワークストレージに、外部からアクセス可能な状態で記憶され、要求に応じて医師が使用するコンピュータにダウンロードされ、インストールされる。あるいは、ＤＶＤ（Digital Versatile Disc）またはＣＤ－ＲＯＭ（Compact Disk Read Only Memory）等の記録媒体に記録されて配布され、その記録媒体からコンピュータにインストールされる。
図２は、コンピュータに関節投影面設定プログラムをインストールして実行することにより実現される、本開示の第１の実施形態による関節投影面設定装置の概略構成を示す図である。図２に示すように、関節投影面設定装置１は、標準的なワークステーションの構成として、ＣＰＵ（Central Processing Unit）１１、メモリ１２およびストレージ１３を備えている。また、関節投影面設定装置１には、表示部１４並びにマウスおよびキーボード等の入力部１５が接続されている。
ストレージ１３には、ネットワーク４を経由して画像保管サーバ３から取得した、被検体の３次元画像、および処理に必要な情報を含む各種情報が記憶されている。なお、本実施形態においては、被検体についての膝関節を診断対象部位とする３次元画像Ｖ０がストレージ１３に記憶されている。
また、メモリ１２には、関節投影面設定プログラムが記憶されている。関節投影面設定プログラムは、ＣＰＵ１１に実行させる処理として、関節を含む３次元画像Ｖ０を取得する画像取得処理、関節の関節面を近似する平面を、３次元画像を投影して後述する投影画像を生成するための投影面に設定する面設定処理、３次元画像Ｖ０の関節にある軟骨を投影面に投影して投影画像を生成する投影処理、投影画像上において、軟骨の定量値を導出する定量化処理、定量化により導出された後述する厚さマップを表示部１４に表示する表示制御処理を規定する。
そして、ＣＰＵ１１がプログラムに従いこれらの処理を実行することで、コンピュータは、画像取得部２１、面設定部２２、投影部２３、定量化部２４および表示制御部２５として機能する。
画像取得部２１は、被検体の膝関節の３次元画像Ｖ０を画像保管サーバ３から取得する。なお、３次元画像Ｖ０が既にストレージ１３に記憶されている場合には、画像取得部２１は、ストレージ１３から３次元画像Ｖ０を取得するようにしてもよい。図３は膝関節の３次元画像Ｖ０を示す図である。図３に示すように、３次元画像Ｖ０には、大腿骨３０および脛骨３１が含まれている。なお、図３においては、説明のために膝蓋骨は省略している。大腿骨３０の脛骨３１と面する部分には軟骨３２が、脛骨３１の大腿骨３０と面する部分には軟骨３３が存在する。また、軟骨３２と軟骨３３との間には半月板３４が存在する。本実施形態においては、３次元画像Ｖ０はＭＲＩ画像であり、骨、軟骨、半月板、並びにこれら以外の筋肉および脂肪等領域のそれぞれにおいて、３次元画像Ｖ０における信号値（ボクセル値）の範囲が異なる。画像取得部２１は、信号値に対するしきい値処理により３次元画像Ｖ０から骨領域および軟骨領域を抽出する。具体的には、３次元画像Ｖ０において、骨の信号値となる範囲の領域を骨領域として抽出する。また、３次元画像Ｖ０において、軟骨の信号値となる範囲の領域を軟骨領域として抽出する。骨領域には大腿骨３０および脛骨３１が含まれ、軟骨領域には軟骨３２および３３が含まれる。
なお、本実施形態においては、画像取得部２１が３次元画像Ｖ０から骨領域および軟骨領域を抽出しているが、これに限定されない。３次元画像Ｖ０から骨領域および軟骨領域を抽出する手段を別途設けるようにしてもよい。なお、本実施形態においては、脛骨３１の軟骨３３を軟骨領域として抽出する。また、３次元画像Ｖ０からの骨領域および軟骨領域の抽出は、しきい値処理に限定されない。例えば、３次元画像Ｖ０から骨領域および軟骨領域を抽出するように、ディープラーニング等により機械学習がなされた判別器を用いてもよい。
一方、３次元画像Ｖ０を３次元画像撮影装置２により取得する際、膝を伸ばした状態、または軽く曲げた状態（１０～２０度）により、膝関節が撮影される。そして、３次元画像Ｖ０においては、図３に示すように、大腿骨３０および脛骨３１が上から下に延びる方向がｚ方向に設定され、膝関節を正面から見た場合の奥から手前に向かう方向がｙ方向に設定され、膝関節を正面から見た場合の左から右に向かう方向がｘ方向に設定される。
面設定部２２は、関節の関節面を近似する平面を、３次元画像を投影して後述する投影画像を生成するための投影面に設定する。以下、面設定部２２が行う処理について説明する。まず、面設定部２２は、３次元画像Ｖ０に含まれる脛骨３１を囲む最小の３次元領域を設定する。最小の３次元領域とは、３次元画像Ｖ０を囲む領域のうち、後述する投影画像を生成するための投影面を設定するために必要最小限の領域をいう。図４は３次元領域を示す図である。図４に示すように面設定部２２は、脛骨３１を囲む最小の３次元領域４０を設定する。なお、３次元領域４０の各辺の方向は、３次元画像Ｖ０のｘ方向、ｙ方向およびｚ方向と一致させる。ここで、ｚ方向が本開示の顆間領域設定方向に対応する。
次に面設定部２２は、３次元領域４０の上面４０Ａにおいて、顆間の領域を除いた領域を関節面領域として抽出する。図５は関節面領域の抽出を説明するための図である。図５に示すように、面設定部２２は、上面４０Ａのｘ方向に延びる辺４１の中点を基準として、辺４１の長さに対して予め定められた顆間領域幅を有する領域を顆間の領域４２に設定する。なお、顆間領域幅としては上面４０Ａのｘ方向における１５～２５％、好ましくは２０％とすることができるが、これに限定されない。そして、面設定部２２は、３次元領域４０のうち、顆間の領域４２よりもｘ方向においてマイナス方向に位置する領域を関節面領域４３に設定する。また、面設定部２２は、３次元領域４０のうち、顆間の領域４２よりもｘ方向においてプラス方向に位置する領域を関節面領域４４に設定する。なお、関節面領域４３が脛骨３１の関節における内側顆側の領域に対応し、関節面領域４４が脛骨３１の関節における外側顆側の領域に対応する。
次に、面設定部２２は、関節面領域４３および４４に含まれる関節面を近似する平面を投影面に設定する。具体的には、関節面領域４３および４４に含まれる関節面を構成する各ボクセルからの距離が最小となる平面を最小自乗法（最小二乗法）により投影面として導出する。図６は投影面の導出を説明するための図である。なお、図６においては説明を簡単にするために、脛骨３１のｙ軸に垂直な断面を示している。また、関節面領域４３および４４上の各ボクセルのうちの一部のボクセルのみを黒丸で示し、ボクセルから平面、すなわち投影面４５までの距離を矢印で示している。図６に示すように、面設定部２２は関節面領域４３および４４に含まれる関節面を構成する各ボクセルからの距離の総和が最小となるような平面を投影面４５として導出する。
次に、面設定部２２は、関節面領域４３および４４における各ボクセルから投影面４５までの距離が、予め定められたボクセル除外距離以上となるボクセルを除外して新たな投影面を設定する。図７はボクセルの除外を説明するための図である。なお、図７はｚｙ平面における脛骨３１の断面を示している。本実施形態においては、面設定部２２は、関節面領域４３および４４に含まれる関節面を構成する各ボクセルからの距離が最小となる平面４５を導出している。しかしながら、関節面には、図７に示すボクセルＰ１のように、平面４５からの距離が他のボクセルと比較して大きくなるボクセルが存在する。このようなボクセルを用いて投影面を設定すると、関節面を近似する平面を精度よく導出することができない。
このため、面設定部２２は、関節面領域４３および４４に含まれる関節面を構成する各ボクセルからの投影面４５までの距離のうち、ボクセル除外距離以上となるボクセルを除外する。そして、面設定部２２は、除外したボクセル以外のボクセルを用いて、関節面領域４３および４４上の各ボクセルからの距離の総和が最小となる平面を、新たな投影面４５Ａとして導出する。ここで、ボクセル除外距離としては、例えば、関節面領域４３および４４に含まれる関節面を構成する各ボクセルからの投影面４５までの距離の最大値の５～２０％、好ましくは１０％とすることができるが、これに限定されない。
本実施形態においては、面設定部２２は、関節面領域４３および４４に含まれる関節面を構成する各ボクセルからの新たな投影面４５Ａまでの距離のうち、さらにボクセル除外距離以上となるボクセルを除外して新たな投影面４５Ｂを設定する処理を繰り返して、最終的な投影面４６を導出する。例えば、新たな投影面を設定する処理を２回繰り返すと、上記ボクセル除外距離を１０％とした場合、関節面領域４３および４４のすべてのボクセルのうち、投影面４６の設定に寄与したボクセルの割合は、８１％となる。なお、本実施形態においては、上記ボクセル除外距離を１０％として、新たな投影面を設定する処理を２回繰り返すとするが、これに限定されない。
投影部２３は、３次元画像Ｖ０の関節にある軟骨を投影面４６に直交する方向に投影して投影画像を生成する。すなわち、図８に示すように、脛骨３１の軟骨３３を、投影面４６に直交する方向４７に投影して、投影画像を生成する。本実施形態においては、投影画像の生成の際に、脛骨３１も投影するが、軟骨３３のみを投影して投影画像を生成してもよい。なお、直交とは完全に直交する場合のみならず、１～２度程度のある程度の誤差を持って直交する場合も含む。
この際、投影部２３は、脛骨３１の内側顆の関節面の重心および外側顆の関節面の重心を結ぶ線が、予め定められた重心方向を向くように投影画像を生成する。図９は重心の導出を説明するための図である。投影面４６の導出に際して使用されるボクセルは、上述した新たな投影面を設定する処理を２回繰り返すことにより、関節面領域４３および４４のすべてのボクセルの８１％に低減されている。投影部２３は、関節面領域４３および４４に含まれる関節面のそれぞれにおいて、投影面４６を設定する処理に使用されたボクセルのみを用いて、脛骨３１の内側顆の関節面の重心および外側顆の関節面の重心を導出する。この際、関節面領域４３内のボクセルを用いて内側顆の関節面の重心Ｇ１が導出され、関節面領域４４内のボクセルを用いて外側顆の関節面の重心Ｇ２がそれぞれ導出される。なお、図９においては、重心Ｇ１およびＧ２の導出に使用されるボクセルが存在する領域を破線で囲んで示している。
投影部２３は、重心Ｇ１およびＧ２を結ぶ線が投影画像において水平となるように、脛骨３１の軟骨３３を投影面４６に直交する方向に投影して、投影画像を生成する。図１０は投影画像を示す図である。図１０に示す投影画像５０においては、重心Ｇ１およびＧ２を結ぶ線４８が水平となっている。
なお、脛骨３１の関節面において赤木ラインを導出し、赤木ラインが投影画像の垂直方向を向くように、脛骨３１の軟骨３３を投影面４６に直交する方向に投影して、投影画像を生成してもよい。赤木ラインとは、後十字靱帯付着部と膝蓋腱付着部内側縁を結んだ線である。図１１は、赤木ラインが投影画像の垂直方向を向くように生成した投影画像を示す図である。なお、図１１に示すように、投影画像５０においては、赤木ライン４９が垂直となっている。
さらに投影部２３は、投影画像５０上に関心領域を設定する。本実施形態においては、関節における軟骨下骨領域に対応する領域を関心領域として設定する。図１２は関心領域の設定を説明するための図である。図１２に示す投影画像５０において、内側顆の関節面の軟骨領域５１Ａおよび外側顆の関節面の軟骨領域５１Ｂのそれぞれにハッチングを付与して示している。ここで、軟骨下骨領域とは、脛骨３１の関節において、大腿骨３０の関節と擦り合う領域である。投影画像５０における軟骨領域５１Ａおよび５１Ｂの周辺部は、大腿骨３０の関節とは擦り合わない。このため、投影部２３は、投影画像５０の軟骨領域５１Ａおよび５１Ｂの縁部から予め定められた軟骨下骨領域抽出範囲にある領域を除外した領域を、軟骨下骨領域として抽出し、抽出した軟骨下骨領域を関心領域５２Ａおよび５２Ｂに設定する。一方、脛骨３１においては、関節において軟骨が存在すべき領域を規定する輪郭が、関節面に凸状部分として含まれる。このため、関節面における凸状部分により囲まれる領域を軟骨領域５１Ａおよび５１Ｂと見なして、関心領域５２Ａおよび５２Ｂを設定してもよい。
定量化部２４は、投影画像５０上における軟骨領域５１Ａおよび５１Ｂの定量値を導出する。具体的には、投影画像５０上の関心領域５２Ａおよび５２Ｂにおいて定量値を導出する。本実施形態においては、定量値の導出は軟骨領域５１Ａおよび５１Ｂの関心領域５２Ａおよび５２Ｂ毎に行うが、これに限定されない。関心領域５２Ａおよび５２Ｂの全領域についての定量値を導出してもよい。
まず、定量化部２４は、関心領域５２Ａおよび５２Ｂの面積、および関心領域５２Ａおよび５２Ｂ内における軟骨領域５１Ａおよび５１Ｂの面積を導出する。なお、投影画像５０においては、１画素当たりの面積が分かる。このため、定量化部２４は、関心領域５２Ａおよび５２Ｂならびに軟骨領域５１Ａおよび５１Ｂの画素数をカウントし、カウントした画素数に１画素当たりの面積を乗算することにより、軟骨領域５１Ａおよび５１Ｂならびに関心領域５２Ａおよび５２Ｂの面積を導出する。軟骨領域５１Ａおよび５１Ｂの面積は定量値の１つである。
また、定量化部２４は、軟骨領域５１Ａおよび５１Ｂの厚さを定量値として導出する。ここで、投影画像５０の各画素位置と、３次元画像Ｖ０における脛骨３１における軟骨３３の位置とは対応づけることができる。定量化部２４は、投影画像５０の関心領域５２Ａおよび５２Ｂにおける各画素位置を、３次元画像Ｖ０における脛骨３１における軟骨３３の各画素位置と対応づける。そして、軟骨３３の画素位置における法線方向において、軟骨３３の表面と脛骨３１の骨部の表面との間の画素数を導出する。そして、導出した画素数に１画素当たりの長さを乗算することにより、軟骨領域５１Ａおよび５１Ｂの厚さを導出する。
また、定量化部２４は、軟骨領域５１Ａおよび５１Ｂの体積を定量値として導出する。軟骨領域５１Ａおよび５１Ｂの体積は、軟骨領域５１Ａおよび５１Ｂの面積に厚さを乗算することにより導出することができる。
また、定量化部２４は、軟骨領域５１Ａおよび５１Ｂの面積および厚さから、他の定量値を導出する。具体的には、関心領域５２Ａおよび５２Ｂ内における軟骨領域５１Ａおよび５１Ｂの被覆率、関心領域５２Ａおよび５２Ｂにおける軟骨領域５１Ａおよび５１Ｂの欠損面積、ならびに、関心領域５２Ａおよび５２Ｂ内の各位置における軟骨領域５１Ａおよび５１Ｂの厚さの代表値を定量値として導出する。
関心領域５２Ａおよび５２Ｂ内における軟骨領域５１Ａおよび５１Ｂの被覆率は、（関心領域５２Ａおよび５２Ｂ内における軟骨領域５１Ａおよび５１Ｂの面積）／（関心領域５２Ａおよび５２Ｂの面積）により導出する。
関心領域５２Ａおよび５２Ｂにおける軟骨領域５１Ａおよび５１Ｂの欠損面積として、定量化部２４は、軟骨領域５１Ａおよび５１Ｂが存在しない領域の面積を導出する。例えば、図１２に示す領域５３Ａは軟骨が存在しない領域であるため、軟骨領域５１Ａにおける軟骨が存在しない領域５３Ａの面積が軟骨領域５１Ａの欠損面積となる。
関心領域５２Ａおよび５２Ｂ内の各位置における軟骨領域５１Ａおよび５１Ｂの厚さの代表値としては、定量化部２４は、関心領域５２Ａおよび５２Ｂ内の軟骨領域５１Ａおよび５１Ｂの厚さの平均値、中間値、最小値または最大値等を導出する。
また、定量化部２４は、関心領域５２Ａおよび５２Ｂ内の各位置における軟骨領域５１Ａおよび５１Ｂの厚さから、厚さマップを生成する。図１３は厚さマップを示す図である。図１３に示すように厚さマップＭ０は、関心領域５２Ａおよび５２Ｂにおける軟骨領域５１Ａおよび５１Ｂの厚さの分布が６段階の色により示されている。厚さマップＭ０においては、色が濃いほど軟骨領域５１Ａおよび５１Ｂが薄いことを示している。なお、図１３においては色の相違をハッチングの相違により示している。また、厚さマップＭ０には、色と厚さとの関係を示すレファレンス６１が含まれる。レファレンス６１を参照することにより、厚さマップＭ０において、関心領域５２Ａおよび５２Ｂにおける軟骨領域５１Ａおよび５１Ｂの厚さの分布を視覚的に容易に認識することができる。
なお、定量化部２４は、投影画像５０において、軟骨領域５１Ａおよび５１Ｂの厚さが予め定められた定量値導出対象厚さ以上となる画素位置のみを用いて定量値を導出してもよい。例えば、投影画像５０において、軟骨領域５１Ａおよび５１Ｂの厚さが０．５ｍｍ以上となる画素位置のみを用いて定量値を導出してもよい。この場合、軟骨の厚さが０．５ｍｍ未満となる画素位置については、定量値の導出から除外する。これにより、厚さが薄く、軟骨として機能していない領域を定量値の導出から外すことができる。なお、欠損面積を定量値とした場合、軟骨の厚さが０．５ｍｍ未満となる画素位置を、軟骨が欠損している画素位置として欠損面積を導出すればよい。
なお、導出された定量値は、患者名、撮影日時、関心領域５２Ａおよび５２Ｂの位置、ならびに投影画像５０等の情報と併せて、３次元画像Ｖ０と対応づけられて画像保管サーバ３に送信され、保存される。
表示制御部２５は、厚さマップＭ０を表示部１４に表示する。
次いで、第１の実施形態において行われる処理について説明する。図１４は第１の実施形態において行われる処理を示すフローチャートである。まず、画像取得部２１が、３次元画像Ｖ０を取得し（ステップＳＴ１）、３次元画像Ｖ０から骨領域および軟骨領域を抽出する（ステップＳＴ２）。次いで、面設定部２２が、脛骨３１の関節の関節面を近似する平面を、３次元画像Ｖ０を投影して投影画像を生成するための投影面４６に設定する（ステップＳＴ３）。
次いで、投影部２３が、３次元画像Ｖ０の関節にある軟骨を投影面４６に直交する方向に投影して投影画像５０を生成する（ステップＳＴ４）。そして、定量化部２４が、投影画像５０上において、軟骨領域５１Ａおよび５１Ｂの定量値を導出し（ステップＳＴ５）、導出された定量値から厚さマップＭ０を生成する（ステップＳＴ６）。さらに、表示制御部２５が厚さマップＭ０を表示部１４に表示し（ステップＳＴ７）、処理を終了する。
このように、本実施形態においては、関節の関節面を近似する平面を、３次元画像Ｖ０を投影して投影画像を生成するための投影面４６に設定するようにした。このため、とくに、脛骨３１のように比較的平坦な関節の関節面に関して、軟骨の欠損部分に軟骨が存在するように投影されることを防止することができる。したがって、本実施形態によれば、関節を含む３次元画像Ｖ０から投影画像５０を生成するための投影面を適切に設定することができる。
次いで、本開示の第２の実施形態について説明する。図１５は本開示の第２の実施形態による関節投影面設定装置の概略構成を示す図である。なお、図１５において図２と同一の構成については同一の参照番号を付与し、詳細な説明は省略する。図１５に示すように第２の実施形態による関節投影面設定装置１Ａは、関節についての機能画像を投影画像にマッピングしたマッピング画像を生成するマッピング部２６をさらに備えた点が第１の実施形態と異なる。なお、図１５においては、定量化部２４を備えているが、第２の実施形態による関節投影面設定装置１Ａは、定量化部２４を備えなくてもよい。
なお、第２の実施形態においては、画像取得部２１は３次元画像Ｖ０に加えて、膝関節の機能画像Ｆ０を取得する。第２の実施形態においては、機能画像Ｆ０としてＴ２マップ画像を用いる。Ｔ２マップ画像はＭＲＩ画像の一種であり、膝関節における水との相関性を信号値として表した画像である。Ｔ２マップ画像においては、例えば、水との相関性が高いほど各画素における信号値が大きくなる。
第２の実施形態においては、機能画像Ｆ０に対しても、第１の実施形態と同様に、面設定部２２が、脛骨３１の関節の関節面を近似する平面を、機能画像Ｆ０を投影して投影画像を生成するための投影面４６に設定する。そして、投影部２３が、機能画像Ｆ０の関節にある軟骨を投影面４６に直交する方向に投影して投影画像を生成する。機能画像Ｆ０の投影画像を機能投影画像７０と称する。ここで、脛骨３１における軟骨３３は厚さを有するため、機能画像Ｆ０においては、軟骨３３の厚さ方向において信号値が異なる。第２の実施形態においては、投影部２３は、軟骨３３の厚さ方向に中間面を設定し、中間面における機能画像Ｆ０の信号値を展開して機能投影画像７０を生成する。なお、軟骨３３の厚さ方向における信号値の平均値、最大値、最小値および中間値等の代表値を用いて機能投影画像７０を生成してもよい。
ここで、３次元画像Ｖ０と機能画像Ｆ０とでは、含まれる部位は同一の膝関節であるが、撮影時間が異なる。このため、マッピング部２６は、投影画像５０と機能投影画像７０との位置合わせを行う。位置合わせの手法としては、剛体位置合わせおよび非剛体位置合わせ等、周知の任意の手法を用いることができる。なお、マッピング部２６は、投影画像５０を機能投影画像７０に一致させるように変形して位置合わせを行ってもよく、機能投影画像７０を投影画像５０に一致させるように変形して位置合わせを行ってもよい。そして、マッピング部２６は、機能投影画像７０を投影画像５０にマッピングしたマッピング画像Ｍ１を生成する。
図１６はマッピング画像Ｍ１を示す図である。図１６に示すように、マッピング画像Ｍ１は、脛骨３１の軟骨３３における水との相関性の大きさが６段階の色により示されている。マッピング画像Ｍ１においては、色が濃いほど水との相関性が低いことを示している。なお、図１６においては色の相違をハッチングの相違により示している。また、マッピング画像Ｍ１には、色と水との相関性との関係を示すレファレンス６２が含まれる。レファレンス６２を参照することにより、軟骨３３の水との相関性を視覚的に容易に認識することができる。なお、図１６に示すマッピング画像Ｍ１においては関心領域５２Ａおよび５２Ｂを設定していないが、第１の実施形態と同様に関心領域５２Ａおよび５２Ｂを設定し、関心領域５２Ａおよび５２Ｂのみにおいてマッピング画像Ｍ１を生成してもよい。
次いで、第２の実施形態において行われる処理について説明する。図１７は第２の実施形態において行われる処理を示すフローチャートである。まず、画像取得部２１が、３次元画像Ｖ０および機能画像Ｆ０を取得し（ステップＳＴ１１）、３次元画像Ｖ０から骨領域および軟骨領域を抽出する（ステップＳＴ１２）。次いで、面設定部２２が、脛骨３１の関節の関節面を近似する平面を、３次元画像Ｖ０を投影して投影画像を生成するための投影面４６に設定する（ステップＳＴ１３）。
次いで、投影部２３が、３次元画像Ｖ０の関節にある軟骨を投影面４６に直交する方向に投影して投影画像５０を生成する（ステップＳＴ１４）。また、投影部２３が、機能画像Ｆ０を投影面４６に直交する方向に投影して機能投影画像７０を生成する（ステップＳＴ１５）。そして、マッピング部２６が、機能投影画像７０を投影画像５０にマッピングしたマッピング画像Ｍ１を生成する（ステップＳＴ１６）。さらに、表示制御部２５がマッピング画像Ｍ１を表示部１４に表示し（ステップＳＴ１７）、処理を終了する。
次いで、本開示の第３の実施形態について説明する。なお、第３の実施形態による関節投影面設定装置の構成は、上記第１および第２の実施形態による関節投影面設定装置のいずれかの構成と同一であり、行われる処理のみが異なるため、ここでは装置についての詳細な説明は省略する。第３の実施形態による関節投影面設定装置においては、投影部２３が生成した投影画像５０を表示制御部２５が表示部１４に表示するようにした点が第１および第２の実施形態と異なる。
次いで、第３の実施形態において行われる処理について説明する。図１８は第３の実施形態において行われる処理を示すフローチャートである。まず、画像取得部２１が、３次元画像Ｖ０を取得し（ステップＳＴ２１）、３次元画像Ｖ０から骨領域および軟骨領域を抽出する（ステップＳＴ２２）。次いで、面設定部２２が、脛骨３１の関節の関節面を近似する平面を、３次元画像Ｖ０を投影して投影画像を生成するための投影面４６に設定する（ステップＳＴ２３）。
次いで、投影部２３が、３次元画像Ｖ０の関節にある軟骨を投影面４６に直交する方向に投影して投影画像５０を生成する（ステップＳＴ２４）。そして、表示制御部２５が投影画像５０を表示部１４に表示し（ステップＳＴ２５）、処理を終了する。
次いで、本開示の第４の実施形態について説明する。上記第１から第３の実施形態においては、２次元の投影画像５０を生成しているが、第４の実施形態においては、３次元の投影画像を生成するようにした点が第１から第３の実施形態と異なる。図１９は、第４の実施形態における３次元の投影画像を示す図である。図１９に示すように、投影画像８０は、第１の実施形態と同様に脛骨３１の軟骨領域８１Ａおよび８１Ｂと、軟骨領域８１Ａおよび８１Ｂの厚さを表す断面図８２Ａおよび８２Ｂを含む。なお，第１の実施形態と同様に、軟骨領域８１Ａおよび８１Ｂにおいて関心領域を設定してもよい。また、断面図８２Ａおよび８２Ｂにおける断面の位置を変更できるようにしてもよい。
このようにして生成した３次元の投影画像８０に対しても、上述したように定量値を導出し、厚さマップＭ０を生成することが可能である。また、第２の実施形態と同様に、３次元の投影画像８０に機能画像Ｆ０を重畳したマッピング画像Ｍ１を生成することも可能である。また、第３の実施形態と同様に、３次元の投影画像８０を、例えばボリュームレンダリング表示等により表示部１４に表示することも可能である。
なお、同一の被検体について、撮影時期が異なる複数の３次元画像を比較して経過観察を行う場合がある。このような場合において、撮影時期が過去の第１の３次元画像Ｖ１において投影画像が生成されている場合、撮影時期が新しい第２の３次元画像Ｖ２においても投影画像を生成して経過観察を行うことが好ましい。なお、第１の３次元画像Ｖ１が本開示の他の３次元画像に対応する。この場合、軟骨は経時によって擦り減ることはあるが、関節の形状が変形することはない。このため、第１の３次元画像Ｖ１について投影画像を生成した際に、投影面４６を表す情報を画像保管サーバ３に保存しておき、第２の３次元画像Ｖ２の投影画像を生成する際には、同一被検体について保存されている投影面の情報を取得し、取得した投影面の情報を用いて投影画像を生成することが好ましい。
これにより、第２の３次元画像Ｖ２の投影画像を生成する際の演算量を低減することができる。また、この場合において、第１の３次元画像Ｖ１の投影画像に関心領域５２Ａおよび５２Ｂが設定されている場合には、第２の３次元画像Ｖ２の投影画像に対して、第１の３次元画像Ｖ１の投影画像と同一の関心領域を設定することが好ましい。これにより、第１の３次元画像Ｖ１の軟骨の状態と第２の３次元画像Ｖ２における軟骨の状態との比較を容易に行うことができる。なお、第１の３次元画像Ｖ１から導出した定量値の導出結果が存在する場合、投影部２３は、第１の３次元画像Ｖ１から定量値を導出した際と同一位置に関心領域を設定してもよい。
また、複数の異なる被検体の３次元画像Ｖ０から投影画像を生成するに際して、同一の投影面を設定するようにしてもよい。この場合、さらに投影画像を生成する際の投影面を複数の異なる被検体の３次元画像Ｖ０において同一とすることが好ましい。これにより、被検体間の軟骨の状態の比較を容易に行うことができる。
また、上記各実施形態においては、関心領域５２Ａおよび５２Ｂの全領域において定量値を導出しているが、これに限定されない。関心領域５２Ａおよび５２Ｂのそれぞれにおいて別個に定量値を導出してもよい。また、図２０に示すように、定量化部２４は、投影画像５０を複数（図２０においては４つ）の領域５０Ａ～５０Ｄに分割することにより関心領域５２Ａおよび５２Ｂを分割し、分割した領域毎に定量値を導出してもよい。
また、上記各実施形態においては、投影画像５０に含まれる関心領域５２Ａおよび５２Ｂにおいて定量値を導出しているが、関心領域５２Ａおよび５２Ｂを設定することなく、投影画像５０に含まれる軟骨領域５１Ａおよび５１Ｂの全体において定量値を導出するようにしてもよい。
また、上記各実施形態においては、軟骨領域５１Ａおよび５１Ｂの面積、関心領域５２Ａおよび５２Ｂにおける軟骨領域５１の面積、軟骨領域５１Ａおよび５１Ｂの厚さ、軟骨領域５１の体積、関心領域５２Ａおよび５２Ｂ内における軟骨領域５１Ａおよび５１Ｂの被覆率、関心領域５２Ａおよび５２Ｂにおける軟骨領域５１Ａおよび５１Ｂの欠損面積、ならびに、関心領域５２Ａおよび５２Ｂ内の各位置における軟骨領域５１Ａおよび５１Ｂの厚さの代表値を定量値として導出している。しかしながら、これらの定量値のうちのいずれか１つ、または任意の組み合わせの定量値を導出するようにしてもよい。
また、上記各実施形態においては、内側顆の関節面および外側顆の関節面の関節面領域４３および４４の双方を用いて投影面を設定しているが、これに限定されない。内側顆の関節面の関節面領域４３のみを用いて投影面を設定してもよい。この場合、設定した投影面を用いて、内側顆の関節面の関節面領域４３のみの投影画像を生成すればよい。一方、外側顆の関節面の関節面領域４４のみを用いて投影面を設定してもよい。この場合、設定した投影面を用いて、外側顆の関節面の関節面領域４４のみの投影画像を生成すればよい。
また，上記実施形態においては、ボクセルを除外しての新たな投影面を設定する処理を繰り返し行っているが、これに限定されない。最初に設定した投影面４５において、関節面領域４３および４４に含まれる関節面を構成する各ボクセルからの投影面４５までの距離のうち、ボクセル除外距離以上となるボクセルを除外して最終的な投影面を設定するようにしてもよい。この場合、ボクセル除外距離としては、ボクセルを除外しての新たな投影面を設定する処理を複数回行う場合と同一であってもよく、ボクセルを除外しての新たな投影面を設定する処理を複数回行う場合よりも多くてもよい。例えば、ボクセル除外距離を最大値の２０％としてもよい。
また、上記実施形態においては、膝関節の軟骨の定量値を導出しているが、これに限定されず、肘関節、股関節、肩関節または椎間関節等の軟骨の定量値を導出する際に本開示を適用できることはもちろんである。
また、上記各実施形態において、例えば、画像取得部２１、面設定部２２、投影部２３、定量化部２４、表示制御部２５およびマッピング部２６といった各種の処理を実行する処理部（Processing Unit）のハードウェア的な構造としては、次に示す各種のプロセッサ（Processor）を用いることができる。上記各種のプロセッサには、上述したように、ソフトウェア（プログラム）を実行して各種の処理部として機能する汎用的なプロセッサであるＣＰＵに加えて、ＦＰＧＡ（Field Programmable Gate Array）等の製造後に回路構成を変更可能なプロセッサであるプログラマブルロジックデバイス（Programmable Logic Device :PLD）、ＡＳＩＣ（Application Specific Integrated Circuit）等の特定の処理を実行させるために専用に設計された回路構成を有するプロセッサである専用電気回路等が含まれる。
１つの処理部は、これらの各種のプロセッサのうちの１つで構成されてもよいし、同種または異種の２つ以上のプロセッサの組み合わせ（例えば、複数のＦＰＧＡの組み合わせまたはＣＰＵとＦＰＧＡとの組み合わせ）で構成されてもよい。また、複数の処理部を１つのプロセッサで構成してもよい。
複数の処理部を１つのプロセッサで構成する例としては、第１に、クライアントおよびサーバ等のコンピュータに代表されるように、１つ以上のＣＰＵとソフトウェアとの組み合わせで１つのプロセッサを構成し、このプロセッサが複数の処理部として機能する形態がある。第２に、システムオンチップ（System On Chip:SoC）等に代表されるように、複数の処理部を含むシステム全体の機能を１つのＩＣ（Integrated Circuit）チップで実現するプロセッサを使用する形態がある。このように、各種の処理部は、ハードウェア的な構造として、上記各種のプロセッサの１つ以上を用いて構成される。
さらに、これらの各種のプロセッサのハードウェア的な構造としては、より具体的には、半導体素子等の回路素子を組み合わせた電気回路（Circuitry）を用いることができる。
１，１Ａ 関節投影面設定装置
２ ３次元画像撮影装置
３ 画像保管サーバ
４ ネットワーク
１１ ＣＰＵ
１２ メモリ
１３ ストレージ
１４ 表示部
１５ 入力部
２１ 画像取得部
２２ 面設定部
２３ 投影部
２４ 定量化部
２５ 表示制御部
２６ マッピング部
３０ 大腿骨
３１ 脛骨
３２，３３ 軟骨
３４ 半月板
４０ ３次元領域
４０Ａ 上面
４１ 辺
４２ 顆間の領域
４３，４４ 関節面領域
４５，４５Ａ，４５Ｂ，４６ 投影面
４７ 投影面に垂直な方向
４８ 重心を結ぶ線
４９ 赤木ライン
５０ 投影画像
５１Ａ，５１Ｂ，８１Ａ，８１Ｂ 軟骨領域
５２Ａ，５２Ｂ 関心領域
５３Ａ 軟骨が存在しない領域
６１，６２ レファレンス
７０ 機能投影画像
８０ ３次元の投影画像
９０ 脛骨
９１ 軟骨
９２ 欠損部分
９３ 投影画像
Ｆ０ 機能画像
Ｇ１，Ｇ２ 重心
Ｍ０ マップ
Ｍ１ マッピング画像
Ｐ１ ボクセル
Ｖ０，Ｖ１，Ｖ２ ３次元画像
Ｘ０ 中心軸
２ ３次元画像撮影装置
３ 画像保管サーバ
４ ネットワーク
１１ ＣＰＵ
１２ メモリ
１３ ストレージ
１４ 表示部
１５ 入力部
２１ 画像取得部
２２ 面設定部
２３ 投影部
２４ 定量化部
２５ 表示制御部
２６ マッピング部
３０ 大腿骨
３１ 脛骨
３２，３３ 軟骨
３４ 半月板
４０ ３次元領域
４０Ａ 上面
４１ 辺
４２ 顆間の領域
４３，４４ 関節面領域
４５，４５Ａ，４５Ｂ，４６ 投影面
４７ 投影面に垂直な方向
４８ 重心を結ぶ線
４９ 赤木ライン
５０ 投影画像
５１Ａ，５１Ｂ，８１Ａ，８１Ｂ 軟骨領域
５２Ａ，５２Ｂ 関心領域
５３Ａ 軟骨が存在しない領域
６１，６２ レファレンス
７０ 機能投影画像
８０ ３次元の投影画像
９０ 脛骨
９１ 軟骨
９２ 欠損部分
９３ 投影画像
Ｆ０ 機能画像
Ｇ１，Ｇ２ 重心
Ｍ０ マップ
Ｍ１ マッピング画像
Ｐ１ ボクセル
Ｖ０，Ｖ１，Ｖ２ ３次元画像
Ｘ０ 中心軸
Claims (29)
- 関節の３次元画像を取得する画像取得部と、
前記関節の関節面を近似する平面を、前記３次元画像を投影して投影画像を生成するための投影面に設定する面設定部であって、前記関節面を構成する複数のボクセルのうち、前記投影面からの距離がボクセル除外距離以上となるボクセルを除外して、新たな投影面を設定する面設定部とを備えた関節投影面設定装置。 - 前記面設定部は、前記ボクセルの除外および前記新たな投影面の設定を複数回繰り返して、前記投影面を設定する請求項１に記載の関節投影面設定装置。
- 前記面設定部は、最小二乗法により前記関節面を近似する平面を前記投影面に設定する請求項１または２に記載の関節投影面設定装置。
- 前記関節は脛骨の関節である請求項１から３のいずれか１項に記載の関節投影面設定装置。
- 前記面設定部は、前記脛骨の関節面を顆間領域設定方向から見た画像において、顆間の領域を除いた領域を関節面領域として抽出し、該関節面領域に含まれる前記関節面を近似する平面を、前記投影面に設定する請求項４に記載の関節投影面設定装置。
- 前記３次元画像の前記関節にある軟骨を前記投影面に直交する方向に投影することにより投影画像を生成する投影部をさらに備えた請求項１から５のいずれか１項に記載の関節投影面設定装置。
- 前記脛骨の内側顆の関節面の重心および外側顆の関節面の重心を結ぶ線が、水平方向を向くように、前記３次元画像の前記関節にある軟骨を前記投影面に直交する方向に投影することにより投影画像を生成する投影部をさらに備えた請求項４または５に記載の関節投影面設定装置。
- 前記投影面に垂直な方向から見た前記脛骨における赤木ラインが、垂直方向を向くように、前記３次元画像の前記関節にある軟骨を前記投影面に直交する方向に投影することにより投影画像を生成する投影部をさらに備えた請求項４または５に記載の関節投影面設定装置。
- 前記投影画像を表示部に表示する表示制御部をさらに備えた請求項６から８のいずれか１項に１に記載の関節投影面設定装置。
- 前記投影画像上において、前記軟骨の定量値を導出する定量化部をさらに備えた請求項６から９のいずれか１項に記載の関節投影面設定装置。
- 前記定量化部は、前記投影画像上の関心領域において、前記定量値を導出する請求項１０に記載の関節投影面設定装置。
- 前記投影部は、前記関節における軟骨が存在すべき領域を規定する輪郭に基づいて前記関心領域を設定する請求項１１に記載の関節投影面設定装置。
- 前記定量化部は、前記関心領域内における前記軟骨の被覆率を、前記定量値として導出する請求項１１または１２に記載の関節投影面設定装置。
- 前記定量化部は、前記関心領域における前記軟骨の欠損面積を、前記定量値として導出する請求項１１から１３のいずれか１項に記載の関節投影面設定装置。
- 前記定量化部は、前記関心領域内の各位置における前記軟骨の厚さの代表値を、前記定量値として導出する請求項１１から１４のいずれか１項に記載の関節投影面設定装置。
- 前記定量化部は、前記関心領域内の各位置における前記軟骨の厚さを、前記定量値として導出する請求項１１から１４のいずれか１項に記載の関節投影面設定装置。
- 前記定量化部は、前記関心領域における前記軟骨の厚さマップを生成する請求項１６に記載の関節投影面設定装置。
- 前記厚さマップを表示部に表示する表示制御部をさらに備えた請求項１７に記載の関節投影面設定装置。
- 前記定量化部は、前記関心領域における前記軟骨の厚さが、定量値導出対象厚さ以上となる領域においてのみ、前記定量値を導出する請求項１１から１８のいずれか１項に記載の関節投影面設定装置。
- 前記定量化部は、前記投影画像上の前記関心領域を分割し、分割により得られる各領域において、前記定量値を導出する請求項１１から１９のいずれか１項に記載の関節投影面設定装置。
- 前記３次元画像を取得した被検体と同一の被検体についての撮影時期が異なる他の３次元画像から導出した定量値の導出結果が存在する場合、前記投影部は、前記他の３次元画像から定量値を導出した際と同一位置に関心領域を設定する請求項１１から２０のいずれか１項に記載の関節投影面設定装置。
- 前記３次元画像を取得した被検体と同一の被検体についての撮影時期が異なる他の３次元画像の投影画像に関心領域が設定されている場合、前記投影部は、前記３次元画像の前記投影画像に対して、前記他の３次元画像の投影画像の関心領域と同一の関心領域を設定する請求項１１から２０のいずれか１項に記載の関節投影面設定装置。
- 前記定量化部は、前記投影画像上の前記軟骨の面積を前記定量値として導出する請求項１０から２２のいずれか１項に記載の関節投影面設定装置。
- 前記定量化部は、前記投影画像上の前記軟骨の体積を前記定量値として導出する請求項１０から２３のいずれか１項に記載の関節投影面設定装置。
- 前記３次元画像を取得した被検体と同一の被検体についての撮影時期が異なる他の３次元画像から生成した他の投影画像が存在する場合、前記面設定部は、前記他の投影画像を生成した投影面を、前記３次元画像を投影して投影画像を生成するための投影面に設定する請求項６から２４のいずれか１項に記載の関節投影面設定装置。
- 前記関節についての機能画像を前記投影画像にマッピングしたマッピング画像を生成するマッピング部をさらに備えた請求項６から２５のいずれか１項に記載の関節投影面設定装置。
- 前記マッピング画像を表示部に表示する表示制御部をさらに備えた請求項２６に記載の関節投影面設定装置。
- 関節の３次元画像を取得し、
前記関節の関節面を近似する平面を、前記３次元画像を投影して投影画像を生成するための投影面に設定するに際し、前記関節面を構成する複数のボクセルのうち、前記投影面からの距離がボクセル除外距離以上となるボクセルを除外して、新たな投影面を設定する関節投影面設定方法。 - 関節の３次元画像を取得する手順と、
前記関節の関節面を近似する平面を、前記３次元画像を投影して投影画像を生成するための投影面に設定するに際し、前記関節面を構成する複数のボクセルのうち、前記投影面からの距離がボクセル除外距離以上となるボクセルを除外して、新たな投影面を設定する手順とをコンピュータに実行させる関節投影面設定プログラム。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2019107116 | 2019-06-07 | ||
| JP2019107116 | 2019-06-07 | ||
| PCT/JP2020/015145 WO2020246124A1 (ja) | 2019-06-07 | 2020-04-02 | 関節投影面設定装置、方法およびプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPWO2020246124A1 JPWO2020246124A1 (ja) | 2020-12-10 |
| JP7242852B2 true JP7242852B2 (ja) | 2023-03-20 |
Family
ID=73651999
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2021524681A Active JP7242852B2 (ja) | 2019-06-07 | 2020-04-02 | 関節投影面設定装置、方法およびプログラム |
Country Status (2)
| Country | Link |
|---|---|
| JP (1) | JP7242852B2 (ja) |
| WO (1) | WO2020246124A1 (ja) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080232661A1 (en) | 2005-08-17 | 2008-09-25 | Koninklijke Philips Electronics, N.V. | Method and Apparatus Featuring Simple Click Style Interactions According To a Clinical Task Workflow |
| US20080312663A1 (en) | 2007-06-15 | 2008-12-18 | Martin Haimerl | Computer-assisted joint analysis using surface projection |
| JP2010000305A (ja) | 2008-06-23 | 2010-01-07 | Aloka Co Ltd | 超音波診断装置 |
| US20100134491A1 (en) | 2007-03-20 | 2010-06-03 | David Borland | Methods, systems, and computer readable media for flexible occlusion rendering |
| JP2011518645A (ja) | 2008-04-29 | 2011-06-30 | オティスメッド コーポレイション | 悪化前の状態を表すと共に関節形成装置の設計及び製造において使用可能であるコンピュータ化骨モデルの生成 |
| JP2013533765A (ja) | 2010-06-16 | 2013-08-29 | エーツー・サージカル | 骨の３ｄ医用画像から幾何学的要素を自動的に判定する方法およびシステム |
| JP2014529314A (ja) | 2011-07-20 | 2014-11-06 | スミス アンド ネフュー インコーポレーテッド | 解剖学的構造へのインプラントの適合を最適化するシステムおよび方法 |
| JP2018042709A (ja) | 2016-09-14 | 2018-03-22 | 富士フイルム株式会社 | 軟骨定量化装置、方法およびプログラム |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7925324B2 (en) * | 2005-11-16 | 2011-04-12 | Brainlab Ag | Measuring the femoral antetorsion angle γ of a human femur in particular on the basis of fluoroscopic images |
-
2020
- 2020-04-02 JP JP2021524681A patent/JP7242852B2/ja active Active
- 2020-04-02 WO PCT/JP2020/015145 patent/WO2020246124A1/ja active Application Filing
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080232661A1 (en) | 2005-08-17 | 2008-09-25 | Koninklijke Philips Electronics, N.V. | Method and Apparatus Featuring Simple Click Style Interactions According To a Clinical Task Workflow |
| US20100134491A1 (en) | 2007-03-20 | 2010-06-03 | David Borland | Methods, systems, and computer readable media for flexible occlusion rendering |
| US20080312663A1 (en) | 2007-06-15 | 2008-12-18 | Martin Haimerl | Computer-assisted joint analysis using surface projection |
| JP2011518645A (ja) | 2008-04-29 | 2011-06-30 | オティスメッド コーポレイション | 悪化前の状態を表すと共に関節形成装置の設計及び製造において使用可能であるコンピュータ化骨モデルの生成 |
| JP2010000305A (ja) | 2008-06-23 | 2010-01-07 | Aloka Co Ltd | 超音波診断装置 |
| JP2013533765A (ja) | 2010-06-16 | 2013-08-29 | エーツー・サージカル | 骨の３ｄ医用画像から幾何学的要素を自動的に判定する方法およびシステム |
| JP2014529314A (ja) | 2011-07-20 | 2014-11-06 | スミス アンド ネフュー インコーポレーテッド | 解剖学的構造へのインプラントの適合を最適化するシステムおよび方法 |
| JP2018042709A (ja) | 2016-09-14 | 2018-03-22 | 富士フイルム株式会社 | 軟骨定量化装置、方法およびプログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2020246124A1 (ja) | 2020-12-10 |
| JPWO2020246124A1 (ja) | 2020-12-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Humbert et al. | 3D-DXA: assessing the femoral shape, the trabecular macrostructure and the cortex in 3D from DXA images | |
| Udupa et al. | A framework for evaluating image segmentation algorithms | |
| US6836557B2 (en) | Method and system for assessment of biomarkers by measurement of response to stimulus | |
| US20200058098A1 (en) | Image processing apparatus, image processing method, and image processing program | |
| US10806392B2 (en) | Cartilage quantification device, cartilage quantification method, and cartilage quantification program | |
| US20040106868A1 (en) | Novel imaging markers in musculoskeletal disease | |
| US20040147830A1 (en) | Method and system for use of biomarkers in diagnostic imaging | |
| Baka et al. | Evaluation of automated statistical shape model based knee kinematics from biplane fluoroscopy | |
| JP2004537357A (ja) | 関節疾患および関節疾患の時間変化を定量的に評価するためのシステムおよび方法 | |
| Al-Saleh et al. | Three-dimensional assessment of temporomandibular joint using MRI-CBCT image registration | |
| JP2004535871A (ja) | 組織構造および組織構造の時間変化を定量化するためのシステムおよび方法 | |
| US8396266B2 (en) | Characterisation of functional medical image scans | |
| JP7242852B2 (ja) | 関節投影面設定装置、方法およびプログラム | |
| Memiş et al. | A novel approach for computerized quantitative image analysis of proximal femur bone shape deformities based on the hip joint symmetry | |
| US20210390764A1 (en) | Joint image unfolding apparatus, joint image unfolding method, and joint image unfolding program | |
| US20060247864A1 (en) | Method and system for assessment of biomarkers by measurement of response to surgical implant | |
| US10580136B2 (en) | Mapping image generation device, method, and program | |
| JP7230240B2 (ja) | 半月板投影面設定装置、方法およびプログラム | |
| JP2023114463A (ja) | 表示装置、方法およびプログラム | |
| Zarychta | Application of fuzzy image concept to medical images matching | |
| Hu et al. | Automatic bone surface restoration for markerless computer-assisted orthopaedic surgery | |
| Morita | Computer-aided Diagnosis Systems Based on Medical Image Registration | |
| Hampali | 3D Shape Reconstruction of Knee Bones from Low Radiation X-ray Images Using Deep Learning | |
| Chen | Automatic segmentation and motion analysis of the knee joint based on MRI and 4DCT images | |
| Sato et al. | Fundamental Technologies for Integration and Pathology in MCA |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20211019 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20221004 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20221122 |
|
| 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: 20230214 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230308 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7242852 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |