画像処理装置、画像処理方法、画像処理プログラム、および、Ｘ線画像撮影装置
【課題】画像のランダムな場所に発生し、比較的幅が狭く、振幅の大きなスジ状の濃度ムラを除去することができるＸ線画像撮影装置を提供する。
【解決手段】画像処理装置は、撮影された被写体および第１の方向に延びるスジ状の濃度ムラを含む元画像に対して、第１の方向と直交する第２の方向にフィルタ処理を行って高周波成分が抽出された第１の画像を生成する第１のフィルタ処理部と、第１の画像に対して、第２の方向に閾値判定を行って閾値よりも小さい画素値の画素が除去された第２の画像を生成する閾値判定部と、第２の画像に対して、第１の方向にフィルタ処理を行って低周波成分が抽出された第３の画像を生成する第２のフィルタ処理部と、元画像から第３の画像を減算することにより、元画像に含まれるスジ状の濃度ムラが補正された補正後の画像を生成する減算部とを備えている。
【解決手段】画像処理装置は、撮影された被写体および第１の方向に延びるスジ状の濃度ムラを含む元画像に対して、第１の方向と直交する第２の方向にフィルタ処理を行って高周波成分が抽出された第１の画像を生成する第１のフィルタ処理部と、第１の画像に対して、第２の方向に閾値判定を行って閾値よりも小さい画素値の画素が除去された第２の画像を生成する閾値判定部と、第２の画像に対して、第１の方向にフィルタ処理を行って低周波成分が抽出された第３の画像を生成する第２のフィルタ処理部と、元画像から第３の画像を減算することにより、元画像に含まれるスジ状の濃度ムラが補正された補正後の画像を生成する減算部とを備えている。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、撮影された画像のランダムな場所に発生するスジ状（直線状）の濃度ムラ（アーチファクト、ノイズ成分）を除去する画像処理装置、画像処理方法、画像処理プログラム、および、前述の画像処理装置を用いたＸ線画像撮影装置に関するものである。
【背景技術】
【０００２】
被写体（被検者）にＸ線を照射し、被写体を透過したＸ線をＸ線検出器で検出して電気信号に変換し、変換した電気信号に基づいて被写体が撮影されたＸ線画像を生成するＸ線画像撮影装置が知られている。
【０００３】
現在では、Ｘ線検出器として、受光したＸ線を電気信号に変換するフラットパネル型検出器（ＦＰＤ（Flat Panel Detector））が利用されている。ＦＰＤを用いたＸ線画像撮影装置では、Ｘ線源から被写体にＸ線を照射し、被写体を透過したＸ線をＦＰＤで検出して電気信号に変換し、ＦＰＤから被写体の画像データに相当する電気信号を読み出してＸ線画像が生成される。
【０００４】
ＦＰＤは、Ｘ線画像の各画素に対応する複数のＸ線検出素子（光電変換素子）をマトリックス状に配列して構成されている。また、ＦＰＤは、図７に示すように、Ｘ線の受光面が行方向に対して複数のサブエリアに分割され、サブエリア毎にデータ読出回路が設けられている。つまり、各々のサブエリアには、Ｘ線画像の行方向の複数の画素から画像データを読み出すための複数の読出ライン（列）が含まれている。
【０００５】
ＦＰＤから画像データの読み出しが行われる場合、ＦＰＤの各行を選択するための選択信号（ゲート信号）が１行ずつアクティブ状態とされるゲートスキャンが行われる。ゲートスキャンにより、ＦＰＤの行が１行ずつ選択され、１行分の画像データが各々の読出ラインおよびデータ読出回路を介して読み出される。また、全ての列について１行分の画像データの読み出しを繰り返すことで１画面分の画像データの読み出しが行われる。
【０００６】
ところで、例えば、低湿環境において、ＦＰＤのデータ読出回路近辺で発生した静電気等に起因して、撮影されたＸ線画像（ＦＰＤから読み出されたＸ線画像の画像データ）に短いスジ状の濃度ムラが重畳される場合がある。
【０００７】
このスジ状の濃度ムラは、ＦＰＤの、ある行をゲートスキャンした瞬間にのみ発生するなど（１ラインの横スジ）、撮影されたＸ線画像のランダムな場所に発生する。また、このランダムな場所に発生するスジ状の濃度ムラは、例えば、１〜３画素程度の幅の狭い高周波成分のものであるが、振幅（画素値）は非常に大きく、画像上では、被写体の関心領域上で真っ白になったり、真っ黒になるなど非常に目立って視認される。
【０００８】
これに対し、フィルタ処理によりスジ状の濃度ムラを抽出して元画像から減算することにより、元画像からスジ状の濃度ムラを除去する処理が知られている。例えば、特許文献１には、グリッドを使用して撮影を行うことにより得られるグリッド像が記録された画像から、グリッドを除去する特性を有するフィルタを用いてグリッド像に起因するノイズを除去する技術が開示されている。
【０００９】
また、特許文献２には、スジ状の濃度ムラの濃さによって除去する／しないを制御することで、単純なフィルタ処理で、発生するスジ状の濃度ムラを除去する処理が提案されている。この処理では、スジ状の濃度ムラと直交する方向の急激な変化を人体パターンとみなし、閾値に応じて、変化を緩やかにする前処理を施してから、フィルタ処理によりスジ成分を抽出し、これを元画像から引くことでスジ状の濃度ムラを除去している。
【００１０】
なお、特許文献２では、閾値の設定によって、スジ状の濃度ムラの成分と線状画像情報成分との分離がなされ、閾値を超えるスジ状の濃度ムラは除去されない。つまり、閾値を超えるものは人体パターンとみなされている。しかし、特許文献２の閾値処理では、振幅の大きなスジ状の濃度ムラのみを除去することができず、無理に除去しようとして閾値を無限大にすると人体パターン（骨梁等）まで除去されてしまう。
【００１１】
また、特許文献１，２が補正対象とするスジ状の濃度ムラは、例えば、１０画素程度の比較的幅の広い低周波成分のものであり、前述のように、１〜３画素程度の比較的幅の狭い高周波成分のものを検出する場合には適していない。
【００１２】
図８は、被写体とスジ状の濃度ムラについて、画素値と周波数との関係を表す概念図である。同図の縦軸は、被写体の画素値（スジ状の濃度ムラの振幅）、横軸は周波数を表す。同図に示すように、被写体は、その大部分が低周波成分（例えば、１０画素以上の幅があるもの）であり、高周波成分（例えば、１〜３画素程度の幅のもの）は少ないと考えられる。
【００１３】
従って、低周波成分のうち、画素値が大きいものの大部分は被写体であり、スジ状の濃度ムラの画素値は小さいのもが多いと考えられる。一方、高周波成分のうち、画素値が大きいものの大部分はスジ状の濃度ムラであり、被写体の画素値は小さいものが多いと考えられる。つまり、周波数が高くなるに従って、被写体の画素値は小さくなり、スジ状の濃度ムラの画素値は大きくなると考えられる。
【００１４】
特許文献２は、このような前提に基づいて、低周波成分および高周波成分の全域にわたって閾値よりも小さいものをスジ状の濃度ムラ（低周波成分）として検出するものであり、高周波成分で高画素値のスジ状の濃度ムラを除去することは考慮されていない。
【００１５】
また、スジ状の濃度ムラの発生箇所が固定されていれば、これを線欠陥としてあらかじめ抽出しておくことも考えられるが、ランダムな場所に発生する振幅の大きなスジ状の濃度ムラの場合は、撮影画像から直接スジ状の濃度ムラを抽出して除去する必要がある。
【先行技術文献】
【特許文献】
【００１６】
【特許文献１】特開平９−４４６４５号公報
【特許文献２】特開２００５−８４９０２号公報
【発明の概要】
【発明が解決しようとする課題】
【００１７】
本発明の目的は、画像のランダムな場所に発生し、比較的幅が狭く、振幅の大きなスジ状の濃度ムラを除去することができる画像処理装置、画像処理方法、画像処理プログラム、および、前述の画像処理装置を用いたＸ線画像撮影装置を提供することにある。
【課題を解決するための手段】
【００１８】
上記目的を達成するために、本発明は、撮影された被写体および第１の方向に延びるスジ状の濃度ムラを含む元画像に対して、前記第１の方向と直交する第２の方向にフィルタ処理を行って高周波成分が抽出された第１の画像を生成する第１のフィルタ処理部と、
前記第１の画像に対して、前記第２の方向に閾値判定を行って該閾値よりも小さい画素値の画素が除去された第２の画像を生成する閾値判定部と、
前記第２の画像に対して、前記第１の方向にフィルタ処理を行って低周波成分が抽出された第３の画像を生成する第２のフィルタ処理部と、
前記元画像から前記第３の画像を減算することにより、前記元画像に含まれるスジ状の濃度ムラが補正された補正後の画像を生成する減算部とを備えていることを特徴とする画像処理装置を提供するものである。
【００１９】
さらに、前記元画像を前記第１の方向に平均縮小した第４の画像を生成する平均縮小部と、
前記第３の画像を前記第１の方向に平均拡大した第５の画像を生成する平均拡大部とを備え、
前記第１のフィルタ処理部は、前記元画像の代わりに、前記平均縮小部によって生成された第４の画像に対してフィルタ処理を行い、かつ、前記減算部は、前記第３の画像の代わりに、前記平均拡大部によって生成された第５の画像を前記元画像から減算するものであることが好ましい。
【００２０】
さらに、前記第１の画像に含まれる高周波成分の振幅を検出し、検出された高周波成分の振幅のうちの最大振幅の半値を前記閾値に設定する閾値設定部を備えていることが好ましい。
【００２１】
もしくは、さらに、撮影条件を変えて撮影された複数の元画像から、それぞれの元画像に撮影された被写体の画素値およびスジ状の濃度ムラの画素値を検出し、前記被写体およびスジ状の濃度ムラの画素値の実績値に基づいて実験的に決定された、前記被写体およびスジ状の濃度ムラを区別することができる閾値が外部から設定される閾値設定部を備えていることが好ましい。
【００２２】
また、前記閾値判定部は、前記閾値を超える点から前記スジ状の濃度ムラの第１の方向の前後で、所定の遡行画素数を超えない範囲で前記スジ状の濃度ムラが単調に変化する塊の画素は残し、それ以外の画素を除去するものであることが好ましい。
【００２３】
また、前記第１のフィルタ処理部は、入力された画像に対して、メディアンフィルタを用いてフィルタ処理を行って第６の画像を生成し、前記入力された画像から前記第６の画像を減算することにより前記第１の画像を生成するものであることが好ましい。
【００２４】
ここで、前記第１のフィルタ処理部は、前記第２の方向に対するスジ状の濃度ムラの画素数がＮである場合、前記入力された画像に対して、（２Ｎ＋１）画素のメディアンフィルタを用いてフィルタ処理を行うものであることが好ましい。
【００２５】
もしくは、前記第１のフィルタ処理部は、入力された画像に対して、多段移動平均フィルタを用いてフィルタ処理を行って第６の画像を生成し、前記入力された画像から前記第６の画像を減算することにより前記第１の画像を生成するものであることが好ましい。
【００２６】
また、前記第２のフィルタ処理部は、前記第２の画像に対して、メディアンフィルタを用いてフィルタ処理を行うものであることが好ましい。
【００２７】
もしくは、前記第２のフィルタ処理部は、前記第２の画像に対して、多段移動平均フィルタを用いてフィルタ処理を行うものであることが好ましい。
【００２８】
また、本発明は、被写体にＸ線を照射し、照射されたＸ線をＦＰＤで検出して前記被写体のＸ線画像を撮影する撮影部と、
前記撮影部から供給されたＸ線画像を前記元画像として前記スジ状の濃度ムラの補正処理を行う、上記のいずれかに記載の画像処理装置とを備えることを特徴とするＸ線画像撮影装置を提供する。
【００２９】
また、本発明は、第１の方向に延びるスジ状の濃度ムラを含む元画像に対して、前記第１の方向と直交する第２の方向にフィルタ処理を行って高周波成分が抽出された第１の画像を生成し、
前記第１の画像に対して、前記第２の方向に閾値判定を行って該閾値よりも小さい画素が除去された第２の画像を生成し、
前記第２の画像に対して、前記第１の方向にフィルタ処理を行って低周波成分が抽出された第３の画像を生成し、
前記元画像から前記第３の画像を減算して前記スジ状の濃度ムラが補正された補正後の画像を生成することを特徴とする画像処理方法を提供する。
【００３０】
ここで、前記元画像の代わりに、前記元画像を前記第１の方向に平均縮小した第４の画像に対してフィルタ処理を行って前記第１の画像を生成し、前記第３の画像の代わりに、前記第３の画像を前記第１の方向に平均拡大した第５の画像を前記元画像から減算して前記補正後の画像を生成することが好ましい。
【００３１】
また、本発明は、第１の方向に延びるスジ状の濃度ムラを含む元画像に対して、前記第１の方向と直交する第２の方向にフィルタ処理を行って高周波成分が抽出された第１の画像を生成するステップと、
前記第１の画像に対して、前記第２の方向に閾値判定を行って該閾値よりも小さい画素が除去された第２の画像を生成し、
前記第２の画像に対して、前記第１の方向にフィルタ処理を行って低周波成分が抽出された第３の画像を生成するステップと、
前記元画像から前記第３の画像を減算して前記スジ状の濃度ムラが補正された補正後の画像を生成するステップとをコンピュータに実行させることを特徴とする画像処理方法を提供する。
【００３２】
ここで、前記元画像の代わりに、前記元画像を前記第１の方向に平均縮小した第４の画像に対してフィルタ処理を行って前記第１の画像を生成するステップと、前記第３の画像の代わりに、前記第３の画像を前記第１の方向に平均拡大した第５の画像を前記元画像から減算して前記補正後の画像を生成するステップとをコンピュータに実行させることが好ましい。
【発明の効果】
【００３３】
本発明によれば、スジ状の濃度ムラと被写体とを区別するための所定の閾値を用いて、スジ状の濃度ムラと被写体との境界を明確化し、閾値を超えるもののみをスジ状の濃度ムラであると判定することにより、被写体を除去することなく、撮影された画像のランダムな場所に発生し、比較的幅が狭く、振幅の大きなスジ状の濃度ムラだけを除去することができる。
【図面の簡単な説明】
【００３４】
【図１】本発明のＸ線画像撮影装置の構成を表す一実施形態のブロック概念図である。
【図２】図１に示す画像処理部の構成を表すブロック図である。
【図３】（Ａ）〜（Ｃ）は、スジ状の濃度ムラの特性（プロファイル）を表すグラフである。
【図４】（Ａ）〜（Ｃ）は、画像処理部の動作を説明する概念図である。
【図５】（Ｄ）〜（Ｆ）は、図４に続く画像処理部の動作を説明する概念図である。
【図６】（Ｇ）〜（Ｉ）は、図５に続く画像処理部の動作を説明する概念図である。
【図７】ＦＰＤの構成を表す概念図である。
【図８】被写体とスジ状の濃度ムラについて、画素値と周波数との関係を表す概念図である。
【発明を実施するための形態】
【００３５】
以下に、添付の図面に示す好適実施形態に基づいて、本発明の画像処理装置、画像処理方法、画像処理プログラム、および、Ｘ線画像撮影装置を詳細に説明する。
【００３６】
図１は、本発明のＸ線画像撮影装置の構成を表す一実施形態のブロック概念図である。同図に示すＸ線画像撮影装置１０は、被写体ＨにＸ線を照射し、被写体Ｈを透過したＸ線をＰＦＤで検出して電気信号に変換し、変換した電気信号に基づいて、被写体Ｈが撮影されたＸ線画像を生成する。撮影装置１０は、撮影部１２と、画像処理部１４と、出力部１６と、指示入力部１８と、制御部２０とによって構成されている。
【００３７】
指示入力部１８は、撮影開始の指示、画像処理条件、出力指示等の各種の指示を入力するためのものである。指示入力部１８からは、撮影開始の指示、画像処理条件、出力指示等を表す指示信号が出力される。
【００３８】
制御部２０は、指示入力部１８から供給された指示信号に応じて、撮影装置１０の動作を制御する部位である。制御部２０は、例えば、撮影部１２における撮影の制御、画像処理部１４における画像処理の制御、出力部１６における出力の制御等を行う。
【００３９】
撮影部１２は、制御部２０の制御に従って、被写体ＨにＸ線を照射し、被写体Ｈを透過したＸ線をＦＰＤ２８で検出することで被写体Ｈの撮影を行う部位である。撮影部１２からは、被写体Ｈが撮影されたＸ線画像の画像データ（デジタル画像データ）が出力される。
【００４０】
画像処理部１４は、本発明の画像処理装置であり、制御部２０の制御に従って、撮影部１２から供給されたＸ線画像の画像データに対して、オフセット補正、残像補正、欠陥画素補正、本発明の画像処理方法に従うスジ状の濃度ムラ補正等を含む画像処理を行う部位である。画像処理部１４は、例えば、コンピュータ上で動作する画像処理プログラム（ソフトウェア）によって構成される。画像処理部１４からは、画像処理後のＸ線画像の画像データが出力される。
【００４１】
出力部１６は、制御部２０の制御に従って、画像処理部１４から供給された画像処理後のＸ線画像の画像データを出力する部位である。出力部１６は、例えば、Ｘ線画像を画面上に表示するモニタ、Ｘ線画像をプリント出力するプリンタ、Ｘ線画像の画像データを記憶する記憶装置等である。
【００４２】
続いて、撮影部１２について説明する。
【００４３】
撮影部１２は、図１に示すように、Ｘ線源２２と、撮影台２４と、Ｘ線検出部２６とを備えている。
【００４４】
Ｘ線源２２は、撮影時に、設定された強度のＸ線を設定された時間だけ照射（曝射）する。Ｘ線源２２から照射されたＸ線は、撮影台２４上の被写体Ｈを透過してＸ線検出部２６に入射される。
【００４５】
Ｘ線検出部２６は、被写体Ｈを透過したＸ線を検出して電気信号（Ｘ線画像の画像データ）に変換するものであり、ＦＰＤ２８、Ａ／Ｄ（アナログ／デジタル）変換器（図示省略）等を備えている。ＦＰＤ２８からは、被写体Ｈが撮影されたＸ線画像のデータ（アナログデータ）が出力され、これがＡ／Ｄ変換器によりＡ／Ｄ変換されて、変換後のＸ線画像の画像データ（デジタルデータ）が出力される。
【００４６】
続いて、画像処理部１４について説明する。
【００４７】
図２は、図１に示す画像処理部の構成を表すブロック図である。同図に示す画像処理部１４は、Ｘ線画像のランダムな場所に発生し、比較的幅が狭く、振幅の大きなスジ状の濃度ムラを補正する部分のみを表したものである。画像処理部１４は、平均縮小部３０と、第１のフィルタ処理部３２と、閾値設定部３４と、閾値判定部３６と、第２のフィルタ処理部３８と、平均拡大部４０と、減算部４２とを備えている。
【００４８】
撮影部１２から画像処理部１４に供給されたＸ線画像（元画像）は平均縮小部３０に入力される。元画像には、撮影された被写体Ｈおよび第１の方向（画像の横方向）に延びるスジ状の濃度ムラが含まれている。
【００４９】
平均縮小部３０は、元画像を横方向に平均縮小した第１の画像を生成する。平均縮小は、例えば、横方向の２画素を平均し、これを新たな画素の画素値とすることにより、元画像を横方向に１／２のサイズに縮小する処理である。
【００５０】
元画像を横方向に平均縮小することにより、画像サイズを数分の１に縮小することができる。そのため、画像処理部１４における処理の計算量を大幅に削減することができ、処理を高速化することができる。また、元画像を横方向に平均縮小することにより、その縮小率に応じて、スジ状の濃度ムラのうちの横方向の長さが補正対象となる所定の長さを満たさない短いスジをスジ候補から除外することができるという利点もある。
【００５１】
続いて、第１のフィルタ処理部３２は、平均縮小部３０によって生成された第１の画像に対して、第１の方向と直交する第２の方向（画像の縦方向）にフィルタ処理を行って高周波成分が抽出された第２の画像を生成する。
【００５２】
例えば、第１の画像に対して、メディアンフィルタを用いてフィルタ処理を行って第６の画像を生成し、第１の画像から第６の画像を減算することにより第２の画像を生成することができる。なお、縦方向に対するスジ状の濃度ムラの画素数がＮである場合、第１の画像に対して、（２Ｎ＋１）画素のメディアンフィルタを用いてフィルタ処理を行うことが望ましい。
【００５３】
また、第１の画像に対して、移動平均の重みを画素ごとに変える多段移動平均フィルタを用いてフィルタ処理を行って第６の画像を生成し、第１の画像から第６の画像を減算することにより第２の画像を生成することもできる。
【００５４】
ここで、スジ状の濃度ムラは、縦方向（幅方向）の画素数が１〜３画素程度であり、被写体Ｈと比べて高周波成分となる。従って、第１のフィルタ処理部３２により縦方向にフィルタ処理を行って、第１の画像から高周波成分を抽出することにより、スジ状の濃度ムラを抽出することができる。なお、この時同時に、被写体Ｈの縦方向の高周波成分も第２の画像に含まれる。
【００５５】
続いて、閾値設定部３４は、元画像に撮影された被写体Ｈと、スジ状の濃度ムラとを区別するための閾値を設定するための部位である。
【００５６】
閾値設定部３４は、例えば、第１の画像に含まれる高周波成分の振幅を検出し、検出された複数の高周波成分の振幅のうちの最大振幅の半値を閾値として設定する。
【００５７】
また、外部から指示入力部１８介して閾値を閾値設定部３４に設定することもできる。この場合、例えば、撮影条件を変えて撮影された複数の元画像から、それぞれの元画像に撮影された被写体Ｈの画素値と、スジ状の濃度ムラの振幅（画素値）とを検出し、両者の画素値の実績値に基づいて実験的に決定された、両者を区別することができる最適な閾値が閾値設定部３４に設定される。
【００５８】
続いて、閾値判定部３６は、第１のフィルタ処理部３２によって生成された第２の画像に対して、閾値設定部３４から供給された閾値を用いて縦方向に閾値判定を行って閾値よりも小さい画素値の画素が除去された第３の画像を生成する。閾値判定により、第２の画像に含まれる被写体Ｈの縦方向の高周波成分が除去されるとともに、スジ状の濃度ムラのうちの閾値よりも小さい成分もノイズ（スジではない）とみなされ、除去される。
【００５９】
ここで、図３（Ａ）は、スジ状の濃度ムラの特性（プロファイル）を表すグラフである。このグラフの縦軸は、スジ状の濃度ムラの幅方向（縦方向）の振幅、横軸は、幅方向の位置（幅の広さ）を表す。このグラフに示す例のスジ状の濃度ムラは、閾値を超える２つのパルス状の波形になっている。同図では、スジ状の濃度ムラの振幅の所定位置に閾値が設定されている。また、スジ状の濃度ムラの閾値を超えない部分は、画素値＝０（ゼロ）を中心としてプラス側ないしマイナス側にわずかに振れている。
【００６０】
続いて、同図（Ｂ）は、特許文献２における前処理後のスジ状の濃度ムラの特性を表すグラフである。このグラフは、同図（Ａ）に示すスジ状の濃度ムラに対して、特許文献２の前処理を施した後のスジ状の濃度ムラを表したものである。このグラフに示すように、特許文献２の前処理後のスジ状の濃度ムラは、振幅の閾値を超える部分が除去（画素値＝０）されている。ただし、ある箇所で閾値を超えた場合、その画素の幅方向の前後で、ある遡行画素数を超えない範囲でスジ状の濃度ムラが単調に変化する塊全部の画素値が０とされる。
【００６１】
同図（Ｃ）は、本実施形態における閾値判定処理後のスジ状の濃度ムラの特性を表すグラフである。このグラフは、同図（Ａ）に示すスジ状の濃度ムラに対して、本実施形態における閾値判定処理を施した後のスジ状の濃度ムラを表したものである。このグラフに示すように、本実施形態の閾値判定処理後のスジ状の濃度ムラは、振幅の閾値よりも小さい部分が除去（画素値＝０）されている。ただし、ある箇所で閾値を超えた場合、その画素の幅方向の前後で、ある遡行画素数を超えない範囲でスジ状の濃度ムラが単調に変化する塊全部の画素値はそのままの状態で残され、それ以外の部分の画素値が０とされる。
【００６２】
被写体Ｈの部分は、高周波成分の高画素値のスジ状の濃度ムラと比べて振幅が小さいため、閾値を超えるもののみをスジ状の濃度ムラであると判定することにより、被写体Ｈを除去することなく、スジ状の濃度ムラだけを除去することができる。
【００６３】
続いて、第２のフィルタ処理部３８は、閾値判定部３６によって生成された第３の画像に対して、例えば、メディアンフィルタまたは多段移動平均フィルタを用いて第１の方向にフィルタ処理を行って低周波成分が抽出された第４の画像を生成する。
【００６４】
前述の平均縮小では、縮小率に応じて、スジ状の濃度ムラのうちの横方向の長さが補正対象となる所定の長さを満たさないものが除去されるが、第３の画像には、横方向の長さが補正対象となる所定の長さを満たさないものが残されている。第２のフィルタ処理部３８は、フィルタのカットオフ周波数に応じて、横方向の長さが補正対象となる所定の長さに満たないものを除去し、補正対象となる所定の長さ以上のものだけを抽出する。
【００６５】
なお、横方向において、被写体Ｈとスジ状の濃度ムラとの境界が急峻な場合、つまり、両者の境界がはっきりしている場合、多段移動平均フィルタを用いてフィルタ処理を行うと両者の境界がぼけてしまうので、メディアンフィルタを用いてフィルタ処理を行うことが望ましい。また、この場合には、同様の理由から、平均縮小の縮小率を小さくする、もしくは平均縮小を行わないことが望ましい。
【００６６】
続いて、平均拡大部４０は、第２のフィルタ処理部３８によって生成された第４の画像を第１の方向に平均拡大した第５の画像を生成する。平均拡大は、平均縮小によって縮小された元画像に対応して、第４の画像を元画像のサイズに戻す処理であり、例えば、平均縮小によって元画像が横方向に１／２のサイズに縮小された場合、各画素を横方向に２倍することにより元画像のサイズに拡大する。
【００６７】
最後に、減算部４２は、元画像から、平均拡大部４０によって生成された第５の画像、つまり、横方向に延びるスジ状の濃度ムラだけを含む画像を減算することにより、元画像に含まれるスジ状の濃度ムラが補正された補正後の画像を生成する。
【００６８】
次に、本発明の画像処理方法に従って、画像処理部１４において、スジ状の濃度ムラを補正する場合の動作を説明する。
【００６９】
撮影部１２から供給されたＸ線画像（元画像）は平均縮小部３０に入力される。
まず、平均縮小部３０により、図４（Ａ）に示す元画像を横方向に１／２のサイズに平均縮小した第１の画像（同図（Ｂ）参照）が生成される。
【００７０】
続いて、第１のフィルタ処理部３２により、同図（Ｂ）に示す第１の画像に対して、縦方向にフィルタ処理を行って高周波成分が抽出され、横方向のスジ状の濃度ムラおよび被写体Ｈの縦方向の高周波成分を含む第２の画像（同図（Ｃ）参照）が生成される。
【００７１】
続いて、閾値判定部３６により、図４（Ｃ）に示す第２の画像に対して、閾値設定部３４から供給された閾値を用いて縦方向に閾値判定を行って閾値よりも小さい画素値の画素が除去された第３の画像（図５（Ｆ）参照）が生成される。
【００７２】
ここで、図５（Ｄ）に示すように、補正対象とする高周波成分の高画素値のスジ状の濃度ムラの振幅（画素値）は、被写体Ｈの画素値よりも大きい。同図（Ｅ）に示すように、閾値判定により、第２の画像に含まれる被写体Ｈの縦方向の高周波成分が除去されるとともに、スジ状の濃度ムラの閾値よりも小さい成分が除去される（図３（Ｃ）参照）。
【００７３】
続いて、第２のフィルタ処理部３８により、図５（Ｆ）に示す第３の画像に対して、横方向にフィルタ処理を行って低周波成分が抽出され、横方向のスジ状の濃度ムラのうちの補正対象となる所定の長さ以上のものだけを含む第４の画像（図６（Ｇ）参照）が生成される。
【００７４】
続いて、平均拡大部４０により、平均縮小部３０によって１／２のサイズに縮小された元画像に対応して、図６（Ｇ）に示す第４の画像を元画像のサイズに戻すために、第４の画像を横方向に２倍に平均拡大した第５の画像（同図（Ｈ）参照）が生成される。なお、平均縮小および平均拡大は必須ではなく、縮小／拡大をせずに上記の処理を行うようにしてもよい。
【００７５】
最後に、減算部４２により、元画像から同図（Ｈ）に示す第５の画像、つまり、横方向に延びるスジ状の濃度ムラだけを含む画像を減算することにより、元画像に含まれるスジ状の濃度ムラが補正された補正後の画像（同図（Ｉ）参照）が生成される。
【００７６】
上記のように、スジ状の濃度ムラと被写体Ｈとを区別するための所定の閾値を用いて、スジ状の濃度ムラと被写体Ｈとの境界を明確化し、閾値を超えるもののみをスジ状の濃度ムラであると判定することにより、被写体Ｈを除去することなく、撮影されたＸ線画像のランダムな場所に発生し、比較的幅が狭く、振幅の大きなスジ状の濃度ムラだけを除去することができる。
【００７７】
次に、撮影装置１０の動作を説明する。
【００７８】
指示入力部１８において、ユーザにより撮影条件が設定され、撮影指示が与えられると、被写体Ｈの撮影が開始される。撮影が開始されると、撮影部１２において、Ｘ線源２２から、撮影条件に応じて設定された強度のＸ線が設定された時間だけ照射される。照射されたＸ線は、撮影台２４上の被写体Ｈを透過してＸ線検出部２６のＦＰＤ２８に入射され、被写体Ｈを透過したＸ線が電気信号（Ｘ線画像の画像データ）に変換される。
【００７９】
続いて、ＦＰＤ２８から画像データ（アナログデータ）が読み出され、Ａ／Ｄ変換器によりＡ／Ｄ変換が行われ、変換後の画像データ（デジタルデータ）が出力される。画像処理部１４では、変換後のＸ線画像の画像データに対して、横方向のスジ状の濃度ムラの補正を含む画像処理が行われる。画像処理後のＸ線画像の画像データは出力部１６に供給され、例えば、モニタ上に表示されたり、プリント出力されたりする。
【００８０】
なお、本発明は、Ｘ線画像だけでなく、どのような画像に対しても適用可能である。また、本発明は、画像のランダムな場所に発生し、比較的幅が狭く、振幅の大きなスジ状の濃度ムラに対して好適なものであるが、これに限定されるものではない。また、横方向のスジ状の濃度ムラに適用する例を挙げて説明したが、縦方向のスジ状の濃度ムラに対しても同様に適用可能である。
【００８１】
本発明は、基本的に以上のようなものである。
以上、本発明について詳細に説明したが、本発明は上記実施形態に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良や変更をしてもよいのはもちろんである。
【符号の説明】
【００８２】
１０ Ｘ線画像撮影装置
１２ 撮影部
１４ 画像処理部
１６ 出力部
１８ 指示入力部
２０ 制御部
２２ Ｘ線源
２４ 撮影台
２６ Ｘ線検出部
２８ ＦＰＤ
３０ 平均縮小部
３２ 第１のフィルタ処理部
３４ 閾値設定部
３６ 閾値判定部
３８ 第２のフィルタ処理部
４０ 平均拡大部
４２ 減算部
【技術分野】
【０００１】
本発明は、撮影された画像のランダムな場所に発生するスジ状（直線状）の濃度ムラ（アーチファクト、ノイズ成分）を除去する画像処理装置、画像処理方法、画像処理プログラム、および、前述の画像処理装置を用いたＸ線画像撮影装置に関するものである。
【背景技術】
【０００２】
被写体（被検者）にＸ線を照射し、被写体を透過したＸ線をＸ線検出器で検出して電気信号に変換し、変換した電気信号に基づいて被写体が撮影されたＸ線画像を生成するＸ線画像撮影装置が知られている。
【０００３】
現在では、Ｘ線検出器として、受光したＸ線を電気信号に変換するフラットパネル型検出器（ＦＰＤ（Flat Panel Detector））が利用されている。ＦＰＤを用いたＸ線画像撮影装置では、Ｘ線源から被写体にＸ線を照射し、被写体を透過したＸ線をＦＰＤで検出して電気信号に変換し、ＦＰＤから被写体の画像データに相当する電気信号を読み出してＸ線画像が生成される。
【０００４】
ＦＰＤは、Ｘ線画像の各画素に対応する複数のＸ線検出素子（光電変換素子）をマトリックス状に配列して構成されている。また、ＦＰＤは、図７に示すように、Ｘ線の受光面が行方向に対して複数のサブエリアに分割され、サブエリア毎にデータ読出回路が設けられている。つまり、各々のサブエリアには、Ｘ線画像の行方向の複数の画素から画像データを読み出すための複数の読出ライン（列）が含まれている。
【０００５】
ＦＰＤから画像データの読み出しが行われる場合、ＦＰＤの各行を選択するための選択信号（ゲート信号）が１行ずつアクティブ状態とされるゲートスキャンが行われる。ゲートスキャンにより、ＦＰＤの行が１行ずつ選択され、１行分の画像データが各々の読出ラインおよびデータ読出回路を介して読み出される。また、全ての列について１行分の画像データの読み出しを繰り返すことで１画面分の画像データの読み出しが行われる。
【０００６】
ところで、例えば、低湿環境において、ＦＰＤのデータ読出回路近辺で発生した静電気等に起因して、撮影されたＸ線画像（ＦＰＤから読み出されたＸ線画像の画像データ）に短いスジ状の濃度ムラが重畳される場合がある。
【０００７】
このスジ状の濃度ムラは、ＦＰＤの、ある行をゲートスキャンした瞬間にのみ発生するなど（１ラインの横スジ）、撮影されたＸ線画像のランダムな場所に発生する。また、このランダムな場所に発生するスジ状の濃度ムラは、例えば、１〜３画素程度の幅の狭い高周波成分のものであるが、振幅（画素値）は非常に大きく、画像上では、被写体の関心領域上で真っ白になったり、真っ黒になるなど非常に目立って視認される。
【０００８】
これに対し、フィルタ処理によりスジ状の濃度ムラを抽出して元画像から減算することにより、元画像からスジ状の濃度ムラを除去する処理が知られている。例えば、特許文献１には、グリッドを使用して撮影を行うことにより得られるグリッド像が記録された画像から、グリッドを除去する特性を有するフィルタを用いてグリッド像に起因するノイズを除去する技術が開示されている。
【０００９】
また、特許文献２には、スジ状の濃度ムラの濃さによって除去する／しないを制御することで、単純なフィルタ処理で、発生するスジ状の濃度ムラを除去する処理が提案されている。この処理では、スジ状の濃度ムラと直交する方向の急激な変化を人体パターンとみなし、閾値に応じて、変化を緩やかにする前処理を施してから、フィルタ処理によりスジ成分を抽出し、これを元画像から引くことでスジ状の濃度ムラを除去している。
【００１０】
なお、特許文献２では、閾値の設定によって、スジ状の濃度ムラの成分と線状画像情報成分との分離がなされ、閾値を超えるスジ状の濃度ムラは除去されない。つまり、閾値を超えるものは人体パターンとみなされている。しかし、特許文献２の閾値処理では、振幅の大きなスジ状の濃度ムラのみを除去することができず、無理に除去しようとして閾値を無限大にすると人体パターン（骨梁等）まで除去されてしまう。
【００１１】
また、特許文献１，２が補正対象とするスジ状の濃度ムラは、例えば、１０画素程度の比較的幅の広い低周波成分のものであり、前述のように、１〜３画素程度の比較的幅の狭い高周波成分のものを検出する場合には適していない。
【００１２】
図８は、被写体とスジ状の濃度ムラについて、画素値と周波数との関係を表す概念図である。同図の縦軸は、被写体の画素値（スジ状の濃度ムラの振幅）、横軸は周波数を表す。同図に示すように、被写体は、その大部分が低周波成分（例えば、１０画素以上の幅があるもの）であり、高周波成分（例えば、１〜３画素程度の幅のもの）は少ないと考えられる。
【００１３】
従って、低周波成分のうち、画素値が大きいものの大部分は被写体であり、スジ状の濃度ムラの画素値は小さいのもが多いと考えられる。一方、高周波成分のうち、画素値が大きいものの大部分はスジ状の濃度ムラであり、被写体の画素値は小さいものが多いと考えられる。つまり、周波数が高くなるに従って、被写体の画素値は小さくなり、スジ状の濃度ムラの画素値は大きくなると考えられる。
【００１４】
特許文献２は、このような前提に基づいて、低周波成分および高周波成分の全域にわたって閾値よりも小さいものをスジ状の濃度ムラ（低周波成分）として検出するものであり、高周波成分で高画素値のスジ状の濃度ムラを除去することは考慮されていない。
【００１５】
また、スジ状の濃度ムラの発生箇所が固定されていれば、これを線欠陥としてあらかじめ抽出しておくことも考えられるが、ランダムな場所に発生する振幅の大きなスジ状の濃度ムラの場合は、撮影画像から直接スジ状の濃度ムラを抽出して除去する必要がある。
【先行技術文献】
【特許文献】
【００１６】
【特許文献１】特開平９−４４６４５号公報
【特許文献２】特開２００５−８４９０２号公報
【発明の概要】
【発明が解決しようとする課題】
【００１７】
本発明の目的は、画像のランダムな場所に発生し、比較的幅が狭く、振幅の大きなスジ状の濃度ムラを除去することができる画像処理装置、画像処理方法、画像処理プログラム、および、前述の画像処理装置を用いたＸ線画像撮影装置を提供することにある。
【課題を解決するための手段】
【００１８】
上記目的を達成するために、本発明は、撮影された被写体および第１の方向に延びるスジ状の濃度ムラを含む元画像に対して、前記第１の方向と直交する第２の方向にフィルタ処理を行って高周波成分が抽出された第１の画像を生成する第１のフィルタ処理部と、
前記第１の画像に対して、前記第２の方向に閾値判定を行って該閾値よりも小さい画素値の画素が除去された第２の画像を生成する閾値判定部と、
前記第２の画像に対して、前記第１の方向にフィルタ処理を行って低周波成分が抽出された第３の画像を生成する第２のフィルタ処理部と、
前記元画像から前記第３の画像を減算することにより、前記元画像に含まれるスジ状の濃度ムラが補正された補正後の画像を生成する減算部とを備えていることを特徴とする画像処理装置を提供するものである。
【００１９】
さらに、前記元画像を前記第１の方向に平均縮小した第４の画像を生成する平均縮小部と、
前記第３の画像を前記第１の方向に平均拡大した第５の画像を生成する平均拡大部とを備え、
前記第１のフィルタ処理部は、前記元画像の代わりに、前記平均縮小部によって生成された第４の画像に対してフィルタ処理を行い、かつ、前記減算部は、前記第３の画像の代わりに、前記平均拡大部によって生成された第５の画像を前記元画像から減算するものであることが好ましい。
【００２０】
さらに、前記第１の画像に含まれる高周波成分の振幅を検出し、検出された高周波成分の振幅のうちの最大振幅の半値を前記閾値に設定する閾値設定部を備えていることが好ましい。
【００２１】
もしくは、さらに、撮影条件を変えて撮影された複数の元画像から、それぞれの元画像に撮影された被写体の画素値およびスジ状の濃度ムラの画素値を検出し、前記被写体およびスジ状の濃度ムラの画素値の実績値に基づいて実験的に決定された、前記被写体およびスジ状の濃度ムラを区別することができる閾値が外部から設定される閾値設定部を備えていることが好ましい。
【００２２】
また、前記閾値判定部は、前記閾値を超える点から前記スジ状の濃度ムラの第１の方向の前後で、所定の遡行画素数を超えない範囲で前記スジ状の濃度ムラが単調に変化する塊の画素は残し、それ以外の画素を除去するものであることが好ましい。
【００２３】
また、前記第１のフィルタ処理部は、入力された画像に対して、メディアンフィルタを用いてフィルタ処理を行って第６の画像を生成し、前記入力された画像から前記第６の画像を減算することにより前記第１の画像を生成するものであることが好ましい。
【００２４】
ここで、前記第１のフィルタ処理部は、前記第２の方向に対するスジ状の濃度ムラの画素数がＮである場合、前記入力された画像に対して、（２Ｎ＋１）画素のメディアンフィルタを用いてフィルタ処理を行うものであることが好ましい。
【００２５】
もしくは、前記第１のフィルタ処理部は、入力された画像に対して、多段移動平均フィルタを用いてフィルタ処理を行って第６の画像を生成し、前記入力された画像から前記第６の画像を減算することにより前記第１の画像を生成するものであることが好ましい。
【００２６】
また、前記第２のフィルタ処理部は、前記第２の画像に対して、メディアンフィルタを用いてフィルタ処理を行うものであることが好ましい。
【００２７】
もしくは、前記第２のフィルタ処理部は、前記第２の画像に対して、多段移動平均フィルタを用いてフィルタ処理を行うものであることが好ましい。
【００２８】
また、本発明は、被写体にＸ線を照射し、照射されたＸ線をＦＰＤで検出して前記被写体のＸ線画像を撮影する撮影部と、
前記撮影部から供給されたＸ線画像を前記元画像として前記スジ状の濃度ムラの補正処理を行う、上記のいずれかに記載の画像処理装置とを備えることを特徴とするＸ線画像撮影装置を提供する。
【００２９】
また、本発明は、第１の方向に延びるスジ状の濃度ムラを含む元画像に対して、前記第１の方向と直交する第２の方向にフィルタ処理を行って高周波成分が抽出された第１の画像を生成し、
前記第１の画像に対して、前記第２の方向に閾値判定を行って該閾値よりも小さい画素が除去された第２の画像を生成し、
前記第２の画像に対して、前記第１の方向にフィルタ処理を行って低周波成分が抽出された第３の画像を生成し、
前記元画像から前記第３の画像を減算して前記スジ状の濃度ムラが補正された補正後の画像を生成することを特徴とする画像処理方法を提供する。
【００３０】
ここで、前記元画像の代わりに、前記元画像を前記第１の方向に平均縮小した第４の画像に対してフィルタ処理を行って前記第１の画像を生成し、前記第３の画像の代わりに、前記第３の画像を前記第１の方向に平均拡大した第５の画像を前記元画像から減算して前記補正後の画像を生成することが好ましい。
【００３１】
また、本発明は、第１の方向に延びるスジ状の濃度ムラを含む元画像に対して、前記第１の方向と直交する第２の方向にフィルタ処理を行って高周波成分が抽出された第１の画像を生成するステップと、
前記第１の画像に対して、前記第２の方向に閾値判定を行って該閾値よりも小さい画素が除去された第２の画像を生成し、
前記第２の画像に対して、前記第１の方向にフィルタ処理を行って低周波成分が抽出された第３の画像を生成するステップと、
前記元画像から前記第３の画像を減算して前記スジ状の濃度ムラが補正された補正後の画像を生成するステップとをコンピュータに実行させることを特徴とする画像処理方法を提供する。
【００３２】
ここで、前記元画像の代わりに、前記元画像を前記第１の方向に平均縮小した第４の画像に対してフィルタ処理を行って前記第１の画像を生成するステップと、前記第３の画像の代わりに、前記第３の画像を前記第１の方向に平均拡大した第５の画像を前記元画像から減算して前記補正後の画像を生成するステップとをコンピュータに実行させることが好ましい。
【発明の効果】
【００３３】
本発明によれば、スジ状の濃度ムラと被写体とを区別するための所定の閾値を用いて、スジ状の濃度ムラと被写体との境界を明確化し、閾値を超えるもののみをスジ状の濃度ムラであると判定することにより、被写体を除去することなく、撮影された画像のランダムな場所に発生し、比較的幅が狭く、振幅の大きなスジ状の濃度ムラだけを除去することができる。
【図面の簡単な説明】
【００３４】
【図１】本発明のＸ線画像撮影装置の構成を表す一実施形態のブロック概念図である。
【図２】図１に示す画像処理部の構成を表すブロック図である。
【図３】（Ａ）〜（Ｃ）は、スジ状の濃度ムラの特性（プロファイル）を表すグラフである。
【図４】（Ａ）〜（Ｃ）は、画像処理部の動作を説明する概念図である。
【図５】（Ｄ）〜（Ｆ）は、図４に続く画像処理部の動作を説明する概念図である。
【図６】（Ｇ）〜（Ｉ）は、図５に続く画像処理部の動作を説明する概念図である。
【図７】ＦＰＤの構成を表す概念図である。
【図８】被写体とスジ状の濃度ムラについて、画素値と周波数との関係を表す概念図である。
【発明を実施するための形態】
【００３５】
以下に、添付の図面に示す好適実施形態に基づいて、本発明の画像処理装置、画像処理方法、画像処理プログラム、および、Ｘ線画像撮影装置を詳細に説明する。
【００３６】
図１は、本発明のＸ線画像撮影装置の構成を表す一実施形態のブロック概念図である。同図に示すＸ線画像撮影装置１０は、被写体ＨにＸ線を照射し、被写体Ｈを透過したＸ線をＰＦＤで検出して電気信号に変換し、変換した電気信号に基づいて、被写体Ｈが撮影されたＸ線画像を生成する。撮影装置１０は、撮影部１２と、画像処理部１４と、出力部１６と、指示入力部１８と、制御部２０とによって構成されている。
【００３７】
指示入力部１８は、撮影開始の指示、画像処理条件、出力指示等の各種の指示を入力するためのものである。指示入力部１８からは、撮影開始の指示、画像処理条件、出力指示等を表す指示信号が出力される。
【００３８】
制御部２０は、指示入力部１８から供給された指示信号に応じて、撮影装置１０の動作を制御する部位である。制御部２０は、例えば、撮影部１２における撮影の制御、画像処理部１４における画像処理の制御、出力部１６における出力の制御等を行う。
【００３９】
撮影部１２は、制御部２０の制御に従って、被写体ＨにＸ線を照射し、被写体Ｈを透過したＸ線をＦＰＤ２８で検出することで被写体Ｈの撮影を行う部位である。撮影部１２からは、被写体Ｈが撮影されたＸ線画像の画像データ（デジタル画像データ）が出力される。
【００４０】
画像処理部１４は、本発明の画像処理装置であり、制御部２０の制御に従って、撮影部１２から供給されたＸ線画像の画像データに対して、オフセット補正、残像補正、欠陥画素補正、本発明の画像処理方法に従うスジ状の濃度ムラ補正等を含む画像処理を行う部位である。画像処理部１４は、例えば、コンピュータ上で動作する画像処理プログラム（ソフトウェア）によって構成される。画像処理部１４からは、画像処理後のＸ線画像の画像データが出力される。
【００４１】
出力部１６は、制御部２０の制御に従って、画像処理部１４から供給された画像処理後のＸ線画像の画像データを出力する部位である。出力部１６は、例えば、Ｘ線画像を画面上に表示するモニタ、Ｘ線画像をプリント出力するプリンタ、Ｘ線画像の画像データを記憶する記憶装置等である。
【００４２】
続いて、撮影部１２について説明する。
【００４３】
撮影部１２は、図１に示すように、Ｘ線源２２と、撮影台２４と、Ｘ線検出部２６とを備えている。
【００４４】
Ｘ線源２２は、撮影時に、設定された強度のＸ線を設定された時間だけ照射（曝射）する。Ｘ線源２２から照射されたＸ線は、撮影台２４上の被写体Ｈを透過してＸ線検出部２６に入射される。
【００４５】
Ｘ線検出部２６は、被写体Ｈを透過したＸ線を検出して電気信号（Ｘ線画像の画像データ）に変換するものであり、ＦＰＤ２８、Ａ／Ｄ（アナログ／デジタル）変換器（図示省略）等を備えている。ＦＰＤ２８からは、被写体Ｈが撮影されたＸ線画像のデータ（アナログデータ）が出力され、これがＡ／Ｄ変換器によりＡ／Ｄ変換されて、変換後のＸ線画像の画像データ（デジタルデータ）が出力される。
【００４６】
続いて、画像処理部１４について説明する。
【００４７】
図２は、図１に示す画像処理部の構成を表すブロック図である。同図に示す画像処理部１４は、Ｘ線画像のランダムな場所に発生し、比較的幅が狭く、振幅の大きなスジ状の濃度ムラを補正する部分のみを表したものである。画像処理部１４は、平均縮小部３０と、第１のフィルタ処理部３２と、閾値設定部３４と、閾値判定部３６と、第２のフィルタ処理部３８と、平均拡大部４０と、減算部４２とを備えている。
【００４８】
撮影部１２から画像処理部１４に供給されたＸ線画像（元画像）は平均縮小部３０に入力される。元画像には、撮影された被写体Ｈおよび第１の方向（画像の横方向）に延びるスジ状の濃度ムラが含まれている。
【００４９】
平均縮小部３０は、元画像を横方向に平均縮小した第１の画像を生成する。平均縮小は、例えば、横方向の２画素を平均し、これを新たな画素の画素値とすることにより、元画像を横方向に１／２のサイズに縮小する処理である。
【００５０】
元画像を横方向に平均縮小することにより、画像サイズを数分の１に縮小することができる。そのため、画像処理部１４における処理の計算量を大幅に削減することができ、処理を高速化することができる。また、元画像を横方向に平均縮小することにより、その縮小率に応じて、スジ状の濃度ムラのうちの横方向の長さが補正対象となる所定の長さを満たさない短いスジをスジ候補から除外することができるという利点もある。
【００５１】
続いて、第１のフィルタ処理部３２は、平均縮小部３０によって生成された第１の画像に対して、第１の方向と直交する第２の方向（画像の縦方向）にフィルタ処理を行って高周波成分が抽出された第２の画像を生成する。
【００５２】
例えば、第１の画像に対して、メディアンフィルタを用いてフィルタ処理を行って第６の画像を生成し、第１の画像から第６の画像を減算することにより第２の画像を生成することができる。なお、縦方向に対するスジ状の濃度ムラの画素数がＮである場合、第１の画像に対して、（２Ｎ＋１）画素のメディアンフィルタを用いてフィルタ処理を行うことが望ましい。
【００５３】
また、第１の画像に対して、移動平均の重みを画素ごとに変える多段移動平均フィルタを用いてフィルタ処理を行って第６の画像を生成し、第１の画像から第６の画像を減算することにより第２の画像を生成することもできる。
【００５４】
ここで、スジ状の濃度ムラは、縦方向（幅方向）の画素数が１〜３画素程度であり、被写体Ｈと比べて高周波成分となる。従って、第１のフィルタ処理部３２により縦方向にフィルタ処理を行って、第１の画像から高周波成分を抽出することにより、スジ状の濃度ムラを抽出することができる。なお、この時同時に、被写体Ｈの縦方向の高周波成分も第２の画像に含まれる。
【００５５】
続いて、閾値設定部３４は、元画像に撮影された被写体Ｈと、スジ状の濃度ムラとを区別するための閾値を設定するための部位である。
【００５６】
閾値設定部３４は、例えば、第１の画像に含まれる高周波成分の振幅を検出し、検出された複数の高周波成分の振幅のうちの最大振幅の半値を閾値として設定する。
【００５７】
また、外部から指示入力部１８介して閾値を閾値設定部３４に設定することもできる。この場合、例えば、撮影条件を変えて撮影された複数の元画像から、それぞれの元画像に撮影された被写体Ｈの画素値と、スジ状の濃度ムラの振幅（画素値）とを検出し、両者の画素値の実績値に基づいて実験的に決定された、両者を区別することができる最適な閾値が閾値設定部３４に設定される。
【００５８】
続いて、閾値判定部３６は、第１のフィルタ処理部３２によって生成された第２の画像に対して、閾値設定部３４から供給された閾値を用いて縦方向に閾値判定を行って閾値よりも小さい画素値の画素が除去された第３の画像を生成する。閾値判定により、第２の画像に含まれる被写体Ｈの縦方向の高周波成分が除去されるとともに、スジ状の濃度ムラのうちの閾値よりも小さい成分もノイズ（スジではない）とみなされ、除去される。
【００５９】
ここで、図３（Ａ）は、スジ状の濃度ムラの特性（プロファイル）を表すグラフである。このグラフの縦軸は、スジ状の濃度ムラの幅方向（縦方向）の振幅、横軸は、幅方向の位置（幅の広さ）を表す。このグラフに示す例のスジ状の濃度ムラは、閾値を超える２つのパルス状の波形になっている。同図では、スジ状の濃度ムラの振幅の所定位置に閾値が設定されている。また、スジ状の濃度ムラの閾値を超えない部分は、画素値＝０（ゼロ）を中心としてプラス側ないしマイナス側にわずかに振れている。
【００６０】
続いて、同図（Ｂ）は、特許文献２における前処理後のスジ状の濃度ムラの特性を表すグラフである。このグラフは、同図（Ａ）に示すスジ状の濃度ムラに対して、特許文献２の前処理を施した後のスジ状の濃度ムラを表したものである。このグラフに示すように、特許文献２の前処理後のスジ状の濃度ムラは、振幅の閾値を超える部分が除去（画素値＝０）されている。ただし、ある箇所で閾値を超えた場合、その画素の幅方向の前後で、ある遡行画素数を超えない範囲でスジ状の濃度ムラが単調に変化する塊全部の画素値が０とされる。
【００６１】
同図（Ｃ）は、本実施形態における閾値判定処理後のスジ状の濃度ムラの特性を表すグラフである。このグラフは、同図（Ａ）に示すスジ状の濃度ムラに対して、本実施形態における閾値判定処理を施した後のスジ状の濃度ムラを表したものである。このグラフに示すように、本実施形態の閾値判定処理後のスジ状の濃度ムラは、振幅の閾値よりも小さい部分が除去（画素値＝０）されている。ただし、ある箇所で閾値を超えた場合、その画素の幅方向の前後で、ある遡行画素数を超えない範囲でスジ状の濃度ムラが単調に変化する塊全部の画素値はそのままの状態で残され、それ以外の部分の画素値が０とされる。
【００６２】
被写体Ｈの部分は、高周波成分の高画素値のスジ状の濃度ムラと比べて振幅が小さいため、閾値を超えるもののみをスジ状の濃度ムラであると判定することにより、被写体Ｈを除去することなく、スジ状の濃度ムラだけを除去することができる。
【００６３】
続いて、第２のフィルタ処理部３８は、閾値判定部３６によって生成された第３の画像に対して、例えば、メディアンフィルタまたは多段移動平均フィルタを用いて第１の方向にフィルタ処理を行って低周波成分が抽出された第４の画像を生成する。
【００６４】
前述の平均縮小では、縮小率に応じて、スジ状の濃度ムラのうちの横方向の長さが補正対象となる所定の長さを満たさないものが除去されるが、第３の画像には、横方向の長さが補正対象となる所定の長さを満たさないものが残されている。第２のフィルタ処理部３８は、フィルタのカットオフ周波数に応じて、横方向の長さが補正対象となる所定の長さに満たないものを除去し、補正対象となる所定の長さ以上のものだけを抽出する。
【００６５】
なお、横方向において、被写体Ｈとスジ状の濃度ムラとの境界が急峻な場合、つまり、両者の境界がはっきりしている場合、多段移動平均フィルタを用いてフィルタ処理を行うと両者の境界がぼけてしまうので、メディアンフィルタを用いてフィルタ処理を行うことが望ましい。また、この場合には、同様の理由から、平均縮小の縮小率を小さくする、もしくは平均縮小を行わないことが望ましい。
【００６６】
続いて、平均拡大部４０は、第２のフィルタ処理部３８によって生成された第４の画像を第１の方向に平均拡大した第５の画像を生成する。平均拡大は、平均縮小によって縮小された元画像に対応して、第４の画像を元画像のサイズに戻す処理であり、例えば、平均縮小によって元画像が横方向に１／２のサイズに縮小された場合、各画素を横方向に２倍することにより元画像のサイズに拡大する。
【００６７】
最後に、減算部４２は、元画像から、平均拡大部４０によって生成された第５の画像、つまり、横方向に延びるスジ状の濃度ムラだけを含む画像を減算することにより、元画像に含まれるスジ状の濃度ムラが補正された補正後の画像を生成する。
【００６８】
次に、本発明の画像処理方法に従って、画像処理部１４において、スジ状の濃度ムラを補正する場合の動作を説明する。
【００６９】
撮影部１２から供給されたＸ線画像（元画像）は平均縮小部３０に入力される。
まず、平均縮小部３０により、図４（Ａ）に示す元画像を横方向に１／２のサイズに平均縮小した第１の画像（同図（Ｂ）参照）が生成される。
【００７０】
続いて、第１のフィルタ処理部３２により、同図（Ｂ）に示す第１の画像に対して、縦方向にフィルタ処理を行って高周波成分が抽出され、横方向のスジ状の濃度ムラおよび被写体Ｈの縦方向の高周波成分を含む第２の画像（同図（Ｃ）参照）が生成される。
【００７１】
続いて、閾値判定部３６により、図４（Ｃ）に示す第２の画像に対して、閾値設定部３４から供給された閾値を用いて縦方向に閾値判定を行って閾値よりも小さい画素値の画素が除去された第３の画像（図５（Ｆ）参照）が生成される。
【００７２】
ここで、図５（Ｄ）に示すように、補正対象とする高周波成分の高画素値のスジ状の濃度ムラの振幅（画素値）は、被写体Ｈの画素値よりも大きい。同図（Ｅ）に示すように、閾値判定により、第２の画像に含まれる被写体Ｈの縦方向の高周波成分が除去されるとともに、スジ状の濃度ムラの閾値よりも小さい成分が除去される（図３（Ｃ）参照）。
【００７３】
続いて、第２のフィルタ処理部３８により、図５（Ｆ）に示す第３の画像に対して、横方向にフィルタ処理を行って低周波成分が抽出され、横方向のスジ状の濃度ムラのうちの補正対象となる所定の長さ以上のものだけを含む第４の画像（図６（Ｇ）参照）が生成される。
【００７４】
続いて、平均拡大部４０により、平均縮小部３０によって１／２のサイズに縮小された元画像に対応して、図６（Ｇ）に示す第４の画像を元画像のサイズに戻すために、第４の画像を横方向に２倍に平均拡大した第５の画像（同図（Ｈ）参照）が生成される。なお、平均縮小および平均拡大は必須ではなく、縮小／拡大をせずに上記の処理を行うようにしてもよい。
【００７５】
最後に、減算部４２により、元画像から同図（Ｈ）に示す第５の画像、つまり、横方向に延びるスジ状の濃度ムラだけを含む画像を減算することにより、元画像に含まれるスジ状の濃度ムラが補正された補正後の画像（同図（Ｉ）参照）が生成される。
【００７６】
上記のように、スジ状の濃度ムラと被写体Ｈとを区別するための所定の閾値を用いて、スジ状の濃度ムラと被写体Ｈとの境界を明確化し、閾値を超えるもののみをスジ状の濃度ムラであると判定することにより、被写体Ｈを除去することなく、撮影されたＸ線画像のランダムな場所に発生し、比較的幅が狭く、振幅の大きなスジ状の濃度ムラだけを除去することができる。
【００７７】
次に、撮影装置１０の動作を説明する。
【００７８】
指示入力部１８において、ユーザにより撮影条件が設定され、撮影指示が与えられると、被写体Ｈの撮影が開始される。撮影が開始されると、撮影部１２において、Ｘ線源２２から、撮影条件に応じて設定された強度のＸ線が設定された時間だけ照射される。照射されたＸ線は、撮影台２４上の被写体Ｈを透過してＸ線検出部２６のＦＰＤ２８に入射され、被写体Ｈを透過したＸ線が電気信号（Ｘ線画像の画像データ）に変換される。
【００７９】
続いて、ＦＰＤ２８から画像データ（アナログデータ）が読み出され、Ａ／Ｄ変換器によりＡ／Ｄ変換が行われ、変換後の画像データ（デジタルデータ）が出力される。画像処理部１４では、変換後のＸ線画像の画像データに対して、横方向のスジ状の濃度ムラの補正を含む画像処理が行われる。画像処理後のＸ線画像の画像データは出力部１６に供給され、例えば、モニタ上に表示されたり、プリント出力されたりする。
【００８０】
なお、本発明は、Ｘ線画像だけでなく、どのような画像に対しても適用可能である。また、本発明は、画像のランダムな場所に発生し、比較的幅が狭く、振幅の大きなスジ状の濃度ムラに対して好適なものであるが、これに限定されるものではない。また、横方向のスジ状の濃度ムラに適用する例を挙げて説明したが、縦方向のスジ状の濃度ムラに対しても同様に適用可能である。
【００８１】
本発明は、基本的に以上のようなものである。
以上、本発明について詳細に説明したが、本発明は上記実施形態に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良や変更をしてもよいのはもちろんである。
【符号の説明】
【００８２】
１０ Ｘ線画像撮影装置
１２ 撮影部
１４ 画像処理部
１６ 出力部
１８ 指示入力部
２０ 制御部
２２ Ｘ線源
２４ 撮影台
２６ Ｘ線検出部
２８ ＦＰＤ
３０ 平均縮小部
３２ 第１のフィルタ処理部
３４ 閾値設定部
３６ 閾値判定部
３８ 第２のフィルタ処理部
４０ 平均拡大部
４２ 減算部
【特許請求の範囲】
【請求項１】
撮影された被写体および第１の方向に延びるスジ状の濃度ムラを含む元画像に対して、前記第１の方向と直交する第２の方向にフィルタ処理を行って高周波成分が抽出された第１の画像を生成する第１のフィルタ処理部と、
前記第１の画像に対して、前記第２の方向に閾値判定を行って該閾値よりも小さい画素値の画素が除去された第２の画像を生成する閾値判定部と、
前記第２の画像に対して、前記第１の方向にフィルタ処理を行って低周波成分が抽出された第３の画像を生成する第２のフィルタ処理部と、
前記元画像から前記第３の画像を減算することにより、前記元画像に含まれるスジ状の濃度ムラが補正された補正後の画像を生成する減算部とを備えていることを特徴とする画像処理装置。
【請求項２】
さらに、前記元画像を前記第１の方向に平均縮小した第４の画像を生成する平均縮小部と、
前記第３の画像を前記第１の方向に平均拡大した第５の画像を生成する平均拡大部とを備え、
前記第１のフィルタ処理部は、前記元画像の代わりに、前記平均縮小部によって生成された第４の画像に対してフィルタ処理を行い、かつ、前記減算部は、前記第３の画像の代わりに、前記平均拡大部によって生成された第５の画像を前記元画像から減算するものであることを特徴とする請求項１に記載の画像処理装置。
【請求項３】
さらに、前記第１の画像に含まれる高周波成分の振幅を検出し、検出された高周波成分の振幅のうちの最大振幅の半値を前記閾値に設定する閾値設定部を備えていることを特徴とする請求項１または２に記載の画像処理装置。
【請求項４】
さらに、撮影条件を変えて撮影された複数の元画像から、それぞれの元画像に撮影された被写体の画素値およびスジ状の濃度ムラの画素値を検出し、前記被写体およびスジ状の濃度ムラの画素値の実績値に基づいて実験的に決定された、前記被写体およびスジ状の濃度ムラを区別することができる閾値が外部から設定される閾値設定部を備えていることを特徴とする請求項１または２に記載の画像処理装置。
【請求項５】
前記閾値判定部は、前記閾値を超える点から前記スジ状の濃度ムラの第１の方向の前後で、所定の遡行画素数を超えない範囲で前記スジ状の濃度ムラが単調に変化する塊の画素は残し、それ以外の画素を除去するものであることを特徴とする請求項１〜４のいずれかに記載の画像処理装置。
【請求項６】
前記第１のフィルタ処理部は、入力された画像に対して、メディアンフィルタを用いてフィルタ処理を行って第６の画像を生成し、前記入力された画像から前記第６の画像を減算することにより前記第１の画像を生成するものであることを特徴とする請求項１〜５のいずれかに記載の画像処理装置。
【請求項７】
前記第１のフィルタ処理部は、前記第２の方向に対するスジ状の濃度ムラの画素数がＮである場合、前記入力された画像に対して、（２Ｎ＋１）画素のメディアンフィルタを用いてフィルタ処理を行うものであることを特徴とする請求項６に記載の画像処理装置。
【請求項８】
前記第１のフィルタ処理部は、入力された画像に対して、多段移動平均フィルタを用いてフィルタ処理を行って第６の画像を生成し、前記入力された画像から前記第６の画像を減算することにより前記第１の画像を生成するものであることを特徴とする請求項１〜５のいずれかに記載の画像処理装置。
【請求項９】
前記第２のフィルタ処理部は、前記第２の画像に対して、メディアンフィルタを用いてフィルタ処理を行うものであることを特徴とする請求項１〜８のいずれかに記載の画像処理装置。
【請求項１０】
前記第２のフィルタ処理部は、前記第２の画像に対して、多段移動平均フィルタを用いてフィルタ処理を行うものであることを特徴とする請求項１〜８のいずれかに記載の画像処理装置。
【請求項１１】
被写体にＸ線を照射し、照射されたＸ線をＦＰＤで検出して前記被写体のＸ線画像を撮影する撮影部と、
前記撮影部から供給されたＸ線画像を前記元画像として前記スジ状の濃度ムラの補正処理を行う、請求項１〜１０のいずれかに記載の画像処理装置とを備えることを特徴とするＸ線画像撮影装置。
【請求項１２】
第１の方向に延びるスジ状の濃度ムラを含む元画像に対して、前記第１の方向と直交する第２の方向にフィルタ処理を行って高周波成分が抽出された第１の画像を生成し、
前記第１の画像に対して、前記第２の方向に閾値判定を行って該閾値よりも小さい画素が除去された第２の画像を生成し、
前記第２の画像に対して、前記第１の方向にフィルタ処理を行って低周波成分が抽出された第３の画像を生成し、
前記元画像から前記第３の画像を減算して前記スジ状の濃度ムラが補正された補正後の画像を生成することを特徴とする画像処理方法。
【請求項１３】
前記元画像の代わりに、前記元画像を前記第１の方向に平均縮小した第４の画像に対してフィルタ処理を行って前記第１の画像を生成し、前記第３の画像の代わりに、前記第３の画像を前記第１の方向に平均拡大した第５の画像を前記元画像から減算して前記補正後の画像を生成することを特徴とする請求項１２に記載の画像処理方法。
【請求項１４】
第１の方向に延びるスジ状の濃度ムラを含む元画像に対して、前記第１の方向と直交する第２の方向にフィルタ処理を行って高周波成分が抽出された第１の画像を生成するステップと、
前記第１の画像に対して、前記第２の方向に閾値判定を行って該閾値よりも小さい画素が除去された第２の画像を生成し、
前記第２の画像に対して、前記第１の方向にフィルタ処理を行って低周波成分が抽出された第３の画像を生成するステップと、
前記元画像から前記第３の画像を減算して前記スジ状の濃度ムラが補正された補正後の画像を生成するステップとをコンピュータに実行させることを特徴とする画像処理プログラム。
【請求項１５】
前記元画像の代わりに、前記元画像を前記第１の方向に平均縮小した第４の画像に対してフィルタ処理を行って前記第１の画像を生成するステップと、前記第３の画像の代わりに、前記第３の画像を前記第１の方向に平均拡大した第５の画像を前記元画像から減算して前記補正後の画像を生成するステップとをコンピュータに実行させることを特徴とする請求項１４に記載の画像処理プログラム。
【請求項１】
撮影された被写体および第１の方向に延びるスジ状の濃度ムラを含む元画像に対して、前記第１の方向と直交する第２の方向にフィルタ処理を行って高周波成分が抽出された第１の画像を生成する第１のフィルタ処理部と、
前記第１の画像に対して、前記第２の方向に閾値判定を行って該閾値よりも小さい画素値の画素が除去された第２の画像を生成する閾値判定部と、
前記第２の画像に対して、前記第１の方向にフィルタ処理を行って低周波成分が抽出された第３の画像を生成する第２のフィルタ処理部と、
前記元画像から前記第３の画像を減算することにより、前記元画像に含まれるスジ状の濃度ムラが補正された補正後の画像を生成する減算部とを備えていることを特徴とする画像処理装置。
【請求項２】
さらに、前記元画像を前記第１の方向に平均縮小した第４の画像を生成する平均縮小部と、
前記第３の画像を前記第１の方向に平均拡大した第５の画像を生成する平均拡大部とを備え、
前記第１のフィルタ処理部は、前記元画像の代わりに、前記平均縮小部によって生成された第４の画像に対してフィルタ処理を行い、かつ、前記減算部は、前記第３の画像の代わりに、前記平均拡大部によって生成された第５の画像を前記元画像から減算するものであることを特徴とする請求項１に記載の画像処理装置。
【請求項３】
さらに、前記第１の画像に含まれる高周波成分の振幅を検出し、検出された高周波成分の振幅のうちの最大振幅の半値を前記閾値に設定する閾値設定部を備えていることを特徴とする請求項１または２に記載の画像処理装置。
【請求項４】
さらに、撮影条件を変えて撮影された複数の元画像から、それぞれの元画像に撮影された被写体の画素値およびスジ状の濃度ムラの画素値を検出し、前記被写体およびスジ状の濃度ムラの画素値の実績値に基づいて実験的に決定された、前記被写体およびスジ状の濃度ムラを区別することができる閾値が外部から設定される閾値設定部を備えていることを特徴とする請求項１または２に記載の画像処理装置。
【請求項５】
前記閾値判定部は、前記閾値を超える点から前記スジ状の濃度ムラの第１の方向の前後で、所定の遡行画素数を超えない範囲で前記スジ状の濃度ムラが単調に変化する塊の画素は残し、それ以外の画素を除去するものであることを特徴とする請求項１〜４のいずれかに記載の画像処理装置。
【請求項６】
前記第１のフィルタ処理部は、入力された画像に対して、メディアンフィルタを用いてフィルタ処理を行って第６の画像を生成し、前記入力された画像から前記第６の画像を減算することにより前記第１の画像を生成するものであることを特徴とする請求項１〜５のいずれかに記載の画像処理装置。
【請求項７】
前記第１のフィルタ処理部は、前記第２の方向に対するスジ状の濃度ムラの画素数がＮである場合、前記入力された画像に対して、（２Ｎ＋１）画素のメディアンフィルタを用いてフィルタ処理を行うものであることを特徴とする請求項６に記載の画像処理装置。
【請求項８】
前記第１のフィルタ処理部は、入力された画像に対して、多段移動平均フィルタを用いてフィルタ処理を行って第６の画像を生成し、前記入力された画像から前記第６の画像を減算することにより前記第１の画像を生成するものであることを特徴とする請求項１〜５のいずれかに記載の画像処理装置。
【請求項９】
前記第２のフィルタ処理部は、前記第２の画像に対して、メディアンフィルタを用いてフィルタ処理を行うものであることを特徴とする請求項１〜８のいずれかに記載の画像処理装置。
【請求項１０】
前記第２のフィルタ処理部は、前記第２の画像に対して、多段移動平均フィルタを用いてフィルタ処理を行うものであることを特徴とする請求項１〜８のいずれかに記載の画像処理装置。
【請求項１１】
被写体にＸ線を照射し、照射されたＸ線をＦＰＤで検出して前記被写体のＸ線画像を撮影する撮影部と、
前記撮影部から供給されたＸ線画像を前記元画像として前記スジ状の濃度ムラの補正処理を行う、請求項１〜１０のいずれかに記載の画像処理装置とを備えることを特徴とするＸ線画像撮影装置。
【請求項１２】
第１の方向に延びるスジ状の濃度ムラを含む元画像に対して、前記第１の方向と直交する第２の方向にフィルタ処理を行って高周波成分が抽出された第１の画像を生成し、
前記第１の画像に対して、前記第２の方向に閾値判定を行って該閾値よりも小さい画素が除去された第２の画像を生成し、
前記第２の画像に対して、前記第１の方向にフィルタ処理を行って低周波成分が抽出された第３の画像を生成し、
前記元画像から前記第３の画像を減算して前記スジ状の濃度ムラが補正された補正後の画像を生成することを特徴とする画像処理方法。
【請求項１３】
前記元画像の代わりに、前記元画像を前記第１の方向に平均縮小した第４の画像に対してフィルタ処理を行って前記第１の画像を生成し、前記第３の画像の代わりに、前記第３の画像を前記第１の方向に平均拡大した第５の画像を前記元画像から減算して前記補正後の画像を生成することを特徴とする請求項１２に記載の画像処理方法。
【請求項１４】
第１の方向に延びるスジ状の濃度ムラを含む元画像に対して、前記第１の方向と直交する第２の方向にフィルタ処理を行って高周波成分が抽出された第１の画像を生成するステップと、
前記第１の画像に対して、前記第２の方向に閾値判定を行って該閾値よりも小さい画素が除去された第２の画像を生成し、
前記第２の画像に対して、前記第１の方向にフィルタ処理を行って低周波成分が抽出された第３の画像を生成するステップと、
前記元画像から前記第３の画像を減算して前記スジ状の濃度ムラが補正された補正後の画像を生成するステップとをコンピュータに実行させることを特徴とする画像処理プログラム。
【請求項１５】
前記元画像の代わりに、前記元画像を前記第１の方向に平均縮小した第４の画像に対してフィルタ処理を行って前記第１の画像を生成するステップと、前記第３の画像の代わりに、前記第３の画像を前記第１の方向に平均拡大した第５の画像を前記元画像から減算して前記補正後の画像を生成するステップとをコンピュータに実行させることを特徴とする請求項１４に記載の画像処理プログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【公開番号】特開２０１１−７６４６７（Ｐ２０１１−７６４６７Ａ）
【公開日】平成２３年４月１４日（２０１１．４．１４）
【国際特許分類】
【出願番号】特願２００９−２２８６９８（Ｐ２００９−２２８６９８）
【出願日】平成２１年９月３０日（２００９．９．３０）
【出願人】（３０６０３７３１１）富士フイルム株式会社 (25,513)
【Ｆターム（参考）】
【公開日】平成２３年４月１４日（２０１１．４．１４）
【国際特許分類】
【出願日】平成２１年９月３０日（２００９．９．３０）
【出願人】（３０６０３７３１１）富士フイルム株式会社 (25,513)
【Ｆターム（参考）】
[ Back to top ]