JP5397873B2 - 膝関節の大腿骨及び脛骨の骨軸自動抽出方法、ならびに骨軸自動抽出プログラム - Google Patents
膝関節の大腿骨及び脛骨の骨軸自動抽出方法、ならびに骨軸自動抽出プログラム Download PDFInfo
- Publication number
- JP5397873B2 JP5397873B2 JP2007183768A JP2007183768A JP5397873B2 JP 5397873 B2 JP5397873 B2 JP 5397873B2 JP 2007183768 A JP2007183768 A JP 2007183768A JP 2007183768 A JP2007183768 A JP 2007183768A JP 5397873 B2 JP5397873 B2 JP 5397873B2
- Authority
- JP
- Japan
- Prior art keywords
- bone
- axis
- image
- luminance
- line
- 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.)
- Expired - Fee Related
Links
- 210000000988 bone and bone Anatomy 0.000 title claims description 233
- 238000000605 extraction Methods 0.000 title claims description 59
- 210000002303 tibia Anatomy 0.000 title claims description 41
- 210000000689 upper leg Anatomy 0.000 title claims description 26
- 210000000629 knee joint Anatomy 0.000 title claims description 16
- 238000000034 method Methods 0.000 claims description 132
- 238000009826 distribution Methods 0.000 claims description 82
- 238000004458 analytical method Methods 0.000 claims description 59
- 210000003127 knee Anatomy 0.000 claims description 49
- 210000003414 extremity Anatomy 0.000 claims description 33
- 238000010521 absorption reaction Methods 0.000 claims description 31
- 210000003141 lower extremity Anatomy 0.000 claims description 20
- 235000005078 Chaenomeles speciosa Nutrition 0.000 claims description 19
- 240000000425 Chaenomeles speciosa Species 0.000 claims description 19
- 238000007906 compression Methods 0.000 claims description 19
- 230000006835 compression Effects 0.000 claims description 18
- 238000012545 processing Methods 0.000 claims description 12
- 210000001519 tissue Anatomy 0.000 claims description 11
- 238000009499 grossing Methods 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 244000286916 Ratibida columnifera Species 0.000 claims description 5
- 235000009413 Ratibida columnifera Nutrition 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 5
- 238000009795 derivation Methods 0.000 claims description 5
- 239000000284 extract Substances 0.000 claims description 5
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 19
- 238000004422 calculation algorithm Methods 0.000 description 10
- 238000003745 diagnosis Methods 0.000 description 7
- 238000002591 computed tomography Methods 0.000 description 6
- 230000002452 interceptive effect Effects 0.000 description 6
- 230000000007 visual effect Effects 0.000 description 6
- 208000003947 Knee Osteoarthritis Diseases 0.000 description 5
- 238000007796 conventional method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 210000002414 leg Anatomy 0.000 description 4
- 238000011883 total knee arthroplasty Methods 0.000 description 4
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 201000008482 osteoarthritis Diseases 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000001356 surgical procedure Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 210000003275 diaphysis Anatomy 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 239000007943 implant Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Description
目視方法は医者がレントゲン写真を目視し、次の手順で骨軸の抽出を手作業で行う方法である(図1)。
A1:膝関節の間隙線を目視して膝間隙ラインの検出し、そのラインから上下に10cm〜15cmの骨軸抽出範囲の設定
A2:上方の骨軸抽出範囲で大腿骨の骨領域の検出し、骨部の左右の骨境界線を抽出して、それらの境界線の中線を大腿骨の骨軸とする。
A3:下方の骨軸抽出範囲で、A2と同様な骨境界線検出を行い、その中線を脛骨の骨軸とする。これにより、大腿骨と脛骨の骨軸が抽出され、さらに診療の重要項目の大腿脛骨外側角(FTA)が容易に算出できることになる。
対話式方法は、医者が介在する対話方式によるCTのデジタルデータによるコンピュータ処理方法である(非特許文献1)。すなわち、コンピュータが困難な処理を医者が行い、それをもとに骨軸の抽出を行う方法である。医者が対話方法で関与する情報は、膝間隙ラインと抽出したい骨の輝度値の2つである。すなわち、上記の目視・手作業方法と比較して、与えられた膝間隙ラインから目視方法のA1ができる。A2とA3については与えられた骨の輝度値の情報をもとに、次の画像処理法で処理が可能となる。
B1:閾値法による骨領域の抽出
骨部のX線透視画像の輝度値は患者によっても撮影条件によっても異なる。よって、対話式に与えられた輝度値を必要とする。この値をもとに、同じ輝度値をもつ骨領域の抽出を行うことが出来る。
B2:エッジ検出法および輪郭連結法
B1で求められた骨領域を2値化してエッジ法で骨境界の検出を行い、輪郭連結法によって、連続した骨境界を抽出して、これらの境界線の中線から、大腿骨と脛骨の骨軸が得られることになる。
(a)下肢透視画像から生体組織の画像領域の切り出しを行うプロセスと、
(b)プロセス(a)で切り出された透視画像から輝度反転変換した画像(以下、吸収画像と称する)を導出し、該吸収画像の輝度分布に基づき膝間隙ラインを抽出するプロセスと、
(c)プロセス(b)で求めた吸収画像から骨組織を強調した画像(以下、強調画像と称する)に変換し、該強調画像の輝度分布と前記膝間隙ラインとに基づき骨境界点を求めるプロセスと、
(d)プロセス(c)で求めた骨境界点と前記膝間隙ラインとに基づき骨軸解析範囲を設定し、該骨軸解析範囲の骨境界点から骨軸を導出するプロセスと、
を有し、プロセス(a)から(d)を順次実行することで膝関節の大腿骨及び脛骨の骨軸が自動抽出され、大腿脛骨外側角(FTA: Femoro Tibial Angle)が自動算出されることを特徴とする。
(1)各i軸座標でj軸方向に輝度分布の平均輝度値IavbM(i)を求める。
(2)IavbM(i)分布においてi軸の両端から輝度値がある輝度値以上に変化したi軸座標をそれぞれIRngA1とIRngA2とし、これらの点から幅(IRngA2−IRngA1)/δ(δ:予め設定された定数)の長さだけ内側にシフトした内部点として、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軸座標の切り出し点とする。
(1)各j軸座標でi軸方向に輝度分布の平均輝度値Javb[j]を求める。
(2)Javb[j]分布においてj軸の両端から輝度値が予め決められた設定値以上に変化した座標をそれぞれJLine1とJLine2とし、これらの座標をj軸方向の切り出し点とする。
ことを特徴とする。
(1)プロセス(a)で切り出された透視画像IBrt0(j,i)を輝度反転変換した吸収画像の輝度分布IBrt1(j,i)を式(1)により求める。
ここでIBmaxとIBminはそれぞれ透視画像の最大輝度値と最小輝度値とする。
(2)吸収画像の輝度分布IBrt1(j,i)を用いて式(2)により輝度差分分布DBrt(j)を求める。
(3)DBrt(j)のピークのj軸座標を膝間隙ラインとして抽出する。
ことを特徴とする。
(1)j軸方向のみにβピクセル(β:予め設定された定数)の移動平均による平滑化を行い、ボケ画像Boke(j,i)を作成する。
(2)ボケ画像Boke(j,i) と吸収画像IBrt1(j,i)とを利用して式(3)による輝度変動画像IBrt2(j,i)を作成する。
(3)プロセス(a)で切り出された透視画像の輝度範囲が(IBmin,IBmax)になるようにIBrt2(j,i)を線形補間し、γピクセル(γ:予め設定された定数)のスケールの移動平均による平滑化を行い、j軸方向の圧縮画像を得る。
ことを特徴とする。ここでβは比較的大きいスケール(例えばβ=50ピクセル)であり、γは比較的小さいスケール(例えばγ=3ピクセル)とする。
(1) プロセス(b)で求めた吸収画像からダイナミックレンジ圧縮画像処理を用いて骨組織の輝度レベルを強調する強調画像を導出する。
(2) 該強調画像の幅方向をi軸、高さ方向をj軸として、i軸方向の輝度分布に対して予め設定された輝度レベル以上の輝度を持つ分布の塊(以下、島と称する)を検出する。
(3) 前記(2)で検出した島が2個の場合、該2島を選択する。前記(2)で検出した島が3個以上の場合、前記膝間隙ラインに基づき大腿骨か脛骨かを判断し、
(3-1) 大腿骨の場合、島の高さの高い方から2島を選択する。
(3-2) 左肢の脛骨の場合、左側にある2島を優先的に選択する。
(3-3) 右肢の脛骨の場合、右側にある2島を優先的に選択する。
(4) 前記(3)で選択した2島からエッジ画像を利用して島の輝度値がピークとなる位置を求め、該位置を骨境界点として抽出する。
ことを特徴とする。
(1)j軸およびi軸の両方向にβピクセル(β:予め設定された定数)の移動平均による平滑化を行い、ボケ画像Boke(j,i)を作成する。
(2)ボケ画像Boke(j,i)と吸収画像IBrt1(j,i)とを利用して式(3)により圧縮画像IBrt2(j,i)を作成する。
(3)圧縮画像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)
ことを特徴とする。
(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)を骨境界線の解析範囲とする。
ことを特徴とする。
さらに骨境界の検出アルゴリズムに、骨組織に反映した輝度分布に基づく骨境界点抽出の解析法を用いることで、単なる輝度値の閾値法に基づく従来法に比べて、精度の信頼性を高める効果がある。
(1)レントゲン写真の外枠
(2)3枚の写真の継ぎ目線と留め金
(3)スケール表示
(4)マーク
C1:下肢透視画像から生体組織の画像領域の切り出しを行うプロセスと、
C2:プロセスC1で切り出された透視画像から輝度反転変換した画像(以下、吸収画像と称する)を導出し、該吸収画像の輝度分布に基づき膝間隙ラインを抽出するプロセスと、
C3:プロセスC2で求めた吸収画像から骨組織を強調した画像(以下、強調画像と称する)に変換し、該強調画像の輝度分布と前記膝間隙ラインとに基づき骨境界点を求めるプロセスと、
C4:プロセスC3で求めた骨境界点と前記膝間隙ラインとに基づき骨軸解析範囲を設定し、該骨軸解析範囲の骨境界点から骨軸を導出するプロセスと、
を有し、コンピュータ上でプロセスC1からC4を順次実行することで膝関節の大腿骨及び脛骨の骨軸が自動抽出され、大腿脛骨外側角(FTA)が自動算出される。以下に、C1からC4の各プロセスについて説明する。
プロセスC1は、生体組織画像以外のノイズ情報を除去するために、膝関節周辺の下肢透視画像の幅方向をi軸、高さ方向をj軸(画像の座標系を図4に示す)として、i軸方向及びj軸方向の輝度分布特性と左肢か右肢かの肢情報とを用いて、生体組織の画像領域だけを含む最小領域の切り出しを行う。図3に切り出しアルゴリズムを示すフローチャート図を示す。本切り出し法は、処理効率の向上のための解析画像サイズの最小化ばかりでなく、後述する骨境界抽出法で支障をもたらす人工的なノイズを除去することにある。人工ノイズは、図5(右図)からわかるように、例えば次のものがある。
(1)レントゲン写真の外枠
(2)3枚の写真の継ぎ目線と留め金
(3)スケール表示
(4)マーク
次にi軸方向の画像切り出し方法(C11)について説明する。図5の下肢透視画像からわかるように、骨のようなX線の吸収率が高いところは、低い輝度値になるために、黒い画像となる。生体組織を1つの塊として分離するために、画像切り出しは初めにi軸方向から行うのが合理的と言える。よって、使用する輝度の平均分布は、各i軸座標でj軸方向の輝度分布の平均値を求め、かつi軸方向に平滑化した分布(図6)とした。また、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)における2つの区間のうち、スケールがある区間を特定することができる。その区間においては、選択する峰は最大輝度値の峰より、IMinIti点により近い峰点とした。ここで峰点の定義はその点を挟む両側の2点、すなわち周りの4点より大きな輝度値を有する点である。
次にj軸方向の画像切り出し方法(C12)について説明する。j軸方向の切り出しはi軸方向の画像切り出し方法(C11)で処理された画像をもとに、次の2つの目的のために行う。
(1)両端の非生体画像領域の除去
(2)レントゲン写真の継ぎ目線(これをラインノイズと呼ぶ)の輝度補正
1)各j軸座標において、i軸方向に輝度分布の平均輝度値Javb[j](図7)を求める。
2) Javb[j]分布においてj軸の両端から輝度値が予め決められた設定値以上に変化した座標をそれぞれJLine1とJLine2とし、これらの座標をj軸方向の切り出し点とする。よって画像のj軸方向の大きさはHeightJ=JLine2−JLine1となる。
次に膝間隙ラインを抽出するプロセスC2について説明する。図9に膝間隙ラインの抽出アルゴリズムを示すフローチャート図を示す。プロセスC2は、C1で切り出された透視画像をもとに、膝間隙ラインの抽出を行う。
(1) 透視画像(輝度値IBrt0(j,i))を式(1)による輝度反転変換を行う。変換された画像(輝度値IBrt1(j,i))を吸収画像と呼ぶ。この変換によって、X線の吸収率が高い骨部は、高い輝度を有することになる。
j=1,…,M i=1,…,N
ここでIBmaxとIBminはそれぞれ透視画像の最大輝度値と最小輝度値である。ここで求める膝間隙ラインは間隙の平均的なj軸座標Jrgkとする。Jrgkの抽出法は次の通りである。
(2) 式(2)により吸収画像の輝度分布IBrt1(j,i)からj軸方向に1階差分の絶対値を求め、かつi軸方向に平均化を行い、輝度差分分布DBrt(j)を求める。
(3) DBrt(j)のピークのj軸座標が求める膝間隙ラインとなる。
(1)
j軸方向にβピクセル(例えばβ=50ピクセル、約2.5cm程度)の移動平均による平滑化を行い、ボケ画像Boke(j,i)を作成する。
(2)
ボケ画像Boke(j,i)を利用して式(3)による輝度変動画像IBrt2(j,i)を作る。
(3)
透視画像の輝度範囲が(IBmin,IBmax)になるようにIBrt2(j,i)を線形補間し、γピクセル(例えばγ=3ピクセル)のスケールの移動平均による平滑化を行い、j軸方向の圧縮画像を得る。
図10に|勾配|分布、DBrt(j)に対してDRPressYフィルター有無での効果を示した。この図から、次のようなDRPressYフィルターの効果が示された。このフィルターは平均値を含む比較的大きなスケールの変動を除去するために、得られた勾配分布は大きさと形において、より安定した分布となり、ピーク位置の抽出が容易となることである。
次に骨境界点を抽出するプロセスC3(骨境界抽出法)について説明する。図12に骨境界点の抽出および骨軸解析のアルゴリズムを示すフローチャート図を示す。プロセスC3は、
(1) プロセス(b)で求めた吸収画像(図14の左図)からダイナミックレンジ圧縮画像処理を用いて骨組織の輝度レベルを強調する強調画像を導出する(C31)。
(2) 該強調画像の幅方向をi軸、高さ方向をj軸として、i軸方向の輝度分布に対して予め設定された輝度レベル以上の輝度を持つ分布の塊(以下、島と称する)を検出する。
(3) 検出した島が2個の場合、該2島を選択する。検出した島が3個以上の場合、前記膝間隙ラインに基づき大腿骨か脛骨かを判断し、
(3-1) 大腿骨の場合、島の高さの高い方から2島を選択する。
(3-2) 左肢の脛骨の場合、左側にある2島を優先的に選択する。
(3-3) 右肢の脛骨の場合、右側にある2島を優先的に選択する。
(4) 前記(3)で選択した2島からエッジ画像を利用して島の輝度値がピークとなる位置を求め、該位置を骨境界点として抽出する(C32:横断輝度分布法)。
次にダイナミックレンジ圧縮強調画像処理(C31)について説明する。求める骨軸は膝間隙から上下に骨軸解析範囲にある、大腿骨と脛骨の骨幹の骨軸である。解析したい骨幹は図13から示されるように、中心部に髄腔という空洞を持ち、幹壁として硬い緻密質(骨)を有する。開発したダイナミックレンジ圧縮強調画像処理は、このような骨組織を反映させる画像処理である。具体的な解析手順は次の通りである。以下で示す画像は吸収画像と輝度範囲が同一になるように規格化される。
(1)ボケ画像Boke(j,i)(図14の中図)の作成
j軸およびi軸の両方の方向にβピクセル(β:予め設定された定数)の移動平均による平滑化を行い、ボケ画像Boke(j,i)を作成する。
(2)j軸およびi軸方向の圧縮画像IBrt2(j,i)(図14の右図)の作成
ボケ画像Boke(j,i) と吸収画像IBrt1(j,i)とを利用して式(3)により圧縮画像IBrt2(j,i)を作成する。
得られた圧縮画像は図13との比較から、骨組織を反映した画像と言える。すなわち、緻密質が最も高い輝度値で、髄腔が低い輝度値の画像である。
(3)強調画像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)は、高い輝度骨部に対して1以上で増幅され、逆に低い輝度の生体組織では1以下で減衰される。従って、強調画像は式(5)で得られ、圧縮画像に比べ、解析する大腿骨と脛骨のみがより高い輝度の画像となる(図15の左図)。このことは代表的な横断輝度分布(図16と図17)からも明白に示される。
I.使用したDR圧縮処理は、輝度レベルIBrtLVL(j)付近の輝度の頻度を集中させる。
II.強調処理は、上記の輝度レベルで規格化されたボケ画像で作られた強調係数を掛けることで、さらに骨部の輝度の集中化を促進する。
上記IIについては、次のピーク強調画像(図15の中図)で明らかに示される。すなわち、ピーク強調画像は、図18のピーク頻度のみの輝度値で作成した2値画像で、骨部が明確に存在することで示される。ピーク強調画像の作成法は、図18の強調画像において、ピークの輝度値の範囲を440〜450とすると、440以下の輝度値は440とし、450以上の輝度値は450とする輝度変換するものである。
次に横断輝度分布法(C32)について説明する。この方法は骨軸解析範囲において、強調画像の横断輝度分布をもとに、エッジ画像を利用して骨境界点を求める。エッジ画像は、図19で示した9×9行列で表されるメキシカンハットフィルター(非特許文献(2))を強調画像に適用して作成し、かつ式(4)のfact(j,i)の処理(式(5))をした(図15の右図)。そして、骨境界点の回帰直線補正により、骨境界点に混在するノイズの検出と補正を行う。
1)強調画像において、正の輝度値を有する分布の塊(今後島と呼ぶ)を検出し、島を高さと位置(島のピーク輝度点)でパラメータ化し、島の高さの高い方から2島を選択する。これらの島が図13の緻密質に対応する。図20に従って、島情報(i軸座標であるAk、Bk、Ck:k=1,2)を決める。中心にある髄腔に対して、添字1が左側で、添字2が右側を示す。
2)エッジ画像は強調画像で示した島、すなわち緻密質の境界点を与える。図20から、ピーク点Dkが骨外境界点に、ゼロ点Ekが骨境界点に対応することになる。
3) 脛骨は図20の下図で示すように大腿骨と異なって、近傍に輝度の高い腓骨を伴う。このために、右肢か左肢かの情報をもとに島の位置情報を利用して、脛骨の島のみを選択するようにした。具体的には、強調画像において、3つ以上の島が存在するとき、右肢の場合腓骨は必ず脛骨の左側にあるために、右側にある島を優先的に選択する。左肢の場合逆に腓骨は右側になるので、左側の島が優先的に選択すればよい。
1)骨境界点に回帰直線近似を行い、標準偏差以上外れたデータはその都度除去し、なくなるまで回帰直線近似を行う。
2)収束した回帰直線を用いて除去されたデータを補正する。
1)島による解析は、2値化法に比べて、輝度値(島の高さ)の情報をより活用する。
2)回帰直線法は骨軸が直線近似で行うことを利用して、ノイズの検出と補正を行い、連結法より効率的である。
次に骨境界点から骨軸を解析するプロセスC4について説明する。図12に骨境界点の抽出および骨軸解析のアルゴリズムを示すフローチャート図を示す。骨軸解析は大腿骨と脛骨に対して、骨幹の左右の骨境界点から、それぞれの骨軸を求めるものである。この解析で、次の2つの解析が重要となる。
(1)骨軸解析の正確な範囲の設定
(2)左右の骨境界点からの骨軸の求め方
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)となる。よって、本発明の骨軸解析範囲は、中線上の点Aから点Bの範囲(j軸座標J15m〜J10mの範囲)で定義される。
上記は、大腿骨の骨軸解析範囲について述べたが、同様に脛骨についても求めることができる。
上記の2つの直線を左右の骨境界線とすれば、式(7)より骨軸の直線が得られることになる。
解析された膝間隙と、大腿骨と脛骨の骨軸をグラフ(図22)と画像(図23)で示した。
これらの2つの骨軸の直線情報(y=b3x+a3とy=b4x+a4とする)から、大腿脛骨外側角(FTA)は容易に次式から計算される。
Claims (14)
- 下肢透視画像をコンピュータ上で自動解析して膝関節の大腿骨及び脛骨の骨軸を自動抽出する方法であって、
(a)下肢透視画像から生体組織の画像領域の切り出しを行うプロセスと、
(b)プロセス(a)で切り出された透視画像から輝度反転変換した画像(以下、吸収画像と称する)を導出し、該吸収画像の輝度分布に基づき膝間隙ラインを抽出するプロセスと、
(c)プロセス(b)で求めた吸収画像から骨組織を強調した画像(以下、強調画像と称する)に変換し、該強調画像の輝度分布と前記膝間隙ラインとに基づき骨境界点を求めるプロセスと、
(d)プロセス(c)で求めた骨境界点と前記膝間隙ラインとに基づき骨軸解析範囲を設定し、該骨軸解析範囲の骨境界点から骨軸を導出するプロセスと、
を有し、プロセス(a)から(d)を順次実行することで膝関節の大腿骨及び脛骨の骨軸が自動抽出され、大腿脛骨外側角が自動算出されることを特徴とする骨軸自動抽出方法。 - プロセス(a)において、下肢透視画像の幅方向をi軸、高さ方向をj軸として、i軸方向及びj軸方向の輝度分布特性と左肢か右肢かの肢情報とを用いて、前記下肢透視画像から生体組織画像以外のノイズ情報を除去した最小領域の切り出しを行うことを特徴とする請求項1に記載の骨軸自動抽出方法。
- 請求項2に記載の切り出しにおいてi軸方向の画像切り出し手順は、
(1)各i軸座標でj軸方向に輝度分布の平均輝度値IavbM(i)を求める。
(2)IavbM(i)分布においてi軸の両端から輝度値がある輝度値以上に変化したi軸座標をそれぞれIRngA1とIRngA2とし、これらの点から幅(IRngA2−IRngA1)/δ(δ:予め設定された定数)の長さだけ内側にシフトした内部点として、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軸座標の切り出し点とする。
ことを特徴とする請求項2に記載の骨軸自動抽出方法。 - 請求項2に記載の切り出しにおいてj軸方向の画像切り出し手順は、
(1)各j軸座標でi軸方向に輝度分布の平均輝度値Javb[j]を求める。
(2)Javb[j]分布においてj軸の両端から輝度値が予め決められた設定値以上に変化した座標をそれぞれJLine1とJLine2とし、これらの座標をj軸方向の切り出し点とする。
ことを特徴とする請求項2に記載の骨軸自動抽出方法。 - 請求項4に記載のj軸方向の画像切り出し手順において、切り出し区間(JLine1,JLine2)で、Javb[j]の分布をj軸方向に平滑化した分布JavbM[j]を求め、変動分布JavbD[j]=Javb[j]−JavbM[j]を用いてレントゲン写真の継ぎ目線(以下、ラインノイズと称する)の区間を特定し、該ラインノイズ区間の輝度分布を区間外の輝度分布から線形補間して補正することによりラインノイズを消去することを特徴とする請求項4に記載の骨軸自動抽出方法。
- 請求項6に記載の膝間隙ラインの抽出において、前処理としてj軸方向にダイナミックレンジ圧縮(以下、DR圧縮と称する)して輝度変動量を増幅させた輝度分布を導出すること、およびDR圧縮のフィルター処理は、
(1)j軸方向のみにβピクセル(β:予め設定された定数)の移動平均による平滑化を行い、ボケ画像Boke(j,i)を作成する。
(2)ボケ画像Boke(j,i) と吸収画像IBrt1(j,i)とを利用して式(3)による輝度変動画像IBrt2(j,i)を作成する。
(3)プロセス(a)で切り出された透視画像の輝度範囲が(IBmin,IBmax)になるようにIBrt2(j,i)を線形補間し、γピクセル(γ:予め設定された定数)のスケールの移動平均による平滑化を行い、j軸方向の圧縮画像を得る。
ことを特徴とする請求項6に記載の骨軸自動抽出方法。 - プロセス(c)において、
(1)
プロセス(b)で求めた吸収画像からダイナミックレンジ圧縮画像処理を用いて骨組織の輝度レベルを強調する強調画像を導出する。
(2)
該強調画像の幅方向をi軸、高さ方向をj軸として、i軸方向の輝度分布に対して予め設定された輝度レベル以上の輝度を持つ分布の塊(以下、島と称する)を検出する。
(3)
前記(2)で検出した島が2個の場合、該2島を選択する。前記(2)で検出した島が3個以上の場合、前記膝間隙ラインに基づき大腿骨か脛骨かを判断し、
(3-1)
大腿骨の場合、島の高さの高い方から2島を選択する。
(3-2)
左肢の脛骨の場合、左側にある2島を優先的に選択する。
(3-3)
右肢の脛骨の場合、右側にある2島を優先的に選択する。
(4)
前記(3)で選択した2島からエッジ画像を利用して島の輝度値がピークとなる位置を求め、該位置を骨境界点として抽出する。
ことを特徴とする請求項1に記載の骨軸自動抽出方法。 - 請求項8に記載した強調画像の導出は、
(1)j軸およびi軸の両方向にβピクセル(β:予め設定された定数)の移動平均による平滑化を行い、ボケ画像Boke(j,i)を作成する。
(2)ボケ画像Boke(j,i)と吸収画像IBrt1(j,i)とを利用して式(3)により圧縮画像IBrt2(j,i)を作成する。
(3)圧縮画像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)
ことを特徴とする請求項8に記載の骨軸自動抽出方法。 - 請求項8に記載の骨境界点の抽出において、前記強調画像にメキシカンハットフィルターを適用してエッジ画像を導出し、該エッジ画像を利用して前記強調画像の輝度分布を基に骨境界点を求めることを特徴とする請求項8に記載の骨軸自動抽出方法。
- 請求項8に記載の骨境界点の抽出において、骨境界点から高さ方向のj軸方向に回帰直線近似法を用いて、骨境界点に混在するノイズの検出と補正を行うことを特徴とする請求項8に記載の骨軸自動抽出方法。
- プロセス(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)を骨境界線の解析範囲とする。
ことを特徴とする請求項1に記載の骨軸自動抽出方法。 - プロセス(d)における骨軸の導出手順は、前記骨軸解析範囲の骨境界点から、回帰直線法を用いて左右の骨境界線を導出し、この左右の骨境界線から等距離にある直線を骨軸の直線として導出することを特徴とする請求項1に記載の骨軸自動抽出方法。
- 請求項1乃至請求項13のいずれかに記載の骨軸自動抽出方法をコンピュータに実行させることを特徴とする骨軸自動抽出プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007183768A JP5397873B2 (ja) | 2007-07-12 | 2007-07-12 | 膝関節の大腿骨及び脛骨の骨軸自動抽出方法、ならびに骨軸自動抽出プログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007183768A JP5397873B2 (ja) | 2007-07-12 | 2007-07-12 | 膝関節の大腿骨及び脛骨の骨軸自動抽出方法、ならびに骨軸自動抽出プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2009018065A JP2009018065A (ja) | 2009-01-29 |
JP5397873B2 true JP5397873B2 (ja) | 2014-01-22 |
Family
ID=40358215
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007183768A Expired - Fee Related JP5397873B2 (ja) | 2007-07-12 | 2007-07-12 | 膝関節の大腿骨及び脛骨の骨軸自動抽出方法、ならびに骨軸自動抽出プログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5397873B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11564652B2 (en) | 2019-10-04 | 2023-01-31 | Konica Minolta, Inc. | Storage medium, dynamic analysis apparatus, and dynamic analysis system |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2926064T3 (es) * | 2015-09-04 | 2022-10-21 | Straxcorp Pty Ltd | Método y aparato para identificar un espacio entre objetos en una imagen |
CN110956622B (zh) * | 2019-11-28 | 2023-12-22 | 天津市天津医院 | 一种从人体x射线图像中自动提取膝关节部分图像的方法 |
CN113177915B (zh) * | 2021-04-20 | 2024-04-09 | 中国科学院高能物理研究所 | 一种胫骨平台后倾角测量方法及装置、存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3639825B2 (ja) * | 2002-04-03 | 2005-04-20 | キヤノン株式会社 | 動画像表示方法、プログラム、コンピュータ可読記憶媒体、及び動画像表示装置 |
GB0504172D0 (en) * | 2005-03-01 | 2005-04-06 | King S College London | Surgical planning |
-
2007
- 2007-07-12 JP JP2007183768A patent/JP5397873B2/ja not_active Expired - Fee Related
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11564652B2 (en) | 2019-10-04 | 2023-01-31 | Konica Minolta, Inc. | Storage medium, dynamic analysis apparatus, and dynamic analysis system |
US11963811B2 (en) | 2019-10-04 | 2024-04-23 | Konica Minolta, Inc. | Storage medium, dynamic analysis apparatus, and dynamic analysis system |
US11998380B2 (en) | 2019-10-04 | 2024-06-04 | Konica Minolta, Inc. | Storage medium, dynamic analysis apparatus, and dynamic analysis system |
Also Published As
Publication number | Publication date |
---|---|
JP2009018065A (ja) | 2009-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5634880B2 (ja) | カスタマイズされた外科器具又はプロテーゼを生産する方法並びに身体部分の患者別モデルを生成する方法及びそのためのプログラムコードを有する媒体 | |
CN108765417B (zh) | 一种基于深度学习和数字重建放射影像的股骨x线片生成***及方法 | |
US8200355B2 (en) | Method for manufacturing artificial implants | |
US8233692B2 (en) | Method of suppressing obscuring features in an image | |
Almeida et al. | Fully automatic segmentation of femurs with medullary canal definition in high and in low resolution CT scans | |
JP2020513869A (ja) | 頭蓋骨を復元する方法 | |
Naranjo et al. | Metal artifact reduction in dental CT images using polar mathematical morphology | |
Chougule et al. | Clinical case study: spine modeling for minimum invasive spine surgeries (MISS) using rapid prototyping | |
CN112509022A (zh) | 一种术前三维影像与术中透视图像的无标定物配准方法 | |
JP5397873B2 (ja) | 膝関節の大腿骨及び脛骨の骨軸自動抽出方法、ならびに骨軸自動抽出プログラム | |
CN103366397A (zh) | 基于c形臂2d投影图像的脊柱3d模型构建方法 | |
Liu et al. | Virtual plate pre-bending for the long bone fracture based on axis pre-alignment | |
CN114897703A (zh) | 一种骨科全长x线片自动拼接***、方法、计算机设备及存储介质 | |
Gelaude et al. | Semi-automated segmentation and visualisation of outer bone cortex from medical images | |
CN111292306A (zh) | 膝关节ct与mr图像融合方法及装置 | |
CN105894508A (zh) | 一种医学图像自动定位质量的评估方法 | |
Rashid et al. | Geometrical model creation methods for human humerus bone and modified cloverleaf plate | |
CN117152126B (zh) | 膝关节状态的确定装置、存储介质及电子设备 | |
CN116277978B (zh) | 多模态骨关节数字化3d打印方法 | |
Laine et al. | Image processing for femoral endosteal anatomy detection: description and testing of a computed tomography based program | |
Stojković et al. | User defined geometric feature for the creation of the femoral neck enveloping surface | |
JP6785237B2 (ja) | チャネル検出を使用する画像データ内の物体の分割 | |
CN111145353B (zh) | 通过图像分割和网格特征点提取算法生成3d点云的方法 | |
CN112842372A (zh) | 一种基于影像扫描的股骨头病变体积的测量方法 | |
Wei et al. | Segmentation of lung lobes in clinical CT images |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20100712 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20120229 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120328 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120516 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20120518 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20121121 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130116 |
|
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: 20130918 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20131016 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5397873 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |