JP4612774B2 - 画像処理装置、画像処理システム、画像処理方法、記憶媒体、及びプログラム - Google Patents
画像処理装置、画像処理システム、画像処理方法、記憶媒体、及びプログラム Download PDFInfo
- Publication number
- JP4612774B2 JP4612774B2 JP2001034121A JP2001034121A JP4612774B2 JP 4612774 B2 JP4612774 B2 JP 4612774B2 JP 2001034121 A JP2001034121 A JP 2001034121A JP 2001034121 A JP2001034121 A JP 2001034121A JP 4612774 B2 JP4612774 B2 JP 4612774B2
- Authority
- JP
- Japan
- Prior art keywords
- frequency component
- image
- original image
- region
- acquiring
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Image Analysis (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、例えば、Ｘ線等の放射線により撮影された画像に対して、階調変換処理、ダイナミックレンジ変更処理、及び鮮鋭化処理等の画像処理を行う装置或いはシステムに用いられる、画像処理装置、画像処理システム、画像処理方法、前記装置もしくはシステムの機能を実現させるためのプログラム、前記方法の処理ステップを実行させるためのプログラム、及びそれらプログラムをコンピュータが読出可能に格納した記憶媒体に関するものである。
【０００２】
【従来の技術】
近年においては、ディジタル技術の進歩により、放射線撮影して得られた画像（放射線画像）をディジタル化してディジタル画像として取得し、当該ディジタル画像に対して画像処理を施し、当該画像処理後のディジタル画像をＣＲＴ等で表示出力する、或いはフィルム出力することが行われている。
【０００３】
ところで、放射線画像については、そのディジタル画像を表示出力或いはフィルム出力に適した画像とするために、上記の画像処理として、階調変換処理やダイナミックレンジ変更処理を行うことがある。
【０００４】
しかしながら、上述のような画像処理を行ったとしても、例えば、Ｘ線により胸部を撮影して得られた画像（Ｘ線胸部画像）は、Ｘ線が透過しやすい肺野の画像、及びＸ線が非常に透過しにくい縦隔部の画像から構成されるため、画素値（ディジタル画像を構成する画素の値）が存在するレンジが非常に広い。このため、肺野及び縦隔部の両方の部分画像を同時に観察することが可能なＸ線胸部画像を得ることは非常に困難であるとされてきた。
【０００５】
そこで、上記の問題を回避する方法として、例えば、「ＳＰＩＥ Ｖｏｌ．６２６ Ｍｅｄｉｃｉｎｅ ＸＩＶ／ＰＡＣＳ（１９８６）」等に記載されている方法（以下、「第１の方法」と言う）がある。
【０００６】
第１の方法は、処理後の画素値ＳD、オリジナル画素値（入力画素値）Ｓorg、オリジナル画像（入力画像）の低周波画像の画素値ＳUS、及び定数Ａ，Ｂ，Ｃ（例えば、Ａ＝３、Ｂ＝０．７）を以って、
【０００７】
【数１】
【０００８】
なる式（１）で表わされるものである。
【０００９】
また、第１の方法は、高周波成分（上記式（１）の第１項）、低周波成分（上記式（２）の第２項）の重み付けを変えることで、具体的には例えば、Ａ＝３、Ｂ＝０．７とすることで、対象画像の高周波成分を強調し、且つ対象画像全体のダイナミックレンジを圧縮する効果が得られるものである。
【００１０】
尚、第１の方法に対して、ある５人の放射線医から、処理なし画像と比較して、診断に有効であるという評価が得られている。
【００１１】
また、例えば、第２５０９５０３特許公報等には、処理後の画素値ＳD、オリジナル画素値（入力画素値）Ｓorg、及びオリジナル画像（入力画像）のＹ方向プロファイルの平均プロファイルＰｙとＸ方向プロファイルの平均プロファイルＰxを以って、
【００１２】
【数２】
【００１３】
なる式（２）で表わされる方法（以下、「第２の方法」と言う）が記載されている。
【００１４】
上記式（２）において、関数Ｆ（ｘ）は、「ｘ＞Ｄth」では、Ｆ（ｘ）が「０」となり、「０≦ｘ≦Ｄth」では、Ｆ（ｘ）が、切片を「Ｅ」、傾き「Ｅ／Ｄth」として単調減少する、という特性を有し、
【００１５】
【数３】
【００１６】
なる式（３）〜式（５）で示される。
【００１７】
上記式（４）及び（５）において、
ｉ＝１〜ｎ
である。“Ｐyi”及び“Ｐxi”は、Ｙ方向及びX方向のプロファイルを示す。また、Ｇ（Ｐx，Ｐy）は例えば、
【００１８】
【数４】
【００１９】
なる式（６）で示される。
【００２０】
第２の方法は、対象画像の低周波部分の画素値で、“Ｄth”以下の濃度レンジが圧縮されるものである。
【００２１】
また、第２の方法と同様の方法として、「日本放射線技術学会雑誌、第４５巻第８号１９８９年８月、１０３０頁、阿南ほか」や、第２６６３１８９特許公報等に記載された方法（以下、「第３の方法」と言う）がある。
【００２２】
第３の方法は、処理後の画素値ＳD、オリジナル画素値（入力画素値）Ｓorg、オリジナル画像（入力画像）をマスクサイズＭ×Ｍ画素で移動平均をとった時の平均画素値SUS、及び単調減少関数ｆ（Ｘ）を以って、
【００２３】
【数５】
【００２４】
なる式（７）及び（８）式で表わされるものである。
【００２５】
また、第３の方法は、上記式（２）で表される第２の方法に対して、低周波画像の作成方法が異なる。すなわち、上記式（２）では１次元データにより低周波画像が作成されるのに対し、上記式（７）及び（８）では、２次元データにより低周波画像が作成される。
【００２６】
第３の方法も、第２の方法と同様に、対象画像の低周波部分の画素値で、“Ｄth”以下の濃度レンジが圧縮されるものである。
【００２７】
上述したような第１〜第３の方法でのダイナミックレンジ圧縮は、低周波画像を変換する関数ｆ１（）と、定数ＡXとを以って、
【００２８】
【数６】
【００２９】
なる式（９）で表される。
【００３０】
また、特開平１０−２７２２８３号等では、次のような方法（以下、「第４の方法」と言う）が記載されている。
【００３１】
第４の方法は、階調変換関数Ｆ１（）、階調変換率ｃ（ｘ、ｙ）、処理済み画像の画素値ｆｄ（ｘ、ｙ）、第１の画像ｆ０（ｘ、ｙ）、第２の画像の画素値ｆ１（ｘ、ｙ）、第２の画像の平滑化（低周波）画像の画素値ｆus（ｘ、ｙ）、及び画像上の座標（ｘ、ｙ）を持って、
【００３２】
【数７】
【００３３】
なる式で示される関係を有し、
【００３４】
【数８】
【００３５】
なる式（１０）で表される。
【００３６】
第４の方法による画像処理は、階調変換処理との組み合わせ等が容易であり、階調変換処理と同時に行える等の効果を有すると共に、高周波成分の調整も同時に行える、という効果を有する。
【００３７】
上述したような第１〜第４の方法による画像処理の他、画像の微細構造を容易に観察できるようにするための、所謂鮮鋭化処理がある。
例えば、鮮鋭化処理としては、原画像ｆ（ｘ、ｙ）、処理後の画像ｆｐ（ｘ、ｙ）、原画像の低周波画像ＳUS２（ｘ、ｙ）、及び定数Cを以って、
【００３８】
【数９】
【００３９】
なる式（１１）で表される処理がある。
【００４０】
しかしながら、上記式（１１）で表される処理では、原画像ｆ（ｘ、ｙ）から、原画像の低周波画像ＳUS２（ｘ、ｙ）を減じて、高周波画像を作成する処理が行われる。この際、原画像ｆ（ｘ、ｙ）のエッジ部分に対応する高周波画像の部分にオーバーシュートやアンダーシュートのもとになる高周波成分が生じる場合がある。この場合、当該高周波成分が、上記の処理において、原画像、又は階調変換後の原画像、又は低周波画像に加算されることになり、この結果、オーバーシュートやアンダーシュートが生じた処理後の画像が得られてしまうことになる。
【００４１】
尚、ここでオーバーシュート及びアンダーシュートは原画像ｆ（ｘ、ｙ）から原画像の低周波画像ＳUS２（ｘ、ｙ）を減じて得られた高周波成分の中で絶対値が過度に大きな成分を意味し、絶対値が正のものをオーバーシュート、絶対値が負のものをアンダーシュートとする。以下、単にオーバーシュートと称した場合、原則としてオーバーシュート及びアンダーシュートを総称しているものとする。
【００４２】
そこで、上記のオーバーシュートを抑制する方法として、例えば、特開平１１−１２４７４９号等では、対象画像の高周波成分をさらに階調変換をすることで、高周波成分の絶対値の大きな成分でありオーバーシュートに対応する部分を抑制する方法が提案されている。また、中間値フィルタやモーフォロジカル・フィルタを用いて低周波成分を作成することで、上記のオーバーシュートを抑制する方法等も提案されている。
【００４３】
【発明が解決しようとする課題】
しかしながら、上述したようなオーバーシュートを抑制する従来の画像処理方法では、例えば、特開平１１−１２４７４９号等に記載の画像処理方法では、振幅の比較的小さなオーバーシュートについては適切に抑制できるが、振幅の大きなオーバーシュートについては、その抑制の際に高周波成分自体の情報量が失われてしまう。このため、オーバーシュートは抑制されるものの、オーバーシュートが抑制された領域の画像の微細構造が見えにくい、という問題があった。
【００４４】
また、中間値フィルタやモーフォロジカル・フィルタを用いた画像処理方法では、当該フィルタにより低周波成分を作成するため、オーバーシュートが発生しにくくなるものの、多くの処理時間（計算時間）を要し、また、低周波成分を作成する領域によっては、微細構造の復元性が悪くなってしまう、という問題があった。
【００４５】
そこで、本発明は、上記の欠点を除去するために成されたもので、オーバーシュートを適切に抑制しつつ、オーバーシュートが抑制された領域の微細構造の情報をも失わないことを実現する構成により、良好な処理後画像を提供できる、画像処理装置、画像処理システム、画像処理方法、前記装置もしくはシステムの機能を実現させるためのプログラム、前記方法の処理ステップを実行させるためのプログラム、及びそれらのプログラムをコンピュータが読出可能に格納した記憶媒体を提供することを目的とする。
【００４６】
【課題を解決するための手段】
本発明の画像処理装置は、原画像から第１の高周波成分を取得する第１の高周波成分取得手段と、前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を前記原画像から取得する第２の高周波成分取得手段と、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とする手段と、前記高周波成分の大きさを所定の規則で変更した後の高周波成分と階調変換した後の前記原画像とを加算して処理後の画像を取得する処理後画像取得手段と、を備えることを特徴とする。
また、本発明の画像処理装置の他の特徴とするところは、原画像から第１の高周波成分を取得する第１の高周波成分取得手段と、前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を前記原画像から取得する第２の高周波成分取得手段と、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とする手段と、前記高周波成分の大きさを所定の規則で変更した後の高周波成分と前記原画像とを加算して処理後の画像を取得する処理後画像取得手段と、を備えることを特徴とする。
また、本発明の画像処理装置のその他の特徴とするところは、原画像から第１の高周波成分を取得する第１の高周波成分取得手段と、前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を前記原画像から取得する第２の高周波成分取得手段と、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とし、前記高周波成分の大きさを所定の規則で変更した後の高周波成分と前記原画像の低周波成分を加算して処理後の画像を取得する処理後画像取得手段と、を備えることを特徴とする。
また、本発明の画像処理装置のその他の特徴とするところは、原画像から第１のフィルタサイズでフィルタ処理して第１の高周波成分を取得する第１の高周波成分取得手段と、前記原画像から前記第１のフィルタサイズよりも小さい第２のフィルタサイズでフィルタ処理して第２の高周波成分を取得する第２の高周波成分取得手段と、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とし、前記高周波成分と前記原画像の所定の低周波成分とを加算して処理後の画像を取得する処理後画像取得手段と、を備えることを特徴とする。
【００６１】
本発明の画像処理システムは、複数の機器が互いに通信可能に接続されてなる画像処理システムであって、前記複数の機器のうち少なくとも１つの機器は、請求項１〜４の何れか１項に記載の画像処理装置の機能を有することを特徴とする。
【００６２】
本発明の画像処理方法は、原画像から第１の高周波成分を取得する第１の高周波成分取得ステップと、前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を前記原画像から取得する第２の高周波成分取得ステップと、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とするステップと、前記高周波成分の大きさを所定の規則で変更した後の高周波成分と階調変換した後の前記原画像とを加算して処理後の画像を取得する処理後画像取得ステップと、を含むことを特徴とする。
また、本発明の画像処理方法の他の特徴とするところは、原画像から第１の高周波成分を取得する第１の高周波成分取得ステップと、前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を原画像から取得する第２の高周波成分取得ステップと、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とするステップと、前記高周波成分の大きさを所定の規則で変更した後の高周波成分と前記原画像とを加算して処理後の画像を取得する処理後画像取得ステップと、を含むことを特徴とする。
また、本発明の画像処理方法のその他の特徴とするところは、原画像から第１の高周波成分を取得する第１の高周波成分取得ステップと、前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を原画像から取得する第２の高周波成分取得ステップと、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とし、前記高周波成分の大きさを所定の規則で変更した後の高周波成分と前記原画像の低周波成分を加算して処理後の画像を取得する処理後画像取得ステップと、を備えることを特徴とする。
また、本発明の画像処理方法のその他の特徴とするところは、原画像から第１のフィルタサイズでフィルタ処理して第１の高周波成分を取得する第１の高周波成分取得ステップと、前記原画像から前記第１のフィルタサイズよりも小さい第２のフィルタサイズでフィルタ処理して第２の高周波成分を取得する第２の高周波成分取得ステップと、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とし、前記高周波成分と前記原画像の所定の低周波成分とを加算して処理後の画像を取得する処理後画像取得ステップと、を備えることを特徴とする。
【００７８】
本発明の記憶媒体は、請求項１〜４の何れか１項に記載の画像処理装置の機能、又は請求項５記載の画像処理システムの機能をコンピュータに実現させるためのプログラムを記録したことを特徴とする。
また、本発明の記憶媒体の他の特徴とするところは、請求項６〜９の何れか１項に記載の画像処理方法の処理ステップをコンピュータに実行させるためのプログラムを記録したことを特徴とする。
【００８０】
本発明のプログラムは、請求項１〜４の何れか１項に記載の画像処理装置の機能、又は請求項５記載の画像処理システムの機能をコンピュータに実現させることを特徴とする。
また、本発明のプログラムの他の特徴とするところは、請求項６〜９の何れか１項に記載の画像処理方法の処理ステップをコンピュータに実行させることを特徴とする。
【００８２】
【発明の実施の形態】
以下、本発明の実施の形態について図面を用いて説明する。
【００８３】
[第１の実施の形態]
本発明は、例えば、図１に示すようなＸ線撮影装置１００に適用される。
本実施の形態のＸ線撮影装置１００は、特に、以下に説明するような構成及び動作により、対象画像の高周波成分を変更することで、オーバーシュート等を適切に防止するようになされている。
【００８４】
＜Ｘ線撮影装置１００の構成＞
Ｘ線撮影装置１００は、上記図１に示すように、Ｘ線を発生するＸ線発生回路（図面では…部）１０１と、被写体１０３を透過したＸ線光が結像される２次元Ｘ線センサ１０４と、２次元Ｘ線センサ１０４から出力される撮像画像を収集するデータ収集回路（図面では…部）１０５と、データ収集回路１０５にて収集された撮像画像に前処理を行う前処理回路（図面では…部）１０６と、前処理回路１０６での前処理後の撮像画像等の各種情報や各種処理実行のための処理プログラムを記憶するメインメモリ１０９と、Ｘ線撮影実行等の指示や各種設定を本装置に対して行うための操作パネル１１０と、前処理回路１０６での前処理後の撮像画像に対して画像処理を施す画像処理回路１１１と、画像処理回路１１１での処理後の画像等を表示出力或いはフィルム出力等する出力回路１１８と、本装置全体の動作制御を司るＣＰＵ１０８とを備えており、データ収集回路１０５、前処理回路１０６、画像処理回路１１１、出力回路１１８と、ＣＰＵ１０８、メインメモリ１０９、及び操作パネル１１０はそれぞれＣＰＵバス１０７に接続され互いにデータ授受できるように構成されている。
【００８５】
画像処理回路１１１は、前処理回路１０６での前処理後の画像から第１の高周波成分を作成する第１の高周波成分作成回路１１２と、第１の高周波成分作成回路１１２で得られた第１の高周波成分の絶対値からオーバーシュートが生じている領域を検出する領域検出回路１１３と、領域検出回路１１３で得られた検出領域に対応する第２の高周波成分を作成する第２の高周波成分作成回路１１４と、第１の高周波成分作成回路１１２で得られた第１の高周波成分のうち領域検出回路１１３で得られた検出領域に対応する高周波成分を第２の高周波成分作成回路１１４で得られた第２の高周波成分に置換する周波成分置換回路１１５と、周波成分置換回路１１５で得られた置換後の高周波成分を原画像、低周波画像、又は階調変換後の画像に加算して処理後画像を取得する処理後画像回路１１６とを備えている。
【００８６】
＜Ｘ線撮影装置１００の動作＞
メインメモリ１０９には、ＣＰＵ１０８での各種制御処理実行に必要なデータや処理プログラム等が予め記憶されると共に、ＣＰＵ１０８の作業用としてのワークメモリを含むものである。メインメモリ１０９に記憶される処理プログラム、特に、画像処理のための処理プログラムとして、ここでは例えば、図２のフローチャートに従った処理プログラムを用いる。
したがって、ＣＰＵ１０８は、上記図２で示される処理プログラム等をメインメモリ１０９から読み出して実行することで、操作パネル１１０からの操作に従った、以下に説明するような本装置１００全体の動作制御を行う。
【００８７】
ステップＳ２００：
先ず、Ｘ線発生回路１０１は、被写体（被検査体）１０３に対してＸ線ビーム１０２を放射する。このＸ線発生回路１０１から放射されたＸ線ビーム１０２は、被検査体１０３を減衰しながら透過して、２次元Ｘ線センサ１０４に到達し、２次元Ｘ線センサ１０４によりＸ線画像として出力される。
２次元Ｘ線センサ１０４から出力されるＸ線画像としては、例えば、胸部画像等の人体部画像が挙げられる。
【００８８】
次に、データ収集回路１０５は、２次元Ｘ線センサ１０４から出力されたＸ線画像を電気信号に変換し、それを前処理回路１０６に供給する。
前処理回路１０６は、データ収集回路１０５からの信号（Ｘ線画像信号）に対して、オフセット補正処理やゲイン補正処理等の前処理を行う。
この前処理回路１０６で前処理が行われたＸ線画像信号は入力画像の情報として、ＣＰＵ１０８の制御により、ＣＰＵバス１１５を介して、メインメモリ１０９及び画像処理回路１１１にそれぞれ転送される。
【００８９】
ステップＳ２０１：
画像処理回路１１１において、不図示の階調変換回路は、前処理回路１０６での処理後の画像ｆ（ｘ、ｙ）に対して、階調変換関数Ｆ１（）を用いた、
【００９０】
【数１０】
【００９１】
なる式（１２）により階調変換処理を施し、原画像（元画像）ｆ０（ｘ、ｙ）を取得する。
尚、上記式（１２）において、“ｘ”及び“ｙ”はそれぞれ、Ｘ座標及びＹ座標を表す。
【００９２】
ステップＳ２０２：
第１の高周波成分作成回路１１２は、
【００９３】
【数１１】
【００９４】
なる式（１３）及び（１４）で示されるように、ステップＳ２０１で得られた原画像ｆ（ｘ、ｙ）から、その低周波画像ｆus（ｘ、ｙ）を減じることで、第１の高周波成分ｈ１（ｘ、ｙ）を取得する。
【００９５】
上記式（１３）における低周波画像ｆus（ｘ、ｙ）は、上記式（１４）に示されるように、元画像ｆ（ｘ、ｙ）を移動平均して作成されるものであり、上記式（１４）における“ｄ１”は、以下単にマスクサイズと称する。
【００９６】
ステップＳ２０３：
領域検出回路１１３は、第１の高周波成分作成回路１１２にて得られた高周波成分ｈ１（ｘ、ｙ）について、その絶対値が一定閾値Ｔｈ以上の領域を検出し、当該検出領域（オーバーシュートが生じやすい領域、以下、単に「オーバーシュート領域」とも言う）の座標（ｘ、ｙ）を、メインメモリ１０９に保存する。
このように、高周波成分ｈ１（ｘ、ｙ）の絶対値が一定閾値Ｔｈ以上の領域をオーバーシュート領域として検出するのは、オーバーシュートを生じさせる高周波成分は、通常の高周波成分よりも絶対値が大きいことによるものである。
【００９７】
ステップＳ２０４：
第２の高周波成分作成回路１１４は、メインメモリ１０９に保存されたオーバーシュート領域の座標に基づいて、
【００９８】
【数１２】
【００９９】
なる式（１５）及び（１６）により、第２の高周波成分ｈ２（ｘ、ｙ）を取得する。
【０１００】
ここで、上記式（１６）における“ｄ２”は、以下単にマスクサイズと称する。マスクサイズｄ２は上記式（１４）でのマスクサイズｄ１に対して、
ｄ２＜ｄ１
なる関係を有する。これにより、第２の高周波成分ｈ２（ｘ、ｙ）は、第１の高周波成分ｈ１（ｘ、ｙ）より高周波成分に周波数成分がかたよっている。また、当該関係により、第２の高周波成分ｈ２（ｘ、ｙ）は、第１の高周波成分ｈ１（ｘ、ｙ）よりも、オーバーシュートを生じさせる成分を減じるものである。
【０１０１】
尚、マスクサイズｄ２を過度に小さく設定すると対象画像における微細構造の復元性が落ちる場合があるため、オーバーシュートが生じにくい範囲で且つ微細構造の復元性をあまり落とさないようにマスクサイズｄ２を決定するのが望ましい。
この場合、中周波数の復元性は通常の大きさのマスクサイズｄ１を用いる場合より劣化するが、微細構造の復元性が落ちることなく、オーバーシュートを抑制できる。
【０１０２】
また、上記式（１５）の“ｆus２（ｘ、ｙ）”は、上記式（１６）により表される低周波画像（平滑化画像）であるが、この低周波画像ｆus２（ｘ、ｙ）の他の取得方法としては、例えば、次のような方法（１）及び（２）が挙げられる。
【０１０３】
方法（１）：モーフォロジカル・フィルタ（モーフォロジ演算）を用いた方法本方法（１）では、次の式（１７）〜（２０）により表されるモーフォロジカル・フィルタ（モーフォロジ演算）により、低周波画像ｆus２（ｘ、ｙ）を生成する。
【０１０４】
【数１３】
【０１０５】
上記式（１７）〜（２０）において、“D（ｘ、ｙ）”は、次の式（２１）により表される円盤状フィルタであり、“ｒ１”を任意の定数として、対象画像に応じて選択されるものである。
【０１０６】
【数１４】
【０１０７】
方法（２）：中間値フィルタを用いた方法
本方法（２）では、中間値フィルタを用いて低周波画像ｆus２（ｘ、ｙ）を生成する。
具体的には、座標（ｘ、ｙ）を中心として、一定距離ｄ２内にある元画像ｆ（ｘ、ｙ）をソートし、中間にある画素値をfus２（ｘ、ｙ）とする。
【０１０８】
上述のような方法（１）或いは方法（２）を用いて得られる低周波画像ｆus２（ｘ、ｙ）は、エッジ構造を保存するものである。そして、上記式（１５）で示したように、この低周波画像ｆus２（ｘ、ｙ）を用いて得られた第２の高周波成分ｈ２（ｘ、ｙ）は、オーバーシュートのもとになる成分がでにくいものである。また、オーバーシュートが生じやすい領域のみに対して、モーフォロジカル・フィルタ或いは中間値フィルタを用いて取得した低周波画像ｆus２（ｘ、ｙ）を生成することで、対象画像全体に対してモーフォロジカル・フィルタ或いは中間値フィルタを施すように構成した場合よりも、処理時間の短縮を図ることができる。
【０１０９】
ステップＳ２０５：
周波数成分置換回路１１５は、ステップＳ２０２において第１の高周波成分作成回路１１２により得られた第１の高周波成分ｈ１（ｘ、ｙ）を、
【０１１０】
【数１５】
【０１１１】
なる関係式（２２）に基づいて、第２の高周波成分ｈ２（ｘ、ｙ）に置換する。
上記関係式（２２）に表されるように、第１の高周波成分ｈ１（ｘ、ｙ）のうち、絶対値が一定閾値Ｔｈ以上の成分について、第２の高周波成分ｈ２（ｘ、ｙ）に置換される。
【０１１２】
ステップＳ２０６：
処理後画像作成回路１１６は、
【０１１３】
【数１６】
【０１１４】
なる式（２３）に従って、周波成分置換回路１１５で得られた置換後の第１の高周波成分ｈ１（ｘ、ｙ）を、例えば、原画像ｆ０（ｘ、ｙ）に加算することで、処理後画像ｆｄ（ｘ、ｙ）を取得する。
【０１１５】
尚、ステップＳ２０６において、周波成分置換回路１１５で得られた置換後の第１の高周波成分ｈ１（ｘ、ｙ）を加算する対象の画像としては、原画像ｆ０（ｘ、ｙ）に限られることはなく、例えば、前処理後の画像や、低周波画像、或いは低周波画像を階調変換した画像を用いるようにしてもよい。
【０１１６】
ステップＳ２０７：
上述のようにして、画像処理回路１１１により得られた処理後画像ｆｄ（ｘ、ｙ）は、出力回路１１８により、表示出力或いはフィルム出力等される。
【０１１７】
上述のように、本実施の形態では、原画像から取得した第１の高周波成分を解析することで、オーバーシュート領域を検出し、そのオーバーシュート領域の第１の高周波成分を、オーバーシュートが生じにくい高周波成分（第２の高周波成分）に置換するように構成したので、オーバーシュートを抑制した良好な処理後画像を提供できる。
【０１１８】
また、オーバーシュート領域に対してのみ、その第１の高周波成分を、より高周波領域にかたよった高周波成分（マスクサイズｄ２を小とした第２の高周波成分）に置き換えるように構成することで、微細構造の復元性を落とすことなくオーバーシュートを抑制できる。
【０１１９】
さらに、第２の高周波成分を生成する方法として、モーフォロジカル・フィルタ或いは中間値フィルタ等を用いた方法を用い、オーバーシュート領域のみに対して当該処理を行うように構成することで、さらなるオーバーシュートの抑制を期待できると共に、対象画像全体にわたってモーフォロジカル・フィルタ或いは中間値フィルタ等を施すように構成した場合に比較して、処理時間の短縮を行える。
【０１２０】
[第２の実施の形態]
本実施の形態では、上記図１に示したＸ線撮影装置１００の動作を、例えば、図３のフローチャートに従った動作とする。
【０１２１】
尚、上記図３のフローチャートにおいて、上記図２のフローチャートと同様の処理実行するステップには同じ符号を付し、その詳細な説明は省略する。
【０１２２】
ステップＳ２００〜ステップＳ２０３：
第１の実施の形態と同様にして、Ｘ線撮影により得られた画像は、前処理回路１０６で処理され、当該処理後の画像ｆ（ｘ、ｙ）は、画像処理回路１１１へ転送される（ステップＳ２００）。
画像処理回路１１１は、画像ｆ（ｘ、ｙ）に対して階調変換処理を施すことで原画像（元画像）ｆ０（ｘ、ｙ）を取得し（ステップＳ２０１）、第１の高周波成分作成回路１１２により、原画像ｆ（ｘ、ｙ）から第１の高周波成分ｈ１（ｘ、ｙ）を取得する（ステップＳ２０２）、領域検出回路１１３により、第１の高周波成分ｈ１（ｘ、ｙ）からオーバーシュート領域を検出する（ステップＳ２０３）。
【０１２３】
尚、ステップＳ２０３は必須ではなく、後述のようにオーバーシュート領域に対応した領域のみ低周波画像を作成する場合には必要であるが、それ以外の場合には必ずしも必要ではない。
【０１２４】
ステップＳ３０１：
第２の高周波成分作成回路１１４は、第１の高周波成分ｈ１（ｘ、ｙ）をダイナミックレンジ圧縮した結果を、第２の高周波成分ｈ２（ｘ、ｙ）として取得する。
【０１２５】
上記のダイナミックレンジ圧縮の処理方法としては、任意の方法を適用可能であるが、本実施の形態では、その一例として、先に説明した特開平１０−２７２２８３号等に記載の方法を用いる。
したがって、ここでは、第１の高周波成分ｈ１（ｘ、ｙ）の絶対値が、一定閾値Ｔｈ以上すなわちオーバーシュート領域のみのダイナミックレンジが圧縮されるように、ダイナミックレンジの圧縮処理を行う。このときの圧縮率ついては任意であるが、本実施の形態では、その一例として、最大の圧縮率とする。これにより、一定閾値Ｔｈ以上の第１の高周波成分ｈ１（ｘ、ｙ）が可能な限り圧縮されることになる。例えば、第１の高周波成分ｈ１（ｘ、ｙ）中の所定の低周波成分が概略一定閾値Ｔｈ以下にクリッピングされるように第１の高周波成分ｈ１（ｘ、ｙ）をダイナミックレンジ圧縮する。
【０１２６】
ここで、ダイナミックレンジ圧縮に用いる低周波画像（平滑化画像）ｆus３（ｘ、ｙ）について、その取得方法としては、例えば、次のような方法（１）〜方法（３）が挙げられる。
【０１２７】
方法（１）：オーバーシュート領域の座標に基づいた方法
本方法（１）では、メインメモリ１０９に保存されたオーバーシュート領域の座標に基づいて、
【０１２８】
【数１７】
【０１２９】
なる式（２４）により、低周波画像ｆus３（ｘ、ｙ）を取得する。
上記式（２４）における“ｄ２”は、マスクサイズであり、上記式（１４）でのマスクサイズｄ１に対して、
ｄ２＜ｄ１
なる関係を有する。
【０１３０】
方法（２）：モーフォロジカル・フィルタを用いた方法
本方法（２）では、次の式（２５）〜（２８）により表されるモーフォロジカル・フィルタにより、低周波画像ｆus３（ｘ、ｙ）を生成する。
【０１３１】
【数１８】
【０１３２】
上記式（２５）〜（２８）において、“D（ｘ、ｙ）”は、次の式（２９）により表される円盤状フィルタであり、“ｒ１”を任意の定数として、対象画像に応じて選択されるものである。
【０１３３】
【数１９】
【０１３４】
方法（３）：中間値フィルタを用いた方法
本方法（３）では、中間値フィルタを用いて低周波画像ｆus３（ｘ、ｙ）を生成する。
具体的には、座標（ｘ、ｙ）を中心として、一定距離ｄ２内にある元画像ｆ（ｘ、ｙ）をソートし、中間にある画素値をfus３（ｘ、ｙ）とする。
【０１３５】
上述のような方法（２）或いは方法（３）を用いて得られる低周波画像ｆus３（ｘ、ｙ）は、エッジ構造を保存するものであり、この低周波画像ｆus３（ｘ、ｙ）を用いて得られた第２の高周波成分ｈ２（ｘ、ｙ）は、オーバーシュートのもとになる成分がでにくいものである。また、オーバーシュートが生じやすい領域のみに対して、モーフォロジカル・フィルタ或いは中間値フィルタを用いて取得した低周波画像ｆus３（ｘ、ｙ）を生成することで、対象画像全体に対してモーフォロジカル・フィルタ或いは中間値フィルタを施すように構成した場合よりも、処理時間の短縮を図ることができる。
【０１３６】
ステップＳ３０２：
処理後画像作成回路１１６は、
【０１３７】
【数２０】
【０１３８】
なる式（３０）〜式（３２）により、処理後画像ｆｄ（ｘ、ｙ）を取得する。
【０１３９】
上記式（３０）〜式（３２）において、“ｈ１（ｘ、ｙ）”は、第２の高周波成分作成回路１１４によりダイナミックレンジ圧縮された後の第１の高周波成分（第２の高周波成分）である。また、上記式（３０）及び（３１）は、鮮鋭化処理を示し、上記式（３２）は、ダイナミックレンジの圧縮処理を示す。また、上記式（３２）において、“Ｆ２（）”は、低周波画像を階調変換する関数を示す。
【０１４０】
ステップＳ２０７：
上述のようにして、画像処理回路１１１により得られた処理後画像ｆｄ（ｘ、ｙ）は、出力回路１１８により、表示出力或いはフィルム出力等される。
【０１４１】
上述のように、本実施の形態では、第２の高周波成分として、第１の高周波成分をダイナミックレンジ圧縮した成分を用いるように構成した。つまり、オーバーシュートが生じやすい領域に対応する第１の高周波成分中の所定の低周波成分を圧縮するように構成したため、微細構造の復元性を損なわずにオーバーシュートを抑制することができる。また、第１の高周波成分の圧縮方法として、所定の低周波成分を圧縮すると共に所定の高周波成分を調節可能なダイナミックレンジ圧縮方法を用いる構成としたので、第１の高周波成分を自在に変更でき、もって微細構造の復元性を損なわずにオーバーシュートを抑制可能な第２の高周波成分を得ることができる。さらに、オーバーシュートが生じやすい領域に対してのみ、モーフォロジカルフィルタまたは中間値フィルタ等を用いてダイナミックレンジ圧縮に用いる低周波画像を作成するように構成すると、効果的にオーバーシュートを抑制することができると共に、第１の高周波成分画像全体に対しモーフォロジカルフィルタまたは中間値フィルタ等を用いる場合に比較して処理時間の短縮を図ることができる。
【０１４２】
また、ダイナミックレンジ圧縮に用いる低周波画像のマスクサイズｄ２を、ｄ１より小さく設定するように構成することで、オーバーシュートを抑制することができる。
【０１４３】
また、第１の高周波成分を、オーバーシュートが生じやすい領域において、より高周波領域にかたよった高周波成分（マスクサイズｄ２を小とした第２の高周波成分）に置き換えるように構成することで、微細構造の復元性を落とすことなくオーバーシュートを抑制できる。
【０１４４】
尚、本発明の目的は、第１及び第２の実施の形態のホスト及び端末の機能を実現するソフトウェアのプログラムコードを記憶した記憶媒体を、システム或いは装置に供給し、そのシステム或いは装置のコンピュータ（又はＣＰＵやＭＰＵ）が記憶媒体に格納されたプログラムコードを読みだして実行することによっても、達成されることは言うまでもない。
この場合、記憶媒体から読み出されたプログラムコード自体が第１及び第２の実施の形態の機能を実現することとなり、そのプログラムコード、及び当該プログラムコードを記憶した記憶媒体は本発明を構成することとなる。
プログラムコードを供給するための記憶媒体としては、ＲＯＭ、フロッピーディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＣＤ−Ｒ、磁気テープ、不揮発性のメモリカード等を用いることができる。
また、コンピュータが読みだしたプログラムコードを実行することにより、第１及び第２の実施の形態の機能が実現されるだけでなく、そのプログラムコードの指示に基づき、コンピュータ上で稼動しているＯＳ等が実際の処理の一部又は全部を行い、その処理によって第１及び第２の実施の形態の機能が実現される場合も含まれることは言うまでもない。
さらに、記憶媒体から読み出されたプログラムコードが、コンピュータに挿入された拡張機能ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに書き込まれた後、そのプログラムコードの指示に基づき、その機能拡張ボードや機能拡張ユニットに備わるＣＰＵなどが実際の処理の一部又は全部を行い、その処理によって第１及び第２の実施の形態の機能が実現される場合も含まれることは言うまでもない。
【０１４５】
【発明の効果】
以上説明したように本発明では、原画像（放射線撮影により得られた画像等）において、第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を前記原画像から取得し、その第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を第１の高周波成分とするように構成した。具体的には例えば、第２の高周波成分を、オーバーシュート領域に基づいて、オーバーシュートが生じにくい高周波成分として取得し、この第２の高周波成分で、オーバーシュート領域に対応する第１の高周波成分を置換する。これにより、オーバーシュートを抑制した良好な対象画像の処理後画像（表示出力やフィルム出力用の画像等）を提供できる。
【０１４６】
また、オーバーシュート領域のみに対して、第２の高周波成分を、第１の高周波成分より高周波領域にかたよった高周波成分（第１の高周波成分から所定の低周波成分を低減した高周波成分等）とするように構成した場合、対象画像における微細構造の復元性を落とすことなく、オーバーシュートを抑制することができる。
【０１４７】
また、第２の高周波成分を、第１の高周波成分をダイナミックレンジ圧縮した成分とするように構成した場合、対象画像における微細構造の復元性を維持した状態で、オーバーシュートを抑制することができる。
また、第１の高周波成分の圧縮として、第１の高周波成分中の所定の低周波成分を圧縮し、所定の高周波成分を調節可能なダイナミックレンジ圧縮を用いるように構成した場合、オーバーシュートを抑制することができると共に、良好な処理後画像を取得することができる。
【０１４９】
また、第２の高周波成分の取得に用いる平滑化画像を、中間値フィルタ（中間値演算）やモーフォロジカル・フィルタ（モーフォロジ演算）等のオーバーシュートの生じにくいフィルタを用いて取得するように構成した場合、オーバーシュートの抑制された良好な処理後画像を取得することができる。
【図面の簡単な説明】
【図１】第１の実施の形態において、本発明を適用したＸ線撮影装置の構成を示すブロック図である。
【図２】上記Ｘ線撮影装置の動作を説明するためのフローチャートである。
【図３】第２の実施の形態における、上記Ｘ線撮影装置の動作を説明するためのフローチャートである。
【符号の説明】
１００ Ｘ線撮影装置
１０１ Ｘ線発生回路
１０２ Ｘ線ビーム
１０３ 被写体
１０４ ２次元Ｘ線センサ
１０５ データ収集回路
１０６ 前処理回路
１０７ ＣＰＵバス
１０８ ＣＰＵ
１０９ メインメモリ
１１０ 操作パネル
１１１ 画像処理回路
１１２ 第１の高周波成分作成回路
１１３ 領域検出回路
１１４ 第２の高周波成分作成回路
１１５ 周波数成分置換回路
１１６ 処理後画像作成回路
１１８ 出力回路
【発明の属する技術分野】
本発明は、例えば、Ｘ線等の放射線により撮影された画像に対して、階調変換処理、ダイナミックレンジ変更処理、及び鮮鋭化処理等の画像処理を行う装置或いはシステムに用いられる、画像処理装置、画像処理システム、画像処理方法、前記装置もしくはシステムの機能を実現させるためのプログラム、前記方法の処理ステップを実行させるためのプログラム、及びそれらプログラムをコンピュータが読出可能に格納した記憶媒体に関するものである。
【０００２】
【従来の技術】
近年においては、ディジタル技術の進歩により、放射線撮影して得られた画像（放射線画像）をディジタル化してディジタル画像として取得し、当該ディジタル画像に対して画像処理を施し、当該画像処理後のディジタル画像をＣＲＴ等で表示出力する、或いはフィルム出力することが行われている。
【０００３】
ところで、放射線画像については、そのディジタル画像を表示出力或いはフィルム出力に適した画像とするために、上記の画像処理として、階調変換処理やダイナミックレンジ変更処理を行うことがある。
【０００４】
しかしながら、上述のような画像処理を行ったとしても、例えば、Ｘ線により胸部を撮影して得られた画像（Ｘ線胸部画像）は、Ｘ線が透過しやすい肺野の画像、及びＸ線が非常に透過しにくい縦隔部の画像から構成されるため、画素値（ディジタル画像を構成する画素の値）が存在するレンジが非常に広い。このため、肺野及び縦隔部の両方の部分画像を同時に観察することが可能なＸ線胸部画像を得ることは非常に困難であるとされてきた。
【０００５】
そこで、上記の問題を回避する方法として、例えば、「ＳＰＩＥ Ｖｏｌ．６２６ Ｍｅｄｉｃｉｎｅ ＸＩＶ／ＰＡＣＳ（１９８６）」等に記載されている方法（以下、「第１の方法」と言う）がある。
【０００６】
第１の方法は、処理後の画素値ＳD、オリジナル画素値（入力画素値）Ｓorg、オリジナル画像（入力画像）の低周波画像の画素値ＳUS、及び定数Ａ，Ｂ，Ｃ（例えば、Ａ＝３、Ｂ＝０．７）を以って、
【０００７】
【数１】
なる式（１）で表わされるものである。
【０００９】
また、第１の方法は、高周波成分（上記式（１）の第１項）、低周波成分（上記式（２）の第２項）の重み付けを変えることで、具体的には例えば、Ａ＝３、Ｂ＝０．７とすることで、対象画像の高周波成分を強調し、且つ対象画像全体のダイナミックレンジを圧縮する効果が得られるものである。
【００１０】
尚、第１の方法に対して、ある５人の放射線医から、処理なし画像と比較して、診断に有効であるという評価が得られている。
【００１１】
また、例えば、第２５０９５０３特許公報等には、処理後の画素値ＳD、オリジナル画素値（入力画素値）Ｓorg、及びオリジナル画像（入力画像）のＹ方向プロファイルの平均プロファイルＰｙとＸ方向プロファイルの平均プロファイルＰxを以って、
【００１２】
【数２】
なる式（２）で表わされる方法（以下、「第２の方法」と言う）が記載されている。
【００１４】
上記式（２）において、関数Ｆ（ｘ）は、「ｘ＞Ｄth」では、Ｆ（ｘ）が「０」となり、「０≦ｘ≦Ｄth」では、Ｆ（ｘ）が、切片を「Ｅ」、傾き「Ｅ／Ｄth」として単調減少する、という特性を有し、
【００１５】
【数３】
なる式（３）〜式（５）で示される。
【００１７】
上記式（４）及び（５）において、
ｉ＝１〜ｎ
である。“Ｐyi”及び“Ｐxi”は、Ｙ方向及びX方向のプロファイルを示す。また、Ｇ（Ｐx，Ｐy）は例えば、
【００１８】
【数４】
なる式（６）で示される。
【００２０】
第２の方法は、対象画像の低周波部分の画素値で、“Ｄth”以下の濃度レンジが圧縮されるものである。
【００２１】
また、第２の方法と同様の方法として、「日本放射線技術学会雑誌、第４５巻第８号１９８９年８月、１０３０頁、阿南ほか」や、第２６６３１８９特許公報等に記載された方法（以下、「第３の方法」と言う）がある。
【００２２】
第３の方法は、処理後の画素値ＳD、オリジナル画素値（入力画素値）Ｓorg、オリジナル画像（入力画像）をマスクサイズＭ×Ｍ画素で移動平均をとった時の平均画素値SUS、及び単調減少関数ｆ（Ｘ）を以って、
【００２３】
【数５】
なる式（７）及び（８）式で表わされるものである。
【００２５】
また、第３の方法は、上記式（２）で表される第２の方法に対して、低周波画像の作成方法が異なる。すなわち、上記式（２）では１次元データにより低周波画像が作成されるのに対し、上記式（７）及び（８）では、２次元データにより低周波画像が作成される。
【００２６】
第３の方法も、第２の方法と同様に、対象画像の低周波部分の画素値で、“Ｄth”以下の濃度レンジが圧縮されるものである。
【００２７】
上述したような第１〜第３の方法でのダイナミックレンジ圧縮は、低周波画像を変換する関数ｆ１（）と、定数ＡXとを以って、
【００２８】
【数６】
なる式（９）で表される。
【００３０】
また、特開平１０−２７２２８３号等では、次のような方法（以下、「第４の方法」と言う）が記載されている。
【００３１】
第４の方法は、階調変換関数Ｆ１（）、階調変換率ｃ（ｘ、ｙ）、処理済み画像の画素値ｆｄ（ｘ、ｙ）、第１の画像ｆ０（ｘ、ｙ）、第２の画像の画素値ｆ１（ｘ、ｙ）、第２の画像の平滑化（低周波）画像の画素値ｆus（ｘ、ｙ）、及び画像上の座標（ｘ、ｙ）を持って、
【００３２】
【数７】
なる式で示される関係を有し、
【００３４】
【数８】
なる式（１０）で表される。
【００３６】
第４の方法による画像処理は、階調変換処理との組み合わせ等が容易であり、階調変換処理と同時に行える等の効果を有すると共に、高周波成分の調整も同時に行える、という効果を有する。
【００３７】
上述したような第１〜第４の方法による画像処理の他、画像の微細構造を容易に観察できるようにするための、所謂鮮鋭化処理がある。
例えば、鮮鋭化処理としては、原画像ｆ（ｘ、ｙ）、処理後の画像ｆｐ（ｘ、ｙ）、原画像の低周波画像ＳUS２（ｘ、ｙ）、及び定数Cを以って、
【００３８】
【数９】
なる式（１１）で表される処理がある。
【００４０】
しかしながら、上記式（１１）で表される処理では、原画像ｆ（ｘ、ｙ）から、原画像の低周波画像ＳUS２（ｘ、ｙ）を減じて、高周波画像を作成する処理が行われる。この際、原画像ｆ（ｘ、ｙ）のエッジ部分に対応する高周波画像の部分にオーバーシュートやアンダーシュートのもとになる高周波成分が生じる場合がある。この場合、当該高周波成分が、上記の処理において、原画像、又は階調変換後の原画像、又は低周波画像に加算されることになり、この結果、オーバーシュートやアンダーシュートが生じた処理後の画像が得られてしまうことになる。
【００４１】
尚、ここでオーバーシュート及びアンダーシュートは原画像ｆ（ｘ、ｙ）から原画像の低周波画像ＳUS２（ｘ、ｙ）を減じて得られた高周波成分の中で絶対値が過度に大きな成分を意味し、絶対値が正のものをオーバーシュート、絶対値が負のものをアンダーシュートとする。以下、単にオーバーシュートと称した場合、原則としてオーバーシュート及びアンダーシュートを総称しているものとする。
【００４２】
そこで、上記のオーバーシュートを抑制する方法として、例えば、特開平１１−１２４７４９号等では、対象画像の高周波成分をさらに階調変換をすることで、高周波成分の絶対値の大きな成分でありオーバーシュートに対応する部分を抑制する方法が提案されている。また、中間値フィルタやモーフォロジカル・フィルタを用いて低周波成分を作成することで、上記のオーバーシュートを抑制する方法等も提案されている。
【００４３】
【発明が解決しようとする課題】
しかしながら、上述したようなオーバーシュートを抑制する従来の画像処理方法では、例えば、特開平１１−１２４７４９号等に記載の画像処理方法では、振幅の比較的小さなオーバーシュートについては適切に抑制できるが、振幅の大きなオーバーシュートについては、その抑制の際に高周波成分自体の情報量が失われてしまう。このため、オーバーシュートは抑制されるものの、オーバーシュートが抑制された領域の画像の微細構造が見えにくい、という問題があった。
【００４４】
また、中間値フィルタやモーフォロジカル・フィルタを用いた画像処理方法では、当該フィルタにより低周波成分を作成するため、オーバーシュートが発生しにくくなるものの、多くの処理時間（計算時間）を要し、また、低周波成分を作成する領域によっては、微細構造の復元性が悪くなってしまう、という問題があった。
【００４５】
そこで、本発明は、上記の欠点を除去するために成されたもので、オーバーシュートを適切に抑制しつつ、オーバーシュートが抑制された領域の微細構造の情報をも失わないことを実現する構成により、良好な処理後画像を提供できる、画像処理装置、画像処理システム、画像処理方法、前記装置もしくはシステムの機能を実現させるためのプログラム、前記方法の処理ステップを実行させるためのプログラム、及びそれらのプログラムをコンピュータが読出可能に格納した記憶媒体を提供することを目的とする。
【００４６】
【課題を解決するための手段】
本発明の画像処理装置は、原画像から第１の高周波成分を取得する第１の高周波成分取得手段と、前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を前記原画像から取得する第２の高周波成分取得手段と、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とする手段と、前記高周波成分の大きさを所定の規則で変更した後の高周波成分と階調変換した後の前記原画像とを加算して処理後の画像を取得する処理後画像取得手段と、を備えることを特徴とする。
また、本発明の画像処理装置の他の特徴とするところは、原画像から第１の高周波成分を取得する第１の高周波成分取得手段と、前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を前記原画像から取得する第２の高周波成分取得手段と、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とする手段と、前記高周波成分の大きさを所定の規則で変更した後の高周波成分と前記原画像とを加算して処理後の画像を取得する処理後画像取得手段と、を備えることを特徴とする。
また、本発明の画像処理装置のその他の特徴とするところは、原画像から第１の高周波成分を取得する第１の高周波成分取得手段と、前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を前記原画像から取得する第２の高周波成分取得手段と、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とし、前記高周波成分の大きさを所定の規則で変更した後の高周波成分と前記原画像の低周波成分を加算して処理後の画像を取得する処理後画像取得手段と、を備えることを特徴とする。
また、本発明の画像処理装置のその他の特徴とするところは、原画像から第１のフィルタサイズでフィルタ処理して第１の高周波成分を取得する第１の高周波成分取得手段と、前記原画像から前記第１のフィルタサイズよりも小さい第２のフィルタサイズでフィルタ処理して第２の高周波成分を取得する第２の高周波成分取得手段と、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とし、前記高周波成分と前記原画像の所定の低周波成分とを加算して処理後の画像を取得する処理後画像取得手段と、を備えることを特徴とする。
【００６１】
本発明の画像処理システムは、複数の機器が互いに通信可能に接続されてなる画像処理システムであって、前記複数の機器のうち少なくとも１つの機器は、請求項１〜４の何れか１項に記載の画像処理装置の機能を有することを特徴とする。
【００６２】
本発明の画像処理方法は、原画像から第１の高周波成分を取得する第１の高周波成分取得ステップと、前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を前記原画像から取得する第２の高周波成分取得ステップと、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とするステップと、前記高周波成分の大きさを所定の規則で変更した後の高周波成分と階調変換した後の前記原画像とを加算して処理後の画像を取得する処理後画像取得ステップと、を含むことを特徴とする。
また、本発明の画像処理方法の他の特徴とするところは、原画像から第１の高周波成分を取得する第１の高周波成分取得ステップと、前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を原画像から取得する第２の高周波成分取得ステップと、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とするステップと、前記高周波成分の大きさを所定の規則で変更した後の高周波成分と前記原画像とを加算して処理後の画像を取得する処理後画像取得ステップと、を含むことを特徴とする。
また、本発明の画像処理方法のその他の特徴とするところは、原画像から第１の高周波成分を取得する第１の高周波成分取得ステップと、前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を原画像から取得する第２の高周波成分取得ステップと、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とし、前記高周波成分の大きさを所定の規則で変更した後の高周波成分と前記原画像の低周波成分を加算して処理後の画像を取得する処理後画像取得ステップと、を備えることを特徴とする。
また、本発明の画像処理方法のその他の特徴とするところは、原画像から第１のフィルタサイズでフィルタ処理して第１の高周波成分を取得する第１の高周波成分取得ステップと、前記原画像から前記第１のフィルタサイズよりも小さい第２のフィルタサイズでフィルタ処理して第２の高周波成分を取得する第２の高周波成分取得ステップと、前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とし、前記高周波成分と前記原画像の所定の低周波成分とを加算して処理後の画像を取得する処理後画像取得ステップと、を備えることを特徴とする。
【００７８】
本発明の記憶媒体は、請求項１〜４の何れか１項に記載の画像処理装置の機能、又は請求項５記載の画像処理システムの機能をコンピュータに実現させるためのプログラムを記録したことを特徴とする。
また、本発明の記憶媒体の他の特徴とするところは、請求項６〜９の何れか１項に記載の画像処理方法の処理ステップをコンピュータに実行させるためのプログラムを記録したことを特徴とする。
【００８０】
本発明のプログラムは、請求項１〜４の何れか１項に記載の画像処理装置の機能、又は請求項５記載の画像処理システムの機能をコンピュータに実現させることを特徴とする。
また、本発明のプログラムの他の特徴とするところは、請求項６〜９の何れか１項に記載の画像処理方法の処理ステップをコンピュータに実行させることを特徴とする。
【００８２】
【発明の実施の形態】
以下、本発明の実施の形態について図面を用いて説明する。
【００８３】
[第１の実施の形態]
本発明は、例えば、図１に示すようなＸ線撮影装置１００に適用される。
本実施の形態のＸ線撮影装置１００は、特に、以下に説明するような構成及び動作により、対象画像の高周波成分を変更することで、オーバーシュート等を適切に防止するようになされている。
【００８４】
＜Ｘ線撮影装置１００の構成＞
Ｘ線撮影装置１００は、上記図１に示すように、Ｘ線を発生するＸ線発生回路（図面では…部）１０１と、被写体１０３を透過したＸ線光が結像される２次元Ｘ線センサ１０４と、２次元Ｘ線センサ１０４から出力される撮像画像を収集するデータ収集回路（図面では…部）１０５と、データ収集回路１０５にて収集された撮像画像に前処理を行う前処理回路（図面では…部）１０６と、前処理回路１０６での前処理後の撮像画像等の各種情報や各種処理実行のための処理プログラムを記憶するメインメモリ１０９と、Ｘ線撮影実行等の指示や各種設定を本装置に対して行うための操作パネル１１０と、前処理回路１０６での前処理後の撮像画像に対して画像処理を施す画像処理回路１１１と、画像処理回路１１１での処理後の画像等を表示出力或いはフィルム出力等する出力回路１１８と、本装置全体の動作制御を司るＣＰＵ１０８とを備えており、データ収集回路１０５、前処理回路１０６、画像処理回路１１１、出力回路１１８と、ＣＰＵ１０８、メインメモリ１０９、及び操作パネル１１０はそれぞれＣＰＵバス１０７に接続され互いにデータ授受できるように構成されている。
【００８５】
画像処理回路１１１は、前処理回路１０６での前処理後の画像から第１の高周波成分を作成する第１の高周波成分作成回路１１２と、第１の高周波成分作成回路１１２で得られた第１の高周波成分の絶対値からオーバーシュートが生じている領域を検出する領域検出回路１１３と、領域検出回路１１３で得られた検出領域に対応する第２の高周波成分を作成する第２の高周波成分作成回路１１４と、第１の高周波成分作成回路１１２で得られた第１の高周波成分のうち領域検出回路１１３で得られた検出領域に対応する高周波成分を第２の高周波成分作成回路１１４で得られた第２の高周波成分に置換する周波成分置換回路１１５と、周波成分置換回路１１５で得られた置換後の高周波成分を原画像、低周波画像、又は階調変換後の画像に加算して処理後画像を取得する処理後画像回路１１６とを備えている。
【００８６】
＜Ｘ線撮影装置１００の動作＞
メインメモリ１０９には、ＣＰＵ１０８での各種制御処理実行に必要なデータや処理プログラム等が予め記憶されると共に、ＣＰＵ１０８の作業用としてのワークメモリを含むものである。メインメモリ１０９に記憶される処理プログラム、特に、画像処理のための処理プログラムとして、ここでは例えば、図２のフローチャートに従った処理プログラムを用いる。
したがって、ＣＰＵ１０８は、上記図２で示される処理プログラム等をメインメモリ１０９から読み出して実行することで、操作パネル１１０からの操作に従った、以下に説明するような本装置１００全体の動作制御を行う。
【００８７】
ステップＳ２００：
先ず、Ｘ線発生回路１０１は、被写体（被検査体）１０３に対してＸ線ビーム１０２を放射する。このＸ線発生回路１０１から放射されたＸ線ビーム１０２は、被検査体１０３を減衰しながら透過して、２次元Ｘ線センサ１０４に到達し、２次元Ｘ線センサ１０４によりＸ線画像として出力される。
２次元Ｘ線センサ１０４から出力されるＸ線画像としては、例えば、胸部画像等の人体部画像が挙げられる。
【００８８】
次に、データ収集回路１０５は、２次元Ｘ線センサ１０４から出力されたＸ線画像を電気信号に変換し、それを前処理回路１０６に供給する。
前処理回路１０６は、データ収集回路１０５からの信号（Ｘ線画像信号）に対して、オフセット補正処理やゲイン補正処理等の前処理を行う。
この前処理回路１０６で前処理が行われたＸ線画像信号は入力画像の情報として、ＣＰＵ１０８の制御により、ＣＰＵバス１１５を介して、メインメモリ１０９及び画像処理回路１１１にそれぞれ転送される。
【００８９】
ステップＳ２０１：
画像処理回路１１１において、不図示の階調変換回路は、前処理回路１０６での処理後の画像ｆ（ｘ、ｙ）に対して、階調変換関数Ｆ１（）を用いた、
【００９０】
【数１０】
なる式（１２）により階調変換処理を施し、原画像（元画像）ｆ０（ｘ、ｙ）を取得する。
尚、上記式（１２）において、“ｘ”及び“ｙ”はそれぞれ、Ｘ座標及びＹ座標を表す。
【００９２】
ステップＳ２０２：
第１の高周波成分作成回路１１２は、
【００９３】
【数１１】
なる式（１３）及び（１４）で示されるように、ステップＳ２０１で得られた原画像ｆ（ｘ、ｙ）から、その低周波画像ｆus（ｘ、ｙ）を減じることで、第１の高周波成分ｈ１（ｘ、ｙ）を取得する。
【００９５】
上記式（１３）における低周波画像ｆus（ｘ、ｙ）は、上記式（１４）に示されるように、元画像ｆ（ｘ、ｙ）を移動平均して作成されるものであり、上記式（１４）における“ｄ１”は、以下単にマスクサイズと称する。
【００９６】
ステップＳ２０３：
領域検出回路１１３は、第１の高周波成分作成回路１１２にて得られた高周波成分ｈ１（ｘ、ｙ）について、その絶対値が一定閾値Ｔｈ以上の領域を検出し、当該検出領域（オーバーシュートが生じやすい領域、以下、単に「オーバーシュート領域」とも言う）の座標（ｘ、ｙ）を、メインメモリ１０９に保存する。
このように、高周波成分ｈ１（ｘ、ｙ）の絶対値が一定閾値Ｔｈ以上の領域をオーバーシュート領域として検出するのは、オーバーシュートを生じさせる高周波成分は、通常の高周波成分よりも絶対値が大きいことによるものである。
【００９７】
ステップＳ２０４：
第２の高周波成分作成回路１１４は、メインメモリ１０９に保存されたオーバーシュート領域の座標に基づいて、
【００９８】
【数１２】
なる式（１５）及び（１６）により、第２の高周波成分ｈ２（ｘ、ｙ）を取得する。
【０１００】
ここで、上記式（１６）における“ｄ２”は、以下単にマスクサイズと称する。マスクサイズｄ２は上記式（１４）でのマスクサイズｄ１に対して、
ｄ２＜ｄ１
なる関係を有する。これにより、第２の高周波成分ｈ２（ｘ、ｙ）は、第１の高周波成分ｈ１（ｘ、ｙ）より高周波成分に周波数成分がかたよっている。また、当該関係により、第２の高周波成分ｈ２（ｘ、ｙ）は、第１の高周波成分ｈ１（ｘ、ｙ）よりも、オーバーシュートを生じさせる成分を減じるものである。
【０１０１】
尚、マスクサイズｄ２を過度に小さく設定すると対象画像における微細構造の復元性が落ちる場合があるため、オーバーシュートが生じにくい範囲で且つ微細構造の復元性をあまり落とさないようにマスクサイズｄ２を決定するのが望ましい。
この場合、中周波数の復元性は通常の大きさのマスクサイズｄ１を用いる場合より劣化するが、微細構造の復元性が落ちることなく、オーバーシュートを抑制できる。
【０１０２】
また、上記式（１５）の“ｆus２（ｘ、ｙ）”は、上記式（１６）により表される低周波画像（平滑化画像）であるが、この低周波画像ｆus２（ｘ、ｙ）の他の取得方法としては、例えば、次のような方法（１）及び（２）が挙げられる。
【０１０３】
方法（１）：モーフォロジカル・フィルタ（モーフォロジ演算）を用いた方法本方法（１）では、次の式（１７）〜（２０）により表されるモーフォロジカル・フィルタ（モーフォロジ演算）により、低周波画像ｆus２（ｘ、ｙ）を生成する。
【０１０４】
【数１３】
上記式（１７）〜（２０）において、“D（ｘ、ｙ）”は、次の式（２１）により表される円盤状フィルタであり、“ｒ１”を任意の定数として、対象画像に応じて選択されるものである。
【０１０６】
【数１４】
方法（２）：中間値フィルタを用いた方法
本方法（２）では、中間値フィルタを用いて低周波画像ｆus２（ｘ、ｙ）を生成する。
具体的には、座標（ｘ、ｙ）を中心として、一定距離ｄ２内にある元画像ｆ（ｘ、ｙ）をソートし、中間にある画素値をfus２（ｘ、ｙ）とする。
【０１０８】
上述のような方法（１）或いは方法（２）を用いて得られる低周波画像ｆus２（ｘ、ｙ）は、エッジ構造を保存するものである。そして、上記式（１５）で示したように、この低周波画像ｆus２（ｘ、ｙ）を用いて得られた第２の高周波成分ｈ２（ｘ、ｙ）は、オーバーシュートのもとになる成分がでにくいものである。また、オーバーシュートが生じやすい領域のみに対して、モーフォロジカル・フィルタ或いは中間値フィルタを用いて取得した低周波画像ｆus２（ｘ、ｙ）を生成することで、対象画像全体に対してモーフォロジカル・フィルタ或いは中間値フィルタを施すように構成した場合よりも、処理時間の短縮を図ることができる。
【０１０９】
ステップＳ２０５：
周波数成分置換回路１１５は、ステップＳ２０２において第１の高周波成分作成回路１１２により得られた第１の高周波成分ｈ１（ｘ、ｙ）を、
【０１１０】
【数１５】
なる関係式（２２）に基づいて、第２の高周波成分ｈ２（ｘ、ｙ）に置換する。
上記関係式（２２）に表されるように、第１の高周波成分ｈ１（ｘ、ｙ）のうち、絶対値が一定閾値Ｔｈ以上の成分について、第２の高周波成分ｈ２（ｘ、ｙ）に置換される。
【０１１２】
ステップＳ２０６：
処理後画像作成回路１１６は、
【０１１３】
【数１６】
なる式（２３）に従って、周波成分置換回路１１５で得られた置換後の第１の高周波成分ｈ１（ｘ、ｙ）を、例えば、原画像ｆ０（ｘ、ｙ）に加算することで、処理後画像ｆｄ（ｘ、ｙ）を取得する。
【０１１５】
尚、ステップＳ２０６において、周波成分置換回路１１５で得られた置換後の第１の高周波成分ｈ１（ｘ、ｙ）を加算する対象の画像としては、原画像ｆ０（ｘ、ｙ）に限られることはなく、例えば、前処理後の画像や、低周波画像、或いは低周波画像を階調変換した画像を用いるようにしてもよい。
【０１１６】
ステップＳ２０７：
上述のようにして、画像処理回路１１１により得られた処理後画像ｆｄ（ｘ、ｙ）は、出力回路１１８により、表示出力或いはフィルム出力等される。
【０１１７】
上述のように、本実施の形態では、原画像から取得した第１の高周波成分を解析することで、オーバーシュート領域を検出し、そのオーバーシュート領域の第１の高周波成分を、オーバーシュートが生じにくい高周波成分（第２の高周波成分）に置換するように構成したので、オーバーシュートを抑制した良好な処理後画像を提供できる。
【０１１８】
また、オーバーシュート領域に対してのみ、その第１の高周波成分を、より高周波領域にかたよった高周波成分（マスクサイズｄ２を小とした第２の高周波成分）に置き換えるように構成することで、微細構造の復元性を落とすことなくオーバーシュートを抑制できる。
【０１１９】
さらに、第２の高周波成分を生成する方法として、モーフォロジカル・フィルタ或いは中間値フィルタ等を用いた方法を用い、オーバーシュート領域のみに対して当該処理を行うように構成することで、さらなるオーバーシュートの抑制を期待できると共に、対象画像全体にわたってモーフォロジカル・フィルタ或いは中間値フィルタ等を施すように構成した場合に比較して、処理時間の短縮を行える。
【０１２０】
[第２の実施の形態]
本実施の形態では、上記図１に示したＸ線撮影装置１００の動作を、例えば、図３のフローチャートに従った動作とする。
【０１２１】
尚、上記図３のフローチャートにおいて、上記図２のフローチャートと同様の処理実行するステップには同じ符号を付し、その詳細な説明は省略する。
【０１２２】
ステップＳ２００〜ステップＳ２０３：
第１の実施の形態と同様にして、Ｘ線撮影により得られた画像は、前処理回路１０６で処理され、当該処理後の画像ｆ（ｘ、ｙ）は、画像処理回路１１１へ転送される（ステップＳ２００）。
画像処理回路１１１は、画像ｆ（ｘ、ｙ）に対して階調変換処理を施すことで原画像（元画像）ｆ０（ｘ、ｙ）を取得し（ステップＳ２０１）、第１の高周波成分作成回路１１２により、原画像ｆ（ｘ、ｙ）から第１の高周波成分ｈ１（ｘ、ｙ）を取得する（ステップＳ２０２）、領域検出回路１１３により、第１の高周波成分ｈ１（ｘ、ｙ）からオーバーシュート領域を検出する（ステップＳ２０３）。
【０１２３】
尚、ステップＳ２０３は必須ではなく、後述のようにオーバーシュート領域に対応した領域のみ低周波画像を作成する場合には必要であるが、それ以外の場合には必ずしも必要ではない。
【０１２４】
ステップＳ３０１：
第２の高周波成分作成回路１１４は、第１の高周波成分ｈ１（ｘ、ｙ）をダイナミックレンジ圧縮した結果を、第２の高周波成分ｈ２（ｘ、ｙ）として取得する。
【０１２５】
上記のダイナミックレンジ圧縮の処理方法としては、任意の方法を適用可能であるが、本実施の形態では、その一例として、先に説明した特開平１０−２７２２８３号等に記載の方法を用いる。
したがって、ここでは、第１の高周波成分ｈ１（ｘ、ｙ）の絶対値が、一定閾値Ｔｈ以上すなわちオーバーシュート領域のみのダイナミックレンジが圧縮されるように、ダイナミックレンジの圧縮処理を行う。このときの圧縮率ついては任意であるが、本実施の形態では、その一例として、最大の圧縮率とする。これにより、一定閾値Ｔｈ以上の第１の高周波成分ｈ１（ｘ、ｙ）が可能な限り圧縮されることになる。例えば、第１の高周波成分ｈ１（ｘ、ｙ）中の所定の低周波成分が概略一定閾値Ｔｈ以下にクリッピングされるように第１の高周波成分ｈ１（ｘ、ｙ）をダイナミックレンジ圧縮する。
【０１２６】
ここで、ダイナミックレンジ圧縮に用いる低周波画像（平滑化画像）ｆus３（ｘ、ｙ）について、その取得方法としては、例えば、次のような方法（１）〜方法（３）が挙げられる。
【０１２７】
方法（１）：オーバーシュート領域の座標に基づいた方法
本方法（１）では、メインメモリ１０９に保存されたオーバーシュート領域の座標に基づいて、
【０１２８】
【数１７】
なる式（２４）により、低周波画像ｆus３（ｘ、ｙ）を取得する。
上記式（２４）における“ｄ２”は、マスクサイズであり、上記式（１４）でのマスクサイズｄ１に対して、
ｄ２＜ｄ１
なる関係を有する。
【０１３０】
方法（２）：モーフォロジカル・フィルタを用いた方法
本方法（２）では、次の式（２５）〜（２８）により表されるモーフォロジカル・フィルタにより、低周波画像ｆus３（ｘ、ｙ）を生成する。
【０１３１】
【数１８】
上記式（２５）〜（２８）において、“D（ｘ、ｙ）”は、次の式（２９）により表される円盤状フィルタであり、“ｒ１”を任意の定数として、対象画像に応じて選択されるものである。
【０１３３】
【数１９】
方法（３）：中間値フィルタを用いた方法
本方法（３）では、中間値フィルタを用いて低周波画像ｆus３（ｘ、ｙ）を生成する。
具体的には、座標（ｘ、ｙ）を中心として、一定距離ｄ２内にある元画像ｆ（ｘ、ｙ）をソートし、中間にある画素値をfus３（ｘ、ｙ）とする。
【０１３５】
上述のような方法（２）或いは方法（３）を用いて得られる低周波画像ｆus３（ｘ、ｙ）は、エッジ構造を保存するものであり、この低周波画像ｆus３（ｘ、ｙ）を用いて得られた第２の高周波成分ｈ２（ｘ、ｙ）は、オーバーシュートのもとになる成分がでにくいものである。また、オーバーシュートが生じやすい領域のみに対して、モーフォロジカル・フィルタ或いは中間値フィルタを用いて取得した低周波画像ｆus３（ｘ、ｙ）を生成することで、対象画像全体に対してモーフォロジカル・フィルタ或いは中間値フィルタを施すように構成した場合よりも、処理時間の短縮を図ることができる。
【０１３６】
ステップＳ３０２：
処理後画像作成回路１１６は、
【０１３７】
【数２０】
なる式（３０）〜式（３２）により、処理後画像ｆｄ（ｘ、ｙ）を取得する。
【０１３９】
上記式（３０）〜式（３２）において、“ｈ１（ｘ、ｙ）”は、第２の高周波成分作成回路１１４によりダイナミックレンジ圧縮された後の第１の高周波成分（第２の高周波成分）である。また、上記式（３０）及び（３１）は、鮮鋭化処理を示し、上記式（３２）は、ダイナミックレンジの圧縮処理を示す。また、上記式（３２）において、“Ｆ２（）”は、低周波画像を階調変換する関数を示す。
【０１４０】
ステップＳ２０７：
上述のようにして、画像処理回路１１１により得られた処理後画像ｆｄ（ｘ、ｙ）は、出力回路１１８により、表示出力或いはフィルム出力等される。
【０１４１】
上述のように、本実施の形態では、第２の高周波成分として、第１の高周波成分をダイナミックレンジ圧縮した成分を用いるように構成した。つまり、オーバーシュートが生じやすい領域に対応する第１の高周波成分中の所定の低周波成分を圧縮するように構成したため、微細構造の復元性を損なわずにオーバーシュートを抑制することができる。また、第１の高周波成分の圧縮方法として、所定の低周波成分を圧縮すると共に所定の高周波成分を調節可能なダイナミックレンジ圧縮方法を用いる構成としたので、第１の高周波成分を自在に変更でき、もって微細構造の復元性を損なわずにオーバーシュートを抑制可能な第２の高周波成分を得ることができる。さらに、オーバーシュートが生じやすい領域に対してのみ、モーフォロジカルフィルタまたは中間値フィルタ等を用いてダイナミックレンジ圧縮に用いる低周波画像を作成するように構成すると、効果的にオーバーシュートを抑制することができると共に、第１の高周波成分画像全体に対しモーフォロジカルフィルタまたは中間値フィルタ等を用いる場合に比較して処理時間の短縮を図ることができる。
【０１４２】
また、ダイナミックレンジ圧縮に用いる低周波画像のマスクサイズｄ２を、ｄ１より小さく設定するように構成することで、オーバーシュートを抑制することができる。
【０１４３】
また、第１の高周波成分を、オーバーシュートが生じやすい領域において、より高周波領域にかたよった高周波成分（マスクサイズｄ２を小とした第２の高周波成分）に置き換えるように構成することで、微細構造の復元性を落とすことなくオーバーシュートを抑制できる。
【０１４４】
尚、本発明の目的は、第１及び第２の実施の形態のホスト及び端末の機能を実現するソフトウェアのプログラムコードを記憶した記憶媒体を、システム或いは装置に供給し、そのシステム或いは装置のコンピュータ（又はＣＰＵやＭＰＵ）が記憶媒体に格納されたプログラムコードを読みだして実行することによっても、達成されることは言うまでもない。
この場合、記憶媒体から読み出されたプログラムコード自体が第１及び第２の実施の形態の機能を実現することとなり、そのプログラムコード、及び当該プログラムコードを記憶した記憶媒体は本発明を構成することとなる。
プログラムコードを供給するための記憶媒体としては、ＲＯＭ、フロッピーディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＣＤ−Ｒ、磁気テープ、不揮発性のメモリカード等を用いることができる。
また、コンピュータが読みだしたプログラムコードを実行することにより、第１及び第２の実施の形態の機能が実現されるだけでなく、そのプログラムコードの指示に基づき、コンピュータ上で稼動しているＯＳ等が実際の処理の一部又は全部を行い、その処理によって第１及び第２の実施の形態の機能が実現される場合も含まれることは言うまでもない。
さらに、記憶媒体から読み出されたプログラムコードが、コンピュータに挿入された拡張機能ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに書き込まれた後、そのプログラムコードの指示に基づき、その機能拡張ボードや機能拡張ユニットに備わるＣＰＵなどが実際の処理の一部又は全部を行い、その処理によって第１及び第２の実施の形態の機能が実現される場合も含まれることは言うまでもない。
【０１４５】
【発明の効果】
以上説明したように本発明では、原画像（放射線撮影により得られた画像等）において、第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を前記原画像から取得し、その第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を第１の高周波成分とするように構成した。具体的には例えば、第２の高周波成分を、オーバーシュート領域に基づいて、オーバーシュートが生じにくい高周波成分として取得し、この第２の高周波成分で、オーバーシュート領域に対応する第１の高周波成分を置換する。これにより、オーバーシュートを抑制した良好な対象画像の処理後画像（表示出力やフィルム出力用の画像等）を提供できる。
【０１４６】
また、オーバーシュート領域のみに対して、第２の高周波成分を、第１の高周波成分より高周波領域にかたよった高周波成分（第１の高周波成分から所定の低周波成分を低減した高周波成分等）とするように構成した場合、対象画像における微細構造の復元性を落とすことなく、オーバーシュートを抑制することができる。
【０１４７】
また、第２の高周波成分を、第１の高周波成分をダイナミックレンジ圧縮した成分とするように構成した場合、対象画像における微細構造の復元性を維持した状態で、オーバーシュートを抑制することができる。
また、第１の高周波成分の圧縮として、第１の高周波成分中の所定の低周波成分を圧縮し、所定の高周波成分を調節可能なダイナミックレンジ圧縮を用いるように構成した場合、オーバーシュートを抑制することができると共に、良好な処理後画像を取得することができる。
【０１４９】
また、第２の高周波成分の取得に用いる平滑化画像を、中間値フィルタ（中間値演算）やモーフォロジカル・フィルタ（モーフォロジ演算）等のオーバーシュートの生じにくいフィルタを用いて取得するように構成した場合、オーバーシュートの抑制された良好な処理後画像を取得することができる。
【図面の簡単な説明】
【図１】第１の実施の形態において、本発明を適用したＸ線撮影装置の構成を示すブロック図である。
【図２】上記Ｘ線撮影装置の動作を説明するためのフローチャートである。
【図３】第２の実施の形態における、上記Ｘ線撮影装置の動作を説明するためのフローチャートである。
【符号の説明】
１００ Ｘ線撮影装置
１０１ Ｘ線発生回路
１０２ Ｘ線ビーム
１０３ 被写体
１０４ ２次元Ｘ線センサ
１０５ データ収集回路
１０６ 前処理回路
１０７ ＣＰＵバス
１０８ ＣＰＵ
１０９ メインメモリ
１１０ 操作パネル
１１１ 画像処理回路
１１２ 第１の高周波成分作成回路
１１３ 領域検出回路
１１４ 第２の高周波成分作成回路
１１５ 周波数成分置換回路
１１６ 処理後画像作成回路
１１８ 出力回路
Claims (13)
- 原画像から第１の高周波成分を取得する第１の高周波成分取得手段と、
前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を前記原画像から取得する第２の高周波成分取得手段と、
前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とする手段と、
前記高周波成分の大きさを所定の規則で変更した後の高周波成分と階調変換した後の前記原画像とを加算して処理後の画像を取得する処理後画像取得手段と、を備えることを特徴とする画像処理装置。 - 原画像から第１の高周波成分を取得する第１の高周波成分取得手段と、
前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を前記原画像から取得する第２の高周波成分取得手段と、
前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とする手段と、
前記高周波成分の大きさを所定の規則で変更した後の高周波成分と前記原画像とを加算して処理後の画像を取得する処理後画像取得手段と、を備えることを特徴とする画像処理装置。 - 原画像から第１の高周波成分を取得する第１の高周波成分取得手段と、
前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を前記原画像から取得する第２の高周波成分取得手段と、
前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とし、
前記高周波成分の大きさを所定の規則で変更した後の高周波成分と前記原画像の低周波成分とを加算して処理後の画像を取得する処理後画像取得手段と、を備えることを特徴とする画像処理装置。 - 原画像から第１のフィルタサイズでフィルタ処理して第１の高周波成分を取得する第１の高周波成分取得手段と、
前記原画像から前記第１のフィルタサイズよりも小さい第２のフィルタサイズでフィルタ処理して第２の高周波成分を取得する第２の高周波成分取得手段と、
前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とし、
前記高周波成分と前記原画像の所定の低周波成分とを加算して処理後の画像を取得する処理後画像取得手段と、を備えることを特徴とする画像処理装置。 - 複数の機器が互いに通信可能に接続されてなる画像処理システムであって、
前記複数の機器のうち少なくとも１つの機器は、請求項１〜４の何れか１項に記載の画像処理装置の機能を有することを特徴とする画像処理システム。 - 原画像から第１の高周波成分を取得する第１の高周波成分取得ステップと、
前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を前記原画像から取得する第２の高周波成分取得ステップと、
前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とするステップと、
前記高周波成分の大きさを所定の規則で変更した後の高周波成分と階調変換した後の前記原画像とを加算して処理後の画像を取得する処理後画像取得ステップと、を含むことを特徴とする画像処理方法。 - 原画像から第１の高周波成分を取得する第１の高周波成分取得ステップと、
前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を原画像から取得する第２の高周波成分取得ステップと、
前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とするステップと、
前記高周波成分の大きさを所定の規則で変更した後の高周波成分と前記原画像とを加算して処理後の画像を取得する処理後画像取得ステップと、を含むことを特徴とする画像処理方法。 - 原画像から第１の高周波成分を取得する第１の高周波成分取得ステップと、
前記第１の高周波成分と比較して低周波帯域の成分が減衰した第２の高周波成分を原画像から取得する第２の高周波成分取得ステップと、
前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とし、
前記高周波成分の大きさを所定の規則で変更した後の高周波成分と前記原画像の低周波成分とを加算して処理後の画像を取得する処理後画像取得ステップと、を備えることを特徴とする画像処理方法。 - 原画像から第１のフィルタサイズでフィルタ処理して第１の高周波成分を取得する第１の高周波成分取得ステップと、
前記原画像から前記第１のフィルタサイズよりも小さい第２のフィルタサイズでフィルタ処理して第２の高周波成分を取得する第２の高周波成分取得ステップと、
前記第１の高周波成分の大きさに基づいて定められる領域に対応する高周波成分を第２の高周波成分とし、前記領域と異なる領域の高周波成分を前記第１の高周波成分とし、
前記高周波成分と前記原画像の所定の低周波成分とを加算して処理後の画像を取得する処理後画像取得ステップと、を備えることを特徴とする画像処理方法。 - 請求項１〜４の何れか１項に記載の画像処理装置の機能、又は請求項５記載の画像処理システムの機能をコンピュータに実現させるためのプログラムを記録したコンピュータ読取可能な記憶媒体。
- 請求項６〜９の何れか１項に記載の画像処理方法の処理ステップをコンピュータに実行させるためのプログラムを記録したコンピュータ読取可能な記憶媒体。
- 請求項１〜４の何れか１項に記載の画像処理装置の機能、又は請求項５記載の画像処理システムの機能をコンピュータに実現させるためのプログラム。
- 請求項６〜９の何れか１項に記載の画像処理方法の処理ステップをコンピュータに実行させるためのプログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001034121A JP4612774B2 (ja) | 2001-02-09 | 2001-02-09 | 画像処理装置、画像処理システム、画像処理方法、記憶媒体、及びプログラム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001034121A JP4612774B2 (ja) | 2001-02-09 | 2001-02-09 | 画像処理装置、画像処理システム、画像処理方法、記憶媒体、及びプログラム |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2010111916A Division JP2010179154A (ja) | 2010-05-14 | 2010-05-14 | 画像処理装置、画像処理システム、画像処理方法、及び記憶媒体 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002237948A JP2002237948A (ja) | 2002-08-23 |
| JP4612774B2 true JP4612774B2 (ja) | 2011-01-12 |
Family
ID=18897777
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001034121A Expired - Fee Related JP4612774B2 (ja) | 2001-02-09 | 2001-02-09 | 画像処理装置、画像処理システム、画像処理方法、記憶媒体、及びプログラム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4612774B2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010179154A (ja) * | 2010-05-14 | 2010-08-19 | Canon Inc | 画像処理装置、画像処理システム、画像処理方法、及び記憶媒体 |
| JP2012033182A (ja) * | 2011-09-21 | 2012-02-16 | Canon Inc | 画像処理装置、画像処理システム、画像処理方法、及び記憶媒体 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4667072B2 (ja) * | 2005-02-25 | 2011-04-06 | 富士フイルム株式会社 | 画像処理装置、画像処理方法およびプログラム |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000307868A (ja) * | 1999-04-23 | 2000-11-02 | Fuji Photo Film Co Ltd | 画像処理方法および装置 |
-
2001
- 2001-02-09 JP JP2001034121A patent/JP4612774B2/ja not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000307868A (ja) * | 1999-04-23 | 2000-11-02 | Fuji Photo Film Co Ltd | 画像処理方法および装置 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010179154A (ja) * | 2010-05-14 | 2010-08-19 | Canon Inc | 画像処理装置、画像処理システム、画像処理方法、及び記憶媒体 |
| JP2012033182A (ja) * | 2011-09-21 | 2012-02-16 | Canon Inc | 画像処理装置、画像処理システム、画像処理方法、及び記憶媒体 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2002237948A (ja) | 2002-08-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3833177B2 (ja) | 画像処理装置、画象処理方法、記憶媒体及びプログラム | |
| US7418122B2 (en) | Image processing apparatus and method | |
| JP2018033745A (ja) | 画像処理装置、方法およびプログラム | |
| US6678400B1 (en) | Image processing apparatus, method and memory medium | |
| JP2006087934A (ja) | 画像処理装置及び画像処理方法 | |
| JP5284434B2 (ja) | Ｘ線画像処理装置、ｘ線撮影装置、ｘ線画像処理方法及びプログラム | |
| JP4612774B2 (ja) | 画像処理装置、画像処理システム、画像処理方法、記憶媒体、及びプログラム | |
| JP4208333B2 (ja) | 画像処理装置、画像処理システム、画像処理方法、及び記憶媒体 | |
| JP4497756B2 (ja) | 画像処理装置、画像処理システム、画像処理方法、記憶媒体、及びプログラム | |
| JP2005109867A (ja) | 画像処理方法および画像処理装置ならびに画像処理プログラム | |
| JP3937660B2 (ja) | 画像処理方法、及び、コンピュータ読み取り可能な記憶媒体 | |
| JP4146958B2 (ja) | 画像処理方法、画像処理装置及び記憶媒体 | |
| JP5539438B2 (ja) | 画像処理装置、画像処理方法及びプログラム | |
| JP4323708B2 (ja) | 画像処理方法、画像処理装置および記録媒体 | |
| JP2013215593A (ja) | 画像処理装置、画像処理システム、画像処理方法、及び記憶媒体 | |
| JP2010179154A (ja) | 画像処理装置、画像処理システム、画像処理方法、及び記憶媒体 | |
| JP3937616B2 (ja) | 画像処理装置、方法及びコンピュータ読み取り可能な記憶媒体 | |
| JPH0738758A (ja) | 放射線画像のダイナミックレンジ圧縮方法 | |
| JP4810002B2 (ja) | 画像処理装置、画像処理システム、画像処理方法、記憶媒体、及びプログラム | |
| JP3780179B2 (ja) | 画像処理装置、画像処理システム、画像処理方法、記憶媒体、及びプログラム | |
| JP4669163B2 (ja) | 画像処理装置、画像処理システム、画像処理方法、記憶媒体、及びプログラム | |
| JP4822607B2 (ja) | 画像処理装置、画像処理システム、記憶媒体、プログラム及び画像処理方法 | |
| JP2002330953A (ja) | 画像処理装置、画像処理システム、画像処理方法、記憶媒体、及びプログラム | |
| JP4447804B2 (ja) | 画像処理装置、画像処理システム、画像処理方法、記憶媒体、及びプログラム | |
| JP2008212493A (ja) | 画像処理装置及び画像処理方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20071218 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20090701 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090929 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20091126 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100316 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100514 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100713 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100909 |
|
| 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: 20101005 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20101016 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131022 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |