JPH0918811A - 性能が制限されたデジタルディスプレイ上への強調表示を行なうためのデジタル画像の補償的前処理方法 - Google Patents

性能が制限されたデジタルディスプレイ上への強調表示を行なうためのデジタル画像の補償的前処理方法

Info

Publication number
JPH0918811A
JPH0918811A JP8055191A JP5519196A JPH0918811A JP H0918811 A JPH0918811 A JP H0918811A JP 8055191 A JP8055191 A JP 8055191A JP 5519196 A JP5519196 A JP 5519196A JP H0918811 A JPH0918811 A JP H0918811A
Authority
JP
Japan
Prior art keywords
image
signal
frequency component
display
high frequency
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.)
Granted
Application number
JP8055191A
Other languages
English (en)
Other versions
JP3681213B2 (ja
Inventor
Edward Muka
エドワード・ムカ
Paul Ho
ポール・シー・ホ
Bruce Raymond Whiting
ブルース・レイモンド・ホワイティング
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Eastman Kodak Co
Original Assignee
Eastman Kodak Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Eastman Kodak Co filed Critical Eastman Kodak Co
Publication of JPH0918811A publication Critical patent/JPH0918811A/ja
Application granted granted Critical
Publication of JP3681213B2 publication Critical patent/JP3681213B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • G06T5/92Dynamic range modification of images or parts thereof based on global image properties
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/40Image enhancement or restoration using histogram techniques

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Radiography Using Non-Light Waves (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Picture Signal Circuits (AREA)
  • Closed-Circuit Television Systems (AREA)
  • Transforming Electric Information Into Light Information (AREA)

Abstract

(57)【要約】 【課題】 高レベルの周囲光が存在する環境下におい
て、観察者による画像理解が容易となるように、市販さ
れているディスプレイ装置上でデジタル画像を表示する
ための補償的な画像処理方法を提供すること。 【解決手段】 ディスプレイ装置上に画像を表示する方
法であって、デジタル画像信号を与える工程と、デジタ
ル画像信号において低い周波数成分の信号を識別する工
程と、デジタル画像信号において高い周波数成分の信号
を識別する工程と、低い周波数成分の信号から減衰され
た低い周波数成分の信号を生成する工程と、高い周波数
成分の信号から増幅された高い周波数成分の信号を生成
する工程と、減衰された低い周波数成分の信号と増幅さ
れた高い周波数成分の信号とを合成する工程とを有して
いる。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、強調的画像処理に
係り、特に、性能が制限されたディスプレイ装置上に表
示されるデジタル画像(放射線撮影画像)の強調処理を
行なうための補償的画像処理に関する。
【0002】
【従来の技術および発明が解決しようとする課題】診断
目的の医学的画像を表示する際には、通常、点灯された
ライトボックス上のフィルムが用いられていた。フィル
ムの画質は優れているが、フィルム自体がかさばるこ
と、保管するのが困難であること、移動が面倒であるこ
と、記憶装置に記録するのが困難であること、複製が容
易でないこと、および全体的な操作に対して費用が多く
かかることから、望ましくない点をも有している。これ
により、フィルムに代えて、デジタル画像およびデジタ
ルディスプレイを使用する傾向が強まっている。フィル
ムに対する主要な改善点として、デジタル画像が、電子
的に転送および記憶可能であることが上げられる。しか
し、現在のデジタルディスプレイは、表示できるデジタ
ル画像の画質に対して、性能的限界がある。(例えば、
CRTディスプレイのようなソフトコピーディスプレイ
は、現在、電子的放射線医学システムにおいて、最も一
般的に使用されている表示手段であり、本願明細書にお
いても、このソフトコピーディスプレイに関すること
が、説明の中心となっている。また、反射型のハードコ
ピー上へのデジタル方式によるプリントも、CRTに対
して記載される制限と同様の制限を受けるとともに、本
発明による手段により、同様にその画質が改善され
る。)
【0003】典型的なデジタルディスプレイに関して、
ライトボックス上で透照されるフィルムと比較すると、
その空間分解能は低く、最大輝度は小さく、そして、輝
度のダイナミックレンジが小さくなっている。これによ
り、デジタルディスプレイの画質は劣り、特に病院内で
高レベルの周囲光が照射される環境においては、画質の
劣化が顕著となる。そして、一般的には、ソフトコピー
ディスプレイに関する制限を解消する唯一の方法は、デ
ィスプレイの空間的分解能および全体的輝度を上げるこ
とであると考えられてきた。この手段はソフトコピーデ
ィスプレイをより有用なものとするが、装置のコストを
増大させるものとなる。それゆえ、高レベルの周囲光が
存在する場合でも、通常の安価なディスプレイを用いて
有用な診断が実施できるように、医学的デジタル画像に
対して補償的画像処理を実施することが望まれている。
【0004】さらに、ICU(集中管理ユニット)内の
CRTディスプレイ上に、胸部の放射線撮影画像を表示
することが強く望まれている。臨床医は、画像が放射線
施設内でプリントされICUに移送されるのをわざわざ
待つことなく、チューブの挿入位置、気胸(pneumothora
x)の欠落、および他の診断事項に関するフィードバック
情報を直ちに得ることができる。また、病院内に設置さ
れるディスプレイのコストは経済的制限を受けることか
ら、購買可能なディスプレイに関して、空間分解能、輝
度、およびダイナミックレンジ等の表示機能は制限を受
けることになる。さらに、通常、ICU内は、照明が明
るく灯されており、ディスプレイ上に高レベルの周囲光
が入射するので、ディスプレイの忠実な画像表示性能が
損なわれることになる。
【0005】本発明により解決される課題を明示するた
めに、フィルムおよびデジタルディスプレイの特性の詳
細な説明とともに、医学的画像の特性と診断者との関係
についても説明がなされる。
【0006】放射線撮影による胸部のフィルム画像は、
通常、感度測定によりlog10単位で約2.5から3.5の輝度
のダイナミックレンジを有する画像を与える。これによ
り、”DR(画像)”、すなわちフィルム画像の輝度の
ダイナミックレンジは、log1 0単位で2.5から3.5の値を
とる。観察用にライトボックス上で透照されたフィルム
画像は、暗い視覚室で測定されたところによると、通
常、最大輝度は2400ニト(nits)(カンテ゛ラ/m2)となり、最小
輝度は0.8ニトとなる。また、空間分解能に関しては、
通常、レーザプリンタは、80ミクロンの画素を有し、1
つのプリントラインには、4000以上の画素が配列されて
いる。
【0007】CRTディスプレイは、ライトボックス上
のフィルムと比較すると、使用可能な輝度のダイナミッ
クレンジが小さくなっている。通常のCRTディスプレ
イに関しては、暗い部屋で測定された場合に、その最大
輝度は160ニトとなり、最小輝度は0.5ニトとなる。通
常、”DR(ディスプレイ)”、すなわちCRTディス
プレイの輝度のダイナミックレンジは、約320、あるい
はlog10単位で2.5の値をとる。さらに、CRTディスプ
レイは、詳細部のコントラストを低減させるように作用
する、高い空間周波数における減衰作用により、空間分
解能が制限される。また、CRTディスプレイでは、増
幅器のスイッチングタイムおよびガラス面内のライトパ
イピング(light piping)に起因する制限により、最も高
い空間周波数において減衰が生じるようになる(Roehri
g, H; Ji, TL; Browne, M; Dall-as, WJ; Blume, Hartw
ig氏による"Signal-to-Noise Ratio and Maximum Infor
m-ation Content of Images Displayed by a CRT", Pro
ceedings SPIE MedicalImaging IV 1990. 1232: pp.115
-133を参照のこと)。現在のCRTディスプレイは、1
つの画像ラインに対して、1000から2000の画素が配列さ
れ、画素サイズは約250ミクロンとなっている。
【0008】さらに、CRTディスプレイでは、振幅信
号あるいは駆動レベルに対して、輝度が非線形的な応答
を示す(E. Muka, H. Blume, S. Daly氏による"Display
ofMedical Images on CRT Soft-Copy Displays", SPI
E Medical Imaging 1995, Vol.2431, Image Display, 1
995を参照のこと)。この輝度の非線形性は、画像表示
を制御するために、補正される必要がある。そして、こ
の非線形性の応答を補償するために、入力画像データ
が、表示される前に参照テーブル(LUT)を通され、
この結果、ディスプレイの非線形性が最小化され、ディ
スプレイ上に画像が輝度の関数として正確に表示され
る。この場合、それぞれのディスプレイが固有の特性を
有しているので、それぞれ所定の輝度応答が得られるよ
うに、固有のLUTが与えられている。この工程を経る
ことにより、透照されたフィルム画像に対して同じレベ
ルの画像をCRTディスプレイ上に表示することが可能
となる。
【0009】しかし、上記の輝度のダイナミックレンジ
の数値は、例えば、放射線撮影画像の読取り用の暗室の
ような周囲光がよく制御された観察用エリアに対して適
用されるものである。これに対して通常の病院内の作業
環境では、CRTディスプレイ上から照射される”グレ
ア輝度(glare luminance)”を生成する高レベルの周囲
光が存在している。例えば、ICU内においては、周囲
光により生成される”グレア輝度”により、通常、最小
輝度が約3.2ニトまで増加される。これにより、CRT
ディスプレイの輝度のダイナミックレンジが、約50あ
るいはlog10単位で1.7まで減少される。また、周囲光
は、場所によって相当に変化する。
【0010】上記の要因から、経験のある放射線技師の
見解によれば、CRTディスプレイ上に表示される画像
は、フィルムと比較すると感覚的にその表示能力が劣
り、デジタル画像はフィルム画像と異なるように見える
とのことである。このような画像表示に関する劣化要因
は、人間の視覚システム(HVS(human visual syste
m))に関する知見に基づいて判断される。HVSの動作
特性は、絶対輝度、空間周波数、画像の詳細部構造、ノ
イズ等のパラメータに基づいたコントラスト感度関数に
より特徴づけられる(A. B. Watson氏編集、S. Daly氏
著の"The Visibledifferences Predictor; Algolism fo
r the Assessment of Image Fidelity",Digital Image
and Human Vision, pp 179-206, Massachusetts: MIT P
ress, 1993を参照のこと)。CRTディスプレイ上で使
用可能な最大輝度は、ライトボックスで使用可能な最大
輝度と比較すると小さいので、観察者が、画像の詳細お
よびコントラストを検知することが困難となる。同様
に、デジタル画像において、フィルムと比較した場合
に、輝度のダイナミックレンジが減少することおよび空
間分解能が低減されることで、観察者の有用な情報を抽
出する能力が低減されることになる。
【0011】表示装置に課せられた機能は、診断用のす
べての画像情報を提供すること、すなわち観察者用に表
示された、現在使用されているフィルムから読取れるす
べての情報を提供することである。しかし、診断用に与
えられるこのような情報が、デジタルディスプレイ装置
の処理性能を越えるという問題が生じている(例えば、
フィルムの輝度のダイナミックレンジがlog10単位で3.2
であるのに対して、CRTディスプレイの輝度のダイナ
ミックレンジはlog10単位で1.7しかない)。この際、画
像処理技術を用いて、デジタル画像の振幅を減衰させ、
CRTディスプレイで利用可能である輝度の小さなダイ
ナミックレンジに画像情報を適合される手法が用いられ
ている。しかし、単純な減衰処理を用いると、画像コン
トラスト、およびエッジやテクスチャのような詳細画像
部の視認性に対して損失が生じることとなる。詳細な画
像特徴部は、原画像においても既に読取り困難であり、
これに輝度に関する要因で詳細な画像特徴部にさらなる
損失が加わると、CRTディスプレイ上の画像が診断に
適さなくなる場合がある。高レベルの周囲光が存在する
環境下のCRTディスプレイ上に、医学的画像を診断に
適応できるレベルで表示するとともに、フィルムと同じ
レベルの詳細画像部識別性能をCRTディスプレイで実
現するためには、より強力な画像処理方法を開発する必
要があることが明らかであろう。
【0012】1994年6月7日にNakazawa氏等に付与された
アメリカ合衆国特許5,319,719号には、オリジナル画像
の低い周波数成分の振幅をアンシャープマスキング(uns
harpmasking)を通して減少させることで、表示される画
像の輝度のダイナミックレンジを圧縮させる、放射線撮
影画像に対する画像処理装置が開示されている。しか
し、この特許には、デジタル放射線撮影画像の高い周波
数成分の増幅(boosti-ng)、あるいは、周囲光またはH
VSに関連する画像補償処理に関しての記載がなされて
いない。
【0013】1990年2月20日にHishinuma氏に付与された
アメリカ合衆国特許4,903,205号には、オリジナルのデ
ジタル放射線撮影画像をディスプレイ表示用に処理す
る、放射線撮影画像表示システムが開示されている。こ
の装置では、超低周波数フィルタに対応するアンシャー
プマスク(unsharp mask)信号が得られ、オリジナル信号
とアンシャープマスク信号との差異信号に強調係数を乗
じて得られた信号を、オリジナル信号に加えることで、
処理信号が生成される。しかし、この特許には、ディス
プレイ装置における高い周波数の画像データの減衰を補
償するための高い周波数領域の強調処理に関しては、開
示されていない。
【0014】1988年5月24日にHishinuma氏等に付与され
たアメリカ合衆国特許4,747,052号には、処理画像信号
Sが以下のオペレーションを用いて得られる、放射線撮
影画像に対する画像処理技術が開示されている。 S = SORG + b(SORG-SUS) ここで、Sは処理画像信号を示し、SORGはオリジナルの
デジタル信号を示し、bは周波数強調係数を示し、SUSは
アンシャープマスク信号を示す。この際、所定のサイズ
のアンシャープマスク信号内のデジタル信号のメジアン
(中央値)がとられる(1982年2月9日に発明者であるKa
to氏等に付与されたアメリカ合衆国特許4,315,318号、1
982年2月23日に発明者であるKato氏等に付与されたアメ
リカ合衆国特許4,317,179号、および1982年8月24日に発
明者であるIshida氏等に付与されたアメリカ合衆国特許
4,346,409号を参照のこと)。
【0015】上記のほとんどの従来技術は、結果を表示
するためのフィルム出力に関するものである。J. H. Ki
m氏等による論文"Improved Visualization of Simulate
dNodules by Adaptive Enhancement of Digital Chest
Radiography; Acad Radio, Vol.1, No.2, October 199
4, pp.93−99"には、CRTディスプレイ上に放射線撮
影による胸部画像を表示するための適応的強調アルゴリ
ズムが開示されている。この手法によれば、輝度のダイ
ナミックレンジが圧縮されるとともに、解剖学的に重要
である特定の部位の画像コントラストが選択的に強調さ
れる。しかし、この論文には、CRTディスプレイの輝
度の較正、あるいは、ディスプレイ装置が設置される視
覚的環境に対するHVSの関連あるいは動作に関して何
ら述べられていない。
【0016】上記のように、医学的画像用のソフトコピ
ーディスプレイに対する画像処理に関しては、種々の問
題点がある。
【0017】
【課題を解決するための手段】本発明により、従来のソ
フトコピーディスプレイに関する問題が解決される。す
なわち、本発明は、最大輝度が小さいこと、輝度のダイ
ナミックレンジが小さいこと、高い周波数領域における
応答性が低いこと等のデジタルディスプレイの特性、お
よび、高レベルの周囲光が存在する環境がディスプレイ
に与える悪影響によって、デジタルディスプレイに生じ
る画質の損失を前もって補償するための、医学的画像に
対する画像処理方法を提供する。この場合、人間の視覚
システムに関する知見に基づいて、観察者により画像処
理に対する適正なパラメータ値が選択されるとの認識が
前提となっている。
【0018】一般的に、本発明は、(例えば、集中管理
ユニット(ICU)内のような)高レベルの周囲光が存
在する環境下において、(例えば、臨床医のような)観
察者による画像理解が容易となるように、(例えば、C
RTディスプレイのような)市販されているディスプレ
イ装置上で、(例えば、放射線撮影画像のような)デジ
タル画像を表示するための補償的な画像処理方法を提供
する。本発明による方法は、単純であるとともに、コン
ピュータを用いた数値演算により実現可能であり、さら
に、視覚的な損失がないという利点を有している。視覚
的損失がないとは、すなわち、詳細部位の視覚的認識に
関して、透照された放射線撮影フィルムと同等の画像
が、CRTディスプレイ上に表示されることを意味す
る。高レベルの周囲光が存在する環境下のICUで実施
された本発明による画像処理方法に関する実験的観測に
よれば、気胸、挿入されたライン、および他の詳細部位
の視覚化に関して、上記の補償的画像処理方法を用いた
ソフトコピーディスプレイ(CRTディスプレイ)上
に、画質的にフィルムと同等の画像が表示される。例え
ば、露光過多のフィルムと比較された場合には、本発明
の画像処理方法を用いたソフトコピーディスプレイの表
示能力は、フィルムより優れていることもある。
【0019】本発明の1つの特徴によれば、医学用画像
処理コンピュータシステム内で、CRTディスプレイ上
の画像表示を強調するために実行される方法は、オリジ
ナルの医学用画像が有する輝度のダイナミックレンジお
よび空間周波数成分を有するデジタル画像データを与え
る工程と、デジタル画像の低い周波数成分の信号を識別
する工程と、デジタル画像の高い周波数成分の信号を識
別する工程と、デジタル画像において、低い周波数成分
の信号から、減衰された低い周波数成分の信号を生成す
る工程と、デジタル画像において、高い周波数成分の信
号から、増幅された高い周波数成分の信号を生成する工
程と、減衰された低い周波数成分の信号と増幅された高
い周波数成分の信号とを合成する工程とを有している。
これにより、ディスプレイ装置上に表示される画像に関
して、輝度のダイナミックレンジおよび空間周波数特性
を向上させるための前もって補償処理が実行された画像
信号が生成される。
【0020】
【発明の実施の形態】図1ないし図6は、本発明を説明
するために用いられる図である。図7は、本発明による
画像処理方法の1つの実施の形態を示す図である。図8
および図9は、定型化された放射線撮影画像を示す図で
ある。図10は、標準的なアルゴリズムによる伝達関数
と修正アンシャープマスクアルゴリズムによる伝達関数
との差異を示す図である。図11および図12は、本発
明を説明するために用いられる図である。図13は、本
発明による方法を組み込んだ装置を示すブロック図であ
る。
【0021】本発明による画像処理方法に関する以下の
記載では、デジタル放射線撮影画像に対する画像処理が
説明されているが、本発明が、他の種類の医学的画像、
およびディスプレイ上に表示可能である任意のデジタル
画像に対して適応可能であることが解されるであろう。
【0022】図13には、本発明が組込まれた装置が開
示されている。しかし、当業者には、本発明の範囲内に
おいて、図13に示された装置以外の装置に対しても本
発明が組込み可能であることが解されるであろう。図に
示されるように、画像ソース(画像源)20によりデジ
タル画像が生成され、この画像データが画像データファ
イル22に記憶される。画像ソース20としては、例え
ば、(X線源、診断用画像源(CT,MRI,US,PET)のよう
な)医学的画像あるいは放射線撮影画像源、蛍光物質リ
ーダ、X線フィルムデジタイザ、(磁気ディスクあるい
は光ディスク等の)画像記録ソース等が上げられる。
(陰極線管(CRT),液晶ディスプレイ,プラズマディス
プレイ,反射型ハードコピーを出力するプリンタ等の)
コンピュータディスプレイ24に表示される前に処理さ
れる画像データを表現する画像データファイルは、コン
ピュータデータメモリ26内に記憶される。また、本発
明による画像処理方法(アルゴリズム)は、(ROM,
ハードディスク等の)メモリ28内に記憶可能であり、
この画像処理プログラムのコピーがコンピュータ命令メ
モリ30内に記憶される。そして、(ユーザ入力装置、
他のコンピュータ等の)画像処理パラメータ入力装置3
3を用いて、画像処理パラメータがメモリ30内に入力
される。
【0023】次に、(例えば、マイクロプロセッサのよ
うな)コンピュータの中央処理装置(CPU)32によ
り、メモリ30に記憶された画像処理アルゴリズムに基
づいて、メモリ26に記憶された画像データが処理され
る。そして、処理された画像データは、コンピュータデ
ィスプレイメモリ34内に記憶され、その後、観察者3
6が見れるようにディスプレイ24上に表示される。
【0024】本発明に関する以下の説明においては、胸
部の放射線撮影画像に対する画像処理が対象となる。す
なわち、フィルム画像と比較した場合の画像の全体的な
変化を最小限に抑えるとともに疑似画像を最小化して、
手ごろな価格の市販のディスプレイ(1K(1024)×1K(1
024)の画素マトリクス、160ニトの最大輝度値)上に医
学上有用な画像が表示されるように、胸部の放射線撮影
画像の画像処理を実行することを意図する。
【0025】初期研究で使用される画像は、中間速度に
より、広いラティテュード(適正露光範囲)のスクリー
ンフィルムに撮影された複数の気胸のフィルムから構成
されている(研究が進むにつれて、他の部位が撮影され
た画像も含まれることになる)。複数の画像は、35cm×
47cmのフィルム上に撮影されたPA(後側−前側)放射
線撮影画像および側面放射線撮影画像から構成されてい
る。そして、20枚のPAフィルムのなかから放射線撮
影画像として特徴的な画像が選択された。これらの画像
のなかには、セントラルカテーテル、ペースメーカ、お
よび胸に挿入されるチューブと同様に、微妙な気胸、張
りのある明らかな気胸、間隙性の肺の病気等が撮影され
た画像が含まれていた。
【0026】[画像処理方法]放射線撮影画像のダイナ
ミックレンジの大部分の領域は、低い周波数成分で占め
られている。しかし、HVSは、これらの低い周波数成
分に対して、比較的感度が小さくなっている。画像にお
ける低い周波数成分は、画像種類、画像内における大ま
かな位置、および(例えば、PA撮影による胸部 VS.
側面撮影による胸部 VS. 腕部、(縦)隔膜 VS. 肺領
域、ソフト組織 VS. 骨、横隔膜下の肺 VS. 横隔膜上の
肺のような)組織種類に関する情報を伝達する。これに
より、有用な画像情報に損失を与えることなしに、低い
周波数成分のダイナミックレンジを減少させることが可
能である。
【0027】また、高い周波数の画像データにより、診
断において最も重要である詳細部のエッジ情報およびテ
クスチャ情報が伝達される。HVSは、これらの高い周
波数成分に対して大きな感度を有している。低い周波数
成分を減衰させる一方で、高い周波数の画像データを現
状維持あるいは強調することで、より有用な画像をデジ
タルディスプレイ上に表示することが可能となる。この
効果は、図1ないし図6に示される画像特性を解析する
ことで明らかにされる。
【0028】図1には、放射線撮影による未処理の胸部
画像の中央部の640,000画素(800画素×800画素の領
域)に対するヒストグラムが示されている。この画像
は、オリジナル画像のすべての周波数成分を有してい
る。曲線の幅は、最も明るい(最も濃度の低い)画素
と、最も暗い(最も濃度の高い)画素との間の差異を表
している。そして、この曲線の幅は、輝度のダイナミッ
クレンジの大きさに比例している。 (濃度)=log10(輝度)
【0029】図2には、一辺が59ピクセルの正方形の
大きさを有するアンシャープマスクを用いて、画像の画
素値(濃度値)が処理あるいはぼかされた後の、図1に
示された放射線撮影による胸部画像の中央部の640,000
画素に対するヒストグラムが示されている。アンシャー
プマスク画像は、低い周波数成分(0.13本/mm以下)の
みを含んでいる。ここで、図2に示される低い周波数成
分のみから構成された画素データの輝度のダイナミック
レンジを表す曲線の幅が、輝度のダイナミックレンジを
表す図1の曲線の幅とほぼ同じであることに留意する必
要がある。これは、画像の低い周波数成分が、全画像の
ダイナミックレンジの大部分を占めることを示してい
る。
【0030】図3には、図1に示された放射線撮影画像
データから、図2に示された放射線撮影画像データを引
いた際の中央部の640,000画素のヒストグラムが示され
ている。この画像は、高い周波数成分(0.13本/mm以
上)のみを有して構成されている。ここで、曲線の幅
が、図1あるいは図2の曲線の幅よりも、かなり小さく
なっていることに留意する必要がある。これは、画像の
高い周波数成分は、非常に小さなダイナミックレンジを
占めるにすぎないことを示している。すなわち、高い周
波数成分を表示するには、オリジナルのダイナミックレ
ンジの比較的小さな領域が要求されるにすぎない。
【0031】図4は、同じ画像に対する高速フーリエ変
換(FFT(Fast Fourier Transfo-rm))を表してい
る。図1に示される放射線撮影による胸部画像の中央部
の512画素×512画素の領域が、垂直方向に平均化され、
512画素の平均スキャンラインが生成される。そして、
この平均スキャンラインにフーリエ変換が適用され、そ
の振幅の大きさが表示されている。ここで、低い周波数
成分の振幅が高い周波数成分の振幅よりも明らかに大き
なことに留意する必要がある。
【0032】図5は、図2に示された放射線撮影による
胸部画像の512画素×512画素の中央領域に対するFFT
が示されている。ここで、低域フィルタ画像において、
高い周波数成分が低減されていることに留意する必要が
ある。
【0033】図6は、図3に示された放射線撮影による
胸部画像の512画素×512画素の中央領域に対するFFT
が示されている。ここで、高域フィルタ画像において、
低い周波数成分が低減されていることに留意する必要が
ある。
【0034】この種の画像処理の1つの問題点は、高い
周波数成分を増幅することで、グラニュラ雑音も強めら
れることである。通常、このグラニュラ雑音は、診断に
用いられる詳細情報の周波数よりも高い周波数で現れ
る。本発明では、ディスプレイにおいて、複雑な処理を
用いることなく、高い周波数成分の減衰作用により、こ
のグラニュラ雑音が減衰される効果を利用している。し
かし、高い周波数成分における減衰作用を有しないディ
スプレイ装置が用いられる場合には、フィルタリングお
よび最終的な低域フィルタに対して、異なるパラメータ
を設定する必要がある。
【0035】[好適な画像処理方法]一般的に、図7に
示される本発明の画像処理方法の好適な実施の形態によ
れば、画像データが、低い空間周波数成分と高い空間周
波数成分とに分離される。低い周波数成分は減衰され、
高い周波数成分は増幅される。低い周波数成分を減衰さ
せることで、画像の詳細部分を実質的に損失させること
なしに、全画像のダイナミックレンジを低減させること
ができる。観察者は、低い周波数成分のコントラストの
変化を認識するが、それでも画像種類および画像内の特
定領域位置を識別することが可能となっている。また、
前もって高い周波数成分を増幅することで、CRTディ
スプレイの特性および周囲光により引き起こされる高い
周波数成分の減衰作用を補償するための画像の前処理的
強調が実行される。さらに、高い周波数成分が増幅され
ることで、画像内に存在するエッジの強調が支援され
る。このような作用が図8および図9に例示されてい
る。
【0036】図8には、放射線撮影画像の1つの領域を
表す定型化されたグラフが示されている。この波形は、
大きな振幅および低い周波数を有する方形波に、小さな
振幅および高い周波数を有する正弦波が重畳されて形成
されている。
【0037】図9には、図8の画像データに対して、本
発明による画像処理を実行した後のフィルタリングされ
た画像データが示されている。ここで、高い周波数の正
弦波の振幅が増加され、低い周波数の方形波の振幅が低
減されていることに留意する必要がある。この際、全体
的な振幅が低減されるとともに、ダイナミックレンジが
小さくなり、これにより、制限されたダイナミックレン
ジを有するディスプレイ装置に対して、画像データが適
合されるようになる。また、同時に、ディスプレイ装置
におけるコントラストの減衰を補償するための、高い周
波数成分の強調が実施される。
【0038】画像を周波数成分に分離するための数値演
算により有効となる1つの方法として、アンシャープマ
スク(unsharp mask)(USM)フィルタアルゴリズムが
上げられる("Pratt, William K. Digital Image Proce
ssing 2nd Ed. John Wiley &Sons,Inc. New York. 199
1. pp.303-305"を参照のこと)。このUSM手法では、
画素の正方形のマスク(カーネル)の均一に重み付けさ
れた空間的平均を用いることで、オリジナルの画像に対
してぼかされた画像、あるいは低域(通過)画像が生成
される。そして、オリジナルの画像データからぼかされ
た(低域)画像データが引かれ、これにより、画像の高
い周波数成分あるいは高域通過成分のみが残される。2
つの周波数成分は、適切にサイズが変更された後に加算
されて、ダイナミックレンジが減少され、鮮鋭さが強調
された画像が生成される。数値演算は、以下のように簡
略化することが可能である。
【0039】 (高域画像) = (オリジナル画像)−(低域画像) (出力画像) = (低域係数)*(低域画像) +(高域係数)*(高域画像) = (低域係数)*(低域画像) +(高域係数)*{(オリジナル画像)−(低域画像)} = {(低域係数)−(高域係数)}*(低域画像) +(高域係数)*(オリジナル画像) この場合、画像を最適化するためのパラメータとして
は、低域係数、高域係数、および周波数成分に分解する
ためのマスクサイズが上げられる。
【0040】本発明による強調手法では、通常のアンシ
ャープマスクと異なり、画像の低域成分および高域成分
に対して、それぞれ独立にサイズ変更が実施される。通
常のアンシャープマスク手法では、画像の高域成分のみ
が増幅され、画像の低域成分は修正されないでそのまま
保持される。これにより、通常のアンシャープマスク手
法では、ダイナミックレンジを減少させることなく、エ
ッジの強調が実行される。
【0041】図10には、高い周波数成分に対する高域
係数2.0を有し、1辺が13画素の正方形として与えられ
る通常のアンシャープマスクフィルタと、低域係数0.7
および高域係数1.2を有し、1辺が59画素の正方形とし
て与えられる修正アンシャープマスクフィルタとに対す
る伝達関数(0.345mm画素)が示されている。ここで、
通常のアンシャープマスク手法と比較すると、修正アン
シャープマスク手法により、低域成分が減衰され、より
低い周波数から強調が開始され、高い周波数におけるゲ
インが低減されていることに留意する必要がある。
【0042】図11には、修正アンシャープマスクアル
ゴリズムが適用された結果が示されている。この図に
は、図1に示された放射線撮影画像を図10に示された
修正アルゴリズムで処理した後の画像の中央部の640,00
0画素のヒストグラムが示されている。ここで、画像の
ダイナミックレンジが減少されていることに留意する必
要がある。
【0043】図12には、図1に示された放射線撮影画
像に修正アンシャープマスクアルゴリズムを適用した結
果が示されており、中央部の512画素×512画素の領域の
画像データがFFTを用いて表現されている。図におい
て、実線で示されるオリジナル画像と、連続した点から
なる点線により示される修正アンシャープマスクアルゴ
リズムにより処理された画像との対比が示されている。
また、破線により、オリジナル画像を低域係数(0.7)で
サイズ変更した際の画像が示されている。
【0044】[パラメータ選択] −低域係数− 本発明による画像処理方法は、画像の輝度のダイナミッ
クレンジ、すなわちDR(画像)に基づいて実施され
る。この画像処理方法は、画像が表示されるCRTディ
スプレイが、DR(画像)より小さな輝度のダイナミッ
クレンジを有しているという事実に基づいて作成されて
いる。これにより、ディスプレイ上の処理された画像の
ダイナミックレンジ、すなわちDR(ディスプレイ)
は、オリジナル画像のダイナミックレンジ、すなわちD
R(画像)よりも小さくなる。そして、画像処理におい
て、当初、{DR(ディスプレイ)/DR(画像)}が
低域係数として設定される。
【0045】好適な実施の形態においては、低域係数値
は、それぞれの適用対象においてダイナミックレンジに
要求されるサイズに応じて設定される。既に述べたよう
に、通常のICU内においては、CRTディスプレイの
利用可能な輝度のダイナミックレンジは、log10単位で
1.7である。また、ICU内においては、通常、高レベ
ルの周囲光が入射するから、透照されたオリジナル画像
の目視可能な輝度のダイナミックレンジは、通常log10
単位の感度測定値で、3.5から約2.5に減少する(R.Boll
en氏およびJ.Vranckx氏による論文"Influence of ambie
nt light on thevisual sensitometric properties of,
and detail perception on, a radiogr-aph", Proc. S
PIE Vol.273 Application of Optical Instrumentation
in Medi-cine IX, pp.57-61, 1981を参照のこと)。こ
れら2つのダイナミックレンジの比(1.7/2.5)は、0.7と
なる。この比は、低域係数値に対して上限値を与えるも
のである。したがって、初期の実験においては、低域係
数値が0.7に設定された。
【0046】−高域係数− しかし、既に述べたように、低域係数を用いて画像コン
トラストを減衰させるのは、望ましいことではない。そ
れゆえ、低域係数により生じた減衰が、画像の高い空間
周波数成分において回復される。図3に示されるよう
に、高い周波数成分に対しては、画像の輝度のダイナミ
ックレンジのわずかの領域が用いられているにすぎな
い。それゆえ、画像処理においては、当初、 (高域係数)=1/(低域係数) として、高域係数値が設定される。
【0047】さらに、CRTディスプレイでは、空間周
波数が増加すると、コントラスト伝達における損失によ
って、画像コントラストに損失が生じることがある。こ
の特性により、ディスプレイ上に表示される画像のさら
なる補償(強調)が要求される。そして、この見地よ
り、高域係数が増加される。
【0048】最終的な高域係数値は、観察者により望ま
れるエッジおよび鮮鋭さの強調の度合により決定され
る。この場合、ユーザが、現在の仕事に対して必要とさ
れる鮮鋭さの度合を選択することもある。詳細な画像を
必要とする仕事に対しては、より詳細な画像、すなわ
ち、より大きな高域係数で処理された画像がディスプレ
イ上に表示されることが望まれる。
【0049】−マスクサイズ− 強調される空間周波数は、以下に示されるようなHVS
モデルおよび通常のソフトコピーディスプレイ特性に基
づいて決定される。HVSは、4(本/度)の空間周波数
において、コントラスト感度関数(CSF(contrast se
nsitivity fun-ction))の最大値を与える。なお、この
空間周波数も、輝度の絶対値の関数となっている。そし
て、医学的画像の表示に対する通常の観察状況において
は、観察者は、ディスプレイ装置から約50cm離間して位
置している。これは、ディスプレイ上では、0.5(本/m
m)の空間周波数においてCSFの最大値が与えられる
ことを示している。少なくとも、画像は、この空間周波
数において強調されるので、HVSが効率的に機能する
ことが可能となる。画像は、この空間周波数の近傍の有
用な範囲において、強調されるのが望ましい。
【0050】医学的画像における有用な画像特徴部は、
その長さが2cmあるいはそれ以下であることが多く、こ
れは画像特徴部の空間周波数が約0.05(本/mm)以上で
あることを示している。これにより、画像強調は、0.05
(本/mm)の空間周波数で開始され、 0.5(本/mm)の
空間周波数を越えて実施されるのが好適である。USM
は、このような高域フィルタとして機能する。2cmのU
SMマスクを選択することで、(低域係数を用いた処理
工程における)輝度のダイナミックレンジの圧縮により
画像に生じた減衰を回復させることが可能となり、これ
により、HVSが、CRT上の画像を効率的に見ること
が可能となる。図10を構成するのに用いられた計算手
法に基づけば、USMにより、{1.4/(マスクサイス゛)}の値
を越えるすべての空間周波数に対して、実質的な強調工
程が実施される。これにより、2cmのマスクサイズが与
えられた場合には、上記の条件を満たす0.07(本/mm)
以上の空間周波数に対して、実質的な強調工程が実施さ
れる。
【0051】通常のCRTディスプレイは、約25cmの画
像高さを有し、この方向の1ラインに約1000個の画素が
配列されるから、画素密度は4(画素/mm)となる。通常
のX線による胸部撮影画像は、CRTディスプレイの可
視領域より大きく、この画像が高さ43cmの大きさを有し
ていることから、画像は縮小され、CRTディスプレイ
に表示される際には小さく見えることになる。それゆ
え、X線撮影画像で2cmの画像特徴部は、CRTディス
プレイ上では、1.2cmに縮小される。これにより、本発
明による修正USM手法のためのマスクを構成するため
には、約49個(マスクサイズは常に奇数である)のC
RT上の画素を必要とする。これに前節に記載された条
件を適用すると、修正USM手法は、0.11(本/mm)の
空間周波数において、最大のゲインを有する。ディスプ
レイのサイズが、実際の画像より幾分小さいので、医学
的画像用のCRTディスプレイには、一般的に、拡大可
能である画像特徴部が含まれている。そして、拡大機能
を設けることで、オペレータがCRTディスプレイ上の
画像特徴部のサイズを増加させることが可能となり、こ
れにより、フィルム上で表示されるのと同じ大きさで、
画像特徴部がCRTディスプレイ上に表示される。この
選択的な拡大工程により、修正USM手法により強調が
生じる空間周波数値が減少される。
【0052】通常、USMのような強調画像処理は、高
い空間周波数において認識可能なノイズを増加させる。
空間周波数のオクターブに基づいたマルチ分解能表示手
法(multiresolution representation)のような複雑な画
像処理を用いて、高い周波数における強調を制限すると
ともに、認識可能な高い周波数成分におけるノイズを低
減することが行われている(S. Ranganath氏による論
文"Image FilteringUsing Multiresolution Represent
ation", IEEE Transactions of Pattern An-alysis an
d Machine Intelligence, Vol.13, No.5, pp. 426-440,
May, 1991を参照のこと)。CRTディスプレイでは、
空間周波数が増加するのに応じて、画像コントラストに
損失が生じるので、CRTディスプレイが、表示される
画像の低域フィルタとして効率的に機能する。この効果
の結果として、本発明により、非常に単純な方法で、修
正USM手法を有効に用いることが可能となる。
【0053】しかし、USM手法は一般的な方法である
ので、マルチ分解能表示アプローチを用いて、修正US
Mアルゴリズムの強調効果を実行することが可能であ
る。また、LUT(参照テーブル)を用いた高域成分の
変更および低域成分の変更に関して、修正USM手法と
同一の機能を実現するために、従来型のUSMフィルタ
を用いることが可能である。
【0054】−ディスプレイ表示に対応した較正− 本発明によれば、輝度が較正されたソフトコピーディス
プレイが与えられ、これにより、画像データが、輝度に
関して正確に表示される。多くの画像処理手法は、この
ような較正工程を有するものではない。HVSが非線形
的な特性を有しているので、認識可能な詳細画像部を生
成するうえで、較正処理は重要な意味を有する。さら
に、フィルムおよびライトボックスを用いた従来の手法
により表示された画像を認識するように訓練された放射
線技師が慣れた画像(例えば、透照されたフィルム画
像)に近い画像が望まれるので、この見地においても、
較正処理は重要な意味を有する。強調処理が完了された
画像空間を表すフィルム濃度値を入力値として受け入れ
るために、LUTを用いてCRTディスプレイに対する
画像データが較正される。
【0055】−サンプル画像に対するテスト結果− 初期研究では、通常のICUの照明条件の下で、通常の
CRTディスプレイ上において、様々なレベルの高域係
数およびマスクサイズに対して、テストが実行された。
このテストにおいては、画像処理に関する既に記載され
たすべての要因に対して、試験サンプルが与えられてい
る。
【0056】研究において使用された胸部画像は、すべ
てlog10単位で2.5のDR(画像)を有していた。また、
周囲光が存在する環境下で動作するディスプレイ装置の
輝度のダイナミックレンジ、DR(ディスプレイ)は、
log10単位で1.7に制限されていた。これにより、低域係
数は、0.7に設定された。そして、高域係数は、当初、1
/0.7で1.4に設定された。この実験条件に対して、肺実
質(lung parenchyma)の表示の変化および高い周波数に
おけるノイズの生成を最小限にすることと、気胸のライ
ンを強調することとのバランスをとるために、高域係数
値として1.2が選択された。
【0057】サンプル画像に対して、高域係数およびマ
スクサイズの最適値は、画像のオリジナルサイズ、およ
びフィルムからデジタル画像データを生成するフィルム
デジタイザの分解能にも影響を受ける。標準サイズの胸
部撮影フィルム(35cm×43cm)は、4k(4096)×5k(5120)
の画素(117画素/cm)にデジタル化される。そして、4k
×4kの正方形領域が抽出され、1k×1kの領域(29画素/c
m、0.345mm/画素)に対する画素の平均化により、ダウ
ンサイジングが実施される。
【0058】また、3から99までの範囲のマスクサイ
ズに対して、テストが実施された。マスクより小さな画
像対象部は、全体が強調され、マスクより大きな画像対
象部は、エッジが強調されるとともに減衰された。小さ
なマスクが使用された場合には、大きな周波数レンジに
わたって画像が減衰され、これにより、輝度のダイナミ
ックレンジがより圧縮される。しかし、2重ラインの疑
似画像を防止するために、マスクは、例えば血管のよう
な画像内に表示される最も標準的な構造体のサイズより
も大きい必要がある。
【0059】最適な画像を生成するマスクサイズは、幅
が約59画素あるいは2.0cmのものであった。このマスク
サイズは、Prokop氏等の研究に関連を有するものとなっ
ている(Prokop氏、Mathias氏、Shaefer氏、Cornelia M
氏、Oestmann氏、Jorg W氏、Galanski氏、Michael氏に
よる"Improved Parameters for Unsharp Mask Filt-eri
ng of Digital Chest Radiographs", Radiology 1993.
97:521-526を参照のこと)。これには、2.5cmのマスク
サイズが、1.4mm、5mm、および7cmのマスクサイズより
も優れているのが判明したと開示されている。
【0060】小さなマスクでは、大きなマスクと比較し
て、自然な画像を得ることができなかった。小さなマス
クにより、”2重ライン”の疑似画像を有するエッジ強
調画像が生成された。この画像は、観察者に対して深さ
を感じさせない平坦な画像として見られた。これは、か
なり高い周波数に到達するまで、強調が行われないこと
に起因する。小さなマスクを使用するすべての画像にお
いて肺部が際だって表示されるので、通常の肺と間隙性
の肺の病気とを識別することが困難となっていた。
【0061】大きなマスクで生成された画像は、階調ス
ケールが良好に維持されている。エッジは少し強調さ
れ、2重ラインの疑似画像は最小化された。1つの画像
では、横隔膜下の小結節が存在していた。大きなマスク
を用いて処理された画像において、小結節は良好に維持
されていた。
【0062】図10に示されるように、修正USMフィ
ルタ(59画素)では、約0.01(本/mm)まで低域係数(0.
7)を用いて減衰が実施され、その後コントラストを回復
(5%のUSMゲイン)させるための高域増幅が実施さ
れる。実質的な回復および強調(最大ゲイン)は、0.07
(本/mm)の空間周波数で開始される。また、図10に
は、2.0のゲインを有する13画素の標準的なUSMフィ
ルタが示されている。通常、市販の放射線撮影画像の数
値処理システムでは、このような小さなマスクが用いら
れている。
【0063】当業者に周知の他の画像処理方法を、アン
シャープマスクによる画像処理方法に代えて使用可能で
あることが解されるであろう。
【0064】−画像解析− 処理画像の低い周波数成分におけるコントラストは、オ
リジナル画像よりも弱くなっている。階調スケールは変
更されるが、画像は容易に認識可能であり、CRTディ
スプレイ上の表示は観察者にとって抵抗のあるものとは
ならない。すなわち、CRTディスプレイ上の画像が、
フィルム画像と類似するものとなる。肺領域のようなオ
リジナル画像で非常に暗い部分は、明るく再生されると
ともに高い周波数のコントラストが強調され、これによ
り、肺実質において詳細画像部がより明確に区別可能と
なる。また、(縦)隔膜のようなオリジナル画像で非常
に明るい部分では、高い周波数のコントラストが強調さ
れ、脊椎骨を容易に識別することが可能となる。また、
気胸ラインのような肺領域内の細かなラインを強調する
弱いエッジ強調が実行されている。肺部は強調される
が、この場合、間隙性の肺の病気を有する画像と、通常
の肺の画像とを区別することが可能となっている。
【0065】−本発明による放射線撮影画像の観察者に
よるテスト− ミズーリ州セントルイスのMIR(Mallinckrodt Instit
ute of Radiology)から出席した3人の胸部放射線技師
が、本発明により処理された画像を観察した。これらの
観察者は、ソフトコピーディスプレイとともに、放射線
画像数値処理システム(CR(computed radiography))
からプリントされたアンシャープマスク画像を観察した
経験を有していた。
【0066】初期研究において使用されたCRTディス
プレイは、イムロギックス・モデル1000(Imlogix model
1000 electronic viewbox)であった。画面サイズは、2
5cm×25cmであり、分解能は、1024画素×1024画素×12
ビットであった。また、ディスプレイは、暗室において
log10単位で2.5のダイナミックレンジを与えるように、
160ニトの最大輝度および0.5ニトの最小輝度を有するよ
うに較正された。ディスプレイは、LUT(参照テーブ
ル)を用いて、log10単位の輝度に対して、個別に線形
化された。
【0067】ICU室内を模擬した照明環境下におい
て、2つのCRTディスプレイが並べて配置された。上
方から光線が照射され、3.2ニトの強度のフレアライト
(flarelight)が生成されるので、CRTディスプレイ上
で利用可能なダイナミックレンジが、log10単位で1.7ま
で減少される。本発明により処理された画像と、処理さ
れない同じ画像とが、並んで配置された2つのCRTデ
ィスプレイ上に表示される。ディスプレイに対する観察
者の好みを削除するために、処理された画像の表示ディ
スプレイが交換された。そして、画像のオリジナルフィ
ルムが、CRTディスプレイ上のライトボックスに設置
された。観察者は、CRTディスプレイから約1m、ラ
イトボックスから約1.5m離れて位置していた。ライト
ボックスは、表示サイズの差異を補償するために、観察
者からより離れて設置された。また、観察者は、最も見
やすい位置に移動することが許されていた。
【0068】本発明による画像処理技術により、胸部の
放射線撮影画像を、高レベルの周囲光が存在する環境下
で動作する1k×1kのCRTディスプレイ上に、医学的に
有効にデジタル表示することが可能となる。気胸および
挿入されたチューブのような詳細画像部の可視化に関し
て、ソフトコピーディスプレイが、フィルムに対して画
質的に同等であることが判明された。また、露光過多の
フィルムに対しては、ソフトコピーディスプレイは、オ
リジナルフィルムよりも画質的に優れていた。さらに、
CRTディスプレイ上においては、すべてのサンプル画
像に対して、本発明による方法で処理された画像は、処
理されない画像よりも、画質的に同等かあるいは優れて
いた。
【0069】1人の放射線技師の意見によれば、本発明
を用いることで、1k×1kのディスプレイ上において、気
胸の可視化に関して満足のいく表示が得られたとのこと
であった。従来では、このような満足のいく表示を得る
ためには、高い分解能のディスプレイが必要であると考
えられていた。
【0070】ここで用いられる一般的な方法は、画像種
類あるいはディスプレイ性能に依存するものではない。
必要であるのは、可視化のために画像データがCRTデ
ィスプレイあるいは類似の装置に表示されることであ
る。画像サイズ、画像種類、ディスプレイ装置の特性、
および観察領域内での周囲光の状態に応じて、パラメー
タ値のみが選択される必要がある。
【0071】高レベルの周囲光が存在する環境下で動作
するソフトコピーディスプレイの制限されたダイナミッ
クレンジを、高画質の反射型プリントのダイナミックレ
ンジに近づけることで、担当医に対して、診断リポート
に付加される医学的画像のペーパコピー(paper copy)を
プリント出力することが可能となる。このようなペーパ
コピーは、写真コピーと比較して、コストおよびスペー
スの面で効率的である。
【0072】図を参照して、本発明が詳細に説明された
が、本発明の範囲内において、開示された方法と同等の
機能を有する変更あるいは修正された方法の実行が可能
であることが解されるであろう。
【図面の簡単な説明】
【図1】未処理の放射線撮影画像における濃度値のヒス
トグラムを示す図である。
【図2】図1の画像をアンシャープマスクを用いて処理
した後の画像における濃度値のヒストグラムを示す図で
ある。
【図3】図1の画像から図2の画像を引いた画像におけ
る濃度値のヒストグラムを示す図である。
【図4】図1に示された画像の高速フーリエ変換を示す
図である。
【図5】図2に示された画像の高速フーリエ変換を示す
図である。
【図6】図3に示された画像の高速フーリエ変換を示す
図である。
【図7】本発明による画像処理方法の1つの実施の形態
を示す図である。
【図8】放射線撮影画像の1つの領域を表す定型化され
た画像データを示す図である。
【図9】図8の画像データに対して、本発明による画像
処理を実行した後のフィルタリングされた画像データを
示す図である。
【図10】標準的なアルゴリズムによる伝達関数と修正
アンシャープマスクアルゴリズムによる伝達関数との差
異を示す図である。
【図11】図1に示された画像に修正アンシャープマス
クアルゴリズムを適用した後の画像における濃度値のヒ
ストグラムを示す図である。
【図12】修正アンシャープマスクアルゴリズムの効果
を示すために、画像データに高速フーリエ変換を適用し
た結果を示す図である。
【図13】本発明による方法を組み込んだ装置を示すブ
ロック図である。
【符号の説明】
20 画像ソース(画像データ) 24 CRTディスプレイ 28 画像処理アルゴリズム格納メモリ
───────────────────────────────────────────────────── フロントページの続き (72)発明者 ポール・シー・ホ アメリカ合衆国・ウィスコンシン・ 53717・マディソン・ポンウッド・サーク ル・66−0 (72)発明者 ブルース・レイモンド・ホワイティング アメリカ合衆国・テキサス・75024・プラ ノ・オールド・ポンド・ドライヴ・4560

Claims (3)

    【特許請求の範囲】
  1. 【請求項1】 周知の特性を有するとともに制御可能に
    動作するディスプレイ装置上に画像を表示する方法であ
    って、 前記ディスプレイ装置上に表示された際に、劣化される
    輝度のダイナミックレンジおよび空間周波数特性を有す
    るデジタル画像信号を与える工程と、 前記デジタル画像信号において、低い周波数成分の信号
    を識別する工程と、 前記デジタル画像信号において、高い周波数成分の信号
    を識別する工程と、 前記低い周波数成分の信号から、前記デジタル画像信号
    に関して減衰された低い周波数成分の信号を生成する工
    程と、 前記高い周波数成分の信号から、前記デジタル画像信号
    に関して増幅された高い周波数成分の信号を生成する工
    程と、 前記減衰された低い周波数成分の信号と前記増幅された
    高い周波数成分の信号とを合成する工程とを有し、 前記ディスプレイ装置上に画像が表示された際に、輝度
    のダイナミックレンジおよび空間周波数特性が改善され
    ている前処理された画像信号を生成することを特徴とす
    るディスプレイ装置上に画像を表示する方法。
  2. 【請求項2】 請求項1記載のディスプレイ装置上に画
    像を表示する方法において、 前記減衰された低い周波数成分の信号と前記増幅された
    高い周波数成分の信号とを合成する前記工程が、 例えば輝度のダイナミックレンジのような表示される
    (オリジナル)画像データに関するパラメータ、例えば
    輝度のダイナミックレンジのような前記ディスプレイ装
    置の特性、前記ディスプレイ装置の動作する環境におけ
    る輝度のダイナミックレンジ、前記ディスプレイ装置の
    空間周波数応答特性、および観察者の視覚システムに関
    する見識のなかの1つあるいは複数の要素に対して関数
    的な処理を行なうことで実施されることを特徴とするデ
    ィスプレイ装置上に画像を表示する方法。
  3. 【請求項3】 周知の特性を有するとともに制御可能に
    動作するディスプレイ装置上に画像を表示する方法であ
    って、 前記ディスプレイ装置上に表示された際に、劣化される
    輝度のダイナミックレンジおよび空間周波数特性を有す
    るデジタル画像信号を与える工程と、 前記デジタル画像信号において、高い周波数成分の信号
    を識別する工程と、 前記高い周波数成分の信号から、前記デジタル画像信号
    に関して振幅が変えられた(増幅された)高い周波数成
    分の信号を生成する工程と、 この振幅が変えられた高い周波数成分の信号と、前記デ
    ジタル画像信号とを合成する工程と、 この合成デジタル画像信号を比例的に減少されること
    で、前記合成デジタル画像信号の減衰信号を生成する工
    程とを有し、 前記ディスプレイ装置上に画像が表示された際に、輝度
    のダイナミックレンジおよび空間周波数特性が改善され
    ている前処理された画像信号を生成することを特徴とす
    るディスプレイ装置上に画像を表示する方法。
JP05519196A 1995-03-14 1996-03-12 性能が制限されたデジタルディスプレイ上への強調表示を行なうためのデジタル画像の補償的前処理方法 Expired - Fee Related JP3681213B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US08/404,399 US5774599A (en) 1995-03-14 1995-03-14 Method for precompensation of digital images for enhanced presentation on digital displays with limited capabilities
US404,399 1995-03-14

Publications (2)

Publication Number Publication Date
JPH0918811A true JPH0918811A (ja) 1997-01-17
JP3681213B2 JP3681213B2 (ja) 2005-08-10

Family

ID=23599445

Family Applications (1)

Application Number Title Priority Date Filing Date
JP05519196A Expired - Fee Related JP3681213B2 (ja) 1995-03-14 1996-03-12 性能が制限されたデジタルディスプレイ上への強調表示を行なうためのデジタル画像の補償的前処理方法

Country Status (4)

Country Link
US (1) US5774599A (ja)
EP (1) EP0732669B1 (ja)
JP (1) JP3681213B2 (ja)
DE (1) DE69625839T2 (ja)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001238870A (ja) * 1999-12-28 2001-09-04 Ge Medical Syst Sa 器官の厚みを補償する方法および装置
JP2001314396A (ja) * 1999-12-30 2001-11-13 General Electric Co <Ge> デジタルx線イメージング検出器のダイナミックレンジ拡張システム
JP2002335447A (ja) * 2001-05-09 2002-11-22 Shimadzu Corp ディジタルサブトラクション装置
WO2006129529A1 (ja) * 2005-06-02 2006-12-07 Konica Minolta Holdings, Inc. 画像処理方法及び画像処理装置
US7958921B2 (en) 2009-03-16 2011-06-14 The Yokohama Rubber Co., Ltd. Pneumatic tire with tread having sub grooves in blocks, center main groove and subsidiary main grooves
JP2020010786A (ja) * 2018-07-17 2020-01-23 コニカミノルタ株式会社 放射線画像表示制御装置、放射線画像解析装置及び放射線撮影システム

Families Citing this family (87)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3706185B2 (ja) * 1995-12-25 2005-10-12 富士写真フイルム株式会社 画像再生方法および装置
JP4035668B2 (ja) * 1996-06-10 2008-01-23 富士ゼロックス株式会社 画像処理装置、画像処理方法、および記憶媒体
EP0818932B1 (en) * 1996-07-09 1999-11-10 STMicroelectronics S.r.l. Circuit for enhancing chrominance transitions in real-time video reception
IL119283A0 (en) * 1996-09-19 1996-12-05 Elscint Ltd Adaptive filtering
US6243095B1 (en) * 1996-12-05 2001-06-05 Peter E. Shile Navigation and display system for digital radiographs
JP3700804B2 (ja) * 1996-12-13 2005-09-28 富士写真フイルム株式会社 画像処理方法および装置
AUPO525897A0 (en) * 1997-02-24 1997-03-20 Redflex Traffic Systems Pty Ltd Digital image processing
US6069979A (en) * 1997-02-25 2000-05-30 Eastman Kodak Company Method for compressing the dynamic range of digital projection radiographic images
US5978518A (en) * 1997-02-25 1999-11-02 Eastman Kodak Company Image enhancement in digital image processing
US6211919B1 (en) * 1997-03-28 2001-04-03 Tektronix, Inc. Transparent embedment of data in a video signal
US6285798B1 (en) * 1998-07-06 2001-09-04 Eastman Kodak Company Automatic tone adjustment by contrast gain-control on edges
US6424730B1 (en) 1998-11-03 2002-07-23 Eastman Kodak Company Medical image enhancement method for hardcopy prints
US6392764B1 (en) * 1998-11-12 2002-05-21 Xerox Corporation Systems and methods for adaptively filtering palettized images
US6421468B1 (en) * 1999-01-06 2002-07-16 Seiko Epson Corporation Method and apparatus for sharpening an image by scaling elements of a frequency-domain representation
US6304277B1 (en) 1999-01-15 2001-10-16 Colorcentric.Com, Inc. Remote modification of digital images using scripts
DE19936615A1 (de) * 1999-08-04 2000-03-02 Werner Stein Verfahren zur Verbesserung der Kontrastwiedergabe bei digitalisierten Bildern
FR2799028B1 (fr) * 1999-09-27 2002-05-03 Ge Medical Syst Sa Procede de reconstitution d'une image tridimentionnelle d'elements de fort contraste
US6690488B1 (en) * 1999-09-30 2004-02-10 Polaroid Corporation Method and apparatus for estimating the spatial frequency response of a digital image acquisition system from the images it produces
EP1257972B1 (en) * 2000-02-02 2006-12-06 Quvis, Inc. System and method for optimizing image resolution using pixelated imaging devices
US6724943B2 (en) * 2000-02-07 2004-04-20 Sony Corporation Device and method for image processing
US6633684B1 (en) 2000-07-07 2003-10-14 Athentech Technologies Corp. Distortion-free image contrast enhancement
US6775407B1 (en) * 2000-08-02 2004-08-10 Eastman Kodak Company Producing a final modified digital image using a source digital image and a difference digital image
WO2002025589A1 (en) * 2000-09-20 2002-03-28 Nik Multimedia, Inc. Digital image sharpening system
US6989859B2 (en) * 2000-12-22 2006-01-24 Eastman Kodak Company Camera having user interface ambient sensor viewer adaptation compensation and method
EP1239415A1 (en) * 2001-03-05 2002-09-11 Ge Medical Systems Sa Method and system of management of the dynamics of a digitized radiological image
US20020154323A1 (en) * 2001-03-10 2002-10-24 Sobol Robert E. Method for variable contrast mapping of digital images
US6807299B2 (en) * 2001-03-10 2004-10-19 Hewlett-Packard Development Company, L.P. Method for contrast mapping of digital images that converges on a solution
FR2823338A1 (fr) * 2001-04-10 2002-10-11 Koninkl Philips Electronics Nv Procede et dispositif de post traitement d'images numeriques
US6804384B2 (en) * 2001-06-12 2004-10-12 Mclean Hospital Corporation Color magnetic resonance imaging
US7003173B2 (en) * 2001-06-12 2006-02-21 Sharp Laboratories Of America, Inc. Filter for combined de-ringing and edge sharpening
US7619585B2 (en) 2001-11-09 2009-11-17 Puredepth Limited Depth fused display
US7064740B2 (en) 2001-11-09 2006-06-20 Sharp Laboratories Of America, Inc. Backlit display with improved dynamic range
US6993200B2 (en) * 2001-11-20 2006-01-31 Sony Corporation System and method for effectively rendering high dynamic range images
CA2412706C (en) * 2001-11-23 2009-12-29 Imaging Dynamics Company Ltd. Correcting geometric distortion in a digitally captured image
NZ517713A (en) 2002-06-25 2005-03-24 Puredepth Ltd Enhanced viewing experience of a display through localised dynamic control of background lighting level
NZ521505A (en) 2002-09-20 2005-05-27 Deep Video Imaging Ltd Multi-view display
US7130485B2 (en) * 2002-10-02 2006-10-31 Eastman Kodak Company Enhancing the tonal and color characteristics of digital images using expansive and compressive tone scale functions
US7116838B2 (en) * 2002-10-25 2006-10-03 Eastman Kodak Company Enhancing the tonal and spatial characteristics of digital images using selective spatial filters
US7444034B1 (en) * 2002-11-06 2008-10-28 Digivision, Inc. Systems and methods for image enhancement in multiple dimensions
US7668390B2 (en) * 2002-11-06 2010-02-23 Digivision, Inc. Systems and methods for image enhancement in multiple dimensions
US7106352B2 (en) * 2003-03-03 2006-09-12 Sun Microsystems, Inc. Automatic gain control, brightness compression, and super-intensity samples
US7742651B2 (en) * 2003-04-04 2010-06-22 Thomson Licensing Method and apparatus for enhancing detail in an image
DE10317367B4 (de) * 2003-04-15 2007-01-11 Siemens Ag Verfahren zur Durchführung der digitalen Subtraktionsangiographie unter Verwendung nativer Volumendatensätze
NZ525956A (en) 2003-05-16 2005-10-28 Deep Video Imaging Ltd Display control system for use with multi-layer displays
US7315656B2 (en) * 2003-05-22 2008-01-01 The Boeing Company Methods and apparatus for enhanced viewing of aerial refueling operations
EP1484707A3 (en) * 2003-06-03 2006-05-17 Konica Minolta Medical & Graphic, Inc. Medical image system, and medical image processing method
US20040263942A1 (en) * 2003-06-24 2004-12-30 Lopez Matthew Grant Display system allowing enhanced dynamic range
CN101686306A (zh) * 2003-09-11 2010-03-31 松下电器产业株式会社 可视处理装置、基于它的装置以及可视处理方法
KR101027825B1 (ko) * 2003-09-11 2011-04-07 파나소닉 주식회사 시각 처리 장치, 시각 처리 방법, 시각 처리 프로그램 및반도체 장치
EP1667065B1 (en) * 2003-09-11 2018-06-06 Panasonic Intellectual Property Corporation of America Visual processing apparatus, visual processing method, visual processing program, and semiconductor device
US7271780B2 (en) * 2003-09-23 2007-09-18 Eastman Kodak Company Display device and system
US7616829B1 (en) * 2003-10-29 2009-11-10 Apple Inc. Reducing undesirable block based image processing artifacts by DC image filtering
US8107760B2 (en) * 2003-11-06 2012-01-31 Z Microsystems Visualization Technologies, Llc Systems and methods for image enhancement in multiple dimensions
WO2005052673A2 (en) * 2003-11-21 2005-06-09 Sharp Laboratories Of America, Inc. Liquid crystal display with adaptive color
US7777714B2 (en) 2004-05-04 2010-08-17 Sharp Laboratories Of America, Inc. Liquid crystal display with adaptive width
US8395577B2 (en) 2004-05-04 2013-03-12 Sharp Laboratories Of America, Inc. Liquid crystal display with illumination control
US7602369B2 (en) * 2004-05-04 2009-10-13 Sharp Laboratories Of America, Inc. Liquid crystal display with colored backlight
US7872631B2 (en) 2004-05-04 2011-01-18 Sharp Laboratories Of America, Inc. Liquid crystal display with temporal black point
US20060074307A1 (en) * 2004-05-10 2006-04-06 Tatsuo Igarashi Body cavity diagnostic system
EA011562B1 (ru) * 2004-07-27 2009-04-28 Дюрр Денталь Гмбх & Ко.Кг Способ и устройство для улучшения распознаваемости различных структур на изображениях, полученных в проходящих лучах
US7426314B2 (en) * 2004-07-30 2008-09-16 Hewlett-Packard Development Company, L.P. Adjusting pixels by desired gains and factors
US7898519B2 (en) 2005-02-17 2011-03-01 Sharp Laboratories Of America, Inc. Method for overdriving a backlit display
EP1650730B1 (en) 2004-10-25 2009-12-30 Barco NV Optical correction for high uniformity panel lights
US8050511B2 (en) 2004-11-16 2011-11-01 Sharp Laboratories Of America, Inc. High dynamic range images from low dynamic range images
US8050512B2 (en) * 2004-11-16 2011-11-01 Sharp Laboratories Of America, Inc. High dynamic range images from low dynamic range images
DE102004061507B4 (de) * 2004-12-21 2007-04-12 Siemens Ag Verfahren zur Korrektur von Inhomogenitäten in einem Bild sowie bildgebende Vorrichtung dazu
KR100640063B1 (ko) * 2005-02-18 2006-10-31 삼성전자주식회사 외부조도를 고려한 영상향상방법 및 장치
JP4538358B2 (ja) * 2005-03-30 2010-09-08 パイオニア株式会社 画像処理装置
US8014034B2 (en) * 2005-04-13 2011-09-06 Acd Systems International Inc. Image contrast enhancement
US7639849B2 (en) 2005-05-17 2009-12-29 Barco N.V. Methods, apparatus, and devices for noise reduction
NZ542843A (en) 2005-10-05 2008-08-29 Pure Depth Ltd Method of manipulating visibility of images on a volumetric display
US9143657B2 (en) 2006-01-24 2015-09-22 Sharp Laboratories Of America, Inc. Color enhancement technique using skin color detection
US8121401B2 (en) 2006-01-24 2012-02-21 Sharp Labortories of America, Inc. Method for reducing enhancement of artifacts and noise in image color enhancement
JP4980625B2 (ja) * 2006-02-21 2012-07-18 富士フイルム株式会社 体腔内観察装置
CN101106563A (zh) * 2006-07-14 2008-01-16 Ge医疗系统环球技术有限公司 服务系统
US8941580B2 (en) 2006-11-30 2015-01-27 Sharp Laboratories Of America, Inc. Liquid crystal display with area adaptive backlight
JP5009377B2 (ja) * 2007-01-15 2012-08-22 ヒューマンアイズ テクノロジーズ リミテッド レンチキュラ印刷のための方法およびシステム
US8520060B2 (en) * 2007-02-25 2013-08-27 Humaneyes Technologies Ltd. Method and a system for calibrating and/or visualizing a multi image display and for reducing ghosting artifacts
US8432411B2 (en) 2007-05-18 2013-04-30 Pure Depth Limited Method and system for improving display quality of a multi-component display
WO2009013744A2 (en) * 2007-07-23 2009-01-29 Humaneyes Technologies Ltd. Multi view displays and methods for producing the same
US8330768B2 (en) 2007-07-27 2012-12-11 Sharp Laboratories Of America, Inc. Apparatus and method for rendering high dynamic range images for standard dynamic range display
US8346010B2 (en) * 2008-09-17 2013-01-01 Panasonic Corporation Image processing device, imaging device, evaluation device, image processing method, and optical system evaluation method
US8249394B2 (en) * 2009-03-12 2012-08-21 Shmuel Peleg Method and system for shift-map image editing
US9524700B2 (en) 2009-05-14 2016-12-20 Pure Depth Limited Method and system for displaying images of various formats on a single display
US8928682B2 (en) 2009-07-07 2015-01-06 Pure Depth Limited Method and system of processing images for improved display
JP5505007B2 (ja) * 2010-03-18 2014-05-28 富士通株式会社 画像処理装置、画像処理方法及び画像処理用コンピュータプログラム
JP6520920B2 (ja) * 2014-03-27 2019-05-29 ノーリツプレシジョン株式会社 画像処理装置

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4315318A (en) * 1978-12-26 1982-02-09 Fuji Photo Film Co., Ltd. Method and apparatus for processing a radiation image
US4317179A (en) * 1978-12-26 1982-02-23 Fuji Photo Film Co., Ltd. Method and apparatus for processing a radiographic image
JPS5691735A (en) * 1979-12-25 1981-07-24 Fuji Photo Film Co Ltd Method and apparatus for treating xxray image
US4674125A (en) * 1983-06-27 1987-06-16 Rca Corporation Real-time hierarchal pyramid signal processing apparatus
US4747052A (en) * 1984-11-14 1988-05-24 Philips Medical Systems, Inc. Radiation image processing
US4697594A (en) * 1985-08-21 1987-10-06 North American Philips Corporation Displaying a single parameter image
US4903205A (en) * 1985-10-15 1990-02-20 Fuji Photo Film Co. Ltd. Method and apparatus for displaying radiation image, and method and apparatus for calculating unsharp mask signal used for the same
DE3629396C2 (de) * 1986-08-29 1993-12-23 Agfa Gevaert Ag Verfahren zur elektronischen Bildverarbeitung
US5012333A (en) * 1989-01-05 1991-04-30 Eastman Kodak Company Interactive dynamic range adjustment system for printing digital images
US5319719A (en) * 1991-05-15 1994-06-07 Konica Corporation Processing apparatus for radiographic image signals
US5387983A (en) * 1991-09-27 1995-02-07 Minolta Camera Kabushiki Kaisha Facsimile apparatus comprising converting means for converting binary image data into multi-value image data and image processing apparatus judging pseudo half-tone image

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001238870A (ja) * 1999-12-28 2001-09-04 Ge Medical Syst Sa 器官の厚みを補償する方法および装置
JP2001314396A (ja) * 1999-12-30 2001-11-13 General Electric Co <Ge> デジタルx線イメージング検出器のダイナミックレンジ拡張システム
JP2002335447A (ja) * 2001-05-09 2002-11-22 Shimadzu Corp ディジタルサブトラクション装置
JP4626086B2 (ja) * 2001-05-09 2011-02-02 株式会社島津製作所 ディジタルサブトラクション装置
WO2006129529A1 (ja) * 2005-06-02 2006-12-07 Konica Minolta Holdings, Inc. 画像処理方法及び画像処理装置
JP4811401B2 (ja) * 2005-06-02 2011-11-09 コニカミノルタホールディングス株式会社 画像処理方法及び画像処理装置
US8131101B2 (en) 2005-06-02 2012-03-06 Konica Minolta Holdings, Inc. Image processing method and image processing apparatus
US7958921B2 (en) 2009-03-16 2011-06-14 The Yokohama Rubber Co., Ltd. Pneumatic tire with tread having sub grooves in blocks, center main groove and subsidiary main grooves
JP2020010786A (ja) * 2018-07-17 2020-01-23 コニカミノルタ株式会社 放射線画像表示制御装置、放射線画像解析装置及び放射線撮影システム

Also Published As

Publication number Publication date
EP0732669B1 (en) 2003-01-22
JP3681213B2 (ja) 2005-08-10
DE69625839T2 (de) 2003-11-06
DE69625839D1 (de) 2003-02-27
US5774599A (en) 1998-06-30
EP0732669A1 (en) 1996-09-18

Similar Documents

Publication Publication Date Title
JP3681213B2 (ja) 性能が制限されたデジタルディスプレイ上への強調表示を行なうためのデジタル画像の補償的前処理方法
US6424730B1 (en) Medical image enhancement method for hardcopy prints
EP1501049B1 (en) Method for rendering digital radiographic images for display based on independent control of fundamental image quality parameters
US5717791A (en) Image contrast enhancing method
US5805721A (en) Method and apparatus for contrast enhancement
US6069979A (en) Method for compressing the dynamic range of digital projection radiographic images
EP0527525B1 (en) Method and apparatus for contrast enhancement
US6813335B2 (en) Image processing apparatus, image processing system, image processing method, program, and storage medium
JP2000050081A (ja) エッジのコントラスト利得制御による自動ト―ン調節方法
JP2012524329A (ja) 多重尺度画像の正規化及び強調
US7155044B2 (en) Multiscale gradation processing method
JP4486784B2 (ja) ノイズ低減方法
Van Metter et al. Enhanced latitude for digital projection radiography
US7626592B1 (en) Multiscale gradation processing method
US8111887B2 (en) Image processing method, image processing apparatus and image processing program
EP1347413A1 (en) Method for enhancing the contrast of an image.
US20020057827A1 (en) Method, apparatus, and recording medium for processing tomographic image
EP1347414A1 (en) Method for enhancing the contrast of an image.
EP1341125A2 (en) Multiscale gradation processing method
JP2004206426A (ja) 画像処理装置、画像処理方法、画像出力方法、プログラム、プログラムを記録した記録媒体及び画像出力システム
Sherrier et al. Image optimization in a computed-radiography/photostimulable-phosphor system
US20030063784A1 (en) Method, apparatus, program and recording media for displaying medical image
Kim et al. Improved visualization of simulated nodules by adaptive enhancement of digital chest radiography
EP4187481A1 (en) Greyscale gradation processing method
US20040252909A1 (en) Image processing method, image processing apparatus, program, and computer storage medium

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20040413

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20040712

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20040715

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20041012

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20041109

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20050202

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20050311

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: 20050419

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050517

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