[go: up one dir, main page]

JPH0991423A - Image processing method and processor - Google Patents

Image processing method and processor

Info

Publication number
JPH0991423A
JPH0991423A JP7250857A JP25085795A JPH0991423A JP H0991423 A JPH0991423 A JP H0991423A JP 7250857 A JP7250857 A JP 7250857A JP 25085795 A JP25085795 A JP 25085795A JP H0991423 A JPH0991423 A JP H0991423A
Authority
JP
Japan
Prior art keywords
function
signal
sorg
sus
difference 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.)
Pending
Application number
JP7250857A
Other languages
Japanese (ja)
Inventor
Wataru Ito
渡 伊藤
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.)
Fujifilm Holdings Corp
Original Assignee
Fuji Photo Film Co Ltd
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 Fuji Photo Film Co Ltd filed Critical Fuji Photo Film Co Ltd
Priority to JP7250857A priority Critical patent/JPH0991423A/en
Publication of JPH0991423A publication Critical patent/JPH0991423A/en
Pending legal-status Critical Current

Links

Landscapes

  • Image Processing (AREA)
  • Facsimile Image Signal Circuits (AREA)

Abstract

PROBLEM TO BE SOLVED: To improve the observed image reading performance of a processed image by setting an absolute value of a function depending on a difference signal based on the difference signal between an original image signal and an unsharpened mask signal. SOLUTION: When an original image signal Sorg obtained by a reader or the like is given to a low pass filter 11, the filter 11 sets a prescribed fog mask to the original image signal Sorg and a fog mask signal Sus is outputted to a 1st arithmetic element 13. The arithmetic element 13 calculates a difference signal between the original image signal Sorg and the fog mask signal Sus and provides an output of the result to a conversion table 12. The conversion table 12 converts the received difference signal into an emphasis signal f(Sorg-Sus) and provides an output of it to a 2nd arithmetic means 14. The 2nd arithmetic element 14 adds the original image signal and the emphasis signal to provide an output of a processed image signal Sproc. As a result, the processed image signal is obtained, in which the sharpness of a desired image part is emphasized.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【発明の属する技術分野】本発明は画像処理方法および
装置に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an image processing method and apparatus.

【0002】[0002]

【従来の技術】従来より、種々の画像取得方法により得
られた画像を表す画像信号に対して、階調処理や周波数
処理等の画像処理を施し、画像の観察読影性能を向上さ
せることが行われている。特に人体を被写体とした放射
線画像のような医用画像の分野においては、医師等の専
門家が、得られた画像に基づいて患者の疾病や傷害の有
無を的確に診断する必要があり、その画像の読影性能を
向上させる画像処理は不可欠なものとなっている。
2. Description of the Related Art Conventionally, an image signal representing an image obtained by various image acquisition methods has been subjected to image processing such as gradation processing and frequency processing to improve the observation / interpretation performance of the image. It is being appreciated. Particularly in the field of medical images such as radiographic images of the human body, a specialist such as a doctor needs to accurately diagnose the presence or absence of a disease or injury of a patient based on the obtained image, and the image Image processing that improves the image interpretation performance of is becoming indispensable.

【0003】この画像処理のうち、いわゆる周波数強調
処理としては、例えば特開昭53−163571号に示されるよ
うに、原画像の濃度値等の画像信号(オリジナル画像信
号という)Sorg を、 Sproc=Sorg +β×(Sorg −Sus) (4) なる画像信号Sprocに変換するものが知られている。こ
こでβはオリジナル画像信号Sorg (一例としては濃度
値など)に依存する周波数強調係数、Susは非鮮鋭マス
ク(いわゆるボケマスク)信号である。このボケマスク
信号Susは、2次元に配置された画素に対してオリジナ
ル画像信号Sorg を中心画素とするN列×N行(Nは奇
数)の画素マトリクスからなるマスク、すなわちボケマ
スクを設定し、 Sus=(ΣSorg )/N2 (5) 等として求められる超低空間周波数成分である。
As the so-called frequency enhancement processing in this image processing, for example, as shown in Japanese Patent Laid-Open No. 163571/1978, an image signal (original image signal) Sorg such as a density value of an original image is represented by Sproc = Sorg + β × (Sorg−Sus) (4) There is known one for converting into an image signal Sproc. Here, β is a frequency enhancement coefficient that depends on the original image signal Sorg (as an example, a density value), and Sus is a non-sharp mask (so-called blur mask) signal. This blur mask signal Sus sets a mask composed of a pixel matrix of N columns × N rows (N is an odd number) with the original image signal Sorg as the central pixel for pixels arranged two-dimensionally, that is, a blur mask, and Sus = It is an ultra-low spatial frequency component obtained as (ΣSorg) / N 2 (5).

【0004】式(4)の第2項括弧内の値(Sorg −S
us)は、オリジナル画像信号から超低空間周波数成分で
あるボケマスク信号を減算したものであるから、オリジ
ナル画像信号のうちの、超低空間周波数成分を除去した
超低空間周波数よりも高い周波数成分Sspを意味する。
この比較的高い周波数成分Sspに周波数強調係数βを乗
じたうえで、オリジナル画像信号Sorg を加算すること
により、オリジナル画像信号Sorg のうちこの比較的高
い周波数成分Ssp、例えば画像信号の変化が急峻なエッ
ジ部などが相対的に強調された処理済画像信号Sprocを
得ることができる。
The value in the second term parenthesis of the equation (4) (Sorg-S
us) is obtained by subtracting the blur mask signal, which is an ultra-low spatial frequency component, from the original image signal, so the frequency component Ssp higher than the ultra-low spatial frequency of the original image signal from which the ultra-low spatial frequency component has been removed. Means
By multiplying the relatively high frequency component Ssp by the frequency enhancement coefficient β and then adding the original image signal Sorg, the relatively high frequency component Ssp of the original image signal Sorg, for example, the image signal changes sharply. It is possible to obtain the processed image signal Sproc in which the edges and the like are relatively emphasized.

【0005】さらに特開昭56−104645号には、下記式
(6)で表される演算をおこなって上記周波数成分Ssp
を強調する技術が開示されている。
Further, in Japanese Patent Application Laid-Open No. 56-104645, the frequency component Ssp is calculated by performing the operation represented by the following equation (6).
The technique which emphasizes is disclosed.

【0006】 Sproc=Sorg +F(x) (6) (ただし、x=Sorg −Sus、F(x)は|x1 |<|
2 |のときF′(x1)≧F′(x2 )≧0であり、
少なくともxのある値x0 (|x1 |<|x0 |<|x
2 |)を境にしてF′(x1 )>F′(x2 )となる単
調増加関数) この技術は、F′(x)すなわちSorg とSusの点の関
数の1次微分が|x|が大きい程小さくなるようにして
いるので|x|が大きいところでの周波数強調の増加率
が抑制され、これによって差信号が大きい程周波数強調
の程度を増大する一方、差信号が大きいところではその
強調の程度の増大を押えて、偽輪郭の発生を防止すると
いうものである。
Sproc = Sorg + F (x) (6) (where x = Sorg−Sus, F (x) is | x 1 | <|
When x 2 |, F ′ (x 1 ) ≧ F ′ (x 2 ) ≧ 0, and
A value x 0 (| x 1 | <| x 0 | <| x
2 |) as a boundary, F ′ (x 1 )> F ′ (x 2 ) is a monotonically increasing function) This technique uses F ′ (x), that is, the first derivative of the function of the points Sorg and Sus is | x Since | is made smaller as | becomes larger, the rate of increase of frequency emphasis at a large | x | is suppressed. As a result, the degree of frequency emphasis is increased as the difference signal is larger, while that at a large difference signal is increased. This is to prevent the generation of false contours by suppressing the increase in the degree of emphasis.

【0007】すなわち、放射線画像において差信号の小
さい部分ではボケマスク処理による周波数強調が普通に
行なわれ、差信号の大きい部分(たとえば骨と筋肉の境
界、軟部とガス部の境界、胃のBa (バリウム)充填部
とその周囲との境界血管造影における血管陰影等)では
周波数強調の程度の増大が押えられて偽輪郭の発生を防
止することができる。
That is, in a radiographic image, a frequency difference is generally emphasized by a blur masking process in a portion having a small difference signal, and a portion having a large difference signal (for example, a boundary between a bone and a muscle, a boundary between a soft portion and a gas portion, a stomach Ba (barium). ) In the boundary angiography between the filling part and its surroundings, blood vessel shadows, etc.), the increase in the degree of frequency enhancement is suppressed, and the occurrence of false contours can be prevented.

【0008】[0008]

【発明が解決しようとする課題】ところで上述の特開昭
56−104645号に開示された技術によると、F(x)は単
調増加関数であるため、オリジナル画像信号Sorg と非
鮮鋭マスク信号Susとの差信号の増大にしたがってF
(x)も増大する。ここで、原画像においてすでにある
程度鮮鋭度の高くコントラストの明瞭な画像部分につい
ては原画像そのままでも観察読影性は十分であるため本
来周波数強調を行う必要はないが、当該画像部分に対応
する差信号は比較的大きい値を示すため、上記式(6)
にしたがって強調処理したのでは当該コントラストの明
瞭な部分も強調され、その強調度合いは、原画像の鮮鋭
度が低いために強調しようとする所望の画像部分の強調
度合いよりも強いものとなり、その結果当該コントラス
トの明瞭な画像部分は過度の強調によりオーバーシュー
トやアンダーシュートを生じて偽輪郭(アーティファク
ト)を生じ、所望の画像部分の観察読影性能を却って低
下させることとなる。
The above-mentioned Japanese Patent Laid-Open No.
According to the technique disclosed in No. 56-104645, since F (x) is a monotonically increasing function, F (x) increases as the difference signal between the original image signal Sorg and the unsharp mask signal Sus increases.
(X) also increases. Here, in the original image, for an image part having a high degree of sharpness and a clear contrast, it is not necessary to originally perform frequency enhancement because the observation and interpretation are sufficient even when the original image is as it is, but the difference signal corresponding to the image part is not necessary. Is a relatively large value, the above equation (6)
According to the emphasizing process according to, the part with clear contrast is also emphasized, and the degree of emphasis becomes stronger than the degree of emphasis of the desired image part to be emphasized because the sharpness of the original image is low. The image portion having a clear contrast causes overshoot or undershoot due to excessive emphasis to generate false contours (artifacts), which deteriorates the observation and interpretation performance of the desired image portion.

【0009】一方、オリジナル画像信号Sorg と非鮮鋭
マスク信号Susとの差信号である高周波成分は、高周波
領域の放射線ノイズ等の高周波ノイズも含まれており、
この高周波ノイズは強調処理しないようにすることが画
像の観察読影性能を向上させるうえで重要である。
On the other hand, the high frequency component, which is the difference signal between the original image signal Sorg and the non-sharp mask signal Sus, contains high frequency noise such as radiation noise in the high frequency region.
It is important not to enhance the high-frequency noise in order to improve the observation / interpretation performance of the image.

【0010】本発明は上記事情に鑑みなされたものであ
って、所望の画像部分の強調を行いつつ、高周波の放射
線ノイズやすでにある程度明瞭なコントラストを有する
画像部分の強調を抑制して観察読影性能を向上した処理
済画像を表す処理済画像信号を得る画像処理方法および
装置を提供することを目的とするものである。
The present invention has been made in view of the above circumstances, and while enhancing a desired image portion, high-frequency radiation noise and enhancement of an image portion having a clear contrast to some extent are suppressed to perform observation and interpretation performance. It is an object of the present invention to provide an image processing method and apparatus for obtaining a processed image signal representing a processed image with improved image quality.

【0011】[0011]

【課題を解決するための手段】本発明の画像処理方法
は、画像を表すオリジナル画像信号Sorg の非鮮鋭マス
ク信号Susを求め、前記オリジナル画像信号Sorg と前
記非鮮鋭マスク信号Susとの差信号を求め、該差信号に
依存する関数f(x)を用いて、下記式(1)にしたが
った強調処理を施すことによって所望の画像部分を周波
数強調する画像処理方法において、 Sproc=Sorg +f(x) (1) (ただし、f(x)は差信号に依存した関数を表す;x
=Sorg −Sus) 前記所望の画像部分に対応する前記差信号の絶対値が第
1の閾値と該第1の閾値より大きい第2の閾値との間の
範囲にあるときに、前記差信号の絶対値が第1の閾値よ
り小さい範囲における前記関数f(x)の値の絶対値
を、前記所望の画像部分に対応する前記関数f(x)の
値より小さく、かつ前記差信号の絶対値が第2の閾値よ
り大きい範囲における前記関数f(x)の値の絶対値
を、前記所望の画像部分に対応する前記関数f(x)の
値より小さく設定することを特徴とするものである。
According to the image processing method of the present invention, an unsharp mask signal Sus of an original image signal Sorg representing an image is obtained, and a difference signal between the original image signal Sorg and the unsharp mask signal Sus is calculated. In the image processing method in which the desired image part is frequency-emphasized by performing the enhancement process according to the following equation (1) using the function f (x) depending on the difference signal, Sproc = Sorg + f (x ) (1) (where f (x) represents a function depending on the difference signal; x
= Sorg-Sus) of the difference signal when the absolute value of the difference signal corresponding to the desired image portion is in a range between a first threshold and a second threshold greater than the first threshold. The absolute value of the value of the function f (x) in the range where the absolute value is smaller than the first threshold is smaller than the value of the function f (x) corresponding to the desired image portion, and the absolute value of the difference signal. Is set to be smaller than the value of the function f (x) corresponding to the desired image portion in the range where is larger than the second threshold value. .

【0012】なお関数fは具体的には、図2に示すよう
に、(Sorg −Sus)≧0においてはf≧0であって、
差信号(Sorg −Sus)が高周波ノイズに対応する第1
の所定値α1 (α1 >0)までの範囲内においては関数
fの値が小さく、差信号(Sorg −Sus)が所望の画像
部分に対応する第2の所定値α2 (α2 >α1 )までは
単調に増加し、第2の所定値α2 以上の範囲においては
単調減少し、(Sorg−Sus)がすでにある程度明瞭な
コントラストを有する画像部分に対応する第3の所定値
α3 (α3 >α2 )以上の範囲は関数fの値が小さくな
るように設定されたものである。
The function f is specifically f ≧ 0 when (Sorg-Sus) ≧ 0, as shown in FIG.
The first difference signal (Sorg-Sus) corresponds to high frequency noise
Within a range up to a predetermined value α 11 > 0) of the function f, the value of the function f is small, and the difference signal (Sorg-Sus) has a second predetermined value α 22 >) corresponding to the desired image portion. up to α 1 ) and monotonically decreases in the range of the second predetermined value α 2 or more, and (Sorg-Sus) already corresponds to the third predetermined value α corresponding to an image part having a clear contrast. The range of 33 > α 2 ) or more is set so that the value of the function f becomes small.

【0013】また関数fを、画像の高周波ノイズを表す
指標値(例えば分散値σ2 等)に応じた補正用指標関数
gを用いて、このノイズに対応した画像部分の強調を一
層抑制するように、下記式(2)にしたがって補正する
ようにしてもよい。
Further, a function f is used as a correction index function g corresponding to an index value (for example, a variance value σ 2 etc.) representing high frequency noise of the image to further suppress the enhancement of the image portion corresponding to this noise. Alternatively, the correction may be performed according to the following equation (2).

【0014】 f′(x)=g(x)×f(x) (2) (ただし、g(x)は0≦g(x)≦1の範囲で単調に
増加する補正用指標関数、f′(x)は式(1)の関数
f(x)に代えて適用される補正後の関数を表す) このように補正することにより式(1)は下記式
(1′)と表すことができる。
F ′ (x) = g (x) × f (x) (2) (where g (x) is a correction index function that monotonically increases in the range 0 ≦ g (x) ≦ 1, f ′ (X) represents a corrected function that is applied in place of the function f (x) in the equation (1)) By performing the correction in this way, the equation (1) may be expressed as the following equation (1 ′). it can.

【0015】 Sproc=Sorg +g(x)×f(x) (1′) なお、上記分散値σ2 は所定の大きさ(例えばm行×n
列)の画素マトリックスからなるマスクを設定し、この
マスク内の全画素のオリジナル画像信号Sorgの平均
値Smeanを求め、下記演算(式(7)参照)により
求めるものであるが、このマスクを非鮮鋭マスクと同一
のものとすれば、Smean=Susとなり演算処理を下記式
(8)に示すように簡単化することができる。
Sproc = Sorg + g (x) × f (x) (1 ′) Note that the variance value σ 2 has a predetermined value (for example, m rows × n).
(Column) pixel mask is set, the average value Smean of the original image signals Sorg of all the pixels in this mask is calculated, and the average value Smean is calculated by the following calculation (see equation (7)). If it is the same as the sharp mask, Smean = Sus and the arithmetic processing can be simplified as shown in the following formula (8).

【0016】 σ2 ={Σ(Sorg −Smean)2 }/(m・n) (7) σ2 =(Σ|Sorg −Sus|2 )/N2 (8) 分散値σ2 を画像のノイズを表す指標値とした補正用の
指標関数gの具体例を図4に示す。図示の補正用指標関
数gは、分散値σ2 が高周波ノイズ領域に対応する範囲
内においては関数gの出力が小さくなるように、分散値
σ2 が信号領域に対応する範囲内においては関数gの出
力が1に近くなるように設定されたものである。
Σ 2 = {Σ (Sorg −Smean) 2 } / (m · n) (7) σ 2 = (Σ | Sorg −Sus | 2 ) / N 2 (8) The variance value σ 2 is the noise of the image. FIG. 4 shows a specific example of the correction index function g that is an index value representing The correction index function g shown in the figure has a function g within the range in which the variance value σ 2 corresponds to the signal region such that the output of the function g is small in the range in which the variance value σ 2 corresponds to the high frequency noise region. Is set so that the output of is close to 1.

【0017】具体的には、補正用の関数gは下記式
(3)で表すことができる。
Specifically, the correction function g can be expressed by the following equation (3).

【0018】 g(x)=G(Σ|Sorg −Sus|) (3) 本発明の画像処理装置は、上記本発明の画像処理方法を
実施するための具体的な装置であって、画像を表すオリ
ジナル画像信号Sorg の非鮮鋭マスク信号Susを求める
非鮮鋭マスク信号演算手段と、前記オリジナル画像信号
Sorg と前記非鮮鋭マスク信号Susとの差信号を求める
差信号演算手段と、該差信号を該差信号に依存する関数
f(x)に変換する変換テーブルと、下記式(1)にし
たがった強調処理を施す強調処理手段とを備え、前記画
像のうち所望の画像部分を周波数強調する画像処理装置
において、 Sproc=Sorg +f(x) (1) (ただし、f(x)は差信号に依存した関数を表す;x
=Sorg −Sus) 前記所望の画像部分に対応する前記差信号の絶対値が第
1の閾値と該第1の閾値より大きい第2の閾値との間の
範囲にあるときに、前記差信号の絶対値が第1の閾値よ
り小さい範囲における前記関数f(x)の値の絶対値
が、前記所望の画像部分に対応する前記関数f(x)の
値より小さく、かつ前記差信号の絶対値が第2の閾値よ
り大きい範囲における前記関数f(x)の値の絶対値
が、前記所望の画像部分に対応する前記関数f(x)の
値より小さく設定されていることを特徴とするものであ
る。
G (x) = G (Σ | Sorg −Sus |) (3) The image processing apparatus of the present invention is a specific apparatus for carrying out the above-described image processing method of the present invention. The non-sharp mask signal calculating means for obtaining the non-sharp mask signal Sus of the original image signal Sorg represented, the difference signal calculating means for obtaining the difference signal between the original image signal Sorg and the non-sharp mask signal Sus, and the difference signal An image processing that includes a conversion table that converts a function f (x) depending on a difference signal, and an emphasis processing unit that performs an emphasis process according to the following expression (1), and frequency-enhances a desired image portion of the image. In the apparatus, Sproc = Sorg + f (x) (1) (where f (x) represents a function dependent on the difference signal; x
= Sorg-Sus) of the difference signal when the absolute value of the difference signal corresponding to the desired image portion is in a range between a first threshold and a second threshold greater than the first threshold. The absolute value of the value of the function f (x) in the range where the absolute value is smaller than the first threshold value is smaller than the value of the function f (x) corresponding to the desired image portion, and the absolute value of the difference signal. Is set to be smaller than a value of the function f (x) corresponding to the desired image portion in a range in which is larger than a second threshold value. Is.

【0019】また、前記差信号を、前記画像のうち該差
信号に依存する高周波ノイズを表す補正用指標関数g
(x)に変換する補正用変換テーブルと、下記式(2)
にしたがって前記関数f(x)をさらに備えた構成を採
用することもできる。
In addition, the difference signal is used as a correction index function g that represents high-frequency noise in the image that depends on the difference signal.
A conversion table for correction to be converted into (x) and the following formula (2)
Accordingly, a configuration further including the function f (x) can be adopted.

【0020】 f′(x)=g(x)×f(x) (2) (ただし、g(x)は0≦g(x)≦1の範囲で単調に
増加する補正用指標関数、f′(x)は式(1)の関数
f(x)に代えて適用される補正後の関数を表す) なお上記強調関数f(x)としては、上記本発明の画像
処理方法における強調関数と同様に、具体的には図2に
示すものなどを適用することができる。
F ′ (x) = g (x) × f (x) (2) (where g (x) is a correction index function that monotonically increases in the range of 0 ≦ g (x) ≦ 1, f ′ (X) represents a corrected function that is applied in place of the function f (x) of the equation (1). Note that the enhancement function f (x) is the enhancement function in the image processing method of the present invention. Similarly, specifically, the one shown in FIG. 2 or the like can be applied.

【0021】さらに、この強調関数f(x)を、オリジ
ナル画像信号のうちノイズを表す分散値σ2 等の指標値
に応じた補正用指標関数g(x)により、このノイズに
対応した画像部分の強調が抑制されるように、下記式
(2)にしたがって補正する補正手段をさらに備えた構
成を採ることもできる。
Further, the enhancement function f (x) is converted into an image portion corresponding to the noise by a correction index function g (x) corresponding to an index value such as a variance value σ 2 representing noise in the original image signal. It is also possible to adopt a configuration further including a correction unit that corrects according to the following formula (2) so that the emphasis of 1 is suppressed.

【0022】 f′(x)=g(x)×f(x) (2) (ただし、g(Sorg −Sus)は指標関数;0≦g(S
org −Sus)≦1、f′(Sorg −Sus)は式(1)の
関数f(Sorg −Sus)に代えて適用される補正後の関
数を表す) 具体的には、指標値に応じた関数gは下記式(3)で表
すことができる。
F ′ (x) = g (x) × f (x) (2) (where g (Sorg-Sus) is an index function; 0 ≦ g (S
org-Sus) ≤ 1, f '(Sorg-Sus) represents a corrected function applied in place of the function f (Sorg-Sus) of the equation (1)) Specifically, depending on the index value The function g can be expressed by the following equation (3).

【0023】 g(x)=G(Σ|Sorg −Sus|2 ) (3)G (x) = G (Σ | Sorg−Sus | 2 ) (3)

【0024】[0024]

【発明の効果】本発明の画像処理方法,装置によれば、
オリジナル画像信号Sorg とその非鮮鋭マスク信号との
差信号が小さい範囲においては、その差信号の中に占め
る高周波ノイズの割合が高いため、これの強調処理を抑
制することによって画像の観察読影性能を向上せしめ、
一方、差信号が所望の画像部分に対応する差信号より大
きい範囲における原画像はある程度明瞭なコントラスト
を有するため、これの強調処理を抑制することによって
当該部分に偽輪郭が生じるのを防止し画像の観察読影性
能を向上させることができる。
According to the image processing method and apparatus of the present invention,
In the range where the difference signal between the original image signal Sorg and its unsharp mask signal is small, the high-frequency noise occupies a large part of the difference signal. Improve
On the other hand, since the original image in a range where the difference signal is larger than the difference signal corresponding to the desired image portion has a clear contrast to some extent, by suppressing the enhancement processing of the original image, it is possible to prevent the false contour from being generated in the portion. It is possible to improve the observation / interpretation performance.

【0025】[0025]

【発明の実施の形態】以下、本発明の画像処理方法を具
体的に実施するための画像処理装置の実施の形態につい
て図面を用いて説明する。
BEST MODE FOR CARRYING OUT THE INVENTION Embodiments of an image processing apparatus for specifically implementing the image processing method of the present invention will be described below with reference to the drawings.

【0026】図1は本発明の画像処理装置の第1の実施
の形態を示す概略ブロック図である。図示の画像処理装
置は、画像を表すオリジナル画像信号(濃度値)Sorg
の、超低空間周波数に対応する非鮮鋭マスク信号Susを
求めるローパスフィルタ11と、このローパスフィルタ11
により得られた非鮮鋭マスク信号Susとオリジナル画像
信号Sorg との差信号(Sorg −Sus)を算出する第1
演算子13と、得られた差信号(Sorg −Sus)が入力さ
れて、その差信号(Sorg −Sus)に対応づけされた所
定の強調信号f(Sorg −Sus)を出力する変換テーブ
ル12と、強調信号f(Sorg −Sus)ともとのオリジナ
ル画像信号Sorg とを加算処理して処理済画像信号Spr
ocを算出する第2演算子14とを備えた構成である。
FIG. 1 is a schematic block diagram showing a first embodiment of the image processing apparatus of the present invention. The illustrated image processing apparatus uses an original image signal (density value) Sorg representing an image.
, A low-pass filter 11 for obtaining an unsharp mask signal Sus corresponding to an ultra-low spatial frequency, and this low-pass filter 11
First, a difference signal (Sorg-Sus) between the unsharp mask signal Sus and the original image signal Sorg obtained by
An operator 13 and a conversion table 12 for inputting the obtained difference signal (Sorg-Sus) and outputting a predetermined emphasis signal f (Sorg-Sus) associated with the difference signal (Sorg-Sus). , The enhancement signal f (Sorg-Sus) and the original image signal Sorg and the processed image signal Spr.
and a second operator 14 for calculating oc.

【0027】ローパスフィルタ11は入力されたオリジナ
ル画像信号Sorg に対して、例えば3列×3行の画素マ
トリクスからなるボケマスクを設定し、下記式(5)
(N=3に設定)にしたがって得られたボケマスク信号
Susを画像強調手段16に出力する。
The low-pass filter 11 sets a blur mask consisting of, for example, a pixel matrix of 3 columns × 3 rows for the input original image signal Sorg, and the following equation (5)
The blur mask signal Sus obtained according to (set to N = 3) is output to the image enhancing means 16.

【0028】 Sus=(ΣSorg )/N2 (5) なお、ボケマスクとしては式(5)に示すようにマスク
内の画素値の単純平均を用いるものの他、中心画素から
の距離に応じてマスク内の画素値の重み付けを変化させ
たものを用いることもできる。
Sus = (ΣS org) / N 2 (5) As the blur mask, a simple average of the pixel values in the mask is used as shown in Expression (5), or the mask in the mask is calculated according to the distance from the center pixel. It is also possible to use one in which the weighting of the pixel value of is changed.

【0029】変換テーブル12は例えば図2に示すような
関数形状にしたがって、差信号(Sorg −Sus)と強調
信号f(Sorg −Sus)とが対応付けされている。
In the conversion table 12, the difference signal (Sorg-Sus) and the enhancement signal f (Sorg-Sus) are associated with each other in accordance with the function shape as shown in FIG.

【0030】この関数形状は、(Sorg −Sus)≧0に
おいてはf≧0であって、(Sorg−Sus)が原画像の
高周波ノイズに対応する第1の所定値α1 (α1 >0)
までの範囲内においては関数fの値が小さく、(Sorg
−Sus)が所望の画像部分に対応する第2の所定値α2
(α2 >α1 )までは単調に増加し、第2の所定値α2
以上の範囲においては単調減少し、(Sorg −Sus)が
すでにある程度明瞭なコントラストを有する画像部分に
対応する第3の所定値α3 (α3 >α2 )以上の範囲は
関数fの値が小さくなるように設定されたものである。
なお、(Sorg−Sus)≦0における関数形状は、(So
rg −Sus)≧0の範囲における関数形状を原点につい
て対称にしたものとなっている。
This function shape is such that f ≧ 0 when (Sorg-Sus) ≧ 0, and (Sorg−Sus) corresponds to the high frequency noise of the original image, the first predetermined value α 11 > 0). )
The value of the function f is small in the range up to (Sorg
-Sus) is the second predetermined value α 2 corresponding to the desired image portion
It increases monotonically until (α 2 > α 1 ) and reaches the second predetermined value α 2
In the range above, the value of the function f is monotonically decreasing, and in the range above the third predetermined value α 33 > α 2 ) where (Sorg-Sus) already corresponds to an image part having a clear contrast, the value of the function f is It is set to be small.
The function shape when (Sorg-Sus) ≦ 0 is (So
The function shape in the range of rg-Sus) ≧ 0 is symmetrical about the origin.

【0031】次に本実施形態の画像処理装置の作用につ
いて説明する。
Next, the operation of the image processing apparatus of this embodiment will be described.

【0032】まず、所定の画像情報読取装置等により得
られたオリジナル画像信号Sorg がローパスフィルタ1
1、第1演算子13、第2演算子14にそれぞれ入力され
る。
First, the original image signal Sorg obtained by a predetermined image information reading device is the low-pass filter 1
1, the first operator 13 and the second operator 14, respectively.

【0033】ローパスフィルタ11は入力されたオリジナ
ル画像信号Sorg に対して、例えば3列×3行の画素マ
トリクスからなるボケマスクを設定し、上記式(5)に
したがって得られたボケマスク信号Susを第1演算子13
に出力する。
The low-pass filter 11 sets, for the input original image signal Sorg, a blur mask consisting of, for example, a pixel matrix of 3 columns × 3 rows, and outputs the blur mask signal Sus obtained according to the above equation (5) as the first blur mask signal Sus. Operator 13
Output to

【0034】なお、ボケマスクとしては式(5)に示す
ようにマスク内の画素値の単純平均を用いるものの他、
前述したように中心画素からの距離に応じてマスク内の
画素値の重み付けを適宜変化させてもよい。
As the blur mask, a simple average of pixel values in the mask is used as shown in equation (5),
As described above, the weighting of pixel values in the mask may be appropriately changed according to the distance from the central pixel.

【0035】第1演算子13は、入力されたオリジナル画
像信号Sorg およびボケマスク信号Susに基づいて、こ
れらの差信号(Sorg −Sus)を算出し、その結果を変
換テーブル12に出力する。
The first operator 13 calculates the difference signal (Sorg-Sus) based on the input original image signal Sorg and blur mask signal Sus, and outputs the result to the conversion table 12.

【0036】変換テーブル12は入力された差信号(Sor
g −Sus)をその差信号に応じた強調信号f(Sorg −
Sus)に変換して第2演算子14に出力する。
The conversion table 12 receives the difference signal (Sor
g-Sus) is an emphasized signal f (Sorg-
Sus) and output to the second operator 14.

【0037】この強調信号f(Sorg −Sus)は、前述
したような関数形状にしたがって差信号(Sorg −Su
s)と強調信号f(Sorg −Sus)とが対応付けされて
いるため、差信号が放射線ノイズなどその大きさが極め
て小さい範囲、すなわち第1の所定値の範囲(−α1
Sorg −Sus<α1 )における強調信号は略ゼロに近い
小さい値に押さえられたものとなり、この差信号が放射
線ノイズによるものより大きく、所望の構造物の空間周
波数に対応する信号成分に応じた範囲、すなわち差信号
が第1の所定値と第3の所定値の範囲(α1 ≦Sorg −
Sus≦α3 または−α3 ≦Sorg −Sus≦−α1 )にお
ける強調信号はある程度大きな値を出力し、所望の構造
物よりも鮮鋭度が高い画像部分に対応する範囲(α3
Sorg −SusまたはSorg −Sus<−α3 )では強調信
号は略ゼロに近い小さい値に押さえられる。
This emphasized signal f (Sorg-Sus) is a difference signal (Sorg-Su) according to the function shape as described above.
s) and the enhancement signal f (Sorg-Sus) are associated with each other, the difference signal has a very small range such as radiation noise, that is, a range of the first predetermined value (-α 1 <
The enhancement signal in Sorg-Sus <α 1 ) is suppressed to a small value close to zero, and this difference signal is larger than that due to radiation noise, and depends on the signal component corresponding to the spatial frequency of the desired structure. The range, that is, the range of the difference signal between the first predetermined value and the third predetermined value (α 1 ≦ Sorg −
The enhancement signal in Sus ≦ α 3 or −α 3 ≦ Sorg −Sus ≦ −α 1 ) outputs a relatively large value, and the range (α 3 <which corresponds to an image portion having a sharpness higher than that of the desired structure is used.
Sorg-Sus or Sorg -Sus <-α 3) the enhancement signal is suppressed to a small value close to approximately zero.

【0038】第2演算子14は入力されたオリジナル画像
信号Sorg と強調信号f(Sorg −Sus)とを加算処理
し、下記式(1)に示す処理済画像信号Sprocを出力す
る。
The second operator 14 performs addition processing on the input original image signal Sorg and the emphasized signal f (Sorg-Sus), and outputs a processed image signal Sproc shown in the following equation (1).

【0039】 Sproc=Sorg +f(Sorg −Sus) (1) この結果、処理済画像信号Sprocは、高周波の放射線ノ
イズおよびすでにある程度明瞭なコントラストを有する
画像部分の強調が抑制され、かつ所望の画像部分の鮮鋭
度が強調されたものとなる。
Sproc = Sorg + f (Sorg−Sus) (1) As a result, in the processed image signal Sproc, high-frequency radiation noise and enhancement of an image portion having a clear contrast to some extent are suppressed, and a desired image portion is suppressed. The sharpness of is emphasized.

【0040】したがって本実施形態の画像処理装置によ
れば、原画像においてすでにある程度明瞭なコントラス
トを有する画像部分については強調が抑制されるため、
当該画像部分の過度の強調によるオーバーシュート、ア
ンダーシュートを防止し、アーティファクトの形成を低
減して、観察読影性能の優れた可視画像を再生するのに
適した処理済画像信号を得ることができる。
Therefore, according to the image processing apparatus of the present embodiment, the enhancement is suppressed for the image portion of the original image which already has a clear contrast to some extent.
It is possible to prevent overshoot and undershoot due to excessive enhancement of the image portion, reduce the formation of artifacts, and obtain a processed image signal suitable for reproducing a visible image with excellent observation and interpretation performance.

【0041】図3は本発明の画像処理装置の第2の実施
の形態を示す概略ブロック図である。図示の画像処理装
置は、図1に示した第1の実施形態の画像処理装置に対
して、強調信号f(Sorg −Sus)を補正する補正手段
である分散値演算手段15と補正値テーブル16と第3演算
子17とをさらに備えた構成である。
FIG. 3 is a schematic block diagram showing a second embodiment of the image processing apparatus of the invention. The image processing apparatus shown in the figure is different from the image processing apparatus of the first embodiment shown in FIG. 1 in that the variance value calculating means 15 and the correction value table 16 are correction means for correcting the enhancement signal f (Sorg-Sus). And a third operator 17 are further provided.

【0042】すなわち、分散値演算手段15は上述の差信
号(Sorg −Sus)に基づいて、画像の高周波ノイズ成
分を含む指標値として分散値σ2 を求め、補正値テーブ
ル16はこの分散値σ2 に応じた補正係数g(0≦g≦
1)を求めるものであり、第3演算子17は得られた補正
係数gを強調信号f(Sorg −Sus)に乗じてその積を
求めるものである。なお分散値σ2 は下記式(8)によ
り求められ、これは差信号(Sorg −Sus)に応じた関
数g(Sorg −Sus)として表すことができる。
That is, the variance value calculating means 15 obtains the variance value σ 2 as an index value including the high frequency noise component of the image based on the above difference signal (Sorg-Sus), and the correction value table 16 uses this variance value σ. correction coefficient corresponding to 2 g (0 ≦ g ≦
1) and the third operator 17 multiplies the enhancement signal f (Sorg-Sus) by the obtained correction coefficient g to obtain the product. The variance value σ 2 is obtained by the following equation (8), which can be expressed as a function g (Sorg-Sus) corresponding to the difference signal (Sorg-Sus).

【0043】 σ2 =(Σ|Sorg −Sus|2 )/N2 (8) ここで補正値テーブル16が出力する補正係数gは、図4
に示すように、分散値σ2 が極小さい範囲においては1
を大幅に下回り、分散値σ2 がある程度以上の大きさの
範囲においては1に略等しい値を採るような関数形状と
して分散値σ2に対応付けされている。すなわち分散値
σ2 が小さい範囲においては、画像信号よりも高周波ノ
イズを含む割合が大きいため、この範囲においては強調
信号を小さく抑制する。
Σ 2 = (Σ | Sorg −Sus | 2 ) / N 2 (8) Here, the correction coefficient g output from the correction value table 16 is as shown in FIG.
As shown in, in the range where the variance value σ 2 is extremely small, 1
Is significantly lower than, and the variance value σ 2 is associated with the variance value σ 2 as a function shape that takes a value substantially equal to 1 in a range of a certain magnitude or more. That is, in the range in which the variance value σ 2 is small, the ratio of high-frequency noise is larger than that in the image signal, so that the emphasis signal is suppressed small in this range.

【0044】次に本実施形態の画像処理装置の作用につ
いて説明する。
Next, the operation of the image processing apparatus of this embodiment will be described.

【0045】まず、所定の画像情報読取装置等により得
られたオリジナル画像信号Sorg がローパスフィルタ1
1、第1演算子13、第2演算子14にそれぞれ入力され
る。
First, the original image signal Sorg obtained by a predetermined image information reading device is the low-pass filter 1
1, the first operator 13 and the second operator 14, respectively.

【0046】ローパスフィルタ11は入力されたオリジナ
ル画像信号Sorg に対して、例えば3列×3行の画素マ
トリクスからなるボケマスクを設定し、上記式(5)に
したがって得られたボケマスク信号Susを第1演算子13
に出力する。
The low-pass filter 11 sets, for the input original image signal Sorg, a blur mask consisting of, for example, a pixel matrix of 3 columns × 3 rows, and outputs the blur mask signal Sus obtained according to the above equation (5) as the first blur mask signal Sus. Operator 13
Output to

【0047】第1演算子13は、入力されたオリジナル画
像信号Sorg およびボケマスク信号Susに基づいて、こ
れらの差信号(Sorg −Sus)を算出し、その結果を変
換テーブル12および分散値演算手段12に出力する。
The first operator 13 calculates a difference signal (Sorg-Sus) based on the input original image signal Sorg and blur mask signal Sus, and the result thereof is converted into a conversion table 12 and a variance value calculating means 12. Output to.

【0048】変換テーブル12は入力された差信号(Sor
g −Sus)をその差信号に応じた強調信号f(Sorg −
Sus)に変換して第2演算子14に出力する。
The conversion table 12 displays the input difference signal (Sor
g-Sus) is an emphasized signal f (Sorg-
Sus) and output to the second operator 14.

【0049】この強調信号f(Sorg −Sus)は、差信
号が放射線ノイズなどその大きさが極小さい範囲(−α
1 <Sorg −Sus<α1 )における強調信号は略ゼロに
近い小さい値に押さえられたものとなり、この差信号が
放射線ノイズによるものより大きく、所望の空間周波数
の構造物による信号成分に対応する範囲(α1 ≦Sorg
−Sus≦α3 または−α3 ≦Sorg −Sus≦−α1 )に
おける強調信号はある程度大きな値を出力し、所望の構
造物よりもコントラストが明瞭な画像部分に対応する部
分(α3 <Sorg −SusまたはSorg −Sus<−α3
では強調信号は略ゼロに近い小さい値に押さえられ、第
3演算子17に入力される。
This emphasized signal f (Sorg-Sus) has a range (-α) in which the difference signal has a very small magnitude such as radiation noise.
The emphasized signal in 1 <Sorg-Sus <α 1 ) is suppressed to a small value close to zero, and this difference signal is larger than that due to radiation noise and corresponds to the signal component due to the structure of the desired spatial frequency. Range (α 1 ≤ Sorg
The enhancement signal in −Sus ≦ α 3 or −α 3 ≦ Sorg −Sus ≦ −α 1 ) outputs a relatively large value, and a portion corresponding to an image portion where the contrast is clearer than that of the desired structure (α 3 <Sorg -Sus or Sorg-Sus <-α 3 )
Then, the emphasis signal is suppressed to a small value close to zero and is input to the third operator 17.

【0050】一方、分散値演算手段15は入力された差信
号(Sorg −Sus)に応じた分散値σ2 を求める。ここ
で、分散値σ2 を求めるための局所的に設定されるマス
クは、演算処理の簡単化のため非鮮鋭マスクと同一のも
のとする。これによって分散値マスク内の平均値Smean
とボケマスク信号Susと等しいものとなり、得られた分
散値σ2 は補正値テーブル16に入力される。
On the other hand, the variance value calculating means 15 obtains the variance value σ 2 according to the input difference signal (Sorg-Sus). Here, the mask that is locally set for obtaining the variance value σ 2 is the same as the non-sharp mask for the sake of simplifying the arithmetic processing. This gives the mean value Smean in the variance mask.
Is equal to the blur mask signal Sus, and the obtained variance value σ 2 is input to the correction value table 16.

【0051】補正値テーブル16は入力された分散値σ2
に対応する補正係数g(Sorg −Sus)を出力する。す
なわち、分散値σ2 が極めて小さい範囲においては1を
大幅に下回る補正係数g(Sorg −Sus)を、分散値σ
2 がある程度以上の大きさの範囲においては1に略等し
い補正係数g(Sorg −Sus)をそれぞれ出力する。
The correction value table 16 stores the input variance value σ 2
The correction coefficient g (Sorg-Sus) corresponding to is output. That is, in the range where the variance value σ 2 is extremely small, the correction coefficient g (Sorg-Sus) that is significantly less than 1 is set to the variance value σ
In the range where 2 is larger than a certain level, the correction coefficient g (Sorg-Sus) which is substantially equal to 1 is output.

【0052】第3演算子17は入力された補正係数g(S
org −Sus)と強調信号f(Sorg−Sus)との積を求
め、得られた積を第2演算子14に出力する。
The third operator 17 inputs the correction coefficient g (S
org-Sus) and the emphasized signal f (Sorg-Sus) are obtained, and the obtained product is output to the second operator 14.

【0053】第2演算子14は入力されたオリジナル画像
信号Sorg と補正係数g(Sorg −Sus)により補正さ
れた強調信号f(Sorg −Sus)・g(Sorg −Sus)
とを加算処理し、下記式(1′)に示す処理済画像信号
Sprocを出力する。
The second operator 14 is the enhancement signal f (Sorg-Sus) .g (Sorg-Sus) corrected by the input original image signal Sorg and the correction coefficient g (Sorg-Sus).
And are subjected to addition processing, and the processed image signal Sproc shown in the following equation (1 ′) is output.

【0054】 Sproc=Sorg +g(Sorg −Sus)×f(Sorg −Sus) (1′) この結果、処理済画像信号Sprocは、高周波の放射線ノ
イズおよびすでにある程度明瞭なコントラストを有する
画像部分の強調が抑制され、かつ所望の画像部分の鮮鋭
度が強調されたものとなる。
Sproc = Sorg + g (Sorg−Sus) × f (Sorg−Sus) (1 ′) As a result, the processed image signal Sproc has high-frequency radiation noise and enhancement of an image portion which already has a clear contrast to some extent. It is suppressed and the sharpness of the desired image portion is emphasized.

【0055】このように本実施形態の画像処理装置によ
れば、原画像においてすでにある程度明瞭なコントラス
トを有する画像部分については強調が抑制されるため、
過度の強調によるオーバーシュート、アンダーシュート
を防止し、アーティファクトの形成を低減して、観察読
影性能の優れた可視画像を再生するのに適した画像信号
を得ることができる。
As described above, according to the image processing apparatus of this embodiment, the enhancement is suppressed for the image portion of the original image which already has a clear contrast to some extent.
It is possible to prevent overshoot and undershoot due to excessive emphasis, reduce the formation of artifacts, and obtain an image signal suitable for reproducing a visible image with excellent observation and interpretation performance.

【図面の簡単な説明】[Brief description of drawings]

【図1】本発明の画像処理装置の第1の実施の形態を示
す概略ブロック図
FIG. 1 is a schematic block diagram showing a first embodiment of an image processing apparatus of the present invention.

【図2】差信号(Sorg −Sus)と強調信号f(Sorg
−Sus)とを対応付ける変換テーブルを示すグラフ
FIG. 2 shows a difference signal (Sorg-Sus) and an emphasis signal f (Sorg).
-Sus) is a graph showing a conversion table in which

【図3】本発明の画像処理装置の第2の実施の形態を示
す概略ブロック図
FIG. 3 is a schematic block diagram showing a second embodiment of the image processing apparatus of the invention.

【図4】分散値σ2 と補正係数g(Sorg −Sus)とを
対応付ける変換テーブルを示すグラフ
FIG. 4 is a graph showing a conversion table that associates a variance value σ 2 with a correction coefficient g (Sorg-Sus).

【符号の説明】[Explanation of symbols]

11 ローパスフィルタ 12 変換テーブル 13,14,17 演算子 15 分散値演算手段 16 補正値テーブル 11 Low-pass filter 12 Conversion table 13, 14, 17 Operator 15 Variance value calculation means 16 Correction value table

Claims (6)

【特許請求の範囲】[Claims] 【請求項1】 画像を表すオリジナル画像信号Sorg の
非鮮鋭マスク信号Susを求め、前記オリジナル画像信号
Sorg と前記非鮮鋭マスク信号Susとの差信号を求め、
該差信号に依存する関数f(x)を用いて、下記式
(1)にしたがった強調処理を施すことによって所望の
画像部分を周波数強調する画像処理方法において、 Sproc=Sorg +f(x) (1) (ただし、f(x)は差信号に依存した関数を表す;x
=Sorg −Sus) 前記所望の画像部分に対応する前記差信号の絶対値が第
1の閾値と該第1の閾値より大きい第2の閾値との間の
範囲にあるときに、 前記差信号の絶対値が第1の閾値より小さい範囲におけ
る前記関数f(x)の値の絶対値を、前記所望の画像部
分に対応する前記関数f(x)の値より小さく、かつ前
記差信号の絶対値が第2の閾値より大きい範囲における
前記関数f(x)の値の絶対値を、前記所望の画像部分
に対応する前記関数f(x)の値より小さく設定するこ
とを特徴とする請求項1記載の画像処理方法。
1. An unsharp mask signal Sus of an original image signal Sorg representing an image is obtained, and a difference signal between the original image signal Sorg and the unsharp mask signal Sus is obtained,
In the image processing method for frequency-enhancing a desired image portion by performing an enhancement process according to the following formula (1) using the function f (x) depending on the difference signal, Sproc = Sorg + f (x) ( 1) (where f (x) represents a function dependent on the difference signal; x
= Sorg-Sus) When the absolute value of the difference signal corresponding to the desired image portion is in a range between a first threshold value and a second threshold value larger than the first threshold value, the difference signal The absolute value of the value of the function f (x) in the range where the absolute value is smaller than the first threshold is smaller than the value of the function f (x) corresponding to the desired image portion, and the absolute value of the difference signal. 2. The absolute value of the value of the function f (x) in the range where is larger than the second threshold value is set smaller than the value of the function f (x) corresponding to the desired image portion. The described image processing method.
【請求項2】 前記差信号に基づく、前記画像のうち高
周波ノイズを表す補正用指標関数g(x)を用いて、下
記式(2)にしたがって前記関数f(x)を補正せしめ
ることを特徴とする請求項1記載の画像処理方法。 f′(x)=g(x)×f(x) (2) (ただし、g(x)は0≦g(x)≦1の範囲で単調に
増加する補正用指標関数、f′(x)は式(1)の関数
f(x)に代えて適用される補正後の関数を表す)
2. The function f (x) is corrected according to the following equation (2) using a correction index function g (x) representing high frequency noise in the image based on the difference signal. The image processing method according to claim 1. f ′ (x) = g (x) × f (x) (2) (where g (x) is a correction index function that monotonically increases in the range of 0 ≦ g (x) ≦ 1, f ′ (x ) Represents a corrected function that is applied in place of the function f (x) in equation (1))
【請求項3】 前記補正用指標関数g(x)は、下記式
(3)に示す前記差信号の分散値に応じた関数G(x)
であることを特徴とする請求項2記載の画像処理方法。 g(x)=G(Σ|Sorg −Sus|2 ) (3)
3. The correction index function g (x) is a function G (x) corresponding to the variance value of the difference signal shown in the following equation (3).
The image processing method according to claim 2, wherein g (x) = G (Σ | Sorg-Sus | 2 ) (3)
【請求項4】 画像を表すオリジナル画像信号Sorg の
非鮮鋭マスク信号Susを求める非鮮鋭マスク信号演算手
段と、前記オリジナル画像信号Sorg と前記非鮮鋭マス
ク信号Susとの差信号を求める差信号演算手段と、該差
信号を該差信号に依存する関数f(x)に変換する変換
テーブルと、下記式(1)にしたがった強調処理を施す
強調処理手段とを備え、前記画像のうち所望の画像部分
を周波数強調する画像処理装置において、 Sproc=Sorg +f(x) (1) (ただし、f(x)は差信号に依存した関数を表す;x
=Sorg −Sus) 前記所望の画像部分に対応する前記差信号の絶対値が第
1の閾値と該第1の閾値より大きい第2の閾値との間の
範囲にあるときに、 前記差信号の絶対値が第1の閾値より小さい範囲におけ
る前記関数f(x)の値の絶対値が、前記所望の画像部
分に対応する前記関数f(x)の値より小さく、かつ前
記差信号の絶対値が第2の閾値より大きい範囲における
前記関数f(x)の値の絶対値が、前記所望の画像部分
に対応する前記関数f(x)の値より小さく設定されて
いることを特徴とする画像処理装置。
4. A non-sharp mask signal calculating means for obtaining a non-sharp mask signal Sus of an original image signal Sorg representing an image, and a difference signal calculating means for obtaining a difference signal between the original image signal Sorg and the non-sharp mask signal Sus. And a conversion table that converts the difference signal into a function f (x) that depends on the difference signal, and an emphasis processing unit that performs an emphasis process according to the following equation (1). In an image processing device for frequency-enhancing a part, Sproc = Sorg + f (x) (1) (where f (x) represents a function depending on the difference signal; x
= Sorg-Sus) When the absolute value of the difference signal corresponding to the desired image portion is in a range between a first threshold value and a second threshold value larger than the first threshold value, the difference signal The absolute value of the value of the function f (x) in the range where the absolute value is smaller than the first threshold value is smaller than the value of the function f (x) corresponding to the desired image portion, and the absolute value of the difference signal. Where the absolute value of the value of the function f (x) in the range where is larger than the second threshold value is set smaller than the value of the function f (x) corresponding to the desired image portion. Processing equipment.
【請求項5】 前記差信号を、前記画像のうち該差信号
に依存する高周波ノイズを表す補正用指標関数g(x)
に変換する補正用変換テーブルと、 下記式(2)にしたがって前記関数f(x)を補正せし
める補正手段とをさらに備えたことを特徴とする請求項
4記載の画像処理装置。 f′(x)=g(x)×f(x) (2) (ただし、g(x)は0≦g(x)≦1の範囲で単調に
増加する補正用指標関数、f′(x)は式(1)の関数
f(x)に代えて適用される補正後の関数を表す)
5. The correction index function g (x) representing the difference signal, which represents high-frequency noise in the image that depends on the difference signal.
5. The image processing apparatus according to claim 4, further comprising: a correction conversion table for converting into the following expression, and a correction unit that corrects the function f (x) according to the following expression (2). f ′ (x) = g (x) × f (x) (2) (where g (x) is a correction index function that monotonically increases in the range of 0 ≦ g (x) ≦ 1, f ′ (x ) Represents a corrected function that is applied in place of the function f (x) in equation (1))
【請求項6】 前記補正用指標関数g(x)は、下記式
(3)に示す前記差信号の分散値に応じた関数G(x)
であることを特徴とする請求項5記載の画像処理装置。 g(x)=G(Σ|Sorg −Sus|2 ) (3)
6. The correction index function g (x) is a function G (x) corresponding to the variance value of the difference signal shown in the following equation (3).
The image processing apparatus according to claim 5, wherein g (x) = G (Σ | Sorg-Sus | 2 ) (3)
JP7250857A 1995-09-28 1995-09-28 Image processing method and processor Pending JPH0991423A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP7250857A JPH0991423A (en) 1995-09-28 1995-09-28 Image processing method and processor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP7250857A JPH0991423A (en) 1995-09-28 1995-09-28 Image processing method and processor

Publications (1)

Publication Number Publication Date
JPH0991423A true JPH0991423A (en) 1997-04-04

Family

ID=17214046

Family Applications (1)

Application Number Title Priority Date Filing Date
JP7250857A Pending JPH0991423A (en) 1995-09-28 1995-09-28 Image processing method and processor

Country Status (1)

Country Link
JP (1) JPH0991423A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002525724A (en) * 1998-09-16 2002-08-13 インテル・コーポレーション General-purpose image enhancement framework
JP2004127064A (en) * 2002-10-04 2004-04-22 Konica Minolta Holdings Inc Image processing method, image processor, image processing program and image recording device
US7248748B2 (en) 2000-11-30 2007-07-24 Canon Kabushiki Kaisha Image processing apparatus, image processing method, storage medium, and program
US7558434B2 (en) 2000-11-30 2009-07-07 Canon Kabushiki Kaisha Image processing apparatus, image processing method, storage medium, and program
JP2012027547A (en) * 2010-07-20 2012-02-09 Mitsubishi Electric Corp Image processing apparatus

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002525724A (en) * 1998-09-16 2002-08-13 インテル・コーポレーション General-purpose image enhancement framework
US7248748B2 (en) 2000-11-30 2007-07-24 Canon Kabushiki Kaisha Image processing apparatus, image processing method, storage medium, and program
US7447376B2 (en) 2000-11-30 2008-11-04 Canon Kabushiki Kaisha Image processing apparatus, image processing method, storage medium, and program
US7558434B2 (en) 2000-11-30 2009-07-07 Canon Kabushiki Kaisha Image processing apparatus, image processing method, storage medium, and program
US7561750B2 (en) 2000-11-30 2009-07-14 Canon Kabushiki Kaisha Image processing apparatus, image processing method, storage medium, and program
JP2004127064A (en) * 2002-10-04 2004-04-22 Konica Minolta Holdings Inc Image processing method, image processor, image processing program and image recording device
JP2012027547A (en) * 2010-07-20 2012-02-09 Mitsubishi Electric Corp Image processing apparatus

Similar Documents

Publication Publication Date Title
US9245323B2 (en) Medical diagnostic device and method of improving image quality of medical diagnostic device
JP3816151B2 (en) Image processing method and apparatus
EP1550981B1 (en) System and method for image noise reduction using a minimal error spatio-temporal recursive filter
JP3690882B2 (en) Image enhancement processing method and apparatus
US5963676A (en) Multiscale adaptive system for enhancement of an image in X-ray angiography
US5150421A (en) System for automated transformation of gray level of image
EP1208536B1 (en) Contrast adjustment by computing a local mean density
US6813335B2 (en) Image processing apparatus, image processing system, image processing method, program, and storage medium
US6760401B2 (en) Apparatus and method for processing of digital images
JPH08329242A (en) Method for enhancement of image quality
JP3700804B2 (en) Image processing method and apparatus
US7149358B2 (en) Method and system for improving contrast using multi-resolution contrast based dynamic range management
JP3738788B2 (en) Image dynamic range compression processing method and apparatus
JP3738791B2 (en) Image processing method and apparatus
JP3700798B2 (en) Image processing method and apparatus
JP3973181B2 (en) Method for processing a series of radiological images of an object
JPH0991423A (en) Image processing method and processor
JP2001283215A (en) Image processor
JP2001344601A (en) Image processing apparatus and image processing program
JP3783116B2 (en) Radiation image enhancement processing method and apparatus
JPH0991421A (en) Image processing method and processor
JP2002125153A (en) Noise suppression processing unit and recording medium
JPH0998293A (en) Image processing method and device
JPH09107479A (en) Method and device for image processing
US6956977B2 (en) Methods for improving contrast based dynamic range management

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20041022

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20041102

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20051115

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20061128

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20061201