JPS63303474A - Digital image processing method - Google Patents

Digital image processing method

Info

Publication number
JPS63303474A
JPS63303474A JP62138873A JP13887387A JPS63303474A JP S63303474 A JPS63303474 A JP S63303474A JP 62138873 A JP62138873 A JP 62138873A JP 13887387 A JP13887387 A JP 13887387A JP S63303474 A JPS63303474 A JP S63303474A
Authority
JP
Japan
Prior art keywords
input image
equation
image
value
maximum value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP62138873A
Other languages
Japanese (ja)
Other versions
JPH0624009B2 (en
Inventor
Shunichi Kohama
小浜 俊一
Hidehiro Tachiki
立木 秀広
Toshio Chiba
千葉 敏夫
Susumu Matsuzaki
進 松崎
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.)
Japan Steel Works Ltd
Oki Electric Industry Co Ltd
Technical Research and Development Institute of Japan Defence Agency
Original Assignee
Japan Steel Works Ltd
Oki Electric Industry Co Ltd
Technical Research and Development Institute of Japan Defence Agency
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 Japan Steel Works Ltd, Oki Electric Industry Co Ltd, Technical Research and Development Institute of Japan Defence Agency filed Critical Japan Steel Works Ltd
Priority to JP62138873A priority Critical patent/JPH0624009B2/en
Publication of JPS63303474A publication Critical patent/JPS63303474A/en
Publication of JPH0624009B2 publication Critical patent/JPH0624009B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Abstract

PURPOSE:To realize sharpness corresponding to the chevron pattern of a light- and-shade level using the maximum value in the neighborhood of each image element, by applying a specific equation on the density value of an input image. CONSTITUTION:To a processor 10, the input image F={fi, j} is inputted. In such a case, fi, j represents the density value in the image element (i, j), and (i) and (j) represent natural numbers. the processor 10 performs a processing shown in equation I on the input F, and outputs an output in which natural sharpness corresponding to the chevron pattern of the input image is applied to a display medium 12. Also, it performs the processing by equation II for the result of the equation I, then, performs the compression of a dynamic range. In such a case Ni, j represents the set of values fi, j in the neighborhood of a point (i, j), and Max(Ni, j) the maximum value of the Ni, j, and (alpha) and (beta) represent parameters possible to be changed from a console 14.

Description

【発明の詳細な説明】 (産業上の利用分野) 本発明は入力画像に対してレベルのダイナミックレンジ
の圧縮、規格化および尖鋭化を同時に行うことかできる
ディジタル画像処理方法に関するちのである。
DETAILED DESCRIPTION OF THE INVENTION (Field of Industrial Application) The present invention relates to a digital image processing method that can simultaneously compress, normalize, and sharpen the level dynamic range of an input image.

(従来の技術) 簡単のため一次元画像の場合について説明する。(Conventional technology) For simplicity, the case of a one-dimensional image will be explained.

従来尖鋭化の技術としては微分処理が広く用いられてい
る。この方法は人力濃度値列(×o) に対し出力濃度
値列(y、)を Yn”Xn−172(Xn−1”Xn+1)     
    (+)Jl=正整数 によって算出する方式である。第2図に人力(X、)第
3図に(1)式の結果(y、、)を示す。
Differential processing has been widely used as a conventional sharpening technique. This method converts the output concentration value sequence (y,) to the human concentration value sequence (xo) to Yn"Xn-172 (Xn-1"Xn+1)
This is a calculation method using (+)Jl=positive integer. Figure 2 shows the human power (X,), and Figure 3 shows the result (y,,) of equation (1).

(発明か解決しようとする問題点) 上述の従来の方法では′c4淡レベルの作る山の形状に
対して1が適切に選択されているところでは、効果があ
るが(図中(() (a) (八) ) iか不適当な
場合不自然な処理となる(図中(ニ)の部分)。
(Problem to be solved by the invention) The above-mentioned conventional method is effective when 1 is appropriately selected for the shape of the mountain created by the 'c4 light level (() in the figure). a) (8)) If i is inappropriate, the processing will be unnatural (part (d) in the figure).

また、この方法のみではレベルのダイナミックレンジの
圧縮効果は不モ分である。
Moreover, the effect of compressing the dynamic range of the level is insufficient if only this method is used.

本発明は上述の欠点を改牌するととも・・・、k、ダイ
ナミックレンジを圧縮するティジタル仏号処理方式を提
供することを目的とする。
It is an object of the present invention to overcome the above-mentioned drawbacks and provide a digital code processing method that compresses the dynamic range.

(問題点を解決するための手段) ト記[目的を達成するだめの本発明の特徴は、入力画像
F・(fij) (fijは画素(i、j)における濃
度値、i=1,2.・・・、k、 j=1.z、・−、
L、 k、Lは自然数)に対し、下記(1)式、又は(
1)式及び(2)式を適用した出力画像を得ることにあ
る。
(Means for Solving the Problems) G [The features of the present invention to achieve the purpose are that the input image F.(fij) (fij is the density value at pixel (i, j), i=1,2 ..., k, j=1.z, ・-,
L, k, L are natural numbers), the following formula (1) or (
The objective is to obtain an output image by applying equations 1) and 2).

0≦α〈l Ni、+は点(i、j)の近傍における値fijの集合
M、、X(N、)はNijの最大値 βΔ十(1−α)M、、(Nij) lk    L O≦β、A=−Σ Σhij kLi=l j=1 (作用) L記(1)式を適用することにより各画素点の近傍にお
ける最大値を用いた?Q 7Aレベルの山の形状に応し
た自然な尖鋭化を図ることかできる。ここで、近傍点N
ijの範囲及びαの値は所望の結果が得られるようにオ
ペレータか装置を操作して調節することができる。さら
・・・、k、(1)式の結果に対して、(2)式を適用
することにより、全域の平均値と局所近傍における最大
値の和で各濃度値を規格化して、ダイナミックレンジの
圧縮をすることができる。
0≦α〈l Ni,+ is the set M of values fij in the vicinity of point (i, j), , X(N,) is the maximum value of Nij βΔ0(1-α)M, , (Nij) lk L O≦β, A=−Σ Σhij kLi=l j=1 (Operation) By applying equation (1) in L, the maximum value in the vicinity of each pixel point is used. Q: Is it possible to create a natural sharpening that corresponds to the shape of the mountain at the 7A level? Here, the neighboring point N
The range of ij and the value of α can be adjusted by the operator or the device to obtain the desired results. Furthermore, k, by applying equation (2) to the result of equation (1), each density value is normalized by the sum of the average value of the entire area and the maximum value in the local vicinity, and the dynamic range is calculated. can be compressed.

(実施例) 画像処理システムの例を第1図に示す。ディジタルプロ
セッサIOは、通常のプログラムされたコンピュータま
たは−Y用のプロセッサ、表示媒体■2はグラフィック
ディスプレイの他、al SA紙等が可能である。ここ
でパラメータα、β+NiJはコンソール14からの人
力Rより適宜変更できる。
(Example) An example of an image processing system is shown in FIG. The digital processor IO can be a normal programmed computer or a processor for -Y, and the display medium 2 can be a graphic display, ALSA paper, etc. Here, the parameters α, β+NiJ can be changed as appropriate using human power R from the console 14.

プロセッサ10には入力画像F=(f、J)か人力され
る。ここでfijは画素(i、j)における濃度値でi
=1、2、−、k、j=1、2、・・・、L、 kとL
は自然数である。
An input image F=(f, J) is input to the processor 10 manually. Here fij is the density value at pixel (i, j)
=1,2,-,k,j=1,2,...,L, k and L
is a natural number.

プロセッサlOに人力1?に対し、(1)式で示される
処理を行って入力画像の山の形状に応じた自然な尖鋭化
を行った出力を表示媒体12に出力する。
Processor IO and 1 human power? Then, the process shown in equation (1) is performed to output an output to the display medium 12 that has been naturally sharpened according to the shape of the mountain in the input image.

([1J−αM、、X(N、J) kQ、・・f iJ−αMaX(N、)<0のとき0≦
αく1 Nijは点い、j)の近傍における値f1Jの集合M、
X(N、、)はNijの最大値 次いで、(1)式の結果に対し、(2)式による処理を
行って、ダイナミックレンジの圧縮を行う。
([1J-αM,,
αku1 Nij is a point, the set M of values f1J in the vicinity of j),
X(N, ,) is the maximum value of Nij.Then, the result of equation (1) is processed by equation (2) to compress the dynamic range.

βΔ十(l−α)M、、(Nij) lk    L O≦β、A=−Σ Σh+j i+Li=1 jl (2)式における分母は、全域の濃度の羽均値と、局所
近傍における最大値との各々に適当な係数をかけた値の
和であり、この和により各濃度値[H,は規格化される
βΔ10(l−α)M,, (Nij) lk L O≦β, A=−Σ Σh+j i+Li=1 jl The denominator in equation (2) is the average value of the concentration over the entire area and the maximum value in the local vicinity. is the sum of the values obtained by multiplying each of them by an appropriate coefficient, and each density value [H, is normalized by this sum.

(1)弐及び(2)式における変数は、コンソール14
から、所望の結果が得られるようにオペレータか人力す
る。従って、変数の組合わせに従って複数の出力画像か
得られる。
The variables in formulas (1) and (2) are
From there, the operator or human labor is used to obtain the desired results. Therefore, multiple output images can be obtained according to the combination of variables.

たとえばα=0.N、Jを全域、β−0とすれば原画像
そのままの表示となり、またNijを隣接要素のみとし
αを1に十分近ずければ極大点のみの表示となる。また
、βを変化させることによりダイナミックレンジの圧縮
特性を制御できる。
For example, α=0. If N and J are the entire area and β-0, the original image will be displayed as is, and if Nij is only the adjacent element and α is sufficiently close to 1, only the maximum point will be displayed. Furthermore, by changing β, the compression characteristics of the dynamic range can be controlled.

本発明を用いて、第2図の人力を処理した結果が第4図
に描かれている。(イ)(ロ)(ハ)については従来の
微分処理と同様の結果かえられている上、すその広い山
である(:)も不自然につぶされることなく忠実に尖鋭
化さねていることかわかる。なお、第4図ではα−0,
75,β=1、N、Jの近傍は前後2点と自己の点の計
3点としてあり、−次元表示なので添字(i、j)の代
わりにnで表示される。
The results of processing the human effort in Figure 2 using the present invention are depicted in Figure 4. Regarding (a), (b), and (c), the results are the same as in conventional differential processing, and the wide mountain at the base (:) is not unnaturally crushed and is faithfully sharpened. I understand that. In addition, in Fig. 4, α-0,
There are three points in the vicinity of 75, β=1, N, J, two points before and after, and the own point, and since it is a -dimensional display, it is indicated by n instead of the subscript (i, j).

以上説明したように本発明によると、すその狭い(濃淡
レベルの)山に対してもすその広い山に対しても自然な
尖鋭化かなされている。つまり、微分処理では、すその
幅が適切でない山は不自然な変形がなされるのに対し、
本方式ではそのような弊害は生じない。
As explained above, according to the present invention, both mountains with narrow bases (gradation level) and mountains with wide bases are naturally sharpened. In other words, in differential processing, a mountain whose base width is not appropriate is deformed unnaturally;
This method does not cause such a problem.

また、第1図に示すシステム例において、?■純なパラ
メータ変化たけて原画像表示や極大点表示なと多様な出
力画像が1itられるのもハードウェア及びソフトウェ
ア上の利点となる。
Moreover, in the system example shown in FIG. - It is also an advantage in terms of hardware and software that various output images such as original image display and local maximum point display can be produced in one unit by simply changing parameters.

(発明の効果) 本発明の方法では、近傍における最大値を用いることか
ら原画像に忠実で、自然な尖鋭化がなされる。また、パ
ラメータα、β、N+jを変化させることにより、ダイ
ナミックレンジの圧縮、コントラストの強調等多様なバ
リエーションか簡単に得られる。従って表示処理を用す
るすべてのシステムに応用できる。
(Effects of the Invention) In the method of the present invention, since the maximum value in the vicinity is used, natural sharpening is achieved that is faithful to the original image. Further, by changing the parameters α, β, and N+j, various variations such as dynamic range compression and contrast enhancement can be easily obtained. Therefore, it can be applied to all systems that use display processing.

【図面の簡単な説明】[Brief explanation of the drawing]

第1図は本発明によるディジタル画像処理システムの構
成図、第2図は入力画像の例、第3図は従来の技術によ
り第2図の入力画像に処理を施した結果を示す図、第4
図は本発明により第2図の入力画像に処理を施した結果
を示す図である。 10:演算器、12:表示媒体、14;コンソール、イ
9ロ、ハ、二二濃淡レベル。
FIG. 1 is a block diagram of a digital image processing system according to the present invention, FIG. 2 is an example of an input image, FIG. 3 is a diagram showing the result of processing the input image of FIG. 2 using a conventional technique, and FIG.
The figure is a diagram showing the result of processing the input image of FIG. 2 according to the present invention. 10: Arithmetic unit, 12: Display medium, 14: Console, I9, B, C, 22 gray level.

Claims (2)

【特許請求の範囲】[Claims] (1)入力画像F=(f_i_j)(f_i_jは画素
(i、j)における濃度値、i=1、2、・・・、k、
j=1、2、・・・、L、k、Lは自然数)に対し次の
処理を施して入力画像の山の形状に応じた自然な尖鋭化
を行った出力画像H=(h_i_j)を得ることを特徴
とする画像処理方法; ▲数式、化学式、表等があります▼(1) 0≦α<1 N_i_jは点(i、j)の近傍における値f_i_j
の集合M_a_x(N_i_j)はN_i_jの最大値
(1) Input image F = (f_i_j) (f_i_j is the density value at pixel (i, j), i = 1, 2, ..., k,
j = 1, 2, ..., L, k, L are natural numbers), the output image H = (h_i_j) is obtained by performing the following processing and natural sharpening according to the shape of the mountain in the input image. Image processing method characterized by obtaining; ▲Mathematical formulas, chemical formulas, tables, etc.▼(1) 0≦α<1 N_i_j is the value f_i_j near the point (i, j)
The set M_a_x(N_i_j) is the maximum value of N_i_j.
(2)入力画像F=(f_i_j)(f_i_jは画素
(i、j)における濃度値、i=1、2、・・・、k、
j=1、2、・・・、L、k、Lは自然数)に対し(1
)式による処理を施して入力画像の山の形状に応じた自
然な尖鋭化を行った出力画像H=(h_i_j)を得、
次いでHに対し(2)式を適用してダイナミックレンジ
を圧縮した出力G(g_i_j)を得ることを特徴とす
る画像処理方法; ▲数式、化学式、表等があります▼(1) 0≦α<1 N_i_jは点(i、j)の近傍における値f_i_j
の集合M_a_x(N_i_j)はN_i_jの最大値
g_i_j=h_i_j/[βA+(1−α)M_a_
x(N_i_j)](2)▲数式、化学式、表等があり
ます▼
(2) Input image F = (f_i_j) (f_i_j is the density value at pixel (i, j), i = 1, 2, ..., k,
j=1, 2, ..., L, k, L is a natural number), (1
) process to obtain an output image H = (h_i_j) that has been naturally sharpened according to the shape of the mountain in the input image,
An image processing method characterized by applying equation (2) to H to obtain an output G (g_i_j) with a compressed dynamic range; ▲Mathematical formulas, chemical formulas, tables, etc.▼(1) 0≦α< 1 N_i_j is the value f_i_j near the point (i, j)
The set M_a_x(N_i_j) is the maximum value of N_i_j g_i_j=h_i_j/[βA+(1-α)M_a_
x(N_i_j)](2)▲There are mathematical formulas, chemical formulas, tables, etc.▼
JP62138873A 1987-06-04 1987-06-04 Digital image processing method Expired - Lifetime JPH0624009B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP62138873A JPH0624009B2 (en) 1987-06-04 1987-06-04 Digital image processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP62138873A JPH0624009B2 (en) 1987-06-04 1987-06-04 Digital image processing method

Publications (2)

Publication Number Publication Date
JPS63303474A true JPS63303474A (en) 1988-12-12
JPH0624009B2 JPH0624009B2 (en) 1994-03-30

Family

ID=15232106

Family Applications (1)

Application Number Title Priority Date Filing Date
JP62138873A Expired - Lifetime JPH0624009B2 (en) 1987-06-04 1987-06-04 Digital image processing method

Country Status (1)

Country Link
JP (1) JPH0624009B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19722761B4 (en) * 1997-06-02 2006-06-01 Hell Gravure Systems Gmbh Process for processing image data

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19722761B4 (en) * 1997-06-02 2006-06-01 Hell Gravure Systems Gmbh Process for processing image data

Also Published As

Publication number Publication date
JPH0624009B2 (en) 1994-03-30

Similar Documents

Publication Publication Date Title
US5781665A (en) Apparatus and method for cropping an image
EP0352016B1 (en) Method and system for enhancement of a digitized image
Pavlidis Algorithms for graphics and image processing
Rahman et al. A comparison of the multiscale retinex with other image enhancement techniques
Rahman et al. Multiscale retinex for color rendition and dynamic range compression
EP1412920B1 (en) A general purpose image enhancement algorithm which augments the visual perception of detail in digital images
US7715656B2 (en) Magnification and pinching of two-dimensional images
US6516087B1 (en) Method for real time correlation of stereo images
Duan et al. Fast tone mapping for high dynamic range images
Rahman Properties of a center/surround Retinex: Part 1: Signal processing design
Dibos et al. Global total variation minimization
JPS60169980A (en) Method and apparatus for correcting deterioration in contrast of image
JPH03102479A (en) Method and apparatus for processing image
JPS6041168A (en) Vide reprojection apparatus and method
JPH09149241A (en) Method and device for enlarging picture
JPS63303474A (en) Digital image processing method
JPH0998293A (en) Image processing method and device
US6453075B1 (en) Method for preserving image detail when adjusting the tone scale of a digital color image
JP3206217B2 (en) Image processing device
JP2973432B2 (en) Image processing method and apparatus
Gao et al. Multiscale decomposition based high dynamic range tone mapping method using guided image filter
JPH02106734A (en) Picture composing device
JPH03265367A (en) Display processing system for reduced picture
JP2916607B2 (en) Image magnifier
Vrscay A hitchhiker’s guide to ‘fractal-based’function approximation and image compression

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080330

Year of fee payment: 14