(57)【特許請求の範囲】
【請求項１】
下肢透視画像をコンピュータ上で自動解析して膝関節の大腿骨及び脛骨の骨軸を自動抽出する方法であって、
（ａ）下肢透視画像から生体組織の画像領域の切り出しを行うプロセスと、
（ｂ）プロセス(a)で切り出された透視画像から輝度反転変換した画像（以下、吸収画像と称する）を導出し、該吸収画像の輝度分布に基づき膝間隙ラインを抽出するプロセスと、
（ｃ）プロセス(b)で求めた吸収画像から骨組織を強調した画像（以下、強調画像と称する）に変換し、該強調画像の輝度分布と前記膝間隙ラインとに基づき骨境界点を求めるプロセスと、
（ｄ）プロセス(c)で求めた骨境界点と前記膝間隙ラインとに基づき骨軸解析範囲を設定し、該骨軸解析範囲の骨境界点から骨軸を導出するプロセスと、
を有し、プロセス(a)から(d)を順次実行することで膝関節の大腿骨及び脛骨の骨軸が自動抽出され、大腿脛骨外側角が自動算出されることを特徴とする骨軸自動抽出方法。
【請求項２】
プロセス(a)において、下肢透視画像の幅方向をi軸、高さ方向をj軸として、i軸方向及びj軸方向の輝度分布特性と左肢か右肢かの肢情報とを用いて、前記下肢透視画像から生体組織画像以外のノイズ情報を除去した最小領域の切り出しを行うことを特徴とする請求項１に記載の骨軸自動抽出方法。
【請求項３】
請求項２に記載の切り出しにおいてi軸方向の画像切り出し手順は、
（１）各i軸座標でj軸方向に輝度分布の平均輝度値IavbM(i)を求める。
（２）IavbM(i)分布においてi軸の両端から輝度値がある輝度値以上に変化したi軸座標をそれぞれIRngA1とIRngA2とし、これらの点から幅(IRngA2−IRngA1)/δ（δ：予め設定された定数）の長さだけ内側にシフトした内部点として、IRngB1とIRngB2とを求める。
（３）区間(IRngB1, IRngB2)における最小輝度値のi軸座標点IMinItiを求める。
（４）左肢の場合、区間(IRngB1,IMinIti)における最大輝度値の峰点をILine1とし、区間(IMinIti,IRngB2)におけるIMinItiにより近い峰点をILine2とし、ILine1とILine2とをi軸座標の切り出し点とする。右肢の場合、区間(IMinIti,IRngB2)における最大輝度値の峰点をILine2とし、区間(IRngB1,IMinIti)におけるIMinItiにより近い峰点をILine1とし、ILine1とILine2とをi軸座標の切り出し点とする。
ことを特徴とする請求項２に記載の骨軸自動抽出方法。
【請求項４】
請求項２に記載の切り出しにおいてj軸方向の画像切り出し手順は、
（１）各j軸座標でi軸方向に輝度分布の平均輝度値Javb[j]を求める。
（２）Javb[j]分布においてj軸の両端から輝度値が予め決められた設定値以上に変化した座標をそれぞれJLine1とJLine2とし、これらの座標をj軸方向の切り出し点とする。
ことを特徴とする請求項２に記載の骨軸自動抽出方法。
【請求項５】
請求項４に記載のj軸方向の画像切り出し手順において、切り出し区間(JLine1,JLine2)で、Javb[j]の分布をj軸方向に平滑化した分布JavbM[j]を求め、変動分布JavbD[j]=Javb[j]−JavbM[j]を用いてレントゲン写真の継ぎ目線（以下、ラインノイズと称する）の区間を特定し、該ラインノイズ区間の輝度分布を区間外の輝度分布から線形補間して補正することによりラインノイズを消去することを特徴とする請求項４に記載の骨軸自動抽出方法。
【請求項６】
プロセス(b)において、画像の幅方向をi軸、高さ方向をj軸として、
（１）プロセス(a)で切り出された透視画像IBrt0(j,i)を輝度反転変換した吸収画像の輝度分布IBrt1(j,i)を式(1)により求める。
0005397873.tif 000012
ここでIBmaxとIBminはそれぞれ透視画像の最大輝度値と最小輝度値とする。
（２）吸収画像の輝度分布IBrt1(j,i)を用いて式(2)により輝度差分分布DBrt(j)を求める。
0005397873.tif 000013
（３）DBrt(j)のピークのj軸座標を膝間隙ラインとして抽出する。
ことを特徴とする請求項１に記載の骨軸自動抽出方法。
【請求項７】
請求項６に記載の膝間隙ラインの抽出において、前処理としてj軸方向にダイナミックレンジ圧縮（以下、DR圧縮と称する）して輝度変動量を増幅させた輝度分布を導出すること、およびDR圧縮のフィルター処理は、
（１）j軸方向のみにβピクセル（β：予め設定された定数）の移動平均による平滑化を行い、ボケ画像Boke(j,i)を作成する。
（２）ボケ画像Boke(j,i) と吸収画像IBrt1(j,i)とを利用して式(3)による輝度変動画像IBrt2(j,i)を作成する。
0005397873.tif 000014
（３）プロセス(a)で切り出された透視画像の輝度範囲が(IBmin,IBmax)になるようにIBrt2(j,i)を線形補間し、γピクセル（γ：予め設定された定数）のスケールの移動平均による平滑化を行い、j軸方向の圧縮画像を得る。
ことを特徴とする請求項６に記載の骨軸自動抽出方法。
【請求項８】
プロセス(c)において、
(1)
プロセス(b)で求めた吸収画像からダイナミックレンジ圧縮画像処理を用いて骨組織の輝度レベルを強調する強調画像を導出する。
(2)
該強調画像の幅方向をi軸、高さ方向をj軸として、i軸方向の輝度分布に対して予め設定された輝度レベル以上の輝度を持つ分布の塊（以下、島と称する）を検出する。
(3)
前記(2)で検出した島が２個の場合、該２島を選択する。前記(2)で検出した島が３個以上の場合、前記膝間隙ラインに基づき大腿骨か脛骨かを判断し、
(3-1)
大腿骨の場合、島の高さの高い方から２島を選択する。
(3-2)
左肢の脛骨の場合、左側にある２島を優先的に選択する。
(3-3)
右肢の脛骨の場合、右側にある２島を優先的に選択する。
(4)
前記(3)で選択した２島からエッジ画像を利用して島の輝度値がピークとなる位置を求め、該位置を骨境界点として抽出する。
ことを特徴とする請求項１に記載の骨軸自動抽出方法。
【請求項９】
請求項８に記載した強調画像の導出は、
（１）j軸およびi軸の両方向にβピクセル（β：予め設定された定数）の移動平均による平滑化を行い、ボケ画像Boke(j,i)を作成する。
（２）ボケ画像Boke(j,i)と吸収画像IBrt1(j,i)とを利用して式(3)により圧縮画像IBrt2(j,i)を作成する。
0005397873.tif 000015
（３）圧縮画像IBrt2(j,i)から、式(4)〜式(6)に示す強調係数fact(j,i)を導入し、強調画像IBrt3(j,i)を作成する。ここで平均的輝度レベルIBrtLVL(j)はボケ画像Boke(j,i)において、各j軸座標に対して、i軸方向の全区間の輝度分布についてその平均値をMean、その標準偏差値をStd、αを予め設定した定数で求める。
fact(j,i)=Boke(j,i)/IBrtLVL(j)    ・・・式(4) 
IBrt3(j,i)= fact(j,i)＊IBrt2(j,i)    ・・・式(5) 
IBrtLVL(j)=Mean+α×Std      ・・・式(6) 
ことを特徴とする請求項８に記載の骨軸自動抽出方法。
【請求項１０】
請求項８に記載の骨境界点の抽出において、前記強調画像にメキシカンハットフィルターを適用してエッジ画像を導出し、該エッジ画像を利用して前記強調画像の輝度分布を基に骨境界点を求めることを特徴とする請求項８に記載の骨軸自動抽出方法。
【請求項１１】
請求項８に記載の骨境界点の抽出において、骨境界点から高さ方向のj軸方向に回帰直線近似法を用いて、骨境界点に混在するノイズの検出と補正を行うことを特徴とする請求項８に記載の骨軸自動抽出方法。
【請求項１２】
プロセス(d)における骨軸解析範囲の設定手順は、
(1)
骨境界点を求める強調画像IBrt3(j,i)の高さ方向のj軸座標範囲(1〜Jmask1)を決める。（Jmask1：膝間隙からの距離が8cmのj軸座標）
(2)
仮の解析範囲(膝間隙ラインから10cm〜15cmの範囲(J15〜J10))で左右の骨境界点とその中点から回帰直線法で、仮の左右の骨境界線と中線を求める。
(3)
中線と膝間隙ラインとの交点Cとする。中線上にBC=10cmとなるB点、AC=15cmとなるA点を求める。
(4) 
A点とB点で中線と直交する直線を求め、上記(2)で得られた仮の左右の骨境界線との交点(A1,A2,B1,B2)を計算し、右骨境界線に対して範囲(A2〜B2)、左骨境界線に対して範囲(A1〜B1)を骨境界線の解析範囲とする。
ことを特徴とする請求項１に記載の骨軸自動抽出方法。
【請求項１３】
プロセス(d)における骨軸の導出手順は、前記骨軸解析範囲の骨境界点から、回帰直線法を用いて左右の骨境界線を導出し、この左右の骨境界線から等距離にある直線を骨軸の直線として導出することを特徴とする請求項１に記載の骨軸自動抽出方法。
【請求項１４】
請求項１乃至請求項１３のいずれかに記載の骨軸自動抽出方法をコンピュータに実行させることを特徴とする骨軸自動抽出プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、TKA(Total Knee Arthroplasty)術前計画および変形性膝関節症などの診断において、下肢透視画像をコンピュータ上で自動解析して膝関節の大腿骨及び脛骨の骨軸を自動抽出する技術に関するものである。
【背景技術】
【０００２】
膝関節の診療は一般にレントゲン写真（Ｘ線透視画像）を基本に行われている。とくに重要な診療項目の１つとして、膝関節に関与している大腿骨と脛骨の骨軸の画像診断がある。また、これらの骨軸情報は、膝関節の治療としての全置換型人工関節手術においても重要な情報となる。
【０００３】
従来法は以下に述べる目視・手作業方法と対話式方法の２つの方法に分けられる。前者はレントゲン写真データで、後者はＣＴ(computed tomography)のデジタルデータに基づく方法である。レントゲン写真は膝関節の診断では下肢の写真となり、これは３枚のレントゲン写真を縦につなげた写真である（図５の右図）。この下肢写真の正確な読み取りとデジタル化の装置はごく最近開発されたのである。一方、ＣＴデジタルデータはレントゲン写真と同じＸ線による撮影であるが、撮影方式は異なり、３次元データで３次元的な情報を提供する。膝関節の診断はレントゲン写真を基本として、近年ＣＴデータを併用して実施される場合がある。ＣＴデータの画像処理は対話方式で行われている。従来法について概略を以下に述べる。
【０００４】
〔目視・手作業による骨軸抽出法〕
目視方法は医者がレントゲン写真を目視し、次の手順で骨軸の抽出を手作業で行う方法である（図１）。
Ａ１：膝関節の間隙線を目視して膝間隙ラインの検出し、そのラインから上下に10cm〜15cmの骨軸抽出範囲の設定
Ａ２：上方の骨軸抽出範囲で大腿骨の骨領域の検出し、骨部の左右の骨境界線を抽出して、それらの境界線の中線を大腿骨の骨軸とする。
Ａ３：下方の骨軸抽出範囲で、Ａ２と同様な骨境界線検出を行い、その中線を脛骨の骨軸とする。これにより、大腿骨と脛骨の骨軸が抽出され、さらに診療の重要項目の大腿脛骨外側角(FTA)が容易に算出できることになる。
【０００５】
〔対話方式による骨軸抽出法〕
対話式方法は、医者が介在する対話方式によるＣＴのデジタルデータによるコンピュータ処理方法である（非特許文献１）。すなわち、コンピュータが困難な処理を医者が行い、それをもとに骨軸の抽出を行う方法である。医者が対話方法で関与する情報は、膝間隙ラインと抽出したい骨の輝度値の２つである。すなわち、上記の目視・手作業方法と比較して、与えられた膝間隙ラインから目視方法のＡ１ができる。Ａ２とＡ３については与えられた骨の輝度値の情報をもとに、次の画像処理法で処理が可能となる。
Ｂ１：閾値法による骨領域の抽出
骨部のＸ線透視画像の輝度値は患者によっても撮影条件によっても異なる。よって、対話式に与えられた輝度値を必要とする。この値をもとに、同じ輝度値をもつ骨領域の抽出を行うことが出来る。
Ｂ２：エッジ検出法および輪郭連結法
Ｂ１で求められた骨領域を２値化してエッジ法で骨境界の検出を行い、輪郭連結法によって、連続した骨境界を抽出して、これらの境界線の中線から、大腿骨と脛骨の骨軸が得られることになる。
【０００６】
【非特許文献１】M. Fadda, D. Bertelli, S. Martelli, M. Marcacci, P. Dario, C. Paggetti, D. Caramella, D. Trippi： "Computer Assisted Planning for Total Knee Arthroplasty", Lecturer Notes In Computer Science N. 1205 Editors J. Troccaz, E. Grimson, R. Mosges "CVRMed - MRCAS 1997", Springer Verlag, London England UK, 1997, pp. 619-628.
【非特許文献２】船水憲一：「これから始まる医用ディジタル画像の世界」  日本放射線技師会出版(2005)
【非特許文献３】伊藤史人，土井章男，羽石秀昭，鈴木昌彦，徳田正幸，“TKA術前計画支援システムの開発”，電子情報通信学会技術研究報告（医用画像研究会），Vol. 106, No. 509, pp. 129-132, 2007/1/26.3)
【発明の開示】
【発明が解決しようとする課題】
【０００７】
しかしながら、従来の骨軸抽出法では医者が介在することで術前計画の負担が大きくなること、さらに術前計画および変形性膝関節症などの診断については医者の技能レベルに依存することで医者による術前計画および変形性膝関節症などの診断のレベルを一定以上に保つことが難しいという問題点がある。
【０００８】
本発明は、上記問題を解決するため、２次元データである下肢写真のデジタルデータ（ダイコムデータ）から、大腿骨および脛骨の骨軸について、コンピュータ上で自動解析して自動抽出するアルゴリズムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記目的を達成するため、請求項１に記載の骨軸自動抽出方法は、下肢透視画像をコンピュータ上で自動解析して膝関節の大腿骨及び脛骨の骨軸を自動抽出する方法であって、
（ａ）下肢透視画像から生体組織の画像領域の切り出しを行うプロセスと、
（ｂ）プロセス(a)で切り出された透視画像から輝度反転変換した画像（以下、吸収画像と称する）を導出し、該吸収画像の輝度分布に基づき膝間隙ラインを抽出するプロセスと、
（ｃ）プロセス(b)で求めた吸収画像から骨組織を強調した画像（以下、強調画像と称する）に変換し、該強調画像の輝度分布と前記膝間隙ラインとに基づき骨境界点を求めるプロセスと、
（ｄ）プロセス(c)で求めた骨境界点と前記膝間隙ラインとに基づき骨軸解析範囲を設定し、該骨軸解析範囲の骨境界点から骨軸を導出するプロセスと、
を有し、プロセス(a)から(d)を順次実行することで膝関節の大腿骨及び脛骨の骨軸が自動抽出され、大腿脛骨外側角(FTA: Femoro Tibial Angle)が自動算出されることを特徴とする。
【００１０】
請求項２に記載の骨軸自動抽出方法は、プロセス(a)において、下肢透視画像の幅方向をi軸、高さ方向をj軸として、i軸方向及びj軸方向の輝度分布特性と左肢か右肢かの肢情報とを用いて、前記下肢透視画像から生体組織画像以外のノイズ情報を除去した最小領域の切り出しを行うことを特徴とする。
【００１１】
請求項３に記載の骨軸自動抽出方法は、請求項２に記載の切り出しにおいてi軸方向の画像切り出し手順は、
（１）各i軸座標でj軸方向に輝度分布の平均輝度値IavbM(i)を求める。
（２）IavbM(i)分布においてi軸の両端から輝度値がある輝度値以上に変化したi軸座標をそれぞれIRngA1とIRngA2とし、これらの点から幅(IRngA2−IRngA1)/δ（δ：予め設定された定数）の長さだけ内側にシフトした内部点として、IRngB1とIRngB2とを求める。
（３）区間(IRngB1, IRngB2)における最小輝度値のi軸座標点IMinItiを求める。
（４）左肢の場合、区間(IRngB1,IMinIti)における最大輝度値の峰点をILine1とし、区間(IMinIti,IRngB2)におけるIMinItiにより近い峰点をILine2とし、ILine1とILine2とをi軸座標の切り出し点とする。右肢の場合、区間(IMinIti,IRngB2)における最大輝度値の峰点をILine2とし、区間(IRngB1,IMinIti)におけるIMinItiにより近い峰点をILine1とし、ILine1とILine2とをi軸座標の切り出し点とする。
【００１２】
請求項４に記載の骨軸自動抽出方法は、請求項２に記載の切り出しにおいてj軸方向の画像切り出し手順は、
（１）各j軸座標でi軸方向に輝度分布の平均輝度値Javb[j]を求める。
（２）Javb[j]分布においてj軸の両端から輝度値が予め決められた設定値以上に変化した座標をそれぞれJLine1とJLine2とし、これらの座標をj軸方向の切り出し点とする。
ことを特徴とする。
【００１３】
請求項５に記載の骨軸自動抽出方法は、請求項４に記載のj軸方向の画像切り出し手順において、切り出し区間(JLine1,JLine2)で、Javb[j]の分布をj軸方向に平滑化した分布JavbM[j]を求め、変動分布JavbD[j]=Javb[j]−JavbM[j]を用いてレントゲン写真の継ぎ目線（以下、ラインノイズと称する）の区間を特定し、該ラインノイズ区間の輝度分布を区間外の輝度分布から線形補間して補正することによりラインノイズを消去することを特徴とする。
【００１４】
請求項６に記載の骨軸自動抽出方法は、プロセス(b)において、画像の幅方向をi軸、高さ方向をj軸として、
（１）プロセス(a)で切り出された透視画像IBrt0(j,i)を輝度反転変換した吸収画像の輝度分布IBrt1(j,i)を式(1)により求める。
0005397873.tif 000002
ここでIBmaxとIBminはそれぞれ透視画像の最大輝度値と最小輝度値とする。
（２）吸収画像の輝度分布IBrt1(j,i)を用いて式(2)により輝度差分分布DBrt(j)を求める。
0005397873.tif 000003
（３）DBrt(j)のピークのj軸座標を膝間隙ラインとして抽出する。
ことを特徴とする。
【００１５】
請求項７に記載の骨軸自動抽出方法は、請求項６に記載の膝間隙ラインの抽出において、前処理としてj軸方向にダイナミックレンジ圧縮（以下、DR圧縮と称する）して輝度変動量を増幅させた輝度分布を導出すること、およびDR圧縮のフィルター処理は、
（１）j軸方向のみにβピクセル（β：予め設定された定数）の移動平均による平滑化を行い、ボケ画像Boke(j,i)を作成する。
（２）ボケ画像Boke(j,i) と吸収画像IBrt1(j,i)とを利用して式(3)による輝度変動画像IBrt2(j,i)を作成する。
0005397873.tif 000004
（３）プロセス(a)で切り出された透視画像の輝度範囲が(IBmin,IBmax)になるようにIBrt2(j,i)を線形補間し、γピクセル（γ：予め設定された定数）のスケールの移動平均による平滑化を行い、j軸方向の圧縮画像を得る。
ことを特徴とする。ここでβは比較的大きいスケール（例えばβ=50ピクセル）であり、γは比較的小さいスケール（例えばγ=3ピクセル）とする。
【００１６】
請求項８に記載の骨軸自動抽出方法は、プロセス(c)において、
(1) プロセス(b)で求めた吸収画像からダイナミックレンジ圧縮画像処理を用いて骨組織の輝度レベルを強調する強調画像を導出する。
(2) 該強調画像の幅方向をi軸、高さ方向をj軸として、i軸方向の輝度分布に対して予め設定された輝度レベル以上の輝度を持つ分布の塊（以下、島と称する）を検出する。
(3) 前記(2)で検出した島が２個の場合、該２島を選択する。前記(2)で検出した島が３個以上の場合、前記膝間隙ラインに基づき大腿骨か脛骨かを判断し、
(3-1) 大腿骨の場合、島の高さの高い方から２島を選択する。
(3-2) 左肢の脛骨の場合、左側にある２島を優先的に選択する。
(3-3) 右肢の脛骨の場合、右側にある２島を優先的に選択する。
(4) 前記(3)で選択した２島からエッジ画像を利用して島の輝度値がピークとなる位置を求め、該位置を骨境界点として抽出する。
ことを特徴とする。
【００１７】
請求項９に記載の骨軸自動抽出方法は、請求項８に記載した強調画像の導出が、
（１）j軸およびi軸の両方向にβピクセル（β：予め設定された定数）の移動平均による平滑化を行い、ボケ画像Boke(j,i)を作成する。
（２）ボケ画像Boke(j,i)と吸収画像IBrt1(j,i)とを利用して式(3)により圧縮画像IBrt2(j,i)を作成する。
0005397873.tif 000005
（３）圧縮画像IBrt2(j,i)から、式(4)〜式(6)に示す強調係数fact(j,i)を導入し、強調画像IBrt3(j,i)を作成する。ここで平均的輝度レベルIBrtLVL(j)はボケ画像Boke(j,i)において、各j軸座標に対して、i軸方向の全区間の輝度分布についてその平均値をMean、その標準偏差値をStd、αを予め設定した定数で求める。
fact(j,i)=Boke(j,i)/IBrtLVL(j)    ・・・式(4) 
IBrt3(j,i)= fact(j,i)＊IBrt2(j,i)    ・・・式(5) 
IBrtLVL(j)=Mean+α×Std      ・・・式(6) 
ことを特徴とする。
【００１８】
請求項１０に記載の骨軸自動抽出方法は、請求項８に記載の骨境界点の抽出において、前記強調画像にメキシカンハットフィルターを適用してエッジ画像を導出し、該エッジ画像を利用して前記強調画像の輝度分布を基に骨境界点を求めることを特徴とする。
【００１９】
請求項１１に記載の骨軸自動抽出方法は、請求項８に記載の骨境界点の抽出において、骨境界点から高さ方向のj軸方向に回帰直線近似法を用いて、骨境界点に混在するノイズの検出と補正を行うことを特徴とする。
【００２０】
請求項１２に記載の骨軸自動抽出方法は、プロセス(d)における骨軸解析範囲の設定手順が、
(1)
骨境界点を求める強調画像IBrt3(j,i)の高さ方向のj軸座標範囲(1〜Jmask1)を決める。（Jmask1：膝間隙からの距離が8cmのj軸座標）
(2)
仮の解析範囲(膝間隙ラインから10cm〜15cmの範囲(J15〜J10))で左右の骨境界点とその中点から回帰直線法で、仮の左右の骨境界線と中線を求める。
(3)
中線と膝間隙ラインとの交点Cとする。中線上にBC=10cmとなるB点、AC=15cmとなるA点を求める。
(4) 
A点とB点で中線と直交する直線を求め、上記(2)で得られた仮の左右の骨境界線との交点(A1,A2,B1,B2)を計算し、右骨境界線に対して範囲(A2〜B2)、左骨境界線に対して範囲(A1〜B1)を骨境界線の解析範囲とする。
ことを特徴とする。
【００２１】
請求項１３に記載の骨軸自動抽出方法は、プロセス(d)における骨軸の導出手順が、前記骨軸解析範囲の骨境界点から、回帰直線法を用いて左右の骨境界線を導出し、この左右の骨境界線から等距離にある直線を骨軸の直線として導出することを特徴とする。
【００２２】
請求項１４に記載の骨軸自動抽出プログラムは、請求項１乃至請求項１３のいずれかに記載の骨軸自動抽出方法をコンピュータに実行させることを特徴とする。このプログラムがコンピュータにより実行されれば、請求項１乃至請求項１３のいずれかに係る方法におけると基本的に同じ原理に従い、同様な作用効果が実現され得る。
【発明の効果】
【００２３】
請求項１に係る発明によれば、膝間隙と大腿骨及び脛骨の骨軸の抽出をコンピュータ上で自動解析して行うことで、医師の術前計画における負担を軽減するとともに、これまで技能に頼ってきた外科医による術前計画および変形性膝関節症などの診断のレベルを一定以上に保つ効果がある。
さらに骨境界の検出アルゴリズムに、骨組織に反映した輝度分布に基づく骨境界点抽出の解析法を用いることで、単なる輝度値の閾値法に基づく従来法に比べて、精度の信頼性を高める効果がある。
【００２４】
請求項２に係る発明によれば、下肢透視画像の輝度分布特性を利用して生体組織の画像領域だけを含む最小領域の切り出しを行うことが、コンピュータ上で自動解析して行うことが可能となり、さらに人工的なノイズを除去することが可能となる。ここで人工的なノイズには、例えば次のものがある。
（１）レントゲン写真の外枠  
（２）３枚の写真の継ぎ目線と留め金  
（３）スケール表示
（４）マーク
【００２５】
請求項３に係る発明によれば、i軸方向の画像切り出し手順を実施することで、解析効率の向上ができ、かつ人工的なノイズ（例えば、両端の非生体画像、スケール表示）を除去することが可能となる。
【００２６】
請求項４に係る発明によれば、j軸方向の画像切り出し手順を実施することで、解析効率の向上ができ、かつ人工的なノイズ（例えば、両端の非生体画像）を除去することが可能となる。
【００２７】
請求項５に係る発明によれば、請求項４に係る発明の効果に加え、人工的なノイズ（例えば、レントゲン写真の継ぎ目線）を除去することが可能となる。
【００２８】
請求項６に係る発明によれば、下肢透視画像を輝度反転変換した画像(吸収画像)を基に輝度分布特性を利用して膝間隙ラインを求めることが、コンピュータ上で自動解析して行うことが可能となり、膝間隙ライン抽出の精度が一定レベル以上に確保されるとともに、医師の負担を軽減する効果がある。
【００２９】
請求項７に係る発明によれば、請求項６に係る発明の効果に加え、前処理としてj軸方向にダイナミックレンジ圧縮して輝度変動量を増幅させた輝度分布を導出することで、膝間隙ライン抽出の精度を向上させる効果がある。
【００３０】
請求項８に係る発明によれば、骨境界点の抽出がコンピュータ上で自動解析して行うことが可能となり、骨境界点の抽出の精度が一定レベル以上に確保されるとともに、医師の負担を軽減する効果がある。
【００３１】
請求項９に係る発明によれば、請求項８に係る発明の効果に加え、ダイナミックレンジ圧縮画像処理を用いて解析対象の骨部（大腿骨と脛骨）の輝度レベルを強調する強調画像を導出することで、骨境界点の抽出の精度を向上させる効果がある。
【００３２】
請求項１０に係る発明によれば、請求項８に係る発明の効果に加え、強調画像にメキシカンハットフィルターを適用したエッジ画像を利用することで、骨境界点の抽出の精度をさらに向上させる効果がある。
【００３３】
請求項１１に係る発明によれば、請求項８に係る発明の効果に加え、骨境界点に混在するノイズの検出と補正を行うことで、骨境界点の抽出の精度をさらに向上させる効果がある。
【００３４】
請求項１２に係る発明によれば、骨軸解析範囲の設定がコンピュータ上で自動解析して行うことが可能となり、骨軸解析範囲の設定が正確になるとともに、医師の負担を軽減する効果がある。
【００３５】
請求項１３に係る発明によれば、骨軸の導出がコンピュータ上で自動解析して行うことが可能となり、骨軸の導出の精度が一定レベル以上に確保されるとともに、医師の負担を軽減する効果がある。
【００３６】
請求項１４に係る発明によれば、請求項１乃至請求項１３に係る発明の効果と同様な作用効果が実現され得る。
【発明を実施するための最良の形態】
【００３７】
次に、本発明の実施の形態に係る骨軸自動抽出方法について図面に基づいて説明する。なお、この実施の形態により本発明が限定されるものではない。
【００３８】
図２は、本発明の実施の形態に係る骨軸自動抽出方法の全体を示すフローチャート図である。本発明の骨軸自動抽出方法は、Ｘ線透視下肢画像（ダイコムデータ）がコンピュータに入力され、右肢か左肢かの情報を与えれば、自動的に大腿骨と脛骨の骨軸が抽出される方法である。図２に示すように、骨軸自動抽出方法は、
C1：下肢透視画像から生体組織の画像領域の切り出しを行うプロセスと、
C2：プロセスC1で切り出された透視画像から輝度反転変換した画像（以下、吸収画像と称する）を導出し、該吸収画像の輝度分布に基づき膝間隙ラインを抽出するプロセスと、
C3：プロセスC2で求めた吸収画像から骨組織を強調した画像（以下、強調画像と称する）に変換し、該強調画像の輝度分布と前記膝間隙ラインとに基づき骨境界点を求めるプロセスと、
C4：プロセスC3で求めた骨境界点と前記膝間隙ラインとに基づき骨軸解析範囲を設定し、該骨軸解析範囲の骨境界点から骨軸を導出するプロセスと、
を有し、コンピュータ上でプロセスC1からC4を順次実行することで膝関節の大腿骨及び脛骨の骨軸が自動抽出され、大腿脛骨外側角(FTA)が自動算出される。以下に、C1からC4の各プロセスについて説明する。
【００３９】
[C1：下肢解析画像領域の切り出しプロセス]
プロセスC1は、生体組織画像以外のノイズ情報を除去するために、膝関節周辺の下肢透視画像の幅方向をi軸、高さ方向をj軸（画像の座標系を図４に示す）として、i軸方向及びj軸方向の輝度分布特性と左肢か右肢かの肢情報とを用いて、生体組織の画像領域だけを含む最小領域の切り出しを行う。図３に切り出しアルゴリズムを示すフローチャート図を示す。本切り出し法は、処理効率の向上のための解析画像サイズの最小化ばかりでなく、後述する骨境界抽出法で支障をもたらす人工的なノイズを除去することにある。人工ノイズは、図５(右図)からわかるように、例えば次のものがある。
（１）レントゲン写真の外枠  
（２）３枚の写真の継ぎ目線と留め金  
（３）スケール表示
（４）マーク
【００４０】
[C11：i軸方向の画像切り出し方法] 
次にi軸方向の画像切り出し方法(C11)について説明する。図５の下肢透視画像からわかるように、骨のようなＸ線の吸収率が高いところは、低い輝度値になるために、黒い画像となる。生体組織を１つの塊として分離するために、画像切り出しは初めにi軸方向から行うのが合理的と言える。よって、使用する輝度の平均分布は、各i軸座標でj軸方向の輝度分布の平均値を求め、かつi軸方向に平滑化した分布（図６）とした。また、i軸方向の切り出しを確実にするために、右肢か左肢かの情報を必要とする。解析手順は次のようになる。
1)各i軸座標において、j軸方向に輝度分布の平均輝度値IavbM(i)を求める。
2)IavbM(i)分布において、i軸の両端から輝度値がある輝度値以上に変化したi軸座標をそれぞれIRngA1とIRngA2とし、これらの点から、写真の外枠の除去のために、幅(IRngA2−IRngA1)/ δ（例えばδ=5）の長さだけ、内側にシフトした内部点として、IRngB1とIRngB2とを求める。
3)
区間(IRngB1,
IRngB2)における最小輝度値のi軸座標点(IMinIti)を求める。
4)
左肢の場合、区間(IRngB1,IMinIti)における最大輝度値の峰点をILine1とし、区間(IMinIti,IRngB2)におけるIMinItiにより近い峰点をILine2とし、ILine1とILine2とをi軸座標の切り出し点とする。右肢の場合、区間(IMinIti,IRngB2)における最大輝度値の峰点をILine2とし、区間(IRngB1,IMinIti)におけるIMinItiにより近い峰点をILine1とし、ILine1とILine2とをi軸座標の切り出し点とする。よって画像のi軸方向の大きさはWidthI=Iline2−ILine1となる。
上記4)において人工的なノイズのうちスケール表示を除去するために、右肢か左肢かの情報を導入した。すなわち、スケール表示は通常調べたい肢の付近の外側に設置される。右肢か左肢かの情報から、4)における２つの区間のうち、スケールがある区間を特定することができる。その区間においては、選択する峰は最大輝度値の峰より、IMinIti点により近い峰点とした。ここで峰点の定義はその点を挟む両側の２点、すなわち周りの４点より大きな輝度値を有する点である。
【００４１】
[C12：j軸方向の画像切り出し方法]
次にj軸方向の画像切り出し方法(C12)について説明する。j軸方向の切り出しはi軸方向の画像切り出し方法(C11)で処理された画像をもとに、次の２つの目的のために行う。
（１）両端の非生体画像領域の除去
（２）レントゲン写真の継ぎ目線（これをラインノイズと呼ぶ）の輝度補正
【００４２】
上記（１）は画像処理の効率化のために、C11と同様な方法で次の手順で行う。
1)各j軸座標において、i軸方向に輝度分布の平均輝度値Javb[j]（図７）を求める。
2) Javb[j]分布においてj軸の両端から輝度値が予め決められた設定値以上に変化した座標をそれぞれJLine1とJLine2とし、これらの座標をj軸方向の切り出し点とする。よって画像のj軸方向の大きさはHeightJ=JLine2−JLine1となる。
【００４３】
上記（２）は切り出し区間(JLine1,JLine2)で、Javb[j]の分布をj軸方向に平滑化した分布JavbM[j]を求め、変動分布JavbD[j]=Javb[j]−JavbM[j]で解析を行う。JavbD[j]は図８に示したように、２ヶ所のラインノイズで最大値をとる（図８のa点とb点）。骨軸解析で問題となるラインノイズはj軸座標の後半にある脛骨を横切るもの（図８のb点）で、消去する必要がある。消去法は次のようにした。ピークの輝度値の1/10の輝度値をもつj軸座標JpstとJpedを、ピーク座標（b点）から両側にそれぞれ求める。ラインノイズである区間(Jpst,Jped)の輝度分布を、区間外の輝度分布から線形補間して補正する。
【００４４】
以上で説明したC1の切り出し法によって、図５(左図)に示すように、下肢透視画像（Height：2140ピクセル、Width：1760ピクセル）から下肢解析画像(HeightJ：1935ピクセル、WidthI：351ピクセル)が得られることになる。
【００４５】
[C2：膝間隙ライン抽出プロセス]
次に膝間隙ラインを抽出するプロセスC2について説明する。図９に膝間隙ラインの抽出アルゴリズムを示すフローチャート図を示す。プロセスC2は、C1で切り出された透視画像をもとに、膝間隙ラインの抽出を行う。
(1) 透視画像（輝度値IBrt0(j,i)）を式(1)による輝度反転変換を行う。変換された画像（輝度値IBrt1(j,i)）を吸収画像と呼ぶ。この変換によって、Ｘ線の吸収率が高い骨部は、高い輝度を有することになる。
0005397873.tif 000006
j=1,…,M   i=1,…,N
ここでIBmaxとIBminはそれぞれ透視画像の最大輝度値と最小輝度値である。ここで求める膝間隙ラインは間隙の平均的なj軸座標Jrgkとする。Jrgkの抽出法は次の通りである。
(2) 式(2)により吸収画像の輝度分布IBrt1(j,i)からj軸方向に１階差分の絶対値を求め、かつi軸方向に平均化を行い、輝度差分分布DBrt(j)を求める。
0005397873.tif 000007
(3) DBrt(j)のピークのj軸座標が求める膝間隙ラインとなる。
【００４６】
次に、膝間隙ライン抽出を確実にするために、前処理として吸収画像IBrt1(j,i)に対してj軸方向にダイナミックレンジ圧縮（DR圧縮と呼び、ここではDRPressYフィルターと呼ぶ）処理を導入した。DRPressYフィルター処理は次の内容である。
(1)
j軸方向にβピクセル（例えばβ=50ピクセル、約2.5cm程度）の移動平均による平滑化を行い、ボケ画像Boke(j,i)を作成する。
(2)
ボケ画像Boke(j,i)を利用して式(3)による輝度変動画像IBrt2(j,i)を作る。
0005397873.tif 000008
(3)
透視画像の輝度範囲が(IBmin,IBmax)になるようにIBrt2(j,i)を線形補間し、γピクセル（例えばγ=3ピクセル）のスケールの移動平均による平滑化を行い、j軸方向の圧縮画像を得る。
図10に｜勾配｜分布、DBrt(j)に対してDRPressYフィルター有無での効果を示した。この図から、次のようなDRPressYフィルターの効果が示された。このフィルターは平均値を含む比較的大きなスケールの変動を除去するために、得られた勾配分布は大きさと形において、より安定した分布となり、ピーク位置の抽出が容易となることである。
【００４７】
以上で説明したC2のプロセスにより、図11で示される骨軸解析用画像（j軸：HeightJJ=801ピクセルi軸：WidthI=351ピクセル）の切り出しができる。この画像の幅方向はC1と同じであるが、高さ方向は膝間隙Jrgkを中心に、上下方向に20cmの範囲となる。この高さ方向の範囲は骨軸解析範囲が、通常膝間隙から上下に10cm〜15cmの範囲を十分にカバーすることから決定した範囲である。
【００４８】
[C3：骨境界点の抽出プロセス]
次に骨境界点を抽出するプロセスC3（骨境界抽出法）について説明する。図12に骨境界点の抽出および骨軸解析のアルゴリズムを示すフローチャート図を示す。プロセスC3は、
(1) プロセス(b)で求めた吸収画像（図14の左図）からダイナミックレンジ圧縮画像処理を用いて骨組織の輝度レベルを強調する強調画像を導出する（C31）。
(2) 該強調画像の幅方向をi軸、高さ方向をj軸として、i軸方向の輝度分布に対して予め設定された輝度レベル以上の輝度を持つ分布の塊（以下、島と称する）を検出する。
(3) 検出した島が２個の場合、該２島を選択する。検出した島が３個以上の場合、前記膝間隙ラインに基づき大腿骨か脛骨かを判断し、
(3-1) 大腿骨の場合、島の高さの高い方から２島を選択する。
(3-2) 左肢の脛骨の場合、左側にある２島を優先的に選択する。
(3-3) 右肢の脛骨の場合、右側にある２島を優先的に選択する。
(4) 前記(3)で選択した２島からエッジ画像を利用して島の輝度値がピークとなる位置を求め、該位置を骨境界点として抽出する（C32：横断輝度分布法）。
【００４９】
[C31：ダイナミックレンジ圧縮強調画像処理] 
次にダイナミックレンジ圧縮強調画像処理(C31)について説明する。求める骨軸は膝間隙から上下に骨軸解析範囲にある、大腿骨と脛骨の骨幹の骨軸である。解析したい骨幹は図13から示されるように、中心部に髄腔という空洞を持ち、幹壁として硬い緻密質（骨）を有する。開発したダイナミックレンジ圧縮強調画像処理は、このような骨組織を反映させる画像処理である。具体的な解析手順は次の通りである。以下で示す画像は吸収画像と輝度範囲が同一になるように規格化される。
（１）ボケ画像Boke(j,i)（図14の中図）の作成
j軸およびi軸の両方の方向にβピクセル（β：予め設定された定数）の移動平均による平滑化を行い、ボケ画像Boke(j,i)を作成する。
（２）j軸およびi軸方向の圧縮画像IBrt2(j,i)（図14の右図）の作成
ボケ画像Boke(j,i) と吸収画像IBrt1(j,i)とを利用して式(3)により圧縮画像IBrt2(j,i)を作成する。
0005397873.tif 000009
得られた圧縮画像は図13との比較から、骨組織を反映した画像と言える。すなわち、緻密質が最も高い輝度値で、髄腔が低い輝度値の画像である。
（３）強調画像IBrt3(j,i)（図15の左図）の作成
上記の圧縮画像は図14の右図から分かるように、解析する大腿骨と脛骨以外にも輝度の高い生体組織があるために、骨境界点の抽出に対して不十分と言える。よって、解析する骨部を増幅させ、反対にそれ以外では減衰させる強調係数を各j軸座標に設定することにした。このために式(4)〜式(6)に示す強調係数fact(j,i)を導入する。
fact(j,i)=Boke(j,i)/IBrtLVL(j)    ・・・式(4) 
IBrt3(j,i)= fact(j,i)*IBrt2(j,i)     ・・・式(5) 
IBrtLVL(j)=Mean+α×Std      ・・・式(6) 
強調係数fact(j,i)は、ボケ画像Boke(j,i)に対して、各j軸座標の平均的輝度レベルIBrtLVL(j)で規格化することで求めた。ここでIBrtLVL(j)はボケ画像の、各j軸座標に対して、i軸方向の全区間（1, WidthI）の輝度分布の平均値をMean、その標準偏差をStd、αを予め設定された定数で式(6)で与えられる。
式(4)より、強調係数fact(j,i)は、高い輝度骨部に対して１以上で増幅され、逆に低い輝度の生体組織では１以下で減衰される。従って、強調画像は式(5)で得られ、圧縮画像に比べ、解析する大腿骨と脛骨のみがより高い輝度の画像となる（図15の左図）。このことは代表的な横断輝度分布（図16と図17）からも明白に示される。
【００５０】
強調画像の特性について、輝度の頻度分布からも検討する。図18は吸収画像、圧縮画像そして強調画像の輝度の頻度分布である。この図から、次のことが示される。
I.使用したＤＲ圧縮処理は、輝度レベルIBrtLVL(j)付近の輝度の頻度を集中させる。
II.強調処理は、上記の輝度レベルで規格化されたボケ画像で作られた強調係数を掛けることで、さらに骨部の輝度の集中化を促進する。
上記IIについては、次のピーク強調画像（図15の中図）で明らかに示される。すなわち、ピーク強調画像は、図18のピーク頻度のみの輝度値で作成した２値画像で、骨部が明確に存在することで示される。ピーク強調画像の作成法は、図18の強調画像において、ピークの輝度値の範囲を440〜450とすると、440以下の輝度値は440とし、450以上の輝度値は450とする輝度変換するものである。
【００５１】
[C32：横断輝度分布法] 
次に横断輝度分布法(C32)について説明する。この方法は骨軸解析範囲において、強調画像の横断輝度分布をもとに、エッジ画像を利用して骨境界点を求める。エッジ画像は、図19で示した9×9行列で表されるメキシカンハットフィルター（非特許文献(2)）を強調画像に適用して作成し、かつ式(4)のfact(j,i)の処理（式(5)）をした（図15の右図）。そして、骨境界点の回帰直線補正により、骨境界点に混在するノイズの検出と補正を行う。
【００５２】
図20の上図は大腿骨における骨軸解析範囲の任意のj軸座標（図ではj=160）における、強調画像とエッジ画像の横断輝度分布を示す。横断輝度分布は基準輝度レベルLevel(j)からの変動輝度分布とする。Level(j)は式(6)で計算される。解析手順は次の通りである。
1)強調画像において、正の輝度値を有する分布の塊（今後島と呼ぶ）を検出し、島を高さと位置（島のピーク輝度点）でパラメータ化し、島の高さの高い方から2島を選択する。これらの島が図13の緻密質に対応する。図20に従って、島情報（i軸座標であるAk、Bk、Ck：k=1,2）を決める。中心にある髄腔に対して、添字１が左側で、添字２が右側を示す。
2)エッジ画像は強調画像で示した島、すなわち緻密質の境界点を与える。図20から、ピーク点Dkが骨外境界点に、ゼロ点Ekが骨境界点に対応することになる。
3) 脛骨は図20の下図で示すように大腿骨と異なって、近傍に輝度の高い腓骨を伴う。このために、右肢か左肢かの情報をもとに島の位置情報を利用して、脛骨の島のみを選択するようにした。具体的には、強調画像において、3つ以上の島が存在するとき、右肢の場合腓骨は必ず脛骨の左側にあるために、右側にある島を優先的に選択する。左肢の場合逆に腓骨は右側になるので、左側の島が優先的に選択すればよい。
【００５３】
次に骨境界点の回帰直線補正について説明する。骨軸は骨軸解析範囲で直線と見なされ、骨境界線も同様に直線と見なされてきた。よって、骨境界線はその範囲で、骨境界点（図20のE1とE2）を縦断方向（j軸方向）に回帰直線近似法で求められた直線で正確に求められることになる。この性質を利用して、骨境界点に混在するノイズの検出と補正を行った。解析手順は次の通りである。
1)骨境界点に回帰直線近似を行い、標準偏差以上外れたデータはその都度除去し、なくなるまで回帰直線近似を行う。
2)収束した回帰直線を用いて除去されたデータを補正する。
【００５４】
以上のことから、横断輝度分布法は主に島による解析法と回帰直線法で構成される。これらの方法は、従来の対話式の2値化法と連結法と比べて、次の点で優れており、最終的に医者の介在を必要としない骨軸の自動解析を可能にした。
1)島による解析は、2値化法に比べて、輝度値（島の高さ）の情報をより活用する。
2)回帰直線法は骨軸が直線近似で行うことを利用して、ノイズの検出と補正を行い、連結法より効率的である。
【００５５】
[C4：骨境界点から骨軸解析] 
次に骨境界点から骨軸を解析するプロセスC4について説明する。図12に骨境界点の抽出および骨軸解析のアルゴリズムを示すフローチャート図を示す。骨軸解析は大腿骨と脛骨に対して、骨幹の左右の骨境界点から、それぞれの骨軸を求めるものである。この解析で、次の２つの解析が重要となる。
(1)骨軸解析の正確な範囲の設定
(2)左右の骨境界点からの骨軸の求め方
【００５６】
上記(1)について、レントゲン写真において、大腿骨と脛骨の角度がある場合、膝間隙から上下10cm〜15cmの範囲の正確な決定法である。骨軸は左右の骨境界線で決まるので、骨軸解析範囲は骨境界線の解析範囲となる。よって、骨境界点から求められる中点よる中線によって、解析範囲を決定することが可能となる。大腿骨の場合に対して、図21を用いて解析手順を示す。
1)骨境界点を求める範囲は、骨境界線の解析範囲より十分広い範囲とした。実際はj軸座標(1〜Jmask1)とした。Jmask1は膝間隙からの距離が8cmのj軸座標である。
2)はじめに仮の解析範囲(膝間隙ラインから10cm〜15cmの範囲(J15〜J10))で左右の骨境界点とその中点から回帰直線法で、仮の左右の骨境界線と中線を求める。
3)中線と膝間隙ラインとの交点Cとする。中線上にBC=10cmとなるB点、AC=15cmとなるA点を求める。
4)A点とB点で中線と直交する直線を求め、2)で得られた仮の左右の骨境界線との交点(A1,A2,B1,B2)を計算する。これらから、正確な骨境界線の解析範囲が決まる。すなわち、右骨境界線に対して範囲(A2〜B2)、左骨境界線に対して範囲(A1〜B1)となる。よって、本発明の骨軸解析範囲は、中線上の点Ａから点Ｂの範囲（j軸座標J15m〜J10mの範囲）で定義される。
上記は、大腿骨の骨軸解析範囲について述べたが、同様に脛骨についても求めることができる。
【００５７】
(2)については、上記の正確な解析範囲の骨境界点から、回帰直線法で骨境界線を計算する。左右の骨境界線から骨軸を求める計算法は、次の数学の公式を使用する。２つの直線(y=b1x+a1とy=b2x+a2)から等距離にある直線は、次式で与えられる。
0005397873.tif 000010
上記の２つの直線を左右の骨境界線とすれば、式(7)より骨軸の直線が得られることになる。
解析された膝間隙と、大腿骨と脛骨の骨軸をグラフ（図22）と画像（図23）で示した。
これらの２つの骨軸の直線情報(y=b3x+a3とy=b4x+a4とする)から、大腿脛骨外側角(FTA)は容易に次式から計算される。
0005397873.tif 000011
【００５８】
次に上記で説明した骨軸自動抽出方法を実施するためにコンピュータ上で実行されるプログラムについて説明する。このプログラムがコンピュータにより実行されれば、上記で説明した骨軸自動抽出方法におけると基本的に同じ原理に従い、同様な作用効果が実現され得る。
【００５９】
図24は全置換型人工膝関節手術の術前計画をコンピュータでシミュレーションするシステムにおいて、本発明の骨軸自動抽出プログラムが組み込まれた実施例である（非特許文献３）。この実施例から、膝間隙ラインおよび大腿骨・脛骨の骨軸が自動計算され、画面上には効率的にテンプレートが配置され、インプラントを設置しやすくしている。さらにマニュアル操作により骨参照点およびテンプレート位置の調整を行い、医師の豊富な経験を反映させることも可能としている。
【００６０】
この実施例からも、膝間隙と大腿骨及び脛骨の骨軸の抽出をコンピュータ上で自動解析して行うことで、医師の術前計画における負担を軽減するとともに、これまで技能に頼ってきた外科医による術前計画および変形性膝関節症などの診断のレベルを一定以上に保つ効果があることがわかる。
【図面の簡単な説明】
【００６１】
【図１】従来法の目視・手作業法を示す図である。
【図２】本発明に係るアルゴリズムの全体を示すフローチャート図である。
【図３】本発明に係る下肢解析画像の切出しアルゴリズムを示すフローチャート図である。
【図４】本発明に係る画像の座標系を示す図である。
【図５】本発明に係る下肢透視画像（右図、j軸2140ピクセル i軸1760ピクセル）と画像切出し後の下肢解析画像（左図、j軸1935ピクセル i軸351ピクセル）を示す図である。
【図６】本発明に係るi軸方向の画像切出の区間（ILine1,ILine2）を示す図である。i軸方向幅1760ピクセル数が区間(ILine2-ILine1)=351ピクセル数に短縮される。
【図７】本発明に係るj軸方向の画像切出の区間（Jpst,Jped）を示す図である。j軸方向高さ2140ピクセル数が区間(Jped-Jpst)=1935ピクセル数に短縮される。
【図８】本発明に係る区間(JLine1,JLine2)における輝度変動分布Javbd[j]を示す図である。
【図９】本発明に係る膝間隙ラインの抽出アルゴリズムを示すフローチャート図である。
【図１０】本発明に係る膝間隙(Jrgk)を求めるj軸方向の勾配分布を示す図である。（DRPressY無は有のピーク値に合わせるようにレベル調整を行っている）
【図１１】本発明に係る膝間隙ライン(Jrgk)を抽出した図である。
【図１２】本発明に係る骨境界点および骨軸の解析アルゴリズムを示すフローチャート図である。
【図１３】長骨の外観（左図）と縦断面（右図）を示す図である。
【図１４】本発明に係る吸収画像（左図）とボケ画像（中図）と圧縮画像（右図）を示す図である。
【図１５】本発明に係る強調画像（左図）とピーク強調画像（中図）とエッジ画像（右図）を示す図である。
【図１６】本発明に係る原画像とボケ画像と圧縮画像と強調画像の横断輝度分布の比較を示す図である。
【図１７】本発明に係る原画像とボケ画像と圧縮画像と強調画像の比較を示す図である。
【図１８】本発明に係る吸収画像と圧縮画像と強調画像の輝度の頻度分布を示す図である。
【図１９】本発明に係る強調画像からエッジ画像を作成する時に用いるメキシカンハットフィルター(9×9)を示す図である。
【図２０】本発明に係る強調画像とエッジ画像の横断輝度分布を示す図である。
【図２１】本発明に係る骨軸解析範囲の求め方（大腿骨の場合）を説明する図である。
【図２２】本発明に係る大腿骨と脛骨の骨軸解析範囲と骨軸を示す図である。
【図２３】本発明に係る膝間隙と大腿骨、脛骨の骨軸の下肢画像を示す図である。
【図２４】本発明に係る全置換型人工膝関節手術の術前計画支援システムの実施例を示す図である。
【図１】
0005397873.tif 000016
【図２】
0005397873.tif 000017
【図３】
0005397873.tif 000018
【図４】
0005397873.tif 000019
【図５】
0005397873.tif 000020
【図６】
0005397873.tif 000021
【図７】
0005397873.tif 000022
【図８】
0005397873.tif 000023
【図９】
0005397873.tif 000024
【図１０】
0005397873.tif 000025
【図１１】
0005397873.tif 000026
【図１２】
0005397873.tif 000027
【図１３】
0005397873.tif 000028
【図１４】
0005397873.tif 000029
【図１５】
0005397873.tif 000030
【図１６】
0005397873.tif 000031
【図１７】
0005397873.tif 000032
【図１８】
0005397873.tif 000033
【図１９】
0005397873.tif 000034
【図２０】
0005397873.tif 000035
【図２１】
0005397873.tif 000036
【図２２】
0005397873.tif 000037
【図２３】
0005397873.tif 000038
【図２４】
0005397873.tif 000039