JP5190942B2 - Energy subtraction processing apparatus and method, and program - Google Patents
Energy subtraction processing apparatus and method, and program Download PDFInfo
- Publication number
- JP5190942B2 JP5190942B2 JP2008145331A JP2008145331A JP5190942B2 JP 5190942 B2 JP5190942 B2 JP 5190942B2 JP 2008145331 A JP2008145331 A JP 2008145331A JP 2008145331 A JP2008145331 A JP 2008145331A JP 5190942 B2 JP5190942 B2 JP 5190942B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- density
- energy image
- low
- energy
- 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.)
- Active
Links
- 238000012545 processing Methods 0.000 title claims description 106
- 238000000034 method Methods 0.000 title description 51
- 238000006243 chemical reaction Methods 0.000 claims description 77
- 230000005855 radiation Effects 0.000 claims description 61
- 238000001514 detection method Methods 0.000 claims description 18
- 238000003672 processing method Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 description 26
- 210000000988 bone and bone Anatomy 0.000 description 20
- 238000007781 pre-processing Methods 0.000 description 10
- 238000001228 spectrum Methods 0.000 description 9
- 230000007423 decrease Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 230000015556 catabolic process Effects 0.000 description 5
- 238000006731 degradation reaction Methods 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 5
- 238000000926 separation method Methods 0.000 description 5
- 239000000126 substance Substances 0.000 description 5
- 230000002238 attenuated effect Effects 0.000 description 4
- 230000003247 decreasing effect Effects 0.000 description 4
- 230000006866 deterioration Effects 0.000 description 3
- 210000000038 chest Anatomy 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 239000003574 free electron Substances 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000001678 irradiating effect Effects 0.000 description 2
- 210000004872 soft tissue Anatomy 0.000 description 2
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 101150030817 HPT1 gene Proteins 0.000 description 1
- 101100462124 Oryza sativa subsp. japonica AHP1 gene Proteins 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 229910001385 heavy metal Inorganic materials 0.000 description 1
- 229910052740 iodine Inorganic materials 0.000 description 1
- 239000011630 iodine Substances 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 210000005075 mammary gland Anatomy 0.000 description 1
- 238000009607 mammography Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Description
本発明は、2回曝射法により取得された低エネルギー画像と高エネルギー画像とを用いて行われるエネルギーサブトラクション処理装置および方法ならびにプログラムに関するものである。 The present invention relates to an energy subtraction processing apparatus, method, and program that are performed using a low energy image and a high energy image acquired by a double exposure method.
従来からエネルギーの異なる2枚の放射線画像を用いてサブトラクション処理を行い、骨陰影を表した骨部画像と軟部組織を表した軟部画像とを生成することが行われている。エネルギーの異なる2枚の放射線画像を取得する方法として、銅板を2枚の放射線検出器で挟み込み、1回の曝射によりエネルギーの異なる2枚の放射線画像を取得する1回曝射法(たとえば特許文献1−3参照)と、高いエネルギーX線で撮影した高エネルギー画像と低いエネルギーX線で撮影した低エネルギー画像とを取得する2回曝射法とが提案されている(たとえば特許文献4参照)。 Conventionally, subtraction processing is performed using two radiation images having different energies to generate a bone part image representing a bone shadow and a soft part image representing a soft tissue. As a method of acquiring two radiation images having different energies, a single exposure method in which a copper plate is sandwiched between two radiation detectors and two radiation images having different energies are acquired by one exposure (for example, a patent) Document 1-3) and a two-time exposure method for acquiring a high-energy image photographed with high-energy X-rays and a low-energy image photographed with low-energy X-rays have been proposed (for example, see Patent Document 4). ).
特に、特許文献1−4には、エネルギーサブトラクション処理におけるパラメータ(荷重減算係数)の様々な設定・調整手法や画像の前処理方法が開示されている。具体的には、特許文献1には、1回曝射法において、高エネルギー画像および低エネルギー画に対して前処理する際に、ビームハードニングの影響を考慮し、高エネルギー画像の低濃度領域に対してコントラストが小さくなるような濃度調整を行い、低エネルギー画像の低濃度領域に対してコントラストが大きくなるような階調調整処理を行うことが開示されている。特許文献2には、1回曝射法において被写体の厚さに応じてパラメータを変更する方法が開示されている。特許文献3には、高エネルギー画像のヒストグラム(信号分布)の幅に基づいてパラメータを変更することが開示されている。
In particular, Patent Documents 1-4 disclose various setting / adjustment methods and image preprocessing methods for parameters (load subtraction coefficients) in energy subtraction processing. Specifically, in
一方、特許文献4には、2曝射法の高エネルギー画像と低エネルギー画像との撮影管電圧の組み合わせによって、予め用意しておいた管電圧組み合わせテーブルに対応するエネルギーサブトラクション係数を仮決定し、仮係数を用いてエネルギーサブトラクションを行い、軟部画像のエッジ強度が最小となるようにパラメータの最適化される方法が開示されている。
ところで、上述した1回曝射法の特許文献1−3に記載されているビームハードニングは2回曝射法においても生じる現象である。そこで、特許文献4に示す2回曝射法によるエネルギーサブトラクション処理を行う際に、前処理として特許文献1に開示されたような階調調整を行うことが考えられる。しかし、ビームハードニングによる影響は1回曝射法と2曝射法とにおいて全く異なるものになり、特許文献1−3に示す1回曝射法のパラメータの最適化をそのまま2回曝射法に適用しても、骨部と軟部との成分分離に失敗するという問題がある。
Incidentally, the beam hardening described in Patent Documents 1-3 of the single exposure method described above is a phenomenon that occurs even in the double exposure method. Therefore, when performing energy subtraction processing by the double exposure method shown in Patent Document 4, it is conceivable to perform gradation adjustment as disclosed in
そこで、本発明は、2回曝射法においてビームハードニングや放射線散乱等によるサブトラクション画像の画質劣化を低減することができるエネルギーサブトラクション処理装置および方法ならびにプログラムを提供することを目的とするものである。 Accordingly, an object of the present invention is to provide an energy subtraction processing apparatus, method, and program capable of reducing image quality degradation of a subtraction image due to beam hardening, radiation scattering, or the like in the double exposure method. .
本発明のエネルギーサブトラクション装置は、放射線源から異なるエネルギーの放射線を被写体に照射し撮影したときの高エネルギー画像と低エネルギー画像とを取得する画像取得手段と、画像取得手段において取得された高エネルギー画像および/または低エネルギー画像に対し濃度値の調整を行う階調処理手段と、階調処理手段により調整された高エネルギー画像および/または低エネルギー画像を用いてエネルギーサブトラクション処理を行うサブトラクション処理手段とを備え、階調処理手段が、高エネルギー画像の低濃度領域に対し濃度勾配を大きくするおよび/または低エネルギー画像の低濃度領域に対し濃度勾配を小さくする階調処理を行うものであることを特徴とするものである。 The energy subtraction device of the present invention includes an image acquisition unit that acquires a high energy image and a low energy image when a subject is irradiated with radiation of different energy from a radiation source, and a high energy image acquired by the image acquisition unit. And / or gradation processing means for adjusting density values for a low energy image, and subtraction processing means for performing energy subtraction processing using a high energy image and / or a low energy image adjusted by the gradation processing means. And the gradation processing means performs gradation processing for increasing the density gradient with respect to the low density area of the high energy image and / or reducing the density gradient with respect to the low density area of the low energy image. It is what.
本発明のエネルギーサブトラクション処理方法は、放射線源から異なるエネルギーの放射線を被写体に照射し撮影したときの高エネルギー画像と低エネルギー画像とを取得し、取得した高エネルギー画像の低濃度領域に対し濃度勾配を小さくする、および/または低エネルギー画像の低濃度領域に対し濃度勾配を大きくする階調処理を行い、調整した高エネルギー画像と低エネルギー画像とを用いてエネルギーサブトラクション処理を行うことを特徴とするものである。 The energy subtraction processing method of the present invention acquires a high energy image and a low energy image when a subject is irradiated with radiation of different energy from a radiation source, and obtains a density gradient with respect to the low density region of the acquired high energy image. And / or energy subtraction processing using the adjusted high-energy image and low-energy image, and performing gradation processing to increase the density gradient in the low-density area of the low-energy image Is.
本発明のエネルギーサブトラクション処理プログラムは、コンピュータに、放射線源から異なるエネルギーの放射線を被写体に照射し撮影したときの高エネルギー画像と低エネルギー画像とを取得し、取得した高エネルギー画像の低濃度領域に対し濃度勾配を大きくする、および/または低エネルギー画像の低濃度領域に対し濃度勾配を小さくする階調処理を行い、調整した高エネルギー画像と低エネルギー画像とを用いてエネルギーサブトラクション処理を行うことを実行させることを特徴とするものである。 The energy subtraction processing program of the present invention acquires a high energy image and a low energy image when a subject is irradiated with radiation of different energy from a radiation source and captures the image in a low density region of the acquired high energy image. In contrast, the gradation process is performed to increase the density gradient and / or to reduce the density gradient in the low density area of the low energy image, and to perform the energy subtraction process using the adjusted high energy image and the low energy image. It is characterized by being executed.
ここで、濃度値が高いもしくは低いとは各画像内の濃度値が高いもしくは低いことを意味するが、被写体の局所的な厚みが小さいときには濃度値が高くなり、被写体の局所的な厚みが大きいときには濃度値が低くなる。このため、濃度値が高いもしくは低いとは、被写体の局所的な厚みが小さいもしくは大きいことを意味する。 Here, a high or low density value means that the density value in each image is high or low, but the density value increases when the local thickness of the subject is small, and the local thickness of the subject is large. Sometimes the density value is low. For this reason, a high or low density value means that the local thickness of the subject is small or large.
また、高エネルギー画像について「低濃度領域に対し濃度勾配を大きくする」とは、中高濃度領域に比べて低濃度領域の濃度勾配を大きくすることであって、非線形な階調変換特性を有することを意味する。同様に、低エネルギー画像について「低濃度領域に対し濃度勾配を大きくする」とは、中高濃度領域に比べて低濃度領域の濃度勾配を小さくすることであって、非線形な階調変換特性を有することを意味する。さらに、「濃度勾配を大きくする」とはコントラストを大きくすることであり、「濃度勾配を小さくする」とはコントラストを小さくすることを意味する。 For high energy images, “increasing the density gradient relative to the low density area” means increasing the density gradient in the low density area as compared to the medium and high density areas, and has non-linear gradation conversion characteristics. Means. Similarly, “increasing the density gradient with respect to the low density area” for the low energy image means reducing the density gradient in the low density area as compared with the middle and high density areas, and has non-linear tone conversion characteristics. Means that. Furthermore, “increasing the density gradient” means increasing the contrast, and “decreasing the density gradient” means decreasing the contrast.
なお、エネルギーサブトラクション処理装置は、被写体の厚みを検出する厚み検出手段をさらに備えたものであってもよい。このとき、階調処理手段は、厚み検出手段により検出された被写体の厚みが大きくなればなるほど、高エネルギー画像に対する濃度勾配をより大きくするおよび/または低エネルギー画像に対する濃度勾配をより小さくする調整を行う。なお、濃度勾配をより大きくするおよび/またはより小さくするとは、予め設定された濃度勾配よりも大きくするおよび/または小さくするとの意味である。 Note that the energy subtraction processing apparatus may further include a thickness detection unit that detects the thickness of the subject. At this time, the gradation processing unit adjusts to increase the density gradient for the high energy image and / or decrease the density gradient for the low energy image as the thickness of the subject detected by the thickness detection unit increases. Do. Note that making the concentration gradient larger and / or smaller means making the concentration gradient larger and / or smaller than a preset concentration gradient.
また、厚み検出手段は、被写体の厚みを検出するものであればどのような方法であっても良く、たとえば高エネルギー画像または低エネルギー画像の信号分布の幅を用いて被写体の厚みを検出するものであってもよい。 Further, the thickness detecting means may be any method as long as it can detect the thickness of the subject. For example, the thickness detecting means detects the thickness of the subject using the width of the signal distribution of the high energy image or the low energy image. It may be.
さらに、エネルギーサブトラクション処理装置は、低エネルギー画像の取得時に放射線源に印加される管電圧を検出する管電圧検出手段をさらに有するものであってもよい。このとき、階調処理手段は、低エネルギー画像の取得時の管電圧毎に、管電圧が大きくなればなるほど濃度勾配が大きい複数の第1濃度変換テーブルおよび/または管電圧が大きくなればなるほど濃度勾配が小さい複数の第2濃度変換テーブルを有するものであり、管電圧検出手段により検出された管電圧に基づいて、複数の第1濃度変換テーブルおよび/または複数の第2濃度変換テーブルのうちいずれかの第1濃度変換テーブルおよび/または第2濃度変換テーブルを選択し、選択した第1濃度変換テーブルおよび/または第2濃度変換テーブルを用いて濃度値の調整を行うものであってもよい。 Furthermore, the energy subtraction processing apparatus may further include a tube voltage detection unit that detects a tube voltage applied to the radiation source when acquiring a low energy image. At this time, for each tube voltage at the time of acquisition of the low energy image, the gradation processing means has a plurality of first density conversion tables and / or a tube voltage that has a larger density gradient as the tube voltage increases and a density that increases as the tube voltage increases. It has a plurality of second density conversion tables with a small gradient, and based on the tube voltage detected by the tube voltage detecting means, any one of the plurality of first density conversion tables and / or the plurality of second density conversion tables. The first density conversion table and / or the second density conversion table may be selected, and the density value may be adjusted using the selected first density conversion table and / or second density conversion table.
本発明のエネルギーサブトラクション処理装置および方法ならびにプログラムによれば、放射線源から異なるエネルギーの放射線を被写体に照射し撮影したときの高エネルギー画像と低エネルギー画像とを取得し、取得した高エネルギー画像の低濃度領域に対し濃度勾配を大きくする、および/または低エネルギー画像の低濃度領域に対し濃度勾配を小さくする階調処理を行い、調整した高エネルギー画像と低エネルギー画像とを用いてエネルギーサブトラクション処理を行うことにより、被写体の局所的な厚みの相違によりビームハードニングの影響度が異なる場合であっても、エネルギーサブトラクション処理の前処理としてビームハードニングの影響を考慮した階調調整を行うことができるため、ビームハードニングの影響によるサブトラクション画像の画質の劣化を低減することができる。 According to the energy subtraction processing apparatus, method, and program of the present invention, a high energy image and a low energy image obtained by irradiating a subject with radiations of different energies from a radiation source are acquired, and the acquired high energy image is reduced. Perform gradation processing to increase the density gradient for the density area and / or reduce the density gradient to the low density area of the low energy image, and perform energy subtraction processing using the adjusted high energy image and low energy image As a result, even when the influence of beam hardening differs due to a difference in the local thickness of the subject, gradation adjustment can be performed in consideration of the influence of beam hardening as pre-processing of the energy subtraction process. Because of the effect of beam hardening It is possible to reduce the deterioration of the image quality of Rakushon image.
なお、被写体の厚みを検出する厚み検出手段をさらに備え、階調処理手段が、厚み検出手段により検出された被写体の厚みが大きくなればなるほど、高エネルギー画像に対する濃度勾配をより大きくするおよび/または低エネルギー画像に対する濃度勾配をより小さくする調整を行うものであるとき、被写体の厚み(体格)によりビームハードニングの影響度が被写体毎に異なる場合であっても、ビームハードニングの影響度に応じた前処理を施して安定したエネルギーサブトラクション結果を得ることができる。 The apparatus further includes a thickness detection unit that detects the thickness of the subject, and the gradation processing unit increases the density gradient with respect to the high-energy image as the thickness of the subject detected by the thickness detection unit increases. Depending on the degree of influence of beam hardening, even when the influence of beam hardening differs depending on the subject depending on the thickness (physique) of the subject when making adjustments to reduce the density gradient for low-energy images A stable energy subtraction result can be obtained by applying a pretreatment.
さらに、低エネルギー画像の取得時に放射線源に印加される管電圧を検出する管電圧検出手段をさらに有し、階調処理手段が、低エネルギー画像の取得時の管電圧毎に、管電圧が大きくなればなるほど濃度勾配が大きい複数の第1濃度変換テーブルおよび/または管電圧が大きくなればなるほど濃度勾配が小さい複数の第2濃度変換テーブルを有するものであり、管電圧検出手段により検出された管電圧に基づいて、複数の第1濃度変換テーブルおよび/または複数の第2濃度変換テーブルのうちいずれかの第1濃度変換テーブルおよび/または第2濃度変換テーブルを選択するものであれば、低エネルギー画像の取得時における管電圧(放射線エネルギー)の変化によりビームハードニングの影響度が撮影毎に異なる場合であって信号分布の幅の比では適正な調整を行うことができない場合であっても、ビームハードニングの影響度に応じた前処理を施して安定したエネルギーサブトラクション結果が得ることができる。 Furthermore, it further has a tube voltage detecting means for detecting a tube voltage applied to the radiation source at the time of acquiring a low energy image, and the gradation processing means increases the tube voltage for each tube voltage at the time of acquiring a low energy image. The tube has a plurality of first concentration conversion tables having a larger concentration gradient and / or a plurality of second concentration conversion tables having a smaller concentration gradient as the tube voltage increases, and the tube detected by the tube voltage detecting means. Low energy as long as one of the plurality of first density conversion tables and / or the plurality of second density conversion tables is selected from the first density conversion table and / or the second density conversion table. Signal distribution is the case where the influence of beam hardening differs for each radiograph due to changes in tube voltage (radiation energy) during image acquisition. Even if it is not possible to perform appropriate adjustment by the ratio of the width, it is possible to stable energy subtraction result is subjected to pretreatment in accordance with the degree of influence of beam hardening is obtained.
以下、図面を参照して本発明のエネルギーサブトラクション処理装置の実施の形態を詳細に説明する。図1は本発明のエネルギーサブトラクション処理装置の好ましい実施の形態を示すブロック図である。なお、図1のようなエネルギーサブトラクション処理装置1の構成は、補助記憶装置に読み込まれたエネルギーサブトラクション処理プログラムをコンピュータ(たとえばパーソナルコンピュータ等)上で実行することにより実現される。このとき、このエネルギーサブトラクション処理プログラムは、CD−ROM等の情報記憶媒体に記憶され、もしくはインターネット等のネットワークを介して配布され、コンピュータにインストールされることになる。
Hereinafter, embodiments of an energy subtraction processing apparatus of the present invention will be described in detail with reference to the drawings. FIG. 1 is a block diagram showing a preferred embodiment of the energy subtraction processing apparatus of the present invention. The configuration of the energy
エネルギーサブトラクション処理装置1は、画像取得手段10、位置合わせ手段20、サブトラクション処理手段30、画像調整手段40等を備えている。画像取得手段10は、放射線検出器により検出された放射線画像を取得するものであって、特に、異なるエネルギーを有する放射線を用いて撮影したときの高エネルギー画像(診断画像)HPと低エネルギー画像(参照画像)SPとを取得するものである。たとえば胸部正面撮影において高エネルギー画像HPは放射線源に管電圧が100〜140kVp程度印加されたときに撮影された画像であり、低エネルギー画像LPは放射線源に管電圧が50〜80kVp程度印加されたときに撮影された画像である。
The energy
階調処理手段15は、高エネルギー画像HPおよび低エネルギー画像LPに対し前処理として濃度値の階調補正を行うものである。ここで、図2Aは高エネルギー画像HPの所定の階調変換特性(入力−出力特性)を有する第1濃度変換テーブルの一例を示すグラフ、図2Bは低エネルギー画像LPの階調変換特性(入力−出力特性)を有する第2濃度変換テーブルの一例を示すグラフである。図2Aに示すように高エネルギー画像HPの低濃度領域LCHに対し濃度調整する際の濃度勾配が大きくなっており、階調処理手段15は高エネルギー画像HPの低濃度領域LCHに対しコントラストが大きくなるように調整する。一方、図2Bに示すように、低エネルギー画像LPの低濃度領域LCLに対し濃度調整する際の濃度勾配が小さくなっており、階調処理手段15は低エネルギー画像LPの低濃度領域LCLに対しコントラストが小さくなるように調整する。 The gradation processing means 15 performs gradation correction of density values as preprocessing for the high energy image HP and the low energy image LP. 2A is a graph showing an example of a first density conversion table having a predetermined gradation conversion characteristic (input-output characteristic) of the high energy image HP, and FIG. 2B is a gradation conversion characteristic (input) of the low energy image LP. 10 is a graph showing an example of a second density conversion table having (output characteristics). As shown in FIG. 2A, the density gradient when the density is adjusted with respect to the low density region LCH of the high energy image HP is large, and the gradation processing means 15 has a large contrast with respect to the low density region LCH of the high energy image HP. Adjust so that On the other hand, as shown in FIG. 2B, the density gradient at the time of adjusting the density with respect to the low density area LCL of the low energy image LP is small, and the gradation processing means 15 has a low density area LCL of the low energy image LP. Adjust the contrast to be small.
位置合わせ手段20は、低エネルギー画像LPを非線形に歪曲変形させる等の公知の位置合わせ手法を用いて高エネルギー画像HPと低エネルギー画像LPとの位置合わせを行うものである。これによりエネルギーサブトラクション処理により軟部組織と骨との分離性能を向上させることができる。
The
サブトラクション処理手段30は、位置合わせ手段20により位置合わせ処理が施された高エネルギー画像HPと低エネルギー画像LPとを用いてエネルギーサブトラクション処理を行うものである。サブトラクション処理手段30は、高エネルギー画像HPと低エネルギー画像LPとを用いて、骨を除去した被写体の軟部を示す軟部画像SPと被写体の骨部を示す骨部画像BPとを生成する機能を有している。一般的に、サブトラクション画像Psubは、第1荷重減算係数Kaを積算した高エネルギー画像HPと第2荷重減算係数Kbを積算した低エネルギー画像LPとの差分により表される。
The
Psub=Ka・HP−Kb・LP+Kc ・・・(1)
なお、式(1)中のKaは第1荷重減算係数、Kbは第2荷重減算係数、Kcは所定のオフセット値である。そして、サブトラクション処理手段30は、所定の係数Ka、Kb、Kcを用いて式(1)を算出することにより軟部画像SPをサブトラクション画像Psubとして生成する。その後、サブトラクション処理手段30は高エネルギー画像HPから軟部画像SPを減算処理することにより骨部画像BPを取得する(BP=HP−SP)。
Psub = Ka · HP−Kb · LP + Kc (1)
In Equation (1), Ka is a first load subtraction coefficient, Kb is a second load subtraction coefficient, and Kc is a predetermined offset value. Then, the
このように、エネルギーサブトラクション処理の前処理として、高エネルギー画像HPの低濃度領域に対し濃度勾配が大きくなるように階調調整され、あるいは低エネルギー画像LPの低濃度領域に対し濃度勾配が小さくなるように階調調整されることにより、ビームハードニングの影響によるコントラストの低下を防止し、サブトラクション画像Psub(軟部画像SPおよび骨部画像BP)の画質劣化を低減することができる。具体的には、2回曝射法の場合、ビームハードニングの影響により被写体の厚みが大きい部分ほど高エネルギー画像HPと低エネルギー画像LPとのコントラスト比は大きくなる。すなわち、放射線が物質を透過するときには放射線の光子数(線量)は減弱するものであり、放射線は低エネルギーなほど減弱係数が大きく、エネルギーが高くなるにつれて減弱係数が小さくなり透過性が増す。放射線の低エネルギー成分の方が高エネルギー成分よりも大きく減弱するため、放射線は物質を透過するほど低エネルギー成分が減弱されて、高エネルギー成分の割合が多くなる(平均エネルギーが高くなる)。この現象はビームハードニングと呼ばれる。図3は骨部と軟部とにおいても放射線のエネルギー(管電圧)が大きくなればなるほど減弱係数が小さくなることを示している。 As described above, as pre-processing of the energy subtraction process, gradation adjustment is performed so that the density gradient is increased with respect to the low density area of the high energy image HP, or the density gradient is decreased with respect to the low density area of the low energy image LP. By adjusting the gradation as described above, it is possible to prevent a decrease in contrast due to the effect of beam hardening and to reduce image quality degradation of the subtraction image Psub (soft part image SP and bone part image BP). Specifically, in the case of the double exposure method, the contrast ratio between the high energy image HP and the low energy image LP increases as the thickness of the subject increases due to the effect of beam hardening. That is, when radiation passes through a substance, the number of photons (dose) of the radiation decreases. The lower the energy, the larger the attenuation coefficient, and the higher the energy, the smaller the attenuation coefficient and the greater the transparency. Since the low energy component of radiation is attenuated more than the high energy component, the lower the energy component of the radiation is attenuated, the higher the ratio of the high energy component (the higher the average energy). This phenomenon is called beam hardening. FIG. 3 shows that the attenuation coefficient decreases as the radiation energy (tube voltage) increases also in the bone part and the soft part.
また、物質の原子番号が大きいほどさらには密度が高いほど、放射線の減弱係数は大きくなる。つまり、放射線の通過経路上に存在する物質の割合に応じてビームハードニングの程度も異なり、同一被写体であっても局所的な厚みが大きい領域ほど強くビームハードニングの影響を受ける。具体的には、図4に示すように、被写体に照射される放射線のエネルギースペクトルに対し、ビームハードニングの影響が大きいときには減弱する割合が大きく、影響が小さいときには減弱する割合が少ない。 In addition, the higher the atomic number of a substance and the higher the density, the greater the attenuation coefficient of radiation. In other words, the degree of beam hardening differs depending on the ratio of the substance present on the radiation passage path, and even in the same subject, the region having a larger local thickness is more strongly affected by beam hardening. Specifically, as shown in FIG. 4, the rate of attenuation is large when the effect of beam hardening is large, and the rate of attenuation is small when the effect is small, with respect to the energy spectrum of radiation irradiated to the subject.
上述したビームハードニングが生じる要因として光電効果とコンプトン効果(コンプトン散乱または非干渉性散乱)とが上げられる。光電効果は、原子核との結合の強い内殻電子を自由電子(光電子)として飛び出させる現象である。このとき、入射した放射線光子はその全エネルギーを失って消滅し、空席ができた内殻電子軌道にはその外側の軌道電子が遷移する。そのため、その前後のエネルギー準位間の差のエネルギーをもった特性X線が二次X線として放射される(蛍光放射線として知られる)。図5に示すように、光電効果は光子エネルギーの比較的低い放射線(およそ50keV以下)で顕著に現れ、物質の密度や原子番号が大きいほど大きくなる。 The photoelectric effect and the Compton effect (Compton scattering or incoherent scattering) are raised as factors causing the beam hardening described above. The photoelectric effect is a phenomenon in which inner-shell electrons having strong bonds with atomic nuclei are ejected as free electrons (photoelectrons). At this time, the incident radiation photon loses all its energy and disappears, and the outer orbital electrons transition to the vacant inner electron trajectory. For this reason, characteristic X-rays having energy between the energy levels before and after that are emitted as secondary X-rays (known as fluorescent radiation). As shown in FIG. 5, the photoelectric effect appears remarkably with radiation having a relatively low photon energy (about 50 keV or less), and becomes larger as the density and atomic number of the substance increases.
コンプトン効果は放射線の粒子としての性質が現れた現象であって、入射した放射線光子が自由電子または原子核との結合の弱い外殻電子に衝突したとき、その電子にエネルギーの一部を与えて反跳電子として軌道外に跳ね飛ばし、自分はその分だけエネルギーの小さい、波長の長い放射線になって方向を変える現象である。コンプトン効果は、光子エネルギーと物質の密度に支配される現象(密度が大きいほど強い)であり、図5に示すように、比較的高エネルギーの放射線領域(50keV程度以上)では、原子番号の特に高い物質(重金属やヨウ素など)を除き、光電効果よりも強く現れる。コンプトン散乱によって画像のコントラストは低下する。 The Compton effect is a phenomenon that appears as a particle of radiation. When an incident radiation photon collides with a free electron or an outer shell electron with a weak bond with a nucleus, it gives a part of energy to the electron and counteracts it. This is a phenomenon that jumps out of the orbit as a jumping electron, and changes its direction by becoming a long wavelength radiation with a smaller energy. The Compton effect is a phenomenon controlled by photon energy and material density (the higher the density, the stronger). As shown in FIG. 5, in the relatively high energy radiation region (about 50 keV or more), the atomic number Except for high substances (such as heavy metals and iodine), it appears stronger than the photoelectric effect. Compton scattering reduces the image contrast.
図6は2回曝射法において、低エネルギー画像LPを60kVp、高エネルギー画像HPを120kVpで撮影し取得したときに、被写体の局所的な厚さが小さい領域(図6A)と被写体の局所的な厚さが大きい領域(図6B)とについてのグラフである。図6Aの被写体の局所的な厚さが小さい場合において、低エネルギー画像LPと高エネルギー画像HPのエネルギースペクトルは近くなる。たとえば低エネルギー画像LPの平均エネルギーが33.9keVであるのに対し、高エネルギー画像HPの平均エネルギーは45.9keVになり、平均エネルギーの差は小さい。従って、高エネルギー画像HPと低エネルギー画像LPのコントラスト比は小さくなる。 FIG. 6 shows a region where the local thickness of the subject is small (FIG. 6A) and the local area of the subject when the low-energy image LP is taken and acquired at 60 kVp and the high-energy image HP is obtained at 120 kVp in the double exposure method. It is a graph about the area | region (FIG. 6B) with a large thickness. When the local thickness of the subject in FIG. 6A is small, the energy spectra of the low energy image LP and the high energy image HP are close to each other. For example, while the average energy of the low energy image LP is 33.9 keV, the average energy of the high energy image HP is 45.9 keV, and the difference in average energy is small. Accordingly, the contrast ratio between the high energy image HP and the low energy image LP becomes small.
一方、図6Bの被写体の局所的な厚さが大きい場合において、被写体でのビームハードニングが強く起こり、高エネルギー画像HPと低エネルギー画像LPとのエネルギースペクトルは離れる。このとき、低エネルギー画像LPの平均エネルギーがたとえば42.9keVであるのに対し、高エネルギー画像HPの平均エネルギーは63.4keVになり、平均エネルギーの差が大きいことがわかる。従って、高エネルギー画像HPと低エネルギー画像LPのコントラスト比は大きくなる。 On the other hand, when the local thickness of the subject in FIG. 6B is large, beam hardening occurs strongly in the subject, and the energy spectra of the high energy image HP and the low energy image LP are separated. At this time, the average energy of the low energy image LP is, for example, 42.9 keV, whereas the average energy of the high energy image HP is 63.4 keV, indicating that the difference in average energy is large. Accordingly, the contrast ratio between the high energy image HP and the low energy image LP is increased.
図6Aおよび図6Bに示すように、被写体の局所的な厚みが小さい部分ではビームハードニングの影響はそれほど起こらないが、被写体の局所的な厚みが大きい部分では強くビームハードニングを起こすことがわかる。つまり、被写体の局所的な厚みが増すにつれて高エネルギー画像HPと低エネルギー画像LPのエネルギースペクトルは離れていき、高エネルギー画像HPと低エネルギー画像LPのコントラスト比は大きくなる。 As shown in FIGS. 6A and 6B, the effect of beam hardening does not occur so much in a portion where the local thickness of the subject is small, but it is strongly caused in a portion where the local thickness of the subject is large. . That is, as the local thickness of the subject increases, the energy spectrum of the high energy image HP and the low energy image LP are separated from each other, and the contrast ratio between the high energy image HP and the low energy image LP increases.
ところで、被写体の局所的な厚みが大きいときには各エネルギー画像LP、HP内の濃度値は低くなり、被写体の局所的な厚みが小さいときには濃度値が高くなる。したがって、上述した厚みによるビームハードニングの影響度は画素濃度値の高低によって判別することができる。そこで、図2Aに示すように階調処理手段15は高エネルギー画像HPの低濃度領域LCHにおいて、コントラストが大きくなるように濃度勾配を大きくする。一方、図2Bに示すように階調処理手段15は低エネルギー画像LPの低濃度領域LCLに対してはコントラストが小さくなるように濃度勾配を小さくする。すると、高エネルギー画像HPは局所的な厚みの大きい低濃度領域LCLが局所的な厚みの小さい中高濃度領域に比べてよりコントラストが大きくなるように調整される。これにより、ビームハードニングの影響によりサブトラクション画像のコントラストの劣化を防止することができる。
By the way, when the local thickness of the subject is large, the density value in each energy image LP, HP is low, and when the local thickness of the subject is small, the density value is high. Therefore, the degree of influence of beam hardening due to the above-described thickness can be determined by the level of the pixel density value. Therefore, as shown in FIG. 2A, the
なお、この関係は1回曝射法とは逆になるものであり、従来1回曝射法において適用された非線形テーブルを用いても成分分離をうまく行うことができない。具体的には、図7は横軸に高エネルギー画像HPの濃度値、縦軸に低エネルギー画像LPの濃度値を取り、その頻度をグラフ化したジョイントヒストグラムであって、図7Aが2回曝射法(2ショット法)によるもの、図7Bが1回曝射法(1ショット法)によるものを示している。なお、濃度は0〜127に正規化されているものであって、高エネルギー画像HPと低エネルギー画像LPのコントラスト比に対応してビームハードニングの影響がない場合は直線になる。 Note that this relationship is opposite to the one-time exposure method, and component separation cannot be performed well even using a nonlinear table applied in the conventional one-time exposure method. Specifically, FIG. 7 is a joint histogram in which the horizontal axis represents the density value of the high energy image HP and the vertical axis represents the density value of the low energy image LP, and the frequency is graphed. FIG. 7B shows a result obtained by a one-shot method (one-shot method). The density is normalized to 0 to 127, and becomes a straight line when there is no influence of beam hardening corresponding to the contrast ratio between the high energy image HP and the low energy image LP.
図7Aの2回曝射法の場合、低エネルギー画像LPと高エネルギー画像HPとのコントラスト比はビームハードニングの影響により濃度値が小さい領域において直線の下側にずれた非線形なものになる。一方、図7Bの1回曝射法の場合、低エネルギー画像LPと高エネルギー画像HPとのコントラスト比はビームハードニングの影響により濃度値が小さい領域において直線の上側にずれた非線形なものになる。このように、1回曝射法と2回曝射法では特性が逆となっていることがわかる。したがって、従来1回曝射法において適用されている非線形テーブルを用いても成分分離をうまく行うことができないことがわかる。 In the case of the double exposure method of FIG. 7A, the contrast ratio between the low energy image LP and the high energy image HP is non-linearly shifted to the lower side of the straight line in the region where the density value is small due to the effect of beam hardening. On the other hand, in the case of the single exposure method of FIG. 7B, the contrast ratio between the low energy image LP and the high energy image HP is non-linearly shifted to the upper side of the straight line in the region where the density value is small due to the effect of beam hardening. . Thus, it can be seen that the characteristics are reversed between the single exposure method and the double exposure method. Therefore, it can be seen that component separation cannot be performed well even using a non-linear table applied in the conventional single exposure method.
図8は本発明のエネルギーサブトラクション処理方法の好ましい実施形態を示すフローチャートであり、図1から図8を参照してエネルギーサブトラクション処理方法について説明する。まず、画像取得手段10において高エネルギー画像HPと低エネルギー画像LPとが取得される(ステップST1)。その後、階調処理手段15により高エネルギー画像HPおよび低エネルギー画像LPの階調が調整される(ステップST2、図2A、図2B参照)。このとき、階調処理手段15において、高エネルギー画像HPの低濃度領域LCHの濃度勾配が大きくなるように調整され、もしくは低エネルギー画像LPの低濃度領域LCLの濃度勾配が小さくなるような調整が行われる。 FIG. 8 is a flowchart showing a preferred embodiment of the energy subtraction processing method of the present invention. The energy subtraction processing method will be described with reference to FIGS. First, the high-energy image HP and the low-energy image LP are acquired by the image acquisition means 10 (step ST1). Thereafter, the gradation processing means 15 adjusts the gradation of the high energy image HP and the low energy image LP (see step ST2, FIG. 2A, FIG. 2B). At this time, the gradation processing means 15 is adjusted so that the density gradient of the low density region LCH of the high energy image HP is increased or the density gradient of the low density region LCL of the low energy image LP is reduced. Done.
さらに、位置合わせ手段20により低エネルギー画像LPが非線形に歪曲変形することにより、高エネルギー画像HPと低エネルギー画像LPとの位置合わせが行われる(ステップST3)。その後、サブトラクション処理手段30において上記式(1)に基づいて軟部画像SPが生成される(ステップST4)。その後、サブトラクション処理手段30により高エネルギー画像HPから軟部画像SPを減算することにより骨部画像BPが生成される(ステップST5)。その後、画像調整手段40により平滑化によるノイズ除去処理等が施され(ステップST6)、軟部画像SPおよび骨部画像BPが表示装置等に表示され、もしくはハードディスク装置等の記憶手段に記憶される。 Further, the low energy image LP is nonlinearly distorted and deformed by the alignment means 20, thereby aligning the high energy image HP and the low energy image LP (step ST3). Thereafter, the soft traction image SP is generated in the subtraction processing means 30 based on the above formula (1) (step ST4). Thereafter, the bone image BP is generated by subtracting the soft part image SP from the high energy image HP by the subtraction processing means 30 (step ST5). Thereafter, a noise removal process by smoothing is performed by the image adjusting unit 40 (step ST6), and the soft part image SP and the bone part image BP are displayed on a display device or the like, or stored in a storage unit such as a hard disk device.
このように、エネルギーサブトラクション処理の前処理として、高エネルギー画像HPの低濃度領域に対し濃度勾配が大きくなるように階調調整され、あるいは低エネルギー画像LPの低濃度領域に対し濃度勾配が小さくなるように階調調整されることにより、ビームハードニングの影響によるコントラストの低下を防止し、サブトラクション画像Psub(軟部画像SPおよび骨部画像BP)の画質劣化を低減することができる。図9は本発明のエネルギーサブトラクション処理装置の第2の実施形態を示すブロック図であり、図9を参照してエネルギーサブトラクション処理装置100について説明する。なお、図9のエネルギーサブトラクション処理装置100において図1のエネルギーサブトラクション処理装置1と同一の構成を有する部位には同一の符号を付してその説明を省略する。図9のエネルギーサブトラクション処理装置100が図1のエネルギーサブトラクション処理装置1と異なる点は、被写体の厚み(体格)に基づいて階調処理手段における濃度変換特性が調整される点である。
As described above, as pre-processing of the energy subtraction process, gradation adjustment is performed so that the density gradient is increased with respect to the low density area of the high energy image HP, or the density gradient is decreased with respect to the low density area of the low energy image LP. By adjusting the gradation as described above, it is possible to prevent a decrease in contrast due to the effect of beam hardening and to reduce image quality degradation of the subtraction image Psub (soft part image SP and bone part image BP). FIG. 9 is a block diagram showing a second embodiment of the energy subtraction processing apparatus of the present invention. The energy
図9のエネルギーサブトラクション処理装置100は、被写体の厚みを検出する厚み検出手段140を有している。この厚み検出手段140は被写体の厚みを検出するものであり、特に、高エネルギー画像HPのスケールファクター(ヒストグラム幅の逆数)GpH値および低エネルギー画像LPのスケールファクター(ヒストグラム幅の逆数)GpL値の比GpH/GpLを厚み情報として検出する機能を有している。
The energy
以下に、比GpH/GpLが厚み情報を示すことについて説明する。図10Aは厚みが小さい(体格が小さい)被写体を撮影したときに取得された低エネルギー画像LPと高エネルギー画像HPのヒストグラムを示すグラフ、図10Bは厚みが大きい(体格が大きい)被写体を撮影したときに取得された低エネルギー画像LPと高エネルギー画像HPのヒストグラムを示すグラフである。図10Aおよび図10Bに示すように、ヒストグラムの幅の変化に着目すると、高エネルギー画像HPは体格の大小に対する変動が小さいが、低エネルギー画像LPは体格の大小に対する変動が大きい。これは被写体の厚みが大きい場合は厚みが薄い場合に比べて放射線検出器へ露光される被写体を透過したX線の放射線量が減少するからである。 Below, it demonstrates that ratio GpH / GpL shows thickness information. FIG. 10A is a graph showing a histogram of a low energy image LP and a high energy image HP acquired when a subject with a small thickness (small physique) is photographed, and FIG. 10B shows a subject with a large thickness (large physique). It is a graph which shows the histogram of the low energy image LP and the high energy image HP which were sometimes acquired. As shown in FIGS. 10A and 10B, focusing on the change in the width of the histogram, the high energy image HP has a small variation with respect to the size of the physique, whereas the low energy image LP has a large variation with respect to the size of the physique. This is because when the subject is thick, the amount of X-ray radiation transmitted through the subject exposed to the radiation detector is reduced as compared to when the subject is thin.
つまり、高エネルギー画像HPのスケールファクター値GpHは被写体の厚さにほとんど依存しないものであるが、低エネルギー画像LPのスケールファクター値GpLは被写体の厚さに依存して大きく変化する。つまり、被写体の厚さが大きいほどヒストグラムの幅が大きくなりその逆数であるGpL値が小さくなりGp値の比GpH/GpLは大きくなる。従って、Gp値の比GpH/GpLは被写体の厚さを示す厚さ情報として用いることができる。 That is, the scale factor value GpH of the high-energy image HP hardly depends on the thickness of the subject, but the scale factor value GpL of the low-energy image LP varies greatly depending on the thickness of the subject. That is, as the thickness of the subject increases, the width of the histogram increases, the GpL value that is the reciprocal thereof decreases, and the Gp value ratio GpH / GpL increases. Therefore, the Gp value ratio GpH / GpL can be used as thickness information indicating the thickness of the subject.
なお、スケールファクターGp値の比GpH/GpLにより体格を推定する場合に限らず、たとえば患者の身長、体重あるいはBMI(Body Mass Index)を用いてもよいし、超音波距離計を用いてX線源と患者の距離を測定することで、被写体の厚みを直接測定するようにしてもよいし、患者の体格を撮影者(医師、あるいは放射線技師)が手入力してもよい。厚み検出手段140がスケールファクターGp値の比GpH/GpL以外の方法で厚み情報を検出したとき、たとえば体格が大きい人、標準的な人、小さい人というように厚みを検出したとき、それぞれ体格に合わせて対応付けされた複数の調整値(GpH/GpL値)を有しており、入力された体格情報によって調整値を選択するようにしてもよい。 In addition, it is not limited to the case of estimating the physique based on the ratio GpH / GpL of the scale factor Gp value. For example, the height, weight or BMI (Body Mass Index) of the patient may be used, By measuring the distance between the source and the patient, the thickness of the subject may be directly measured, or the physique of the patient may be manually input by the photographer (doctor or radiographer). When the thickness detection means 140 detects thickness information by a method other than the ratio GpH / GpL of the scale factor Gp value, for example, when the thickness is detected such as a large physique, a standard person, or a small person, A plurality of adjustment values (GpH / GpL values) associated with each other may be included, and the adjustment value may be selected based on the inputted physique information.
そして、階調処理手段115は、厚み検出手段140により検出された厚み情報に基づいて階調調整する際の濃度勾配を調整する機能を有している。具体的には、図11Aに示すように、階調処理手段115は、たとえば予め設定されている濃度勾配(図2A参照)にGpH/GpL値を乗算して、高エネルギー画像HPに対する濃度変換特性の全領域について濃度勾配を厚みの大きさに比例してより大きくするように調整する。つまり、被写体の厚みが大きい場合には勾配が大きくなり、被写体の厚みが小さい場合には勾配が小さくなる。
The
同様に、図11Bに示すように、階調処理手段115は、たとえば予め設定されている濃度勾配(図2B参照)にGpH/GpL値を乗算して、低エネルギー画像LPに対する濃度変換特性の全領域の勾配を厚みの大きさに比例してより小さくするように調整する。つまり、被写体の厚みが大きい場合には勾配が小さくなり、被写体の厚みが小さい場合には勾配が大きくなる。 Similarly, as shown in FIG. 11B, the gradation processing means 115 multiplies the GpH / GpL value by, for example, a preset density gradient (see FIG. 2B) to obtain all density conversion characteristics for the low energy image LP. The gradient of the area is adjusted to be smaller in proportion to the thickness. That is, when the subject thickness is large, the gradient is small, and when the subject thickness is small, the gradient is large.
このように、被写体の厚み(体格)に応じて低濃度領域側の濃度勾配を調整することにより、体格差によりビームハードニングによる影響の大小が生じることによるサブトラクション画像のコントラストの低下を低減することができる。 In this way, by adjusting the density gradient on the low density area side according to the thickness (physique) of the subject, the reduction in the contrast of the subtraction image due to the influence of beam hardening caused by the physique difference is reduced. Can do.
図12は本発明のエネルギーサブトラクション処理装置の第3の実施形態を示すブロック図であり、図12を参照してエネルギーサブトラクション処理装置について説明する。なお、図12のエネルギーサブトラクション処理装置200において図9のエネルギーサブトラクション処理装置100と同一の構成を有する部位には同一の符号を付してその説明を省略する。図12のエネルギーサブトラクション処理装置200が図9のエネルギーサブトラクション処理装置100と異なる点は、放射線源に印加された管電圧値に応じて異なる濃度変換テーブルを用いる点である。
FIG. 12 is a block diagram showing a third embodiment of the energy subtraction processing apparatus of the present invention. The energy subtraction processing apparatus will be described with reference to FIG. In the energy
図12のエネルギーサブトラクション処理装置200は、管電圧を検出する管電圧検出手段240を備えている。管電圧検出手段240は、たとえば放射線源を制御する図示しない放射線制御部から放射線源に印加した管電圧値を取得するものである。
The energy
一方、テーブルデータベースには異なる管電圧毎に複数の第1濃度変換テーブルHPT11〜HPT13、複数の第2濃度変換テーブルLPT11〜LPT13が記憶されており、階調処理手段215は、管電圧検出手段240により検出された低エネルギー画像LPの取得時における管電圧値に基づいて、いずれかの第1濃度変換テーブルHPT11〜HPT13および/または第2濃度変換テーブルLPT11〜LPT13を選択する。そして、階調処理手段215は、選択した各濃度変換テーブルHPT11〜HPT13、LPT11〜LPT13を用いて階調調整を行う。
On the other hand, a plurality of first density conversion tables HPT11 to HPT13 and a plurality of second density conversion tables LPT11 to LPT13 are stored for each different tube voltage in the table database, and the gradation processing means 215 is a tube voltage detection means 240. One of the first density conversion tables HPT11 to HPT13 and / or the second density conversion tables LPT11 to LPT13 is selected based on the tube voltage value at the time of acquiring the low energy image LP detected by the above. The
具体的には、高エネルギー画像HPに対し、図13A〜Cに示すようにたとえば60、70、80kVpに対応して3つの第1濃度変換テーブルHPT11〜HPT13が記憶されている。3つの第1濃度変換テーブルHPT11〜HPT13において、曲線の勾配はほぼ同一のものであるが、低エネルギー画像の取得時の管電圧が高い第1濃度変換テーブルHPT11〜HPT13ほど濃度勾配が大きくなるように設定される。 Specifically, for the high energy image HP, as shown in FIGS. 13A to 13C, for example, three first density conversion tables HPT11 to HPT13 are stored corresponding to 60, 70, and 80 kVp. In the three first density conversion tables HPT11 to HPT13, the gradients of the curves are substantially the same, but the first gradient conversion tables HPT11 to HPT13 having a higher tube voltage at the time of acquiring a low energy image have a larger density gradient. Set to
一方、低エネルギー画像LPに対し、図14A〜Cに示すようにたとえば60、70、80kVpに対応して3つの第2濃度変換テーブルLPT11〜LPT13が記憶されている。3つの第2濃度変換テーブルLPT11〜LPT13において、曲線の勾配はほぼ同一のものであるが、低エネルギー画像の取得時の管電圧が高い第2濃度変換テーブルLPT11〜LPT13ほど濃度勾配が小さくなるように設定される。 On the other hand, for the low energy image LP, as shown in FIGS. 14A to 14C, for example, three second density conversion tables LPT11 to LPT13 corresponding to 60, 70, and 80 kVp are stored. In the three second density conversion tables LPT11 to LPT13, the slopes of the curves are substantially the same, but the density gradients become smaller as the second density conversion tables LPT11 to LPT13 having a higher tube voltage when acquiring a low energy image. Set to
そして、階調処理手段215は、低エネルギー画像LPの管電圧によって第1濃度変換テーブルHPT11〜HPT13もしくは第2濃度変換テーブルLPT11〜LPT13のいずれを使用するかを決定する。その後、階調処理手段215は決定したいずれかの第1濃度変換テーブルHPT11〜HPT13もしくは第2濃度変換テーブルLPT11〜LPT13を用いて濃度調整処理を行う。さらに、階調処理手段215は、厚み検出手段140により検出した厚みを用いて、選択した各濃度変換テーブルHPT11〜HPT13、LPT11〜LPT13の濃度勾配を調整し、調整した各濃度変換テーブルHPT11〜HPT13、LPT11〜LPT13を用いて階調の調整を行う(図9〜図11参照)。
Then, the
このように、低エネルギー画像LPの取得時における管電圧により複数の第1濃度変換テーブルHPT11〜HPT13もしくは複数の第2濃度変換テーブルLPT11〜LPT13から選択して階調調整を行うことにより、低エネルギー画像LPにおけるビームハードニングの影響によるサブトラクション画像の画質劣化を低減することができる。 As described above, the tone voltage is selected from the plurality of first density conversion tables HPT11 to HPT13 or the plurality of second density conversion tables LPT11 to LPT13 according to the tube voltage at the time of acquisition of the low energy image LP, thereby reducing the low energy. It is possible to reduce image quality degradation of the subtraction image due to the effect of beam hardening in the image LP.
以下に、複数の濃度変換テーブルを用意した理由について説明する。はじめに、高エネルギー画像HPは患者の体格や撮影目的に応じて管電圧を変更して撮影されるものであって、管電圧が高くなると透過性の高いエネルギ―成分の割合が増える。このため、メインヒストグラムの対数露光量は増加するとともにヒストグラム幅は狭くなり、高エネルギー画像HPのコントラストは小さくなる。つまりヒストグラム幅の逆数であるスケールファクターGp値は大きくなる。 The reason why a plurality of density conversion tables are prepared will be described below. First, the high-energy image HP is imaged by changing the tube voltage according to the patient's physique and the imaging purpose. As the tube voltage increases, the ratio of highly transmissive energy components increases. For this reason, the logarithmic exposure amount of the main histogram increases, the histogram width becomes narrower, and the contrast of the high energy image HP becomes smaller. That is, the scale factor Gp value, which is the reciprocal of the histogram width, increases.
図15は高エネルギー画像HPを取得する際の管電圧を100〜140kVpに変化させて同一被写体を撮影した場合のヒストグラムを示すグラフである。図15に示すように、管電圧が高くなるほど高エネルギー成分のX線の増加によって透過性が増す。このため、放射線画像のヒストグラムのうち直接X線が照射された領域(素抜け部分)のヒストグラムの山と被写体を透過した部分のメインヒストグラムの山が近づく。よって、ヒストグラムの幅が小さくなりスケールファクターGp値は大きくなる。 FIG. 15 is a graph showing a histogram when the same subject is photographed while changing the tube voltage when acquiring the high energy image HP to 100 to 140 kVp. As shown in FIG. 15, the higher the tube voltage, the higher the transmission due to the increase in the X-rays of high energy components. For this reason, in the histogram of the radiographic image, the histogram peak of the region directly irradiated with X-rays (the missing portion) approaches the peak of the main histogram of the portion that has passed through the subject. Therefore, the width of the histogram is reduced and the scale factor Gp value is increased.
ここで、低エネルギー画像LPを取得する際の管電圧を60kVpに固定させた場合、Gp値の比GpH/GpLは、高エネルギー画像HPの管電圧が高いほど大きくなる。但し、この場合には上述したようにGp値の比GpH/GpLにより調整することにより、管電圧変化による分離性能の劣化を低減することができる(図9〜図11参照)。つまり、高エネルギー画像HPの管電圧変化に対しては、高エネルギー画像HPのGp値の比GpH/GpLを用いて調整することができる。 Here, when the tube voltage at the time of acquiring the low energy image LP is fixed to 60 kVp, the Gp value ratio GpH / GpL increases as the tube voltage of the high energy image HP increases. However, in this case, as described above, the adjustment of the Gp value ratio GpH / GpL can reduce the degradation of the separation performance due to the change in tube voltage (see FIGS. 9 to 11). That is, the tube voltage change of the high energy image HP can be adjusted using the ratio GpH / GpL of the Gp values of the high energy image HP.
ところで、低エネルギー画像LPにおいても高エネルギー画像HPにおいて散乱線の影響が生じるのと同様に散乱線の影響が生じる。具体的には、低エネルギー画像LPも管電圧が高くなるほどヒストグラム幅が狭くなり、GpL値が大きくなる。つまり、管電圧が大きくなると低エネルギー画像LPのコントラストが小さくなるため、その分高エネルギー画像HPのコントラスト(特に低濃度領域LCH)を小さく調整する必要がある。このとき、低エネルギー画像LPの管電圧の変化についても上述したGp値の比GpH/GpLで調整することが考えられる。 By the way, in the low energy image LP, the influence of the scattered radiation is caused in the same manner as the influence of the scattered radiation in the high energy image HP. Specifically, the histogram width of the low energy image LP becomes narrower and the GpL value becomes larger as the tube voltage becomes higher. That is, as the tube voltage increases, the contrast of the low energy image LP decreases, and accordingly, the contrast of the high energy image HP (particularly, the low density region LCH) needs to be adjusted to be small. At this time, it is conceivable to adjust the change in the tube voltage of the low energy image LP with the Gp value ratio GpH / GpL described above.
しかし、低エネルギー画像LPの管電圧の変化をGp値の比GpH/GpLで調整した場合、軟部画像に骨が残ってしまって軟部と骨部との分離が十分ではなく、さらに低エネルギー画像LPの管電圧が大きくなればなる程その傾向が顕著に表れることがわかった。 However, when the change in the tube voltage of the low energy image LP is adjusted by the ratio GpH / GpL of the Gp value, the bone remains in the soft part image, and the soft part and the bone part are not sufficiently separated, and the low energy image LP It has been found that this tendency becomes more pronounced as the tube voltage increases.
ここで、図16は低エネルギー画像LPを取得する際の管電圧を60〜80kVpに変化させて同一被写体を撮影した場合のエネルギースペクトルを示すグラフである。図16に示すように、低エネルギー画像LPの管電圧を上昇させた場合、50keV以上のエネルギー成分が含まれる割合は、60kVpで9%, 70kVpで22%,80kVpでは32%になる。これは、上述のビームハードニングの要因である光電効果およびコンプトン効果のうち、低エネルギー画像LPの管電圧上昇に応じて光電効果の影響よりもコンプトン散乱の影響がより強くなるためである。言い換えれば、低エネルギー画像LPの取得時に用いる管電圧値は、光電効果とコンプトン効果との双方が混在し管電圧値によって双方の影響度が異なるため(図3参照)、GpL値を用いただけでは適正な補正を行うことができないことがわかった。 Here, FIG. 16 is a graph showing an energy spectrum when the same subject is photographed while changing the tube voltage when acquiring the low energy image LP to 60 to 80 kVp. As shown in FIG. 16, when the tube voltage of the low energy image LP is increased, the ratio of energy components of 50 keV or more is 9% at 60 kVp, 22% at 70 kVp, and 32% at 80 kVp. This is because, among the photoelectric effect and the Compton effect, which are the factors of the beam hardening described above, the influence of Compton scattering becomes stronger than the influence of the photoelectric effect in accordance with the increase in the tube voltage of the low energy image LP. In other words, the tube voltage value used when acquiring the low-energy image LP is a mixture of both the photoelectric effect and the Compton effect, and the degree of influence differs depending on the tube voltage value (see FIG. 3). It was found that proper correction cannot be performed.
そこで、低エネルギー画像LPの管電圧変化による係数の調整をGp値の比GpH/GpLによる調整ではなく(図13参照)、別途複数の第1濃度変換テーブルHPT11〜HPT13、第2濃度変換テーブルLPT11〜LPT13から各管電圧に対応した適切な濃度変換テーブルHPT11〜LPT13、LPT11〜LPT13を選択して階調調整が行われる。これにより、確実にサブトラクション画像の画質劣化を低減することができる。 Therefore, the adjustment of the coefficient due to the change in the tube voltage of the low energy image LP is not the adjustment by the Gp value ratio GpH / GpL (see FIG. 13), but a plurality of first density conversion tables HPT11 to HPT13 and a second density conversion table LPT11. The gradation adjustment is performed by selecting appropriate density conversion tables HPT11 to LPT13 and LPT11 to LPT13 corresponding to each tube voltage from .about.LPT13. Thereby, it is possible to reliably reduce the image quality deterioration of the subtraction image.
なお、従来、高エネルギー画像HPと低エネルギー画像LPの管電圧に応じてパラメータ(画像内で一様のスカラー量)を、ルックアップテーブルを用いて選択しているが、胸部正面撮影において高エネルギー画像HPは診断画像であり、そのコントラストが大きく変動する要因は、管電圧よりもむしろ個人の体格変動によるものであると考えられ、上記図12〜16とは異なるものである。 Conventionally, a parameter (a uniform scalar amount in the image) is selected using a look-up table according to the tube voltage of the high energy image HP and the low energy image LP. The image HP is a diagnostic image, and it is considered that the factor that the contrast greatly fluctuates is due to the individual physique fluctuation rather than the tube voltage, and is different from the above-described FIGS.
また、図12において厚み検出手段140が設けられている場合について例示しているが、厚みに基づいた調整を行うことなく、低エネルギー画像LPの取得時の管電圧により各濃度変換テーブルHPT11〜HPT13、LPT11〜HPT13を選択して、第1の実施形態に示すような階調調整を行うようにしてもよい(図1参照)。
12 illustrates the case where the
上記各実施の形態によれば、放射線源から異なるエネルギーの放射線を被写体に照射し撮影したときの高エネルギー画像HPと低エネルギー画像とを取得し、取得した高エネルギー画像HPの低濃度領域LCHに対し濃度勾配を大きくする、および/または低エネルギー画像の低濃度領域LCLに対し濃度勾配を小さくする階調処理を行い、調整した高エネルギー画像HPと低エネルギー画像LPとを用いてエネルギーサブトラクション処理を行うことにより、被写体の局所的な厚みの相違によりビームハードニングの影響度が異なる場合であっても、エネルギーサブトラクション処理の前処理としてビームハードニングの影響を考慮した階調調整を行うことができるため、ビームハードニングの影響によるサブトラクション画像の画質の劣化を低減することができる。 According to each of the embodiments described above, a high energy image HP and a low energy image obtained by irradiating a subject with radiation of different energy from a radiation source are acquired, and the acquired low energy region LCH of the high energy image HP is acquired. On the other hand, gradation processing for increasing the density gradient and / or reducing the density gradient for the low density region LCL of the low energy image is performed, and energy subtraction processing is performed using the adjusted high energy image HP and the low energy image LP. As a result, even when the influence of beam hardening differs due to a difference in the local thickness of the subject, gradation adjustment can be performed in consideration of the influence of beam hardening as pre-processing of the energy subtraction process. Therefore, the image quality of the subtraction image is inferior due to the effect of beam hardening. It is possible to reduce the.
また、図9に示すように、被写体の厚みを検出する厚み検出手段140をさらに備え、階調処理手段115が、厚み検出手段140により検出された被写体の厚みが大きくなればなるほど、高エネルギー画像HPの低濃度領域LCHに対する濃度勾配をより大きくするおよび/または低エネルギー画像LPの低濃度領域LCLに対する濃度勾配をより小さくする調整を行うものであるとき、被写体の厚み(体格)によりビームハードニングの影響度が被写体毎に異なる場合であっても、ビームハードニングの影響度に応じた前処理を施して安定したエネルギーサブトラクション結果を得ることができる。
In addition, as shown in FIG. 9, the image processing apparatus further includes a
さらに、図12に示すように、低エネルギー画像LPの取得時に放射線源に印加される管電圧を検出する管電圧検出手段240をさらに有し、階調処理手段215が、低エネルギー画像LPの取得時の管電圧毎に、管電圧が大きくなればなるほど濃度勾配が大きい複数の第1濃度変換テーブルHPT11〜13および/または管電圧が大きくなればなるほど濃度勾配が小さい複数の第2濃度変換テーブルLPT11〜13を有するものであり、管電圧検出手段240により検出された管電圧に基づいて、複数の第1濃度変換テーブルHPT11〜13および/または複数の第2濃度変換テーブルLPT11〜13のうちいずれかの第1濃度変換テーブルHPT11〜13および/または第2濃度変換テーブルLPT11〜13を選択し、選択した第1濃度変換テーブルHPT11〜13および/または第2濃度変換テーブルLPT11〜13を用いて濃度値を調整するものであれば、低エネルギー画像LPの取得時における管電圧(放射線エネルギー)の変化によりビームハードニングの影響度が撮影毎に異なる場合であって信号分布の幅の比では適正な調整を行うことができない場合であっても、ビームハードニングの影響度に応じた前処理を施して安定したエネルギーサブトラクション結果が得ることができる。 Furthermore, as shown in FIG. 12, it further has a tube voltage detecting means 240 for detecting a tube voltage applied to the radiation source when acquiring the low energy image LP, and the gradation processing means 215 acquires the low energy image LP. For each tube voltage, a plurality of first concentration conversion tables HPT11 to 13 having a larger concentration gradient as the tube voltage increases and / or a plurality of second concentration conversion tables LPT11 having a smaller concentration gradient as the tube voltage increases. To 13 and based on the tube voltage detected by the tube voltage detecting means 240, one of the plurality of first concentration conversion tables HPT11 to 13 and / or the plurality of second concentration conversion tables LPT11 to LPT13. The first density conversion tables HPT 11 to 13 and / or the second density conversion tables LPT 11 to 13 are selected and selected. If the density value is adjusted using the first density conversion tables HPT11 to 13 and / or the second density conversion tables LPT11 to LP13, the change in tube voltage (radiation energy) at the time of acquiring the low energy image LP Even if the degree of influence of beam hardening varies from one shot to the next, and even when the ratio of the signal distribution width cannot be adjusted properly, pre-processing is performed according to the degree of influence of beam hardening. Stable energy subtraction results can be obtained.
本発明の実施形態は、上記実施形態に限定されない。たとえば、FPDを用いた撮影コンソールに搭載される画像処理について例示しているが、撮影システムにネットワークを介して接続された画像処理ワークステーションや、画像保管通信システム(PACS:Picture Archiving and Communication System)などに転送された画像データに対してのエネルギーサブトラクション処理にも適用することが可能である。 The embodiment of the present invention is not limited to the above embodiment. For example, the image processing installed in the imaging console using FPD is illustrated, but the image processing workstation connected to the imaging system via the network and the image archiving communication system (PACS: Picture Archiving and Communication System) It is also possible to apply to energy subtraction processing for image data transferred to the above.
また、撮影モダリティもFPDに限定されず、CRシステムを用いた2ショットエネサブにも適用可能である。また、撮影コンソールにCRとFPDが同時に接続されているような、CR-DR混在環境の下である場合、CR(1ショット法)用のエネサブ係数パラメータセットと、DR(2ショット)用のエネサブ係数パラメータセットを撮影されたモダリティに応じて切り替えてもよい。 In addition, the imaging modality is not limited to FPD, and can be applied to a two-shot energy sub using a CR system. Also, in a CR-DR mixed environment where CR and FPD are connected to the shooting console at the same time, the energy sub coefficient parameter set for CR (1 shot method) and the energy sub for DR (2 shot) The coefficient parameter set may be switched according to the captured modality.
さらに、サブトラクション処理により生成された軟部画像SPの平滑化を行い、ノイズが抑制された軟部画像SPを高エネルギー画像HPから減算することで、ノイズが抑制された骨部画像BPを得ることもできる。また、ノイズが抑制された骨部画像BPを高エネルギー画像HPもしくは高エネルギー画像HPと低エネルギー画像LPとを加算した加算画像から減算することにより、ノイズが抑制された軟部画像SPを得ることができ、ノイズ抑制処理と診断画像からの減算処理を繰り返すことで、段階的にノイズ抑制することも可能である。 Furthermore, by performing smoothing of the soft part image SP generated by the subtraction process and subtracting the soft part image SP in which noise is suppressed from the high energy image HP, a bone part image BP in which noise is suppressed can be obtained. . Further, by subtracting the bone image BP in which noise is suppressed from the high energy image HP or an addition image obtained by adding the high energy image HP and the low energy image LP, a soft part image SP in which noise is suppressed can be obtained. It is also possible to suppress noise in stages by repeating the noise suppression process and the subtraction process from the diagnostic image.
また、上記実施の形態において、サブトラクション処理手段30は胸部撮影したときの軟部画像SPと骨部画像BPとをサブトラクション画像Psubとして生成する場合について例示しているが、マンモグラフィ装置により撮影した人体の乳房を被写体とした場合の乳腺が強調された画像と悪性腫瘍が強調された画像をサブトラクション画像Psubとして生成するものであってもよい。
Further, in the above embodiment, the
1、100、200 エネルギーサブトラクション処理装置
10 画像取得手段
15、115、215 階調処理手段
20 位置合わせ手段
30 サブトラクション処理手段
40 画像調整手段
140 厚み検出手段
BP 骨部画像
SP 軟部画像
GpH 高エネルギー画像のスケールファクター値
GpL 低エネルギー画像のスケールファクター値
HP 高エネルギー画像
LP 低エネルギー画像
Ka 第1荷重減算係数
Kb 第2荷重減算係数
HPT1、HPT11〜HPT13 第1濃度変換テーブル
LPT11〜LPT13 第2濃度変換テーブル
1, 100, 200 Energy
Claims (6)
該画像取得手段において取得された前記高エネルギー画像および/または前記低エネルギー画像に対し濃度値の調整を行う階調処理手段と、
該階調処理手段により調整された前記高エネルギー画像および/または前記低エネルギー画像を用いてエネルギーサブトラクション処理を行うサブトラクション処理手段と
を備え、
前記階調処理手段が、前記高エネルギー画像の低濃度領域に対し濃度勾配を大きくするおよび/または前記低エネルギー画像の低濃度領域に対し濃度勾配を小さくするような濃度値の調整を行うものであることを特徴とするエネルギーサブトラクション処理装置。 Image acquisition means for acquiring a high energy image and a low energy image when a subject is irradiated with radiation of different energy from a radiation source, and
Gradation processing means for adjusting a density value for the high energy image and / or the low energy image acquired by the image acquisition means;
Subtraction processing means for performing energy subtraction processing using the high energy image and / or the low energy image adjusted by the gradation processing means,
The gradation processing means adjusts the density value so as to increase the density gradient with respect to the low density area of the high energy image and / or reduce the density gradient with respect to the low density area of the low energy image. An energy subtraction processing device characterized by being.
前記階調処理手段が、前記厚み検出手段により検出された前記被写体の厚みが大きくなればなるほど、前記高エネルギー画像に対する濃度勾配をより大きくするおよび/または前記低エネルギー画像に対する濃度勾配をより小さくする調整を行うものであることを特徴とする請求項1記載のエネルギーサブトラクション処理装置。 A thickness detecting means for detecting the thickness of the subject;
The gradation processing means increases the density gradient for the high energy image and / or reduces the density gradient for the low energy image as the thickness of the subject detected by the thickness detection means increases. The energy subtraction processing apparatus according to claim 1, wherein adjustment is performed.
前記階調処理手段が、前記低エネルギー画像の取得時の管電圧毎に、該管電圧が大きくなればなるほど前記濃度勾配が大きく設定されている複数の前記第1濃度変換テーブルおよび/または前記管電圧が大きくなればなるほど前記濃度勾配が小さく設定されている複数の前記第2濃度変換テーブルを有するものであり、前記管電圧検出手段により検出された前記管電圧に基づいて、前記複数の第1濃度変換テーブルおよび/または前記複数の第2濃度変換テーブルのうちいずれかの前記第1濃度変換テーブルおよび/または前記第2濃度変換テーブルを選択し、選択した前記第1濃度変換テーブルおよび/または前記第2濃度変換テーブルを用いて濃度値の調整を行うものであることを特徴とする請求項1記載のエネルギーサブトラクション処理装置。 Tube voltage detecting means for detecting a tube voltage applied to the radiation source at the time of acquiring the low energy image;
For each tube voltage at the time of acquisition of the low energy image, the gradation processing means has a plurality of the first density conversion tables and / or the tubes in which the density gradient is set to be larger as the tube voltage becomes higher. The plurality of second concentration conversion tables in which the concentration gradient is set to be smaller as the voltage is increased, and the plurality of first voltages are based on the tube voltage detected by the tube voltage detecting means. One of the first density conversion table and / or the second density conversion table is selected from the density conversion table and / or the plurality of second density conversion tables, and the selected first density conversion table and / or the selected density conversion table is selected. 2. The energy subtraction according to claim 1, wherein the density value is adjusted using the second density conversion table. Management apparatus.
取得した前記高エネルギー画像の低濃度領域に対し濃度勾配を大きくするおよび/または前記低エネルギー画像の低濃度領域に対し濃度勾配を小さくする階調処理を行い、
調整した前記高エネルギー画像と前記低エネルギー画像とを用いてエネルギーサブトラクション処理を行う
ことを特徴とするエネルギーサブトラクション処理方法。 Acquire a high energy image and a low energy image when the subject is irradiated with radiation of different energy from the radiation source,
Performing gradation processing to increase the density gradient for the low density region of the acquired high energy image and / or reduce the density gradient to the low density region of the low energy image;
An energy subtraction processing method, wherein energy subtraction processing is performed using the adjusted high energy image and the low energy image.
放射線源から異なるエネルギーの放射線を被写体に照射し撮影したときの高エネルギー画像と低エネルギー画像とを取得し、
取得した前記高エネルギー画像の低濃度領域に対し濃度勾配を大きくする、および/または前記低エネルギー画像の低濃度領域に対し濃度勾配を小さくする階調処理を行い、
調整した前記高エネルギー画像と前記低エネルギー画像とを用いてエネルギーサブトラクション処理を行う
ことを実行させるためのエネルギーサブトラクション処理プログラム。 On the computer,
Acquire a high energy image and a low energy image when the subject is irradiated with radiation of different energy from the radiation source,
Performing gradation processing to increase the density gradient with respect to the low density area of the acquired high energy image and / or to reduce the density gradient with respect to the low density area of the low energy image;
An energy subtraction processing program for executing energy subtraction processing using the adjusted high energy image and the low energy image.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008145331A JP5190942B2 (en) | 2008-06-03 | 2008-06-03 | Energy subtraction processing apparatus and method, and program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008145331A JP5190942B2 (en) | 2008-06-03 | 2008-06-03 | Energy subtraction processing apparatus and method, and program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2009291271A JP2009291271A (en) | 2009-12-17 |
JP5190942B2 true JP5190942B2 (en) | 2013-04-24 |
Family
ID=41540039
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008145331A Active JP5190942B2 (en) | 2008-06-03 | 2008-06-03 | Energy subtraction processing apparatus and method, and program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5190942B2 (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6458540B2 (en) * | 2015-02-23 | 2019-01-30 | コニカミノルタ株式会社 | Radiographic imaging system, radiographic imaging apparatus, and body thickness estimation method |
KR101848177B1 (en) * | 2016-12-07 | 2018-04-11 | 건양대학교산학협력단 | Method and system for weighted-image subtraction based on material characteristics for high-contrast decomposition |
JP6906479B2 (en) * | 2018-05-25 | 2021-07-21 | 富士フイルム株式会社 | Bone mineral information acquisition device, method and program |
JP7221981B2 (en) * | 2018-09-18 | 2023-02-14 | 富士フイルム株式会社 | Image processing device, image processing method, and image processing program |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6244224A (en) * | 1985-08-23 | 1987-02-26 | 富士写真フイルム株式会社 | Image processing method and apparatus |
JP2981693B2 (en) * | 1992-06-30 | 1999-11-22 | 富士写真フイルム株式会社 | Energy subtraction image generation method |
JP3165538B2 (en) * | 1992-12-10 | 2001-05-14 | 富士写真フイルム株式会社 | Radiation image alignment method and apparatus |
JP2000121579A (en) * | 1998-10-13 | 2000-04-28 | Nagoya Denki Kogyo Kk | Method and apparatus for x-ray image pickup as well as medium with recorded x-ray image pickup control program |
JP2002143137A (en) * | 2000-11-15 | 2002-05-21 | Hitachi Medical Corp | X-ray radiography apparatus |
JP4245297B2 (en) * | 2002-02-18 | 2009-03-25 | 東芝医用システムエンジニアリング株式会社 | X-ray diagnostic imaging equipment |
JP2007236766A (en) * | 2006-03-10 | 2007-09-20 | Shimadzu Corp | X-ray photographing apparatus |
-
2008
- 2008-06-03 JP JP2008145331A patent/JP5190942B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2009291271A (en) | 2009-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4549513B2 (en) | Imaging system including radiation filter for X-ray imaging | |
Graham et al. | Compensators for dose and scatter management in cone‐beam computed tomography | |
US8417048B2 (en) | Method for correction of digital images | |
US9345443B2 (en) | Calibration free dual energy radiography method | |
US8165379B2 (en) | Method of processing radiological images, and, in particular, mammographic images | |
JP2024045424A (en) | Method and system for determining virtual outputs for multi-energy x-ray imaging apparatus | |
JP6678541B2 (en) | Image processing apparatus, method and program | |
Crotty et al. | Experimental spectral measurements of heavy K-edge filtered beams for x-ray computed mammotomography | |
WO2019026409A1 (en) | Radiation image pickup device and radiation image pickup system | |
JP5190942B2 (en) | Energy subtraction processing apparatus and method, and program | |
Baldelli et al. | A prototype of a quasi-monochromatic system for mammography applications | |
JP6083990B2 (en) | Radiation imaging apparatus, control method thereof, and program | |
JP2008167948A (en) | Radiographic image processing method and apparatus, and program | |
US8121372B2 (en) | Method for reducing image noise in the context of capturing an image using two different radiation spectra | |
Stecke et al. | Alternative X-ray filters for an intra-oral digital radiographic system | |
JP5210048B2 (en) | Energy subtraction processing apparatus and method, and program | |
JP5534699B2 (en) | X-ray diagnostic apparatus and image processing apparatus | |
Ducote et al. | Dual-energy cardiac imaging: an image quality and dose comparison for a flat-panel detector and x-ray image intensifier | |
JP5044457B2 (en) | Radiation imaging apparatus and imaging condition setting method | |
WO2015045379A1 (en) | Image processing device and method | |
JP2007268033A (en) | Radiography system and radiography method | |
JP7531337B2 (en) | Radiography system, control device, control method, and program | |
Nickoloff et al. | Pediatric high KV/filtered airway radiographs: comparison of CR and film-screen systems | |
Wong et al. | Effect of area x-ray beam equalization on image quality and dose in digital mammography | |
US20100061654A1 (en) | Scatter estimation and reduction method and apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20110125 |
|
RD15 | Notification of revocation of power of sub attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7435 Effective date: 20110208 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20120814 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120821 |
|
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: 20130115 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20130124 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5190942 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: 20160208 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 |