JP3937660B2

JP3937660B2 - 画像処理方法、及び、コンピュータ読み取り可能な記憶媒体

Info

Publication number: JP3937660B2
Application number: JP25418499A
Authority: JP
Inventors: 弘之新畠
Original assignee: Canon Inc
Current assignee: Canon Inc
Priority date: 1998-09-18
Filing date: 1999-09-08
Publication date: 2007-06-27
Anticipated expiration: 2019-09-08
Also published as: JP2000155838A

Description

【０００１】
【発明の属する技術分野】
本発明は、Ｘ線画像等の画像のダイナミックレンジ圧縮処理機能を有する画像処理方法及びコンピュータ読み取り可能な記憶媒体に関するものである。
【０００２】
【従来の技術】
例えば、Ｘ線胸部画像は、Ｘ線が透過しやすい肺野の画像及びＸ線が非常に透過しにくい縦隔部の画像より構成されるため、画素値の存在するレンジが非常に広い。このため、肺野及び縦隔部の両方を同時に観察することが可能なＸ線胸部画像を得ることは困難であるとされてきた。
【０００３】
そこで、この問題を回避する方法として、「ＳＰＩＥＶｏｌ．６２６ＭｅｄｉｃｉｎｅＸＩＶ／ＰＡＣＳＩＶ（１９８６）」に記載された方法（以下、「方法１」と言う）がある。この方法１は、処理後の画素値Ｓ_D、オリジナル画素値（入力画素値）Ｓ_org、オリジナル画像（入力画像）の低周波数画像の画素値Ｓ_US、定数Ａ，Ｂ，Ｃ（例えば、Ａ＝２，Ｂ＝０．７）を持って、
Ｓ_D＝Ａ［Ｓ_org−Ｓ_US］＋Ｂ［Ｓ_US］＋Ｃ・・・・（１）
なる式（１）で表されるものである。
【０００４】
上述のような方法１では、高周波成分（第１項）、低周波成分（第２項）の重み付けを変えることが可能で、例えば、定数Ａ，ＢをＡ＝２，Ｂ＝０．７とした場合、高周波成分を強調し、且つ全体のダイナミックレンジを圧縮する効果が得られる。これは、ある放射線医師達によれば、本処理なしの画像と比較して、診断に有効であるという評価が得られている。
【０００５】
また、第２５０９５０３特許公報等には、処理後の画素値Ｓ_D、オリジナル画素値（入力画素値）Ｓ_org、オリジナル画像（入力画像）のＹ方向のプロファイルの平均プロファイルＰｙとＸ方向のプロファイルの平均プロファイルＰｘを持って、
Ｓ_D＝Ｓ_org＋Ｆ［Ｇ（Ｐｘ，Ｐｙ）］・・・・（２）
なる式（２）で表される方法（以下、「方法２」と言う）が記載されている。
【０００６】
ここで、式（２）での関数ｆ（ｘ）が有する特性について説明すると、まず、「ｘ＞Ｄｔｈ」ではｆ（ｘ）が「０」となり、「０≦ｘ≦Ｄｔｈ」ではｆ（ｘ）が切片を「Ｅ」、傾きを「Ｅ／Ｄｔｈ」として単調減少するものとなり、次の式（３）で表される。
ｉｆ（ｘ≦Ｄｔｈ）ｆ（ｘ）＝Ｅ−（Ｅ／Ｄｔｈ）ｘ
ｅｌｓｅｆ（ｘ）＝０・・・・（３）
また、”Ｐｘ”及び”Ｐｙ”は、
Ｐｙ＝（ΣＰｙｉ）／ｎ・・・・（４）
Ｐｘ＝（ΣＰｘｉ）／ｎ・・・・（５）
なる式（４）及び（５）で表される。
【０００７】
式（４）及び（５）において、”ｉ”は”ｉ＝１〜ｎ”であり、”Ｐｙｉ”及び”Ｐｘｉ”はプロファイルであり、また、式（２）での”Ｇ（Ｐｘ，Ｐｙ）”は、例えば、
Ｇ（Ｐｘ，Ｐｙ）＝ｍａｘ（Ｐｘ，Ｐｙ）・・・・（６）
で示される。
【０００８】
上述のような方法２では、低周波画像の画素値で”Ｄｔｈ”以下の濃度レンジが圧縮される。
【０００９】
また、方法２（第２５０９５０３特許公報等に記載の方法）と同様の方法（以下、「方法３」と言う）が「日本放射線技術学会雑誌第４５巻第８号１９８９年８月１０３０頁阿部ほか」に記載されている。この方法３は、処理後の画素値Ｓ_D、オリジナル画素値（入力画素値）Ｓ_org、オリジナル画像（入力画像）をマスクサイズＭ×Ｍ画素で移動平均をとった時の平均画素値Ｓ_US、単調減少関数ｆ（ｘ）を持って、
Ｓ_D＝Ｓ_org＋ｆ（Ｓ_US）・・・・（７）
Ｓ_US＝ΣＳ_org／Ｍ² ・・・・（８）
なる式（７）及び（８）で表されるものである。
【００１０】
方法３は、式（２）によって表される方法２とは低周波画像の作成方法が異なり、方法２では１次元データで低周波画像を作成していたのに対し、２次元データで低周波画像を作成する方法である。このような方法３によっても、低周波画像の画素値で”Ｄｔｈ”以下の濃度レンジが圧縮される。
【００１１】
また、第２６６３１８９特許公報等には、処理後の画素値Ｓ_D、オリジナル画素値（入力画素値）Ｓ_org、オリジナル画像（入力画像）をマスクサイズＭ×Ｍ画素で移動平均をとった時の平均画素値Ｓ_US、単調増加関数ｆ１（ｘ）を持って、
Ｓ_D＝Ｓ_org＋ｆ１（Ｓ_US）・・・・（７）
Ｓ_US＝ΣＳ_org／Ｍ² ・・・・（８）
なる式（７）及び（８）で表される方法（以下、「方法４」と言う）が記載されている。
【００１２】
ここで、式（７）での関数ｆ１（ｘ）が有する特性について説明すると、まず、「ｘ＜Ｄｔｈ」ではｆ（ｘ）が「０」となり、「Ｄｔｈ≦ｘ」ではｆ（ｘ）が切片を「Ｅ」、傾きを「Ｅ／Ｄｔｈ」として単調減少するものとなり、次の式（１１）で表される。
ｆ１［ｘ］＝Ｅ−（Ｅ／Ｄｔｈ）Ｘ・・・・（１１）
【００１３】
上述のような方法４では、低周波画像の画素値Ｄｔｈ以上の画素値（濃度値）が圧縮される。このときのアルゴリズムについては、方法３（「日本放射線技術学会雑誌第４５巻第８号１９８９年８月１０３０頁阿部ほか」等に記載の方法）と同様である。
【００１４】
【発明が解決しようとする課題】
しかしながら、上述したような従来の方法では、取得画像（入力画像）を階調変換してフィルムやＣＲＴに対して出力した場合、視認確認できない濃度値もしくは輝度値の範囲に階調変換後の注目領域（被写体領域等）の濃度値もしくは輝度値がかかるのを防ぐことを目的としていた。
【００１５】
このため、従来の方法では、圧縮開始位置や圧縮率が固定であり、ダイナミックレンジ圧縮を施した取得画像を階調変換した場合、注目領域の濃度値もしくは輝度値が視認領域に入らないことがあった。これは診断等に支障をきたす問題につながる。
【００１６】
また、取得画像によっては、被写体部分の圧縮が必要以上に行なわれてしまい、本来被写体部分が保持している情報量が無意味に損出される場合があった。これは、診断能等が低下する問題につながる。
【００１７】
さらに、ダイナミックレンジを拡大するという発想が無かったので、フィルムやＣＲＴ上に画像を出力する場合に、注目領域のダイナミックレンジを有効に出力することができないという問題があった。すなわち、注目領域のダイナミックレンジを所定の領域に変換するという思想がなかった。
【００１８】
本発明は上記のような問題を解決するためになされたもので、高周波成分の振幅を保った状態で、注目領域のダイナミックレンジを所定の濃度領域に変換できるようにすることを目的とする。
【００１９】
【課題を解決するための手段】
本発明の画像処理方法は、原画像データに対し諧調変換処理及び周波数処理を施す画像処理方法であって、前記原画像データにおける被写体内の領域から、階調変換曲線を生成するための値を抽出する工程と、前記原画像データに基づく高周波成分を生成する工程と、前記高周波成分を、前記階調変換曲線の傾きに基づいて変更する工程と、前記変更した高周波波成分に基づき画像データを生成する工程とを備えることを特徴とする。
本発明のコンピュータ読み取り可能な記憶媒体は、原画像データに対し諧調変換処理及び周波数処理を施す画像処理方法をコンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記憶媒体であって、前記原画像データにおける被写体内の領域から、階調変換曲線を生成するための値を抽出する工程と、前記原画像データに基づく高周波成分を生成する工程と、前記高周波成分を、前記階調変換曲線の傾きに基づいて変更する工程と、前記変更した高周波波成分に基づき画像データを生成する工程とをコンピュータに実行させるためのプログラムを記憶したことを特徴とする。
【００４３】
【発明の実施の形態】
以下、本発明の実施の形態を図面と共に説明する。
【００４４】
（第１の実施の形態）
図１は、本発明の実施の形態による画像処理装置１００を示す。
この画像処理装置１００は、濃度値変換機能を有するＸ線画像の画像処理装置であり、前処理回路１０６、濃度値変換部１１３，ＣＰＵ１０８、メインメモリ１０９、操作パネル１１０、画像表示装置１１１を備えており、ＣＰＵバス１０７を介して互いにデータ授受されるようになされている。
【００４５】
また、画像処理装置１００は、前処理回路１０６に接続されたデータ収集回路１０５と、データ収集回路１０５に接続された２次元Ｘ線センサ１０４及びＸ線発生回路１０１とを備えており、これらの各回路はＣＰＵバス１０７にも接続されている。
【００４６】
図２は画像処理装置の処理の流れを示すフローチャートである。
図３（ａ）は入力画像を示し、図３（ｂ）は上記入力画像からＸ線の照射領域とす抜け領域（Ｘ線が素通りした領域）とを削除した図である。
図４は階調変換曲線上で視認できる濃度域と注目領域の濃度域との関係を示す図であり、横軸が入力画像の濃度値、縦軸が出力画像の濃度値である。
【００４７】
上記構成による画像処理装置１００において、メインメモリ１０９は、ＣＰＵ１０８の処理に必要な各種のデータなどが記憶されると共に、ＣＰＵ１０８の作業用としてのワークメモリを含む。ＣＰＵ１０８は、メインメモリ１０９を用いて、操作パネル１１０からの操作に従って装置全体の動作制御等を行う。これにより画像処理装置１００は、以下のように動作する。
【００４８】
まず、Ｘ線発生回路１０１は、被検査体１０３に対してＸ線ビーム１０２を放射する。このＸ線ビーム１０２は、被検査体１０３を減衰しながら透過して、２次元Ｘ線センサ１０４に到達し、２次元Ｘ線センサ１０４よりＸ線画像として出力される。ここでは、２次元Ｘ線センサ１０４から出力されるＸ線画像を、例えば図３（ａ）のような膝、胸椎画像とする。
【００４９】
データ収集回路１０５は、２次元Ｘ線センサ１０４から出力されたＸ線画像を電気信号に変換して前処理回路１０６に供給する。前処理回路１０６は、データ収集回路１０５からの信号（Ｘ線画像信号）に対して、オフセット補正処理やゲイン補正処理等の前処理を行う。前処理が行われたＸ線画像信号は入力画像として、ＣＰＵ１０８の制御により、ＣＰＵバス１０７を介してメインメモリ１０９、照射領域抽出回路１１２及び濃度値変換部１１３のす抜け削除回路１１３ａに転送される。
【００５０】
濃度値変換部１１３において、１１３ａはす抜け領域とす抜け領域と一定幅で接する体領域とを削除するす抜け削除回路、１１３ｂはす抜け削除回路１１３ａで削除されなかった領域から、濃度値の最大値と最小値を算出する最高値最小値抽出回路、１１３ｃは入力画像の平滑化画像を作成する平滑化画像作成回路、１１３ｄは入力画像と平滑化画像との差分から高周波成分を作成する高周波成分作成回路、１１３ｅは最大値最小値抽出回路１１３ｂで抽出された最大値と最小値に基づき平滑化画像の濃度値を変換する濃度値変換回路、１１３ｆは濃度値変換回路１１３ｅで変換された画像に高周波成分作成回路１１３ｄで作成した高周波を足し込む高周波成分足し込み回路である。
【００５１】
次に濃度値変換部１１３の動作について図２のフローチャートにより説明する。
ＣＰＵバス１０７を介して前処理回路１０６で処理された入力画像をＣＰＵ１０８の制御により受信した照射領域抽出回路１１２は、入力画像中の照射領域を抽出する（ステップＳ２０１）。また、同時に入力画像を受信したす抜け削除回路１１３ａは照射領域外及び照射領域内のす抜け領域とす抜け領域と一定間隔内で接する体領域とを例えば０画素で置き換える（ステップＳ２０２）。具体的には次の式（１２）で表されるような画像の変換を行う。
【００５２】
【数１】

【００５３】
ここで、ｆ（ｘ，ｙ）は画像データを示し、ｆ１（ｘ，ｙ）はす抜け領域及びす抜け領域と一定間隔内で接する体領域を削除した後の画像を示す。Ｓｇｎ（ｘ，ｙ）は以下の式（４）のように表される。Ｔｈ１は実験により定められる定数、ｄ１，ｄ２は体領域を削除する幅を決める定数である。
ｓｇｎ（ｘ，ｙ）＝０ｆ（ｘ，ｙ）≧Ｔｈ１のとき
ｓｇｎ（ｘ，ｙ）＝１その他・・・・（１３）。
図３（ｂ）は、入力画像の照射領域外及びす抜け領域を０で置き換えた場合の画像である。
【００５４】
次に、最大最小値抽出回路１１３ｂでは、画像濃度値の最大値（図４の”ｄ２”）と最小値（図４の”ｄ３”）を算出する（ステップＳ２０３）。ここで、最大値、最小値を算出するのに、照射領域及びす抜け領域を除去した領域の平滑化画像から抽出してもよい。次に平滑化画像作成回路１１３ｃは次式に従い平滑化画像を作成する（ステップＳ２０４）。ｆｕｓ（ｘ，ｙ）は入力画像ｆ０（ｘ，ｙ）の平滑化（低周波）画像の画素値であり、例えば、式（１４）〜（１８）、又は式（１９）で示される。
【００５５】

【００５６】
”Ｄ（ｘ，ｙ）”は、円盤状フィルタであり、ｒ１を任意の定数とし、入力画像に応じて選択される。

【００５７】
ここで得られたｆｕｓ（ｘ，ｙ）のプロファイルはエッジ構造を保存しているものであり、従来ダイナミックレンジ圧縮の欠点であるオーバーシュート、アンダーシュートが起きないものである。
【００５８】
平滑化画像は、例えば式（１７）で示される平均濃度を用いてもよいし、次の式（１９）を用いてもよい。ここでｄは定数である。
【００５９】
【数２】

【００６０】
また、例えばＥｒｏｓｉｏｎ、Ｄａｉ１ａｔｉｏｎ、Ｏｐｅｎｉｎｇ、Ｃ１ｏｓｉｎｎｇ等のモルフォロジカルフィルタを用い平滑化画像を作成してもよい。
【００６１】
次に、高周波成分作成回路１１３ｃでは、入力画像ｆ０（ｘ，ｙ）と平滑化画像ｆｕｓ（ｘ，ｙ）とから次式に従い高周波画像ｆｈ（ｘ，ｙ）を作成する（ステップＳ２０５）。次に濃度値変換回路１１３ｅは、最大最小値値抽出回路１１３ｂで抽出した最大値（ｄ３）、最小値（ｄ２）、そして階調変換後の視認濃度値で決まる濃度値ｄ１，ｄ４に基づき次式により濃度値を変換し、濃度値変換した平滑化画像Ｓｕｓ０×（ｘ，ｙ）を作成する（ステップＳ２０６）。
【００６２】
図４において、ＳｉｇｈｔｍａｘとＳｉｇｈｔｍｉｎが視認できる最大最小濃度値であり、ｄ１，ｄ４に対応する。平滑化画像の濃度値がｍｉｎ（平滑化画像が取り得る最小値）からｄ１以下の場合、

【００６３】
ｄ２＜Ｓｕｓ（ｘ，ｙ）≦ｄ３の場合、

【００６４】
ｄ３＜Ｓｕｓ（ｘ，ｙ）≦ｍａｘの場合、

【００６５】
そして、濃度値変換した画像Ｓｕｓｕｏ（ｘ，ｙ）に高周波画像ｆｈ（ｘ，ｙ）を足し込み、最終画像ｆｅ（ｘ，ｙ）を得る（ステップＳ２０７）。
ｆｅ（ｘ，ｙ）＝Ｓｕｓｕ０（ｘ，ｙ）十ｆｈ（ｘ，ｙ）・・・（２３）
最終的には、得られた画像ｆｅ（ｘ，ｙ）を階調変換回路１１４で階調変換してフィルム出力又は画像表示装置１１１に表示する。
【００６６】
尚、本実施の形態では、す抜け領域がある画像に対する場合について説明したが、す抜け領域がない場合は、す抜け削除回路１１３ａを経由せず、照射領域内の画像から、最大最小値抽出回路１１３ｂにより最大値と最小値を抽出すればよい。
【００６７】
本実施の形態によれば、高周波成分を保持したまま、注目領域の画像を一定幅に圧縮、拡大することが可能であり、画像情報を有効に使える効果がある。さらに、上記一定幅をフィルムなどの視認領域とすれば、フィルム上の視認領域の幅に注目領域を拡大縮小することができる効果が有る。また、最大値、最小値を平滑化画像から抽出する場合は、濃度値変換部１１３での変換をより精度よく行える効果が有る。
【００６８】
（第２の実施の形態）
図５は、本発明の実施の形態による画像処理装置を適用したＸ線撮影装置６００の構成を示す。
【００６９】
Ｘ線撮影装置６００は、画像処理機能を有するものであり、前処理回路６０６、ＣＰＵ６０８、メインメモリ６０９、操作パネル６１０、画像表示器６１１、及び画像処理回路６１２を備えており、これらの各回路はＣＰＵバス６０７を介して互いにデータ授受されるようになされている。
【００７０】
また、Ｘ線撮影装置６００は、前処理回路６０６に接続されたデータ収集回路６０５と、データ収集回路６０５に接続された２次元Ｘ線センサ６０４及びＸ線発生回路６０１とを備えており、これらの各回路もＣＰＵバス６０７に接続されている。
【００７１】
画像処理回路６１２は、階調変換のための特徴量を算出する特徴抽出回路６１２ａと、後述する第１の階調変換回路６１２ｄの階調変換曲線からダイナミックレンジを変更する範囲と変更量を算出する制御回路６１２ｂと、制御回路６１２ｂで算出されたダイナミックレンジを変更する範囲と変更量から原画像のダイナミックレンジを変更するＤＲＣ回路６１２ｃと、ＤＲＣ回路６１２ｃにてダイナミックレンジが変更された原画像の階調変換を行なう第１の階調変換回路６１２ｄとを備えている。
【００７２】
ＤＲＣ回路６１２ｃは、原画像の階調変換を行なう第２の階調変換回路６１３と、第２の階調変換回路６１３にて階調変換された画像に対して原画像の高周波成分を足し込む高周波成分調整回路とを備えている。
【００７３】
図６は、例えば、頸椎側面で撮影を行なった場合に得られる画像（Ｘ線撮影画像）において、特徴量を抽出する領域ａを示したものである。
【００７４】
図７は、第１の階調変換回路６１２ｄの階調変換曲線を示したものであり、この図７において、横軸は画素値を示し、縦軸は濃度値を示す。”Ｓｌ”は、濃度上で視認識限界の下限濃度値に対応する画素値を示し、”Ｓｈ”は、濃度上で視認識限界の上限濃度値に対応する画素値を示す。”Ｓｍｉｎ”及び”Ｓｍａｘ”は、注目領域の最小画素値及び最大画素値を示す。”Ｓａ”は、特徴抽出回路６１２ａにて算出された特徴量に対応する画素値を示し、”ｄｌ”及び”ｄｈ”はそれぞれ、ダイナミックレンジの変更しない範囲の下限濃度値及び上限濃度値に対応する画素値を示す。
【００７５】
図８は、第２の階調変換回路６１３の階調変換曲線を示したものであり、この図８において、横軸は入力画素値を示し、縦軸は出力画素値を示す。”ｄｌ”、”ｄｈ”、及び”Ｓａ”は、図７での画素値と同じ画素値を示す。
図８に示す階調変換曲線ａは、直線成分のみから構成され、同図に示す階調変換曲線ｂは、階調変換曲線ａの微分値が連続となるようななだらかな曲線から構成されている。
【００７６】
ここで、上述のようなＸ線撮影装置６００において、メインメモリ６０９は、ＣＰＵ６０８による本装置全体の動作制御のための処理に必要な処理プログラムや各種のデータ等が記憶されるものであると共に、ＣＰＵ６０８の作業用メモリとしてのワークメモリを含むものである。
したがって、ＣＰＵ６０８は、メインメモリ６０９及びメインメモリ６０９に記憶された情報を用いて、操作パネル６１０からの操作に従った本装置全体の動作制御等を行なう。これにより、Ｘ線撮影装置６００は次のように動作する。
【００７７】
先ず、Ｘ線発生回路６０１が被写体６０３に対してＸ線ビーム６０２を放射する。Ｘ線発生回路６０１から放射されたＸ線ビーム６０２は、被写体６０３を減衰しながら透過して、２次元Ｘ線センサ６０４に到達し、２次元Ｘ線センサ６０４によりＸ線画像として出力される。この２次元Ｘ線センサ６０４から出力されるＸ線画像は、例えば、人体部画像である。
【００７８】
データ収集回路６０５は、２次元Ｘ線センサ６０４から出力されたＸ線画像を電気信号に変換して前処理回路６０６に供給する。前処理回路６０６は、データ収集回路６０５からの信号（Ｘ線画像信号）に対して、オフセット補正処理やゲイン補正処理等の前処理を施す。この前処理回路６０６にて前処理が施されたＸ線画像信号は入力画像として、ＣＰＵ６０８の制御により、ＣＰＵバス６０７を介して、メインメモリ６０９及び画像処理回路６１２にそれぞれ転送される。
【００７９】
画像処理回路６１２は、例えば、図９に示すフローチャートに従って次のように動作する。
【００８０】
先ず、ＣＰＵバス６０７を介して前処理回路６０６から供給された入力画像（前処理が施された原画像）ｆ０（ｘ，ｙ）を、ＣＰＵ６０８の制御により受信した特徴抽出回路６１２ａは、階調変換のための特徴量Ｓ_aを算出する（ステップＳ７０１）。
ここでの特徴量の抽出方法は、入力画像ｆ０（ｘ，ｙ）での被写体の部位毎に異なり、本出願人から複数の方法が提案されている。例えば、入力画像ｆ０（ｘ，ｙ）が、頸椎の撮影により得られた画像である場合には、特願平１０−２７２２８４号に記載された方法を用いる。この方法は、特徴量Ｓａとして、図６の”ａ”に示した領域内の画素値の平均値を算出する方法である。
【００８１】
次に、第１の階調変換回路６１２ｄでの階調変換曲線を、特徴抽出回路６１２ａにて得られた特徴量Ｓａを用いて規定する（ステップＳ７０２）。
例えば、図７に示すように、特徴量Ｓａの画素値が濃度値１．０に変更されるような階調変換曲線を規定する。
【００８２】
次に、制御回路６１２ｂは、ステップＳ７０２にて規定された階調変換曲線から入力画像ｆ０（ｘ，ｙ）のダイナミックレンジを変更する範囲と変更量を算出する（ステップＳ７０３）。この算出方法としての一例を次に述べる。
【００８３】
例えば、通常の日本人は、画像での視認識できるの濃度値が、濃度値０．２から３．０の範囲となっている。
そこで、先ず、ステップＳ７０２にて規定された階調変換曲線に基づき、濃度値０．２、３．０に階調変換される画素値Ｓｌ、Ｓｈを算出する。
次に、ダイナミックレンジを変更しない幅を、
ｄｌ＝Ｓａ−ｈｌ・・・・（２４）
ｄｈ＝Ｓａ＋ｈｈ・・・・（２５）
なる式（２４）及び（２５）に従って算出する。
式（２４）及び（２５）において、”ｈｌ”及び”ｈｈ”は定数であり、”ｄｌ”及び”ｄｈ”はダイナミックレンジを変更する起点となる画素である。
【００８４】
そして、”ｄｌ”以下の画素値、”ｄｌ”以上の画素値の変更量Ｒｌ，Ｒｈを、
Ｒｌ＝（Ｓｌ−ｄｌ）／（Ｓｍｉｎ−ｄｌ）・・・・（２６）
Ｒｈ＝（Ｓｈ−ｄｈ）／（Ｓｍａｘ−ｄｈ）・・・・（２７）
なる式（２６）及び（２７）に従って算出する。
式（２６）及び（２７）において、”Ｓｍｉｎ”及び”Ｓｍａｘ”は、被写体領域（注目領域）の画素値のうちの最小値及び最大値を示す。すなわち、ｄｌ以下の画素値のダイナミックレンジをＲｌ倍すると、画素値Ｓｍｉｎは画素値Ｓｌに変更される。同様に、画素値Ｓｍａｘは画素値Ｓｈに変更される。換言すれば、起点をｄｌ，ｄｈとし、変更量をＲｌ，Ｒｈとすると、階調変換後の被写体領域の濃度値は、視野領域（ここでは例えば、濃度値０．２〜３．０）の幅と一致する。
【００８５】
上述のように、変更量Ｒｌ，Ｒｈを求めることにより、ダイナミックレンジを圧縮しても、撮影画像の注目領域の全てを出力画像上に再現することができる。したがって、特徴抽出の不安定性を吸収することができ、診断能力の向上を図ることができる。
【００８６】
また、変更量Ｒｌ，Ｒｈに対する拘束条件として、
Ｒｌ＜Ｃｌ・・・・（２８）
Ｒｈ＜Ｃｈ・・・・（２９）
なる式（２８）及び（２９）に示すような条件を設けるようにしてもよい。
式（２８）及び（２９）において、”Ｃｌ”及び”Ｃｈ”はそれぞれ定数を示す。これらの定数Ｃｌ，Ｃｈ、及び上述した定数ｈｌ，ｈｈは、例えば、撮影対象となる被写体の部位毎に実験的に決まる定数である。
このように、拘束条件を設けることによって、圧縮のしすぎを防ぐことができるため、圧縮のしすぎによる診断能力の低下を確実に防止することができる。
【００８７】
上述のようにして、ダイナミックレンジを変更する起点となる画素ｄｌ，ｄｈ、及び変更量Ｒｌ，Ｒｈが求まると、これらのｄｌ，ｄｈ，Ｒｌ，Ｒｈに基づき、図８に示したような第２の階調変換回路６１３の階調変換曲線を規定する。
図８に示す階調変換曲線ａでの画素値ｄｌ以下の傾きは”Ｒｌ”であり、画素値ｄｌから画素値ｄｈの範囲の傾きは”１”であり、画素値ｄｈ以上の傾きは”Ｒ_h”である。一方、同じく図８に示す階調変換曲線ｂは、特願昭１１−７６８８２号等に記載されている技術によって、階調変換曲線ａの微分値が連続となるようになだらかにした曲線である。第２の階調変換回路６１３の階調変換曲線として、階調変換曲線ａを用いる場合には、微分不連続点で偽輪郭が発生することがあるが、階調変換曲線ｂを用いる場合には、当該偽輪郭が発生することはない。ここでは一例として、階調変換曲線ｂを第２の階調変換回路６１３で用いるものとする。
【００８８】
次に、ＤＲＣ回路６１２ｃは、以降の処理で、高周波成分を保持するか否かを選択する（ステップＳ７０４）。
高周波成分を保持せずに処理実行した場合に得られる結果画像は、低周波成分が重要である軟部組織等を診断する場合に用いて有効である。一方、高周波成分を保持して処理実行した場合に得られる結果画像は、高周波成分が重要である骨部、肺野等を診断する場合に用いて有効である。すなわち、処理対象の画像が軟部組織等の画像であれば、高周波成分を保持しないと選択し、処理対象の画像が骨部、肺野等の画像であれば、高周波成分を保持すると選択することによって、軟部組織や、骨部、肺野等のそれぞれの診断に有効な結果画像が得られる。
このステップＳ７０４での選択は、撮影部位（軟部組織、骨部、肺野等）に基づき自動選択する構成としてもよいし、ユーザがマニュアル選択する構成としてもよい。
【００８９】
ステップＳ７０４での選択の結果、高周波成分を保持しない場合、第２の階調変換回路６１３により、入力画像（原画像）に対して階調変換曲線ｂ（図８参照）を用いた階調変換を施す（ステップＳ７０８）。
そして、第１の階調変換回路６１２ｄにより、第２の階調変換回路６１３での階調変換後の入力画像に対して階調変換を施す（ステップＳ７０９）。
【００９０】
ここで、入力画像（原画像）の座標を”ｘ，ｙ”、原画像の画素値を”ｆ０（ｘ，ｙ）”、第２の階調変換回路６１３の階調変換曲線ｂを”Ｆ（ｘ）”、階調変換曲線ｂを用いた階調変換後の画像の画素値を”ｆ１（ｘ，ｙ）”、第１の階調変換回路６１２ｄの階調変換曲線を”Ｆ１（ｘ）”とすると、結果的に得られる画像（最終画像）ｆ２（ｘ，ｙ）は、
ｆ２（ｘ，ｙ）＝Ｆ１（Ｆ（ｆ０（ｘ，ｙ））・・・・（３０）
なる式（３０）で表される。
この場合、”Ｆ１（Ｆ（ｘ））を、１つの階調変換曲線と考えてよい。
【００９１】
一方、ステップＳ７０４での選択の結果、高周波成分を保持する場合、先ず、高周波成分調整回路６１４により、原画像ｆ０（ｘ，ｙ）の高周波成分ｆｈ（ｘ，ｙ）を、
ｆｈ（ｘ，ｙ）＝ｆ０（ｘ，ｙ）−ｆｕｓ（ｘ，ｙ）・・・・（３１）
なる式（３１）により算出する（ステップＳ７０５）。
式（３１）において、”ｆｕｓ（ｘ，ｙ）”は、平滑化画像であり、マスクサイズの大きさを示す定数ｄを持って、
【００９２】
【数３】

【００９３】
なる式（３２）により算出される。
【００９４】
次に、第２の階調変換回路６１３により、原画像ｆ０（ｘ，ｙ）に対して、
ｆ１（ｘ，ｙ）＝Ｆ（ｆ０（ｘ，ｙ））・・・・（３３）
なる式（３３）で示される階調変換を施し、当該階調変換後の画像ｆ１（ｘ，ｙ）を取得する（ステップＳ７０６）。
【００９５】
次に、高周波成分調整回路６１４により、第２の階調変換回路６１３の階調変換曲線ｂの微係数から係数ｃ（ｘ）を、
【００９６】
【数４】

【００９７】
なる式（３４）により算出する。この式（３４）に示されるように、係数ｃ（ｘ）は、階調変換曲線ｂの傾きから”１”を引いた値である。
そして、係数ｃ（ｘ）に基づいて、ステップＳ７０６にて取得された第２の階調変換回路６１３での階調変換後の画像ｆ１（ｘ，ｙ）に対して、ステップＳ７０６にて取得された原画像ｆ０（ｘ，ｙ）の高周波成分ｆｈ（ｘ，ｙ）を足し込み、その処理後画像ｆ３（ｘ，ｙ）を得る。すなわち、処理後画像ｆ３（ｘ，ｙ）は、

なる式（３５）で表される（ステップＳ７０７）。
【００９８】
そして、第１の階調変換回路６１２ｄにより、ステップＳ７０７にて得られた処理後画像ｆ３（ｘ，ｙ）に対して、
ｆ２（ｘ，ｙ）＝Ｆ１（ｆ３（ｘ，ｙ））・・・・（３６）
なる式（３６）に従った階調変換を行い、この結果得られた画像ｆ２（ｘ，ｙ）を最終画像として得る（ステップＳ７０９）。
【００９９】
尚、ステップＳ７０５の式（３１）における平滑化画像ｆｕｓ（ｘ，ｙ）を、例えば、次の式（３７）〜（４１）で示されるモルフォジ演算を用いて算出するようにしてもよい。

式（３７）〜（４０）において、”Ｄ（ｘ１，ｙ１）”は、円盤状フィルタであり、入力画像に応じて選択される任意の定数ｒ１を持って、

なる式（４１）で表される。
【０１００】
式（３７）〜（４１）により得られる平滑化画像ｆｕｓ（ｘ，ｙ）のプロファイルは、エッジ構造を保存しているものであり、従来ダイナミックレンジ圧縮の欠点であるオーバーシュート、アンダーシュートが起きないものである。
【０１０１】
本実施の形態によれば、入力画像の任意の階調領域の濃度分布幅を圧縮、伸張することができ、且つ階調変換後の高周波成分の振幅を階調変換前の画像の高周波成分の振幅を保持できる効果がある。
また、階調変換曲線で決まる画素値に基づき、画像のダイナミックレンジを変更するので、階調変換後の画像の濃度値を想定して、ダイナミックレンジを変更できるので、階調変換後の濃度値の範囲を一意に調整できる効果がある。
【０１０２】
また、階調変換曲線から視認限界で決まる画素値を求め、その画素値の値に基づきダイナミックレンジを変更する範囲、量を決定できるので、注目領域の画素値の濃度値の範囲を視認範囲にすることができる効果がある。
また、注目領域が視認領域にあることから、注目領域全体を観察することが可能であり、診断能が上がる効果がある。
また、注目領域を視認領域の幅に一致することが可能であり、注目領域の情報量を濃度値として最大限に表現できる。
また、画像のダイナミックレンジを変更しない範囲を有するので、診断等で重要となる画素範囲をダイナミックレンジを変更しない領域とし、従来通りの濃度値で表現し、さらに、従来では視認できなかった領域も濃度値として観察できるため、診断能が上がる効果がある。
また、ダイナミックレンジの変更手段において高周波成分を保持するため、高周波成分の情報量を落とすことがなく、さらに従来では濃度値上で視認できなかった領域の濃度値も観察できるため、診断能が上がる効果がある。
【０１０３】
（他の実施の形態）
上述の実施の形態は、ハード的に実現してもよく、また、ソフト的に構成してもよい。例えば、ＣＰＵやメインメモリ等のメモリ等からなるコンピュータシステムを用いて、図２又は図９に示したフローチャートによる処理を実行するためのプログラムを実現することにより、上述の実施の形態における構成を実現してもよい。
【０１０４】
上記プログラムを記憶する記憶媒体としては、ＲＯＭ、ＲＡＭ等の半導体メモリ、光ディスク、光磁気ディスク、磁気記憶媒体等を用いてよく、これらをＣＤ−ＲＯＭ、ＦＤ、磁気カード、磁気テープ、不揮発性メモリカード等を用いることができる。
【０１０５】
また、コンピュータ上で稼働しているＯＳ等が処理の一部又は全部を行う場合、あるいは記憶媒体から読み出されたプログラムコードが、コンピュータに挿入された拡張機能ボードやコンピュータに接続された拡張機能ユニットに備わるメモリに書き込まれた後、そのプログラムコードの指示に基づいて、上記拡張機能ボードや拡張機能ユニットに備わるＣＰＵ等が処理の一部又は全部を行う場合にも、上記実施の形態と同等の機能を実現できると共に、同等の効果を得ることができ、本発明に含まれることは言うまでもない。
【０１０６】
【発明の効果】
以上説明したように、本発明によれば、画像の微細な部分である高周波成分を保持したまま、注目領域の画像を一定幅に圧縮、拡大することが可能であり、これにより、画像情報を有効に使うことができる効果がある。また、上記一定幅をフィルムなどの視認領域とすれば、フィルム上の視認領域の幅に注目領域を拡大縮小することができる効果がある。
【０１０７】
また、画像にす抜け領域がある場合に、す抜け領域を除去することにより、注目領域をさらに精度よく抽出することが可能であり、濃度値変換の精度が上がる効果がある。
【０１０８】
さらに、平滑化画像から最大値、最小値を抽出することにより、濃度値変換をより精度よく行うことができる効果がある。
【図面の簡単な説明】
【図１】第１の実施の形態に係る画像処理装置の構成を示すブロック図である。
【図２】第１の実施の形態での処理手順を示すフローチャートである。
【図３】入力画像と照射領域とす抜け領域とを除去した画像を示す構成図である。
【図４】第１の実施の形態での動作を説明するための特性図である。
【図５】第２の実施の形態に係る画像処理装置を適用したＸ線撮影装置の構成を示すブロック図である。
【図６】特徴量を抽出する領域を説明するための図である。
【図７】上記画像処理装置の第１の階調変換曲線の階調変換曲線を説明するための図である。
【図８】上記画像処理装置の第２の階調変換曲線の階調変換曲線を説明するための図である。
【図９】第２の実施の形態での処理手順を示すフローチャートである。
【符号の説明】
１０５データ収集回路
１０８ＣＰＵ
１０９メインメモリ
１１０操作パネル
１１２照射領域抽出回路
１１３濃度値変換回路
１１３ａす抜け削除回路
１１３ｂ最大最小値抽出回路
１１３ｃ平滑化画像作成回路
１１３ｄ高周波成分作成回路
１１３ｅ濃度値変換回路
１１２ｆ高周波成分足し込み回路

Claims

原画像データに対し諧調変換処理及び周波数処理を施す画像処理方法であって、
前記原画像データにおける被写体内の領域から、階調変換曲線を生成するための値を抽出する工程と、
前記原画像データに基づく高周波成分を生成する工程と、
前記高周波成分を、前記階調変換曲線の傾きに基づいて変更する工程と、
前記変更した高周波波成分に基づき画像データを生成する工程とを備えることを特徴とする画像処理方法。
前記階調変換曲線を生成するための値は、前記被写体内の領域から抽出された最大値、最小値及び特徴量の少なくともいずれか一つであることを特徴とする請求項１に記載の画像処理方法。
前記諧調変換曲線に基づいて前記原画像データを階調変換処理する工程を備えることを特徴とする請求項１に記載の画像処理方法。
前記被写体内の領域は、直接X線が照射されている領域であるす抜け領域及び該す抜け領域と一定幅で接する領域に基づき抽出されることを特徴とする請求項１乃至３のいずれか１項に記載の画像処理方法。
原画像データに対し諧調変換処理及び周波数処理を施す画像処理方法をコンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記憶媒体であって、
前記原画像データにおける被写体内の領域から、階調変換曲線を生成するための値を抽出する工程と、
前記原画像データに基づく高周波成分を生成する工程と、
前記高周波成分を、前記階調変換曲線の傾きに基づいて変更する工程と、
前記変更した高周波波成分に基づき画像データを生成する工程とをコンピュータに実行させるためのプログラムを記憶したコンピュータ読み取り可能な記憶媒体。