JPS63303474A - Digital image processing method - Google Patents
Digital image processing methodInfo
- 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
Links
- 238000003672 processing method Methods 0.000 title claims description 5
- 238000000034 method Methods 0.000 claims description 7
- 230000006835 compression Effects 0.000 abstract description 3
- 238000007906 compression Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
Abstract
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.
第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)
(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.
(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.▼
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)
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 |
-
1987
- 1987-06-04 JP JP62138873A patent/JPH0624009B2/en not_active Expired - Lifetime
Cited By (1)
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 |