JP4458309B2 - 3D image display device - Google Patents

3D image display device Download PDF

Info

Publication number
JP4458309B2
JP4458309B2 JP26427899A JP26427899A JP4458309B2 JP 4458309 B2 JP4458309 B2 JP 4458309B2 JP 26427899 A JP26427899 A JP 26427899A JP 26427899 A JP26427899 A JP 26427899A JP 4458309 B2 JP4458309 B2 JP 4458309B2
Authority
JP
Japan
Prior art keywords
region
curved surface
image
curved
voxel
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 - Lifetime
Application number
JP26427899A
Other languages
Japanese (ja)
Other versions
JP2001061834A (en
Inventor
元章 齊藤
Original Assignee
テラリコン・インコーポレイテッド
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by テラリコン・インコーポレイテッド filed Critical テラリコン・インコーポレイテッド
Priority to JP26427899A priority Critical patent/JP4458309B2/en
Publication of JP2001061834A publication Critical patent/JP2001061834A/en
Application granted granted Critical
Publication of JP4458309B2 publication Critical patent/JP4458309B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Description

【0001】
【発明の属する技術】
本発明は,被検体の物理的性質の空間的分布,特に上顎または下顎,上下顎などの湾曲した領域の性質を画像として表示する三次元画像表示装置に関する。
【0002】
【従来の技術】
歯科領域では,上顎または下顎,上下顎を連続した画像として1枚のフィルムに描出する撮影法であるパノラマX線撮影が行われている。顎関節および上下顎骨,歯列の状態が一枚のフィルム上で観察できるので,歯科領域では広く使用されている。
【0003】
X線CT装置によって人体の物理的性質をあらわす正確な断面画像データが得られるようになって以来,異なる断面位置で撮影した複数の断面画像データを使用して三次元立体画像を再構成することが行われてきた。特に最近,ヘリカルスキャンX線CT装置やコーンビームX線CT装置が実用化されて,より精密な三次元立体画像を再構成することが可能になったので,歯科領域でもX線CT装置の利用が進んでいる。
【0004】
X線CT装置で収集したデータを使用して歯列の状態を観察する方法として,歯科用曲面再構成と呼ばれる画像処理が行われている。これはX線CT装置で撮影したデータから作成した上顎または下顎,上下顎を含む三次元ボクセルデータを使用する。歯列に沿った複数の曲面とその曲面に直交する多数の平面を設定し,歯列に沿った複数の曲面のパノラマ状断面像とその曲面に直交する多数の平面の断面像を作成,表示する方法である。この歯科用曲面再構成と呼ばれる画像処理法を使用することによって,顎関節および上下顎骨,歯列の断面を展開して,パノラマX線撮影に類似した様式で観察できるようになった。
【0005】
図7はこの従来の歯科用曲面再構成を説明するブロック図である。61はX線CT装置で撮影した上顎または下顎,上下顎の画像データを使用して構築したCT三次元ボクセルデータ66と,断面像の断面位置を設定する設定面67を示す。64は歯列に沿った複数の曲面とその曲面に直交する多数の平面を設定する処理を示す。62は設定面における上顎または下顎,上下顎の領域に歯列に沿って設定した複数の曲面68とこれに直交する多数の平面69を示している。65は歯列に平行な曲面の断面像とこれに直交する平面の断面像を作成する処理を示す。63は歯列に平行な複数の曲面の断面像を平面に展開した画像70とこれに直交する多数の平面の断面像71を模式的に示す。歯科用曲面再構成では,歯列に平行な曲面の断面像は1mm程度の間隔で作成し,またこれに直交する平面の断面像も1mm程度の間隔で作成するので,実際にはこれよりも多数の画像が作成される。
【0006】
パノラマX線撮影で得られる画像はX線の透過像が一枚であるが,歯科用曲面再構成で得られる画像は歯列に沿って設定した複数枚の曲面の断面像とその曲面に直交する多数枚の断面像である。従って,パノラマX線撮影では一枚の画像を観察することによって歯列全体の情報を把握することが可能であるが,歯科用曲面再構成では歯列全体の情報を把握するためには歯列に沿って設定した複数枚の曲面の断面像とその曲面に直交する多数枚の断面像を同時に観察することが必要になる。歯科用曲面再構成で作成された画像を全て一枚のフィルムにプリントする場合には歯列全体の情報の把握は可能であるが,ワークステーションの表示画面上で対話式に処理を行いながら作成した画像を表示する場合には歯列全体を観察するためには多数の断面像を順次観察することになるので立体構造を直感的に把握することが難しい欠点がある。
【0007】
【発明が解決しようとする課題】
従来の歯科用曲面再構成法においては,情報に漏れがないようにするためには歯列に沿って設定する曲面とこれに直交する平面の設定を慎重に行う必要があり,また間隔を密に設定することが必要になるので断面像の数が膨大になる。従来の歯科用曲面再構成法においては,歯列全体の情報を把握するためには歯列に沿って設定した複数枚の曲面の断面像とその曲面に直交する多数枚の断面像を観察する必要があるために,対話式の操作が難しく,また立体構造を直感的に把握することが難しい。上顎または下顎,上下顎のような湾曲した厚板状解剖学的領域の情報を全て一枚に含んだ平面展開画像を表示すると共に,対話式の操作と立体構造の直感的な把握を可能にする画像表示装置が必要である。
【0008】
【課題を解決するための手段】
本発明は,従来の歯科用曲面再構成法における欠点を改善するために考案したものである。請求項1の発明は,図2において,上顎または下顎,上下顎を含む三次元ボクセルデータ41に歯列に沿った基準曲面46と,歯列の関心オブジェクトを含む湾曲した厚板状領域43とを設定する手段28と,この基準曲面46の垂線方向48における湾曲した厚板状領域43での光の減衰をボクセルの画素値から計算してその結果を基準曲面46に投影することによって,基準曲面46上に関心オブジェクトの投影画像データ49を作成する手段29と,基準曲面46を平面に展開する関係式を作成してこの関係式に従って基準曲面上の画像データ49を平面に展開する手段30とを備えることによって,パノラマ平面画像25を作成して表示するようにした。従来の歯科用曲面再構成法においては断面像を使用しているために体積データの画像化には多数の断面像を必要とするが,本発明では歯列に沿って設定した基準曲面に垂直な方向における関心オブジェクトでの光の減衰を計算して画像を作成するので,関心オブジェクトの全ての情報を一枚の画像によって,パノラマX線撮影による画像と同様に,画像化することが可能である。
【0009】
請求項3の発明は,図2において,上顎または下顎,上下顎を含む三次元ボクセルデータ41に歯列に沿った基準曲面46と,歯列の関心オブジェクトを含む湾曲した厚板状領域43とを設定する手段28と,新しい三次元ボクセル空間24と基準平面55を用意し,基準曲面46を基準平面55に変換する関係式を作成し,この関係式に従って基準曲面に基づいた湾曲した厚板状領域43のボクセルデータを,基準平面に基づいた新しい三次元ボクセル空間のボクセルに写像することによって変換された厚板状の関心オブジェクト52を作成する手段31と,基準平面55に垂直な方向における非湾曲厚板領域52の関心オブジェクトでの光の減衰をボクセルの画素値から計算してその結果を基準平面55上に投影する手段32とを備えることによって,パノラマ状の平面画像25を作成して表示するようにした。
【0010】
本発明では,CTデータのボクセル空間に歯列に沿った基準曲面46と関心オブジェクトを含む湾曲した厚板状領域43を設定する。請求項1の発明は,基準曲面46の垂線方向48における湾曲した厚板状領域43での光の減衰を計算して基準曲面46に投影することで関心対象物の投影画像データ49を作成した後,これを平面に展開することで,パノラマ平面画像25を作成する。一方,請求項3の発明は,新しいボクセル空間24に基準平面55を設定し,基準曲面46が基準平面55になるように湾曲した厚板状領域43を新しいボクセル空間24の非湾曲した厚板状領域52に写像した後,基準平面55の垂線方向57における厚板領域52での光の減衰を計算して基準平面55に投影することで,パノラマ平面画像25を作成する。(以降,CTデータのボクセル空間をCTボクセル空間と略称する。また,新しいボクセル空間をパノラマボクセル空間と略称する。)
【0011】
すなわち請求項3の発明は,CTボクセル空間において湾曲している基準曲面が,パノラマボクセル空間において平面になるように変換を行う。従って,CTボクセル空間における基準曲面に垂直方向の関心オブジェクトでの光の減衰の基準曲面への投影は,パノラマボクセル空間においては基準平面に垂直方向の関心オブジェクトでの光の減衰の基準平面への投影になる。
【0012】
請求項1の発明と請求項3の発明は得られる結果は同一であるが,請求項1の発明は関心オブジェクトでの光の減衰の基準曲面への投影を計算するときに基準曲面の垂直方向が変化するが,請求項3の発明は関心オブジェクトでの光の減衰の基準平面への投影を計算するときに基準平面の垂直方向が変化しない。このため,関心オブジェクトでの光の減衰の投影を求める計算は請求項3の発明の方が容易である。従って,対話的な操作を行う場合には請求項3の発明の方が優れている。
【0013】
請求項3の発明はパノラマボクセル空間に変換した後の関心オブジェクトを取り扱っているので,基準平面だけでなく他の任意の平面への光の減衰の投影処理によってもパノラマ状に展開した画像を得ることができる。また,パノラマボクセル空間に変換した後の関心オブジェクトに,断面変換処理,ボリュームレンダリング処理などの他の三次元画像処理を行うことによってもパノラマ状に展開した画像を得ることができる。
【0014】
【発明の実施の形態】
以下,本発明による三次元画像表示装置について説明する。図1は本発明の実施例を説明するブロック図である。この実施例ではX線CT装置を使用する場合について説明しているが,MR装置など他の装置の場合にも同様に機能する。データ収集部1は,X線CT装置においては被検体の周囲からX線を照射し,透過するX線の線量を測定することによって投影データを収集する部分である。
【0015】
この実施例では電子ビーム・スキャン方式のX線CT装置を例に示している。電子銃12から放射した電子ビーム13を制御して,被検体の周囲に環状に配置されたX線ターゲット11上を走査する。X線ターゲットで発生したX線は,寝台16上に横たわる被検体の横断面を透過した後,X線検出器14で電流に変換される。このX線検出器の出力はデータ収集回路15でデジタルデータに変換される。寝台16を移動させながら電子ビームの走査を繰り返して行うことによって,被検体の複数横断面のデータを収集する。この実施例では電子ビーム・スキャン方式のX線CT装置について説明したが,他の方式のX線CT装置においても同様に機能する。
【0016】
再構成演算処理部2は,収集したデータに前処理を行って投影データを作成する。次に,この投影データに再構成処理,後処理を行って画像データを作成する。この画像データは被検体のX線透過係数に相当するCT値の空間的な分布を表している。再構成演算処理部2は投影データを再構成する機能をもつ高速演算装置を含むデータ処理装置である。データ蓄積装置17は再構成演算処理部2で得た投影データおよび画像データを蓄積する機能を備えている。
【0017】
三次元画像処理部3は三次元画像を再構成する機能をもつデータ処理装置である。三次元画像処理部3は再構成演算処理部2と一体化されている場合もあり,分離されている場合もある。三次元画像処理部3は再構成演算処理部2によって作成された投影データまたは画像データを使用して三次元画像を再構成する。再構成演算処理部2で得た投影データまたは画像データを直接使用する場合もあり,データ蓄積装置17に蓄積した投影データまたは画像データを使用する場合もあり,またオフライン媒体に保存した投影データまたは画像データを使用する場合もある。
【0018】
画像表示部4は,三次元ボクセルデータにパラメタ設定面を設定するCTボクセル表示領域6,パラメタ設定面に関心オブジェクトの湾曲した厚板状領域と基準曲線とを設定する関心領域・基準曲面表示領域7,関心オブジェクトの湾曲した厚板状領域を変換した非湾曲厚板領域を表示するパノラマボクセル表示領域8,得られたパノラマ画像を表示するパノラマ画像表示領域9,三次元画像を表示する三次元画像表示領域10を持ち,画像処理パラメタを設定する機能を備えている。
【0019】
図3は本発明の実施例1を説明するブロック図であり,請求項1に相当する。21はCTボクセル空間で,X線CT装置で撮影した上顎または下顎,上下顎の投影データまたは画像データを使用して構築したCT三次元ボクセルデータ41と関心オブジェクトの領域と基準曲面の決定を行う設定面42を示している。
【0020】
28は関心オブジェクトの領域と基準曲面の決定を行う処理である。CT三次元ボクセルデータ41の適当な横断面を設定面42として選択して,この設定面42上で関心オブジェクトの湾曲領域と画像を投影する基準曲線とを設定する。次に,CT三次元ボクセルデータの異なる適当な横断面を設定面として選択して,この設定面上で関心オブジェクトの湾曲領域と画像を投影する基準曲線を設定する。この設定面の数は複数面であっても,2面であっても,1面であってもよい。
【0021】
幾つかの設定面における湾曲領域を連結することによって関心オブジェクトの湾曲した厚板状領域を決定する。また幾つかの設定面における基準曲線を連結することによって基準曲面を決定する。設定面の数が1面の場合には,例えば,関心オブジェクトの湾曲領域と基準曲線を設定面に平行な他の面に写像することによって関心オブジェクトの湾曲した厚板状領域および基準曲面を作成することができる。
【0022】
22はCTボクセル空間で,43は上顎または下顎,上下顎を含む湾曲した厚板状領域を示している。ここで44は関心オブジェクトの湾曲した厚板状領域の外側曲面,45は関心オブジェクトの湾曲した厚板状領域の内側曲面,46は基準曲面,48は基準曲面41上のある点47の法線ベクトルを模式的に示している。
【0023】
29は湾曲した厚板状領域43の光の減衰の基準曲面46への投影処理を示す。基準曲面46上のある点47の法線ベクトル48に沿って曲面44と曲面45で囲まれた湾曲した厚板状領域を透過する光が受ける減衰をボクセル値から計算して,これを基準曲面46上の点47における値とする。この処理を基準曲面46上の全ての点に対して行う。
【0024】
23はCTボクセル空間で,49は投影処理29によって得られた湾曲した厚板状領域を透過する光が受ける減衰を値として持つ基準曲面画像データを示している。この基準曲面上の点50はその点の法線ベクトル51の方向における湾曲した厚板状領域を透過する光が受ける減衰を値として持つ。同様にこの基準曲面上の全ての点はその点の法線ベクトルの方向における湾曲した厚板状領域を透過する光が受ける減衰を値として持つ。
【0025】
30は湾曲した厚板状領域を透過する光が受ける減衰を値としてもつ基準曲面画像データ49をパノラマ平面画像に展開する処理を示している。
【0026】
25は展開処理30によって得られたパノラマ平面画像を示している。これはCTボクセルデータ41の上顎または下顎,上下顎を含む湾曲した厚板状領域を透過する光が受ける減衰を値としてもつ画像である。これはX線パノラマ撮影に類似した画像である。
【0027】
図4は本発明の実施例2を説明するブロック図であり,請求項2に相当する。図3と同一の番号は図3と同一の機能をもつ。21はX線CT装置で撮影した上顎,下顎または上下顎の投影データまたは画像データを使用して構築したCT三次元ボクセルデータ41と設定面42を示している。
【0028】
28は関心オブジェクトと基準曲面の決定を行う処理である。
【0029】
22は決定された上顎,下顎または上下顎を含む湾曲した厚板状領域43と基準曲面46を示す。
【0030】
29は湾曲した厚板状領域43における光の減衰の基準曲面46への投影処理を示す。
【0031】
23は投影処理29によって得られた湾曲した厚板状領域を透過する光が受ける減衰を値として持つ基準曲面画像データ49を示す。
【0032】
30は湾曲した厚板状領域を透過する光が受ける減衰を値としてもつ基準曲面画像データ49をパノラマ平面画像に展開する処理を示している。
【0033】
25は展開処理30によって得られたパノラマ平面画像を示している。
【0034】
34は湾曲した厚板状領域43への三次元画像処理を示す。三次元画像処理としてはボリュームレンダリング,MPRと呼ばれる断面変換処理などがある。
【0035】
27は湾曲した厚板状領域43の三次元画像処理34によって得られた画像の一例を模式的に示している。
【0036】
35はパノラマ平面画像25上のROIとCTボクセル空間の画像27上のROIとの連携を示す。パノラマ平面画像25上に設定したROIの場所がCTボクセル空間の画像27上の対応する場所にROIとして表示される。また,CTボクセル空間の画像27上に設定したROIの場所がパノラマ平面画像25上の対応する場所にROIとして表示される。これによって,パノラマ平面画像とCTボクセル空間の三次元画像上の位置関係を把握することが容易になった。
【0037】
図5は本発明の実施例3を説明するブロック図であり,請求項3に相当する。図3および図4と同一の番号は図3および図4と同一の機能をもつ。21はX線CT装置で撮影した上顎,下顎または上下顎の投影データまたは画像データを使用して構築したCT三次元ボクセルデータ41と設定面42を示している。
【0038】
28は関心オブジェクトと基準曲面の決定を行う処理である。
【0039】
22は決定された上顎または下顎,上下顎を含む湾曲した厚板状領域43と基準曲面46を示す。
【0040】
31はCTボクセル空間における湾曲した厚板状領域43を,パノラマボクセル空間の非湾曲厚板領域52に変換する処理を示す。CTボクセル空間における基準曲面46をパノラマボクセル空間における基準平面55に写像する関係式を作成する。この関係式に従って,CTボクセル空間における基準曲面46をパノラマボクセル空間における基準平面55に写像する。またこの関係式に従って,CTボクセル空間における湾曲した厚板状領域43を,パノラマボクセル空間の非湾曲厚板領域52に変換する。
【0041】
24はパノラマボクセル空間で,52はCTボクセル空間における湾曲厚板領域43をパノラマボクセル空間に変換した後の上顎または下顎,上下顎を含む非湾曲厚板領域を示す。ここで53は関心オブジェクトの湾曲厚板領域の外側曲面44を変換した曲面または平面,54は関心オブジェクトの湾曲厚板領域の内側曲面45を変換した曲面または平面,55は基準曲面46を変換した基準平面,57は基準平面上のある点56の法線ベクトルを模式的に示している。
【0042】
32は非湾曲厚板領域の基準平面への投影処理を示す。基準平面55上のある点56の法線ベクトル57に沿って曲面または平面53と曲面または平面54で囲まれた厚板領域を透過する光が受ける減衰をボクセル値から計算して,これを基準平面55上の点56における値にする。これを基準平面55上の全ての点に対して行う。
【0043】
25は投影処理32によって得られた非湾曲厚板領域52を透過する光が受ける減衰を値として持つパノラマ平面画像である。この画像はパノラマボクセル空間において,基準平面上の点の法線ベクトル方向における非湾曲厚板領域52を透過する光が受ける減衰を値として持つ画像であるが,これは同時にCTボクセル空間において,基準曲面46上の点の法線ベクトル48方向における湾曲した厚板状領域43を透過する光が受ける減衰を値として持つ画像である。すなわち,これはCTボクセルデータ41の上顎,下顎または上下顎を含む湾曲した厚板状領域を透過する光が受ける減衰を値としてもつ画像である。これはX線パノラマ撮影に類似した画像である。
【0044】
33は厚板領域24に対するMPRなどの断面変換処理,ボリュームレンダリング処理などの三次元画像処理を示す。
【0045】
26は非湾曲厚板領域52に対するMPRなどの断面変換処理,ボリュームレンダリング処理などの三次元画像処理33によって得られた画像を示す。これはパノラマボクセル空間において非湾曲厚板領域52に対して三次元画像処理を行った結果の画像であるが,これはCTボクセル空間において湾曲した厚板状領域43に対して三次元画像処理を行った結果の画像をパノラマ状に展開した画像と同一の画像である。たとえば,非湾曲厚板領域52における基準平面55に平行な複数の平面とその平面に直交する多数の平面を設定してその断面像を得ることによって,その歯列に沿った複数の曲面におけるパノラマ状断面像とその曲面に直交する多数の平面の断面像を作成,表示することが可能になった。この方法で,X線CT画像データを使用して歯列の状態を観察する方法である歯科用曲面再構成で得られる画像と同一の画像が得られる。この方法では,CTボクセル空間における湾曲した厚板状領域43をパノラマボクセル空間における非湾曲厚板柱状領域52に変換した後で三次元画像処理を行うので,対話的にパラメタを変えながらパノラマ状に展開した画像を得ることが可能になった。
【0046】
図6は本発明の実施例4を説明するブロック図であり,請求項4に相当する。図3および図4,図5と同一の番号は図3および図4,図5と同一の機能をもつ。21はX線CT装置で撮影した上顎または下顎,上下顎の投影データまたは画像データを使用して構築したCT三次元ボクセルデータ41と設定面42を示している。
【0047】
28は関心オブジェクトと基準曲面の決定を行う処理である。
【0048】
22は決定された上顎または下顎,上下顎を含む湾曲した厚板状領域43と基準曲面46を示す。
【0049】
31はCTボクセル空間における湾曲した厚板状領域43を,パノラマボクセル空間の非湾曲厚板領域52に変換する処理を示す。
【0050】
24はパノラマボクセル空間で,52はCTボクセル空間における湾曲した厚板状領域43をパノラマボクセル空間に変換した後の非湾曲厚板領域を示す。
【0051】
32は非湾曲厚板領域の基準平面への投影処理を示す。
【0052】
25は投影処理32によって得られた非湾曲厚板領域52を透過する光が受ける減衰を値として持つパノラマ平面画像である。
【0053】
33は厚板領域24に対するMPRなどの断面変換処理,ボリュームレンダリング処理などの三次元画像処理を示す。
【0054】
26は非湾曲厚板領域52に対するMPRなどの断面変換処理,ボリュームレンダリング処理などの三次元画像処理33によって得られた画像を示す。
【0055】
34は湾曲した厚板状領域43への三次元画像処理を示す。三次元画像処理としてはボリュームレンダリング,MPRと呼ばれる断面変換処理などがある。
【0056】
27は湾曲した厚板状領域43の三次元画像処理34によって得られた画像の一例を模式的に示している。
【0057】
35はパノラマ平面画像25上のROIとCTボクセル空間の画像27上のROIとの連携を示す。これによって,パノラマ平面画像とCTボクセル空間の三次元画像上の位置関係を把握することが容易になった。
【0058】
【実施例】
この発明の実施例はX線CT装置を使用してデータを収集する場合について説明したが,MR装置など他の装置を使用してデータを収集する場合にも同様に機能する。
【0059】
この発明の実施例は上顎または下顎,上下顎を含む領域について説明したが,他の湾曲している解剖学的領域に対しても同様に適用することができる。
【0060】
【発明の効果】
本発明の請求項1は,CTボクセル空間に関心オブジェクトを含む湾曲した厚板状領域とこれにほぼ平行する基準曲面とを設定する。そして基準曲面に垂直な光が湾曲した厚板状領域を透過する際に受ける減衰を基準曲面の値として投影した後この基準曲面を平面に展開することによって,湾曲した厚板状領域をもつ関心オブジェクトの性質を一枚の画像によって観察可能にするパノラマ状画像を作成する方法を示した。従来の歯科用曲面再構成法においては断面像を使用しているために体積データの画像化には多数の断面像を必要とするが,本発明では歯列に沿って設定した基準曲面に垂直な方向における関心オブジェクトでの光の減衰を計算して画像を作成するので,関心オブジェクトの厚さ方向の全ての情報を一枚の画像によって,パノラマX線撮影による画像と同様に,画像化することが可能になった。
【0061】
本発明の請求項2は,CTボクセル空間に関心オブジェクトを含む湾曲した厚板状領域とこれにほぼ平行する基準曲面とを設定する。そして基準曲面に垂直な光が湾曲した厚板状領域を透過する際に受ける減衰を基準曲面の値として投影した後この基準曲面を平面に展開することによって,湾曲した厚板状領域をもつ関心オブジェクトの性質を一枚の画像によって観察可能にするパノラマ状画像を作成すると共に,CTボクセル空間の湾曲した厚板状領域に三次元画像処理を行って三次元画像を作成する。そして,このパノラマ状画像上のROIと三次元画像上のROIを連携させることによって解剖学的構造の位置関係を明確にする方法を示した。これによって請求項1の効果のほかに,解剖学的構造の位置関係をより明確にすること可能になった。
【0062】
本発明の請求項3は,CTボクセル空間に関心オブジェクトを含む湾曲した厚板状領域とこれにほぼ平行する基準曲面とを設定する。そして新しいボクセル空間と基準平面を設け,基準曲面を基準平面に写像する変換式を作成する。次に,この変換式に従って,CTボクセル空間の湾曲した厚板状領域を新しいボクセル空間の非湾曲厚板領域に写像する。そして基準平面に垂直な光が非湾曲厚板領域を透過する際に受ける減衰を基準平面の値として投影することによって,湾曲した厚板領域をもつ関心オブジェクトの性質を一枚の画像によって観察可能にするパノラマ状画像を作成する方法を示した。これによって請求項1の効果のほかに,パノラマ状画像のパラメタの設定とその変更を,対話的に行うことが可能になった。
【0063】
本発明の請求項4は,CTボクセル空間に関心オブジェクトを含む湾曲した厚板状領域とこれにほぼ平行する基準曲面とを設定する。そして新しいボクセル空間と基準平面を設け,基準曲面を基準平面に写像する変換式を作成する。次に,この変換式に従って,CTボクセル空間の湾曲した厚板状領域を新しいボクセル空間の非湾曲厚板領域に写像する。そして基準平面に垂直な光が非湾曲厚板領域を透過する際に受ける減衰を基準平面の値として投影することによって,湾曲した厚板領域をもつ関心オブジェクトの性質を一枚の画像によって観察可能にするパノラマ状画像を作成すると共に,CTボクセル空間の湾曲した厚板状領域に三次元画像処理を行って三次元画像を作成する。そして,このパノラマ状画像上のROIと三次元画像上のROIを連携させることによって解剖学的構造の位置関係を明確にする方法を示した。
【0064】
本発明によって,上顎,下顎,上下顎などの湾曲した関心オブジェクトの情報を全て一枚に含んだ画像を平面に展開して表示することが可能になった。従来の歯科用曲面再構成法において多数の断面画像を同時に観察しなければ関心オブジェクトの全体像が把握できないのとは異なり,本発明は一枚の画像によって関心オブジェクトの全体像を把握することができる。従って,対話式の操作と立体構造の直感的な把握を可能にする画像表示装置を提供することができた。
【図面の簡単な説明】
【図1】本発明の実施例の構成を示す概略図。
【図2】本発明を説明するブロック図。
【図3】本発明の実施例1を説明するブロック図。
【図4】本発明の実施例2を説明するブロック図。
【図5】本発明の実施例3を説明するブロック図。
【図6】本発明の実施例4を説明するブロック図。
【図7】従来の歯科用曲面再構成法を説明するブロック図。
【符号の説明】
1 データ収集部
2 再構成演算処理部
3 三次元画像処理部
4 画像表示部
5 操作部
6 CTボクセル表示領域
7 関心領域・基準曲面表示領域
8 パノラマボクセル表示領域
9 パノラマ画像表示領域
10 三次元画像表示領域
11 X線ターゲット
12 電子銃
13 電子ビーム
14 X線検出器
15 データ収集回路
16 寝台
17 データ蓄積装置
21 CTボクセル空間におけるCTボクセルデータと設定面
22 CTボクセル空間における湾曲した厚板状領域と基準曲面
23 CTボクセル空間における基準曲面画像データ
24 パノラマボクセル空間における非湾曲厚板領域と基準平面
25 パノラマ平面画像
26 パノラマ状断面画像
27 三次元画像
28 湾曲した厚板状領域と基準曲面の決定
29 湾曲した厚板状領域の基準曲面への投影
30 基準曲面画像データのパノラマ平面画像への展開
31 湾曲した厚板状領域の非湾曲厚板領域への写像
32 非湾曲厚板領域のパノラマ平面画像への投影
33 パノラマボクセル空間の三次元画像処理
34 CTボクセル空間の三次元画像処理
35 パノラマ平面画像ROIとCTボクセル空間画像ROIの連携
41 CTボクセルデータ
42 設定面
43 湾曲した厚板状領域
44 湾曲した厚板状領域の外側曲面
45 湾曲した厚板状領域の内側曲面
46 基準曲面
47 基準曲面上の点
48 基準曲面の法線ベクトル
49 基準曲面画像データ
50 基準曲面画像データ上の点
51 基準曲面画像データの法線ベクトル
52 パノラマボクセル空間における非湾曲厚板領域
53 パノラマボクセル空間における非湾曲厚板領域の外側曲面
54 パノラマボクセル空間における非湾曲厚板領域の内側曲面
55 パノラマボクセル空間における基準平面
56 パノラマボクセル空間における基準平面上の点
57 パノラマボクセル空間における基準平面の法線ベクトル
61 CTボクセルデータと設定面
62 歯列に平行な曲面と直交する平面の設定
63 歯列に平行な曲面と直交する平面の断面像
64 歯列に平行な曲面と直交する平面を設定する処理
65 歯列に平行な曲面と直交する平面の断面像を作成する処理
66 CTボクセルデータ
67 設定面
68 歯列に平行な曲面
69 歯列に直交する平面
70 歯列に平行な断面像
71 歯列に直交する断面像
[0001]
[Technology to which the invention belongs]
The present invention relates to a three-dimensional image display device that displays the spatial distribution of physical properties of a subject, particularly the properties of curved regions such as the upper or lower jaw and the upper and lower jaws, as an image.
[0002]
[Prior art]
In the dental area, panoramic X-ray photography is performed, which is a photographing method in which the upper jaw, lower jaw, and upper and lower jaws are drawn as a continuous image on a single film. The temporomandibular joint, maxilla and mandible, and dentition can be observed on a single film, so it is widely used in the dental field.
[0003]
Since accurate cross-sectional image data representing the physical properties of the human body can be obtained by the X-ray CT apparatus, a three-dimensional stereoscopic image is reconstructed using a plurality of cross-sectional image data taken at different cross-sectional positions. Has been done. In particular, helical scan X-ray CT devices and cone beam X-ray CT devices have recently been put into practical use, and it has become possible to reconstruct more precise three-dimensional stereoscopic images. Is progressing.
[0004]
Image processing called dental curved surface reconstruction is performed as a method for observing the state of a dentition using data collected by an X-ray CT apparatus. This uses three-dimensional voxel data including the upper jaw, lower jaw, and upper and lower jaws created from data photographed by an X-ray CT apparatus. Set multiple curved surfaces along the dentition and multiple planes orthogonal to the curved surface, create and display panoramic cross-sectional images of multiple curved surfaces along the dentition and multiple planes orthogonal to the curved surface It is a method to do. By using an image processing method called dental curved surface reconstruction, the temporomandibular joint, maxilla and mandible, and dentition sections can be developed and observed in a manner similar to panoramic radiography.
[0005]
FIG. 7 is a block diagram for explaining this conventional dental curved surface reconstruction. Reference numeral 61 denotes CT three-dimensional voxel data 66 constructed using image data of the upper jaw, lower jaw, and upper and lower jaws taken by the X-ray CT apparatus, and a setting surface 67 for setting the sectional position of the sectional image. Reference numeral 64 denotes a process of setting a plurality of curved surfaces along the dentition and a large number of planes orthogonal to the curved surfaces. Reference numeral 62 denotes a plurality of curved surfaces 68 set along the dentition in the upper or lower jaw and upper and lower jaw regions on the setting surface, and a large number of planes 69 orthogonal thereto. Reference numeral 65 denotes processing for creating a cross-sectional image of a curved surface parallel to the dentition and a cross-sectional image of a plane orthogonal to the curved surface. 63 schematically shows an image 70 in which cross-sectional images of a plurality of curved surfaces parallel to the dentition are developed on a plane and cross-sectional images 71 of a number of planes orthogonal to the image 70. In dental curved surface reconstruction, cross-sectional images of curved surfaces parallel to the dentition are created at intervals of about 1 mm, and cross-sectional images of planes perpendicular to this are also created at intervals of about 1 mm. A number of images are created.
[0006]
An image obtained by panoramic X-ray imaging is a single X-ray transmission image, but an image obtained by dental curved surface reconstruction is a cross-sectional image of a plurality of curved surfaces set along the dentition and orthogonal to the curved surface. It is the cross-sectional image of many sheets. Therefore, in panoramic X-ray imaging, it is possible to grasp the information of the entire dentition by observing a single image, but in the dental curved surface reconstruction, in order to grasp the information of the entire dentition, the dentition It is necessary to simultaneously observe a plurality of cross-sectional images of a curved surface set along the line and a plurality of cross-sectional images orthogonal to the curved surface. When all images created by dental curved surface reconstruction are printed on a single film, it is possible to grasp the information of the entire dentition, but it is created while processing interactively on the display screen of the workstation. In the case of displaying the obtained image, in order to observe the entire dentition, a large number of cross-sectional images are sequentially observed, so that it is difficult to intuitively grasp the three-dimensional structure.
[0007]
[Problems to be solved by the invention]
In the conventional curved surface reconstruction method, it is necessary to carefully set the curved surface set along the dentition and the plane perpendicular to this in order to prevent leakage of information, and the interval is narrow. Therefore, the number of cross-sectional images becomes enormous. In the conventional dental curved surface reconstruction method, in order to grasp the information of the entire dentition, the cross-sectional images of a plurality of curved surfaces set along the dentition and a plurality of cross-sectional images orthogonal to the curved surface are observed. Because it is necessary, interactive operation is difficult, and it is difficult to intuitively grasp the three-dimensional structure. Displays a flat development image that includes all information of curved plank-like anatomical regions such as the upper jaw, lower jaw, and upper and lower jaws, and enables interactive operation and intuitive grasp of the three-dimensional structure. An image display device is required.
[0008]
[Means for Solving the Problems]
  The present invention has been devised to remedy the drawbacks of conventional curved surface reconstruction methods. In FIG. 2, the invention of claim 1 includes a three-dimensional voxel data 41 including the upper jaw, lower jaw, and upper and lower jaws, a reference curved surface 46 along the dentition, and a curved plank-like region 43 including the object of interest in the dentition. Is calculated from the voxel pixel value, and the result is projected onto the reference curved surface 46, by calculating the light attenuation in the curved thick plate region 43 in the perpendicular direction 48 of the reference curved surface 46. Means 29 for creating the projection image data 49 of the object of interest on the curved surface 46, and means 30 for creating a relational expression for developing the reference curved surface 46 on a plane and developing the image data 49 on the reference curved surface on the plane according to this relational expression The panorama plane image 25 is created and displayed. The conventional dental curved surface reconstruction method uses a cross-sectional image, so that a large number of cross-sectional images are required for imaging volume data. In the present invention, however, the vertical surface is perpendicular to the reference curved surface set along the dentition. Since the image is created by calculating the attenuation of light at the object of interest in any direction, it is possible to image all the information of the object of interest with a single image, as with panoramic X-ray images. is there.
[0009]
  In the invention of claim 3, in FIG. 2, three-dimensional voxel data 41 including the upper jaw or lower jaw, upper and lower jaws, a reference curved surface 46 along the dentition, and a curved plank region 43 including the object of interest of the dentition, Is prepared, a new three-dimensional voxel space 24 and a reference plane 55 are prepared, a relational expression for converting the reference curved surface 46 to the reference plane 55 is created, and a curved thick plate based on the reference curved surface according to this relational expression Means 31 for creating a thick plate-like object of interest 52 converted by mapping the voxel data of the region 43 into voxels in a new three-dimensional voxel space based on the reference plane, and in a direction perpendicular to the reference plane 55 Means 32 for calculating the attenuation of light at the object of interest in the non-curved plank region 52 from the voxel pixel values and projecting the result onto the reference plane 55. And by, and to display to create a panorama-like planar image 25.
[0010]
  In the present invention, the reference surface 46 along the dentition and the object of interest are included in the voxel space of the CT data.Curved plank-like region43 is set. In the first aspect of the present invention, the reference curved surface 46 is perpendicular to the vertical direction 48.Curved plank-like regionThe projection image data 49 of the object of interest is created by calculating the light attenuation at 43 and projecting it onto the reference curved surface 46, and then developing this onto a plane to create the panorama plane image 25. On the other hand, the invention of claim 3 sets the reference plane 55 in the new voxel space 24 so that the reference curved surface 46 becomes the reference plane 55.Curved plank-like region43 is mapped onto the non-curved thick plate-like region 52 of the new voxel space 24, and then the attenuation of light in the thick plate region 52 in the normal direction 57 of the reference plane 55 is calculated and projected onto the reference plane 55, A panoramic plane image 25 is created. (Hereafter, the voxel space of CT data is abbreviated as CT voxel space. The new voxel space is abbreviated as panoramic voxel space.)
[0011]
That is, the invention according to claim 3 performs conversion so that the reference curved surface curved in the CT voxel space becomes a plane in the panoramic voxel space. Therefore, the projection of the light attenuation of the object of interest perpendicular to the reference curved surface in the CT voxel space onto the reference curved surface in the panoramic voxel space is directed to the reference plane of the light attenuation of the object of interest perpendicular to the reference plane. Become a projection.
[0012]
The invention of claim 1 and the invention of claim 3 obtain the same result, but the invention of claim 1 is the vertical direction of the reference curved surface when calculating the projection of the attenuation of light on the object of interest onto the reference curved surface. However, the invention of claim 3 does not change the vertical direction of the reference plane when calculating the projection of the attenuation of light on the object of interest onto the reference plane. Therefore, the calculation for obtaining the projection of the attenuation of light on the object of interest is easier in the invention of claim 3. Therefore, the invention of claim 3 is superior when performing interactive operation.
[0013]
Since the invention of claim 3 handles the object of interest after conversion to the panoramic voxel space, an image developed in a panoramic shape is obtained not only by the reference plane but also by projection processing of attenuation of light on any other plane. be able to. A panoramic image can also be obtained by performing other three-dimensional image processing such as cross-section conversion processing and volume rendering processing on the object of interest after conversion to the panoramic voxel space.
[0014]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, a three-dimensional image display apparatus according to the present invention will be described. FIG. 1 is a block diagram illustrating an embodiment of the present invention. In this embodiment, the case where an X-ray CT apparatus is used has been described, but the same functions in the case of another apparatus such as an MR apparatus. In the X-ray CT apparatus, the data collection unit 1 is a part that collects projection data by irradiating X-rays from around the subject and measuring the dose of transmitted X-rays.
[0015]
In this embodiment, an electron beam scanning type X-ray CT apparatus is shown as an example. The electron beam 13 emitted from the electron gun 12 is controlled to scan the X-ray target 11 arranged in a ring around the subject. X-rays generated by the X-ray target pass through the cross section of the subject lying on the bed 16 and are then converted into current by the X-ray detector 14. The output of the X-ray detector is converted into digital data by the data acquisition circuit 15. By scanning the electron beam repeatedly while moving the bed 16, data of a plurality of cross sections of the subject is collected. In this embodiment, the electron beam scanning type X-ray CT apparatus has been described. However, other types of X-ray CT apparatuses function in the same manner.
[0016]
The reconstruction calculation processing unit 2 preprocesses the collected data to create projection data. Next, image data is created by performing reconstruction processing and post-processing on the projection data. This image data represents a spatial distribution of CT values corresponding to the X-ray transmission coefficient of the subject. The reconstruction calculation processing unit 2 is a data processing device including a high-speed calculation device having a function of reconstructing projection data. The data storage device 17 has a function of storing projection data and image data obtained by the reconstruction calculation processing unit 2.
[0017]
The three-dimensional image processing unit 3 is a data processing device having a function of reconstructing a three-dimensional image. The three-dimensional image processing unit 3 may be integrated with the reconstruction calculation processing unit 2 or may be separated. The 3D image processing unit 3 reconstructs a 3D image using the projection data or image data created by the reconstruction calculation processing unit 2. The projection data or image data obtained by the reconstruction calculation processing unit 2 may be used directly, the projection data or image data stored in the data storage device 17 may be used, or the projection data stored in the offline medium or Sometimes image data is used.
[0018]
  The image display unit 4 includes a CT voxel display area 6 for setting a parameter setting surface in the three-dimensional voxel data, and an object of interest on the parameter setting surface.Curved plank-like regionRegion of interest and reference curved surface display area 7 for setting the object and reference curve,Curved plank-like regionA panoramic voxel display area 8 for displaying the non-curved thick plate area converted from the above, a panoramic image display area 9 for displaying the obtained panoramic image, and a three-dimensional image display area 10 for displaying a three-dimensional image. It has a function to set.
[0019]
FIG. 3 is a block diagram for explaining the first embodiment of the present invention and corresponds to claim 1. Reference numeral 21 denotes a CT voxel space, which determines CT three-dimensional voxel data 41 constructed by using projection data or image data of the upper jaw or lower jaw, upper and lower jaws taken by an X-ray CT apparatus, an object region of interest, and a reference curved surface. A setting surface 42 is shown.
[0020]
Reference numeral 28 denotes a process for determining the region of interest and the reference curved surface. An appropriate cross section of the CT three-dimensional voxel data 41 is selected as the setting plane 42, and a curved region of the object of interest and a reference curve for projecting an image are set on the setting plane 42. Next, an appropriate cross section having different CT three-dimensional voxel data is selected as a setting plane, and a curved region of the object of interest and a reference curve for projecting an image are set on the setting plane. The number of setting surfaces may be a plurality of surfaces, two surfaces, or one surface.
[0021]
  By connecting curved areas in several setting planes,Curved plank-like regionTo decide. The reference curved surface is determined by connecting the reference curves on several setting surfaces. When the number of setting planes is one, for example, the curved area of the object of interest and the reference curve are mapped to another plane parallel to the setting plane.Curved plank-like regionAnd a reference curved surface can be created.
[0022]
  22 is CT voxel space, 43 includes maxilla or mandible and maxillaCurved plank-like regionIs shown. Where 44 is the object of interestCurved plank-like region, 45 is the object of interestCurved plank-like region, 46 is a reference curved surface, and 48 is a normal vector of a point 47 on the reference curved surface 41.
[0023]
  29 isCurved plank-like regionThe projection processing of 43 light attenuation onto the reference curved surface 46 is shown. Surrounded by a curved surface 44 and a curved surface 45 along a normal vector 48 of a point 47 on the reference curved surface 46Curved plank-like regionIs calculated from the voxel value, and this is set as a value at a point 47 on the reference curved surface 46. This process is performed for all points on the reference curved surface 46.
[0024]
  23 is a CT voxel space, 49 is obtained by the projection process 29Curved plank-like regionThe reference curved surface image data having attenuation as a value received by light passing through is shown. The point 50 on the reference curved surface is in the direction of the normal vector 51 of the point.Curved plank-like regionAttenuation received by the light passing through is taken as a value. Similarly, all points on this reference surface are in the direction of their normal vectorsCurved plank-like regionAttenuation received by the light passing through is taken as a value.
[0025]
  30 isCurved plank-like regionThe reference curved surface image data 49 having, as a value, the attenuation received by the light passing through is developed into a panoramic plane image.
[0026]
  Reference numeral 25 denotes a panoramic plane image obtained by the expansion process 30. This includes the upper or lower jaw, upper and lower jaws of CT voxel data 41Curved plank-like regionIt is an image which has the attenuation which the light which permeate | transmits receives. This is an image similar to X-ray panoramic imaging.
[0027]
FIG. 4 is a block diagram for explaining a second embodiment of the present invention, and corresponds to claim 2. The same numbers as those in FIG. 3 have the same functions as those in FIG. Reference numeral 21 denotes CT three-dimensional voxel data 41 and a setting surface 42 constructed using projection data or image data of the upper jaw, lower jaw or upper and lower jaws taken by an X-ray CT apparatus.
[0028]
Reference numeral 28 denotes a process for determining an object of interest and a reference curved surface.
[0029]
  22 includes determined maxilla, mandible or maxillaCurved plank-like region43 and a reference curved surface 46 are shown.
[0030]
  29 isCurved plank-like regionThe projection processing of the light attenuation at 43 onto the reference curved surface 46 is shown.
[0031]
  23 was obtained by the projection process 29Curved plank-like regionThe reference curved surface image data 49 having attenuation as a value received by light passing through is shown.
[0032]
  30 isCurved plank-like regionThe reference curved surface image data 49 having, as a value, the attenuation received by the light passing through is developed into a panoramic plane image.
[0033]
Reference numeral 25 denotes a panoramic plane image obtained by the expansion process 30.
[0034]
  34 isCurved plank-like region3D image processing to 43 is shown. Three-dimensional image processing includes volume rendering, cross-sectional conversion processing called MPR, and the like.
[0035]
  27 isCurved plank-like regionAn example of an image obtained by 43 three-dimensional image processing 34 is schematically shown.
[0036]
Reference numeral 35 denotes the linkage between the ROI on the panorama plane image 25 and the ROI on the image 27 in the CT voxel space. The ROI location set on the panorama plane image 25 is displayed as the ROI at the corresponding location on the image 27 in the CT voxel space. Further, the location of the ROI set on the image 27 in the CT voxel space is displayed as the ROI at the corresponding location on the panorama plane image 25. This makes it easy to grasp the positional relationship between the panoramic planar image and the CT voxel space on the three-dimensional image.
[0037]
FIG. 5 is a block diagram for explaining a third embodiment of the present invention, and corresponds to claim 3. The same numbers as those in FIGS. 3 and 4 have the same functions as those in FIGS. Reference numeral 21 denotes CT three-dimensional voxel data 41 and a setting surface 42 constructed using projection data or image data of the upper jaw, lower jaw or upper and lower jaws taken by an X-ray CT apparatus.
[0038]
Reference numeral 28 denotes a process for determining an object of interest and a reference curved surface.
[0039]
  22 includes determined upper or lower jaw, upper and lower jawCurved plank-like region43 and a reference curved surface 46 are shown.
[0040]
  31 in CT voxel spaceCurved plank-like region43 shows a process of converting 43 into a non-curved thick plate region 52 of a panoramic voxel space. A relational expression for mapping the reference curved surface 46 in the CT voxel space to the reference plane 55 in the panoramic voxel space is created. According to this relational expression, the reference curved surface 46 in the CT voxel space is mapped to the reference plane 55 in the panoramic voxel space. Also, according to this relational expression, in CT voxel spaceCurved plank-like region43 is converted into a non-curved thick plate region 52 of the panoramic voxel space.
[0041]
Reference numeral 24 denotes a panoramic voxel space, and 52 denotes a non-curved slab area including the upper jaw, the lower jaw, and the upper and lower jaws after the curved slab area 43 in the CT voxel space is converted into the panoramic voxel space. Here, 53 is a curved surface or plane obtained by converting the outer curved surface 44 of the curved plank region of the object of interest, 54 is a curved surface or plane transformed by the inner curved surface 45 of the curved plank region of the object of interest, and 55 is a converted standard curved surface 46. Reference plane 57, schematically shows a normal vector of a certain point 56 on the reference plane.
[0042]
Reference numeral 32 denotes a process of projecting the non-curved thick plate region onto the reference plane. Attenuation received by light passing through a curved plate or plane 53 and a plank region surrounded by the curved surface or plane 54 along a normal vector 57 of a certain point 56 on the reference plane 55 is calculated from the voxel value, and this is used as a reference. The value at the point 56 on the plane 55 is set. This is performed for all points on the reference plane 55.
[0043]
  Reference numeral 25 denotes a panoramic plane image having the attenuation received by the light transmitted through the non-curved thick plate region 52 obtained by the projection processing 32 as a value. In the panoramic voxel space, this image has an attenuation value received by the light transmitted through the non-curved slab region 52 in the normal vector direction of the point on the reference plane. The point on the curved surface 46 in the direction of the normal vector 48Curved plank-like region43 is an image having attenuation as a value received by light transmitted through 43. That is, this includes the upper jaw, lower jaw or upper and lower jaws of CT voxel data 41Curved plank-like regionIt is an image which has the attenuation which the light which permeate | transmits receives. This is an image similar to X-ray panoramic imaging.
[0044]
Reference numeral 33 denotes three-dimensional image processing such as MPR or the like cross-section conversion processing or volume rendering processing for the thick plate region 24.
[0045]
  Reference numeral 26 denotes an image obtained by three-dimensional image processing 33 such as MPR or the like cross-section conversion processing or volume rendering processing for the non-curved thick plate region 52. This is an image obtained as a result of performing the three-dimensional image processing on the non-curved plank region 52 in the panoramic voxel space.Curved plank-like region43 is the same image as the image obtained by performing the three-dimensional image processing on 43 in a panoramic shape. For example, by setting a plurality of planes parallel to the reference plane 55 in the non-curved plank region 52 and a number of planes orthogonal to the planes, and obtaining cross-sectional images thereof, panoramas on the plurality of curved surfaces along the dentition are obtained. It has become possible to create and display cross-sectional images of a plane and a number of planes orthogonal to the curved surface. By this method, the same image as that obtained by dental curved surface reconstruction, which is a method of observing the state of the dentition using X-ray CT image data, is obtained. In this method, CT voxel spaceCurved plank-like regionSince the three-dimensional image processing is performed after converting 43 into the non-curved thick plate columnar region 52 in the panoramic voxel space, it is possible to obtain a panoramic image while changing parameters interactively.
[0046]
FIG. 6 is a block diagram for explaining a fourth embodiment of the present invention, and corresponds to claim 4. The same reference numerals as those in FIGS. 3, 4, and 5 have the same functions as those in FIGS. 3, 4, and 5. Reference numeral 21 denotes CT three-dimensional voxel data 41 and a setting surface 42 constructed using projection data or image data of the upper jaw or the lower jaw and the upper and lower jaws taken by the X-ray CT apparatus.
[0047]
Reference numeral 28 denotes a process for determining an object of interest and a reference curved surface.
[0048]
  22 includes determined upper or lower jaw, upper and lower jawCurved plank-like region43 and a reference curved surface 46 are shown.
[0049]
  31 in CT voxel spaceCurved plank-like region43 shows a process of converting 43 into a non-curved thick plate region 52 of a panoramic voxel space.
[0050]
  24 is a panoramic voxel space and 52 is a CT voxel space.Curved plank-like regionThe non-curved slab area after converting 43 into panoramic voxel space is shown.
[0051]
Reference numeral 32 denotes a process of projecting the non-curved thick plate region onto the reference plane.
[0052]
Reference numeral 25 denotes a panoramic plane image having the attenuation received by the light transmitted through the non-curved thick plate region 52 obtained by the projection processing 32 as a value.
[0053]
Reference numeral 33 denotes three-dimensional image processing such as MPR or the like cross-section conversion processing or volume rendering processing for the thick plate region 24.
[0054]
Reference numeral 26 denotes an image obtained by three-dimensional image processing 33 such as MPR or the like cross-section conversion processing or volume rendering processing for the non-curved thick plate region 52.
[0055]
  34 isCurved plank-like region3D image processing to 43 is shown. Three-dimensional image processing includes volume rendering, cross-sectional conversion processing called MPR, and the like.
[0056]
  27 isCurved plank-like regionAn example of an image obtained by 43 three-dimensional image processing 34 is schematically shown.
[0057]
Reference numeral 35 denotes the linkage between the ROI on the panorama plane image 25 and the ROI on the image 27 in the CT voxel space. This makes it easy to grasp the positional relationship between the panoramic planar image and the CT voxel space on the three-dimensional image.
[0058]
【Example】
Although the embodiment of the present invention has been described with respect to the case where data is collected using an X-ray CT apparatus, it functions similarly when data is collected using another apparatus such as an MR apparatus.
[0059]
Although the embodiment of the present invention has been described for the region including the upper jaw, the lower jaw, and the upper and lower jaws, it can be similarly applied to other curved anatomical regions.
[0060]
【The invention's effect】
  Claim 1 of the present invention includes an object of interest in a CT voxel space.Curved plank-like regionAnd a reference curved surface substantially parallel to this. And light perpendicular to the reference curved surfaceCurved plank-like regionBy projecting the attenuation received when passing through as the value of the reference surface, and expanding the reference surface to a plane,Curved plank-like regionA method for creating a panoramic image that allows observation of the properties of an object of interest with a single image is presented. The conventional dental curved surface reconstruction method uses a cross-sectional image, so that a large number of cross-sectional images are required for imaging volume data. In the present invention, however, the vertical surface is perpendicular to the reference curved surface set along the dentition. Since the image is created by calculating the attenuation of light at the object of interest in any direction, all the information in the thickness direction of the object of interest is imaged as a panoramic X-ray image with a single image. It became possible.
[0061]
  Claim 2 of the present invention includes an object of interest in the CT voxel space.Curved plank-like regionAnd a reference curved surface substantially parallel to this. And light perpendicular to the reference curved surfaceCurved plank-like regionBy projecting the attenuation received when passing through as the value of the reference surface, and expanding the reference surface to a plane,Curved plank-like regionCreate a panoramic image that allows observation of the nature of the object of interest with a single image, and the CT voxel spaceCurved plank-like region3D image processing is performed to create a 3D image. And the method of clarifying the positional relationship of the anatomical structure by linking the ROI on the panoramic image and the ROI on the three-dimensional image was shown. This makes it possible to clarify the positional relationship of the anatomical structure in addition to the effect of the first aspect.
[0062]
  Claim 3 of the present invention includes an object of interest in the CT voxel space.Curved plank-like regionAnd a reference curved surface substantially parallel to this. Then, a new voxel space and a reference plane are provided, and a conversion formula for mapping the reference curved surface to the reference plane is created. Next, according to this conversion formula, the CT voxel spaceCurved plank-like regionIs mapped to the non-curved slab area of the new voxel space. And by projecting the attenuation received when light perpendicular to the reference plane passes through the non-curved plank region as the value of the reference plane,Curved plank areaA method for creating a panoramic image that allows observation of the properties of an object of interest with a single image is presented. This makes it possible to interactively set and change the parameters of the panoramic image in addition to the effect of claim 1.
[0063]
  Claim 4 of the present invention includes an object of interest in the CT voxel space.Curved plank-like regionAnd a reference curved surface substantially parallel to this. Then, a new voxel space and a reference plane are provided, and a conversion formula for mapping the reference curved surface to the reference plane is created. Next, according to this conversion formula, the CT voxel spaceCurved plank-like regionIs mapped to the non-curved slab area of the new voxel space. And by projecting the attenuation received when light perpendicular to the reference plane passes through the non-curved plank region as the value of the reference plane,Curved plank areaCreate a panoramic image that allows observation of the nature of the object of interest with a single image, and the CT voxel spaceCurved plank-like region3D image processing is performed to create a 3D image. And the method of clarifying the positional relationship of the anatomical structure by linking the ROI on the panoramic image and the ROI on the three-dimensional image was shown.
[0064]
According to the present invention, it is possible to display an image including all information of curved objects of interest such as the upper jaw, the lower jaw, and the upper and lower jaws on a flat surface. Unlike the conventional method for reconstructing a curved surface of a dental surface, the entire image of an object of interest cannot be grasped unless a large number of cross-sectional images are observed at the same time. it can. Therefore, it has been possible to provide an image display device that enables interactive operation and intuitive grasp of the three-dimensional structure.
[Brief description of the drawings]
FIG. 1 is a schematic diagram showing the configuration of an embodiment of the present invention.
FIG. 2 is a block diagram illustrating the present invention.
FIG. 3 is a block diagram illustrating Example 1 of the present invention.
FIG. 4 is a block diagram illustrating Embodiment 2 of the present invention.
FIG. 5 is a block diagram illustrating Example 3 of the present invention.
FIG. 6 is a block diagram illustrating Embodiment 4 of the present invention.
FIG. 7 is a block diagram for explaining a conventional dental curved surface reconstruction method.
[Explanation of symbols]
1 Data collection department
2 Reconstruction calculation processing unit
3 3D image processor
4 Image display area
5 Operation part
6 CT voxel display area
7 Area of interest / Reference curved surface display area
8 Panorama voxel display area
9 Panorama image display area
10 3D image display area
11 X-ray target
12 electron gun
13 Electron beam
14 X-ray detector
15 Data collection circuit
16 sleeper
17 Data storage device
21 CT voxel data and setting surface in CT voxel space
22 In CT voxel spaceCurved plank-like regionAnd reference curved surface
23 Reference curved surface image data in CT voxel space
24 Non-curved slab area and reference plane in panoramic voxel space
25 panoramic plane image
26 Panoramic cross-sectional image
27 3D image
28Curved plank-like regionAnd reference surface determination
29Curved plank areaProjection of a region onto a reference curved surface
30 Development of reference curved surface image data to panorama plane image
31Curved plank-like regionTo the non-curved plank region
32 Projection of non-curved plank area to panoramic plane image
33 3D image processing of panoramic voxel space
34 Three-dimensional image processing of CT voxel space
35 Cooperation between panoramic planar image ROI and CT voxel space image ROI
41 CT voxel data
42 Setting surface
43Curved plank-like region
44Curved plank-like regionOutside curved surface
45Curved plank-like regionInside curved surface
46 Reference curved surface
47 Points on the reference curved surface
48 Normal vector of the reference surface
49 Reference curved surface image data
50 Points on the reference curved surface image data
51 Normal vector of reference curved surface image data
52 Uncurved slab area in panoramic voxel space
53 Outside curved surface of non-curved thick plate area in panoramic voxel space
54 Inside curved surface of non-curved thick plate area in panoramic voxel space
55 Reference plane in panoramic voxel space
56 Points on the reference plane in panoramic voxel space
57 Normal vector of the reference plane in panoramic voxel space
61 CT voxel data and setting surface
62 Setting a plane perpendicular to the curved surface parallel to the dentition
63 Cross-sectional image of a plane perpendicular to the curved surface parallel to the dentition
64 Processing to set a plane orthogonal to the curved surface parallel to the dentition
65 Processing to create a cross-sectional image of a plane perpendicular to the curved surface parallel to the dentition
66 CT voxel data
67 Setting surface
68 Curved surface parallel to dentition
69 Plane perpendicular to dentition
70 Cross-sectional image parallel to the dentition
71 Cross-sectional image orthogonal to the dentition

Claims (4)

被検体の物理的性質に対応した画素値を持つ三次元ボクセルを使用して三次元画像を表示する装置において,上顎または下顎,上下顎を含む三次元ボクセルデータに歯列に沿った基準曲面と歯列の関心オブジェクトを含む湾曲した厚板状領域とを設定する手段と,この基準曲面の垂直方向における関心オブジェクトでの光の減衰をボクセルの画素値から計算してその結果を基準曲面上に投影することによって基準曲面上に関心オブジェクトの投影画像データを得る手段と,基準曲面上の画像データを平面に展開する手段とを備えることによって,上顎または下顎,上下顎における歯列の状態をパノラマ状に画像表示することを可能にした三次元画像表示装置。  In a device that displays 3D images using 3D voxels with pixel values corresponding to the physical properties of the subject, 3D voxel data including the maxilla, mandible, and maxilla is used as a reference curved surface along the dentition. A means for setting a curved slab-like region containing the object of interest in the dentition, and calculating the light attenuation at the object of interest in the vertical direction of this reference curved surface from the voxel pixel values and putting the result on the reference curved surface By providing means for obtaining projection image data of an object of interest on the reference curved surface by projecting and means for expanding the image data on the reference curved surface into a plane, the state of the dentition in the upper or lower jaw and the upper and lower jaws is panoramic. A three-dimensional image display device that can display an image in a shape. 被検体の物理的性質に対応した画素値を持つ三次元ボクセルを使用して三次元画像を表示する装置において,上顎または下顎,上下顎を含む三次元ボクセルデータに歯列に沿った基準曲面と歯列の関心オブジェクトを含む湾曲した厚板状領域とを設定する手段と,この基準曲面の垂直方向における関心オブジェクトでの光の減衰をボクセルの画素値から計算してその結果を基準曲面上に投影することによって基準曲面上に関心オブジェクトの投影画像データを得る手段と,基準曲面上の画像データを平面に展開する手段とを備えることによって,上顎または下顎,上下顎における歯列の状態をパノラマ状に第一の表示領域に画像表示すると共に,湾曲した厚板状領域に対してボリュウムレンダリングの三次元画像表示法による処理を行った三次元画像を第二の表示領域に表示する手段と,この第一表示領域に設定した関心領域に対応する領域を第二の表示領域に表示するかまたはこの第二表示領域に設定した関心領域に対応する領域を第一表示領域に表示する手段とを備えることによって,パノラマ状画像表示と三次元画像表示の相互関係を明らかにする機能を持たせた三次元表示装置。  In a device that displays 3D images using 3D voxels with pixel values corresponding to the physical properties of the subject, 3D voxel data including the maxilla, mandible, and maxilla is used as a reference curved surface along the dentition. A means for setting a curved slab-like region containing the object of interest in the dentition, and calculating the light attenuation at the object of interest in the vertical direction of this reference curved surface from the voxel pixel values and putting the result on the reference curved surface By providing means for obtaining projection image data of an object of interest on the reference curved surface by projecting and means for expanding the image data on the reference curved surface into a plane, the state of the dentition in the upper or lower jaw and the upper and lower jaws is panoramic. The image is displayed in the first display area, and the curved thick plate area is processed by the volume rendering three-dimensional image display method. A means for displaying the original image in the second display area, and an area corresponding to the area of interest set in the first display area is displayed in the second display area or in the area of interest set in the second display area. A three-dimensional display device having a function of clarifying the interrelation between panoramic image display and three-dimensional image display by providing means for displaying a corresponding region in a first display region. 被検体の物理的性質に対応した画素値を持つ三次元ボクセルを使用して三次元画像を表示する装置において,上顎または下顎,上下顎を含む三次元ボクセルデータに歯列に沿った基準曲面と歯列の関心オブジェクトを含む湾曲した厚板状領域とを設定する手段と,基準曲面が平面になるように湾曲した厚板状領域のボクセルを新しい三次元ボクセル空間の厚板領域のボクセルに写像する手段と,基準平面に垂直な方向におけるこの厚板領域の関心オブジェクトでの光の減衰をボクセルの画素値から計算してその結果を基準平面上に投影する手段とを備えることによって,上顎または下顎,上下顎における歯列の状態をパノラマ状に画像表示することを可能にした三次元画像表示装置。  In a device that displays 3D images using 3D voxels with pixel values corresponding to the physical properties of the subject, 3D voxel data including the maxilla, mandible, and maxilla is used as a reference curved surface along the dentition. A means for setting a curved plank-like region containing the object of interest in the dentition and a voxel in the plank-like region curved so that the reference curved surface becomes a plane is mapped to a voxel in the plank region of the new three-dimensional voxel space Means for calculating the light attenuation at the object of interest in this plank region in a direction perpendicular to the reference plane from the voxel pixel values and projecting the result onto the reference plane, A three-dimensional image display device that can display panoramic images of the dentition in the lower and upper jaws. 被検体の物理的性質に対応した画素値を持つ三次元ボクセルを使用して三次元画像を表示する装置において,上顎または下顎,上下顎を含む三次元ボクセルデータに歯列に沿った基準曲面と歯列の関心オブジェクトを含む湾曲した厚板状領域とを設定する手段と,基準曲面が平面になるように湾曲した厚板状領域のボクセルを新しい三次元ボクセル空間の厚板領域のボクセルに写像する手段と,基準平面に垂直な方向におけるこの厚板領域の関心オブジェクトでの光の減衰をボクセルの画素値から計算してその結果を基準平面上に投影する手段とを備えることによって,上顎または下顎,上下顎における歯列の状態をパノラマ状に第一の表示領域に画像表示すると共に,湾曲した厚板状領域に対してボリュウムレンダリングの三次元画像表示法による処理を行った三次元画像を第二の表示領域に表示する手段と,この第一表示領域に設定した関心領域に対応する領域を第二の表示領域に表示するかまたはこの第二表示領域に設定した関心領域に対応する領域を第一表示領域に表示する手段とを備えることによって,パノラマ状画像表示と三次元画像表示の相互関係を明らかにする機能を持たせた三次元表示装置。  In a device that displays 3D images using 3D voxels with pixel values corresponding to the physical properties of the subject, 3D voxel data including the maxilla, mandible, and maxilla is used as a reference curved surface along the dentition. A means for setting a curved plank-like region containing the object of interest in the dentition and a voxel in the plank-like region curved so that the reference curved surface becomes a plane is mapped to a voxel in the plank region of the new three-dimensional voxel space Means for calculating the light attenuation at the object of interest in this plank region in a direction perpendicular to the reference plane from the voxel pixel values and projecting the result onto the reference plane, The state of the dentition in the lower and upper jaws is displayed in a panoramic image in the first display area, and volume rendering 3D image display for the curved thick plate-like area Means for displaying in the second display area the 3D image that has undergone the processing according to the above, and displaying the area corresponding to the region of interest set in the first display area in the second display area or in the second display area A three-dimensional display device having a function of clarifying the interrelationship between the panoramic image display and the three-dimensional image display by including means for displaying a region corresponding to the region of interest set in the first display region.
JP26427899A 1999-08-16 1999-08-16 3D image display device Expired - Lifetime JP4458309B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP26427899A JP4458309B2 (en) 1999-08-16 1999-08-16 3D image display device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP26427899A JP4458309B2 (en) 1999-08-16 1999-08-16 3D image display device

Publications (2)

Publication Number Publication Date
JP2001061834A JP2001061834A (en) 2001-03-13
JP4458309B2 true JP4458309B2 (en) 2010-04-28

Family

ID=17400956

Family Applications (1)

Application Number Title Priority Date Filing Date
JP26427899A Expired - Lifetime JP4458309B2 (en) 1999-08-16 1999-08-16 3D image display device

Country Status (1)

Country Link
JP (1) JP4458309B2 (en)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100440492B1 (en) * 2002-03-18 2004-07-15 신종계 Method for expanding a curved surface to a plane
US7873192B1 (en) 2005-05-20 2011-01-18 Imaging Sciences International Llc Method for providing a smoothed image
JP4594844B2 (en) * 2005-11-04 2010-12-08 朝日レントゲン工業株式会社 Cone beam X-ray CT imaging system for head and neck
JP4823780B2 (en) * 2006-06-19 2011-11-24 株式会社日立メディコ Panoramic tomographic image generation apparatus and panoramic tomographic image generation program
DE102008009643B4 (en) 2007-02-22 2019-01-24 J. Morita Manufacturing Corporation Image processing method, image processing apparatus and X-ray imaging apparatus for generating a panoramic tomogram.
JP2008229322A (en) * 2007-02-22 2008-10-02 Morita Mfg Co Ltd Image processing method, image displaying method, image processing program, storage medium, image processor, and x-ray imaging device
JP4786685B2 (en) * 2007-11-16 2011-10-05 株式会社モリタ製作所 X-ray image display method, X-ray imaging apparatus, and X-ray image display apparatus
DE102008008733A1 (en) 2008-02-12 2009-08-13 Sirona Dental Systems Gmbh Method for creating a tomogram
JP5396657B2 (en) * 2008-07-17 2014-01-22 株式会社吉田製作所 Composite X-ray imaging method and imaging apparatus
JP5308739B2 (en) * 2008-08-22 2013-10-09 株式会社モリタ製作所 X-ray image display method and X-ray imaging apparatus
JP5337472B2 (en) * 2008-12-24 2013-11-06 イマグノーシス株式会社 Method for displaying cut surface image in medical three-dimensional image, medical three-dimensional image display device, and medical image display program
DK2459070T3 (en) * 2009-07-31 2020-06-02 Dental Imaging Technologies Corp Panoramic tooth imaging using segmentation and a master arc
JP5881192B2 (en) * 2014-06-24 2016-03-09 学校法人日本大学 Panoramic X-ray imaging apparatus and panoramic X-ray imaging method
CN114143602B (en) * 2021-12-06 2023-07-07 深圳创维数字技术有限公司 Panoramic picture playing method and device, panoramic playing server and storage medium

Also Published As

Publication number Publication date
JP2001061834A (en) 2001-03-13

Similar Documents

Publication Publication Date Title
Dobbins III et al. Digital x-ray tomosynthesis: current state of the art and clinical potential
JP4079632B2 (en) Method and apparatus for selecting and displaying medical image data
Scarfe et al. Maxillofacial cone beam computed tomography: essence, elements and steps to interpretation
JP6007386B2 (en) Data processing device for radiation imaging
Ayoub et al. Towards building a photo-realistic virtual human face for craniomaxillofacial diagnosis and treatment planning
US5872828A (en) Tomosynthesis system for breast imaging
JP4495926B2 (en) X-ray stereoscopic reconstruction processing apparatus, X-ray imaging apparatus, X-ray stereoscopic reconstruction processing method, and X-ray stereoscopic imaging auxiliary tool
JP4340334B2 (en) Method and apparatus for determining a conversion between an object and a three-dimensional representation of the object
JP4458309B2 (en) 3D image display device
CN104066376B (en) Equipment for digital radiation photograph and method
JP2008532612A (en) Image processing system and method for alignment of two-dimensional information with three-dimensional volume data during medical treatment
JPH1057365A (en) X-ray imaging method
JPH04144548A (en) Digital panorama x-ray photographing device
US11024061B2 (en) Apparatus and method for scattered radiation correction
JP2004320771A (en) Method for performing digital subtraction angiography
JP3897925B2 (en) Cone beam CT system
JP2006187608A (en) Method and configuration for multiresolutive reconstruction for medical x-ray imaging
Baumrind The road to three-dimensional imaging in orthodontics
JPH119583A (en) X-ray ct scanner
US4769756A (en) Systematic method for matching existing radiographic projections with radiographs to be produced from a specified region of interest in cancellous bone
JP2007289569A (en) Medical image processor and medical image processing method
KR20180004134A (en) Method for improving image data from a tooth image generation system
Li et al. X-ray digital intra-oral tomosynthesis for quasi-three-dimensional imaging: system, reconstruction algorithm, and experiments
Montúfar et al. Perspective and orthogonal CBCT/CT digitally reconstructed radiographs compared to conventional cephalograms
JPH08126634A (en) Image display of three-dimensional (stereoscopic) or four-dimensional (three dimensions having time axis) data in stereoscopic imaging centering around blood vessel angiography utilizing x-rays

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20060721

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20060721

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090217

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090420

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20091104

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20091228

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: 20100122

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20100203

R150 Certificate of patent or registration of utility model

Ref document number: 4458309

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130219

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140219

Year of fee payment: 4

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

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

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term