JP2000268171A - Picture processing method and picture processor and record medium - Google Patents

Picture processing method and picture processor and record medium

Info

Publication number
JP2000268171A
JP2000268171A JP11363655A JP36365599A JP2000268171A JP 2000268171 A JP2000268171 A JP 2000268171A JP 11363655 A JP11363655 A JP 11363655A JP 36365599 A JP36365599 A JP 36365599A JP 2000268171 A JP2000268171 A JP 2000268171A
Authority
JP
Japan
Prior art keywords
image signal
pixel density
conversion
conversion function
original image
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
JP11363655A
Other languages
Japanese (ja)
Inventor
Masahiko Yamada
雅彦 山田
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 JP11363655A priority Critical patent/JP2000268171A/en
Publication of JP2000268171A publication Critical patent/JP2000268171A/en
Pending legal-status Critical Current

Links

Landscapes

  • Image Processing (AREA)
  • Studio Circuits (AREA)

Abstract

PROBLEM TO BE SOLVED: To obtain the same frequency responding characteristics regardless of the pixel density of an original picture signal in a picture processing method for preparing plural picture signals each having different band limiting characteristics from an original picture signal, and for preparing plural converted picture signals by converting each picture signal based on plural conversion functions, and for obtaining an already processed picture signal to which processing such as the emphasis of the specific frequency components of the original picture signal is applied from each converted picture signal. SOLUTION: The pixel density information M of an original picture signal Sorg is inputted by a pixel density information inputting means 5, and the definition parameter of a conversion function in a conversion function defining means 4 is decided based on the pixel density information M. Then, processing is executed by a non-linear processing means 3 by using the function defined by the decided parameter.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、画像信号に対して
所定の周波数成分を強調するなどの処理を行うための画
像処理方法および装置並びに画像処理方法をコンピュー
タに実行させるためのプログラムを記録したコンピュー
タ読取り可能な記録媒体に関するものである。
The present invention relates to an image processing method and apparatus for performing processing such as emphasizing a predetermined frequency component on an image signal, and a program for causing a computer to execute the image processing method. The present invention relates to a computer-readable recording medium.

【0002】[0002]

【従来の技術】従来、本出願人により、非鮮鋭マスク画
像信号(以下、ボケ画像信号という)を用いて周波数強
調処理あるいはダイナミックレンジ圧縮処理などを行っ
て放射線画像の診断性能を向上させる数々の画像処理方
法および装置が提案されている(特開昭55-163472号、
同55-87953号、特開平3-222577号、特開平10-75395号、
同10-171983号など)。例えば周波数強調処理は、原画
像信号Sorgからボケ画像信号Susを引いたものに強調
係数βを乗じたものを、原画像信号Sorgに加算するこ
とにより、原画像信号の所定の空間周波数成分を強調す
るものである。これを式で表すと下記の式(1)のよう
になる。
2. Description of the Related Art Heretofore, a number of techniques have been proposed by the present applicant to improve the diagnostic performance of radiation images by performing frequency enhancement processing or dynamic range compression processing using an unsharp mask image signal (hereinafter referred to as a blurred image signal). An image processing method and apparatus have been proposed (JP-A-55-163472,
No. 55-87953, JP-A-3-222577, JP-A-10-75395,
No. 10-171983). For example, in the frequency emphasis processing, a predetermined spatial frequency component of the original image signal is enhanced by adding a value obtained by subtracting the blurred image signal Sus from the original image signal Sorg and an enhancement coefficient β to the original image signal Sorg. Is what you do. This is expressed by the following equation (1).

【0003】 Sproc=Sorg+β×(Sorg−Sus) …(1) (Sproc:周波数強調処理された信号、Sorg:原画像
信号、Sus:ボケ画像信号、β:強調係数) また、特開平10-75395号には、原画像信号に加算する加
算信号の周波数応答特性を調整することにより、周波数
強調処理された信号に対してアーチファクトが発生する
ことを防止する方法が提案されている。この方法とは、
まず鮮鋭度の異なる、すなわち周波数応答特性の異なる
複数のボケ画像信号を作成し、そのボケ画像信号および
原画像信号の中の2つの信号の差分をとることにより、
原画像信号の、ある限られた周波数帯域の周波数成分を
表す複数の帯域制限画像信号を作成し、さらにその帯域
制限画像信号をそれぞれ異なる変換関数によって所望の
大きさとなるように変換してから、その複数の抑制され
た帯域制限画像信号を積算することにより上記加算信号
を作成するものである。この処理は例えば下記の式
(2)により表すことができる。
Sproc = Sorg + β × (Sorg−Sus) (1) (Sproc: frequency-enhanced signal, Sorg: original image signal, Sus: blurred image signal, β: enhancement coefficient) Japanese Patent Application Laid-Open No. H11-163873 proposes a method of adjusting the frequency response characteristics of an addition signal to be added to an original image signal, thereby preventing the occurrence of an artifact in a signal subjected to frequency emphasis processing. This method is
First, a plurality of blurred image signals having different sharpness, that is, different frequency response characteristics are created, and a difference between the two signals in the blurred image signal and the original image signal is obtained.
From the original image signal, create a plurality of band-limited image signals representing frequency components of a limited frequency band, and further convert the band-limited image signals to have desired sizes by different conversion functions, The addition signal is created by integrating the plurality of suppressed band-limited image signals. This processing can be represented, for example, by the following equation (2).

【0004】 Sproc=Sorg+β(Sorg)×Fusm(Sorg,Sus1,Sus2,…SusN) Fusm(Sorg,Sus1,Sus2,…SusN) =f(Sorg−Sus1)+f(Sus1−Sus2)+… +f(Susk-1−Susk)+…+f(SusN-1−SusN) …(2 ) (但し、Sproc:処理済み画像信号 Sorg :原画像信号 Susk(k=1〜N):ボケ画像信号 f(k=1〜N):各帯域制限画像信号を変換する変換関
数 β(Sorg):原画像信号に基づいて定められる強調係
数) さらに、特開平10-171983号には、周波数強調処理とダ
イナミックレンジ圧縮処理とを同時に施す場合におい
て、処理が施された信号に対してアーチファクトが発生
することを防止する方法が提案されている。この方法と
は、上記特開平10-75395号に記載されたように、複数の
帯域制限画像信号を作成し、この帯域制限画像信号に基
づいて原画像信号の高周波成分に関する信号(高周波成
分信号)および低周波成分に関する信号(低周波成分信
号)を得、原画像信号にこれらの高周波成分に関する信
号および低周波成分に関する信号を加算することにより
周波数強調処理を施すとともにダイナミックレンジ圧縮
処理を施すようにしたものである。この処理は例えば下
記の式(3)により表すことができる。
Sproc = Sorg + β (Sorg) × Fusm (Sorg, Sus1, Sus2,... SusN) Fusm (Sorg, Sus1, Sus2,... SusN) = f 1 (Sorg−Sus1) + f 2 (Sus1−Sus2) +. k (Susk-1−Susk) +... + f N (SusN−1−SusN) (2) (where, Sproc: processed image signal Sorg: original image signal Susk (k = 1 to N): blurred image signal f k (k = 1 to N): a conversion function for converting each band-limited image signal β (Sorg): an enhancement coefficient determined based on the original image signal) Further, Japanese Patent Application Laid-Open No. 10-171983 discloses a frequency enhancement process. There has been proposed a method for preventing the occurrence of artifacts in a processed signal when the dynamic range compression processing is performed simultaneously. As described in Japanese Patent Application Laid-Open No. H10-75395, this method creates a plurality of band-limited image signals, and based on the band-limited image signals, signals related to high-frequency components of the original image signal (high-frequency component signals). And a signal relating to a low-frequency component (low-frequency component signal) is obtained, and a signal relating to the high-frequency component and a signal relating to the low-frequency component are added to the original image signal to perform the frequency emphasis processing and the dynamic range compression processing. It was done. This processing can be represented, for example, by the following equation (3).

【0005】 Sproc=Sorg +β(Sorg)・Fusm(Sorg,Sus1,Sus2,…SusN) +D(Sorg−Fdrc(Sorg,Sus1,Sus2,…SusN))…(3 ) Fusm(Sorg,Sus1,Sus2,…SusN) ={fu1(Sorg−Sus1)+fu2(Sus1 −Sus2)+… +fuk(Susk-1−Susk )+…+fuN(SusN-1−SusN)} Fdrc(Sorg,Sus1,Sus2,…SusN) ={fd1(Sorg−Sus1)+fd2(Sus1 −Sus2)+… +fdk(Susk-1−Susk )+…+fdN(SusN-1−SusN)} (但し、Sproc:処理済み画像信号 Sorg:原画像信号 Susk(k=1〜N):ボケ画像信号 fuk(k=1〜N):高周波成分信号を得るために使用す
る変換関数 fdk(k=1〜N):低周波成分信号を得るために使用す
る変換関数 β(Sorg):原画像信号に基づいて定められる強調係
数 D(Sorg−Fdrc):低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数(DはSorg−Fdrcを
変換する関数)) これらの周波数強調処理やダイナミックレンジ圧縮処理
(以下変換処理とする)においては、帯域制限画像信号
を変換する変換関数などの定義を変更することによって
原画像信号に加算する加算信号の周波数応答特性を調整
することができる。このため、各変換関数の定義次第
で、アーチファクトの発生防止など所望の周波数応答特
性を有する処理済み画像信号を得ることができる。しか
しながら、実際に変換関数を定義しようとした場合、各
変換関数をどのように定義すれば所望の結果が得られる
かは、上記(2)および(3)のような式から簡単に分
かるというものではない。このため、処理済み画像信号
の所望とする周波数応答特性を指定し、この指定された
特性に基づいて変換関数の定義パラメータを決定するこ
とにより、所望の周波数応答特性を有する処理済み画像
信号を簡易に得ることができる方法が提案されている
(特開平10-63838号)。
Sproc = Sorg + β (Sorg) · Fusm (Sorg, Sus1, Sus2,... SusN) + D (Sorg−Fdrc (Sorg, Sus1, Sus2,... SusN)) (3) Fusm (Sorg, Sus1, Sus2, ... SusN) = {f u1 ( Sorg-Sus1) + f u2 (Sus1 -Sus2) + ... + f uk (Susk-1-Susk) + ... + f uN (SusN-1-SusN)} Fdrc (Sorg, Sus1, Sus2, ... SusN) = {f d1 (Sorg−Sus1) + f d2 (Sus1−Sus2) +... + F dk (Susk−1−Susk) +... + F dN (SusN−1−SusN)} (where Sproc: processed image signal Sorg: original image signal Susk (k = 1~N): unsharp image signals f uk (k = 1~N): transformation function f dk (k = 1~N) used to obtain a high frequency component signal: low A conversion function β (Sorg) used to obtain a frequency component signal: an enhancement coefficient determined based on the original image signal D (Sorg−Fdrc): determined based on the low frequency component signal Dynamic range compression coefficient (D is a function for converting Sorg-Fdrc) In these frequency emphasizing processing and dynamic range compression processing (hereinafter referred to as conversion processing), a conversion function for converting a band-limited image signal is defined. By changing the frequency, the frequency response characteristic of the added signal to be added to the original image signal can be adjusted. For this reason, a processed image signal having desired frequency response characteristics such as prevention of generation of an artifact can be obtained depending on the definition of each conversion function. However, when actually trying to define a conversion function, it is easy to determine how to define each conversion function to obtain a desired result from the expressions such as the above (2) and (3). is not. For this reason, by specifying a desired frequency response characteristic of the processed image signal and determining the definition parameter of the conversion function based on the specified characteristic, the processed image signal having the desired frequency response characteristic can be simplified. (Japanese Patent Laid-Open No. 10-63838).

【0006】一方、上述した変換処理に用いられるボケ
画像信号は、原画像信号の画素に対して所定間隔毎に所
定のフィルタリング処理を施すことによって画素を間引
きし、その後間引きした数分の画素を所定の補間方法に
よって補間することにより作成される。このフィルタリ
ング処理としては、ローパスフィルタにより原画像信号
の高周波成分を取り除くような処理、具体的にはフィル
タ内の画素値の平均値や、荷重平均値を求めるといった
処理が行われている。上記特開平10-75395号などにおい
てボケ画像信号を得るために行われるフィルタリング処
理では、原画像信号に対してフィルタリング処理を施し
て得られた画素数の少ない信号に対してさらにこのフィ
ルタリング処理を施して、フィルタリングの各段階で得
られた画素数の少ない画像信号をそれぞれ原画像信号と
同じ画素数となるように補間することにより、複数の異
なるボケ画像信号を作成している。
On the other hand, the blurred image signal used in the above-described conversion processing is thinned out by subjecting the pixels of the original image signal to predetermined filtering processing at predetermined intervals, and then the pixels corresponding to the thinned number are obtained. It is created by interpolating according to a predetermined interpolation method. As the filtering process, a process of removing a high-frequency component of an original image signal by a low-pass filter, specifically, a process of obtaining an average value of pixel values in the filter and a weighted average value is performed. In the filtering process performed to obtain a blurred image signal in Japanese Patent Application Laid-Open No. 10-75395, the filtering process is further performed on a signal having a small number of pixels obtained by performing a filtering process on an original image signal. Then, a plurality of different blurred image signals are created by interpolating the image signals having a small number of pixels obtained at each stage of the filtering so as to have the same number of pixels as the original image signal.

【0007】ボケ画像信号は、上記のように原画像信号
に基づいて作成されるものであるが、この原画像信号
は、読取装置により原画像信号を所定の読取密度で読み
取って所定の画素密度を有する画像を再現可能なデジタ
ル信号としたものである。一般に、デジタル化された画
像信号を例えばプリント出力として再生する場合、画素
密度によって決定されるある周波数(ナイキスト周波
数)以下の周波数成分は正しく再生されることが知られ
ている。すなわち、読取密度すなわち画素密度は再生時
に必要となる画質レベルを考慮して決定されるため、必
ずしも一定であるとは限らない。
The blurred image signal is generated based on the original image signal as described above. The original image signal is read by a reading device at a predetermined reading density to obtain a predetermined pixel density. Is a digital signal that can be reproduced. In general, when a digitized image signal is reproduced as, for example, a print output, it is known that a frequency component equal to or lower than a certain frequency (Nyquist frequency) determined by the pixel density is reproduced correctly. That is, the reading density, that is, the pixel density is determined in consideration of the image quality level required at the time of reproduction, and is not always constant.

【0008】例えば、放射線画像読取再生システムで
は、蓄積性蛍光体シートに記録された人体の放射線画像
をレーザビーム走査によりデジタル画像信号として読み
取っているが、読取密度すなわち画素密度は蓄積性蛍光
体シートの大きさによって異なったり、ユーザの設定に
より任意の値に変更することができるものである。
For example, in a radiation image reading / reproducing system, a radiation image of a human body recorded on a stimulable phosphor sheet is read as a digital image signal by laser beam scanning. May be different depending on the size of the data, or can be changed to an arbitrary value by a user setting.

【0009】ここで、画素密度が異なる画像信号、すな
わちナイキスト周波数が異なる画像信号に対し、同一の
ローパスフィルタによるフィルタリング処理、および同
一の補間方法による補間処理を施した場合、得られる帯
域制限画像信号の周波数特性、具体的には帯域制限画像
信号の周波数帯域は画素密度によって異なることとな
る。これにより、例えば1つの原画像を2種類の読取密
度で読み取って2種類の画素密度を有する原画像信号を
得た場合、同じボケ画像信号を用いて帯域制限画像信号
を得て周波数強調処理あるいはダイナミックレンジ圧縮
処理を行っても、強調される周波数帯域あるいは圧縮さ
れる周波数帯域は2種類の原画像信号でそれぞれ異なっ
てしまうという問題が生ずる。
Here, when filtering processing using the same low-pass filter and interpolation processing using the same interpolation method are performed on image signals having different pixel densities, that is, image signals having different Nyquist frequencies, a band-limited image signal obtained is obtained. , Specifically, the frequency band of the band-limited image signal depends on the pixel density. Thus, for example, when one original image is read at two reading densities and an original image signal having two pixel densities is obtained, a band-limited image signal is obtained using the same blurred image signal, and a frequency emphasis process or Even if the dynamic range compression processing is performed, there is a problem that the emphasized frequency band or the compressed frequency band differs between the two types of original image signals.

【0010】このため、原画像信号の画素密度に関する
情報を得、その情報に基づいて複数のフィルタ係数リス
トからフィルタ係数を選択し、選択されたフィルタ係数
のフィルタによって原画像信号に対してフィルタリング
処理を行ってボケ画像信号を得るようにした方法が提案
されている(特開平10-63837号)。ここで例えば読取密
度が5本/mmと6.7 本/mmの原画像信号に対して同
一のローパスフィルタによりフィルタリング処理を施す
ことにより得られる帯域制限画像信号はその周波数帯域
が異なるものとなるが、この方法によれば各原画像信号
に対してそれぞれ異なるローパスフィルタ(実空間フィ
ルタ)によりフィルタリング処理を施すことにより、各
原画像信号から得られる各帯域制限画像信号の周波数帯
域のピーク位置を読取密度に依存することなく略一致さ
せることができ、このピーク位置が揃った帯域制限画像
信号を用いて周波数強調処理などを行うと、読取密度に
ある程度依存しない強調特性にすることができる。
For this reason, information on the pixel density of the original image signal is obtained, a filter coefficient is selected from a plurality of filter coefficient lists based on the information, and a filtering process is performed on the original image signal by a filter of the selected filter coefficient. Has been proposed to obtain a blurred image signal (Japanese Patent Laid-Open No. 10-63837). Here, for example, the band-limited image signals obtained by performing the filtering process with the same low-pass filter on the original image signals having the reading densities of 5 lines / mm and 6.7 lines / mm have different frequency bands. According to this method, each original image signal is subjected to a filtering process using a different low-pass filter (real space filter), so that the peak position of the frequency band of each band-limited image signal obtained from each original image signal can be read. Can be substantially matched without depending on the image density. If the frequency emphasizing process or the like is performed using the band-limited image signal having the aligned peak positions, the emphasizing characteristic can be made somewhat independent of the reading density.

【0011】[0011]

【発明が解決しようとする課題】しかしながら、上記特
開平10-63837号に記載された方法では、元々読取密度が
異なっており、帯域制限画像信号により表される帯域制
限画像が持っているエネルギ、すなわち帯域制限画像信
号の大きさは画素密度に応じて異なるため、上述のよう
に実空間フィルタによるフィルタリング処理を施して
も、各帯域制限画像信号の周波数帯域のピーク位置や形
状を完全に基準特性に一致させるということは困難(特
性上多少の差が生じる)で粗調のみができるに留まると
いう問題があり、各帯域制限画像信号に対して同じパラ
メータ(変換関数)を用いて特定の周波数帯域の帯域制
限画像信号を強調するような処理を施しても、画素密度
が異なると処理済画像信号の周波数応答特性が画素密度
に応じて微妙に異なり、同じ特性の画像にはならない。
However, according to the method described in Japanese Patent Application Laid-Open No. H10-63837, the reading density is originally different, and the energy of the band-limited image represented by the band-limited image signal has That is, since the size of the band-limited image signal differs according to the pixel density, even if the filtering process using the real space filter is performed as described above, the peak position and the shape of the frequency band of each band-limited image signal are completely determined by the reference characteristic. There is a problem that it is difficult to perform the rough adjustment only because it is difficult (there is a slight difference in characteristics), and the same parameter (conversion function) is used for each band-limited image signal. Even if a process that emphasizes the band-limited image signal is performed, if the pixel density is different, the frequency response characteristics of the processed image signal will differ slightly depending on the pixel density. Not be the image of the same characteristics.

【0012】本発明は上記事情に鑑みなされたものであ
り、原画像の画素密度に拘わらず所望の変換処理、例え
ば特定の周波数成分を強調するような処理を、常に同じ
ように行うことができる画像処理方法および装置並びに
画像処理方法をコンピュータに実行させるためのプログ
ラムを記録したコンピュータ読取り可能な記録媒体を提
供することを目的とするものである。
The present invention has been made in view of the above circumstances, and a desired conversion process, for example, a process for emphasizing a specific frequency component, can be always performed in the same manner regardless of the pixel density of an original image. An object of the present invention is to provide an image processing method and apparatus, and a computer-readable recording medium in which a program for causing a computer to execute the image processing method is recorded.

【0013】[0013]

【課題を解決するための手段】本発明による画像処理方
法は、所定の画素密度を有する原画像信号から、互いに
異なる帯域制限特性を有する複数の画像信号を作成し、
該各画像信号に対して所定の変換関数に基づいて変換処
理を施して変換画像信号を得、該変換画像信号から処理
済み画像信号を得る画像処理方法において、前記画素密
度に基づいて前記変換関数の定義パラメータを決定して
該変換関数を定義することを特徴とするものである。
According to the image processing method of the present invention, a plurality of image signals having different band limiting characteristics are created from an original image signal having a predetermined pixel density.
In an image processing method of performing a conversion process on each of the image signals based on a predetermined conversion function to obtain a converted image signal and obtaining a processed image signal from the converted image signal, the conversion function based on the pixel density Are defined, and the conversion function is defined.

【0014】ここで、「帯域制限特性」とは、原画像信
号における複数の周波数帯域毎の周波数応答特性のこと
をいう。
Here, the "band limiting characteristic" refers to a frequency response characteristic of each of a plurality of frequency bands in the original image signal.

【0015】また、「変換処理」とは、具体的には、上
記式(2)に示すような原画像信号に含まれる特定の周
波数成分を強調する周波数強調処理や、上記式(3)に
示すような原画像の最高濃度と最低濃度との差すなわち
ダイナミックレンジを狭めるように高濃度域もしくは低
濃度域あるいは高濃度域および低濃度域の双方のコント
ラストを下げるダイナミックレンジ圧縮処理などが挙げ
られる。
The "conversion processing" is specifically a frequency emphasis processing for emphasizing a specific frequency component included in the original image signal as shown in the above equation (2), or a frequency emphasis processing in the above equation (3). A dynamic range compression process that lowers the contrast in the high-density region or low-density region or in both the high-density region and the low-density region so as to narrow the difference between the maximum density and the minimum density of the original image, that is, the dynamic range, is exemplified. .

【0016】さらに、「所定の画素密度」は、オペレー
タが入力してもよく、原画像信号に画素密度に関する情
報を付与しておき、原画像信号を処理する際に画素密度
に関する情報を自動的に取得するようにしてもよい。ま
た、「画素密度」とは、例えば蓄積性蛍光体シートに記
録された放射線画像を読み取る際の読取密度のみなら
ず、例えば原画像のサイズと原画像信号を得る際のサン
プリングの関係を表す解像度の値(例えばdpi)であ
ってもよい。
Further, the "predetermined pixel density" may be input by an operator. Information relating to the pixel density is added to the original image signal, and the information relating to the pixel density is automatically processed when the original image signal is processed. May be obtained. The “pixel density” refers to not only the reading density at the time of reading a radiation image recorded on a stimulable phosphor sheet, for example, but also the resolution representing the relationship between the size of the original image and the sampling at the time of obtaining the original image signal. (For example, dpi).

【0017】さらにまた、「所定の画素密度に基づいて
前記変換関数の定義パラメータを決定してその変換関数
を定義する」とは、最終的に得られる処理済画像信号の
周波数応答特性が、画素密度の影響を受けないように、
画素密度に基づいて変換関数を定義することを意味す
る。定義パラメータを決定するに際しては、例えば、前
記複数の画像信号、所望の周波数応答特性、および変換
関数の定義パラメータの関係式を、複数の画像信号と周
波数応答特性を既知の値とし、変換関数の定義パラメー
タを変数とする連立方程式として解くことによりパラメ
ータを求めるなどの方法や、処理された画像を観察しな
がら、変換関数を定義するパラメータを少しずつ調整し
て、それぞれの変換関数を定義する方法を用いたり、あ
るいは、予め設定された少なくとも2以上の画素密度ご
とに定義パラメータを用意し、この用意した定義パラメ
ータのうち、予め設定された2以上の画素密度のうちの
所定の画素密度に最も近い画素密度に対応する定義パラ
メータを選択することにより変換関数の定義パラメータ
を決定して変換関数を定義する方法などを用いることも
できる。なお、変換関数としては、線形関数、非線型関
数、定数など種々の態様のものを用いることができる。
Further, "determining the definition parameter of the conversion function based on a predetermined pixel density and defining the conversion function" means that the frequency response characteristic of the finally obtained processed image signal is To be unaffected by density,
This means that a conversion function is defined based on the pixel density. When determining the definition parameters, for example, the relationship between the plurality of image signals, the desired frequency response characteristics, and the definition parameters of the conversion function, the plurality of image signals and the frequency response characteristics are known values, A method such as finding a parameter by solving it as a simultaneous equation with the definition parameter as a variable, or a method of defining each conversion function by gradually adjusting the parameters defining the conversion function while observing the processed image Alternatively, a definition parameter is prepared for each of at least two or more pixel densities that are set in advance. By selecting the definition parameter corresponding to the close pixel density, the definition parameter of the conversion function is determined and the conversion function is determined. And a method of definition may also be used. As the conversion function, various functions such as a linear function, a nonlinear function, and a constant can be used.

【0018】本発明による画像処理方法においては、決
定した定義パラメータを原画像信号に関連づけて保存す
ることが望ましい。
In the image processing method according to the present invention, it is desirable to store the determined definition parameters in association with the original image signal.

【0019】ここで「決定した定義パラメータを原画像
信号に関連づけて保存する」とは、画像を再出力させる
ときに、前回決定した定義パラメータを用いた処理済み
画像信号を用いることができるように、変換処理後に原
画像信号を保存する際に、決定した定義パラメータと原
画像とを関連付けて(対応付けて)保存することを意味
する。対応付けがとられている限りにおいては、定義パ
ラメータは、原画像信号と一緒に同じ記憶装置(記憶媒
体)に保存してもよいし、両者を別個の記憶装置に保存
してもよい。
Here, "storing the determined definition parameters in association with the original image signal" means that when the image is re-output, a processed image signal using the previously determined definition parameters can be used. This means that when the original image signal is stored after the conversion processing, the determined definition parameters and the original image are stored in association with each other. As long as the association is established, the definition parameters may be stored in the same storage device (storage medium) together with the original image signal, or both may be stored in separate storage devices.

【0020】なお、本発明による画像処理方法において
は、前記原画像信号に対して前記画素密度に基づいて定
められたフィルタ係数を有するフィルタによりフィルタ
リング処理を施して互いに周波数応答特性が異なる複数
のボケ画像信号を作成し、前記原画像信号および前記複
数のボケ画像信号に基づいて、前記原画像信号の複数の
周波数帯域毎の信号を表す複数の帯域制限画像信号を作
成し、この複数の帯域制限画像信号を複数の画像信号と
することが好ましい。
In the image processing method according to the present invention, the original image signal is subjected to a filtering process using a filter having a filter coefficient determined based on the pixel density, so that a plurality of blurs having different frequency response characteristics are obtained. Creating an image signal, creating a plurality of band-limited image signals representing signals for each of a plurality of frequency bands of the original image signal, based on the original image signal and the plurality of blurred image signals, Preferably, the image signal is a plurality of image signals.

【0021】ここで、「ボケ画像信号」とは、画素数は
原画像信号と同じであるが原画像信号よりも鮮鋭度が低
い画像を表す画像信号である。ボケ画像信号は、まず原
画像信号の画素に対して所定間隔毎に所定のフィルタリ
ング処理を施すことによって画素を間引きし、そのよう
にして得た画像信号に対して同様のフィルタリング処理
を繰り返してさらに画素数を少なくした画像信号を複数
作成し、そのそれぞれに対して、所定の補間方法により
原画像と画素数が同じになるように補間処理を施すこと
により作成する。ここで、フィルタリング処理は画素密
度に基づいて定められたフィルタ係数を有するフィルタ
により行われるものであり、その具体的手法については
上述した特開平10-63837号に記載されている。
Here, the "blurred image signal" is an image signal representing an image having the same number of pixels as the original image signal but lower in sharpness than the original image signal. The blurred image signal is first thinned out by performing a predetermined filtering process at predetermined intervals on the pixels of the original image signal, and the same filtering process is repeatedly performed on the image signal thus obtained. A plurality of image signals with a reduced number of pixels are created, and each of them is created by performing an interpolation process by a predetermined interpolation method so that the number of pixels is the same as that of the original image. Here, the filtering process is performed by a filter having a filter coefficient determined based on the pixel density, and a specific method is described in the above-mentioned Japanese Patent Application Laid-Open No. 10-63837.

【0022】また、「原画像信号の複数の周波数帯域毎
の信号を表す複数の帯域制限画像信号」は、例えば隣接
する周波数帯域のボケ画像信号同士で差分をとって作成
してもよいし、原画像信号と各ボケ画像信号の差分をと
って作成してもよい。あるいは原画像信号とボケ画像信
号の他の組み合わせで差分をとって作成することもでき
る。
The "plurality of band-limited image signals representing signals of a plurality of frequency bands of the original image signal" may be created by taking a difference between blurred image signals of adjacent frequency bands, for example. It may be created by taking the difference between the original image signal and each blurred image signal. Alternatively, it can be created by taking a difference with another combination of the original image signal and the blurred image signal.

【0023】本発明による画像処理装置は、所定の画素
密度を有する原画像信号から、互いに異なる帯域制限特
性を有する複数の画像信号を作成する画像信号作成手段
と、該各画像信号に対して所定の変換関数に基づいて変
換処理を施して変換画像信号を得、該変換画像信号から
処理済み画像信号を得る変換手段とを備えた画像処理装
置において、前記画素密度に基づいて前記変換関数の定
義パラメータを決定して該変換関数を定義する変換関数
定義手段を備えたことを特徴とするものである。
An image processing apparatus according to the present invention comprises: image signal creating means for creating a plurality of image signals having different band limiting characteristics from an original image signal having a predetermined pixel density; And a conversion means for performing a conversion process based on the conversion function to obtain a converted image signal and obtaining a processed image signal from the converted image signal, wherein the conversion function is defined based on the pixel density. A conversion function defining means for determining a parameter and defining the conversion function is provided.

【0024】本発明による画像処理装置においては、変
換関数定義手段を、予め設定された少なくとも2以上の
画素密度ごとに対応する定義パラメータを格納したテー
ブルを有し、このテーブルに格納された定義パラメータ
のうち、予め設定された2以上の画素密度のうちの所定
の画素密度に最も近い画素密度に対応する定義パラメー
タを読み出すことにより定義パラメータの決定を行なう
ものとすることができる。
In the image processing apparatus according to the present invention, the conversion function defining means has a table storing definition parameters corresponding to at least two or more predetermined pixel densities, and the definition parameter stored in this table is stored. The definition parameter may be determined by reading out the definition parameter corresponding to the pixel density closest to the predetermined pixel density among the two or more pixel densities set in advance.

【0025】本発明による画像処理装置においては、決
定した定義パラメータを原画像信号に関連づけて保存す
る記憶手段を備えたものとするのが望ましい。
It is desirable that the image processing apparatus according to the present invention includes a storage unit for storing the determined definition parameters in association with the original image signal.

【0026】なお、前記画像信号作成手段は、前記原画
像信号に対して前記画素密度に基づいて定められたフィ
ルタ係数を有するフィルタによりフィルタリング処理を
施して互いに周波数応答特性が異なる複数のボケ画像信
号を作成するボケ画像信号作成手段と、前記原画像信号
および前記複数のボケ画像信号に基づいて、前記原画像
信号の複数の周波数帯域毎の信号を表す複数の帯域制限
画像信号を前記複数の画像信号として作成する帯域制限
画像信号作成手段とからなることが好ましい。
The image signal creating means performs a filtering process on the original image signal with a filter having a filter coefficient determined based on the pixel density, and a plurality of blurred image signals having different frequency response characteristics from each other. A plurality of band-limited image signals representing signals of a plurality of frequency bands of the original image signal based on the original image signal and the plurality of blur image signals. It is preferable to include a band-limited image signal creating means for creating a signal.

【0027】また、本発明による画像処理方法をコンピ
ュータに実行させるためのプログラムとして、コンピュ
ータ読取り可能な記録媒体に記録して提供してもよい。
Further, a program for causing a computer to execute the image processing method according to the present invention may be provided by being recorded on a computer-readable recording medium.

【0028】[0028]

【発明の効果】本発明は、周波数強調処理などの変換処
理を行う場合に、原画像信号の画素密度に基づいて処理
済画像信号の周波数応答特性が画素密度の影響を受けな
いような、変換処理で用いられる変換関数の定義パラメ
ータを求め、そのパラメータにより定義される変換関数
に基づいて変換処理を行うようにしたものである。この
ため、読取密度が異なっても、互いに帯域制限特性が異
なるような複数の原画像信号に対する処理画像(強調処
理などが施された画像)の周波数応答特性を一致させる
ことができる、つまり、原画像信号の画素密度に影響さ
れない一定の周波数応答特性を有する画像を再現可能な
処理済み画像信号を得ることができる。
According to the present invention, when a conversion process such as a frequency emphasis process is performed, the frequency response characteristic of the processed image signal is not influenced by the pixel density based on the pixel density of the original image signal. A conversion function definition parameter used in the processing is obtained, and the conversion processing is performed based on the conversion function defined by the parameter. For this reason, even if the reading density differs, the frequency response characteristics of the processed images (images subjected to enhancement processing and the like) with respect to a plurality of original image signals having different band limiting characteristics can be matched. It is possible to obtain a processed image signal capable of reproducing an image having a constant frequency response characteristic that is not affected by the pixel density of the image signal.

【0029】また、予め設定された少なくとも2以上の
画素密度ごとに定義パラメータを用意し、この用意した
定義パラメータのうち、所定の画素密度に最も近い画素
密度に対応する画素密度の定義パラメータを選択するよ
うにすれば、変換関数の定義パラメータの決定を簡単に
行なうことができる。
Further, a definition parameter is prepared for each of at least two or more pixel densities set in advance, and a definition parameter of a pixel density corresponding to a pixel density closest to a predetermined pixel density is selected from the prepared definition parameters. By doing so, it is possible to easily determine the definition parameters of the conversion function.

【0030】また、決定した定義パラメータを原画像に
関連づけて保存するようにすれば、画像をフィルムやモ
ニターに再出力するとき、前回出力した画像に使用した
定義パラメータを用いて処理済み画像信号が得られるの
で、前回と同等画質の画像を得ることができる。
Further, if the determined definition parameters are stored in association with the original image, when the image is re-output to a film or a monitor, the processed image signal is processed using the definition parameters used for the previously output image. Therefore, an image having the same image quality as the previous image can be obtained.

【0031】[0031]

【発明の実施の形態】以下、本発明の画像処理方法およ
び装置の一実施形態について、図面を参照して詳細に説
明する。以下に示す画像処理装置は、蓄積性蛍光体シー
トに記録された人体の放射線画像を読み取って得た画像
信号に対して、その画像が診断に適した画像となるよう
に、ボケ画像信号を使用して周波数強調処理を施すもの
であり、処理された画像信号は主としてフィルムに記録
され、診断に用いられる。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Hereinafter, an embodiment of an image processing method and apparatus according to the present invention will be described in detail with reference to the drawings. The image processing apparatus described below uses a blurred image signal for an image signal obtained by reading a radiation image of a human body recorded on a stimulable phosphor sheet so that the image becomes an image suitable for diagnosis. Then, the processed image signal is mainly recorded on a film and used for diagnosis.

【0032】図1は本実施形態による画像処理装置の構
成を示す概略ブロック図である。画像処理装置1は、読
取装置などにおいて得られた所定の画素密度を有する原
画像信号Sorgからボケ画像信号を作成するボケ画像信
号作成手段2と、特定の周波数を強調するための周波数
強調処理を行って処理済み画像信号Sprocを得る非線形
処理手段3とを有する。さらに、この画像処理装置1
は、変換関数定義手段4と画素密度情報入力手段5とを
備える。変換関数定義手段4は、非線形処理手段3が変
換処理に使用する変換関数を定義する手段であり、例え
ば関数の傾き、非線形の程度などのパラメータを画素密
度情報入力手段5から入力された画素密度情報Mに基づ
いて決定して変換関数を定義する。画素密度情報入力手
段5は原画像信号Sorgの画素密度情報Mを得るための
手段である。画素密度情報入力手段5による入力は、ユ
ーザがキーボードから数値として入力してもよいし、操
作画面に数種類の密度を表示してユーザに選択させると
いった形態でもよい。あるいは、読取装置側で、画素密
度情報Mを原画像信号Sorgに添付しておき、画像処理
装置1が入力された画像信号毎に添付された情報Mを認
識するようにしてもよく、変換関数定義手段4が画素密
度を認識することができれば、どのような形態であって
もよい。
FIG. 1 is a schematic block diagram showing the configuration of the image processing apparatus according to the present embodiment. The image processing apparatus 1 includes a blurred image signal generating unit 2 that generates a blurred image signal from an original image signal Sorg having a predetermined pixel density obtained by a reading device and the like, and a frequency emphasis process for emphasizing a specific frequency. And a non-linear processing means 3 for obtaining a processed image signal Sproc. Further, the image processing apparatus 1
Comprises a conversion function definition means 4 and a pixel density information input means 5. The conversion function definition means 4 is a means for defining a conversion function used by the non-linear processing means 3 for the conversion processing. For example, parameters such as the slope of the function and the degree of non-linearity are input to the pixel density information input means 5. A conversion function is determined based on the information M. The pixel density information input means 5 is a means for obtaining the pixel density information M of the original image signal Sorg. The input by the pixel density information input unit 5 may be input by the user as a numerical value from a keyboard, or by displaying several types of densities on the operation screen and allowing the user to select. Alternatively, the reading device may attach the pixel density information M to the original image signal Sorg, and the image processing device 1 may recognize the information M attached to each input image signal. Any form may be used as long as the definition unit 4 can recognize the pixel density.

【0033】ここで、まずボケ画像信号の作成処理につ
いて詳細に説明する。図2はボケ画像信号作成処理の概
要を示すブロック図である。図2に示されるように、図
1のボケ画像信号作成手段2は、まずフィルタリング処
理手段10により、原画像信号Sorgに対し、原画像の
画素のx方向およびy方向に対してフィルタリング処理
を施して原画像信号Sorgよりも解像度が低い画像信号
(以下、低解像度画像信号という)を作成し、次に
この低解像度画像信号Bに対して同様のフィルタリン
グ処理を施してこの低解像度画像信号Bよりもさらに
解像度が低い低解像度画像信号Bを作成し、以降順次
同様のフィルタリング処理を重ねていくものである。そ
して、補間処理手段11により、このフィルタリング処
理の各段において得られる低解像度画像信号Bに対し
て、それぞれ補間拡大処理を施して、鮮鋭度の異なる複
数のボケ画像信号Sus1〜SusN(以下Susk(k=1〜
N)で代表させる)を得るものである。
Here, the process of creating a blurred image signal will be described in detail. FIG. 2 is a block diagram showing an outline of the blurred image signal creation processing. As shown in FIG. 2, the blurred image signal creating means 2 of FIG. 1 first performs filtering processing on the original image signal Sorg in the x direction and the y direction of the pixels of the original image by the filtering processing means 10. To generate an image signal B 1 (hereinafter referred to as a low-resolution image signal) having a lower resolution than the original image signal Sorg, and then perform the same filtering processing on the low-resolution image signal B 1 to perform the low-resolution image signal processing. create a low-resolution image signal B 2 further has lower resolution than the signal B 1, those to superimpose successive similar filtering process later. Then, the interpolation processing unit 11, the low-resolution image signal B k obtained at each stage of this filtering process, respectively subjected to interpolation enlargement processing, the sharpness of different blurred image signal SUS1 to SUSn (hereinafter Susk (K = 1 ~
N).

【0034】本実施形態においては、上記フィルタリン
グ処理のフィルタとして、1次元ガウス分布に略対応し
たフィルタを使用する。すなわちフィルタのフィルタ係
数を、ガウス信号に関する下記の式(4)にしたがって
定める。
In this embodiment, a filter substantially corresponding to a one-dimensional Gaussian distribution is used as a filter for the filtering process. That is, the filter coefficient of the filter is determined according to the following equation (4) regarding the Gaussian signal.

【0035】[0035]

【数1】 これは、ガウス信号は周波数空間および実空間の双方に
おいて局在性がよいためであり、例えば上記(4)式に
おいてσ=1とした場合の7×1の1次元フィルタは図
3(a)に示すようなものとなる。
(Equation 1) This is because the Gaussian signal has good localization in both the frequency space and the real space. For example, a 7 × 1 one-dimensional filter when σ = 1 in the above equation (4) is shown in FIG. It is as shown in

【0036】フィルタリング処理は、図4に示すよう
に、原画像信号Sorgに対して、あるいは低解像度画像
信号に対して1画素おきに行う。このような1画素おき
のフィルタリング処理をx方向、y方向に行うことによ
り、低解像度画像信号Bの画素数は原画像の1/4と
なり、フィルタリング処理により得られる低解像度画像
信号に対して繰り返しこのフィルタリング処理を施すこ
とにより、得られるn個の低解像度画像信号B(k=
1〜n)は、それぞれ、画素数が原画像信号Sorgの1
/22kの画像信号となる。
The filtering process is performed on the original image signal Sorg or every other pixel on the low-resolution image signal as shown in FIG. By performing such a filtering process for every other pixel in the x direction and the y direction, the number of pixels of the low resolution image signal B 1 becomes 4 of the original image, and the low resolution image signal obtained by the filtering process is reduced. By repeatedly performing this filtering process, n low-resolution image signals B k (k =
1 to n) each have a pixel number of 1 of the original image signal Sorg.
/ 22k image signal.

【0037】次に、このようにして得られた低解像度画
像信号Bに対して施される補間拡大処理について説明
する。補間演算の方法としては、Bスプラインによる方
法など種々の方法が挙げられるが、本実施形態において
は、上記フィルタリング処理においてガウス信号に基づ
くローパスフィルタを用いているため、補間演算につい
てもガウス信号を用いるものとする。具体的には、下記
の式(5)において、σ=2k−1 と近似したものを
用いる。
Next, a description will be given interpolation enlargement process performed on the low-resolution image signal B k obtained in this manner. As the method of the interpolation operation, various methods such as a method using a B-spline can be cited. In the present embodiment, a Gaussian signal is also used for the interpolation operation because a low-pass filter based on a Gaussian signal is used in the filtering processing. Shall be. Specifically, in equation (5) below, a value approximated to σ = 2 k−1 is used.

【0038】[0038]

【数2】 画像信号Bを補間する際には、k=1であるためσ=
1となる。この場合、補間処理を行うためのフィルタ
は、図5に示すように5×1の1次元フィルタとなる。
この補間処理は、まず低解像度画像信号Bに対して1
画素おきに値が0の画素を1つずつ補間することにより
低解像度画像信号Bを原画像と同一のサイズに拡大
し、次に、この補間された低解像度画像信号Bに対し
て上述した図5に示す1次元フィルタによりフィルタリ
ング処理を施すことにより行われる。
(Equation 2) When interpolating the image signals B 1 represents, for a k = 1 σ =
It becomes 1. In this case, the filter for performing the interpolation processing is a 5 × 1 one-dimensional filter as shown in FIG.
The interpolation process, first the low resolution image signal B 1 1
Expanding a low-resolution image signal B 1 in the original image and the same size by the value in every other pixel is interpolated by one pixel of 0, then above for the interpolated low resolution image signal B 1 The filtering is performed by the one-dimensional filter shown in FIG.

【0039】同様に、この補間拡大処理を全ての低解像
度画像信号Bに対して行う。低解像度画像信号B
補間する際には、上記式(5)に基づいて、3×2
1の長さのフィルタを作成し、画像信号Bの各画素の
間に値が0の画素を2−1個ずつ補間することによ
り、原画像と同一サイズに拡大し、この値が0の画素が
補間された画像信号Bに対して3×2−1の長さの
フィルタにより、フィルタリング処理を施すことにより
補間拡大する。
Similarly, this interpolation enlargement processing is performed on all the low-resolution image signals Bk . When interpolating the low-resolution image signal B k , 3 × 2 k
By creating a filter having a length of 1 and interpolating 2 k -1 pixels each having a value of 0 between each pixel of the image signal B k , the pixel is enlarged to the same size as the original image. Are subjected to a filtering process with a filter having a length of 3 × 2 k −1 on the image signal B k in which the pixels have been interpolated.

【0040】ここで、周波数強調処理では、原画像信号
Sorgが取得される際の画素密度は常に一定であること
が前提となっており、処理済み画像信号Sprocを作成す
るにあたり、この画素密度については何ら考慮がなされ
ていなかった。しかしながら、実際の画像処理システム
においては、種々の画像信号が入力されるため、全ての
画像信号の画素密度が同じであるとは限らない。例え
ば、本実施形態の放射線画像処理システムでは、この画
素密度は蓄積性蛍光体シートの大きさによって異なって
おり、半切/大角サイズで5本/mm、4つ切サイズで
6.7本/mm、6つ切サイズで10本/mmとなってい
る。また、この画素密度はユーザが任意に設定すること
も可能である。
Here, in the frequency emphasizing process, it is assumed that the pixel density at the time when the original image signal Sorg is obtained is always constant, and when creating the processed image signal Sproc, this pixel density Was not considered at all. However, in an actual image processing system, since various image signals are input, not all image signals have the same pixel density. For example, in the radiation image processing system according to the present embodiment, the pixel density differs depending on the size of the stimulable phosphor sheet, and is 5 / mm in half-cut / large-angle size, and in 4-cut size.
It is 6.7 lines / mm, and 10 lines / mm in a 6-cut size. Further, the pixel density can be arbitrarily set by the user.

【0041】しかしながら、画素密度が異なる画像信
号、すなわちナイキスト周波数が異なる画像信号に対
し、同一のローパスフィルタによるフィルタリング処
理、および同一の補間方法による補間処理を施した場
合、後述するようにしてボケ画像信号Suskから得られ
る帯域制限画像信号の周波数応答特性は画素密度によっ
て異なるものとなる。例えば、読取密度が10本/mm、
5本/mm、6.7本/mmのそれぞれの原画像信号に対
して、図3(a)に示す1次元フィルタを使用したフィ
ルタリング処理を施して得られたボケ画像信号Suskか
ら帯域制限画像信号を作成した場合、画素密度10本/m
mと5本/mmの場合においてはナイキスト周波数が1
/2の関係になっているため、その周波数帯域は等しく
なるものの、画素密度6.7本/mmの場合は画素密度10
本/mmと5本/mmの場合と全く異なるものとなる。
このため、本実施形態においては上記特開平10-63837号
に記載された方法のように、画素密度10本/mmおよび
5本/mmの場合には図3(a)に示すフィルタを、画
素密度6.7本/mmの場合には図3(b)に示すフィル
タを使用してフィルタリング処理を行い、いずれの画素
密度においても得られる帯域制限画像信号の周波数帯域
が略同一のものとなるようにしている。なお、ボケ画像
信号作成手段2におけるフィルタの切り替えは、画素密
度情報入力手段5から入力される画素密度情報Mに基づ
いて行われる。
However, when image signals having different pixel densities, that is, image signals having different Nyquist frequencies are subjected to the filtering process using the same low-pass filter and the interpolation process using the same interpolation method, the blurred image will be described later. The frequency response characteristics of the band-limited image signal obtained from the signal Susk differ depending on the pixel density. For example, the reading density is 10 lines / mm,
A band-limited image signal is obtained from a blur image signal Susk obtained by performing a filtering process using a one-dimensional filter shown in FIG. 3A on each of the original image signals of 5 lines / mm and 6.7 lines / mm. If created, pixel density 10 lines / m
In the case of m and 5 lines / mm, the Nyquist frequency is 1
/ 2, the frequency bands are equal, but when the pixel density is 6.7 lines / mm, the pixel density is 10
This is completely different from the case of 5 lines / mm and 5 lines / mm.
For this reason, in the present embodiment, as shown in the method described in the above-mentioned Japanese Patent Application Laid-Open No. 10-63837, when the pixel density is 10 lines / mm and 5 lines / mm, the filter shown in FIG. In the case of a density of 6.7 lines / mm, filtering is performed using the filter shown in FIG. 3B so that the frequency band of the band-limited image signal obtained at any pixel density is substantially the same. ing. Note that the filter switching in the blurred image signal creating means 2 is performed based on the pixel density information M input from the pixel density information input means 5.

【0042】次に、上記のようにして作成されたボケ画
像信号Suskを用いて行われる非線形処理について説明
する。図6は非線形処理の一例としての周波数強調処理
を行う装置の構成をボケ画像信号作成手段2とともに示
す概略ブロック図である。図6に示すように、原画像信
号Sorgと、フィルタリング処理手段10および補間処
理手段11により作成されたボケ画像信号Suskについ
て、減算器21によりそれらの信号の差分がとられ、原
画像信号Sorgの限られた周波数帯域の成分である帯域
制限画像信号(Sorg-Sus1、Sus1-Sus2など)が作成
される。
Next, the non-linear processing performed using the blurred image signal Susk generated as described above will be described. FIG. 6 is a schematic block diagram showing a configuration of an apparatus for performing a frequency emphasizing process as an example of the non-linear process together with the blurred image signal creating means 2. As shown in FIG. 6, a subtractor 21 calculates a difference between the original image signal Sorg and the blurred image signal Susk created by the filtering processing unit 10 and the interpolation processing unit 11, and calculates a difference between the signals. Band limited image signals (Sorg-Sus1, Sus1-Sus2, etc.), which are components of a limited frequency band, are created.

【0043】そして、この帯域制限画像信号が変換器2
2においてそれぞれ異なる変換関数f〜fにより所望
の大きさとなるように抑制され、さらに下記の式(2)
にしたがって、その複数の抑制された帯域制限画像信号
が演算器23において積算され、さらに原画像信号に加
算されて、処理済み画像信号Sprocが生成される。これ
により所望とする周波数成分を目的に応じた度合いで強
調した処理済み画像信号Sprocを得ることができる。
Then, the band-limited image signal is converted by the converter 2
2 are suppressed so as to have a desired size by different conversion functions f 1 to f N , respectively, and the following equation (2)
Accordingly, the plurality of suppressed band-limited image signals are integrated in the arithmetic unit 23 and further added to the original image signal to generate a processed image signal Sproc. As a result, it is possible to obtain a processed image signal Sproc in which a desired frequency component is emphasized to a degree according to the purpose.

【0044】 Sproc=Sorg +β(Sorg )×Fusm (Sorg,Sus1,Sus2,…SusN) Fusm (Sorg,Sus1,Sus2,…SusN) =f(Sorg −Sus1)+f(Sus1 −Sus2)+… +f(Susk-1−Susk )+…+f(SusN-1−SusN)…(2 ) (但し、Sproc:高周波成分が強調された画像信号 Sorg :原画像信号 Susk(k=1〜N):ボケ画像信号 f(k=1〜N):各帯域制限画像信号を変換する変換関
数 β(Sorg):原画像信号に基づいて定められる強調係
数) ここで、変換関数fとしては、下記の式(6)に示す
非線形関数が用いられる。
[0044] Sproc = Sorg + β (Sorg) × Fusm (Sorg, Sus1, Sus2, ... SusN) Fusm (Sorg, Sus1, Sus2, ... SusN) = f 1 (Sorg -Sus1) + f 2 (Sus1 -Sus2) + ... + F k (Susk-1−Susk) +... + F N (SusN−1−SusN)... (2) (where, Sproc: an image signal in which high-frequency components are emphasized Sorg: an original image signal Susk (k = 1 to N) : Blurred image signal f k (k = 1 to N): conversion function for converting each band-limited image signal β (Sorg): enhancement coefficient determined based on the original image signal) Here, the conversion function f k is: A non-linear function expressed by the following equation (6) is used.

【0045】 f(Sin)=Sout=Sin×Y×(exp(X/Sin)-1)/(exp(X/Sin)+1) …(6) 式(6)において、Sinは入力信号、Soutは出力信
号、Xは非線形の程度を決定するパラメータ、Yは関数
全体の傾きをコントロールするパラメータである。そし
て、パラメータX,Yを調整することにより帯域制限画
像信号の周波数応答特性を変更することができる。
F (Sin) = Sout = Sin × Y × (exp (X / Sin) −1) / (exp (X / Sin) +1) (6) In equation (6), Sin is an input signal, Sout is an output signal, X is a parameter for determining the degree of nonlinearity, and Y is a parameter for controlling the slope of the entire function. By adjusting the parameters X and Y, the frequency response characteristics of the band-limited image signal can be changed.

【0046】以上、処理済み画像信号Sprocの作成につ
いて説明したが、次に、本発明が解決しようとする問題
点およびその解決手段について例を示して説明する。上
記特開平10-63837号に記載された方法のように、画素密
度毎に異なるフィルタを用いて原画像信号Sorgに対し
てフィルタリングを行うことによりボケ画像信号Susk
を得、このボケ画像信号Suskから帯域制限画像信号を
得れば、画素密度に拘わらず帯域制限画像信号を同一の
周波数帯域のものとすることができる。しかしながら、
同一の周波数帯域とすることができるものの、帯域制限
画像信号の周波数応答特性は画素密度によって異なるも
のとなる。
The generation of the processed image signal Sproc has been described above. Next, problems to be solved by the present invention and means for solving the problems will be described with reference to examples. By filtering the original image signal Sorg using a different filter for each pixel density as in the method described in the above-mentioned Japanese Patent Application Laid-Open No. 10-63837, the blur image signal Susk
If the band-limited image signal is obtained from the blurred image signal Susk, the band-limited image signal can be in the same frequency band regardless of the pixel density. However,
Although the same frequency band can be set, the frequency response characteristics of the band-limited image signal differ depending on the pixel density.

【0047】ここで、読取密度が10本/mm、5本/m
m、6.7 本/mmのそれぞれの原画像信号に対して、図
3(a)、(b)に示す1次元フィルタを使用したフィ
ルタリング処理を施して作成した帯域制限画像信号の周
波数応答特性を、読取密度が10本/mmの場合を図7
に、6.7 本/mmの場合を図8に、5本/mmの場合を
図9にそれぞれ示す。図7から図9に示すように、帯域
制限画像信号のうち最高周波数帯域については画素密度
によって定められ、それ以外の周波数帯域の帯域制限画
像信号のピーク周波数(cycle/mm)は、各画素密度に
おいて、1(10本/mmのみ)、0.5、0.25、
0.125、0.06、0.03(6.7 本/mmおよび
5本/mmのみ)のように、k段目の帯域制限画像信号
のピーク周波数はk−1段目の帯域制限画像信号のピー
ク周波数の1/2になる。一方、各帯域制限画像信号
は、略同じピーク周波数であっても、各ピーク周波数で
のレスポンスの大きさ、すなわち帯域制限画像信号が有
するエネルギや、レスポンスの形状は、画素密度毎に若
干異なるものとなっている。このため、式(6)に示す
変換関数において画素密度に拘わらず同一のパラメータ
X,Yを使用したのでは、処理済み画像信号Sprocの周
波数応答特性が画素密度毎に若干異なるものとなってし
まう。
Here, the reading density is 10 lines / mm and 5 lines / m
The frequency response characteristics of the band-limited image signal created by performing the filtering process using the one-dimensional filter shown in FIGS. FIG. 7 shows the case where the reading density is 10 lines / mm.
FIG. 8 shows the case of 6.7 lines / mm, and FIG. 9 shows the case of 5 lines / mm. As shown in FIGS. 7 to 9, the highest frequency band of the band-limited image signal is determined by the pixel density, and the peak frequency (cycle / mm) of the band-limited image signal of the other frequency band is determined by each pixel density. In 1 (only 10 / mm), 0.5, 0.25,
As in 0.125, 0.06, and 0.03 (only 6.7 lines / mm and 5 lines / mm), the peak frequency of the band-limited image signal of the k-th stage is the peak frequency of the band-limited image signal of the (k-1) -th stage. It becomes 1/2 of the peak frequency. On the other hand, even if each band-limited image signal has substantially the same peak frequency, the magnitude of the response at each peak frequency, that is, the energy of the band-limited image signal and the shape of the response are slightly different for each pixel density. It has become. Therefore, if the same parameters X and Y are used regardless of the pixel density in the conversion function shown in Expression (6), the frequency response characteristics of the processed image signal Sproc will be slightly different for each pixel density. .

【0048】そこで、本発明では、画素密度情報入力手
段5から入力された原画像信号Sorgの画素密度に関す
る情報Mに基づいて、変換関数定義手段4において変換
関数のパラメータX,Yを設定して、画素密度毎に異な
るパラメータを有する変換関数により帯域制限画像信号
を変換するようにしたものである。変換関数定義手段4
の詳細な構成を図10に示す。図10に示すように、変
換関数定義手段4は、パラメータ設定手段7と画素密度
毎のパラメータX,Yの組を記憶したメモリ8とからな
る。
Therefore, in the present invention, the parameters X and Y of the conversion function are set in the conversion function definition means 4 based on the information M on the pixel density of the original image signal Sorg input from the pixel density information input means 5. The band-limited image signal is converted by a conversion function having a different parameter for each pixel density. Conversion function definition means 4
FIG. 10 shows a detailed configuration of the first embodiment. As shown in FIG. 10, the conversion function defining means 4 comprises a parameter setting means 7 and a memory 8 storing a set of parameters X and Y for each pixel density.

【0049】パラメータX,Yの組の例を下記の表1に
示す。表1に示すパラメータX,Yの組は10本/mmの
画素密度を基準とし、画素密度が6.7 本/mmおよび5
本/mmのパラメータX,Yにより処理済み画像信号S
procを得た際に、10本/mmのパラメータX,Yにより
得られた処理済み画像信号Sprocと略同じ周波数応答特
性が得られるように設定されている。なお、表1におい
ては、最上段がナイキスト周波数までの間にある帯域制
限画像信号に対するもの、最下段が0.06cycle/mm
にピークを持つ帯域制限画像信号に対するものであり、
上の行ほど高周波帯域の帯域制限画像信号に対するパラ
メータとなっている。また、本実施形態では、基準とす
る10本/mmの画素密度で読み取られたものにおいては
0.03cycle/mmにピークを有する帯域制限画像信号
が得られないため、6.7 本/mmおよび5本/mmの画
素密度で読み取られたものにおいて取得される0.03
cycle/mmにピークを有する帯域制限画像信号に対して
は強調処理が施されないように、パラメータの組が、最
高周波数帯域および0.5〜0.06までの5組に設定
されており、画素密度に拘わらず、0.06cycle/mm
にピークを持つ帯域制限画像信号からナイキスト周波数
までの間にある帯域制限画像信号を用いて強調処理が施
されることとなる。
An example of a set of parameters X and Y is shown in Table 1 below. The set of parameters X and Y shown in Table 1 is based on a pixel density of 10 lines / mm, and the pixel density is 6.7 lines / mm and 5 pixels / mm.
Image signal S processed by parameters X and Y
When the proc is obtained, the frequency response characteristics are set to be substantially the same as the processed image signal Sproc obtained by the parameters X and Y of 10 lines / mm. In Table 1, the uppermost row corresponds to the band-limited image signal up to the Nyquist frequency, and the lowermost row represents 0.06 cycle / mm.
For a band-limited image signal having a peak at
The upper row is a parameter for a band-limited image signal in a high frequency band. Further, in the present embodiment, since a band-limited image signal having a peak at 0.03 cycle / mm cannot be obtained for a pixel read at a pixel density of 10 pixels / mm as a reference, 6.7 pixels / mm and 5 pixels / 3 obtained at a pixel density of / 3
The set of parameters is set to the highest frequency band and five sets from 0.5 to 0.06 so that the enhancement processing is not performed on the band-limited image signal having a peak at cycle / mm. 0.06 cycle / mm regardless of density
The emphasis processing is performed using the band-limited image signal between the band-limited image signal having the peak at the Nyquist frequency and the band-limited image signal.

【0050】[0050]

【表1】 そして、画素密度情報入力手段5から入力された画素密
度情報Mに基づいてパラメータ選択手段7がメモリ8に
記憶されたパラメータの組を選択し、この選択されたパ
ラメータの組が変換関数定義手段4に入力されて変換関
数が定義され、さらに定義された変換関数に基づいて非
線形処理手段3において帯域制限画像信号に対して上記
式(2)に基づく周波数強調処理が行われる。
[Table 1] Then, the parameter selecting means 7 selects a set of parameters stored in the memory 8 based on the pixel density information M input from the pixel density information input means 5, and the selected set of parameters is stored in the conversion function defining means 4. To define a conversion function, and based on the defined conversion function, the nonlinear processing means 3 performs the frequency emphasis processing on the band-limited image signal based on the above equation (2).

【0051】上記パラメータX,Yにより設定された変
換関数による画素密度毎の周波数応答特性を図11から
図13に示す。図11から図13に示すように、表1に
示すパラメータX,Yにより設定された変換関数に基づ
いて周波数強調処理を行うと、画素密度に拘わらず周波
数応答特性が略同一のものとなる。
FIGS. 11 to 13 show frequency response characteristics for each pixel density based on the conversion function set by the parameters X and Y. As shown in FIGS. 11 to 13, when the frequency enhancement processing is performed based on the conversion functions set by the parameters X and Y shown in Table 1, the frequency response characteristics become substantially the same regardless of the pixel density.

【0052】なお、このパラメータX,Yは例えば上記
特開平10-63838号に記載された方法により設定すればよ
い。以下その設定の方法について説明する。なお、パラ
メータX,Yを設定するための装置としては、ディスプ
レイなどの表示装置を介してユーザに周波数応答特性の
入力を促して、その入力からユーザが所望する周波数応
答特性を認識するものであり、具体的には、そのような
処理を行うソフトウェア、入力装置、表示装置などを意
味する。
The parameters X and Y may be set, for example, by the method described in the above-mentioned Japanese Patent Application Laid-Open No. 10-63838. The setting method will be described below. An apparatus for setting the parameters X and Y prompts a user to input a frequency response characteristic via a display device such as a display, and recognizes a desired frequency response characteristic from the input. Specifically, it means software, an input device, a display device, and the like for performing such processing.

【0053】パラメータX,Yを設定する装置において
は、図14に示されるような操作画面がディスプレイに
表示される。このとき画面上には、その時点で設定され
ている変換関数(初期値)に基づいて帯域制限画像信号
を変換した場合に得られる信号の周波数応答特性曲線が
表示されている(図14の(1))。その特性曲線上に
は、移動可能な6個の指示点が表示されている。すなわ
ち、その画像処理がn個の帯域制限画像信号を作成して
非線形処理を行うものである場合、n個の指示点が表示
されるものとする。
In the apparatus for setting the parameters X and Y, an operation screen as shown in FIG. 14 is displayed on the display. At this time, a frequency response characteristic curve of a signal obtained when the band-limited image signal is converted based on the conversion function (initial value) set at that time is displayed on the screen ((FIG. 14) 1)). Six movable indication points are displayed on the characteristic curve. That is, if the image processing is to generate n band-limited image signals and perform nonlinear processing, it is assumed that n designated points are displayed.

【0054】ユーザは、各指示点を、例えばマウスなど
のポインティングデバイスを用いて所望の周波数応答に
対応する点に移動する(図14の(2))。この移動
は、画素密度に拘わらず各周波数帯域のレスポンスが同
一のものとなるように行われる。但し、表示されている
ままの周波数応答でよければ、必ずしも移動する必要は
ない。移動操作後の点(移動されなかった点を含む)
は、その点が示す周波数およびその周波数における所望
の周波数応答として認識され、これに基づいて周波数応
答特性曲線が計算しなおされ、移動操作後の点を通るよ
うな新しい特性曲線が画面上に表示される(図14の
(3))。ユーザは、確認画面において、新たに表示さ
れた特性曲線が所望のものであれば、パラメータ設定を
指示し、所望のものでなければ、再度指示点の移動操作
を行うことができる(図14の(4))。なお、パラメ
ータの設定はこのような形態に限られるものではなく、
例えばいくつかの周波数を順次画面に示し、ユーザがそ
の周波数に対する所望の周波数応答特性を数値として順
次入力するような形態など、他の様々な形態が考えられ
る。
The user moves each designated point to a point corresponding to a desired frequency response using a pointing device such as a mouse (FIG. 14 (2)). This movement is performed so that the response in each frequency band becomes the same regardless of the pixel density. However, if the frequency response as displayed is sufficient, it is not always necessary to move. Points after the move operation (including points that were not moved)
Is recognized as the frequency indicated by the point and the desired frequency response at that frequency, the frequency response characteristic curve is recalculated based on this, and a new characteristic curve passing through the point after the movement operation is displayed on the screen (3 in FIG. 14). If the newly displayed characteristic curve is a desired one on the confirmation screen, the user instructs parameter setting. If not, the user can perform the operation of moving the indicated point again (see FIG. 14). (4)). The parameter setting is not limited to such a form.
For example, various other forms are conceivable, such as a form in which several frequencies are sequentially displayed on a screen, and a user sequentially inputs desired frequency response characteristics for the frequencies as numerical values.

【0055】パラメータ設定の指示がなされると、例え
ば6つの帯域制限画像信号を用いて処理を行う場合に
は、入力された所望の周波数応答特性は、以下のような
連立方程式(7) a1 = Y1×S11+Y2×S21+Y3×S31+Y4×S41+Y5×S51+Y6×S61 a2 = Y1×S12+Y2×S22+Y3×S32+Y4×S42+Y5×S52+Y6×S62 a3 = Y1×S13+Y2×S23+Y3×S33+Y4×S43+Y5×S53+Y6×S63 a4 = Y1×S14+Y2×S24+Y3×S34+Y4×S44+Y5×S54+Y6×S64 a5 = Y1×S15+Y2×S25+Y3×S35+Y4×S45+Y5×S55+Y6×S65 a6 = Y1×S16+Y2×S26+Y3×S36+Y4×S46+Y5×S56+Y6×S66 Y1〜Y6 ≧ 0 …(7) (但し、a1〜a6 : 所望の値として指定された周波数応
答 S11〜S66 : 原画像信号に基づいて作成された帯域制限
画像信号の周波数応答特性 Y1〜Y6 : 変換関数の0付近の傾きを表すパラメー
タ) に、a1〜a6として当てはめられる。但し、方程式の次元
は、帯域制限画像信号の数、すなわち指定された値の数
に基づいて決まるものであり、上記連立方程式は一例に
すぎない。
When the parameter setting is instructed, for example, when processing is performed using six band-limited image signals, the input desired frequency response characteristic is calculated by the following simultaneous equation (7) a1 = Y1 x S11 + Y2 x S21 + Y3 x S31 + Y4 x S41 + Y5 x S51 + Y6 x S61 a2 = Y1 x S12 + Y2 x S22 + Y3 x S32 + Y4 x S42 + Y5 x S52 + Y6 x S62 a3 = Y1 x S13 + Y2 x S43 + Y3 x S3 + Y3 x S3 + Y3 x S4 + Y3 x S4 + Y3 x S4 + Y3 x S4 + Y3 x S4 + Y3 x S4 + Y3 x S4 + Y3 x S3 + Y3 x S4 + Y3 x S4 + Y3 x S4 + Y3 x S4 .Times.S34 + Y4.times.S44 + Y5.times.S54 + Y6.times.S64 a5 = Y1.times.S15 + Y2.times.S25 + Y3.times.S35 + Y4.times.S45 + Y5.times.S55 + Y6.times.S65 a6 = Y1.times.S16 + Y2.times.S26 + Y3.times.S36 + Y4.times.S46 + Y5.times.S6 + Y6.times.S6 + Y6. a1 to a6: Frequency response specified as desired value S11 to S66: Frequency response characteristic of band-limited image signal created based on original image signal Y1 to Y6: Parameter representing slope of conversion function near 0) , A1 ~ Fits as a6. However, the dimension of the equation is determined based on the number of band-limited image signals, that is, the number of designated values, and the above simultaneous equations are merely examples.

【0056】そして、式(7)に示す連立方程式を、指
定された周波数応答特性a1〜a6および帯域制限画像信号
S11〜S66を既定値とし、パラメータY1〜Y6を変数とし
て解くことによりこのパラメータY1〜Y6を求め、こ
れにより変換関数を定義するパラメータを設定する。な
お、非線形の程度を表すパラメータXは、パラメータY
1〜Y6に対応させて任意の値に設定すればよいが、表
1に示すように画素密度毎に一定の値を用いればよい。
そして、画素密度毎に設定されたパラメータX,Yの組
がメモリ8に記憶されることとなる。
Then, the simultaneous equations shown in the equation (7) are converted into the designated frequency response characteristics a1 to a6 and the band-limited image signal.
The parameters Y1 to Y6 are obtained by solving S11 to S66 as default values and the parameters Y1 to Y6 as variables, thereby setting the parameters defining the conversion function. Note that the parameter X representing the degree of nonlinearity is a parameter Y
Any value may be set corresponding to 1 to Y6, but a constant value may be used for each pixel density as shown in Table 1.
Then, a set of parameters X and Y set for each pixel density is stored in the memory 8.

【0057】また、上記式(7)において、帯域制限画
像信号の周波数応答特性S11〜S66として表されている値
は、S11〜S16(S21〜S26,…S61〜S66も同様)が1つの
帯域制限画像信号の所定の6つの周波数の応答特性に対
応するものであり、これらはボケ画像信号作成において
用いられるフィルタのフィルタ係数と、補間係数から、
フーリエ変換など、周波数解析において通常用いられて
いる変換方法により求めることができる。
In the above equation (7), the values expressed as the frequency response characteristics S11 to S66 of the band-limited image signal are such that S11 to S16 (S21 to S26,. These correspond to the response characteristics of the predetermined six frequencies of the restricted image signal, and these correspond to the filter coefficient of the filter used in the blurred image signal creation and the interpolation coefficient.
It can be obtained by a conversion method generally used in frequency analysis, such as Fourier transform.

【0058】次いで、本実施形態の動作について説明す
る。図15は本実施形態の動作を示すフローチャートで
ある。まず、読取装置などから原画像信号Sorgが画像
処理装置1に入力される(ステップS1)。原画像信号
Sorgはボケ画像信号作成手段2に入力されてここで原
画像信号Sorgの画素密度に応じてボケ画像信号Suskが
作成される(ステップS2)。一方、画素密度情報入力
手段5には、原画像信号Sorgの画素密度情報Mが入力
され(ステップS3)、入力された画素密度情報Mに基
づいて変換関数定義手段4において上述したように画素
密度に応じたパラメータの組が選択され、選択されたパ
ラメータの組により変換関数が定義される(ステップS
4)。なお、ステップS1,2の処理よりもステップS
3,4の処理を先に行ってもよいが、これらの処理を並
列して行うことにより演算時間を短縮することができ
る。非線形処理手段3においては、ボケ画像信号Susk
に基づいて帯域制限画像信号が作成され、さらに変換関
数定義手段4において定義された変換関数に基づいて上
記式(2)に示す非線形処理(周波数強調処理)が行わ
れ処理済み画像信号Sprocが得られる(ステップS
5)。
Next, the operation of this embodiment will be described. FIG. 15 is a flowchart showing the operation of the present embodiment. First, an original image signal Sorg is input to the image processing device 1 from a reading device or the like (step S1). The original image signal Sorg is input to the blurred image signal creating means 2, where the blurred image signal Susk is created according to the pixel density of the original image signal Sorg (step S2). On the other hand, pixel density information M of the original image signal Sorg is input to the pixel density information input means 5 (step S3), and based on the input pixel density information M, the pixel density information Is selected, and a conversion function is defined by the selected set of parameters (step S
4). It should be noted that the processing in step S1 is more than the processing in step S1
Although the processes 3 and 4 may be performed first, the processing time can be reduced by performing these processes in parallel. In the nonlinear processing means 3, the blurred image signal Susk
, A non-linear processing (frequency emphasis processing) represented by the above equation (2) is performed based on the conversion function defined by the conversion function definition means 4, and a processed image signal Sproc is obtained. (Step S
5).

【0059】この周波数強調処理は、ボケ画像信号を用
いて、強調のために原画像信号に加算する加算信号を作
成するものであり、周波数強調処理を行うことによって
アーチファクトが発生しないようにするために、この加
算信号を構成する各周波数帯域の信号がそれぞれ所望の
信号となるよう、周波数帯域毎に異なる変換処理を施す
ものである。所望の信号を作成するためには、原画像信
号Sorgの画素密度に拘わらず所望の周波数応答特性を
有する信号となることが好ましい。本実施形態によれ
ば、画素密度に拘わらず処理済み画像信号Sprocの周波
数応答特性が略同一のものとなるように変換関数のパラ
メータを定義するようにしたため、原画像信号Sorgの
画素密度に拘わらず略同一の周波数応答特性を有する処
理済み画像信号Sprocを得ることができる。
This frequency emphasis processing is to generate an addition signal to be added to the original image signal for emphasis using the blurred image signal. In order to prevent an artifact from occurring by performing the frequency emphasis processing. Then, different conversion processing is performed for each frequency band so that the signals of each frequency band constituting the addition signal become desired signals. In order to create a desired signal, the signal preferably has a desired frequency response characteristic regardless of the pixel density of the original image signal Sorg. According to the present embodiment, since the parameters of the conversion function are defined so that the frequency response characteristics of the processed image signal Sproc are substantially the same regardless of the pixel density, regardless of the pixel density of the original image signal Sorg. First, a processed image signal Sproc having substantially the same frequency response characteristics can be obtained.

【0060】上記説明では、パラメータX,Yを特開平
10-63838号 に記載の方法により設定するものとして説
明したが、これに限らず、試行錯誤により設定してもよ
い。
In the above description, the parameters X and Y are
Although the method is described as being set according to the method described in No. 10-63838, the present invention is not limited to this, and may be set by trial and error.

【0061】また、上記表1に示すような基本解像度
(a,b,c;a<b<c:cが最高解像度、上記例で
は、a=5,b=6.7 ,c=10)の定義パラメータのテ
ーブルを用意しておき、処理対象となる原画像解像度R
sが基本解像度(a,b,c;a<b<c:cが最高解
像度)と異なる場合には、基本解像度a,b,cの内、
最も近い基本解像度の定義パラメータを以下のルールに
したがって選択してもよい。
The definition of the basic resolution (a, b, c; a <b <c: c is the highest resolution, a = 5, b = 6.7, c = 10 in the above example) as shown in Table 1 above. A parameter table is prepared, and the resolution R of the original image to be processed
When s is different from the basic resolution (a, b, c; a <b <c: c is the highest resolution), among the basic resolutions a, b, c,
The definition parameter of the closest basic resolution may be selected according to the following rules.

【0062】 Rs>b+(c−b)/2のとき、c用のパラメータ b+(c−b)/2≧Rs>a+(b−a)/2のと
き、b用のパラメータ b+(c−b)/2≧Rsのときa用のパラメータ この方法を用いると、変換関数の定義パラメータの決定
を簡単に行なうことができる。
When Rs> b + (c−b) / 2, the parameter for c b + (c−b) / 2 ≧ Rs> a + (ba) / 2, the parameter for b b + (c− b) Parameter for a when / 2 ≧ Rs This method makes it easy to determine the definition parameter of the conversion function.

【0063】なお、記憶装置を設けて、処理が終了した
後に原画像信号Sorg を該記憶装置に保存することもで
きるが、この場合には、定義パラメータの決定方法の種
類を問わず、決定した定義パラメータ並びに基本解像度
情報および画像解像度情報を原画像信号Sorg の付帯情
報として、原画像に関連づけて該記憶装置に保存すると
よい。この付帯情報の保存先は、例えば原画像信号Sor
g を保存するものと同じ記憶装置であってもよいし、別
の記憶装置であってもよい。これにより、画像をフィル
ムやモニターに再出力するとき、前回出力した画像の原
画像信号Sorgに対応する定義パラメータを記憶装置か
ら読み出し、この読み出した定義パラメータを使用して
処理済み画像信号procを生成して、前回出力した画像と
同等画質の画像を得ることができる。また、定義パラメ
ータを決定する処理を省略できるので、短時間で画像が
出力される。
It should be noted that a storage device may be provided to store the original image signal Sorg after the processing is completed. In this case, however, the definition parameter is determined regardless of the method of determining the definition parameter. The definition parameter, the basic resolution information and the image resolution information may be stored in the storage device in association with the original image as supplementary information of the original image signal Sorg. The storage destination of the supplementary information is, for example, the original image signal Sor.
It may be the same storage device that stores g, or may be another storage device. Thereby, when the image is re-output to the film or the monitor, the definition parameter corresponding to the original image signal Sorg of the previously output image is read from the storage device, and the processed image signal proc is generated using the read definition parameter. Thus, an image having the same image quality as the previously output image can be obtained. Further, since the process of determining the definition parameters can be omitted, an image is output in a short time.

【0064】なお、上記実施形態においては、非線形処
理として上記式(2)に示す周波数強調処理を行ってい
るが、周波数強調処理とダイナミックレンジ圧縮処理と
を同時に行う処理であってもよい。以下この処理を行う
装置について説明する。図16は、非線形処理の一例と
しての周波数強調処理およびダイナミックレンジ圧縮処
理を行う装置の構成をボケ画像信号作成手段2とともに
示す概略ブロック図である。図16に示すように、原画
像信号Sorgと、フィルタリング処理手段10および補
間拡大手段11により作成されたボケ画像信号Suskに
ついて、減算器21によりそれらの信号の差分がとら
れ、原画像信号の、限られた周波数帯域の成分である帯
域制限画像信号(Sorg-Sus1、Sus1-Sus2など)が作
成される。ここで、フィルタリング処理手段10におい
て使用されるフィルタは上記実施形態と同様に画素密度
に応じて定められる。このようにして求められた帯域制
限画像信号は、図16に示すように第1の変換手段3a
および第2の変換手段3bにそれぞれ入力され、各変換
手段の変換器22a,22bにより処理される。
In the above embodiment, the frequency emphasis processing shown in the above equation (2) is performed as the non-linear processing. However, the frequency emphasis processing and the dynamic range compression processing may be performed simultaneously. Hereinafter, an apparatus for performing this processing will be described. FIG. 16 is a schematic block diagram showing a configuration of an apparatus for performing a frequency emphasis process and a dynamic range compression process as an example of the non-linear process together with the blurred image signal creating means 2. As shown in FIG. 16, a difference between the original image signal Sorg and the blurred image signal Susk created by the filtering processing unit 10 and the interpolation and enlarging unit 11 is calculated by a subtracter 21, and Band limited image signals (Sorg-Sus1, Sus1-Sus2, etc.), which are components of a limited frequency band, are created. Here, the filter used in the filtering processing means 10 is determined according to the pixel density as in the above embodiment. The band-limited image signal obtained in this manner is, as shown in FIG.
And the second conversion means 3b, respectively, and are processed by the converters 22a and 22b of each conversion means.

【0065】第1の変換手段3aの変換器22aによる
変換は、上述したように原画像信号Sorgの画素密度情
報Mに基づいて定義された変換関数を使用して行う。こ
こで、変換器22aにおいて変換は、例えば図17、図
18に示される変換関数、あるいはこれらを組み合わせ
た関数を使用して行う。なお、これらの変換関数は、原
画像信号Sorgの画素密度に拘わらず処理済み画像信号
Sprocの周波数応答特性が略同一のものとなるように、
画素密度毎にパラメータが設定されている。
The conversion by the converter 22a of the first conversion means 3a is performed using the conversion function defined based on the pixel density information M of the original image signal Sorg as described above. Here, the conversion in the converter 22a is performed using, for example, the conversion functions shown in FIGS. Note that these conversion functions are set so that the frequency response characteristics of the processed image signal Sproc are substantially the same regardless of the pixel density of the original image signal Sorg.
Parameters are set for each pixel density.

【0066】ここで、図17に示す変換関数は、振幅の
大きな帯域制限画像信号を抑制するような変換を行うも
のであり、周波数帯域の高い帯域制限画像信号の抑制の
度合いを、周波数帯域の低い帯域制限画像信号よりも強
くするものであるが、これは実際の放射線画像のエッジ
に含まれている高周波成分が、低周波成分に比べてその
振幅が小さいということを考慮したものである。実際の
放射線画像においては、かなり急峻なエッジでさえも正
確な階段状にはなっておらず、高周波成分になるほどそ
の振幅が小さくなっていることが多い。このため、各周
波数成分の振幅に合わせて、周波数の高い帯域制限画像
信号ほど小さい振幅から抑制を行うことが望ましく、本
関数によりそれを実現することができる。
Here, the conversion function shown in FIG. 17 performs a conversion to suppress a band-limited image signal having a large amplitude, and determines the degree of suppression of the band-limited image signal having a high frequency band. This is made stronger than the low band-limited image signal in consideration of the fact that the amplitude of the high-frequency component contained in the edge of the actual radiation image is smaller than that of the low-frequency component. In an actual radiographic image, even a very steep edge does not have an accurate step shape, and the higher the frequency component, the smaller its amplitude in many cases. For this reason, it is desirable to suppress the band-limited image signal with a higher frequency from the smaller amplitude in accordance with the amplitude of each frequency component, and this can be realized by this function.

【0067】また図18の関数は、帯域制限画像信号
を、帯域制限画像信号の絶対値に基づいて決まる、その
絶対値以下の値となるように変換を行うもので、この関
数が低周波帯域を処理する関数であるほど、帯域制限画
像信号の絶対値が0近傍の所定の範囲内の値である帯域
制限画像信号を変換した際に得られる変換画像信号の絶
対値が小さい値であることを特徴とするものである。言
い換えれば、これらの関数はそれぞれ、原点を通り、関
数の傾きがその関数により処理される値に拘わらず1以
下であり、その関数の0近傍における傾きが、低周波帯
域を処理する関数であるほど小さいことを特徴とするも
のである。これらの関数は、変換画像信号を積算して得
られる信号を、原画像信号Sorgに加えた場合、原画像
信号Sorgと加算された信号とのつなぎ目、すなわち信
号の立ち上がりをより自然なものとするという効果があ
る。
The function of FIG. 18 converts a band-limited image signal to a value that is determined based on the absolute value of the band-limited image signal and is equal to or less than the absolute value. The smaller the absolute value of the converted image signal obtained when the band-limited image signal whose absolute value is within a predetermined range near 0 is converted, the smaller the value of the function is It is characterized by the following. In other words, each of these functions passes through the origin and the slope of the function is less than or equal to 1 irrespective of the value processed by the function, and the slope of the function near 0 is a function for processing the low frequency band. It is characterized by being smaller. When these functions add a signal obtained by integrating the converted image signal to the original image signal Sorg, the joint between the original image signal Sorg and the added signal, that is, the rise of the signal becomes more natural. This has the effect.

【0068】同様に、第2の変換手段3bの変換器22
bによる変換は、例えば図19に示される関数、前述の
図18の関数、あるいは図19と図18の関数を組み合
わせたものを使用して行う。上記各変換手段3a,3b
により出力された変換画像信号は、それぞれ演算器23
a,23bに入力される。ここで演算器23aは周波数
強調処理に必要な信号を作成するための演算を行い、演
算器23bはダイナミックレンジ圧縮処理に必要な信号
を作成するための演算を行うものである。
Similarly, the converter 22 of the second converter 3b
The conversion by b is performed using, for example, the function shown in FIG. 19, the function of FIG. 18 described above, or a combination of the functions of FIG. 19 and FIG. Each of the conversion means 3a, 3b
The converted image signals output by
a and 23b. Here, the arithmetic unit 23a performs an operation for generating a signal necessary for the frequency emphasizing process, and the arithmetic unit 23b performs an operation for generating a signal necessary for the dynamic range compression process.

【0069】演算器23aは、上述した図6に示す実施
形態と同様の周波数強調処理を行うものである。この演
算器23aにおいては次のような処理が行われる。ま
ず、上述したように第1の変換手段3aにより変換され
た帯域制限画像信号は積算される。そしてこの積算信号
が求められると、さらに原画像信号Sorg の値に応じた
強調度βが乗じられる。
The arithmetic unit 23a performs the same frequency emphasizing processing as in the embodiment shown in FIG. The following processing is performed in the arithmetic unit 23a. First, as described above, the band-limited image signals converted by the first conversion unit 3a are integrated. When the integrated signal is obtained, the integrated signal is further multiplied by the enhancement degree β according to the value of the original image signal Sorg.

【0070】一方、演算器23bは、ダイナミックレン
ジ圧縮処理を行うものである。この演算器23bにおい
ては次のような処理が行われる。まず、上述したように
第2の変換手段3bにより変換された帯域制限画像信号
は積算される。そして求められた積算信号は原画像信号
Sorgから減じられる。さらにその減算により作成され
た差分信号が変換関数により変換され、ダイナミックレ
ンジ圧縮係数が得られる。
On the other hand, the arithmetic unit 23b performs a dynamic range compression process. The following processing is performed in the arithmetic unit 23b. First, as described above, the band-limited image signals converted by the second conversion unit 3b are integrated. Then, the obtained integrated signal is subtracted from the original image signal Sorg. Further, the difference signal created by the subtraction is converted by a conversion function to obtain a dynamic range compression coefficient.

【0071】演算器23a,23bにより求められた信
号は、それぞれ加算手段28により原画像信号Sorgに
加算され処理済み画像信号Sprocが得られる。
The signals obtained by the arithmetic units 23a and 23b are respectively added to the original image signal Sorg by the adding means 28 to obtain a processed image signal Sproc.

【0072】以上の処理は、下記の式(3)により表さ
れる。
The above processing is represented by the following equation (3).

【0073】 Sproc=Sorg +β(Sorg)・Fusm(Sorg,Sus1,Sus2,…SusN) +D(Sorg−Fdrc(Sorg,Sus1,Sus2,…SusN))…(3 ) Fusm(Sorg,Sus1,Sus2,…SusN) ={fu1(Sorg−Sus1)+fu2(Sus1 −Sus2)+… +fuk(Susk-1−Susk )+…+fuN(SusN-1−SusN)} Fdrc(Sorg,Sus1,Sus2,…SusN) ={fd1(Sorg−Sus1)+fd2(Sus1 −Sus2)+… +fdk(Susk-1−Susk )+…+fdN(SusN-1−SusN)} (但し、Sproc:処理済み画像信号 Sorg :原画像信号 Susk(k=1〜N):ボケ画像信号 fuk(k=1〜N):第1の変換処理に使用する変換関数 fdk(k=1〜N):第2の変換処理に使用する変換関数 β(Sorg):原画像信号に基づいて定められる強調係
数 D(Sorg−Fdrc):低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数(DはSorg−Fdrcを
変換する関数)) ここで、図16に示す実施形態においては、周波数強調
処理に使用する信号とダイナミックレンジ圧縮処理に使
用する信号をそれぞれ原画像信号Sorgに基づいて作成
して最後に加算を行うものであるが、原画像信号Sorg
に対して周波数強調処理あるいはダイナミックレンジ圧
縮処理のいずれか一方の処理を先に行い、その結果得ら
れた信号に対して他の一方の処理を行うようにする場合
も含むものとする。但し、一般に放射線画像の低濃度部
は撮影時の撮影線量が少ない部分であるためノイズが比
較的多く含まれている。これに対し、周波数強調処理で
は高濃度部ほど強調の度合いが高くなる濃度依存の強調
処理を行っているため、原画像信号に対して直接周波数
強調処理を行えばこの低濃度部は強調されず、ノイズが
強調されることがない。一方、例えば原画像信号に対し
ダイナミックレンジ圧縮処理を施してから周波数強調処
理を行うと、低濃度部はダイナミックレンジ圧縮処理に
より濃度が高められることになる。つまり、強調処理に
おいてその部分が強調されることとなり、同時にその部
分に含まれるノイズが強調されてしまうことになる。し
たがって、望ましくは図16に示す実施形態のように原
画像信号Sorgに基づいて各信号を作成するのがよい。
また処理時間の短縮という意味でも、2つの処理は並列
に行うことが望ましい。
Sproc = Sorg + β (Sorg) · Fusm (Sorg, Sus1, Sus2,... SusN) + D (Sorg−Fdrc (Sorg, Sus1, Sus2,... SusN)) (3) Fusm (Sorg, Sus1, Sus2, ... SusN) = {f u1 ( Sorg-Sus1) + f u2 (Sus1 -Sus2) + ... + f uk (Susk-1-Susk) + ... + f uN (SusN-1-SusN)} Fdrc (Sorg, Sus1, Sus2, ... SusN) = {f d1 (Sorg−Sus1) + f d2 (Sus1−Sus2) +... + F dk (Susk−1−Susk) +... + F dN (SusN−1−SusN)} (where Sproc: processed image signal Sorg: original image signal Susk (k = 1~N): unsharp image signals f uk (k = 1~N): transformation function f dk used for the first conversion processing (k = 1~N): second A conversion function β (Sorg) used in the conversion processing of: an enhancement coefficient determined based on the original image signal D (Sorg-Fdrc): a dynamic lens determined based on the low-frequency component signal Compression coefficient (D is a function for converting Sorg-Fdrc) Here, in the embodiment shown in FIG. 16, a signal used for frequency emphasis processing and a signal used for dynamic range compression processing are respectively based on the original image signal Sorg. The original image signal Sorg
In this case, one of the frequency emphasis processing and the dynamic range compression processing is performed first, and the other signal is processed for the other signal. However, a low-density portion of a radiographic image is generally a portion where the imaging dose at the time of imaging is small, and therefore contains a relatively large amount of noise. On the other hand, in the frequency emphasizing process, since the density-dependent emphasizing process in which the degree of emphasis is higher in the higher density portion is performed, if the frequency emphasizing process is directly performed on the original image signal, the lower density portion is not emphasized. , Noise is not emphasized. On the other hand, if, for example, the dynamic range compression process is performed on the original image signal and then the frequency emphasis process is performed, the density of the low density portion is increased by the dynamic range compression process. That is, the part is emphasized in the emphasis processing, and at the same time, the noise included in the part is emphasized. Therefore, it is desirable to generate each signal based on the original image signal Sorg as in the embodiment shown in FIG.
Also, in terms of shortening the processing time, it is desirable that the two processes be performed in parallel.

【0074】このように、ダイナミックレンジ圧縮処理
を行う場合においても、原画像信号Sorgの画素密度に
拘わらず処理済み画像信号Sprocの周波数応答特性が略
同一のものとなるように変換関数のパラメータを定義す
ることにより、原画像信号Sorgの画素密度に拘わらず
略同一の周波数応答特性を有する処理済み画像信号Spr
ocを得ることができる。
As described above, even when the dynamic range compression processing is performed, the parameters of the conversion function are set so that the frequency response characteristics of the processed image signal Sproc become substantially the same regardless of the pixel density of the original image signal Sorg. By definition, the processed image signal Spr having substantially the same frequency response characteristics regardless of the pixel density of the original image signal Sorg
You can get oc.

【0075】なお、上記実施形態においては、10本/m
mの画素密度の場合に、下記の表2に示すような変換関
数のパラメータを設定した場合、処理済み画像信号Spr
ocの周波数応答特性は図20に示すものとなる。ここ
で、図7から図9に示すように10本/mmの画素密度に
おいては1cycle/mmに相当する帯域制限画像信号が存
在するが、6.7本/mmおよび5本/mmの画素密度に
おいては1cycle/mmの帯域制限画像信号は存在しない
ため、全ての周波数帯域に亘って処理済み画像信号Spr
ocの周波数応答特性を画素密度に拘わらず一致させるこ
とができない。一方、処理済み画像信号Sprocを再生す
ることにより得られる画像においては高周波成分よりも
低周波成分の方が視覚的に目立つものとなる。したがっ
て、このような場合には高周波成分についてはある程度
まで一致させ、1cycle/mm以下の低周波成分において
周波数応答特性を一致させるように、6.7本/mmおよ
び5本/mmの画素密度における帯域制限画像信号に対
しては、表2に示すようなパラメータX,Yにより変換
関数を定義すればよい。6.7本/mmおよび5本/mm
の画素密度の帯域制限画像信号を表2に示すパラメータ
により定義された変換関数により変換して得られた処理
済み画像信号Sprocの周波数応答特性を図21および図
22にそれぞれ示す。図21および図22に示すよう
に、表2に示すパラメータにより定義された変換関数に
より帯域制限画像信号を変換すると、高周波成分につい
ては一致の程度は低いものの、視覚的に目立つ1cycle/
mm以下の低周波成分についてはその周波数応答特性は
一致したものとなる。したがって、画素密度に拘わらず
視覚的に目立つ低周波成分において周波数応答特性が一
致した画像を再現可能な処理済み画像信号Sprocを得る
ことができる。
In the above embodiment, 10 wires / m
If the parameters of the conversion function as shown in Table 2 below are set for a pixel density of m, the processed image signal Spr
The frequency response characteristic of oc is as shown in FIG. Here, as shown in FIGS. 7 to 9, at a pixel density of 10 lines / mm, there is a band-limited image signal corresponding to 1 cycle / mm, but at a pixel density of 6.7 lines / mm and 5 lines / mm, Since there is no band-limited image signal of 1 cycle / mm, the processed image signal Spr is processed over all frequency bands.
The frequency response characteristics of oc cannot be matched regardless of the pixel density. On the other hand, in an image obtained by reproducing the processed image signal Sproc, the low-frequency component is visually more noticeable than the high-frequency component. Therefore, in such a case, the band limitation at a pixel density of 6.7 lines / mm and 5 lines / mm is performed so that the high-frequency components match to some extent and the low-frequency components of 1 cycle / mm or less match the frequency response characteristics. For an image signal, a conversion function may be defined by parameters X and Y as shown in Table 2. 6.7 / mm and 5 / mm
21 and 22 show the frequency response characteristics of the processed image signal Sproc obtained by converting the band-limited image signal having the pixel density described above with the conversion function defined by the parameters shown in Table 2. As shown in FIGS. 21 and 22, when the band-limited image signal is converted by the conversion function defined by the parameters shown in Table 2, the degree of coincidence is low for high-frequency components, but 1 cycle / cycle that is visually conspicuous.
The frequency response characteristics of low-frequency components of less than mm are the same. Therefore, it is possible to obtain a processed image signal Sproc capable of reproducing an image whose frequency response characteristics match in a low-frequency component that is visually noticeable regardless of the pixel density.

【0076】[0076]

【表2】 また、上記実施形態においては、蓄積性蛍光体シートを
10本/mm、6.7本/mmおよび5本/mmの画素密度
により読み取ることにより得られた原画像信号Sorgに
対して、画素密度に拘わらず一定の周波数応答特性が得
られるように非線形処理を施しているが、原画像信号S
orgとしては蓄積性蛍光体シートから読み取られたもの
に限定されるものではなく、例えば被写体のサイズとサ
ンプリング間隔との関係(dpi)などのその画像信号
により表される画像の解像度が予め分かっていれば、い
かなる画像信号に対しても本発明を適用することができ
る。
[Table 2] Further, in the above embodiment, the stimulable phosphor sheet is
Non-linear processing is performed on the original image signal Sorg obtained by reading at a pixel density of 10 lines / mm, 6.7 lines / mm and 5 lines / mm so that a constant frequency response characteristic is obtained regardless of the pixel density. The original image signal S
The org is not limited to the one read from the stimulable phosphor sheet. For example, the resolution of the image represented by the image signal such as the relationship between the size of the subject and the sampling interval (dpi) is known in advance. Then, the present invention can be applied to any image signal.

【0077】さらに、上記実施形態においては、変換関
数を非線形関数として帯域制限画像信号に対して非線形
処理を施しているが、非線形処理に限定されるものでは
なく、変換関数を線形関数や定数としてもよい。この場
合、原画像信号の画素密度に応じて線形関数のパラメー
タあるいは定数の値が設定される。
Further, in the above-described embodiment, the band-limited image signal is subjected to the non-linear processing by using the conversion function as a non-linear function. However, the present invention is not limited to the non-linear processing. Is also good. In this case, a parameter or a constant value of a linear function is set according to the pixel density of the original image signal.

【0078】また、上記実施形態においては、原画像信
号からフィルタリングおよび補間拡大によりボケ画像信
号を求め、原画像信号およびボケ画像信号から帯域制限
画像信号を作成しているが、これに限定されるものでは
なく、例えば原画像信号をウェーブレット変換やラプラ
シアンピラミッドなどの手法により多重解像度に変換
し、変換された各解像度毎の画像信号からボケ画像信号
を作成して帯域制限画像信号を作成してもよい。
In the above embodiment, the blurred image signal is obtained from the original image signal by filtering and interpolation and enlargement, and the band-limited image signal is created from the original image signal and the blurred image signal. However, the present invention is not limited to this. For example, even if an original image signal is converted to multiple resolutions by a method such as wavelet transform or Laplacian pyramid, and a blurred image signal is created from the converted image signals for each resolution, a band-limited image signal is created. Good.

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

【図1】本発明の実施形態による画像処理装置の構成を
示す概略ブロック図
FIG. 1 is a schematic block diagram illustrating a configuration of an image processing apparatus according to an embodiment of the present invention.

【図2】ボケ画像信号作成手段の構成を示す概略ブロッ
ク図
FIG. 2 is a schematic block diagram illustrating a configuration of a blurred image signal creating unit.

【図3】フィルタリング処理に使用されるフィルタの一
例を示す図
FIG. 3 is a diagram illustrating an example of a filter used in a filtering process;

【図4】低解像度画像信号作成処理の詳細を示す図FIG. 4 is a diagram showing details of a low-resolution image signal creation process;

【図5】補間拡大処理に使用されるフィルタの一例を示
す図
FIG. 5 is a diagram illustrating an example of a filter used for interpolation enlargement processing;

【図6】周波数強調処理装置の一例を示す図FIG. 6 is a diagram illustrating an example of a frequency emphasis processing device.

【図7】帯域制限画像信号の周波数応答特性を示す図
(画素密度10本/mm)
FIG. 7 is a diagram showing a frequency response characteristic of a band-limited image signal (pixel density 10 lines / mm).

【図8】帯域制限画像信号の周波数応答特性を示す図
(画素密度6.7本/mm)
FIG. 8 is a diagram showing a frequency response characteristic of a band-limited image signal (pixel density: 6.7 lines / mm);

【図9】帯域制限画像信号の周波数応答特性を示す図
(画素密度5本/mm)
FIG. 9 is a diagram illustrating a frequency response characteristic of a band-limited image signal (pixel density: 5 lines / mm);

【図10】変換関数定義手段の構成を示す概略ブロック
FIG. 10 is a schematic block diagram illustrating a configuration of a conversion function definition unit.

【図11】本実施形態により処理がなされた処理済み画
像信号の周波数応答特性を示す図(画素密度10本/m
m)
FIG. 11 is a diagram illustrating a frequency response characteristic of a processed image signal processed according to the present embodiment (pixel density 10 lines / m);
m)

【図12】本実施形態により処理がなされた処理済み画
像信号の周波数応答特性を示す図(画素密度6.7本/m
m)
FIG. 12 is a diagram illustrating a frequency response characteristic of a processed image signal processed according to the present embodiment (pixel density: 6.7 lines / m);
m)

【図13】本実施形態により処理がなされた処理済み画
像信号の周波数応答特性を示す図(画素密度5本/m
m)
FIG. 13 is a diagram illustrating a frequency response characteristic of a processed image signal processed according to the present embodiment (pixel density 5 lines / m);
m)

【図14】周波数応答特性を指定するためのユーザイン
タフェースの一例を示す図
FIG. 14 is a diagram showing an example of a user interface for designating frequency response characteristics.

【図15】本実施形態の動作を示すフローチャートFIG. 15 is a flowchart showing the operation of the present embodiment.

【図16】ダイナミックレンジ圧縮処理装置の一例を示
す図
FIG. 16 is a diagram illustrating an example of a dynamic range compression processing device.

【図17】第1の変換手段における変換関数の一例を表
す図
FIG. 17 is a diagram illustrating an example of a conversion function in a first conversion unit.

【図18】第1の変換手段における変換関数の他の例を
表す図
FIG. 18 is a diagram illustrating another example of the conversion function in the first conversion unit.

【図19】第2の変換手段における変換関数の一例を表
す図
FIG. 19 is a diagram illustrating an example of a conversion function in a second conversion unit.

【図20】本実施形態により処理がなされた処理済み画
像信号の周波数応答特性を示す図(画素密度10本/m
m)
FIG. 20 is a diagram illustrating a frequency response characteristic of a processed image signal processed according to the present embodiment (pixel density 10 lines / m);
m)

【図21】本実施形態により処理がなされた処理済み画
像信号の周波数応答特性を示す図(画素密度6.7本/m
m)
FIG. 21 is a diagram illustrating a frequency response characteristic of a processed image signal processed according to the present embodiment (pixel density: 6.7 lines / m);
m)

【図22】本実施形態により処理がなされた処理済み画
像信号の周波数応答特性を示す図(画素密度5本/m
m)
FIG. 22 is a diagram illustrating a frequency response characteristic of a processed image signal processed according to the present embodiment (pixel density 5 lines / m);
m)

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

1 画像処理装置 2 ボケ画像信号作成手段 3 非線形処理手段 4 変換関数定義手段 5 画素密度情報入力手段 10 フィルタリング処理手段 11 補間処理手段 21 減算器 22 変換器 23 演算器 DESCRIPTION OF SYMBOLS 1 Image processing apparatus 2 Blur image signal creation means 3 Non-linear processing means 4 Conversion function definition means 5 Pixel density information input means 10 Filtering processing means 11 Interpolation processing means 21 Subtractor 22 Transformer 23 Computing unit

Claims (21)

【特許請求の範囲】[Claims] 【請求項1】 所定の画素密度を有する原画像信号か
ら、互いに異なる帯域制限特性を有する複数の画像信号
を作成し、該各画像信号に対して所定の変換関数に基づ
いて変換処理を施して変換画像信号を得、該変換画像信
号から処理済み画像信号を得る画像処理方法において、 前記所定の画素密度に基づいて前記変換関数の定義パラ
メータを決定して該変換関数を定義することを特徴とす
る画像処理方法。
An image signal having different band limiting characteristics is created from an original image signal having a predetermined pixel density, and a conversion process is performed on each of the image signals based on a predetermined conversion function. An image processing method for obtaining a converted image signal and obtaining a processed image signal from the converted image signal, wherein a definition parameter of the conversion function is determined based on the predetermined pixel density, and the conversion function is defined. Image processing method.
【請求項2】 前記変換関数の定義パラメータの決定
を、予め設定された少なくとも2以上の画素密度ごとに
定義パラメータを用意し、該用意した定義パラメータの
うち、前記予め設定された2以上の画素密度のうちの前
記所定の画素密度に最も近い画素密度に対応する定義パ
ラメータを選択することにより行なうことを特徴とする
請求項1記載の画像処理方法。
2. A method for determining a definition parameter of the conversion function, comprising: defining a definition parameter for each of at least two or more preset pixel densities; 2. The image processing method according to claim 1, wherein the step is performed by selecting a definition parameter corresponding to a pixel density closest to the predetermined pixel density among the densities.
【請求項3】 前記決定した定義パラメータを前記原
画像信号に関連づけて保存することを特徴とする請求項
1または2記載の画像処理方法。
3. The image processing method according to claim 1, wherein the determined definition parameters are stored in association with the original image signal.
【請求項4】 前記原画像信号に対して前記画素密度
に基づいて定められたフィルタ係数を有するフィルタに
よりフィルタリング処理を施して互いに周波数応答特性
が異なる複数のボケ画像信号を作成し、前記原画像信号
および前記複数のボケ画像信号に基づいて前記原画像信
号の複数の周波数帯域毎の信号を表す複数の帯域制限画
像信号を作成し、該複数の帯域制限画像信号を前記複数
の画像信号とすることを特徴とする請求項1から3いず
れか1項記載の画像処理方法。
4. The original image signal is subjected to a filtering process using a filter having a filter coefficient determined based on the pixel density to generate a plurality of blurred image signals having different frequency response characteristics from each other. Creating a plurality of band-limited image signals representing signals for each of a plurality of frequency bands of the original image signal based on the signals and the plurality of blurred image signals, and using the plurality of band-limited image signals as the plurality of image signals The image processing method according to any one of claims 1 to 3, wherein:
【請求項5】 前記所定の変換関数は、非線形関数で
あることを特徴とする請求項1から4いずれか1項記載
の画像処理方法。
5. The image processing method according to claim 1, wherein the predetermined conversion function is a non-linear function.
【請求項6】 前記変換処理は、周波数強調処理であ
ることを特徴とする請求項1から5いずれか1項記載の
画像処理方法。
6. The image processing method according to claim 1, wherein the conversion process is a frequency emphasis process.
【請求項7】 前記変換処理は、ダイナミックレンジ
圧縮処理であることを特徴とする請求項1から5いずれ
か1項記載の画像処理方法。
7. The image processing method according to claim 1, wherein the conversion processing is a dynamic range compression processing.
【請求項8】 所定の画素密度を有する原画像信号か
ら、互いに異なる帯域制限特性を有する複数の画像信号
を作成する画像信号作成手段と、該各画像信号に対して
所定の変換関数に基づいて変換処理を施して変換画像信
号を得、該変換画像信号から処理済み画像信号を得る変
換手段とを備えた画像処理装置において、 前記画素密度に基づいて前記変換関数の定義パラメータ
を決定して該変換関数を定義する変換関数定義手段を備
えたことを特徴とする画像処理装置。
8. An image signal generating means for generating a plurality of image signals having band limiting characteristics different from each other from an original image signal having a predetermined pixel density, based on a predetermined conversion function for each of the image signals. A conversion means for performing a conversion process to obtain a converted image signal and obtaining a processed image signal from the converted image signal, wherein a definition parameter of the conversion function is determined based on the pixel density. An image processing apparatus comprising a conversion function defining means for defining a conversion function.
【請求項9】 前記変換関数定義手段が、予め設定さ
れた少なくとも2以上の画素密度ごとに対応する定義パ
ラメータを格納したテーブルを有し、該テーブルに格納
された定義パラメータのうち、前記予め設定された2以
上の画素密度のうちの前記所定の画素密度に最も近い画
素密度に対応する定義パラメータを読み出すことにより
前記定義パラメータの決定を行なうものであることを特
徴とする請求項8記載の画像処理装置。
9. The conversion function defining means has a table storing definition parameters corresponding to at least two or more pixel densities set in advance, and among the definition parameters stored in the table, the preset function includes 9. The image according to claim 8, wherein the definition parameter is determined by reading a definition parameter corresponding to a pixel density closest to the predetermined pixel density among the two or more pixel densities obtained. Processing equipment.
【請求項10】 前記決定した定義パラメータを前記
原画像信号に関連づけて保存する記憶手段を備えたこと
を特徴とする請求項8または9記載の画像処理装置。
10. The image processing apparatus according to claim 8, further comprising a storage unit configured to store the determined definition parameter in association with the original image signal.
【請求項11】 前記画像信号作成手段は、前記原画
像信号に対して前記画素密度に基づいて定められたフィ
ルタ係数を有するフィルタによりフィルタリング処理を
施して互いに周波数応答特性が異なる複数のボケ画像信
号を作成するボケ画像信号作成手段と、前記原画像信号
および前記複数のボケ画像信号に基づいて、前記原画像
信号の複数の周波数帯域毎の信号を表す複数の帯域制限
画像信号を前記複数の画像信号として作成する帯域制限
画像信号作成手段とからなることを特徴とする請求項8
から10いずれか1項記載の画像処理装置。
11. The image signal generating means performs a filtering process on the original image signal with a filter having a filter coefficient determined based on the pixel density, and outputs a plurality of blurred image signals having different frequency response characteristics from each other. A plurality of band-limited image signals representing signals of a plurality of frequency bands of the original image signal based on the original image signal and the plurality of blur image signals. 9. A band-limited image signal creating means for creating a signal.
The image processing apparatus according to any one of claims 1 to 10.
【請求項12】 前記所定の変換関数は、非線形関数
であることを特徴とする請求項8から11いずれか1項
記載の画像処理装置。
12. The image processing apparatus according to claim 8, wherein the predetermined conversion function is a non-linear function.
【請求項13】 前記変換手段における変換処理は、
周波数強調処理であることを特徴とする請求項8から1
2いずれか1項記載の画像処理装置。
13. The conversion processing in the conversion means,
8. The method according to claim 8, wherein the frequency emphasis processing is performed.
3. The image processing apparatus according to claim 1.
【請求項14】 前記変換手段における変換処理は、
ダイナミックレンジ圧縮処理であることを特徴とする請
求項8から12いずれか1項記載の画像処理装置。
14. The conversion processing in the conversion means,
13. The image processing device according to claim 8, wherein the image processing is a dynamic range compression process.
【請求項15】 所定の画素密度を有する原画像信号
から、互いに異なる帯域制限特性を有する複数の画像信
号を作成し、該各画像信号に対して所定の変換関数に基
づいて変換処理を施して変換画像信号を得、該変換画像
信号から処理済み画像信号を得る画像処理方法をコンピ
ュータに実行させるためのプログラムを記録したコンピ
ュータ読取り可能な記録媒体において、 前記プログラムは、前記画素密度に基づいて前記変換関
数の定義パラメータを決定して該変換関数を定義する手
順を有することを特徴とするコンピュータ読取り可能な
記録媒体。
15. A plurality of image signals having different band limiting characteristics from an original image signal having a predetermined pixel density, and performing a conversion process on each of the image signals based on a predetermined conversion function. A computer-readable recording medium storing a program for causing a computer to execute an image processing method for obtaining a converted image signal and obtaining a processed image signal from the converted image signal, wherein the program is based on the pixel density A computer-readable recording medium having a procedure for determining a definition parameter of a conversion function and defining the conversion function.
【請求項16】 前記変換関数の定義パラメータの決
定して該変換関数を定義する手順は、予め設定された少
なくとも2以上の画素密度ごとに定義パラメータを用意
し、該用意した定義パラメータのうち、前記予め設定さ
れた2以上の画素密度のうちの前記所定の画素密度に最
も近い画素密度に対応する定義パラメータを選択するこ
とにより行なう手順であることを特徴とする請求項11
記載のコンピュータ読取り可能な記録媒体。
16. The step of determining the definition parameter of the conversion function and defining the conversion function includes preparing a definition parameter for each of at least two or more pixel densities set in advance, and 12. A procedure performed by selecting a definition parameter corresponding to a pixel density closest to the predetermined pixel density among the two or more predetermined pixel densities.
The computer-readable recording medium according to claim 1.
【請求項17】 前記決定した定義パラメータを前記
原画像信号に関連づけて保存する手順を有することを特
徴とする請求項11または12記載のコンピュータ読取
り可能な記録媒体。
17. The computer-readable recording medium according to claim 11, further comprising a step of storing the determined definition parameter in association with the original image signal.
【請求項18】 前記複数の画像信号を作成する手順
は、前記原画像信号に対して前記画素密度に基づいて定
められたフィルタ係数を有するフィルタによりフィルタ
リング処理を施して互いに周波数応答特性が異なる複数
のボケ画像信号を作成する手順と、前記原画像信号およ
び前記複数のボケ画像信号に基づいて、前記原画像信号
の複数の周波数帯域毎の信号を表す複数の帯域制限画像
信号を前記複数の画像信号として作成する手順とからな
ることを特徴とする請求項15から17いずれか1項記
載のコンピュータ読取り可能な記録媒体。11記載のコ
ンピュータ読取り可能な記録媒体。
18. The method of producing a plurality of image signals, wherein the original image signal is subjected to a filtering process using a filter having a filter coefficient determined based on the pixel density, and the plurality of image signals have different frequency response characteristics. And generating a plurality of band-limited image signals representing signals for a plurality of frequency bands of the original image signal based on the original image signal and the plurality of blur image signals. 18. The computer-readable recording medium according to claim 15, wherein the recording medium is a signal. 12. The computer-readable recording medium according to item 11.
【請求項19】 前記所定の変換関数は、非線形関数
であることを特徴とする請求項15から18いずれか1
項記載のコンピュータ読取り可能な記録媒体。
19. The method according to claim 15, wherein the predetermined conversion function is a non-linear function.
A computer-readable recording medium according to claim 1.
【請求項20】 前記変換処理は、周波数強調処理で
あることを特徴とする請求項15から19いずれか1項
記載のコンピュータ読取り可能な記録媒体。
20. The computer-readable recording medium according to claim 15, wherein said conversion process is a frequency emphasis process.
【請求項21】 前記変換処理は、ダイナミックレン
ジ圧縮処理であることを特徴とする請求項15から19
いずれか1項記載のコンピュータ読取り可能な記録媒
体。
21. The method according to claim 15, wherein the conversion processing is a dynamic range compression processing.
The computer-readable recording medium according to claim 1.
JP11363655A 1999-01-14 1999-12-22 Picture processing method and picture processor and record medium Pending JP2000268171A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP11363655A JP2000268171A (en) 1999-01-14 1999-12-22 Picture processing method and picture processor and record medium

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP11-7293 1999-01-14
JP729399 1999-01-14
JP11363655A JP2000268171A (en) 1999-01-14 1999-12-22 Picture processing method and picture processor and record medium

Publications (1)

Publication Number Publication Date
JP2000268171A true JP2000268171A (en) 2000-09-29

Family

ID=26341570

Family Applications (1)

Application Number Title Priority Date Filing Date
JP11363655A Pending JP2000268171A (en) 1999-01-14 1999-12-22 Picture processing method and picture processor and record medium

Country Status (1)

Country Link
JP (1) JP2000268171A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011151985A (en) * 2010-01-22 2011-08-04 Mitsubishi Electric Corp Device for setting of setting value, and protective relay device
CN109564678A (en) * 2016-07-04 2019-04-02 Hoya株式会社 Image processing apparatus and electronic endoscope system

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61112285A (en) * 1984-11-07 1986-05-30 Hitachi Ltd Conversational image emphasis system
JPH01106578A (en) * 1987-10-19 1989-04-24 Fuji Photo Film Co Ltd Method for setting picture processing data
JPH05143727A (en) * 1991-11-19 1993-06-11 Fuji Photo Film Co Ltd Method for processing image data and device therefor
JPH06292009A (en) * 1993-04-02 1994-10-18 Konica Corp Method for processing radiation picture
JPH0998293A (en) * 1995-09-29 1997-04-08 Fuji Photo Film Co Ltd Image processing method and device
JPH1063838A (en) * 1996-08-19 1998-03-06 Fuji Photo Film Co Ltd Method for processing picture and device therefor
JPH1063837A (en) * 1996-08-19 1998-03-06 Fuji Photo Film Co Ltd Method for processing picture and device therefor
JPH1075364A (en) * 1995-09-29 1998-03-17 Fuji Photo Film Co Ltd Dynamic range compression processing method and device for image
JPH1075395A (en) * 1995-09-29 1998-03-17 Fuji Photo Film Co Ltd Image processing method and device
JPH10105701A (en) * 1996-09-25 1998-04-24 Fuji Photo Film Co Ltd Method and device for radio graph emphasis processing
JPH10171983A (en) * 1996-12-13 1998-06-26 Fuji Photo Film Co Ltd Method and device for processing picture
JPH10200756A (en) * 1996-11-18 1998-07-31 Seiko Epson Corp Image processing unit, image processing method, and medium recording image processing program

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61112285A (en) * 1984-11-07 1986-05-30 Hitachi Ltd Conversational image emphasis system
JPH01106578A (en) * 1987-10-19 1989-04-24 Fuji Photo Film Co Ltd Method for setting picture processing data
JPH05143727A (en) * 1991-11-19 1993-06-11 Fuji Photo Film Co Ltd Method for processing image data and device therefor
JPH06292009A (en) * 1993-04-02 1994-10-18 Konica Corp Method for processing radiation picture
JPH0998293A (en) * 1995-09-29 1997-04-08 Fuji Photo Film Co Ltd Image processing method and device
JPH1075364A (en) * 1995-09-29 1998-03-17 Fuji Photo Film Co Ltd Dynamic range compression processing method and device for image
JPH1075395A (en) * 1995-09-29 1998-03-17 Fuji Photo Film Co Ltd Image processing method and device
JPH1063838A (en) * 1996-08-19 1998-03-06 Fuji Photo Film Co Ltd Method for processing picture and device therefor
JPH1063837A (en) * 1996-08-19 1998-03-06 Fuji Photo Film Co Ltd Method for processing picture and device therefor
JPH10105701A (en) * 1996-09-25 1998-04-24 Fuji Photo Film Co Ltd Method and device for radio graph emphasis processing
JPH10200756A (en) * 1996-11-18 1998-07-31 Seiko Epson Corp Image processing unit, image processing method, and medium recording image processing program
JPH10171983A (en) * 1996-12-13 1998-06-26 Fuji Photo Film Co Ltd Method and device for processing picture

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011151985A (en) * 2010-01-22 2011-08-04 Mitsubishi Electric Corp Device for setting of setting value, and protective relay device
CN109564678A (en) * 2016-07-04 2019-04-02 Hoya株式会社 Image processing apparatus and electronic endoscope system

Similar Documents

Publication Publication Date Title
US7221806B2 (en) Method of and system for image processing and recording medium for carrying out the method
JP3816151B2 (en) Image processing method and apparatus
US20070086671A1 (en) Image processing apparatus
JP3700804B2 (en) Image processing method and apparatus
JP3738788B2 (en) Image dynamic range compression processing method and apparatus
JP3738791B2 (en) Image processing method and apparatus
JP3999432B2 (en) Image processing method and apparatus, and recording medium
JP3696339B2 (en) Image processing method and apparatus
JP2007180851A (en) Gray scale transformation device, program and method for raw image, and electronic camera
JP6174903B2 (en) Spatial / gradation super-resolution device and program
JPH09149241A (en) Method and device for enlarging picture
JPH06325165A (en) Display method for part of radiation picture
JP4598115B2 (en) Image processing method and apparatus, and recording medium
JP2000268171A (en) Picture processing method and picture processor and record medium
JP2002183727A (en) Image processing device
JP3783116B2 (en) Radiation image enhancement processing method and apparatus
JP3965302B2 (en) Spatial filtering method
JP2002074356A (en) Method and device for processing picture and recording medium
JP4267159B2 (en) Image processing method and apparatus, and recording medium
JP2008211840A (en) Image processing apparatus
JP2011049696A (en) Image processing device, and image processing method
JP3669452B2 (en) Image processing method and apparatus
JP4213337B2 (en) Image processing method and apparatus, and recording medium
JP3557673B2 (en) Image processing device
JPH10150572A (en) Image-processing unit, image-processing method, and medium recording image-processing program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040823

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20061201

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20070713

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070724

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070925

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20071211