JP6132938B2 - Image processing apparatus, radiation imaging system, image processing method, and recording medium - Google Patents
Image processing apparatus, radiation imaging system, image processing method, and recording medium Download PDFInfo
- Publication number
- JP6132938B2 JP6132938B2 JP2016002063A JP2016002063A JP6132938B2 JP 6132938 B2 JP6132938 B2 JP 6132938B2 JP 2016002063 A JP2016002063 A JP 2016002063A JP 2016002063 A JP2016002063 A JP 2016002063A JP 6132938 B2 JP6132938 B2 JP 6132938B2
- Authority
- JP
- Japan
- Prior art keywords
- range
- pixel values
- image
- gradation conversion
- image processing
- 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 105
- 238000003384 imaging method Methods 0.000 title claims description 23
- 230000005855 radiation Effects 0.000 title claims description 11
- 238000003672 processing method Methods 0.000 title claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 108
- 238000012937 correction Methods 0.000 claims description 37
- 238000000034 method Methods 0.000 claims description 28
- 230000008569 process Effects 0.000 claims description 21
- 229920006395 saturated elastomer Polymers 0.000 claims description 5
- 230000006870 function Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 11
- 238000010191 image analysis Methods 0.000 description 10
- 230000002708 enhancing effect Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 4
- FTGYKWAHGPIJIT-UHFFFAOYSA-N hydron;1-[2-[(2-hydroxy-3-phenoxypropyl)-methylamino]ethyl-methylamino]-3-phenoxypropan-2-ol;dichloride Chemical compound Cl.Cl.C=1C=CC=CC=1OCC(O)CN(C)CCN(C)CC(O)COC1=CC=CC=C1 FTGYKWAHGPIJIT-UHFFFAOYSA-N 0.000 description 3
- 239000013307 optical fiber Substances 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000002601 radiography Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Images
Description
本発明は、イメージセンサにより撮像された画像に対する画像処理に関する。 The present invention relates to image processing on an image captured by an image sensor.
通常の撮影システムでは、イメージセンサが受ける光の量を調節したり、受光量に対する出力の特性を調整したりすることで、画像の黒つぶれや白とびを低減する。また出力された画像に処理を施すことで画像の画質の水準を保っている。X線撮影システムを例にすると、線量制御機構がX線源を制御することでイメージセンサに到達する線量が調節される。到達した線量に対する出力値をイメージセンサのゲインを変更することにより調整する。調整された出力値に応じて得られるX線画像に階調数の圧縮・拡大処理や、画像の特定範囲を強調するコントラスト強調処理などの画像処理が施される。 In a normal photographing system, the amount of light received by an image sensor is adjusted, or the characteristics of output with respect to the amount of received light are adjusted, thereby reducing blackout and overexposure of an image. In addition, the level of image quality is maintained by processing the output image. Taking an X-ray imaging system as an example, the dose reaching the image sensor is adjusted by the dose control mechanism controlling the X-ray source. The output value for the reached dose is adjusted by changing the gain of the image sensor. The X-ray image obtained according to the adjusted output value is subjected to image processing such as gradation compression / enlargement processing and contrast enhancement processing for enhancing a specific range of the image.
線量や入出力の特性に応じて画像処理を行う技術として、動画撮影において線量が頻繁に変更されている場合、線量が安定するまで画像処理の応答を弱め、画素値のふらつきによる影響を抑える技術がある(特許文献1)。また、被写体の輝度範囲に合わせてセンサのゲインを変更してセンサの感度領域内に収まるようにし、ゲイン毎に設定された関数で階調変換する方法が開示されている(特許文献2)。 Technology that performs image processing according to the dose and input / output characteristics. When the dose is frequently changed during movie shooting, this technology reduces the response of the image processing until the dose stabilizes and suppresses the effects of pixel value fluctuations. (Patent Document 1). In addition, a method is disclosed in which the gain of the sensor is changed in accordance with the luminance range of the subject so as to be within the sensitivity region of the sensor, and gradation conversion is performed using a function set for each gain (Patent Document 2).
しかしながら、出力特性が良好な範囲から大きく外れた画像が得られた場合に、ノイズを強調した画像が出力されてしまうという問題がある。例えば受光量が不足した画像では光量不足に起因するノイズが多く重畳することとなる。この場合に適切な受光量の画像と同様にコントラスト強調処理をしてしまうと、このノイズを強調した画像が出力されてしまう。 However, there is a problem in that an image with enhanced noise is output when an image whose output characteristics deviate greatly from a good range is obtained. For example, in an image where the amount of received light is insufficient, a large amount of noise due to insufficient light amount is superimposed. In this case, if contrast enhancement processing is performed in the same manner as an image having an appropriate amount of received light, an image in which this noise is enhanced is output.
そこで本発明は、階調変換の特性を取得する取得手段と、X線センサのダイナミックレンジに対応する画素値の範囲を記憶する記憶手段と、前記階調変換の特性により階調が設定される画素値の範囲が、前記記憶された画素値の範囲に含まれるか否かを判定する判定手段と、前記判定の結果に応じて前記階調変換の特性を補正する補正手段と、を有することを特徴とする。 Therefore, the present invention sets the gradation by the acquisition means for acquiring the characteristics of the gradation conversion , the storage means for storing the pixel value range corresponding to the dynamic range of the X-ray sensor, and the characteristics of the gradation conversion. Determining means for determining whether a range of pixel values is included in the stored range of pixel values; and correcting means for correcting the characteristics of the gradation conversion in accordance with the result of the determination. It is characterized by.
かかる本発明によれば、階調変換の特性により階調が設定される画素値の範囲がX線センサのダイナミックレンジに対応する画素値の範囲に応じて階調変換の特性を補正するため、ノイズの発生を抑えた画像を生成することができる。 According to the present invention, since the range of pixel values in which gradation is set by the characteristics of gradation conversion corrects the characteristics of gradation conversion according to the range of pixel values corresponding to the dynamic range of the X-ray sensor , An image with reduced noise generation can be generated.
以下、適宜図面を参照しながら本発明を実施するための形態を説明する。 DESCRIPTION OF EMBODIMENTS Hereinafter, embodiments for carrying out the present invention will be described with reference to the drawings as appropriate.
本発明を放射線撮影システムで実施する例を図1乃至図5を参照しながら説明する。 An example in which the present invention is implemented by a radiation imaging system will be described with reference to FIGS.
図1は本実施例の画像処理装置100の構成を示す。画像処理装置100は画像入力部101、設定部102、特性取得部103、判定部104、補正部105、階調処理部106、表示制御部107を有する。画像入力部101は、イメージセンサで取得された画像に対してセンサ特性が補正された画像を取得する。この画像処理装置100による処理がされる前の画像をオリジナル画像と呼ぶことにする。
FIG. 1 shows the configuration of an image processing apparatus 100 of this embodiment. The image processing apparatus 100 includes an
設定部102は入力された画像についてコントラストを強調する画素値の範囲を設定する。本実施例では設定部102が有する画像解析部108と階調変換作成部109によりこの処理が実現される。画像解析部108は、オリジナル画像について階調処理に必要な情報をオリジナル画像から計算し解析結果を出力する。階調変換作成部109は、解析結果に基づいてオリジナル画像を適切な画像にするための階調変換曲線を作成する。階調変換曲線の傾きによって強調される画素値の範囲が設定される。なお階調変換作成部109では階調変換曲線に限らず入力画素値に対する出力画素値を示すLUTを作成してもよい。さらには画像において強調する領域を設定すればよく、階調変換を作成しなくてもよい。
The
特性取得部103はオリジナル画像を撮像する際に用いられたイメージセンサのダイナミックレンジを示す値を記憶部から取得する。ここでダイナミックレンジとは、イメージセンサの線量に対する出力特性が良好な線量の範囲またはこれに対応する画素値の範囲をいう。
The
判定部104では、このダイナミックレンジに、設定部102にて設定された強調される画素値の範囲が包含されるか否かを判定する。包含されないと判定された場合には、補正部105は階調変換曲線を入力とし、ダイナミックレンジに包含されるように画像において強調される画素値の範囲を補正する。補正済みの階調変換曲線を出力する。一方、判定部104にて包含されると判定された場合には階調変換曲線は補正されない。階調処理部106は補正されなかった階調変換曲線または補正済みの階調変換曲線とオリジナル画像とを入力とし、階調処理した画像を出力する。表示制御部107は階調処理後の画像を表示部へと出力する。
The
上述の画像処理装置100では回路で実装した場合を示したが、画像処理装置100をソフトウェアで実装することもできる。その場合、先述の各部はソフトウェアとハードウェアにより実現される。 In the above-described image processing apparatus 100, the case where the image processing apparatus 100 is implemented by a circuit has been shown. In that case, the above-described units are realized by software and hardware.
図2には画像処理装置100をハードウェアとソフトウェアにより実装する場合のハードウェア構成を示している。画像処理装置100であるコントロールPC201はCPU203、RAM204、ROM205、記憶部206を有している。これらはバス209を介して接続されている。さらに外部装置として表示部207、入力部208が接続される。後述する図3に示す処理をコントロールPC201に実行させるプログラムがソフトウェアモジュールとしてROM205または記憶部206に格納されている。このプログラムはRAM204に読み込まれ、CPU203に実行される。CPU203はプログラムされた命令に応じてコントロールPC201の各ハードウェア資源を制御し、全体として図1に示す各機能が実現される。かかるプログラムは撮影モード毎に分かれたソフトウェアモジュールとして格納されている。
FIG. 2 shows a hardware configuration when the image processing apparatus 100 is implemented by hardware and software. A
図2はまた、画像処理装置100を含む放射線撮影システム200のハードウェア構成を示している。画像処理装置100であるコントロールPC201とイメージセンサ202、X線発生装置211が信号線である光ファイバー214でつながっている。信号線は光ファイバーでなくてもCAN(Controller Area Network)やギガビットイーサなどでもよい。光ファイバー214には、表示部210、データベース212、ネットワークインタフェース部213が接続されている。
FIG. 2 also shows a hardware configuration of the radiation imaging system 200 including the image processing apparatus 100. The
コントロールPC201はイメージセンサ202や放射線源であるX線発生装置211などにコマンドを送る。イメージセンサ202はX線を光に変換する蛍光体と、蛍光体が発する光を電気信号に変換する半導体基板と、電気信号を読み出すための駆動回路と、読み出された電気信号からデジタル画像データを作成する回路とを含む。表示部210はX線発生装置211や撮影に応じて得られる動画像をユーザがリアルタイムで確認するためのGUIの表示に主に用いられる。一方表示部207はユーザが撮影条件の設定を行うためのGUIの表示に主に用いられる。
The
本システムでは動画撮影と静止画撮影が可能となっている。動画撮影の場合、イメージセンサ202が受光して得られる出力値に基づいて、放射線の線量が制御される。これによりイメージセンサ202が出力するX線画像の一定画質を保持することができる。
In this system, video shooting and still image shooting are possible. In the case of moving image shooting, the radiation dose is controlled based on an output value obtained by receiving light by the
この制御においてはCPU203が出力値を定期的にチェックすることでイメージセンサ202の決められた範囲に到達する線量を監視し、その値が一定になるように入射線量を制御する。CPU203は到達線量が基準より少ないと判定した場合に入射線量を多くするようにX線発生装置211に命令する。到達線量が基準を超えていると判定すれば入射線量を少なくするようにX線発生装置211に対して命令する。X線発生装置211は命令に応じて適切な線量のX線をX線管から射出する。また、表示部210のGUIや不図示の操作ボタン、または入力部208を介して手動での線量の制御が可能である。
In this control, the
このように透視撮影やシネ撮影などの動画撮影では線量制御は自動で行われるため、線量の監視に用いる範囲が適切でないと必ずしも最適な線量になるとは限らない。また動画撮影を開始した直後や撮影領域が移動する場合は、フィードバック制御が収束するために時間がかかることから、その間イメージセンサからは適切な線量でない画像が出力される恐れがある。また、線量制御は正しく動作する場合でも、安全性の観点から線量に上限があるため、必要な量だけの線量を曝射することができない場合もある。この場合もノイズだけを強調してしまった画像を出力することになる。 Thus, in moving image photography such as fluoroscopic photography and cine photography, dose control is automatically performed. Therefore, if the range used for dose monitoring is not appropriate, the optimum dose is not always obtained. In addition, immediately after the start of moving image shooting or when the shooting area moves, it takes time for feedback control to converge, and during that time, an image with an appropriate dose may be output from the image sensor. Even when dose control operates correctly, there is an upper limit on the dose from the viewpoint of safety, and there is a case where only a necessary dose cannot be exposed. In this case, an image in which only noise is emphasized is output.
画像処理装置100により行われる処理の概要を図3に基づいて説明する。図3は、イメージセンサ202に検出される線量に応じた画素値と、出力特性が良好なダイナミックレンジの範囲と、イメージセンサ202が生成した画像の画素値の範囲との関係を示す図である。
An outline of processing performed by the image processing apparatus 100 will be described with reference to FIG. FIG. 3 is a diagram illustrating a relationship between a pixel value corresponding to a dose detected by the
図3(a)で、出力特性が良好な範囲を示すダイナミックレンジ300が示されている。これに対して、線量(受光量)に応じた画像が生成されることとなる。範囲301はある画像の画素値の範囲を示している。ここでは、画像の画素値のヒストグラムを取ったときに出現回数が所定の閾値よりも高い値の下限と上限で画される範囲を示している。設定範囲302は、設定部102が画像の関心領域の画素値に基づいて設定した画像において強調される画素値の範囲である。
FIG. 3A shows a
範囲301を画素値の範囲とする画像では、設定範囲302がダイナミックレンジ300に一部包含されない。この場合、補正部105は包含されない範囲分だけ設定範囲302を高画素値側にシフトする。ここでシフトとは、画素値の範囲の大きさを変えずに下限値と上限値を変えることを意味する。範囲303はこのシフト補正された範囲を示している。これにより、画像の本来強調されるべき領域が強調されず潰れたような画像となる。
よって作業者に対して、この画像を見て画像の線量が不足していることを直感的に理解させることができる。
In the image having the range 301 as the pixel value range, the
Therefore, the operator can intuitively understand that the dose of the image is insufficient by viewing this image.
また別の画像の画素値の範囲が範囲311で示されている。この画像において強調される画素値の範囲が312で示されている。この場合、設定範囲312はダイナミックレンジ300の範囲を超えている部分があるため、補正部105が強調する範囲を範囲313にシフトする。これにより、画像の本来強調されるべき領域よりも低画素値側が強調されるため、白とびしたような画像を得ることができる。よって作業者に対して、この画像を見て画像の線量が不足していることを直感的に理解させることができる。
A range of pixel values of another image is indicated by a range 311. A range of pixel values to be emphasized in this image is indicated by 312. In this case, since the
図3(b)では、ある画像の画素値の範囲が範囲331で示されており、その強調されるべき範囲が設定範囲332で示されている。この場合、画像の画素値の範囲も、設定範囲332もダイナミックレンジ300内に包含されているため、補正部105による補正は必要ない。本来強調されるべき設定範囲332が強調された画像を得ることができる。
In FIG. 3B, the range of pixel values of an image is indicated by a range 331, and the range to be emphasized is indicated by a
また別の画像における画素値の範囲が範囲341で示され、強調される範囲が設定範囲342で示されている。この場合、画像の画素値の範囲341は一部ダイナミックレンジ300に包含されていないが、設定範囲342は包含されている。この場合にも、ノイズやアーチファクトを強調することがないため、設定範囲342を強調した画像を得ることができる。
Further, a range of pixel values in another image is indicated by a range 341, and a range to be emphasized is indicated by a
上述の放射線撮影システム200により実現される処理の流れを図4に従い説明する。
図4は放射線撮影システム200が実行する動画撮影処理の流れを示している。
The flow of processing realized by the radiation imaging system 200 will be described with reference to FIG.
FIG. 4 shows the flow of the moving image capturing process executed by the radiation imaging system 200.
まずステップS401の前段階として、X線発生装置211はX線管を制御し断続的にX線を射出させる。ここで射出するX線の線量制御は、イメージセンサ202の所定領域が受光して得られる出力値を所定のタイミングで取得することで監視し、それに応じてX線発生装置211により制御される。
First, as a step before step S401, the
X線の射出と同期して、イメージセンサ202はX線を受光しフレーム画像を連続して出力する。この画像に対してイメージセンサ202内の回路でセンサ特性の補正がされる。ここでセンサ特性の補正とは、オフセット補正、Log変換、ゲイン補正、欠陥画素補正などである。この補正された画像の生成と画像処理装置100の処理が並行して行われる。
In synchronization with the emission of X-rays, the
ステップS401では、画像処理に必要な情報を外部から取得する。画像入力部101は、センサ特性を補正したオリジナル画像を取得する。また特性取得部103は、予め計測され記憶部206に格納された、イメージセンサ202のダイナミックレンジを示す値を取得する。
In step S401, information necessary for image processing is acquired from the outside. The
以下ステップS402乃至S404で、設定部102が画像に基づいて強調する画素値の範囲を設定する。画像解析部108は、オリジナル画像を解析し、階調変換カーブを作成するためのパラメータを計算する(ステップS402)。本実施例では階調処理方法として、固定の階調変換カーブをシフトさせて使用する方法を行う。シフトさせる量はオリジナル画像上での関心領域に最適な輝度が割り当てられるように決める。そこで画像解析部108は関心領域内の平均値を計算し出力する。
In steps S402 to S404, the
関心領域の設定はROI設定部として機能する設定部102により行われる。関心領域は画像の中心から数インチの正方形領域とする。これは、診断上重要と思われる部位を中心に撮影することが多いことに対応したものである。画像の中心付近を関心領域とすることで結果的に診断上注目すべき部位を関心領域とすることができる。関心領域を設定する別の方法としては、画像を解析することで撮影部位毎に適切な領域を設定することができる。例えば胸部の画像であれば肺野を抽出し、肺野の領域を関心領域とする。この技術は、必ずしも注目したい領域が画像の中央に位置しない場合があることに対応するものである。そのような場合としては例えば、撮影部位や目的によっては画像の複数の部位を注目する必要がある場合や、または予め撮影フォーマットが決まっている場合がある。これに対応して、撮影部位の情報に基づいて画像を解析し関心領域を設定することで、注目したい部位を基準に画像の階調を設定することができる。更に別の方法として、関心領域を線量の監視に用いる先述の所定領域に合わせて同じ領域に設定する方法がある。これにより、線量の制御の基準となる領域と画像処理の基準となる領域が揃うこととなるため、画像の品質をより適切に向上させることができる。
The region of interest is set by a
次に階調変換作成部109は解析により得られた関心領域内の平均値と、階調変換カーブ用の外部パラメータを取得し、階調変換カーブ(第1の階調変換関数)を作成する(ステップS403)。本実施例では階調変換カーブとして式1に示す画素値wのシグモイド関数S(w)を利用する。
Next, the gradation
式1においてpxSmin、pxSmaxはそれぞれS(w)の最小値、最大値である。sPrmCはS(w)の傾きで、sX0は、S(w)のw方向のシフト量である。
In Equation 1, pxSmin and pxSmax are the minimum value and maximum value of S (w), respectively. sPrmC is the slope of S (w), and sX0 is the shift amount of S (w) in the w direction.
次に階調変換作成部109はこの階調変換カーブS(w)のコントラストを示す値として微分値S’(w)を計算する(ステップS404)。この微分値がある入力画素値において所定の値以上である場合、例えば1より大きい場合には、オリジナル画像よりもコントラストが強調される。逆に1より小さい場合には、オリジナル画像よりも画素値の差が小さく扱われることとなるため、コントラストが弱められる。階調変換作成部109はこの微分値が所定の閾値より大きい画素値の範囲を、画像において強調する画素値の範囲として設定する。この閾値は上述のとおり1としてもよいが、1より大きい値としてもよく、画像や撮影部位などの情報に応じて決められる。以上ステップS402乃至S404の処理により設定部102は強調する画素値の範囲を設定する。
Next, the gradation
ステップS405では判定部104が階調変換カーブへの入力値と微分値S’(w)との間に指定の条件が成り立っているかを判定する(ステップS405)。指定の条件とは、ノイズが支配的な画素値の範囲を高コントラストに変換しないという条件である。ノイズが支配的な画素に高コントラストを与えたとしてもノイズだけが強調されて情報がなく、診断上不適切な画像となる。これを防ぐことを目的に、設定部102により設定された強調する画素値の範囲がイメージセンサの出力特性が良好な画素値の範囲を意味するダイナミックレンジに含まれるか否かを判定する。含まれると判定された場合(ステップS406でYes)には、画像の関心領域に基づいて設定された強調範囲がノイズやアーチファクトを強調することがない。よって設定部102で設定された範囲がそのまま強調される(ステップS408)。
In step S405, the
含まれないと判定された場合には(ステップS406でNo)、ノイズやアーチファクトを強調してしまい、画像として不自然なものとなる。さらには、作業者は画像を見ても適切な線量が照射されているのかいないのかを判断することができなくなってしまう。例えば、本来低線量であれば画像は暗くなるのに対して、低画素値のノイズの多い画素値範囲をコントラスト強調してしまうと、画像全体がざらついた画像となってしまう。これでは作業者は、適切な線量が照射されていないことは分かっても、線量を上げればいいのか下げればいいのかが分からなくなってしまう。そこで、関心領域に応じたコントラストの強調を放棄し、基準を超えた低線量であることを示すような画像となるように階調を変更する。そのために、ステップS407で補正部105により強調する画素値の範囲が補正される。
If it is determined that the image is not included (No in step S406), noise and artifacts are emphasized, resulting in an unnatural image. Furthermore, it becomes impossible for the operator to determine whether or not an appropriate dose is being applied even when viewing the image. For example, if the dose is originally low, the image becomes dark, but if the pixel value range with a low pixel value and noisy is contrast-enhanced, the entire image becomes rough. This makes it difficult for the operator to know that the appropriate dose has not been applied, but to increase or decrease the dose. Therefore, the contrast enhancement according to the region of interest is abandoned, and the gradation is changed so that an image showing a low dose exceeding the standard is obtained. For this purpose, the range of pixel values to be emphasized by the
ステップS406での判定部104の判定処理と、ステップS307での補正部105による補正処理を図5及び図6を参照しながら詳述する。まず、図5はイメージセンサの特性を示す図に基づいてイメージセンサ202のダイナミックレンジについて説明する。
横軸にイメージセンサへの到達線量(受光量)、縦軸はセンサ出力値とし、実線が信号成分、破線がノイズ成分を表している。イメージセンサの出力値は、この信号値とノイズ値を合わせた値となる。線量(受光量)E1まではノイズ成分の方が支配的であるが、そこからは信号成分が支配的となる。そして線量(受光量)E2に達すると信号成分は飽和する。このような関係の中で、E1からE2までの範囲をセンサのダイナミックレンジと定義する。
The determination process of the
The horizontal axis represents the arrival dose (light reception amount) to the image sensor, the vertical axis represents the sensor output value, the solid line represents the signal component, and the broken line represents the noise component. The output value of the image sensor is a value obtained by combining the signal value and the noise value. The noise component is dominant up to the dose (received light amount) E1, but the signal component is dominant from there. When the dose (received light amount) E2 is reached, the signal component is saturated. In such a relationship, the range from E1 to E2 is defined as the dynamic range of the sensor.
定性的に言えば、イメージセンサの信号値及びノイズ値の出力特性が良好な画素値または線量(受光量)の範囲である。また別の観点では、ダイナミックレンジの中ではSN比が閾値より大きい範囲である。また別の観点では、信号成分が飽和しない、または信号成分の出力の線型性が良好に保たれている画素値の範囲である。 Qualitatively speaking, it is a range of pixel values or doses (amount of received light) in which the output characteristics of signal values and noise values of the image sensor are good. From another viewpoint, the SN ratio is larger than the threshold in the dynamic range. From another viewpoint, the pixel value range is such that the signal component is not saturated or the linearity of the output of the signal component is kept good.
各到達線量に対応するセンサ出力値がOutDRMinとOutDRMaxである。OutDRMin以下の出力値となる場合、ノイズが支配的な値であることがいえる。一方、OutDRMax以上の出力値となる場合、センサ出力値が飽和し、センサ特有のアーチファクトが発生する。したがって、OutDRMin以下、OutDRMax以上となる値は、センサ出力値として有効ではないことがいえる。 The sensor output values corresponding to each arrival dose are OutDRMin and OutDRMax. When the output value is less than or equal to OutDRMin, it can be said that noise is a dominant value. On the other hand, when the output value is greater than or equal to OutDRMax, the sensor output value is saturated and artifacts peculiar to the sensor occur. Therefore, it can be said that values that are equal to or lower than OutDRMin and equal to or higher than OutDRMax are not effective as sensor output values.
次に図6に基づいて階調変換カーブの補正処理を説明する。入力画素値は、センサ出力値に修正を施して得られるオリジナル画像の画素値である。図5の線量E1に対応する入力画素値はwE1、線量E2に対応する入力画素値はwE2である。wE1とwE2の範囲をグレーで示している。これがイメージセンサ202のダイナミックレンジに対応する。すなわちセンサ出力値として有効な範囲はグレーゾーン内だけと判断できる。補正部105はこの範囲に階調変換カーブの高コントラストな領域があたるように制御する(ステップ407)。シグモイド関数S(w)の微分値S’(w)について、低コントラストから高コントラストへと変わる閾値をTh1とし、高コントラストから低コントラストへと変わる閾値をTh2とする。また、Th1とTh2に対応する入力画素値wをそれぞれwth1とwth2とする。この範囲が画像において強調される領域である。
Next, gradation correction curve correction processing will be described with reference to FIG. The input pixel value is a pixel value of the original image obtained by correcting the sensor output value. The input pixel value corresponding to the dose E1 in FIG. 5 is wE1, and the input pixel value corresponding to the dose E2 is wE2. The range of wE1 and wE2 is shown in gray. This corresponds to the dynamic range of the
判定部104は、これら2つの閾値wth1、wth2で定義される範囲がwE1からwE2の範囲内に包含されるか否かを判定する。ここで包含されないと判定された場合には、階調変換カーブを式(2)のように補正する。この処理は、階調変換カーブを低画素値シフトする処理である。
The
図5に示すようにwth1がwE1以下である場合は、ノイズだけを強調することになるため、階調変換カーブを高線量域側(高画素値側)へとシフトする。一方、wth2がwE2以上である場合は、センサのアーチファクトだけを強調することになるため、階調変換カーブを低線量域側(低画素値側)へとシフトする。
As shown in FIG. 5, when wth1 is equal to or less than wE1, only noise is emphasized, so the gradation conversion curve is shifted to the high dose range side (high pixel value side). On the other hand, when wth2 is greater than or equal to wE2, only the sensor artifact is emphasized, so the gradation conversion curve is shifted to the low dose region side (low pixel value side).
補正部105により決められるシフト量は、強調される画素値の範囲がダイナミックレンジの範囲内に含まれるように決められる。強調される画素値の範囲がダイナミックレンジにぎりぎり含まれるようにシフトする。例えば、ダイナミックレンジよりも低画素値側を強調するよう設定されていた場合には、強調する範囲の下限がダイナミックレンジの下限と一致するようにシフト量を決める。このようにすると、シフト量は最小に抑えられることとなる。シフト量を最小に抑えることで画像の見えやすさを極力保つことができる。
The shift amount determined by the
以上の処理により、階調変換カーブが補正されることに応じて強調する画素値の範囲も補正される。このようにして、強調する範囲がダイナミックレンジに範囲内に含まれない場合には強調する範囲をシフトする処理が行われる。 Through the above processing, the pixel value range to be emphasized is also corrected in accordance with the correction of the gradation conversion curve. In this way, when the range to be emphasized is not included in the dynamic range, processing for shifting the range to be emphasized is performed.
画像のダイナミックレンジの範囲が十分に大きい場合には、画素値の範囲をシフトさせるだけではノイズ等の強調を防ぐことができない場合がある。信号出力の線型性のよい画素値の範囲に収めることよりもSN比が閾値より大きい画素値の範囲に収めることを優先させてシフトさせる。こうすることで強調される画素値の範囲が高画素値側にシフトされ白とびした画像となるため、手動で線量をより高く設定されてしまうのを防ぐことができる。 When the range of the dynamic range of the image is sufficiently large, it may not be possible to prevent enhancement of noise or the like only by shifting the pixel value range. The shift is performed with priority given to the range of pixel values with a signal-to-noise ratio larger than the threshold value, rather than the range of pixel values with good linearity of signal output. By doing so, the range of pixel values to be emphasized is shifted to the high pixel value side and an overexposed image is obtained, so that it is possible to prevent the dose from being set higher manually.
図4のステップS408にて、階調処理部106は、オリジナル画像に階調変換を施す。設定部102により設定された画素値の範囲がダイナミックレンジに包含される場合には、設定された範囲を強調する階調変換カーブS(w)(第1の階調変換関数)を画像に適用する。一方でダイナミックレンジに包含されない範囲がある場合には、包含されるように画素値の範囲がシフトされた階調変換カーブSnew(w)(第2の階調変換関数)を適用する。
In step S408 of FIG. 4, the
ここで、オリジナル画像において特定の二次元領域が決められていた場合は、その特定の領域のみをSnew(w)にて変換し、特定の領域外はS(w)にて変換する。これにより特定の領域だけはノイズやアーチファクトを強調しない、線量の過不足が明らかな部分画像として表示される。 Here, when a specific two-dimensional area is determined in the original image, only the specific area is converted by Snew (w), and outside the specific area is converted by S (w). As a result, only a specific region is displayed as a partial image in which noise and artifacts are not emphasized and the excess or deficiency of the dose is clear.
階調変換した画像は、表示制御部107が表示部210に階調変換された画像を表示させる。その他画像処理はメモリ保存、PACS転送などを行い、ユーザに提供する(ステップS409)。
The gradation-converted image is displayed by the
上述の処理により、撮影線量が足りずノイズだけが強調されてしまった画像が、ノイズが強調されることのない画像となり、線量が不足していることが明らかな暗い画像が出力される。撮影線量が多すぎてセンサの特性に起因するアーチファクトが強調されてしまったような画像が、センサの特性が強調されない画像となり、線量が多すぎであることが明らかな明るい画像が出力される。このように、線量が不足または過多で画像処理による画質の改善が望めない場合には、線量が不足または過多であることを示す画像処理を画像に施す。これによりノイズやアーチファクトを強調した不自然な画像が出力されてしまう状況を減らすことができる。また、作業者が画像を見ることで線量が不足しているか過多であるのかが一目瞭然となるため、線量の制御を適切に行うことができる。 By the above-described processing, an image in which only the noise is emphasized due to insufficient imaging dose becomes an image in which noise is not enhanced, and a dark image that clearly shows that the dose is insufficient is output. An image in which artifacts due to the characteristics of the sensor are emphasized due to an excessive imaging dose becomes an image in which the characteristics of the sensor are not enhanced, and a bright image that clearly shows that the dose is excessive is output. As described above, when the image quality cannot be improved by image processing due to insufficient or excessive dose, image processing indicating that the dose is insufficient or excessive is performed on the image. As a result, it is possible to reduce the situation in which an unnatural image in which noise and artifacts are emphasized is output. Moreover, since it becomes clear at a glance whether the dose is insufficient or excessive when the operator looks at the image, it is possible to appropriately control the dose.
なお関心領域の画素値などの情報から、線量が不足した画像であることが明らかである場合には、特性取得部103はダイナミックレンジの情報として下限値E1またはwE1の値のみを取得する。この場合ダイナミックレンジの情報としてはSN比が閾値よりも大きい画素値の範囲を示す情報となる。
If it is clear from the information such as the pixel value of the region of interest that the dose is insufficient, the
また別の例として、イメージセンサのダイナミックレンジが十分に広い場合や、イメージセンサの信号成分の出力について線型性が保たれる範囲が十分に広い場合も考えられる。この場合にも特性取得部103はダイナミックレンジの情報として下限値E1またはwE1の値のみを取得することとしてもよい。
As another example, there may be a case where the dynamic range of the image sensor is sufficiently wide, or a case where the linearity of the output of the signal component of the image sensor is sufficiently wide. Also in this case, the
実施例2では、強調される画素値の範囲が前記ダイナミックレンジに包含されるように補正するため、階調変換曲線をシフトさせる代わりに曲線の傾きを変更する。ハードウェア構成は実施例1とほぼ同様であるが、補正部105による処理が異なっている。実施例2を図7を参照しながら説明する。
In the second embodiment, in order to correct the emphasized pixel value range to be included in the dynamic range, the slope of the curve is changed instead of shifting the gradation conversion curve. The hardware configuration is almost the same as in the first embodiment, but the processing by the
図1の構成図と図4の全体フローに基づいて実施例2の処理の流れを説明する。ステップS401からステップS405までは実施例1と同様であるため省略する。図7に階調変換カーブでの変換を示した。入力画素値は、センサ出力値に修正を施した値である。図5の線量E1に対応する入力画素値はwE1、線量E2に対応する入力画素値はwE2である。wE1とwE2の範囲をグレーで示している。すなわちセンサ出力値として有効な範囲はグレーゾーン内だけと判断できる。補正部105はこの範囲に階調変換カーブの高コントラストな領域があたるように制御する(ステップ406)。微分値S’(w)上にて、低コントラストと高コントラストとの閾値をそれぞれTh1とTh2とするとする。
入力画素値wE1とwE2での傾きがTh1、Th2となるようなS(w)の傾きsPrmCth1、sPrmCth2に設定することで制御ができる。すなわち階調変換カーブを構成するf(w)を式(3)のように補正する。
The processing flow of the second embodiment will be described based on the configuration diagram of FIG. 1 and the overall flow of FIG. Steps S401 to S405 are the same as those in the first embodiment, and will not be described. FIG. 7 shows the conversion using the gradation conversion curve. The input pixel value is a value obtained by correcting the sensor output value. The input pixel value corresponding to the dose E1 in FIG. 5 is wE1, and the input pixel value corresponding to the dose E2 is wE2. The range of wE1 and wE2 is shown in gray. That is, it can be determined that the effective range as the sensor output value is only in the gray zone. The
Control can be performed by setting the slopes sPrmCth1 and sPrmCth2 of S (w) such that the slopes at the input pixel values wE1 and wE2 are Th1 and Th2. That is, f (w) constituting the gradation conversion curve is corrected as shown in Expression (3).
図7に示すようにwE1におけるコントラストが高い場合は、ノイズだけを強調することになるため、階調変換カーブの傾きを小さくする。一方、wE2におけるコントラストが高い場合は、センサのアーチファクトだけを強調することになるため、階調変換カーブの傾きを小さくする。式(3)では全ての入力画素に対し一律で傾きを変えているが、S(w)の中心を境に傾きの異なる階調変換カーブなどへと変換してもよい。(ステップS407)。ステップS408以降は実施例1と同様であるため省略する。
As shown in FIG. 7, when the contrast at wE1 is high, only noise is emphasized, so that the gradient of the gradation conversion curve is reduced. On the other hand, when the contrast at wE2 is high, only the sensor artifact is emphasized, so that the gradient of the gradation conversion curve is reduced. In equation (3), the gradient is uniformly changed for all input pixels, but it may be converted into a gradation conversion curve or the like having a different gradient from the center of S (w). (Step S407). Since step S408 and subsequent steps are the same as those in the first embodiment, description thereof is omitted.
また別の例として、wth1とwth2の範囲に収まらないときは、wth1とwth2の範囲に階調変換カーブを圧縮(もしくは入力画素値を圧縮)する。 As another example, when it does not fall within the range of wth1 and wth2, the gradation conversion curve is compressed (or the input pixel value is compressed) to the range of wth1 and wth2.
このように階調変換曲線をシフトさせるだけでなく傾きを変更したり画素値方向に圧縮したりすることでも、ノイズやアーチファクトを強調した不自然な画像を出力させないようにすることができる。 Thus, not only shifting the gradation conversion curve but also changing the inclination or compressing in the pixel value direction can prevent an unnatural image with enhanced noise and artifacts from being output.
実施例3では、線量制御に応じたコントラスト強調処理を行うか否かを切り替える機能を有している。ハードウェア構成は実施例1とほぼ同様であるが、階調処理部106は表示部207に表示されたGUIへの入力に応じて処理を行う点が異なっている。以下本実施例を図8乃至図10を参照しながら説明する。
The third embodiment has a function of switching whether or not to perform contrast enhancement processing according to dose control. The hardware configuration is almost the same as that of the first embodiment, except that the
ユーザ入力からの入力を得るための調整画面800を図8に基づいて説明する。制御モードの選択GUI802は例えば図8のようにGUIのチェックボックスを含む。選択GUIでの選択に応じて、判定部104が上述の実施例に記載された画像処理を行うか否かを選択する。ユーザがこのチェックボックスをオンにすると、図3のフローチャートに示す画像処理が行われるが、オフにすると図3の処理は行われない。輝度を調整するための調整GUI803はユーザが輝度を調整するための図8のようなGUIを含む。この調整GUI803への入力に応じて補正部105は階調変換カーブをシフトされる。コントラストを調整する調整GUI804はユーザがコントラストを調整するための図9のようなGUIを含む。この調整GUI804への入力に応じて補正部105により階調変換カーブが変更され、コントラストが調整される。
An adjustment screen 800 for obtaining input from user input will be described with reference to FIG. The control mode selection GUI 802 includes a GUI check box as shown in FIG. 8, for example. In accordance with the selection on the selection GUI, the
調整GUI803の詳細を図9に基づいて説明する。調整GUI803には図9に図示するように輝度調整ステップが中央(Step0)に近い位置からStep0、Step2、Step−2、Step4、・・・、StepTHと表示されている。また、各輝度調整ステップに対応して輝度シフト量が中央に近い位置から0、5、−5、10、−10、・・・・と5刻みで表示されている。例えば、Step2は輝度シフト量10に対応し、StepTHは輝度シフト量−30に対応している。ユーザが輝度調整ステップを選択すると、補正部105は対応した輝度シフト量分だけ階調変換カーブがシフトされ輝度が調整される。
Details of the adjustment GUI 803 will be described with reference to FIG. In the adjustment GUI 803, as shown in FIG. 9, the brightness adjustment step is displayed as
調整GUI803では、あるステップを超えると階調変換カーブが指定のコントラスト内に収まらなくなる境界となるステップをStepTHとしている。StepTHよりも(絶対値が)大きく輝度を変更する場合、診断上ふさわしくない画像が作成される。 In the adjustment GUI 803, Step TH is defined as a step which becomes a boundary where the gradation conversion curve does not fall within the designated contrast when a certain step is exceeded. When the luminance is changed to be larger than that of StepTH (absolute value), an image that is not suitable for diagnosis is created.
ここで、任意のステップが選択された場合に、輝度シフト量の最小単位を小さくするため、輝度シフト量と輝度調整ステップとの再割当てが実行される。これは、表示制御部107により行われる。図9の下段には、輝度調整ステップStep0が選択された場合の再割当てを示している。図9ではステップが負の方向にいく場合に指定のコントラスト以上になる。したがってStepTHからStep0までの間を3刻みで割り当てる。この結果、StepTHからStep0までの間で詳細に調整することができる。コントラスト調整部904においても同様の調整機構を採用することが可能である。
Here, when an arbitrary step is selected, reallocation of the luminance shift amount and the luminance adjustment step is executed in order to reduce the minimum unit of the luminance shift amount. This is performed by the
上述の調整機能を利用する際の画像処理装置100が実行する処理の流れを図10のフローチャートに従い説明する。まず、階調変換作成部109がユーザによって輝度・コントラスト調整された量を取得する(ステップS1001)。次に階調変換作成部109は取得した輝度調整量を使って階調変換カーブを更新する。そして判定部104は、センサダイナミックレンジ外の範囲にある画素値が指定のコントラストの範囲内であるかを計算する(ステップS1002)。ここで判定部104が、指定のコントラストの範囲内であるか否かを判定する(ステップS1003)。この処理は上述の実施例と同様である。ここで指定のコントラストの範囲内であると判定された場合には(ステップS1003でYes)、ステップS806へと進む。一方、範囲内ではないと判定すると(ステップS1003でNo)、引き続き判定部104は現在の階調制御モードがONであるかOFFであるかを判定する(ステップS1004)。ここでの判定は、調整画面制御モードの選択GUI802がONの時は、センサダイナミックレンジの範囲外にある入力画素が指定のコントラスト以上である場合に、以下となるように制御する。実施例の1や2と同様に階調変換カーブを修正する(ステップS1005)。OFFの時は制御をしない。つまりステップS1002で作成した階調変換カーブで処理をする(ステップS1006)。そして表示制御部107は階調処理済みの画像を画像表示GUI801に表示させる(ステップS1007)。またステップS1008で表示制御部107は上述の図9に示したとおり、調整GUI803及び804のシフト量と調整ステップとの再割当て処理を行う。
The flow of processing executed by the image processing apparatus 100 when using the adjustment function described above will be described with reference to the flowchart of FIG. First, the gradation
以上のとおり本実施例では、イメージセンサのダイナミックレンジに応じたコントラスト強調処理を行うか行わないかを選択可能とし、ユーザが調整可能とする。これにより適切な画像処理を画像に施すことができる。 As described above, in this embodiment, it is possible to select whether or not to perform the contrast enhancement processing according to the dynamic range of the image sensor, and the user can adjust. Accordingly, appropriate image processing can be performed on the image.
実施例4では、先述の実施例において、強調される画素値の範囲がダイナミックレンジの範囲内に含まれていないと判定された場合には、表示制御部107で表示部210に警告を表示させる。システムの構成及び処理の流れについて、先述の実施例と重複する部分については説明を省略する。
In the fourth embodiment, when it is determined in the above-described embodiment that the range of pixel values to be emphasized is not included in the range of the dynamic range, the
関心領域に基づいて設定された画像において強調される画素値の範囲がダイナミックレンジの範囲内に含まれないと判定部104で判定された場合、適切な階調変換が行われない旨の警告を表示する。
If the
または、設定部102で設定された画素値の範囲がダイナミックレンジよりも低画素値側の範囲を包含する場合には、線量が不足している旨の表示をする。この表示は、「線量が不足している」旨を文字情報で知らせる表示でも、線量不足であることを示すアイコンでもその両方でもよい。また表示部210には警告の表示と並べて、実施例1に示す処理により階調処理部106にて処理された画像が表示されることとなる。この画像は設定部102にて設定された範囲がダイナミックレンジを外れている場合には、線量不足であることが明らかな階調変換が施された画像である。そのためX線画像と、文字情報及び/またはアイコンとの両方によって線量不足であることをユーザに分かりやすく提示することができる。
Alternatively, when the pixel value range set by the
同様に設定部102で設定された画素値の範囲がダイナミックレンジよりも高画素値側の範囲を包含する場合には、線量が過多となっている旨の警告表示をする。この場合でも線量過多であることが強調された画像と警告表示が並べて表示されるため、線量が過多であることをユーザに分かりやすく提示することができる。
Similarly, when the range of pixel values set by the
また実施例3のステップS1003の処理で、指定のコントラストの範囲に含まれない場合は、線量不足または線量過多である旨の警告を出してもよい。警告を出し、ユーザに階調変換カーブを修正せずに処理をするか否かを選択させるメッセージを表示してもよい。もしくは警告を出すだけで表示を行ってしまってもよい。 Further, in the process of step S1003 of the third embodiment, if the specified contrast range is not included, a warning indicating that the dose is insufficient or excessive may be issued. A warning may be issued and a message may be displayed that allows the user to select whether or not to process without correcting the tone conversion curve. Alternatively, it may be displayed only by issuing a warning.
また別の例として、手動での線量制御を支援する方法に適用することができる。動画撮影をする際に、最初ダイナミックレンジを外れた低線量で撮影を開始する。この場合、表示部210には線量が不足していることが明らかな画像と警告表示とが表示部210に表示される。そこで作業者は徐々に線量を上昇させ、ダイナミックレンジにぎりぎり含まれるようななるべく少ない線量に調節していく。線量が適切となったところで線量の変更を終了する。特に放射線撮影においては被写体の被曝を抑えるためなるべく低線量で撮影したいという要求がある。本発明は係る手動での線量制御を支援することができる。
As another example, the present invention can be applied to a method for supporting manual dose control. When shooting a movie, start shooting at a low dose outside the dynamic range. In this case, the
上述の実施例では、強調する範囲がイメージセンサのダイナミックレンジを外れてしまった場合にもコントラストを強調する処理を行っていたが、実施例5では当該場合に強調処理を省略する。図11を参照しながら本実施例を説明する。なお図11にて図1と同様の名称を付されたブロックは図1と同等であるため説明を省略する。画像処理装置1100を含むシステムは図2と同様であるため説明を省略する。 In the above-described embodiment, the process of enhancing the contrast is performed even when the emphasized range is out of the dynamic range of the image sensor. However, in the fifth embodiment, the enhancement process is omitted in this case. The present embodiment will be described with reference to FIG. 11 with the same names as those in FIG. 1 are the same as those in FIG. The system including the image processing apparatus 1100 is the same as that shown in FIG.
画像処理装置1100の制御部1110は判定部1104の判定結果を取得する。ここで、制御部1110は階調処理部1106に処理を実行させるか否かを制御する。制御部1110は設定されたコントラストを強調する範囲がイメージセンサのダイナミックレンジに包含されない場合には、階調処理部1106に対して所定の信号を送る。一方、コントラストを強調する範囲がイメージセンサのダイナミックレンジに包含されない場合には、当該所定の信号を送信しないよう制御する。この所定の信号を入力した場合、階調処理部1106は入力された画像に対して階調変換を施さず表示制御部1107に送る。所定の信号が入力されない場合には、先述の実施例と同様に階調処理部1106は設定された範囲のコントラストを強調する処理を実行する。このようにして階調処理部1106に入力された画像はコントラスト強調処理が施されないまま表示制御部1107に入力され、表示部207に表示される。
The
これにより画像の線量が不足している場合や、線量が多すぎる場合にはコントラストを強調する処理が施されないこととなる。これにより、画像から線量不足であるか、線量過多であるかを示す情報をシステムの操作者に直感的に確認させることができる。システムの操作者は線量の調節やイメージセンサのゲイン切り替え、または撮影の中止など必要に応じた対処を行うことができる。 As a result, when the dose of the image is insufficient or when the dose is too high, the processing for enhancing the contrast is not performed. As a result, the system operator can intuitively confirm information indicating whether the dose is insufficient or excessive from the image. The system operator can take necessary measures such as adjusting the dose, switching the gain of the image sensor, or stopping the imaging.
本実施例では、あらかじめ作成した第1の階調変換特性を補正して第2の階調変換特性処理を得ることをせず、画像から直接に階調変換特性を決定する。 In this embodiment, the gradation conversion characteristic is determined directly from the image without correcting the first gradation conversion characteristic created in advance and obtaining the second gradation conversion characteristic process.
図12に画像処理装置1200の構成を示す。画像入力部1201、階調処理部1206、表示制御部1207は上述の実施例と同様である。図2に示すハードウェアと、ROM205または記憶部206に記憶されたソフトウェアプログラムとが共同して各部の機能を実現する。以下上述の実施例と異なる点を説明する。
FIG. 12 shows the configuration of the image processing apparatus 1200. The
本実施例では、画像解析部1208はX線画像を解析して階調変換特性を決めるための基準の画素値を取得する。基準の画素値の一例は被写体領域の平均値または中央値等の統計値を用いる。被写体領域は、画像解析部1208がX線画像から周知のアルゴリズムにより被写体領域を抽出することで得られる。被写体領域を基準に階調変換を行うことが可能となるため、被写体領域を診断しやすい画像を得ることができる。
In this embodiment, the
また別の例としては、画像中心を含む部分領域内の統計値を用いても良い。X線撮影においては写損を防ぐ観点から被写体の位置決めを行ってから撮影を行うことが多いため画像の中心付近に診断上重要な部位がくることが多い。本例はこれを利用したものであり、被写体抽出処理が不要となるため、解析処理の負荷を低減することができる。 As another example, a statistical value in a partial region including the center of the image may be used. In X-ray imaging, from the viewpoint of preventing image loss, imaging is often performed after the subject is positioned, and thus there are many diagnostically important parts near the center of the image. This example uses this, and the subject extraction process becomes unnecessary, so the load of the analysis process can be reduced.
記憶部206には、イメージセンサ202の出荷前等のタイミングで予め取得されたダイナミックレンジの範囲が記憶される。なお、特性の変化などに対応して、ダイナミックレンジを取得する特性取得部1203を設けることとしてもよい。この場合特性取得部1203は、にイメージセンサ202から得られる画像についてSN比が閾値より大きくなる線量に対応する画素値の範囲を取得し、記憶部206に記憶する。また、ダイナミックレンジの情報として、先述の図5に示すようにイメージセンサ202により得られる画像のSN比が閾値より大きい画素値の範囲でありかつ信号成分の出力の線型性が保たれる画素値の範囲を記憶部206に記憶することとしても良い。これには、たとえば所定の規則で線量透過率が空間的に変化するファントムを撮影した画像について、SN比及びセンサに到達した入射線量に対する出力画素値の変化を計測することで得ることができる。
The
判定部1204は、画像解析部1208により得られた階調変換の基準の画素値が、記憶部206に記憶されたダイナミックレンジに含まれるか否かを判定する。基準の画素値がダイナミックレンジに含まれる場合には、決定部として機能する階調変換作成部1209は、画像解析部1208により得られた基準の画素値を基準に階調変換特性が決定する。この場合、階調変換により階調が広げられる画素値の範囲の中心が画像解析部1208により得られた基準の画素値をとなる。
The
画像解析部1208により得られた基準の画素値が前記記憶された画素値の範囲に含まれると判定されなかった場合には、階調変換作成部1209はダイナミックレンジに含まれる特定の画素値を新たな基準の画素値として選ぶ。この新たな基準の画素値が、階調変換により階調が広げられる画素値の範囲の中心となるように階調変換の特性が決定される。
When it is not determined that the reference pixel value obtained by the
上述の機能により実現される処理の例を図13、図14に基づき説明する。図13(a)、図13(a)は本実施例の処理を示す図である。範囲1301はX線画像の画素値の範囲を示す。ダイナミックレンジ300よりも範囲1301が大きくなっているのは、X線画像の場合被写体を通過せず直接X線がイメージセンサ202に到達する場合があるからである。画像特徴量の値は、たとえば被写体領域に基づいて得られた階調変換の基準となる画素値の値である。特に低線量の撮影を行う場合、この特徴量の値がダイナミックレンジ300よりも下回る場合があり、その場合には、階調カーブの基準点をダイナミックレンジ300に含まれる特定の画素値とする。この基準の値が、階調カーブにより階調が広げられる画素値の範囲の中心位置となる。図13(a)の例では、特定の画素値はダイナミックレンジ300の下限値となっている。図14(a)が、図13(a)の処理で決定される階調カーブの例である。これに限らず、ダイナミックレンジ300の下限値よりも所定の画素値だけ大きい画素値を基準の値としてもよい。
An example of processing realized by the above function will be described with reference to FIGS. FIG. 13A and FIG. 13A are diagrams showing the processing of this embodiment. A range 1301 indicates a range of pixel values of the X-ray image. The reason why the range 1301 is larger than the
上述の処理により画像特徴量の値を基準に決められる階調カーブを用いるよりも、ノイズの領域を抑圧することができる。 The noise region can be suppressed rather than using the gradation curve determined based on the value of the image feature amount by the above-described processing.
図13(b)は画像特徴量の値がダイナミックレンジ300に包含される場合を示したものである。この場合画像特徴量の値がそのまま階調変換の基準の値となる。この場合は、被写体領域の階調を適切に広げ、診断しやすい画像を提供することができる。図14(b)は図13(b)の場合に対応する階調カーブの例である。
FIG. 13B shows a case where the image feature value is included in the
本実施例の処理により、被写体の領域を反映した階調処理を行うことができる。また、適正な線量を下回る線量で撮影した場合には、ノイズの抑圧を優先し被写体領域の少なくとも一部の階調を狭める。これにより、不自然にノイズが強調された画像が表示されないようにすることができる。 By the processing of this embodiment, gradation processing that reflects the area of the subject can be performed. In addition, when imaging is performed at a dose lower than an appropriate dose, priority is given to noise suppression, and at least part of the gradation of the subject area is narrowed. Thereby, it is possible to prevent an image with unnaturally emphasized noise from being displayed.
(その他の実施例) 本実施例では強調される画素値の範囲を階調変換曲線によって設定したが、LUTによって行っても、また階調変換曲線やLUT等の階調変換によらずともよい。 Other Embodiments In this embodiment, the range of pixel values to be emphasized is set by a gradation conversion curve. However, it may be performed by an LUT or not by gradation conversion such as a gradation conversion curve or LUT. .
また、予めイメージセンサの出力特性が良好な画素値の範囲を抑圧するコントラスト変更処理を行うこととしてもよい。また、関心領域の画素値の範囲に基づいて設定された強調される画素値の範囲がイメージセンサの出力特性が良好な範囲よりも低線量側の領域を含む場合には、素抜け領域を除いてコントラストを抑圧することで、低線量であることが診断者にとって分かりやすい画像を生成することができる。 In addition, a contrast change process for suppressing a range of pixel values with good output characteristics of the image sensor may be performed in advance. In addition, if the range of pixel values to be emphasized set based on the range of pixel values in the region of interest includes a region on the lower dose side than the range where the output characteristics of the image sensor are good, exclude the missing region. By suppressing the contrast, it is possible to generate an image that is easy for the diagnostician to understand that the dose is low.
画像処理装置100または画像処理装置1100で行われている処理を複数の装置で分散させ画像処理システムとして実現してもよい。また図1または図10、図12に示す1つの機能ブロックとしてまとめられている処理を複数の回路または機能ブロックで分散させて実現してもよい。またイメージセンサと、上述の実施例における画像処理装置及び表示装置の機能とを組み込んだFPDまたは撮像装置により本実施例のシステムを実現することとしてもよい。 The processing performed in the image processing apparatus 100 or the image processing apparatus 1100 may be distributed by a plurality of apparatuses and realized as an image processing system. Further, the processing grouped as one functional block shown in FIG. 1, FIG. 10, or FIG. 12 may be realized by being distributed by a plurality of circuits or functional blocks. Further, the system of the present embodiment may be realized by an FPD or an imaging device incorporating the image sensor and the functions of the image processing device and the display device in the above-described embodiment.
図2で示すように本実施例をコンピュータのハードウェアとソフトウェアの協働により実現する場合、複数のCPUで分散させて本実施例を実現することとしてもよい。ここでコンピュータ上で稼動しているオペレーティングシステム(OS)などが実際の処理の一部または全部を行い、その処理によって上述の機能が実現してもよい。またソフトウェアであるプログラムまたはプログラムコードを記録した記録媒体もまた本実施例を構成する。ここでいう記録媒体は、non−transitoryなものであり、電波等記録対象が固定されないは様な媒体は含まないが、キャッシュや揮発メモリ等は含む概念である。 As shown in FIG. 2, when the present embodiment is realized by the cooperation of computer hardware and software, the present embodiment may be realized by being distributed by a plurality of CPUs. Here, an operating system (OS) or the like running on the computer may perform part or all of the actual processing, and the above-described functions may be realized by the processing. A recording medium on which a program or program code which is software is recorded also constitutes the present embodiment. The recording medium here is a non-transitory medium and does not include a medium in which a recording target such as a radio wave is not fixed, but includes a cache, a volatile memory, and the like.
上述した本実施の形態における記述は、実施例の好適な画像処理装置の一例であり、実施の形態はこれに限定されるものではない。 The description in this embodiment described above is an example of a suitable image processing apparatus of the embodiment, and the embodiment is not limited to this.
以上、上述した実施例では、被写体に照射されるX線の線量が線量である場合に得られたX線画像について、X線画像における被写体領域の階調を狭めることでノイズを抑圧することができる。このように被写体領域の階調を広げることよりもノイズを抑圧することを優先した階調処理を行うことで、不自然にノイズが強調された画像が表示されないようにすることができる。 As described above, in the above-described embodiment, noise can be suppressed by narrowing the gradation of the subject region in the X-ray image for the X-ray image obtained when the X-ray dose irradiated to the subject is a dose. it can. Thus, by performing gradation processing that prioritizes suppression of noise over widening the gradation of the subject area, it is possible to prevent an image with unnaturally emphasized noise from being displayed.
100 画像処理装置
101 画像入力部
102 設定部
103 特性取得部
104 判定部
105 補正部
106 階調処理部
107 表示制御部
200 放射線撮影システム
DESCRIPTION OF SYMBOLS 100
Claims (14)
X線センサのダイナミックレンジに対応し、信号成分の出力が飽和しない画素値の範囲を記憶する記憶手段と、
前記階調変換が行われる画素値の範囲が、前記記憶された画素値の範囲に含まれるか否かを判定する判定手段と、
前記判定の結果に応じて階調変換の特性を補正する補正手段と、
を有することを特徴とする画像処理装置。 Setting means for setting a range of pixel values to be subjected to gradation conversion based on an X-ray image captured by an X-ray sensor ;
Storage means for storing a range of pixel values corresponding to the dynamic range of the X-ray sensor and in which the output of the signal component is not saturated ;
Determination means for determining whether or not a range of pixel values to which the gradation conversion is performed is included in the range of stored pixel values;
Correction means for correcting the characteristics of gradation conversion according to the result of the determination;
An image processing apparatus comprising:
放射線源と、
前記放射線源の射出するX線を検出するX線センサと、
前記X線センサにより得られ前記画像処理装置で処理された画像を表示する表示手段と、
を有することを特徴とする放射線撮影システム。 An image processing apparatus according to any one of claims 1 to 9,
A radiation source;
An X-ray sensor for detecting X-rays emitted from the radiation source;
Display means for displaying an image obtained by the X-ray sensor and processed by the image processing device;
A radiation imaging system comprising:
X線センサのダイナミックレンジに対応し、信号成分の出力が飽和しない画素値の範囲を取得する取得ステップと、
前記階調変換が行われる画素値の範囲が、前記X線センサのダイナミックレンジに対応する画素値の範囲に含まれるか否かを判定する判定ステップと、
前記判定の結果に応じて階調変換の特性を補正する補正ステップと、
を有することを特徴とする画像処理方法。 A setting step for setting a range of pixel values to be subjected to gradation conversion based on an X-ray image captured by an X-ray sensor ;
An acquisition step for acquiring a range of pixel values corresponding to the dynamic range of the X-ray sensor and in which the output of the signal component is not saturated ;
A determination step of determining whether a range of pixel values to which the gradation conversion is performed is included in a range of pixel values corresponding to a dynamic range of the X-ray sensor;
A correction step of correcting the characteristics of gradation conversion according to the result of the determination ;
An image processing method comprising:
X線センサのダイナミックレンジに対応し、X線センサの特性によりSN比が閾値より大きくなる画素値の範囲を記憶する記憶手段と、
前記階調変換が行われる画素値の範囲が、前記記憶された画素値の範囲に含まれるか否かを判定する判定手段と、
前記判定の結果に応じて階調変換の特性を補正する補正手段と、
を有することを特徴とする画像処理装置。 Setting means for setting a range of pixel values to be subjected to gradation conversion based on an X-ray image captured by an X-ray sensor;
Storage means for storing a range of pixel values corresponding to a dynamic range of the X-ray sensor and having an SN ratio larger than a threshold value due to the characteristics of the X-ray sensor;
Determination means for determining whether or not a range of pixel values to which the gradation conversion is performed is included in the range of stored pixel values;
Correction means for correcting the characteristics of gradation conversion according to the result of the determination;
An image processing apparatus comprising:
X線センサのダイナミックレンジに対応し、X線センサの特性によりSN比が閾値より大きくなる画素値の範囲を取得する取得ステップと、
前記階調変換が行われる画素値の範囲が、前記X線センサのダイナミックレンジに対応する画素値の範囲に含まれるか否かを判定する判定ステップと、
前記判定の結果に応じて階調変換の特性を補正する補正ステップと、
を有することを特徴とする画像処理方法。 A setting step for setting a range of pixel values to be subjected to gradation conversion based on an X-ray image captured by an X-ray sensor;
An obtaining step for obtaining a range of pixel values corresponding to a dynamic range of the X-ray sensor and having an SN ratio larger than a threshold value due to the characteristics of the X-ray sensor;
A determination step of determining whether a range of pixel values to which the gradation conversion is performed is included in a range of pixel values corresponding to a dynamic range of the X-ray sensor;
A correction step of correcting the characteristics of gradation conversion according to the result of the determination;
An image processing method comprising:
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010275145 | 2010-12-09 | ||
JP2010275145 | 2010-12-09 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011237570A Division JP5868119B2 (en) | 2010-12-09 | 2011-10-28 | Image processing apparatus, radiation imaging system, image processing method, and recording medium |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2016104196A JP2016104196A (en) | 2016-06-09 |
JP6132938B2 true JP6132938B2 (en) | 2017-05-24 |
Family
ID=56101968
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016002063A Active JP6132938B2 (en) | 2010-12-09 | 2016-01-07 | Image processing apparatus, radiation imaging system, image processing method, and recording medium |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6132938B2 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108852401B (en) * | 2018-07-11 | 2022-04-15 | 上海联影医疗科技股份有限公司 | Correction table generation method and device, CT system and storage medium |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2952519B2 (en) * | 1990-11-21 | 1999-09-27 | コニカ株式会社 | Radiation image gradation converter |
JP3312281B2 (en) * | 1992-09-03 | 2002-08-05 | コニカ株式会社 | Radiation image processing device |
JP3375237B2 (en) * | 1995-08-29 | 2003-02-10 | 株式会社モリタ製作所 | Automatic density correction method for X-ray equipment |
JP3814864B2 (en) * | 1996-04-02 | 2006-08-30 | コニカミノルタホールディングス株式会社 | Radiographic image processing condition determining apparatus, image processing apparatus, radiographic image processing condition determining method, and image processing method |
JPH1032756A (en) * | 1996-07-17 | 1998-02-03 | Fuji Photo Film Co Ltd | Method and device for processing gradation of image information |
JP3690882B2 (en) * | 1996-08-16 | 2005-08-31 | 富士写真フイルム株式会社 | Image enhancement processing method and apparatus |
JP2000030046A (en) * | 1998-04-20 | 2000-01-28 | Konica Corp | Radiation image detecting and processing apparatus |
JP3805107B2 (en) * | 1998-06-24 | 2006-08-02 | 株式会社日立メディコ | Measuring device using X-ray |
JP3800892B2 (en) * | 1999-10-29 | 2006-07-26 | コニカミノルタホールディングス株式会社 | Radiation image processing device |
JP2002314881A (en) * | 2001-04-13 | 2002-10-25 | Toshiba Corp | X-ray diagnostic system |
FR2847698B1 (en) * | 2002-11-27 | 2005-05-06 | Ge Med Sys Global Tech Co Llc | METHOD FOR MANAGING THE DYNAMICS OF A DIGITAL RADIOLOGICAL IMAGE |
JP4626138B2 (en) * | 2003-09-30 | 2011-02-02 | コニカミノルタエムジー株式会社 | Image processing apparatus, image processing method, and program |
JP2006079590A (en) * | 2004-08-10 | 2006-03-23 | Konica Minolta Medical & Graphic Inc | Image processor |
JP2006102044A (en) * | 2004-10-04 | 2006-04-20 | Konica Minolta Medical & Graphic Inc | Medical image processing system and medical image processing method |
EP1820447B1 (en) * | 2004-12-10 | 2012-06-06 | Konica Minolta Medical & Graphic, Inc. | Image processing device, image processing method, and image processing program |
JP2006247204A (en) * | 2005-03-11 | 2006-09-21 | Fuji Photo Film Co Ltd | Medical image processing apparatus and method |
JP2008212493A (en) * | 2007-03-06 | 2008-09-18 | Canon Inc | Image processor and image processing method |
JP5275668B2 (en) * | 2008-04-04 | 2013-08-28 | 富士フイルム株式会社 | Image processing apparatus and image processing method |
-
2016
- 2016-01-07 JP JP2016002063A patent/JP6132938B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2016104196A (en) | 2016-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5868119B2 (en) | Image processing apparatus, radiation imaging system, image processing method, and recording medium | |
JP5921180B2 (en) | Image processing apparatus, image processing method, and program | |
KR101367747B1 (en) | Radiographic apparatus and control method for the same | |
JP5713752B2 (en) | Image processing apparatus and control method thereof | |
JP5719418B2 (en) | High dynamic range image exposure time control method | |
JP6041502B2 (en) | Image processing apparatus and control method | |
KR20140039939A (en) | Photograph image generating method, apparatus therof, and medium storing program source thereof | |
JP2007082181A (en) | Imaging apparatus and image processing method | |
JP6491545B2 (en) | Image processing apparatus, radiation imaging apparatus, image processing method, program, and storage medium | |
JP2014179756A (en) | Image processing device, control method therefor and control program | |
JP2017012445A (en) | Image processing device, image processing method and image processing program | |
JP4863700B2 (en) | X-ray inspection equipment | |
JP5238175B2 (en) | X-ray diagnostic imaging equipment | |
JP6132938B2 (en) | Image processing apparatus, radiation imaging system, image processing method, and recording medium | |
JP6105903B2 (en) | Image processing apparatus, image processing method, radiation imaging system, and program | |
JP5932392B2 (en) | Image processing apparatus and image processing method | |
JP2005149762A (en) | X-ray inspection device, and tube voltage/tube current control method of the same | |
JPWO2015174206A1 (en) | Image diagnostic apparatus and gradation information setting method | |
JP2011229012A (en) | Image processing system, image processing method and program thereof | |
JP5693099B2 (en) | Image processing apparatus, image processing method, radiation imaging system, and program | |
JP2016201634A (en) | Imaging device, imaging method and program | |
JP2012044559A (en) | Imaging apparatus and control method of the same | |
US9959614B2 (en) | Image processing device and method | |
JP2013059511A (en) | Image processing apparatus | |
JP2015093013A (en) | Radiation image processing device, radiographic apparatus, and control methods, and programs thereof |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20161130 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20161206 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20170130 |
|
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: 20170321 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170418 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 6132938 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |