JP2005160544A - 画像処理装置 - Google Patents
画像処理装置 Download PDFInfo
- Publication number
- JP2005160544A JP2005160544A JP2003400083A JP2003400083A JP2005160544A JP 2005160544 A JP2005160544 A JP 2005160544A JP 2003400083 A JP2003400083 A JP 2003400083A JP 2003400083 A JP2003400083 A JP 2003400083A JP 2005160544 A JP2005160544 A JP 2005160544A
- Authority
- JP
- Japan
- Prior art keywords
- image
- ray
- filter
- edge
- edge intensity
- 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.)
- Pending
Links
- 238000009499 grossing Methods 0.000 claims abstract description 22
- 238000001914 filtration Methods 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims description 31
- 238000012545 processing Methods 0.000 claims description 30
- 230000000694 effects Effects 0.000 claims description 13
- 238000002591 computed tomography Methods 0.000 claims description 10
- 230000009467 reduction Effects 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 22
- 230000008569 process Effects 0.000 description 11
- 238000003384 imaging method Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 8
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000007781 pre-processing Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 239000000470 constituent Substances 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 210000001015 abdomen Anatomy 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 239000006096 absorbing agent Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Transforming Light Signals Into Electric Signals (AREA)
Abstract
【課題】本発明の目的は、エッジ強度に加えて、スキャン条件や画像再構成条件に応じた特性でフィルタをかけることでX線CT画像の画質を向上することにある。
【解決手段】X線CT画像の局所領域におけるエッジ強度を算出するエッジ強度算出部120と、エッジ強度に基づいて、局所領域ごとに平滑化程度に関するフィルタ特性を決定するフィルタ要素算出部123と、局所領域ごとに決定したフィルタ特性に従ってX線CT画像にフィルタをかけるガウシアンフィルタ処理部122と、X線CT画像に対応する画像再構成条件とエッジ強度の算出方法との対応を記憶する記憶部118とを具備する。
【選択図】 図1
【解決手段】X線CT画像の局所領域におけるエッジ強度を算出するエッジ強度算出部120と、エッジ強度に基づいて、局所領域ごとに平滑化程度に関するフィルタ特性を決定するフィルタ要素算出部123と、局所領域ごとに決定したフィルタ特性に従ってX線CT画像にフィルタをかけるガウシアンフィルタ処理部122と、X線CT画像に対応する画像再構成条件とエッジ強度の算出方法との対応を記憶する記憶部118とを具備する。
【選択図】 図1
Description
本発明は、X線CT画像の画質を改善する画像処理装置およびX線CT装置に関する。
画像のエッジを保持、あるいは鮮鋭化しながら、ノイズを抑制することで画質を改善する手法として、画像の局所的なエッジ信号を抽出し、エッジ信号の強度が強いほど鮮鋭化効果が、エッジ信号の強度が弱いほど平滑化効果が高くなるように局所的にフィルタの特性を変更して処理する方法が一般的に知られている。空間フィルタとして、2次微分フィルタの演算結果を用いてフィルタ特性を変える方法と、ウェーブレットフィルタとして、周波数空間で非線形フィルタ処理を行う方法がある。
X線CT画像は、管電圧やスキャン速度等のスキャン条件、さらに再構成関数、拡大率(再構成FOV)等の画像再構成条件により画質(粒状性・ノイズ・コントラスト・鮮鋭度・アーチファクトなど)が大きく変動するため、上述した従来の技術を用いても予定した画質向上効果が得られないことがある。
本発明の目的は、エッジ強度に加えて、スキャン条件や画像再構成条件に応じた特性でフィルタをかけることでX線CT画像の画質を向上することにある。
本発明の第1局面は、X線CT画像の局所領域におけるエッジ強度を算出するエッジ強度算出手段と、前記エッジ強度に基づいて、局所領域ごとに平滑化程度に関するフィルタ特性を決定するフィルタ特性決定手段と、前記局所領域ごとに決定したフィルタ特性に従って、前記X線CT画像にフィルタをかける手段と、前記X線CT画像に対応する画像再構成条件に従って前記エッジ強度の算出方法を決定するエッジ強度算出方法決定手段とを具備する。
本発明の第2局面は、X線CT画像の局所領域におけるエッジ強度を算出するエッジ強度算出手段と、前記エッジ強度と、前記X線CT画像に対応する画像再構成条件とに基づいて、局所領域ごとに平滑化程度に関するフィルタ特性を決定するフィルタ特性決定手段と、前記決定したフィルタ特性に従って前記X線CT画像にフィルタをかけるフィルタ手段とを具備する。
本発明の第2局面は、X線CT画像の局所領域におけるエッジ強度を算出するエッジ強度算出手段と、前記エッジ強度と、前記X線CT画像に対応する画像再構成条件とに基づいて、局所領域ごとに平滑化程度に関するフィルタ特性を決定するフィルタ特性決定手段と、前記決定したフィルタ特性に従って前記X線CT画像にフィルタをかけるフィルタ手段とを具備する。
本発明によれば、エッジ強度に加えて、スキャン条件や画像再構成条件に応じた特性でフィルタをかけることによりX線CT画像の画質を向上することができる。
以下、図面を参照して本発明による画像処理装置の実施形態を説明する。本発明による画像処理装置は単独で又はX線コンピュータ断層撮影装置若しくはPACSに組み込まれて構成される。ここでは本発明による画像処理装置はX線コンピュータ断層撮影装置に組み込まれたものとして説明する。周知のとおり、X線コンピュータ断層撮影装置には、X線管と放射線検出器とが1体として被検体の周囲を回転する回転/回転型と、リング状に多数の検出素子がアレイされ、X線管のみが被検体の周囲を回転する固定/回転型等様々なタイプがあり、いずれの型でも適用可能である。ここでは、現在、主流を占めている回転/回転型として説明する。また、1スライスの断層画像データを再構成するには、被検体の周囲1周、約360°分の投影データが、またハーフスキャン法でも180°+ビュー角分の投影データが必要とされる。いずれの再構成方式にも適用可能である。ここでは、前者を例に説明する。また、入射X線を電荷に変換するメカニズムは、シンチレータ等の蛍光体でX線を光に変換し更にその光をフォトダイオード等の光電変換素子で電荷に変換する間接変換形と、X線による半導体内の電子正孔対の生成及びその電極への移動すなわち光導電現象を利用した直接変換形とが主流である。X線検出素子としては、それらのいずれの方式を採用してもよいが、ここでは、前者の間接変換形として説明する。また、近年では、X線管とX線検出器との複数のペアを回転フレームに搭載したいわゆる多管球型のX線CT装置の製品化が進み、その周辺技術の開発が進んでいる。本発明では、従来からの一管球型のX線CT装置であっても、多管球型のX線CT装置であってもいずれにも適用可能である。ここでは、一管球型として説明する。
図1に、本実施形態に係るX線コンピュータ断層撮影装置の構成をブロック図により示している。架台部100は、架台駆動部107により回転駆動される円環形状の回転フレーム102を有する。回転フレーム102には、X線管101とX線検出器103とが搭載される。X線検出器103は、撮影時に被検体が挿入される空洞(撮影領域)Sを挟んでX線管101に対向する。X線管101には、高電圧発生装置109から連続的又は周期的に発生される高電圧がスリップリング108を経由して印加される。それによりコーンビーム状のX線が発生される。X線検出器103は、複数の検出素子を有する。複数の検出素子は、チャンネル方向とスライス方向との直交2方向に沿って配列される。このようなX線検出器103は、一般的に、マルチスライス形(2次元アレイ形ともいう)と呼ばれる。
X線検出器103の出力には、一般的にDAS(data acquisition system) と呼ばれているデータ収集回路104が接続されている。このデータ収集回路104には、X線検出器103の各チャンネルの電流信号を電圧に変換するI−V変換器と、この電圧信号をX線の曝射周期に同期して周期的に積分する積分器と、この積分器の出力信号を増幅するアンプと、このプリアンプの出力信号をディジタル信号に変換するアナログ・ディジタル・コンバータとが、チャンネルごとに設けられている。前処理装置106は、データ収集回路104の出力(ディジタル信号(純生データと呼ばれる))を光又は磁気を媒介して非接触のデータ伝送を実現する非接触データ伝送装置105を経由して受け取り、生データに対して、チャンネル間の感度不均一を補正したり、またX線強吸収体、主に金属部による極端な信号強度の低下又は信号脱落を補正する等の前処理を実行する。前処理装置106から出力されるデータは、投影データと呼ばれる。投影データは補助記憶装置112に記憶され、ホストコントローラ110の制御に従って再構成装置114に読み出される。再構成装置114は、記憶装置112からの投影データのセットに基づいてフェルドカンプ法又は他の再構成法に従ってスライス又はボリュームに関する画像データを再構成する。この画像データは記憶装置112に記憶され、また表示装置116に可視画像として表示される。
ホストコントローラ110には、X線CT装置としての基本的な構成要素としての上記記憶装置112、再構成装置114、表示装置116、キーボードやポインティングデバイス等を有する入力装置115の他に、フィルタリング制御部117、再構成FOV/サブ領域サイズ対応表記憶部118、再構成FOV/フィルタサイズ対応表記憶部119、エッジ強度算出部120、エッジ強度/フィルタパラメータ対応表記憶部121、ガウシアンフィルタ処理部122、フィルタ要素算出部123が接続される。フィルタリング制御部117は、再構成FOV/サブ領域サイズ対応表記憶部118、再構成FOV/フィルタサイズ対応表記憶部119、エッジ強度算出部120、エッジ強度/フィルタパラメータ対応表記憶部121、ガウシアンフィルタ処理部122、フィルタ要素算出部123とともに画像処理装置を構成する。
これら構成要素の機能を本実施形態による画像処理動作と共に説明する。本動作前に、被検体の任意の部位を含む範囲が、予め設定されたスキャン条件のもとでスキャンされる。スキャン条件には、周知の通り、管電圧、管電流、撮影FOV(撮影領域径)、撮影スライス厚、撮影スライスピッチ、ヘリカルピッチ等が含まれる。当該スキャンにより収集されたデータに基づいて、予め設定された再構成条件のもとで再構成装置114で画像(X線CT画像)データが再構成され、記憶装置112に記憶される。再構成条件には、再構成再構成関数、再構成FOV(再構成領域径)、画像スライス厚、画像スライスピッチ等が含まれる。撮影FOV内の指定した一部分だけを装置固有の最高空間解像度で再構成するいわゆる拡大再構成(ズーミング再構成)では、指定した再構成FOVに応じて画像拡大率が決まる。例えば撮影FOVが512mmと仮定すると、256mmの再構成FOVでは、画像拡大率は2倍になる。
画像処理装置は、再構成されたX線CT画像データに対して画像処理、特にエッジを残しながら、ノイズを抑圧するいわゆる量子フィルタ処理を施すために設けられている。この量子フィルタ処理には、平滑化フィルタ処理をベースとしたものと、鮮鋭化フィルタ処理をベースとしたものとがある。平滑化と鮮鋭化とは通常、反比例の関係にある。平滑化を強くかけることで、高いノイズ抑圧効果を発揮できるが、その反面、臓器等の輪郭を強調するエッジ鮮鋭化効果は低減する(エッジが鈍る)。逆に、弱く平滑化を掛けることで、ノイズ抑圧効果はあまり発揮できないが、その反面、臓器等の輪郭を強調するエッジ鮮鋭化効果を発揮できる。ここでは平滑化ベースのフィルタとしてガウシアンフィルタ処理を一例として、エッジを保持し、画像をぼかすことなく、ノイズを低減することを目指すフィルタについて示す。
量子フィルタ処理においては、先ず、図2に示すように、エッジ強度算出部120において入力画像Iin(i,j)の各画素(i,j)におけるエッジ強度IE(i,j)が算出される。比較的高いエッジ強度IE(i,j)はその画素(i,j)がエッジ(輪郭)上の点である可能性が高いことを表しており、比較的低いエッジ強度IE(i,j)はその注目画素(i,j)がエッジ(輪郭)上の点である可能性が低いことを表している。エッジ強度算出手段として、2次微分フィルタが使用される。一般的に2次微分フィルタは、
で表され、2次微分フィルタh1と入力画像Iin(i,j)を積和演算することで、画素(i,j)における画素値の変化量、すなわちエッジ強度を算出することができる。
となる。
ところが、2次微分フィルタh1は、画素ごとの画素値変化量からエッジ強度を算出するため、再構成条件としての例えば画像拡大率に応じてエッジ強度の算出方法を変更する必要がある。ここで、図3に示すとおり、フィルタ処理を実施する対象となる注目画素(i,j)を含むn×m×l画素のサブ領域0と、サブ領域0の周囲にn×m×l画素のサブ領域を8個設定し、各サブ領域のCT値の平均に対し、2次微分フィルタをかける。すなわち、サブ領域kの平均値をave(k)とした場合、エッジ強度は、
となる。画像の拡大率に比例し、サブ領域のサイズn×m×lを変更することで、画像の拡大率に依存しないで、安定したノイズ除去効果を実現することが可能となる。例えば、図4に示すように、画像の拡大率に依存しないで、安定したノイズ除去効果を実現するように再構成FOV(拡大率)とサブ領域のサイズとの適当な対応を決定し、予め再構成FOV/サブ領域サイズ対応表記憶部118に記憶させておく。エッジ強度算出処理にあたっては、フィルタリング制御部117の制御のもとで、予め再構成条件として設定されている再構成FOVに対応するサブ領域サイズがその記憶部118から選択的にエッジ強度算出部120に選択的に読み出される。図4の例では、再構成FOV(画像拡大率)が256mm(2倍)までがサブ領域サイズが3×3に設定され、再構成FOV(画像拡大率)が256mm(2倍)を超えたときサブ領域サイズが5×5に設定されるように2段階で示しているが、再構成FOV(画像拡大率)が高くなるに応じてサブ領域サイズが大きくなる対応関係である限りにおいて、再構成FOV(画像拡大率)に対するサブ領域サイズの対応は図4の例に限定されるものではない。なお、中心のサブ領域に対して周囲のサブ領域のサイズを違えるようにしてもよいし、周囲のサブ領域どうしでも違えるようにしても良い。また、使用する平均値ave(k)は、単純平均に限定されず、重み付き平均でも構わない。
なお、上述では、画像の拡大率に応じてサブ領域のサイズを決定したが、それに代えて又はそれとともに、画像の拡大率に応じて2次微分フィルタのフィルタサイズN×Mを変更するようにしても良い。例えば、図5に示すように、再構成FOV(拡大率)と2次微分フィルタのフィルタサイズN×Mとの適当な対応を予め決定し、再構成FOV/フィルタサイズ対応表記憶部119に記憶させておく。エッジ強度算出処理にあたっては、フィルタリング制御部117の制御のもとで、予め再構成条件として設定されている再構成FOVに対応するフィルタサイズがその記憶部119から選択的にエッジ強度算出部120に選択的に読み出される。図5の例では、再構成FOV(画像拡大率)が256mm(2倍)までがフィルタサイズが3×3に設定され、256mm(2倍)から128mm(4倍)までがフィルタサイズが5×5に設定され、128mm(4倍)を超過したときフィルタサイズは7×7に設定されるように3段階で示しているが、再構成FOV(画像拡大率)が高くなるに応じてフィルタサイズが大きくなる対応関係である限りにおいて、再構成FOV(画像拡大率)に対するフィルタサイズの対応は図5の例に限定されるものではない。
として表され、フィルタパラメータσによって、平滑化の度合いを制御することができる。フィルタパラメータσは0より大きく、1以下の範囲から選択され、0に接近するほど、平滑化の度合いは低くなる。よって、図6に示すように、エッジ強度IEが高くなるほどフィルタパラメータσが0に接近して平滑化の度合いが小さくなり、逆に鮮鋭化の度合いが大きくなるような特性をもつ関数σ=f(IE)を設定することで、エッジ部分の信号を鈍らせること無く、平坦部のノイズを平滑化により低減することができる。図6に示すように、関数σ=f(IE)を再構成FOV(画像の拡大率)の相違に応じてここでは3種類用意して、記憶部121に記憶する。フィルタ要素の算出処理にあたっては、フィルタリング制御部117の制御のもとで、予め再構成条件として設定されている再構成FOVと、計算したエッジ強度とに対応するフィルタパラメータσが画素ごとに記憶部121から選択的にフィルタ要素計算部123に選択的に読み出される。フィルタ要素計算部123では、フィルタパラメータσに基づいて式(3)の各フィルタ要素が計算される。計算されたフィルタ要素からなるガウシアンフィルタによりX線CT画像がフィルタ処理に供される。なお、上述では、フィルタパラメータσが再構成FOVと計算したエッジ強度とに対して対応付けられて記憶されているが、再構成FOVは頭部、胸部、腹部等の部位に応じてほぼ決まることから、フィルタパラメータσが頭部、胸部、腹部等の部位と計算したエッジ強度とに対して対応付けられて記憶されるようにしてもよい。
上述では、再構成条件の画像拡大率に対してサブ領域のサイズ、フィルタサイズ、フィルタパラメータを対応付けるようにしたが、再構成条件の画像拡大率に代えて、スキャン条件の管電圧、スキャン条件の管電流、スキャン条件の撮影スライス厚、スキャン条件の撮影スライスピッチ、スキャン条件のヘリカルピッチ、再構成条件の再構成関数、再構成条件の画像スライス厚、または再構成条件の画像スライスピッチに対してサブ領域のサイズ、フィルタサイズ、フィルタパラメータを対応付けるようにしてもよい。
なお、撮像時の管電圧や撮像部位等により、画像のノイズ強度が変わる。それにより平坦部分でもエッジ強度が高くなるため、ノイズ強度が大きくなっているにも関わらず、ノイズ低減効果が低くなるため画質が改善されにくくなるという問題が生じる。そこで、2次微分フィルタでエッジ強度を計算する前に、ノイズ強度に比例した平滑化効果を有するフィルタにより、ノイズを低減した後、2次微分フィルタを実行することで、ノイズによる影響を軽減することができる。ノイズ強度の計算方法は従来からあるノイズ強度計算方法から任意の方法を選択的に採用すれば良い。上記は、平滑化フィルタと2次微分フィルタを順番に処理してもよいし、平滑化フィルタと2次微分フィルタとの効果を併せ持つ1個のフィルタにより1回で処理することも可能である。
また、エッジ強度の大きい画像では、ガウシアンフィルタの平滑化効果を高めなくては、ノイズが低減されず画質が改善されない。よって、エッジIEからガウシアンフィルタのパラメータσを算出する関数σ=f(IE)は、ノイズ強度cも考慮した関数σ=f’(IE,c)とし、ノイズ強度cが大きくなるほどσが大きくなる特性を与えるようにしても良い。つまり、フィルタパラメータσがエッジ強度とノイズ強度とに対応付けられて記憶され、フィルタ要素の算出に際しては、エッジ強度とノイズ強度とに対応するフィルタパラメータσが記憶部121から選択的にフィルタ要素計算部123に選択的に読み出される。
以上、画像を平滑化するガウシアンフィルタへの適用例を述べたが、画像の先鋭化するウェーブレットフィルタであるMRAフィルタにも適用可能である。図7に示すとおり、MRAでは、画像を低周波成分と高周波成分とに分解し、その低周波成分をさらに低周波成分と高周波成分とに分解し、その分解を階層的にnレベルまで繰り返し、各レベルの高周波成分に対してそのレベルに応じてゲインをかけ、各レベルの低周波成分はそのまま維持して、合成処理することによりエッジを強調する鮮鋭化処理を行う。
画像のエッジを保持しながら、ノイズを除去するフィルタを構成する場合、図8(a)、図8(b)に示すとおり、低周波成分はそのまま通過させ、所定の分解レベルL未満では、強度があるしきい値th未満の高周波成分は除去し、しきい値th以上の高周波成分はゲインを与えてエッジとして通過するような非線形処理を行う。スキャン条件又は画像の拡大率等の再構成条件に応じて分解レベルLを変更し、またノイズ強度に応じてしきい値thを変更することで、スキャン条件や再構成条件が変化した場合でも、一定のフィルタ効果を得ることができる。
以上、2次元画像についての説明したが、3次元以上の多次元画像についても同様に適用可能である。
(変形例)
本発明は、上述した実施形態に限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で種々変形して実施することが可能である。さらに、上記実施形態には種々の段階が含まれており、開示される複数の構成要件における適宜な組み合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件から幾つかの構成要件が削除されてもよい。
本発明は、上述した実施形態に限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で種々変形して実施することが可能である。さらに、上記実施形態には種々の段階が含まれており、開示される複数の構成要件における適宜な組み合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件から幾つかの構成要件が削除されてもよい。
101…X線管、102…回転フレーム、103…2次元検出器、104…DAS、105…非接触データ電装装置、106…前処理装置、107…架台駆動部、108…スリップリング、109…高電圧発生装置、110…ホストコントローラ、112…記憶装置、114…再構成装置、115…入力装置、116…表示装置、117…フィルタリング制御部、118…再構成FOV/サブ領域サイズ対応表記憶部、119…再構成FOV/フィルタサイズ対応表記憶部、120…エッジ強度算出部、121…エッジ強度/フィルタパラメータ対応表記憶部、122…ガウシアンフィルタ処理部、123…フィルタ要素算出部。
Claims (12)
- X線CT画像の局所領域におけるエッジ強度を算出するエッジ強度算出手段と、
前記エッジ強度に基づいて、局所領域ごとに平滑化程度に関するフィルタ特性を決定するフィルタ特性決定手段と、
前記局所領域ごとに決定したフィルタ特性に従って、前記X線CT画像にフィルタをかける手段と、
前記X線CT画像に対応する画像再構成条件に従って前記エッジ強度の算出方法を決定するエッジ強度算出方法決定手段とを具備することを特徴とする画像処理装置。 - 前記エッジ強度算出方法決定手段は、n×m×lサイズを有するサブ領域の平均値に対して2次微分フィルタをかけることにより前記エッジ強度を算出するとともに、前記サブ領域のn×m×lサイズを前記画像再構成条件に含まれる画像拡大率に応じて変更することを特徴とする請求項1記載の画像処理装置。
- 前記エッジ強度算出方法決定手段は、前記画像拡大率に対する前記サブ領域のサイズの対応関係を記憶する記憶する手段を有することを特徴とする請求項2記載の画像処理装置。
- 前記エッジ強度算出方法決定手段は、前記X線CT画像に対して2次微分フィルタをかけることにより前記エッジ強度を算出するとともに、前記2次微分フィルタのサイズを前記画像再構成条件に含まれる画像拡大率に応じて変更することを特徴とする請求項1記載の画像処理装置。
- 前記エッジ強度の算出前に前記X線CT画像に対してノイズ低減処理を施す手段を更に備えることを特徴とする請求項1記載の画像処理装置。
- 前記エッジ強度算出手段は、エッジ強度算出とともにノイズ低減効果を有する2次微分フィルタを前記X線CT画像に対して施すことを特徴とする請求項1記載の画像処理装置。
- X線CT画像の局所領域におけるエッジ強度を算出するエッジ強度算出手段と、
前記エッジ強度と、前記X線CT画像に対応する画像再構成条件とに基づいて、局所領域ごとに平滑化程度に関するフィルタ特性を決定するフィルタ特性決定手段と、
前記決定したフィルタ特性に従って前記X線CT画像にフィルタをかけるフィルタ手段とを具備することを特徴とする画像処理装置。 - 前記フィルタ特性決定手段は、前記エッジ強度と、前記画像再構成条件に含まれる画像拡大率とに基づいて前記フィルタ特性を決定することを特徴とする請求項7記載の画像処理装置。
- 前記フィルタ特性決定手段は、前記エッジ強度と、前記画像再構成条件に含まれる画像拡大率と、前記X線CT画像のノイズ強度とに基づいて前記フィルタ特性を決定することを特徴とする請求項7記載の画像処理装置。
- 前記フィルタ特性決定手段は、前記エッジ強度と、前記画像再構成条件に含まれる画像拡大率と、前記X線CT画像に撮像された部位とに基づいて前記フィルタ特性を決定することを特徴とする請求項7記載の画像処理装置。
- 被検体をスキャンして収集したデータに基づいて前記被検体に関するX線CT画像を再構成するX線コンピュータ断層撮影装置において、
前記X線CT画像の局所領域におけるエッジ強度を算出するエッジ強度算出手段と、
前記エッジ強度に基づいて、局所領域ごとに平滑化程度に関するフィルタ特性を決定するフィルタ特性決定手段と、
前記局所領域ごとに決定したフィルタ特性に従って、前記X線CT画像にフィルタをかける手段と、
前記X線CT画像に対応する画像再構成条件に従って前記エッジ強度の算出方法を決定するエッジ強度算出方法決定手段とを具備することを特徴とするX線コンピュータ断層撮影装置。 - 被検体をスキャンして収集したデータに基づいて前記被検体に関するX線CT画像を再構成するX線コンピュータ断層撮影装置において、
X線CT画像の局所領域におけるエッジ強度を算出するエッジ強度算出手段と、
前記エッジ強度と、前記X線CT画像に対応する画像再構成条件とに基づいて、局所領域ごとに平滑化程度に関するフィルタ特性を決定するフィルタ特性決定手段と、
前記決定したフィルタ特性に従って前記X線CT画像にフィルタをかけるフィルタ手段とを具備することを特徴とするX線コンピュータ断層撮影装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003400083A JP2005160544A (ja) | 2003-11-28 | 2003-11-28 | 画像処理装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003400083A JP2005160544A (ja) | 2003-11-28 | 2003-11-28 | 画像処理装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2005160544A true JP2005160544A (ja) | 2005-06-23 |
Family
ID=34724445
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003400083A Pending JP2005160544A (ja) | 2003-11-28 | 2003-11-28 | 画像処理装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2005160544A (ja) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011158893A1 (ja) * | 2010-06-17 | 2011-12-22 | 株式会社 日立メディコ | X線ct装置、及びその制御方法 |
JP2012175992A (ja) * | 2011-02-25 | 2012-09-13 | Tanita Corp | 体組成計 |
US8374460B2 (en) | 2008-07-29 | 2013-02-12 | Ricoh Company, Ltd. | Image processing unit, noise reduction method, program and storage medium |
JP2014236810A (ja) * | 2013-06-06 | 2014-12-18 | キヤノン株式会社 | 画像処理装置、断層撮影装置、画像処理方法およびプログラム |
JP2016526962A (ja) * | 2013-06-10 | 2016-09-08 | ユニバーシティー オブ ミシシッピ メディカル センター | 医用画像処理法 |
JP2016202246A (ja) * | 2015-04-15 | 2016-12-08 | 東芝メディカルシステムズ株式会社 | X線診断装置 |
JP2019213960A (ja) * | 2019-09-25 | 2019-12-19 | キヤノンメディカルシステムズ株式会社 | X線診断装置 |
CN111795980A (zh) * | 2020-08-04 | 2020-10-20 | 合肥工业大学 | 一种基于逐像素高斯函数拟合法的x射线边界照明成像方法 |
JP2021021600A (ja) * | 2019-07-25 | 2021-02-18 | 株式会社島津製作所 | X線イメージング装置およびx線イメージング方法 |
-
2003
- 2003-11-28 JP JP2003400083A patent/JP2005160544A/ja active Pending
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8374460B2 (en) | 2008-07-29 | 2013-02-12 | Ricoh Company, Ltd. | Image processing unit, noise reduction method, program and storage medium |
US9025848B2 (en) | 2010-06-17 | 2015-05-05 | Hitachi Medical Corporation | X-ray CT apparatus and control method for reconstructing X-ray CT images from filtered projection data |
CN102946807A (zh) * | 2010-06-17 | 2013-02-27 | 株式会社日立医疗器械 | X射线ct装置及其控制方法 |
WO2011158893A1 (ja) * | 2010-06-17 | 2011-12-22 | 株式会社 日立メディコ | X線ct装置、及びその制御方法 |
JP5878119B2 (ja) * | 2010-06-17 | 2016-03-08 | 株式会社日立メディコ | X線ct装置、及びその制御方法 |
JP2012175992A (ja) * | 2011-02-25 | 2012-09-13 | Tanita Corp | 体組成計 |
US9418417B2 (en) | 2013-06-06 | 2016-08-16 | Canon Kabushiki Kaisha | Image processing apparatus, tomography apparatus, image processing method, and storage medium |
EP2824638A1 (en) | 2013-06-06 | 2015-01-14 | Canon Kabushiki Kaisha | Image processing apparatus and image processing method |
JP2014236810A (ja) * | 2013-06-06 | 2014-12-18 | キヤノン株式会社 | 画像処理装置、断層撮影装置、画像処理方法およびプログラム |
JP2016526962A (ja) * | 2013-06-10 | 2016-09-08 | ユニバーシティー オブ ミシシッピ メディカル センター | 医用画像処理法 |
JP2019088810A (ja) * | 2013-06-10 | 2019-06-13 | ユニバーシティー オブ ミシシッピ メディカル センター | 医用画像処理法 |
JP2016202246A (ja) * | 2015-04-15 | 2016-12-08 | 東芝メディカルシステムズ株式会社 | X線診断装置 |
JP2021021600A (ja) * | 2019-07-25 | 2021-02-18 | 株式会社島津製作所 | X線イメージング装置およびx線イメージング方法 |
JP7180566B2 (ja) | 2019-07-25 | 2022-11-30 | 株式会社島津製作所 | X線イメージング装置およびx線イメージング方法 |
JP2019213960A (ja) * | 2019-09-25 | 2019-12-19 | キヤノンメディカルシステムズ株式会社 | X線診断装置 |
CN111795980A (zh) * | 2020-08-04 | 2020-10-20 | 合肥工业大学 | 一种基于逐像素高斯函数拟合法的x射线边界照明成像方法 |
CN111795980B (zh) * | 2020-08-04 | 2022-04-26 | 合肥工业大学 | 一种基于逐像素高斯函数拟合法的x射线边界照明成像方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7865006B2 (en) | Methods and apparatus for artifact reduction | |
JP5180181B2 (ja) | コンピュータ断層撮影データ収集装置及び方法 | |
US8965078B2 (en) | Projection-space denoising with bilateral filtering in computed tomography | |
JP4820582B2 (ja) | ヘリカルマルチスライスctのための回復ノイズを伴うヘリカルウィンドミルアーチファクトを低減する方法 | |
JP4020451B2 (ja) | もとの像データを強化するシステム | |
US10789738B2 (en) | Method and apparatus to reduce artifacts in a computed-tomography (CT) image by iterative reconstruction (IR) using a cost function with a de-emphasis operator | |
WO2010028027A1 (en) | Method for reconstruction in dual energy, dual source helical computed tomography | |
JP2009131464A (ja) | X線ct装置 | |
US10702234B2 (en) | Image combining using images with different focal-spot sizes | |
US8284893B2 (en) | X-ray computer tomography apparatus and image processing apparatus | |
JP2007307417A (ja) | 像データの処理方法及び像データの処理装置 | |
JP4519434B2 (ja) | 超解像処理装置及び医用画像診断装置 | |
US8731267B2 (en) | X-ray computed tomography apparatus and image processing apparatus | |
KR101076321B1 (ko) | 콘빔 ct 장치에서의 3차원 영상 획득 방법 및 이를 적용한 콘빔 ct 장치 | |
US5963614A (en) | Filters for single slice helical image reconstruction in a computed tomography system | |
JP2004000632A (ja) | 物体の画像を再構成する方法及び装置 | |
CN111065335A (zh) | 医用图像处理装置和医用图像处理方法 | |
JP2005160544A (ja) | 画像処理装置 | |
JP5329204B2 (ja) | X線ct装置 | |
WO2016132880A1 (ja) | 演算装置、x線ct装置、及び画像再構成方法 | |
US8385620B2 (en) | Method and system for multi-detector-row computed tomography utilizing projection data up-sampling with shift | |
US6647084B1 (en) | Method and apparatus for filtering projection data of a helical scan | |
JP6615531B2 (ja) | X線コンピュータ断層撮影装置及び医用画像処理装置 | |
US6091840A (en) | Methods and apparatus for single slice helical image reconstruction in a computed tomography system | |
JP2005095329A (ja) | 超解像処理装置及び医用画像診断装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20061121 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20091020 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20100302 |