JP4535795B2 - Image processing apparatus and X-ray CT system - Google Patents
Image processing apparatus and X-ray CT system Download PDFInfo
- Publication number
- JP4535795B2 JP4535795B2 JP2004205079A JP2004205079A JP4535795B2 JP 4535795 B2 JP4535795 B2 JP 4535795B2 JP 2004205079 A JP2004205079 A JP 2004205079A JP 2004205079 A JP2004205079 A JP 2004205079A JP 4535795 B2 JP4535795 B2 JP 4535795B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- comparison
- image quality
- processing apparatus
- reconstruction
- 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
Links
- 238000012545 processing Methods 0.000 title claims description 39
- 238000000034 method Methods 0.000 claims description 90
- 230000006872 improvement Effects 0.000 claims description 41
- 238000005259 measurement Methods 0.000 claims description 17
- 238000007781 pre-processing Methods 0.000 claims description 10
- 230000000052 comparative effect Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 description 13
- 230000008569 process Effects 0.000 description 12
- 238000012805 post-processing Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 8
- 238000001514 detection method Methods 0.000 description 7
- 238000003672 processing method Methods 0.000 description 5
- 238000013480 data collection Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 210000004185 liver Anatomy 0.000 description 2
- 230000032258 transport Effects 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
本発明は、スキャンされた被検者のX線投影データからX線断層像を再構成する画像処理装置及び画像処理方法並びにX線CTシステムに関する。 The present invention relates to an image processing apparatus, an image processing method, and an X-ray CT system for reconstructing an X-ray tomographic image from scanned X-ray projection data of a subject.
X線CT(Computerized Tomography)システムは、ドーナツ状の空洞部を有してX線を被検者(被検体)に照射して投影データを計測する計測装置(以下、「ガントリ」と呼ぶ。)と、ガントリに対して各種制御信号を与えるとともに、ガントリにおいて計測されたX線の投影データに基づいてX線断層像を再構成して画面表示等する操作コンソールと、被検者をガントリの空洞部内に固定・支持するために、被検者を載せたスライド可能なテーブル部を当該空洞部に向けて搬送するためのテーブル装置とで構成される。 An X-ray CT (Computerized Tomography) system has a doughnut-shaped cavity, and measures a projection data by irradiating a subject (subject) with X-rays (hereinafter referred to as “gantry”). In addition to providing various control signals to the gantry, an operation console for reconstructing an X-ray tomogram based on the X-ray projection data measured in the gantry and displaying it on the screen, etc. In order to fix and support the inside of the unit, the table unit is configured to transport a slidable table unit on which a subject is placed toward the cavity.
このように、X線CTシステムは、被検者を透過した投影データを収集して、当該データから被検者の特定の部位のX線断層像を再構成することを主目的とする。一般に、X線CTシステムでは、異なる角度から被検者に向けてX線を照射し、各角度で被検者を透過したX線の投影データを検出する。一般に、上述した一連のX線の照射・検出動作を「スキャン」と呼ぶ。そして、ガントリで検出された投影データは操作コンソールで受信され、前処理された後にX線断層像が再構成される。通常、この場合の再構成方法として、フィルタ補正逆投影(FBP:Filter Back Projection)法が用いられている。また、再構成画像に対して逐次近似法を用いて後処理を行うことにより、所定部位の再構成画像の画質を改善することが知られている(例えば、特許文献1参照)。 Thus, the main purpose of the X-ray CT system is to collect projection data transmitted through the subject and reconstruct an X-ray tomographic image of a specific part of the subject from the data. In general, in an X-ray CT system, X-rays are irradiated toward a subject from different angles, and X-ray projection data transmitted through the subject at each angle is detected. In general, the series of X-ray irradiation / detection operations described above is called “scanning”. Then, the projection data detected by the gantry is received by the operation console, and after pre-processing, the X-ray tomogram is reconstructed. Normally, a filter back projection (FBP) method is used as a reconstruction method in this case. It is also known to improve the image quality of a reconstructed image at a predetermined site by performing post-processing on the reconstructed image using a successive approximation method (see, for example, Patent Document 1).
また、逐次近似法を用いた画像再構成法として、期待値最大化(EM:Expectation Maximization)アルゴリズムを用いた最尤推定(ML:Maximum Likelihood)法(MLEM法)が知られており、これを改善したものとしてサブセット化による期待値最大化(OSEM:Ordered Subset Expectaion Maximization)法も知られている。
しかしながら、上記逐次近似法において、MLEM法は計算速度が遅く、OSEM法は繰り返し演算が多いため、FBP法のような解析的方法と比べて計算時間が多くかかってしまうという問題があった。そのため、スキャンによって取得された投影データはFBP法を用いて再構成された後、特に詳細に再構成する必要があるような部位については、さらにOSEM法を用いた後処理が行われていた。しかしながら、OSEM法を用いた従来の後処理の演算においては、OSEM法を用いてどの程度演算すればよいかといった演算の終了条件が定まっていない。そのため、従来は検査者等の経験によって決められているにすぎず、後処理をした後の再構成画像の定量的な評価をすることができなかった。 However, in the above-described successive approximation method, the MLEM method has a slow calculation speed, and the OSEM method has many repetitive operations, so that there is a problem that it takes a lot of calculation time compared to an analytical method such as the FBP method. Therefore, after the projection data acquired by the scan is reconstructed using the FBP method, a post-process using the OSEM method is further performed on a part that needs to be reconstructed in detail. However, in the conventional post-processing calculation using the OSEM method, the calculation end condition such as how much calculation should be performed using the OSEM method is not determined. For this reason, conventionally, it is determined only by the experience of the examiner and the like, and the reconstructed image after the post-processing cannot be quantitatively evaluated.
本発明は、このような事情を考慮してなされたものであり、CT画像の特定部位について、エッジ部分のシャープさを保ったまま、それ以外の部分をより滑らかに補正することができるとともに、再構成画像の評価を定量的に行うことができる画像処理装置及び画像処理方法並びにX線CTシステムを提供することを目的とする。 The present invention has been made in consideration of such circumstances, and for the specific part of the CT image, while maintaining the sharpness of the edge part, other parts can be corrected more smoothly, An object of the present invention is to provide an image processing apparatus, an image processing method, and an X-ray CT system capable of quantitatively evaluating a reconstructed image.
上記課題を解決するため、本発明は、被検体のX線投影データからCT画像を画像再構成する画像処理装置であって、
前記X線投影データの前処理を行う前処理手段と、
前処理された前記X線投影データからCT画像の画像再構成を行う再構成手段と、
前記再構成手段によって画像再構成された前記CT画像上に所定のROIを設定して、該ROI内の画像計測値を算出し、その画像計測値に基づいて前記CT画像の画質改善が必要か否かを判定する判定手段と、
前記判定手段によって前記CT画像再構成のパラメータ更新が必要と判定された場合、再構成パラメータを変更して画像再構成することにより、前記CT画像を画質改善する画質改善手段と、
前記判定手段によって前記CT画像の画質改善が必要でないと判定された場合、該CT画像を出力する画像出力手段と
を備えることを特徴とする。
In order to solve the above problems, the present invention is an image processing apparatus for reconstructing a CT image from X-ray projection data of a subject,
Preprocessing means for performing preprocessing of the X-ray projection data;
Reconstruction means for performing image reconstruction of a CT image from the preprocessed X-ray projection data;
Is it necessary to set a predetermined ROI on the CT image reconstructed by the reconstruction means, calculate an image measurement value in the ROI, and improve the image quality of the CT image based on the image measurement value? Determining means for determining whether or not;
When it is determined by the determination means that the CT image reconstruction parameter needs to be updated, image quality improvement means for improving the image quality of the CT image by changing the reconstruction parameter and reconstructing the image;
And an image output unit that outputs the CT image when it is determined by the determining unit that the image quality of the CT image does not need to be improved.
また、本発明は、被検体のX線投影データからCT画像を画像再構成する画像処理装置であって、
前記X線投影データからフィルタ補正逆投影法を用いてCT画像の再構成を行う再構成手段と、
前記CT画像と所定の比較画像とを用いた逐次近似法によって、該比較画像を画質改善する画質改善手段と、
前記画質改善手段によって画質改善された前記比較画像上に所定のROIを設定して、該ROI内の標準偏差を算出する算出手段と、
前記標準偏差に基づいて前記比較画像の画質改善がさらに必要か否かを判定する判定手段と、
前記判定手段によって前記比較画像の画質改善がさらに必要であると判定された場合、前記画質改善と該判定手段によって判定された前記比較画像とを用いた逐次近似法によって、該比較画像を画質再改善する画質再改善手段と、
前記判定手段によって前記比較画像の画質改善が必要でないと判定された場合、該比較画像をCT画像として出力する出力手段と
を備えることを特徴とする。
Further, the present invention is an image processing apparatus for reconstructing a CT image from X-ray projection data of a subject,
Reconstruction means for reconstructing a CT image from the X-ray projection data using a filtered back projection method;
An image quality improvement means for improving the image quality of the comparison image by a successive approximation method using the CT image and the predetermined comparison image;
Calculating means for setting a predetermined ROI on the comparative image whose image quality has been improved by the image quality improving means, and calculating a standard deviation in the ROI;
Determining means for determining whether or not further improvement in image quality of the comparison image is necessary based on the standard deviation;
If it is determined by the determining means that further improvement in image quality of the comparative image is necessary, the comparison image is re-imaged by a successive approximation method using the image quality improvement and the comparative image determined by the determining means. Image re-improvement means to improve,
And an output unit that outputs the comparison image as a CT image when it is determined by the determination unit that the image quality improvement of the comparison image is not necessary.
さらに、本発明は、被検体のX線投影データからCT画像を画像再構成する画像処理方法であって、
前記X線投影データの前処理を行う前処理工程と、
前処理された前記X線投影データからCT画像の画像再構成を行う再構成工程と、
前記再構成工程によって画像再構成された前記CT画像上に所定のROIを設定して、該ROI内の画像計測値を算出する画像計測値算出工程と、
前記画像計測値に基づいて前記CT画像の画質改善が必要か否かを判定する判定工程と、
前記判定工程によって前記CT画像の画質改善が必要と判定された場合、再構成パラメータを変更して画像再構成することにより、前記CT画像を画質改善する画質改善工程と、
前記判定工程によって前記CT画像の画質改善が必要でないと判定された場合、該CT画像を出力する画像出力工程と
を有することを特徴とする。
Furthermore, the present invention is an image processing method for reconstructing a CT image from X-ray projection data of a subject,
A preprocessing step for preprocessing the X-ray projection data;
A reconstruction step of reconstructing a CT image from the preprocessed X-ray projection data;
An image measurement value calculation step of setting a predetermined ROI on the CT image reconstructed by the reconstruction step and calculating an image measurement value in the ROI;
A determination step of determining whether image quality improvement of the CT image is necessary based on the image measurement value;
When it is determined that the image quality improvement of the CT image is necessary by the determination step, the image quality improvement step of improving the image quality of the CT image by changing the reconstruction parameter and reconstructing the image,
An image output step of outputting the CT image when it is determined that the image quality improvement of the CT image is not necessary in the determination step.
さらにまた、本発明は、被検体のX線投影データからCT画像を画像再構成する画像処理方法であって、
前記X線投影データからフィルタ補正逆投影法を用いてCT画像の画像再構成を行う再構成工程と、
前記CT画像と所定の比較画像とを用いた逐次近似法によって、該比較画像を画質改善する画質改善工程と、
前記画質改善工程によって画質改善された前記比較画像上に所定のROIを設定して、該ROI内の標準偏差を算出する算出工程と、
前記標準偏差に基づいて前記比較画像の画質改善がさらに必要か否かを判定する判定工程と、
前記判定工程によって前記比較画像の画質改善がさらに必要であると判定された場合、前記CT画像と該判定工程によって判定された前記比較画像とを用いた逐次近似法によって、該比較画像を画質再改善する画質再改善工程と、
前記判定工程によって前記比較画像の画質改善が必要でないと判定された場合、該比較画像をCT画像として出力する出力工程と
を有することを特徴とする画像処理方法。
Furthermore, the present invention is an image processing method for reconstructing a CT image from X-ray projection data of a subject,
A reconstruction step of performing image reconstruction of a CT image using a filtered back projection method from the X-ray projection data;
An image quality improvement step of improving the image quality of the comparison image by a successive approximation method using the CT image and a predetermined comparison image;
A calculation step of setting a predetermined ROI on the comparative image whose image quality has been improved by the image quality improvement step, and calculating a standard deviation in the ROI;
A determination step of determining whether or not further improvement in image quality of the comparison image is necessary based on the standard deviation;
When it is determined that the image quality of the comparison image needs to be further improved by the determination step, the comparison image is re-imaged by a successive approximation method using the CT image and the comparison image determined by the determination step. Image quality re-improvement process to improve,
An image processing method comprising: an output step of outputting the comparison image as a CT image when it is determined that the image quality improvement of the comparison image is not necessary in the determination step.
さらにまた、本発明は、被検体のX線投影データを収集するガントリと、前記X線投影データからCT画像を画像再構成する画像処理装置と、前記ガントリの空洞部に対して前記被検体を搬送するテーブル装置とを備えるX線CTシステムであって、
前記画像処理装置が、
前記X線投影データの前処理を行う前処理手段と、
前処理された前記X線投影データからCT画像の画像再構成を行う再構成手段と、
前記再構成手段によって画像再構成された前記CT画像上に所定のROIを設定して、該ROI内の画像計測値を算出する画像計測値算出手段と、
前記画像計測値に基づいて前記CT画像の画質改善が必要か否かを判定する判定手段と、
前記判定手段によって前記CT画像の画質改善が必要と判定された場合、再構成パラメータを変更して画像再構成することにより、前記CT画像を画質改善する画質改善手段と、
前記判定手段によって前記CT画像の画質改善が必要でないと判定された場合、該CT画像を出力する画像出力手段と
を備えることを特徴とする。
Furthermore, the present invention relates to a gantry that collects X-ray projection data of a subject, an image processing device that reconstructs a CT image from the X-ray projection data, and the subject to the cavity of the gantry. An X-ray CT system comprising a transporting table device,
The image processing apparatus is
Preprocessing means for performing preprocessing of the X-ray projection data;
Reconstruction means for performing image reconstruction of a CT image from the preprocessed X-ray projection data;
Image measurement value calculation means for setting a predetermined ROI on the CT image reconstructed by the reconstruction means and calculating an image measurement value in the ROI;
Determination means for determining whether or not image quality improvement of the CT image is necessary based on the image measurement value;
When it is determined by the determination means that the image quality of the CT image needs to be improved, the image quality improvement means for improving the image quality of the CT image by changing the reconstruction parameter and reconstructing the image;
And an image output means for outputting the CT image when it is determined by the determining means that no improvement in the image quality of the CT image is necessary.
さらにまた、本発明は、被検体のX線投影データを収集するガントリと、前記X線投影データからCT画像を再構成する画像処理装置と、前記ガントリの空洞部に対して前記被検体を搬送するテーブル装置とを備えるX線CTシステムであって、
前記画像処理装置が、
前記X線投影データからフィルタ補正逆投影法を用いてCT画像の画像再構成を行う再構成手段と、
前記CT画像と所定の比較画像とを用いた逐次近似法によって、該比較画像を画質改善する画質改善手段と、
前記更新手段によって更新された前記比較画像上に所定のROIを設定して、該ROI内の標準偏差を算出する算出手段と、
前記標準偏差に基づいて前記比較画像の画質改善がさらに必要か否かを判定する判定手段と、
前記判定手段によって前記比較画像の更新がさらに必要であると判定された場合、前記CT画像と該判定手段によって判定された前記比較画像とを用いた逐次近似法によって、該比較画像を画質再改善する画質再改善手段と、
前記判定手段によって前記比較画像の画質改善が必要でないと判定された場合、該比較画像をCT画像として出力する出力手段と
を備えることを特徴とする。
Furthermore, the present invention provides a gantry that collects X-ray projection data of a subject, an image processing device that reconstructs a CT image from the X-ray projection data, and transports the subject to a cavity of the gantry. An X-ray CT system comprising a table device for performing
The image processing apparatus is
Reconstruction means for performing image reconstruction of a CT image from the X-ray projection data using a filtered back projection method;
An image quality improvement means for improving the image quality of the comparison image by a successive approximation method using the CT image and the predetermined comparison image;
Calculating means for setting a predetermined ROI on the comparison image updated by the updating means and calculating a standard deviation in the ROI;
Determining means for determining whether or not further improvement in image quality of the comparison image is necessary based on the standard deviation;
If it is determined by the determination means that the comparison image needs further updating, the image quality of the comparison image is improved again by a successive approximation method using the CT image and the comparison image determined by the determination means. Image quality re-improvement means,
And an output unit that outputs the comparison image as a CT image when it is determined by the determination unit that the image quality improvement of the comparison image is not necessary.
本発明によれば、CT画像の特定部位について、エッジ部分のシャープさを保ったまま、それ以外の部分をより滑らかに補正することができるとともに、再構成画像の評価を定量的に行うことができる。 According to the present invention, it is possible to more smoothly correct other portions of a specific portion of a CT image while maintaining the sharpness of the edge portion, and to quantitatively evaluate the reconstructed image. it can.
以下、図面を参照して、本発明の一実施形態について詳細に説明する。 Hereinafter, an embodiment of the present invention will be described in detail with reference to the drawings.
図1は、本発明の一実施形態に係るX線CTシステムの外観図である。図1に示すように、本実施形態に係るX線CTシステムは、ガントリ100とテーブル装置200と操作コンソール300とから構成される。
FIG. 1 is an external view of an X-ray CT system according to an embodiment of the present invention. As shown in FIG. 1, the X-ray CT system according to this embodiment includes a
ガントリ100は、空洞部101を有するドーナツ状の回転体102を備えている。回転体102内には、X線の発生源であるX線管と、当該X線管から発生したX線を検出するX線検出素子群で構成されるX線検出部とが、空洞部101を挟んだ対向する位置にそれぞれ収容されており、当該X線管及びX線検出器が互いの位置関係を保ったまま、空洞部101の周りを回転する構造になっている。従って、この空洞部101内に被検者を位置させた状態で、X線管及びX線検出器を備える回転体102の回転を行いながら、X線管の駆動及びX線検出器の検出動作を行うことによって、異なる投影角度からのX線照射及び被検者を透過したX線の検出を行うことが可能である。
The
また、ガントリ100には、操作パネル103が備えられており、この操作パネル103を用いて、オペレータはテーブル装置200のテーブル等のスライド操作や投影データの取得の実行等を行うことが可能である。
Further, the
また、操作コンソール300は、ガントリ100と通信する機能や、ガントリ100から転送されてきた投影データを受信して、X線断層像を再構成する機能等を有する。
The
図1に示すようなX線CTシステムを用いて実際にスキャンを行う場合、オペレータは、操作コンソール300において、スキャン条件等を設定してスキャン計画を立てる。そして、当該スキャン計画に基づいて、メインスキャンが行われる。すなわち、初期位置が設定されたテーブル装置200を用いて、被検者を載せたテーブル等がガントリ100の空洞部101に挿入され、回転体102を回転させると共にX線管を駆動しつつ、X線検出部で被検体を透過したX線を検出する。そして、回転体102の各回転角における透過X線の投影データを操作コンソール300に転送する。そして、操作コンソール300において、算術的にX線断層像を再構成する。
When an actual scan is performed using the X-ray CT system as shown in FIG. 1, the operator sets a scan condition and the like on the
図2は、本実施形態に係るX線CTシステムの細部構成を示すブロック図である。 FIG. 2 is a block diagram showing a detailed configuration of the X-ray CT system according to the present embodiment.
ガントリ100は、図2に示すように、全体の制御を行うメインコントローラ104をはじめ、以下の構成を備えている。
As shown in FIG. 2, the
105aは操作コンソール300との通信を行うためのインタフェース、105bはテーブル装置200との通信を行うためのインタフェースである。102はテーブル装置200上に横になった被検者を位置させるための空洞部101を有する回転体である。回転体102の内部には、X線管コントローラ106により駆動制御されるX線発生源であるX線管105と、X線の照射範囲を画定するためのスリットを有するコリメータ107が設けられている。
105a is an interface for communicating with the
また、回転体102には、被検者を透過したX線を検出するX線検出部108と、X線検出部108で検出された透過X線より得られる投影データを収集するデータ収集部109も備わっている。尚、X線管105及びコリメータ107とX線検出部108とは、回転体102において互いに空洞部101を挟んで反対位置、すなわち、被検者を挟んで対向する位置に設けられ、その関係が維持された状態で空洞部101の周りを回転するようになっている。回転体102の回転は、モータコントローラ110からの駆動信号により駆動される回転モータ111によって行われる。
The
112は、各種プログラム等を記録しているROMである。113は、ワークエリアとして使用されるとともに、テーブル装置200から送られてくるテーブル部の位置情報等を保持するための領域を確保したRAMである。
メインコントローラ104は、インタフェース105aを介して受信した各種コマンドの解析を行い、それに基づいて上記のX線管コントローラ106、モータコントローラ110、データ収集部109、及びテーブル装置200に対して、各種制御信号を出力する。また、メインコントローラ104は、電源投入後に、ROM112に格納された後述する図3のフローチャートに係るプログラムをRAM113にロードし、そのプログラムを実行する。さらに、メインコントローラ104は、データ収集部109で収集された投影データ及びRAM113に記憶されたデータを、インタフェース105aを介して操作コンソール300に送出する処理を行う。
The
操作コンソール300は、いわゆるワークステーション等で実現される画像処理装置であって、図2に示すように、装置全体の制御を行うCPU51、ブートプログラム等を記憶するROM52、主記憶装置として機能するRAM53をはじめとして、以下の構成を備える。
The
HDD54は、ハードディスク装置であって、OSの他に、スキャン計画を行うためのプログラム、ガントリ100に各種指示を与えたりガントリ100より受信したデータに基づいてX線断層像を再構成するための診断プログラム等が格納されている。また、VRAM55は、表示しようとするイメージデータを展開するためのメモリであり、ここにイメージデータ等を展開することでCRT56に表示させることができる。57及び58は、それぞれ各種設定を行うためのキーボード及びマウスである。また、59は、ガントリ100と通信を行うためのインタフェースである。
The HDD 54 is a hard disk device, in addition to the OS, a program for performing a scan plan, a diagnosis for giving various instructions to the
図3は、本実施形態に係るX線CTシステムにおける操作コンソール300による画像再構成処理を説明するためのフローチャートである。先に説明したとおり、図3のフローチャートに係るプログラムはROM52に格納されていたものであり、電源投入後、RAM53にロードされ、CPU51によって実行されるものである。
FIG. 3 is a flowchart for explaining image reconstruction processing by the
まず、所定角度ごとに被検体に対して照射されたX線の投影データ(scan data)がインタフェース59を介してガントリ100から取得される(ステップS1)。このとき、取得されるデータは、被検体に対して180度〜360度回転時のそれぞれの角度における投影データであって、例えば984枚の投影データが取得される。
First, projection data (scan data) of X-rays irradiated on the subject at every predetermined angle is acquired from the
次に、それぞれの投影データに対して、ノイズ除去等の前処理が行われる(ステップS2)。そして、全ての投影データを用いてFBP法を用いた画像の再構成処理が行われる(ステップS3)。そして、特定の部位に相当する関心領域(以下、「ROI」と称す。)についてCT画像(再構成画像)の画質を改善するために、後処理としてOSEM法を用いた画質改善処理が行われる(ステップS4)。尚、本実施形態における当該OSEM法を用いたROIの画質改善処理の詳細については後述する。そして、画質改善後のCT画像が取得され、CRT56上に表示される(ステップS5)。 Next, pre-processing such as noise removal is performed on each projection data (step S2). Then, image reconstruction processing using the FBP method is performed using all projection data (step S3). Then, in order to improve the image quality of the CT image (reconstructed image) for a region of interest (hereinafter referred to as “ROI”) corresponding to a specific part, image quality improvement processing using the OSEM method is performed as post-processing. (Step S4). Details of ROI image quality improvement processing using the OSEM method in the present embodiment will be described later. Then, a CT image after image quality improvement is acquired and displayed on the CRT 56 (step S5).
図4は、本実施形態で再構成されたCT画像のROIの画質改善処理手順を詳細に説明するためのフローチャートである。 FIG. 4 is a flowchart for explaining in detail the ROI image quality improvement processing procedure of the CT image reconstructed in the present embodiment.
まず、上述したステップS3におけるFBP法を用いて再構成されたCT画像Iを再投影して、再投影画像inを生成する(ステップS41)。ここで、iはFBP法を用いて作成されたCT画像Iの再投影画像を表し、nは当該CT画像に対するビュー(すなわち、角度又は位相に対応する。)を識別する番号であり、例えば再投影が100ビューの場合は、1から100までのいずれかの値を持つ。尚、このnの値は、経験値又は改善の程度に応じた任意の値を用いることができる。 First, re-projecting the reconstructed CT image I using the FBP method in step S3 described above, to produce a re-projected image i n (step S41). Here, i represents a reprojected image of the CT image I created using the FBP method, and n is a number for identifying a view (that is, corresponding to an angle or a phase) with respect to the CT image. When the projection is 100 views, it has any value from 1 to 100. The value of n can be an empirical value or an arbitrary value depending on the degree of improvement.
次に、再構成されたCT画像Iの比較対象となる比較画像I’(m)を作成する(ステップS42)。ここで、mは比較画像I’の更新回数を表し、最初は0である(ステップS43)。本実施形態では、比較画像I’(m)の初期画像I’(0)として、全画素が1の値を有するCT画像と同じ大きさの画像を作成するものとする。尚、この画素値は0以外の値であれば、それぞれの画素値がばらばらであってもかまわない。そして、FBP法を用いて再構成したCT画像Iの再投影の場合と同一のビューを用いて当該初期画像I’(0)を投影して、投影画像i’n(0)を生成する(ステップS44)。ここで、i’は、比較画像I’の投影画像を表すものとする。 Next, a comparison image I ′ (m) to be compared with the reconstructed CT image I is created (step S42). Here, m represents the number of updates of the comparison image I ′, and is initially 0 (step S43). In the present embodiment, as the initial image I ′ (0) of the comparison image I ′ (m), an image having the same size as the CT image having a value of 1 for all the pixels is created. In addition, if this pixel value is a value other than 0, each pixel value may vary. Then, the initial image I ′ (0) is projected using the same view as in the case of reprojection of the CT image I reconstructed using the FBP method to generate a projection image i ′ n (0) ( Step S44). Here, i ′ represents a projected image of the comparative image I ′.
そして、再投影画像inと投影画像i’nとを同一の角度(すなわち、nが同一のビュー)で比較することによって比較画像I’(m)を更新し、更新された比較画像I’(m+1)を生成する(ステップS45)。すなわち、画像I’(0)をI’(1)に更新する。尚、本実施形態における逐次近似法を用いた更新方法は、公知のOSEM(サブセット化による期待値最大化)法を用いて行うため、詳細な説明は省略する。次に、更新された比較画像I’(m+1)上に設定されたROI内の画素の標準偏差SDを算出する(ステップS46)。 Then, the comparison image I ′ (m) is updated by comparing the reprojection image i n and the projection image i ′ n with the same angle (that is, the view with the same n), and the updated comparison image I ′. (M + 1) is generated (step S45). That is, the image I ′ (0) is updated to I ′ (1). Note that the update method using the successive approximation method in the present embodiment is performed using a known OSEM (expected value maximization by subsetting) method, and thus detailed description thereof is omitted. Next, the standard deviation SD of the pixels in the ROI set on the updated comparison image I ′ (m + 1) is calculated (step S46).
そして、算出された標準偏差SDが一定の範囲内の値になっているかどうかを判断する(ステップS47)。その結果、標準偏差SDが一定の範囲内にない場合(No)、比較画像の更新がさらに必要であると判定され、ステップS44に戻って更新された比較画像I’(m+1)を同一ビューを用いて再投影して上述したステップS45〜S47の処理を実行する。すなわち、比較画像の更新が必要でないと判定されるまで、当該ループが繰り返され、そのたびに比較画像が更新される。一方、標準偏差SDが一定の範囲内にある場合(Yes)、終了条件を満たしたと判定され、本後処理を終了してステップS5に進む。そして、ステップS5では、終了条件を満たしたときに最終的に更新された比較画像I’が、画質の改善されたCT画像としてCRT56に表示される。
Then, it is determined whether or not the calculated standard deviation SD is a value within a certain range (step S47). As a result, if the standard deviation SD is not within a certain range (No), it is determined that the comparison image needs to be updated again, and the process returns to step S44 to display the updated comparison image I ′ (m + 1) in the same view. Then, the above-described steps S45 to S47 are executed after reprojection. That is, the loop is repeated until it is determined that it is not necessary to update the comparison image, and the comparison image is updated each time. On the other hand, when the standard deviation SD is within a certain range (Yes), it is determined that the termination condition is satisfied, the post-processing is terminated, and the process proceeds to step S5. In step S5, the comparison image I ′ that is finally updated when the end condition is satisfied is displayed on the
上述したように、操作コンソール300は、被検体のX線投影データからCT画像を再構成する画像処理装置として機能する。すなわち、当該画像処理装置は、ガントリ100でスキャンされたX線投影データからFBP法を用いてCT画像の再構成を行い、当該CT画像と所定の比較画像とを用いた逐次近似法によって、その比較画像を更新する。次いで、更新された比較画像上に所定のROIを設定して、当該ROI内の標準偏差を算出する。ここで、算出された標準偏差に基づいて比較画像の更新がさらに必要か否かが判定される。その結果、比較画像の更新がさらに必要であると判定された場合は、CT画像と比較画像とを用いた逐次近似法によって、上記比較画像を再更新する。一方、比較画像の更新が必要でないと判定された場合、当該比較画像をCT画像として出力することを特徴とする。
As described above, the
尚、本実施形態では、後処理の終了条件として標準偏差SDを用いているが、分散を用いた場合であっても同様に処理を行うことができる。すなわち、比較画像上に設定されたROI内の分散を算出し、当該分散に基づいて比較画像の更新が必要であるか否かを判定させる。 In the present embodiment, the standard deviation SD is used as the post-processing termination condition, but the same processing can be performed even when variance is used. That is, the variance within the ROI set on the comparison image is calculated, and it is determined whether or not the comparison image needs to be updated based on the variance.
また、従来から知られている反復計算回数の最大値を予め決めておく方法、或いはOSEM法で用いられる期待値の変化率が所定の範囲内に収まった場合に終了する方法と、上述した本実施形態における標準偏差SDを用いる方法とを組み合わせて上記終了条件を定めても良い。例えば、上記比較画像の更新回数を計測する計測手段をさらに備えておき、標準偏差に基づいて比較画像の更新が必要であるか否かの判定結果にかかわらず、所定の更新回数に達した場合には、比較画像の更新処理を終了させ、そのときの比較画像をCT画像として出力させる。 In addition, a method of previously determining the maximum value of the number of iterations known in the past, a method of ending when the rate of change of the expected value used in the OSEM method falls within a predetermined range, and the book described above The end condition may be determined in combination with the method using the standard deviation SD in the embodiment. For example, when a measurement means for measuring the number of updates of the comparison image is further provided and the predetermined number of updates is reached regardless of the determination result of whether or not the comparison image needs to be updated based on the standard deviation In this case, the comparison image update process is terminated, and the comparison image at that time is output as a CT image.
ここで、上記ステップS46における標準偏差SDの算出処理について詳細に説明する。本実施形態では、ステップS3で再構成したCT画像IにROIを設けて、当該ROI内の画素値について標準偏差SDを算出するものとする。そこで、ROIの設定方法として、ステップS3で生成された再構成画像から自動的に設定された固定のROIを用いる方法と、ステップS4における後処理後の画像に対してその都度設定される動的なROIを用いる方法等がある。 Here, the process of calculating the standard deviation SD in step S46 will be described in detail. In the present embodiment, it is assumed that an ROI is provided for the CT image I reconstructed in step S3, and the standard deviation SD is calculated for the pixel values in the ROI. Therefore, as a ROI setting method, a method that uses a fixed ROI that is automatically set from the reconstructed image generated in step S3 and a dynamic that is set each time for the post-processed image in step S4. There is a method using a simple ROI.
まず、上記固定のROIを用いる方法は、ステップS3で再構成された画像中においてROIを自動検出し、最後まで当該ROIを比較画像I’に対して適用して標準偏差SDを求める方法である。図5は、標準偏差SDを算出するためのROIを自動検出するための詳細な処理手順を説明するためのフローチャートである。まず、前述した図3のステップS3において再構成された画像中においてROIとして注目する領域に対して期待されるCT値の範囲で当該CT画像の2値化を行う(ステップS51)。例えば、撮影部位が腹部であって、ROIとして注目する領域が肝臓の場合、CT画像I(再構成画像)中における当該肝臓領域に期待されるCT値の範囲でCT画像Iを2値化する。 First, the method using the fixed ROI is a method in which the ROI is automatically detected in the image reconstructed in step S3, and the standard deviation SD is obtained by applying the ROI to the comparison image I ′ until the end. . FIG. 5 is a flowchart for explaining a detailed processing procedure for automatically detecting the ROI for calculating the standard deviation SD. First, the CT image is binarized in the range of CT values expected for the region of interest as ROI in the image reconstructed in step S3 of FIG. 3 (step S51). For example, when the imaging region is the abdomen and the region of interest as the ROI is the liver, the CT image I is binarized within the CT value range expected for the liver region in the CT image I (reconstructed image). .
次に、ステップS51で得られた2値画像をラベリングすることにより、2値画像中に存在するそれぞれの独立した領域をラベル付けする(ステップS52)。そして、ラベル付けされたそれぞれの領域のうち、一定以上の面積(一定数以上の画素数)を有する領域のみを抽出する(ステップS53)。次いで、抽出された領域に対して論理フィルタリングを行うことによって、当該領域の収縮処理を行ってエッジ部分を除外する(ステップS54)。これによって、エッジ部分に囲まれた内部の領域だけの画像計測値で画像ノイズを定量的に表現できる標準偏差SDや分散を求めることが可能となる。そして、収縮された領域に対応する位置の比較画像I’の領域をROIとして、当該ROIに対してステップS4のOSEM法を用いた後処理を行う。尚、この場合、上述したように、更新された比較画像I’におけるROIの標準偏差SDや分散が所定範囲内になった場合、その後処理を終了するようにする。 Next, by labeling the binary image obtained in step S51, each independent region existing in the binary image is labeled (step S52). And only the area | region which has an area more than fixed (the number of pixels more than a fixed number) among each labeled area | region is extracted (step S53). Next, by performing logical filtering on the extracted region, the region is contracted to remove the edge portion (step S54). As a result, it is possible to obtain the standard deviation SD and the variance that can quantitatively represent the image noise with the image measurement values of only the inner area surrounded by the edge portion. Then, the region of the comparison image I ′ at the position corresponding to the contracted region is set as the ROI, and post-processing using the OSEM method in step S4 is performed on the ROI. In this case, as described above, when the standard deviation SD or variance of the ROI in the updated comparison image I ′ falls within a predetermined range, the process is terminated thereafter.
また、ステップS4における後処理した後の比較画像I’に対してその都度ROIを設定する方法では、前述した図4のフローチャートにおけるステップS44の画像の更新によって更新された比較画像I’に対して、その都度、図5に示した上記手順によってROIを設定する。 Further, in the method of setting the ROI for each comparison image I ′ after the post-processing in step S4, the comparison image I ′ updated by updating the image in step S44 in the flowchart of FIG. 4 described above. Each time, the ROI is set by the above procedure shown in FIG.
上述したように、本実施形態では、CT画像(再構成画像)の全体ではなく特定のROIについてのみ標準偏差SDや分散を算出することによって、注目する部位の画像再構成を好適に行うことが可能となる。また、上述したように、ROIをCT画I(再構成画像)又は更新された比較画像I’から自動抽出して、当該ROIの標準偏差SDや分散の値に基づいて後処理である更新処理の終了条件を設定することにより、更新された比較画像、ステップS5ですなわち最終的に得られるCT画像の定量的な評価が可能となる。 As described above, in this embodiment, it is possible to suitably perform image reconstruction of a region of interest by calculating the standard deviation SD and variance only for a specific ROI, not the entire CT image (reconstructed image). It becomes possible. In addition, as described above, the ROI is automatically extracted from the CT image I (reconstructed image) or the updated comparison image I ′, and the update process is post-processing based on the standard deviation SD and variance value of the ROI. By setting the end condition, it is possible to quantitatively evaluate the updated comparison image, that is, the CT image finally obtained in step S5.
尚、ROIの設定は、CT画像Iや比較画像I’から自動検出するのではなく、あらかじめ決められた範囲の領域をROIとして用いるものであってもよい。すなわち、再構成される画像の種類に関わらず、例えば、当初から画像中の所定の領域をROIとして設定しておき、その範囲内の画素値の標準偏差SDや分散を求める方法である。この方法によれば、エッジ部の保存には適さないが、ROIを設定するための演算の必要がないので、高速に同様の終了条件を用いて更新処理を行うことが可能である。 The ROI setting may not be automatically detected from the CT image I or the comparison image I ′, but may use an area in a predetermined range as the ROI. That is, regardless of the type of image to be reconstructed, for example, a predetermined region in the image is set as the ROI from the beginning, and the standard deviation SD and variance of the pixel values within the range are obtained. According to this method, although it is not suitable for storing the edge portion, it is not necessary to perform an operation for setting the ROI, so that the update process can be performed at high speed using the same end condition.
<その他の実施形態>
OSEM法を用いた演算の反復ごとに、サブセット内の投影数を変化することにより、比較的少ない計算回数で再構成画像を改善することができることが知られている。そこで、本実施形態では、例えば、図4のステップS46で算出された標準偏差SDや分散の大きさに応じて、次のループのステップS45におけるOSEM法による画像I’の更新時のサブセット内の投影数を変化させる例について説明する。
<Other embodiments>
It is known that the reconstructed image can be improved with a relatively small number of calculations by changing the number of projections in the subset for each iteration of the operation using the OSEM method. Therefore, in the present embodiment, for example, in accordance with the standard deviation SD calculated in step S46 in FIG. 4 and the size of the variance, in the subset when the image I ′ is updated by the OSEM method in step S45 of the next loop. An example of changing the number of projections will be described.
図6は、反復演算の各段階で算出される標準偏差SDや分散によってサブセット内の投影数を変化させて画像I’を更新する場合の例を示す図である。例えば、全投影データが6(6ビュー)の場合、1回目の演算では(A)に示すように各画素は1つの投影数、6つのサブセットで投影し、画像I’を更新してROIの標準偏差SD1を求める。そして、2回目は一回目の標準偏差SD1に基づいて、各画素を2つの投影数、3つのサブセットで投影し、画像I’をさらに更新してROIの標準偏差SD2を求める。さらに、3回目は2回目の標準偏差SD2に基づいて、各画素を6つの投影数、1つのサブセットで投影し、画像I’をさらに更新してROIの標準偏差SD3を求めていく。 FIG. 6 is a diagram illustrating an example in which the image I ′ is updated by changing the number of projections in the subset based on the standard deviation SD and the variance calculated at each stage of the iterative calculation. For example, if the total projection data is 6 (6 views), in the first calculation, each pixel is projected with one projection number and 6 subsets as shown in (A), and the image I ′ is updated to update the ROI. A standard deviation SD1 is obtained. Then, in the second time, each pixel is projected with two projection numbers and three subsets based on the first standard deviation SD1, and the image I ′ is further updated to obtain the ROI standard deviation SD2. Further, in the third time, each pixel is projected with six projection numbers and one subset based on the second standard deviation SD2, and the image I 'is further updated to obtain the ROI standard deviation SD3.
上述したように、上述した実施形態のOSEM法を用いた演算の反復ごとに、サブセット内の投影数を変化することにより、高周波成分の収束を早め、少ない計算回数で良好な画像を得ることが可能となる。 As described above, by changing the number of projections in the subset for each iteration of the operation using the OSEM method of the above-described embodiment, it is possible to accelerate the convergence of high-frequency components and obtain a good image with a small number of calculations. It becomes possible.
また、OSEM法の特性として、サブセット内の投影数が少ない場合はエッジ成分が強調されてあまりノイズを軽減することができないが、サブセット内の投影数が多い場合はエッジが滑らかになる一方でノイズを軽減する程度は少ない。従って、目的に応じてサブセット内の投影数を変更して、上記実施形態で説明した標準偏差SDを用いた終了条件を用いて更新処理を行うようにしてもよい。 Further, as a characteristic of the OSEM method, when the number of projections in the subset is small, the edge component is emphasized and noise cannot be reduced much. However, when the number of projections in the subset is large, the edge becomes smooth while the noise is smoothed. There is little degree to reduce. Accordingly, the number of projections in the subset may be changed according to the purpose, and the update process may be performed using the termination condition using the standard deviation SD described in the above embodiment.
また、逐次近似アルゴリズムとして、MLEM法やOSEM法の他に、画像の先験的知識としての事前確率を組み込んだ最大事後確率(MAP−EM:Maximum A Posteriori)推定法が知られている。MAP−EM法は、尤度関数と事前確率により得られる事後確率を最大にするように画像を推定するため、MLEM法に対してノイズに対して安定である。そこで、上記実施形態に2つのパラメータを用いるMAP−EM法を併用することによって、最適なパラメータの組み合わせを決定し、エッジ部のシャープさとROI内の画像の滑らかさの両立を図ることもできる。 In addition to the MLEM method and the OSEM method, a maximum a posteriori (MAP-EM: Maximum A Posteriori) estimation method incorporating prior probabilities as image a priori knowledge is known as a successive approximation algorithm. The MAP-EM method is more stable against noise than the MLEM method because the image is estimated so as to maximize the posterior probability obtained by the likelihood function and the prior probability. Therefore, by combining the MAP-EM method using two parameters with the above-described embodiment, it is possible to determine an optimal combination of parameters and achieve both the sharpness of the edge portion and the smoothness of the image in the ROI.
また、フィルタ補正逆投影法においては、ノイズに関わる画像再構成パラメータとしては、再構成関数が考えられ、これの各周波数帯の強弱の調整で画像ノイズを同様に調整できる。 In the filtered back projection method, a reconstruction function can be considered as an image reconstruction parameter related to noise, and image noise can be similarly adjusted by adjusting the strength of each frequency band.
尚、本発明は、前述した実施形態の機能をコンピュータで実現するためのプログラムを、コンピュータ読み取り可能な記録媒体(記憶媒体)に記録して、当該記録媒体に記録されたプログラムを上記システム内のコンピュータ等が読み出して実行するようにしてもよい。また、コンピュータ読み取り可能な記録媒体には、フレキシブルディスク、CR−ROM、DVD−ROM等の可搬記録媒体、システム又はコンピュータに内蔵されているハードディスク、揮発性メモリ(RAM)等の記憶装置(記録装置)を含む。また、上記プログラムは、インターネット等のネットワークや電話回線等の通信回線を介して上記システムやコンピュータ等にダウンロードするようにしてもよい。 In the present invention, a program for realizing the functions of the above-described embodiments by a computer is recorded on a computer-readable recording medium (storage medium), and the program recorded on the recording medium is stored in the system. A computer or the like may read and execute. Computer-readable recording media include portable recording media such as a flexible disk, CR-ROM, and DVD-ROM, a storage device (recording device) such as a hard disk or volatile memory (RAM) built in the system or computer. Device). The program may be downloaded to the system, computer, or the like via a network such as the Internet or a communication line such as a telephone line.
そして、上記システムやコンピュータにおいて、読み出されたプログラムを実行することにより、前述した実施形態の機能が実現されることとなる。この場合、上記記録媒体等は、前述したフローチャートに対応するプログラムを格納している。 Then, the functions of the above-described embodiments are realized by executing the read program in the system or the computer. In this case, the recording medium or the like stores a program corresponding to the flowchart described above.
Claims (15)
前記X線投影データの前処理を行う前処理手段と、
前処理された前記X線投影データからCT画像の画像再構成を行う再構成手段と、
前記再構成手段によって画像再構成された前記CT画像上に所定のROIを設定して、該ROI内の画像計測値を算出し、その画像計測値に基づいて前記CT画像の画質改善が必要か否かを判定する判定手段と、
前記判定手段によって前記CT画像の画質改善が必要と判定された場合、逐次近似法のサブセット内の投影数、尤度関数および事後確率の少なくとも1つの再構成パラメータを変更して画像再構成することにより、前記CT画像を画質改善する画質改善手段と、
前記判定手段によって前記CT画像の画質改善が必要でないと判定された場合、該CT画像を出力する画像出力手段とを備えることを特徴とする画像処理装置。 An image processing apparatus for reconstructing a CT image from X-ray projection data of a subject,
Preprocessing means for performing preprocessing of the X-ray projection data;
Reconstruction means for performing image reconstruction of a CT image from the preprocessed X-ray projection data;
Is it necessary to set a predetermined ROI on the CT image reconstructed by the reconstruction means, calculate an image measurement value in the ROI, and improve the image quality of the CT image based on the image measurement value? Determining means for determining whether or not;
When it is determined by the determination means that the image quality of the CT image needs to be improved, the image reconstruction is performed by changing at least one reconstruction parameter of the number of projections, likelihood function, and posterior probability in the subset of the successive approximation method. An image quality improvement means for improving the image quality of the CT image,
An image processing apparatus comprising: an image output unit that outputs the CT image when it is determined by the determination unit that the image quality of the CT image is not required to be improved.
前記X線投影データからフィルタ補正逆投影法を用いてCT画像の画像再構成を行う再構成手段と、
前記CT画像と所定の比較画像とを用いた逐次近似法によって、該比較画像を画質改善する画質改善手段と、
前記画質改善手段によって画質改善された前記比較画像上に所定のROIを設定して、該ROI内の標準偏差を算出する算出手段と、
前記標準偏差に基づいて前記比較画像の画質改善がさらに必要か否かを判定する判定手段と、
前記判定手段によって前記比較画像の画質改善がさらに必要であると判定された場合、前記CT画像と前記画質改善がさらに必要であると該判定手段によって判定された前記比較画像とを用いた逐次近似法によって、該比較画像を画質再改善する画質再改善手段と、
前記判定手段によって前記比較画像の画質改善が必要でないと判定された場合、該比較画像をCT画像として出力する出力手段とを備えており、
前記画質改善手段及び前記画質再改善手段が、前記逐次近似法として、サブセット化による期待値最大化法を使用することを特徴とする画像処理装置。 An image processing apparatus for reconstructing a CT image from X-ray projection data of a subject,
Reconstruction means for performing image reconstruction of a CT image from the X-ray projection data using a filtered back projection method;
An image quality improvement means for improving the image quality of the comparison image by a successive approximation method using the CT image and the predetermined comparison image;
Calculating means for setting a predetermined ROI on the comparative image whose image quality has been improved by the image quality improving means, and calculating a standard deviation in the ROI;
Determining means for determining whether or not further improvement in image quality of the comparison image is necessary based on the standard deviation;
When it is determined that the image quality improvement of the comparison image is further required by the determination means, the successive approximation using the CT image and the comparison image determined by the determination means that the image quality improvement is further required Image quality re-improvement means for re-improvement of the image quality of the comparison image by the method,
An output unit that outputs the comparison image as a CT image when it is determined by the determination unit that image quality improvement of the comparison image is not necessary ;
The image quality improvement means and the image quality again improved means, wherein a successive approximation, the image processing apparatus characterized that you use the expectation maximization method by subsetting.
前記X線投影データからフィルタ補正逆投影法を用いてCT画像の画像再構成を行う再構成手段と、
前記CT画像と所定の比較画像とを用いた逐次近似法によって、該比較画像を画質改善する画質改善手段と、
前記画質改善手段によって画質改善された前記比較画像上に所定のROIを設定して、該ROI内の標準偏差を算出する算出手段と、
前記標準偏差に基づいて前記比較画像の画質改善がさらに必要か否かを判定する判定手段と、
前記判定手段によって前記比較画像の画質改善がさらに必要であると判定された場合、前記CT画像と前記画質改善がさらに必要であると該判定手段によって判定された前記比較画像とを用いた逐次近似法によって、該比較画像を画質再改善する画質再改善手段と、
前記判定手段によって前記比較画像の画質改善が必要でないと判定された場合、該比較画像をCT画像として出力する出力手段と、
前記標準偏差の大きさに基づいて、前記逐次近似法における前記CT画像及び前記比較画像の投影数及びサブセット数を決定するパラメータ決定手段をさらに備えることを特徴とする画像処理装置。 An image processing apparatus for reconstructing a CT image from X-ray projection data of a subject,
Reconstruction means for performing image reconstruction of a CT image from the X-ray projection data using a filtered back projection method;
An image quality improvement means for improving the image quality of the comparison image by a successive approximation method using the CT image and the predetermined comparison image;
Calculating means for setting a predetermined ROI on the comparative image whose image quality has been improved by the image quality improving means, and calculating a standard deviation in the ROI;
Determining means for determining whether or not further improvement in image quality of the comparison image is necessary based on the standard deviation;
When it is determined that the image quality improvement of the comparison image is further required by the determination means, the successive approximation using the CT image and the comparison image determined by the determination means that the image quality improvement is further required Image quality re-improvement means for re-improvement of the image quality of the comparison image by the method,
An output means for outputting the comparison image as a CT image when it is determined by the determination means that image quality improvement of the comparison image is not necessary ;
An image processing apparatus , further comprising: parameter determination means for determining the number of projections and the number of subsets of the CT image and the comparison image in the successive approximation method based on the size of the standard deviation .
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004205079A JP4535795B2 (en) | 2004-07-12 | 2004-07-12 | Image processing apparatus and X-ray CT system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004205079A JP4535795B2 (en) | 2004-07-12 | 2004-07-12 | Image processing apparatus and X-ray CT system |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2006025868A JP2006025868A (en) | 2006-02-02 |
JP4535795B2 true JP4535795B2 (en) | 2010-09-01 |
Family
ID=35892903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004205079A Expired - Lifetime JP4535795B2 (en) | 2004-07-12 | 2004-07-12 | Image processing apparatus and X-ray CT system |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4535795B2 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10210633B2 (en) | 2015-06-12 | 2019-02-19 | Hitachi, Ltd. | X-ray CT device and sequential correction parameter determination method |
Families Citing this family (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4509971B2 (en) * | 2006-06-09 | 2010-07-21 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X-ray CT system |
US8135186B2 (en) * | 2008-01-25 | 2012-03-13 | Purdue Research Foundation | Method and system for image reconstruction |
JP5406063B2 (en) * | 2010-01-27 | 2014-02-05 | 株式会社日立メディコ | Reconstruction calculation device, reconstruction calculation method, and X-ray CT apparatus |
JP5759159B2 (en) | 2010-12-15 | 2015-08-05 | 富士フイルム株式会社 | Radiation tomographic image generation method and apparatus |
WO2012127761A1 (en) * | 2011-03-22 | 2012-09-27 | 株式会社日立メディコ | X-ray ct device and method |
WO2013008702A1 (en) * | 2011-07-08 | 2013-01-17 | 株式会社 日立メディコ | Image reconstruction device and image reconstruction method |
US8712134B2 (en) * | 2011-10-18 | 2014-04-29 | Kabushiki Kaisha Toshiba | Method and system for expanding axial coverage in iterative reconstruction in computer tomography (CT) |
US8625870B2 (en) * | 2011-10-19 | 2014-01-07 | Kabushiki Kaisha Toshiba | Method and system for supplementing detail image in successive multi-scale reconstruction |
JP6000539B2 (en) * | 2011-12-12 | 2016-09-28 | 株式会社根本杏林堂 | Medical image processing system |
US8837797B2 (en) * | 2012-01-10 | 2014-09-16 | Kabushiki Kaisha Toshiba | Spatial resolution improvement in computer tomography (CT) using iterative reconstruction |
US10304217B2 (en) * | 2012-07-30 | 2019-05-28 | Toshiba Medical Systems Corporation | Method and system for generating image using filtered backprojection with noise weighting and or prior in |
US9420986B2 (en) | 2012-09-13 | 2016-08-23 | Hitachi, Ltd. | X-ray CT apparatus and X-ray CT image processing method |
JP6479317B2 (en) * | 2013-01-25 | 2019-03-06 | キヤノンメディカルシステムズ株式会社 | Medical imaging device |
JP6329490B2 (en) * | 2013-02-05 | 2018-05-23 | 株式会社日立製作所 | X-ray CT apparatus and image reconstruction method |
CN105188543B (en) * | 2013-04-08 | 2017-11-14 | 株式会社日立制作所 | X ray CT device, restructing operation device and restructing operation method |
JP6162241B2 (en) * | 2013-07-25 | 2017-07-12 | 株式会社日立製作所 | X-ray CT system |
JP6494218B2 (en) * | 2013-12-05 | 2019-04-03 | キヤノン株式会社 | Underdetermined linear equation solver, subject information acquisition device |
JP6215449B2 (en) | 2014-03-14 | 2017-10-18 | 株式会社日立製作所 | X-ray CT apparatus and processing apparatus |
US10019798B2 (en) | 2014-06-05 | 2018-07-10 | Hitachi, Ltd. | Image processing apparatus and reconstruction condition setting method |
JP6486100B2 (en) * | 2014-12-22 | 2019-03-20 | キヤノン株式会社 | Image processing apparatus, image processing method, and program |
JP6491011B2 (en) * | 2015-03-30 | 2019-03-27 | キヤノンメディカルシステムズ株式会社 | Medical image processing apparatus and medical image processing method |
WO2017154217A1 (en) * | 2016-03-11 | 2017-09-14 | 株式会社島津製作所 | Image reconstruction method, image reconstruction program, and tomography device equipped with same |
US20200175732A1 (en) * | 2017-06-02 | 2020-06-04 | Koninklijke Philips N.V. | Systems and methods to provide confidence values as a measure of quantitative assurance for iteratively reconstructed images in emission tomography |
JP7300811B2 (en) * | 2018-06-11 | 2023-06-30 | キヤノンメディカルシステムズ株式会社 | Medical information processing apparatus, medical information processing method, and program |
CN109884090B (en) * | 2019-03-07 | 2021-06-29 | 重庆大学 | CT spatial resolution measurement method for improving disk card method |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS58206726A (en) * | 1982-05-28 | 1983-12-02 | 株式会社日立製作所 | Image treating apparatus |
JPS6017568A (en) * | 1983-07-11 | 1985-01-29 | Hitachi Ltd | Method and device for processing image |
JPS61225640A (en) * | 1985-03-30 | 1986-10-07 | Toshiba Corp | Tomographic camera |
JPH05237094A (en) * | 1991-07-26 | 1993-09-17 | Yokogawa Medical Syst Ltd | Image processing method and apparatus for reducing streak-like false image |
JPH09238936A (en) * | 1996-03-07 | 1997-09-16 | Toshiba Medical Eng Co Ltd | Abnormal shadow detection system |
JPH10127622A (en) * | 1996-10-31 | 1998-05-19 | Toshiba Corp | X-ray computer tomography device |
JP2003000579A (en) * | 2001-05-10 | 2003-01-07 | Agfa Gevaert Nv | Method for retrospectively correcting inhomogeneity in radiograph |
JP2003135450A (en) * | 2001-10-31 | 2003-05-13 | Yoshihiko Nomura | Method for reducing artifacts in x-ray ct reconstructed image |
JP2004000508A (en) * | 2002-03-15 | 2004-01-08 | Koninkl Philips Electronics Nv | Method for determining objective function |
JP2004097665A (en) * | 2002-09-12 | 2004-04-02 | Hitachi Medical Corp | Blood flow dynamics analysis apparatus |
-
2004
- 2004-07-12 JP JP2004205079A patent/JP4535795B2/en not_active Expired - Lifetime
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS58206726A (en) * | 1982-05-28 | 1983-12-02 | 株式会社日立製作所 | Image treating apparatus |
JPS6017568A (en) * | 1983-07-11 | 1985-01-29 | Hitachi Ltd | Method and device for processing image |
JPS61225640A (en) * | 1985-03-30 | 1986-10-07 | Toshiba Corp | Tomographic camera |
JPH05237094A (en) * | 1991-07-26 | 1993-09-17 | Yokogawa Medical Syst Ltd | Image processing method and apparatus for reducing streak-like false image |
JPH09238936A (en) * | 1996-03-07 | 1997-09-16 | Toshiba Medical Eng Co Ltd | Abnormal shadow detection system |
JPH10127622A (en) * | 1996-10-31 | 1998-05-19 | Toshiba Corp | X-ray computer tomography device |
JP2003000579A (en) * | 2001-05-10 | 2003-01-07 | Agfa Gevaert Nv | Method for retrospectively correcting inhomogeneity in radiograph |
JP2003135450A (en) * | 2001-10-31 | 2003-05-13 | Yoshihiko Nomura | Method for reducing artifacts in x-ray ct reconstructed image |
JP2004000508A (en) * | 2002-03-15 | 2004-01-08 | Koninkl Philips Electronics Nv | Method for determining objective function |
JP2004097665A (en) * | 2002-09-12 | 2004-04-02 | Hitachi Medical Corp | Blood flow dynamics analysis apparatus |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10210633B2 (en) | 2015-06-12 | 2019-02-19 | Hitachi, Ltd. | X-ray CT device and sequential correction parameter determination method |
Also Published As
Publication number | Publication date |
---|---|
JP2006025868A (en) | 2006-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4535795B2 (en) | Image processing apparatus and X-ray CT system | |
JP4438053B2 (en) | Radiation imaging apparatus, image processing method, and computer program | |
JP5142664B2 (en) | X-ray computed tomography system | |
US6266388B1 (en) | Methods and apparatus for two-pass cone beam image reconstruction | |
US8311305B2 (en) | X-ray diagnostic apparatus, image processing apparatus, and method of calculating filter coefficients used for image formation processing in x-ray diagnostic apparatus and the like | |
US8788012B2 (en) | Methods and apparatus for automatically registering lesions between examinations | |
US20110044559A1 (en) | Image artifact reduction | |
JPH09253079A (en) | X-ray tomography system | |
KR101844513B1 (en) | Apparatus and method for processing medical image, and computer readable recording medium related to the method | |
JPH10262960A (en) | Partial volume artifact reducing method and its system | |
JP2001224586A (en) | Image mapping method and system | |
US20030161433A1 (en) | Backprojection method and X-ray CT apparatus | |
US7050527B2 (en) | Methods and apparatus for artifact reduction in cone beam CT image reconstruction | |
US8194945B2 (en) | Computer aided image acquisition and diagnosis system | |
JP2001299739A (en) | Method and apparatus for selecting retroactive reconstitution parameter | |
JP2006034785A (en) | X-ray ct image processing method, and x-ray ct apparatus | |
JP4509255B2 (en) | Perspective image creation method and apparatus | |
JP4463960B2 (en) | X-ray CT system, operation console, control method therefor, and storage medium | |
JP6321405B2 (en) | Image generating apparatus, radiation tomography apparatus, image generating method, and program | |
US6647084B1 (en) | Method and apparatus for filtering projection data of a helical scan | |
US7209580B2 (en) | Fast computed tomography method | |
US7769126B2 (en) | Computed tomography system | |
JP2009542282A (en) | Local motion compensated reconstruction of stenosis | |
US20050018889A1 (en) | Systems and methods for filtering images | |
JP4509497B2 (en) | X-ray CT system and control method thereof |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A625 | Written request for application examination (by other person) |
Free format text: JAPANESE INTERMEDIATE CODE: A625 Effective date: 20070508 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20100106 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100119 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100419 |
|
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: 20100518 |
|
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: 20100615 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130625 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 4535795 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130625 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130625 Year of fee payment: 3 |
|
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 |
|
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 |