JP3696339B2 - 画像処理方法および装置 - Google Patents
画像処理方法および装置 Download PDFInfo
- Publication number
- JP3696339B2 JP3696339B2 JP21733496A JP21733496A JP3696339B2 JP 3696339 B2 JP3696339 B2 JP 3696339B2 JP 21733496 A JP21733496 A JP 21733496A JP 21733496 A JP21733496 A JP 21733496A JP 3696339 B2 JP3696339 B2 JP 3696339B2
- Authority
- JP
- Japan
- Prior art keywords
- image signal
- original image
- signals
- frequency response
- signal
- 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
- 238000003672 processing method Methods 0.000 title claims description 11
- 230000006870 function Effects 0.000 claims description 74
- 238000006243 chemical reaction Methods 0.000 claims description 66
- 230000004044 response Effects 0.000 claims description 60
- 238000000034 method Methods 0.000 claims description 36
- 230000002708 enhancing effect Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 description 23
- 238000001914 filtration Methods 0.000 description 18
- 238000010586 diagram Methods 0.000 description 17
- 238000004364 calculation method Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 101000920618 Homo sapiens Transcription and mRNA export factor ENY2 Proteins 0.000 description 1
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 102100031954 Transcription and mRNA export factor ENY2 Human genes 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
- G06T5/75—Unsharp masking
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Transforming Light Signals Into Electric Signals (AREA)
- Image Analysis (AREA)
Description
【発明の属する技術分野】
本発明は、画像信号の所定の周波数成分を強調するための画像処理方法および装置に関し、特に詳しくは所望の画像処理結果を得るための画像処理パラメータの自動設定に関する。
【0002】
【従来の技術】
従来、本出願人により、非鮮鋭マスク画像信号(以下、ボケ画像信号という)を用いて周波数強調処理を行って放射線画像の診断性能を向上させる数々の画像処理方法および装置が提案されている(特開昭55-163472 号、同55-87953号など)。ここで、ボケ画像信号とは、画素数は原画像信号と同じであるが原画像信号よりも鮮鋭度が低い画像を表す画像信号であり、原画像信号の所定の周波数以上の高周波成分が除去された周波数応答特性を有する信号である。
【0003】
上記周波数強調処理は、原画像信号Sorg からボケ画像信号Susを引いたものに強調係数βを乗じたものを、原画像信号Sorgに加算することにより、原画像信号の所定の空間周波数成分を強調するものである。これを式で表すと下記の式(1)のようになる。
【0004】
Sproc=Sorg+β×(Sorg−Sus) …(1)
(Sproc:周波数強調処理された信号、Sorg :原画像信号、Sus:ボケ画像信号、β:強調係数)
この際、上記処理では信号を加算することによりアーチファクトが発生してしまうことがあるが、これは、原画像信号Sorgに加算する加算信号の周波数応答特性を調整することにより解決でき、そのための具体的な調整方法として、次のような方法が提案されている。(特願平8-182155号など)。
【0005】
この方法とは、まず鮮鋭度の異なる、すなわち周波数応答特性の異なる複数のボケ画像信号を作成し、そのボケ画像信号および原画像信号の中の2つの信号の差分をとることにより、原画像信号の、ある限られた周波数帯域の周波数成分を表す複数の帯域制限画像信号(以下、バンドパス信号という)を作成し、さらにそのバンドパス信号をそれぞれ異なる変換関数によって所望の大きさとなるように抑制してから、その複数の抑制されたバンドパス信号を積算することにより上記加算信号を作成するものである。この処理を式として表すと、例えば、下記の式(2)
(但し、Sproc:高周波成分が強調された画像信号
Sorg :原画像信号
Susk(k=1〜N):ボケ画像信号
fk(k=1〜N):各バンドパス信号を変換する変換関数
β(Sorg):原画像信号に基づいて定められる強調係数)
のようになる。
【0006】
この方法によれば、バンドパス信号を変換する変換関数f1〜fNの定義を変えることによって原画像信号Sorgに加算する加算信号の周波数応答特性を調整することができる。したがって、アーチファクトを防止するためには、変換関数f1〜fNを、アーチファクトが発生しないような周波数応答特性を有する画像信号が作成されるような関数として定義すればよい。さらには、アーチファクトの防止という目的に限らず、各変換関数の定義次第で、所望の周波数応答特性を有する処理済画像信号を得ることができる。
【0007】
【発明が解決しようとする課題】
しかしながら、上記方法に基づいて、実際に変換関数を定義しようとした場合、各変換関数をどのように定義すれば所望の結果が得られるかは、上記(2)のような式から簡単にわかるというものではない。そのため、従来は、強調処理された画像を観察しながら、変換関数を定義するパラメータ(例えば傾きなど)を少しずつ調整して、それぞれの変換関数を定義していた。これは、経験の少ないユーザにとっては非常に困難な作業であり、また熟練したユーザにとっても、時間のかかる作業であった。
【0008】
本発明は上記問題に鑑みて、経験の少ないユーザでも簡単に、上記変換関数の定義パラメータを調整して、所望の周波数応答特性を有する処理済画像信号を得ることができる画像処理方法および装置を提供することを目的とする。
【0009】
【課題を解決するための手段】
本発明の画像処理方法は、原画像を表す原画像信号に基づいて互いに周波数応答特性が異なる複数のボケ画像信号を作成し、原画像信号およびその複数のボケ画像信号に基づいて、あるいは複数のボケ画像信号に基づいて原画像信号の複数の周波数帯域ごとの信号を表す複数のバンドパス信号を作成し、その各バンドパス信号を複数の変換関数に基づいて変換して複数の変換画像信号を作成し、その各変換画像信号を積算して得た積算信号を原画像信号に加算することにより原画像信号の所定の周波数成分が強調された処理済画像信号を得る画像処理方法であって、前記処理済画像信号の所望の周波数応答特性を指定し、指定された周波数応答特性に基づいて前記変換関数の定義パラメータを決定してその変換関数を定義することを特徴とするものである。
【0010】
また、本発明の画像処理装置は、上記画像処理方法に基づいて画像処理を行う装置であって、原画像を表す原画像信号に基づいて互いに周波数応答特性が異なる複数のボケ画像信号を作成するボケ画像信号作成手段と、原画像信号およびその複数のボケ画像信号に基づいて、あるいは複数のボケ画像信号に基づいて原画像信号の複数の周波数帯域ごとの信号を表す複数のバンドパス信号を作成し、その各バンドパス信号を複数の変換関数に基づいて変換して複数の変換画像信号を作成し、その各変換画像信号を積算して得た積算信号を原画像信号に加算することにより原画像信号の所定の周波数成分が強調された処理済画像信号を得る処理手段とからなる画像処理装置であり、前記処理済画像信号の所望の周波数応答特性を指定する周波数応答特性入力手段と、指定された周波数応答特性に基づいて前記変換関数の定義パラメータを決定してその変換関数を定義する変換関数定義手段とを有することを特徴とするものである。
【0011】
ここで、「ボケ画像信号」とは、画素数は原画像信号と同じであるが原画像信号よりも鮮鋭度が低い画像を表す画像信号である。ボケ画像信号は、まず原画像信号の画素に対して所定間隔ごとに所定のフィルタリング処理を施すことによって画素を間引きし、そのようにして得た画像信号に対して同様のフィルタリング処理を繰り返してさらに画素数を少なくした画像信号を複数作成し、そのそれぞれに対して、所定の補間方法により原画像と画素数が同じになるように補間処理を施すことにより作成する。ここで、フィルタリングや補間については、一般に広く使用されている種々の方法を適用することができる。
【0012】
また、「原画像信号の複数の周波数帯域ごとの信号を表す複数のバンドパス信号」は、例えば隣接する周波数帯域のボケ画像信号同士で差分をとって作成してもよいし、原画像信号と各ボケ画像信号の差分をとって作成してもよい。あるいは原画像信号とボケ画像信号の他の組み合わせで差分をとって作成することもできる。なお、変換画像信号の作成、積算信号の作成、積算信号の原画像信号への加算は、例えば上記(2)式で表すことができる。
【0013】
また、「所定の周波数成分を強調する」とは、例えば画像のエッジ部を強調するために、高周波成分を強調するといったことを意味する。
【0014】
また、「処理済画像信号の所望の周波数応答特性を指定する」方法としては、例えば、操作画面上でユーザに対していくつかの周波数を示し、そのそれぞれに対して所望の周波数応答特性を数値として入力させてもよいし、あるいは周波数強調された画像信号の周波数応答特性のグラフを画面上に表示し、ユーザがマウス等のポインティングデバイスを用いてそのグラフを所望の形状となるように変形するようにしてもよい。
【0015】
また、「指定された周波数応答特性に基づいて前記変換関数の定義パラメータを決定してその変換関数を定義する」とは、例えば、前記バンドパス信号、入力された所望の周波数応答特性、および変換関数の定義パラメータの関係式を、バンドパス信号と周波数応答特性を既知の値とし、変換関数の定義パラメータを変数とする連立方程式として解くことによりパラメータを求めるなどの方法を意味する。
【0016】
【発明の効果】
本発明の画像処理方法および装置は、指定された周波数応答特性から、画像処理で用いられる変換関数の定義パラメータを逆算して求め、そのパラメータにより定義される変換関数を用いて画像処理を行うため、所望の周波数応答特性を指定するだけで、指定したとおりの周波数応答特性の画像信号を得る画像処理を行うことができ、パラメータ決定のための試行錯誤が不要となるばかりでなく、経験の少ないユーザでも簡単に所望の画像信号が得られるような画像処理を行うことができる。
【0017】
【発明の実施の形態】
以下、本発明の画像処理方法および装置の一実施の形態について、図面を参照して詳細に説明する。以下に示す画像処理装置は、蓄積性蛍光体シートに記録された人体の放射線画像を読み取って得た画像信号に対して、その画像が診断に適した画像となるように、ボケ画像信号を使用して周波数強調処理を施すものであり、処理された画像信号は主としてフィルムに記録され、診断に用いられる。
【0018】
図1はこの画像処理装置の概略を表す図である。画像処理装置1は、ボケ画像信号を作成するボケ画像信号作成手段2と、特定の周波数を強調するための周波数強調処理を行う非線形処理手段3を有する。さらに、この画像処理装置1は、変換関数定義手段4と、周波数応答特性入力手段5とを備える。変換関数定義手段4は、非線形処理手段3が変換処理に使用する変換関数を定義する手段であり、例えば関数の傾きなどのパラメータを決定して変換関数を定義する。周波数応答特性入力手段5は、強調処理された画像信号の所望の周波数応答特性をユーザが指定できるようにするための手段であり、変換関数定義手段4は周波数応答特性入力手段5によりユーザが指定した周波数応答特性に基づいて、パラメータを決定する。
【0019】
ここで、まずボケ画像信号の作成処理について詳細に説明する。図2はボケ画像信号作成処理の概要を示すブロック図である。図2に示されるように、図1のボケ画像信号作成手段2は、まずフィルタリング処理手段10により、原画像信号Sorgに対し、原画像の画素のx方向およびy方向に対してフィルタリング処理を施して原画像信号よりも解像度が低い画像信号B1(以下、低解像度画像信号という)を作成し、次にこの低解像度画像信号B1に対して同様のフィルタリング処理を施してこの低解像度画像信号B1よりもさらに解像度が低い低解像度画像信号B2を作成し、以降順次同様のフィルタリング処理を重ねていくものである。そして、補間処理手段11により、このフィルタリング処理の各段において得られる低解像度画像信号Bkに対して、それぞれ補間拡大処理を施して、鮮鋭度の異なる複数のボケ画像信号Sus1〜SusNを得るものである。
【0020】
本実施の形態においては、上記フィルタリング処理のフィルタとして、一次元ガウス分布に略対応したフィルタを使用する。すなわちフィルタのフィルタ係数を、ガウス信号に関する下記の式(3)
【0021】
【数1】
【0022】
にしたがって定める。これは、ガウス信号は周波数空間および実空間の双方において、局在性がよいためであり、例えば上記(3)式においてσ=1とした場合の5×1の1次元フィルタは図3に示すようなものとなる。
【0023】
フィルタリング処理は、図4に示すように、原画像信号Sorg に対して、あるいは低解像度画像信号に対して1画素おきに行う。このような1画素おきのフィルタリング処理をx方向、y方向に行うことにより、低解像度画像信号B1の画素数は原画像の1/4となり、フィルタリング処理により得られる低解像度画像信号に対して繰り返しこのフィルタリング処理を施すことにより、得られるn個の低解像度画像信号Bk(k=1〜n)は、それぞれ、画素数が原画像信号の1/22kの画像信号となる。
【0024】
次に、このようにして得られた低解像度画像信号Bkに対して施される補間拡大処理について説明する。補間演算の方法としては、Bスプラインによる方法など種々の方法が挙げられるが、本実施の形態においては、上記フィルタリング処理においてガウス信号に基づくローパスフィルタを用いているため、補間演算についてもガウス信号を用いるものとする。具体的には、下記の式(4)
【0025】
【数2】
【0026】
において、σ=2k-1 と近似したものを用いる。
【0027】
画像信号B1を補間する際には、k=1であるためσ=1となる。この場合、補間処理を行うためのフィルタは、図5に示すように5×1の一次元フィルタとなる。この補間処理は、まず低解像度画像信号B1に対して1画素おきに値が0の画素を1つずつ補間することにより低解像度画像信号B1を原画像と同一のサイズに拡大し、次に、この補間された低解像度画像信号B1に対して上述した図5に示す一次元フィルタによりフィルタリング処理を施すことにより行われる。
【0028】
同様に、この補間拡大処理を全ての低解像度画像信号Bkに対して行う。低解像度画像信号Bkを補間する際には、上記式(4)に基づいて、3×2k−1の長さのフィルタを作成し、画像信号Bkの各画素の間に値が0の画素を2k−1個ずつ補間することにより、原画像と同一サイズに拡大し、この値が0の画素が補間された画像信号Bkに対して3×2k−1の長さのフィルタにより、フィルタリング処理を施すことにより補間拡大する。
【0029】
次に、上記のようにして作成されたボケ画像信号を用いて行われる非線形処理について説明する。以下に示す処理は、特定の周波数成分を強調するための処理であり、原画像信号およびボケ画像信号は、上記(2)式にしたがって処理されるものとする。
【0030】
図6はこのような処理を行う周波数強調処理装置の一例を示す図である。この周波数強調処理装置では、減算器21として示されるように、原画像信号Sorgとボケ画像信号Suskのうち隣接する2つの信号に基づいてバンドパス信号が作成される。図7は、このようにして作成されたバンドパス信号の周波数応答特性の一例を示す図であり、6個のバンドパス信号の特性が示されている。各バンドパス信号は、変換器22においてそれぞれ異なる変換関数f1〜fNにより所望の大きさとなるように抑制され、さらに上記(2)式にしたがって、その複数の抑制されたバンドパス信号が演算器23において積算され、さらに原画像信号に加算されて、処理済画像信号Sprocが生成される。
【0031】
ここで、2種類の変換関数と、それにより得られる処理済画像信号の例を示す。図8は変換関数の一例を示す図であり、No.1が図6のバンドパス信号(Sorg−Sus1)を変換する変換関数、No.2が図6のバンドパス信号(Sus1−Sus2)を変換する変換関数、以下No.が大きい変換関数ほど低い周波数帯域のバンドパス信号を変換する関数を表している。この6つの変換関数により、図7に示される6個のバンドパス信号をそれぞれ変換し、前記(2)式にしたがって処理した結果得られる画像信号の入力値が小さな時の周波数特性は、図9に示されるようなものとなる。これに対し、変換関数として、図10に示されるような関数を用いた場合(但し図8同様、No.が大きい変換関数ほど低い周波数帯域のバンドパス信号を変換する関数を表している)、図7のバンドパス信号はそれぞれ図11に示されるような周波数応答特性の信号に変換され、最終的な処理済画像信号としては、図12に示されるな周波数応答特性の画像信号が得られる。
【0032】
この例から明らかであるように、変換関数を調整することにより所望の異なる周波数応答特性の画像信号を得ることができる。しかし、図8と図9、そして図10と図12を比べてわかるように、図8の変換関数を図10の変換関数に変えたら、得られる画像信号の周波数応答特性が図9のようなものから図12のようなものになるということは、必ずしも容易にわかることではない。すなわち、従来は、試行錯誤しながら変換関数の傾きなどのパラメータを調整することによって、所望の周波数応答特性が得られるような関数を見つけなければならなかった。
【0033】
本発明の画像処理方法および装置は、この問題を解決するものであるが、以下そのための手段である、図1の周波数応答特性指定手段5と変換関数定義手段4について説明する。ここで、周波数応答特性指定手段5は、ディスプレイなどの表示装置を介してユーザに周波数応答特性の入力を促して、その入力からユーザが所望する周波数応答特性を認識するものであり、具体的には、そのような処理を行うソフトウェア、入力装置、表示装置などを意味する。
【0034】
本実施の形態では、図13に示されるような操作画面が装置に付属するディスプレイに表示される。このとき画面上には、その時点で設定されている変換関数に基づいて入力信号を処理した場合に得られる信号の周波数応答特性曲線が表示されている(図13の(1))。その特性曲線上には、移動可能な6個の指示点が表示されている。すなわち、その画像処理がn個のバンドパス信号を作成して非線形処理を行うものである場合、n個の指示点が表示されるものとする。
【0035】
ユーザは、各指示点を、例えばマウスなどのポインティングデバイスを用いて所望の周波数応答に対応する点に移動する(図13の(2))。但し、表示されているままの周波数応答でよければ、必ずしも移動する必要はない。移動操作後の点(移動されなかった点を含む)は、その点が示す周波数およびその周波数における所望の周波数応答として認識され、これに基づいて特性曲線が計算しなおされ、移動操作後の点を通るような新しい特性曲線が画面上に表示される(図13の(3))。ユーザは、確認画面において、新たに表示された特性曲線が所望のものであれば、画像処理を行うことを指示し、所望のものでなければ、再度指示点の移動操作を行うことができる(図13の(4))。
【0036】
但し、周波数応答特性指定手段5の形態は上記実施の形態に限られるものではなく、例えばいくつかの周波数を順次画面に示し、ユーザがその周波数に対する所望の周波数応答特性を数値として順次入力するような形態など、他の様々な形態が考えられる。
【0037】
次に、変換関数定義手段4により行われる処理について説明する。本実施の形態では、上記のように6つのバンドパス信号を用いて処理を行っているため、上記のようなユーザインタフェースにより入力された所望の周波数応答特性は、以下のような連立方程式(5)
a1 = p1×S11+p2×S21+p3×S31+p4×S41+p5×S51+p6×S61
a2 = p1×S12+p2×S22+p3×S32+p4×S42+p5×S52+p6×S62
a3 = p1×S13+p2×S23+p3×S33+p4×S43+p5×S53+p6×S63
a4 = p1×S14+p2×S24+p3×S34+p4×S44+p5×S54+p6×S64
a5 = p1×S15+p2×S25+p3×S35+p4×S45+p5×S55+p6×S65
a6 = p1×S16+p2×S26+p3×S36+p4×S46+p5×S56+p6×S66
p1〜p6 ≧ 0 …(5)
(但し、a1〜a6 : 所望の値として指定された周波数応答
S11〜S66 : 原画像信号に基づいて作成されたバンドパス信号の周波数応答特性
p1〜p6 : 変換関数の0付近の傾きを表すパラメータ)
に、a1〜a6として当てはめらる。但し、方程式の次元は、バンドパス信号の数、すなわち指定された値の数に基づいて決まるものであり、上記連立方程式は一例にすぎない。
【0038】
変換関数定義手段4は、上記(5)式のような連立方程式を、指定された周波数応答特性a1〜a6およびバンドパス信号S11〜S66を既定値とし、パラメータp1〜p6を変数として解くことによりこのパラメータを求め、これにより変換関数を定義するものである。なお、上記(5)式において、バンドパス信号の周波数応答特性S11〜S66として表されている値は、S11〜S16(S21〜S26,…S61〜S66も同様)が1つのバンドパス信号の所定の6つの周波数の応答特性に対応するものであり、これらはボケ画像信号作成において用いられるフィルタのフィルタ係数と、補間係数から、フーリエ変換など、周波数解析において通常用いられている変換方法により求めることができる。このための計算は、上記連立方程式を解く度に行ってもよいし、あらかじめ計算しておいてその結果をメモリに記憶しておいてもよい。
【0039】
また、本実施の形態では、変換関数のパラメータp1〜p6は、それぞれ各バンドパス信号を変換する関数の0付近の傾きを表すものとしている。これは、実際に変換処理に使用される範囲が0付近であるからである。なお、その関数が非線形関数であっても、0付近の所定の範囲を線形関数として近似し、その傾きをパラメータとして上記方程式に当てはめればよい。ただし、これは一例にすぎず、変換関数を定義するためのパラメータは必ずしも傾きに限られるものではない。
【0040】
以上のようにして定義された変換関数に基づいて、非線形処理手段3により画像信号を処理すれば、ユーザが指定したとおりの周波数応答特性を有する処理済画像信号を、直ちに得ることができる。つまり、本発明により、変換関数を定義するための従来のような試行錯誤が必要なくなり、経験の少ないユーザでも、画像処理装置の機能を十分に生かして所望の画像処理を行い、所望の処理済画像信号を得ることができる。これは、経験の少ないユーザばかりでなく、経験豊富なユーザにとっても作業時間の短縮という利益をもたらすものであり、実用上の効果は極めて大きい。
【図面の簡単な説明】
【図1】本発明の画像処理装置の概略を示す図
【図2】ボケ画像信号作成処理の概要を示すブロック図
【図3】フィルタリング処理に使用されるフィルタの一例を示す図
【図4】低解像度画像信号作成処理の詳細を示す図
【図5】補間拡大処理に使用されるフィルタの一例を示す図
【図6】周波数強調処理装置の一例を示す図
【図7】バンドパス信号の周波数応答特性の一例を示す図
【図8】変換関数の一例を示す図
【図9】図8の変換関数を変換に使用した場合に得られる画像信号の周波数応答特性を示す図
【図10】変換関数の他の例を示す図
【図11】図10の変換関数により変換された図6のバンドパス信号を示す図
【図12】図10の変換関数を変換に使用した場合に得られる画像信号の周波数応答特性を示す図
【図13】周波数応答特性を指定するためのユーザインタフェースの一例を示す図
【符号の説明】
1 画像処理装置
2 ボケ画像信号(非鮮鋭マスク画像信号)作成手段
3 非線形処理手段
4 変換関数定義手段
5 周波数応答特性入力手段
10 フィルタリング処理手段
11 補間処理手段
21 減算器
22 変換器
23 演算器
Claims (2)
- 原画像を表す原画像信号に基づいて互いに周波数応答特性が異なる複数の非鮮鋭マスク画像信号を作成し、前記原画像信号および前記複数の非鮮鋭マスク画像信号に基づいて、または前記複数の非鮮鋭マスク画像信号に基づいて前記原画像信号の複数の周波数帯域ごとの信号を表す複数の帯域制限画像信号を作成し、該各帯域制限画像信号を複数の変換関数に基づいて変換して複数の変換画像信号を作成し、該各変換画像信号を積算して得た積算信号を前記原画像信号に加算することにより該原画像信号の所定の周波数成分が強調された処理済画像信号を得る画像処理方法において、
前記処理済画像信号の所望の周波数応答特性を指定し、
該指定された周波数応答特性に基づいて前記変換関数の定義パラメータを決定して該変換関数を定義することを特徴とする画像処理方法。 - 原画像を表す原画像信号に基づいて互いに周波数応答特性が異なる複数の非鮮鋭マスク画像信号を作成する非鮮鋭マスク画像信号作成手段と、前記原画像信号および前記複数の非鮮鋭マスク画像信号に基づいて、または前記複数の非鮮鋭マスク画像信号に基づいて前記原画像信号の複数の周波数帯域ごとの信号を表す複数の帯域制限画像信号を作成し、該各帯域制限原画像信号を複数の変換関数に基づいて変換して複数の変換画像信号を作成し、該各変換画像信号を積算して得た積算信号を前記原画像信号に加算することにより該原画像信号の所定の周波数成分が強調された処理済画像信号を得る非線形処理手段とからなる画像処理装置において、
前記処理済画像信号の所望の周波数応答特性を指定する周波数応答特性入力手段と、
前記指定された周波数応答特性に基づいて前記変換関数の定義パラメータを決定して該変換関数を定義する変換関数定義手段とを有することを特徴とする画像処理装置。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP21733496A JP3696339B2 (ja) | 1996-08-19 | 1996-08-19 | 画像処理方法および装置 |
US08/912,847 US6035071A (en) | 1996-08-19 | 1997-08-19 | Image processing method and apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP21733496A JP3696339B2 (ja) | 1996-08-19 | 1996-08-19 | 画像処理方法および装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH1063838A JPH1063838A (ja) | 1998-03-06 |
JP3696339B2 true JP3696339B2 (ja) | 2005-09-14 |
Family
ID=16702558
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP21733496A Expired - Fee Related JP3696339B2 (ja) | 1996-08-19 | 1996-08-19 | 画像処理方法および装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US6035071A (ja) |
JP (1) | JP3696339B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9031302B2 (en) | 2011-09-26 | 2015-05-12 | Canon Kabushiki Kaisha | Image processing apparatus configured to change a parameter determined based on user input so that an amplitude response of a processed image with respect to an original image is to be zero or more and related image processing method |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE60022237T2 (de) | 1999-01-14 | 2006-07-06 | Fuji Photo Film Co., Ltd., Minami-Ashigara | Bidverarbeitungsverfahren und -System und Aufzeichnungsmedium zur Durchfürung des Verfahrens |
JP2000268171A (ja) * | 1999-01-14 | 2000-09-29 | Fuji Photo Film Co Ltd | 画像処理方法および装置並びに記録媒体 |
JP4598115B2 (ja) * | 1999-01-26 | 2010-12-15 | 富士フイルム株式会社 | 画像処理方法および装置並びに記録媒体 |
EP1217826A3 (en) * | 2000-12-19 | 2005-05-25 | Konica Corporation | Image processing apparatus |
US8206301B2 (en) * | 2005-02-01 | 2012-06-26 | Fujifilm Corporation | Ultrasonic imaging apparatus and ultrasonic image processing apparatus, method and program |
JP2006287636A (ja) * | 2005-03-31 | 2006-10-19 | Pioneer Electronic Corp | 画質調整装置、画質調整方法およびディスプレイ装置 |
KR100790980B1 (ko) * | 2006-02-09 | 2008-01-02 | 삼성전자주식회사 | 영상 신호의 주파수 성분에 따라 후처리를 달리하는 영상신호의 후처리 회로 |
JP2010005373A (ja) * | 2008-02-14 | 2010-01-14 | Fujifilm Corp | 放射線画像補正方法および装置ならびにプログラム |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2445536B1 (fr) * | 1978-12-26 | 1989-12-22 | Fuji Photo Film Co Ltd | Procede et appareil pour traiter une image radiographique |
US4315318A (en) * | 1978-12-26 | 1982-02-09 | Fuji Photo Film Co., Ltd. | Method and apparatus for processing a radiation image |
US4819188A (en) * | 1984-11-22 | 1989-04-04 | Kabushiki Kaisha Toshiba | Image processing apparatus |
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 |
CA1316591C (en) * | 1987-10-20 | 1993-04-20 | Kazuhiro Hishinuma | Method and apparatus for radiation image processing and x-ray image processing |
US5485534A (en) * | 1990-03-28 | 1996-01-16 | Fuji Photo Film Co., Ltd. | Method and apparatus for emphasizing sharpness of image by detecting the edge portions of the image |
JPH09107479A (ja) * | 1995-09-29 | 1997-04-22 | Fuji Photo Film Co Ltd | 画像処理方法および装置 |
-
1996
- 1996-08-19 JP JP21733496A patent/JP3696339B2/ja not_active Expired - Fee Related
-
1997
- 1997-08-19 US US08/912,847 patent/US6035071A/en not_active Expired - Lifetime
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9031302B2 (en) | 2011-09-26 | 2015-05-12 | Canon Kabushiki Kaisha | Image processing apparatus configured to change a parameter determined based on user input so that an amplitude response of a processed image with respect to an original image is to be zero or more and related image processing method |
Also Published As
Publication number | Publication date |
---|---|
US6035071A (en) | 2000-03-07 |
JPH1063838A (ja) | 1998-03-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3816151B2 (ja) | 画像処理方法および装置 | |
JP3700804B2 (ja) | 画像処理方法および装置 | |
US7221806B2 (en) | Method of and system for image processing and recording medium for carrying out the method | |
US20070086671A1 (en) | Image processing apparatus | |
JP3754933B2 (ja) | 画像処理装置、画像処理システム、画像処理方法、プログラム及び記憶媒体 | |
JP3738788B2 (ja) | 画像のダイナミックレンジ圧縮処理方法および装置 | |
JP3696339B2 (ja) | 画像処理方法および装置 | |
JP3738791B2 (ja) | 画像処理方法および装置 | |
JP3814491B2 (ja) | 画像処理装置、画像処理方法、記憶媒体及びプログラム | |
JP3700798B2 (ja) | 画像処理方法および装置 | |
JP6174903B2 (ja) | 空間・階調超解像装置及びプログラム | |
JP3783116B2 (ja) | 放射線画像強調処理方法および装置 | |
JP2002183727A (ja) | 画像処理装置 | |
JP2002074356A (ja) | 画像処理方法および装置並びに記録媒体 | |
JP4598115B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
JP4014671B2 (ja) | 多重解像度変換方法および装置 | |
TW200534185A (en) | Image processing apparatus, image processing method, program, and storage medium | |
JP3669452B2 (ja) | 画像処理方法および装置 | |
JP4267159B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
JP4030655B2 (ja) | 階調変換方法および装置 | |
JP4235622B2 (ja) | 画像処理装置、画像処理システム、画像処理方法、プログラム及び記憶媒体 | |
JP2000268171A (ja) | 画像処理方法および装置並びに記録媒体 | |
JP4213337B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
JP2002183726A (ja) | 画像処理装置 | |
JP2004173078A (ja) | 画像処理装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050620 |
|
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: 20050628 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050629 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080708 Year of fee payment: 3 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080708 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090708 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090708 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100708 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110708 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110708 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120708 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120708 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130708 Year of fee payment: 8 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |