JPS60150183A - Corrector for digital picture distortion - Google Patents

Corrector for digital picture distortion

Info

Publication number
JPS60150183A
JPS60150183A JP59264268A JP26426884A JPS60150183A JP S60150183 A JPS60150183 A JP S60150183A JP 59264268 A JP59264268 A JP 59264268A JP 26426884 A JP26426884 A JP 26426884A JP S60150183 A JPS60150183 A JP S60150183A
Authority
JP
Japan
Prior art keywords
image
picture
point
coefficient
digital
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
JP59264268A
Other languages
Japanese (ja)
Other versions
JPH0127467B2 (en
Inventor
Mitsutaka Kosaka
小坂 満隆
Koichi Honma
弘一 本間
Fuminobu Furumura
文伸 古村
Akira Tsuboi
坪井 晃
Koichi Ihara
井原 広一
Yutaka Kubo
裕 久保
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.)
Hitachi Ltd
Original Assignee
Hitachi 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 Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP59264268A priority Critical patent/JPS60150183A/en
Publication of JPS60150183A publication Critical patent/JPS60150183A/en
Publication of JPH0127467B2 publication Critical patent/JPH0127467B2/ja
Granted legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/40Image enhancement or restoration by the use of histogram techniques
    • G06T5/80

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)

Abstract

PURPOSE:To correct the density variance of picture data with high accuracy by obtaining a coefficient at each grid point of a split picture for linear approximation and calculating the coefficient at each point from said coefficient at each grid point and the relative coordinates within a block to perform corrections of the digital picture distortion. CONSTITUTION:A picture 21 free from density variance and a picture 23 having density variance are supplied from digital picture input devices 22 and 24 by the control signal given from a processor 25. In this case, a peripheral histogram of a reference point P set on a picture having density variance and a peripheral histogram of a point Q on a picture free from density variance and opposite to the point P are calculated previously by means of the devices 22 and 24. The picture having the density variance is transferred to a picture memory 26 as soon as the histogram is calculated. The processor 25 sends the peripheral histograms of P and Q at grid points that are divided by (mXn) on a picture to a corrector 27 to obtain the coincidence between both histograms. Thus the density variance of pictures can be removed.

Description

【発明の詳細な説明】 【発明の利用分野〕 本発明は、ディジタル画像歪補正装置に関し、詳しくは
、濃淡むらの少ないセンサで撮った画像データを利用し
、あるいはあらかじめ陰面線近似の係数を記憶すること
により、入力画像のd(淡歪を補正する装置に関するも
のである。 〔発明の背景〕 ディジタルil+[機を用いて画像を解析・処理するた
めには、先ず画像の各点の強度を示すアナログ信号を一
?゛イジタル信号に変換してd1算機に入力することが
必要である。 第1図は、従来の画像入力装置の概略図である。 第1図において、1は被写体52はセンサ、3はアナロ
グ・ディジタル変換器、4は記憶装置である。 被写体lの放射強度は、センサ2により2次元的強度分
布として画像8のように形成さJし、その各点の強度を
示すアナログ信号がアナログ・ティシタル変換器3によ
り標本化、i子化されたディジタル(11号に変換され
て、記憶装置4に読み込ま11、る。 第1図に示ずセンサ2では、第2図のように、被写体か
らの光線11が、レンズ等の光学系12により、焦点面
にあるビジコン等の感光面13に電荷とし・て像を゛結
ぶ。読出装置]5は、電子ビーム14によりこの444
Bを読み出し、電圧に変えて信号線1Gを介し゛C送出
する。 ところで、第1図に示す入力装置により取り込]:れた
画像は、センサ2の感光面13の各点の感度むI)に起
因した濃淡むらを含んでいる。つまり、被写(に1に均
一な明るさ分布を持つ光源を使用しても、センサ出力画
像8の濃淡強度分布、あるいは記1.α装買4に読み込
まれた画像9の濃淡強度分布は1面的に均一とならずに
、むらを生ずる。 そこで、画像データを清書な解析、処理に利用するため
には、この濃淡むらを取り除くことが必要となる。 従来1画像の濃淡むらを取り除くためには、第1図に示
すように、変換装置6.記憶装置7および変換表をυ、
える入力装置5よりなる補正装置を、入力装置の後段に
接続する。 変換装置6は、記憶装置4にl(′1えられた各画像デ
ータ、つまり標本化された画像単位の強度を、変換表に
したがって濃淡むら、を取り除いた強度に変換する。2
変換表は、均一な明るさを持つ光源の被写体1に対して
、画像のa淡強度分布を一様にするだめのもので、この
変換表にしたがって変換された画素データは、記憶装置
7に損えられる。 この変換表は、第3図(a>(b) (c)の順序で作
成される。 第3図(11)では、センサ2が取り込むことのできる
入力M刺強度の最大値をXよ(単位を例えばcd)とし
、既タイ1の均一光源を被写体として、入力放射強度X
に勾するセンサ出力信号線16に現われる出カフ6圧を
y(単位を例えばV)とすると、現像の各画素に列J、
bするセンサ面上の位置で(x 、’ y’)の関係が
得られる。すなわち、第3図(’a)の各点17は、放
射強度Xを0からx3まで変化させて。 い<−)かの光源を用いて測定した場合の1光源の各画
素に文1応するS定データである。 第31”I (b )は、アナログ・ディジタル変換器
の特性をモデル化したものである。センサ出力電圧の範
1/Jlを0から)・いまで変化した場合の、変換後の
ディジタル・データのレベル数を1−とすれは、センサ
11シカ電圧yに対するアナログ・ディジタル変換!I
:+の出力Zは、次式で与えられれる。 Z = (1,、X −3・・(1) yIlツ ここで、(、〕は、非負実数aを超えない最大整数であ
る。 補正装置において、濃淡むらを取り除くために人力装置
5からLy、えられる変数表は、次のようにしてめられ
る。 第3図(、)の各点17について、濃淡むらを取り除い
た補正強度Uを、次式で表すものとする。 u = (L X −3・・(2) X 暉 第3図(c)は、点j7の各々に対応して、アナログ・
ディジタル変換器出力2と補正強度Uとの関係を示す点
+9.および点19を結んだ曲線20を示す十ノのであ
る。 第3図(c)の各点19を真綿補間して得られる関数2
0は、次式でノjえられる。 u = f(7,) (:s) また、変換表は、変数ZのレヘルO,I・・・・l、−
1に苅し“C1次式で表わされる。 1」= (r (7,):l ・・(4)通i;tのヒ
シコ〉・センサでは、−に記(3)式の値、3十分ノ1
′:精度で線形関数により近似することかできろ1.こ
のJ、うにして、画像の一各画素について、以にの1・
°六Ii”r”A+うことにより所望の変換表が得ら4
し る 。 変換表?r用いて人力画像の、11ル淡むらを補正する
従来の方法は、センサ2の澹淡むら特性が完全にU+り
定でき、かツi*+記(11)のu=(f(Z))にお
(づる関数1が完全にモデル化できるという前提のもと
に成立する方法である。 しかし、現実に得られる画像データでは、a淡むら4、
〒性を表わす関数fが測定できないものや、未知なもの
、あるいは一旦関数fを測定しても、センサ2の劣化に
より実際に得られた画像の濃淡むらの特性と、事前に測
定した濃淡むらの特性が一致しでいないようなものが多
い。 このような画像に対しては、従来の濃淡むら補正方法で
は、濃淡むらを取り除くことが難しいという欠点がある
。 また、従来の、農淡むら補正方法では、画面の各画素ご
とに変換表を用意する必要があり、変換表を記憶するた
めメモリ容貝が大きくなる欠点もある。 (発明の目的) 本発明の目的は、これらの欠点を解消するため、fit
正閏数fが明確でないような謔(淡むらのある画像デー
タの濾淡歪を補正することができ、また記憶すべき変換
表の容置を減少して、高速高精度に瀦淡むらを補正でき
るディジタル画像歪補正装置を提供することにある。 〔発明の概要〕 本発明によるディジタル画像歪補正装置は、同一・被写
体′Iと濃淡むらの少ないセンサで撮った画像データ、
あるいはすでに濃淡むらの補正が行われた画像データを
利用するか、または常道のビジコン・カメラでは濃淡歪
の画面的変化がなめらかであるどいフ点を利用して、濃
淡也を補正することを特徴としている。 〔発明の実施例〕 先す、 1J”J淡むらの少ない画像データを利用する
装置(第1の装置C’i )について説明した後、次に
そのような画像データが利用できない場合に、陰面綿近
似の係数を人力することにより、有効に1農淡歪を11
11正する装置(第2の装置)について説明する。 第1の:(11置においては、a淡むらのある画像上で
、濃淡むら補正のために基準点P(Q、p)を設定し、
濃淡むらのない画像上で1)に対するのをめ pと(a
の周辺の画像のヒストグラムを一致させることにより、
画像の′a、淡むらを取り除く方法を用いる。。 いま、濃淡むらのある画像の画素強度をXとし、濃淡む
らのない画像の画素強度をyとすると、一般に、濃淡む
ら補正関数[は線形関数により−1分な精度で近似でき
るから、次式が成立する。 y −a :< + l+ ・ (5)ここで、a、1
〕は、画イg!の馬所(Q、p)の関数=1 Co、 
+ IIL b、(p、 、p)で表わすことが−(t
! ルuそこで、13周辺のヒス[−グラムをとり、平
均をl’X+分散をCX2とし、0周辺のヒス(−グラ
ムをとり、平均をμン2分散をσン2とず、flば、点
P、l二でa(1)、、p)、lJ(+1.p)は次の
各人でlj、えられる。 h (A I 1+)= t、t ン −t+ (R,
+ p) μ x −(7)」二組(6)、 (7)式
が導か、1℃るまての詳細を記述する。濃淡むらのある
画像の画素Xと4!淡むらのない画像の画素yについて
、平均強度μ7.μノを中心に分散をとると、第4図(
、)および第4図0〕)1;示す度数曲線となる。第4
図(b)におりる度数の平均値■Σ(y)を、Jξめる
と、前記(5)式よすy = a x + b ・・・
(5)IE (y) =li (ax+b) IE ((y py )2=E (((ax+b) (
aμx +b)) 2)’ −E (F12(X tL
×)2)=a2E (CX p×)2 ) −18)こ
こで、 σy 2−(y−pJ)2 ゜、2 =(、X’4 つ )2 ・・ (9ンである
から、(9)式を(8)式に代入すると、σy 2−n
 2σX2”’(10) となり、前記(6)式が導かれる。 また、q)41ン1(a)(b)より、平均強度1’J
+uxの間には、1)II記(5)式より次の関係があ
る。 71 、y : F1μ、 + b 、 −(11)し
たがって、」二式(11)式より前記(7)式が導かれ
る。 次の、前記(f;)(7)式で与えられる処理を、濃淡
歪のある画像で設定された基準点P(Q;、pJ)上で
すへでについて行って、a 、(0,□、PJ)+ b
(Q□、■・J)を計算し、これらの値を用いて画面、
]−の任意の点(Q、p)における関数”(QI P)
+b(Q’+p)を空間的な補間により君(算する。 このa(Q、pL b(Q、p)を用いて、濃淡むらの
ある画像の画素強度x(A、p)から、濃淡むらをとっ
た画素強度”(Q + P )を、次のHIEによりめ
る。 ”(Q、p)=a(D、、p)・x(11,p)+b(
A、p)(12) 第5図は、本発明の原理説明図であり、第6図は本発明
の第1実施例を示すティジタル画像の歪補正装置のブロ
ック図である。 第6図における濃淡むらのある画像23と#淡むらのな
い画像21は、それぞれ第5図の(a)と(b)に対応
する。第5図(a)は、nXmブロックに分割された濃
淡むらのある画像であり、ブロックごとに周辺のヒスト
ダラムを使用して、第5図(b)の画像に対応させるこ
とにより濃淡歪を補正する。 a淡むらのない画像21と濃淡むらのある画像23は、
処理装置25からの制御信号により、ティジタル画像人
力装置22、および24か゛ら入力される。 このとき、第5図(a)(b)に示すように、あらかじ
め′a淡むらのある画ftR上で設定されたj、4; 
f1点1)(Qi、p、>)の周辺ヒストグラムと1つ
に対応する濃淡むらのンJ′:い画像」二のQの周辺ヒ
ストグラムとを、i:lli像人力IRIに、画像入力
装置22,2/Iを用いて剖1つする。 濃淡も四)のある画像は、ヒスI−グラムかnI算され
ると同11,1に、rljij像記庶装置23に転送さ
れる。 画伶テータタ人力した後、処理装置25は画像人力装置
l”t +! 2+ 24か1)ヒストグラムを取り出
し、ノ、Cr(14点1’ (+!□、pJ)でのヒス
トグラ11の平均μm(i、j)と、ヒスI−グラムの
分散σX”NIJ)と、それに7・1応するQの平均μ
:V(++JL分散’7:V2N+J)を81算し、さ
らにa(i、jL b(i、j)を、11tつする。 次に、基準点の周辺のヒストグラムについて、第7図、
第8図により説明する。 第7図に示すように、画素強度Xのむらのある画像中で
基iiQ点l〕をとると、場所によっては緩慢な変化の
関数となり、また他の場所によっては急激な変化の関数
となるので、近似をとるためには、緩慢な変化のときP
の領域を広く、急激な変化のとき狭く設定する1、領域
の格子点にお()る関数y(y+〜y、4)をa!算し
て他の点は、空間的な補間を行う。 !/ l :Il l X+” l + / 2 ” 
、12 χ→−b2y :l = :r 3:< +)
+ :、I 、 y I= a 、1x + b 、1
(+3) a(i、j)、1)(i、j)は、格子点上の係% D
 H。 C2+ ’l :J + ’] −1+ b 1 r 
b 2 + ” :4 + bdを用いて、陰面線関数
として石質する。 各8,1)を画(Q<の濃淡レベルで立体的に記述すれ
は、第8図(+1)(1))に示す形状となる。 a(1)、、p)=ag +al 1−1− C2P+
a3 L・Pl、+(12,P)=βυ+β、L十β2
P+β3L・、P・・・(14) ここで、(L、P)は、ブロック内相対位置である。 十記(I4)式で、”+ pの陰面線関数、QIPの多
項式関数どして割算することにより、係数a(A+ r
’)+ h(’l+ I))をめ、こオしらを濃淡むら
補正装置Vj27に転送し、カシ)画像記憶装置2Gの
画(etテデーを(11!!、淡むら補正装置27に転
送することによって、補正画像28を1′)るごとがで
きる。 第1]図(J、関数a CQ + P L b (Q 
+ I’ ):/l’陪直陰面数である場合のニガ1(
談むらFl’l!正シ七置のブ装4ツク1ソ[である。 処理装Jli(z25かt)送I)れてくる信r4II
 Oは1.1(Qr P)l l) C1)、 、 p
)の係数が、あに)いは、農1び5む1ら浦iLの開り
ri <、1号てあり、いづ’;hの信じであるかを判
定回路29で判定する。。 テークかa<ay pL b(QI P)の係数であ5
1(ば、陰面線モ六゛ル・テーブル33にl1lJ′i
次データを転送し、テーブル33の内容%、 j、ff
1j淡むら油止ができるように七)1−する、1 a淡むら補正の開始(rJ号であ4(げ、カウンタ:3
0に信号が送られる。カウンタ3oけ、傭シ3を受け取
ると、カウンタの内容を0にクリアし1、順次カウンタ
内容をカラン1−・アップするが、カウンタ内容後1カ
ウン1−だけアップする間に、諾淡むら補止装賀は、次
の動作を行う。 先す、カウンタの内容Kが陰面線モテル・テーブル′3
:3のア1−レス81tつ回路32に送られるので、ア
ドレス81算回路32は次式の関係を使、ノて画作上の
画素位置(jl 、 p)および該当するブロック(i
。 j)を石質する。 ■(=(ρ−] )X P L 十p −(15)ここ
で、Pl−ば1ラインの画素数である。 さらに、」ユ記(15)式し二より絹質した(Q、P)
から、その点がどのブロック(i、j)に含まれるかを
β1算する。 屁し1て、アドレス削9回路32は、該当するテーブル
の71〜レスADと、プロッタ内の相対ライン1、を次
式によりH1算する。 Δu=(i −1)Xm+j −116)ここで、mは
イ黄方向ブロック数である。 L = Q、 −(i−1) X L B−1・・・(
17)ここで、t、 Bはブロックのライン数である。 H1算した結果、テーブル・ア1−レスΔ、、 1.i
陰面線モデル・テーブル33のアドレス・レジスタ31
に、また、相+4ラインl、ば;寅算)1(3・]よン
よび:35に送られる。 ア1−レス・レジスタ31の内容にしたかって陰面餘モ
テル(y)係数(t1+ (’21CL3 + (C4
+ (’L +(’2+ β・2.θ、か陰面線モデル
・テーブル33から読み出され、;寅算器34.3(に
転送されろ。 第1 CI 1′!l Ink、第9図における演勢5
器34 、 :35の構成1ift eある。 r寅%゛回j、18 、I 11では、cl+a2Lま
た(まtノ1−1−β2Lのl寅算を、またC寅算回路
45て1ま、C11J十α、1■7または03+βIL
、の演初″を行う。演算回路44.45で31算された
値は、係数+1に関ずろα14−+z21.+ α3+
α1■−7と、係数I)に関するβ1 (−〇、+−+
 β3+β、1■−であり、各4・が’t’r’(W)
)跨36 + 37に送られる。 検出器/l 3 It、、アドレス言1算凹略32の出
力である陰面線モデル・テーブル33のア1−レスAD
と、相対うイン番53. Lが前回のそれと同一か否か
を検出【11、いずれか一方でも前回と異なっていれば
、演算i!1;3 Fl 1 37に信号−′制°′を
送り、異なっていかければ(0号” o ”を送る。検
出8343の目的は、1ライン内のブロックの変り目の
検出である。 演21 冊3 C;はa(Q+p)を;:1算し、演算
器371:):a(Q、p)を+tl算する。 第1J図は、第!11!、jlにおける演算!iB3 
に 、37の内部構成図である。 検出!if;43からの信号4Gが1“であれば、選択
回路/+ 91J、演算+iH3’I (あるいは35
)の出力であるαl−1−ri21、(あるいはβ1+
β2【−)を選択し、選択回路50はOを選択する。 このとき、加算H’rf 52は、ブロック内の(11
対ピクセルPの値がF)=0である場合に対応する係数
を出力する1、また、選択回路51(J、ゲインに対応
するレジスタ53に(χ3+tr 41− (j5るい
はβ3+β11.)を格納する。 1句、′1冒(侶4:玉からの信号4〔1が″(1”の
場合には、選V(回路49と50は、加算))踵52の
l IiF、I i)口の出力内容を格納したレジスタ
5ノ1の内容と、ゲインを表わすレジスタ5:3の内容
を選択し、これらを加算)!:+ 52に入力する。加
算器52の出力r) U Tが請求める係数a(あるい
はb)である。 111り1ビクセ/lN〕は整数であるから、その回数
だけ加勢゛冊52−r加算を繰り返すことしこより、a
 (Q、 、p)=(+z’o 4−α1L)−1−1
)(ζλ2→−αr+L)の値がめら1しる。 第91!1に才イいて、?寅算器3G、37の出力を来
算器38、加い器39に転送し、画像記憶μ!1置26
からあしす、出された濃淡むらのある画像の画イ((〕
強度Xを用いて、平算器38でaメ、加算器31)でa
 x −1−bのrtij ’t3’、を行うことによ
り、次式か−1くめられる。 仝= a x 十b −(12) ’ 前述のように、仝はi11τ1−さ、bた画素強度であ
る7このJ、うにして、第1の装置においては、#淡も
・らのある画像の強度を、l農法むらのない基7j9の
画像の強度に簡単に重ね合わせることにより、濃淡むら
を補正できる。 次に、第2の装置について、説明する。 第2の装置は、第1の装置の場合と異なり、濃淡むらの
ない画像が得られない場合に、陰面線のモテルによる係
数をあらかじめ記憶しておくことにより、メモリ8凧を
増大させることなく、高精度にl!A淡むらを補正する
ものである。 第12図は、本発明の第2実施例を示すデ2rジタル画
像歪補正装置のブロック図である。 被写体1からの光をセンサ2により電圧に変換し、補正
f」アナログ・ディジタル変換器63を介して、記憶装
置4に記憶する。補正伺アナログ・ディジタル変換器6
3は、第13図に示すような変換特性57を有し、人力
装置64より陰面線近似の係数を人力して、アナログ・
ディジタル変換と同時に濃淡むら補正を行う。第13図
の変換特性57は、前記(1,)、(2)、(3)、(
4)式を用いて、u=(f(Z)) l111x = (g(y):] ・・(18) として−ケえられる。 ここで、前述のように、通常のビジコン・センサでは、
[(Z)は線形近似可能であることにもとづいて、rr
 (y )も線形近似可能どなる。 いま、上記(18)式のg(y)を線形近似したものを
、次式で表わすものとする。 g (y )= a y 十b ・(19)濃淡むらを
持つ人力画像55の各画素について、上1!g(18)
式による濃淡強度の変換を行った後、記憶装置4に読み
込まれた画像の濃淡は、5Gに示すように均一となる。 次に、センサ出力信号の標本化方式について説明する。 センサ感光面の電何として呪られだ画像データは、第1
4図に示すように、電子ビーl、によって画面の左]二
から順次lラインずつ合計Nラインが読み出される。読
み出されたセンサ出力(i t’359は、第15図に
示すように、時刻とともに変動する電圧となる。これを
Δt1の周1期で標本化し、1ライン当りM個の画素デ
ータを得る。ここで。 ライン間の画像データのない時間をΔt2とする。 第16図は、第12図における補正(ツアナログ・ディ
ジタル変換器の構成図である。 変換特性発生装置60は、画面の各画素について前記(
I9)式の係数a、bを作成するものであり、サンプラ
65はタイミング発生回路67からのタイミンク信号(
Δt1ごとのタイミング信号、ただしライン間のΔL2
には信号はない)により、八t1ごとにセン勺2の出力
信号を標本化する。 71−レス発生口’J868は、タイミング発生回路6
7からのタイミング信号により、標本化されている画素
に対応するア1−レス(ライン番号Qど、ライン上の画
素番号l〕)を変換特性発生装rti、60に送る1、
変換特性発生装置60は、上記アドレスに一対応する(
19)式の係数a、1)をめて変換回路ら6に送る。 第17図は、第16図にお(づる変換回路6Gの内部構
成図である。 変換特性発生装置60から送らAした係数のうち、(1
9)式の1〕は、ディジタル・アナログ変換’A’+’
r 661でアナログ信号化されて演算増幅器663に
、また、(19)式のdはディジタル・アナログ変換器
662でアナログ信号化されて演算増幅?!B 664
に、それぞ」し送られる。その後、サンプラ65で標本
化されたll!II像テータはデー勢増幅器664で8
倍され−Cから、演算増幅器663でbだけ加算される
。その後、第3図(I))と同一の特性を持つアナ1」
グ・ディジタル変換器665を通ることにより、i?i
l記(18)式に相当する変換を受り、記憶装置4に送
出される。このようにして、タイミング信号のI’ik
 1f(I+により、各画素ごとに濃淡歪h1i正と標
本化・−5tf−化が順次、高速に行われる。 第18図および第19図は、第16図に才言)る変4J
i% 45性発生装置60の変換特性作成原理を示すI
S〈1で、らる、。 第133図のように、画面を縦横それぞれr】個。 111個のゾI」7りに分割する。分割されたプロッタ
1うの各洛「・点で、前述した方法により(19)式の
係数d、1〕がめられているものとする。 このとき、1つのブロックB内の画素70(座標を(Q
ap)とする)のa、bを、次のようにしてめる。。 第15〕図に示すように、1ブロツクの大きさをり、1
3ラインXPT3画素とする。このとき、画素70の対
応するブロック番号(i、j)は、次式でめられる。 また、ブロノグ内の相対座標(u、■)は、次式%式% このとき、(+、1.V)におりるaとbを次式で表わ
す。 a (u、v)= tto+tz 1 u+a2 v+
a3 uv (22)1) (IIIV) =βo−1
−βI IJ+β2v−1−β:(uv (23)この
陪直線近似は、分割ブロック数m、+1を1・分多く(
JO個以−に)設定すれば、通16のビジコン・カメラ
にj@川用能である。。 ここで、前記(22)、 (23)式の係数α□、β1
(1=O,i、’2.3)は、ブロック分割の格子点で
測定されたa、l)から次のようにめられる。 す1.1′、わイノ、第19図の各格子点のa□、b4
(i= I) 、 I 、 二! 、 3 )から、次
式により得・られる。 +zg 引鳥 0 1 第20図は、第16図にお番プる変換特性発生装置60
の内部構成図である。 人力;iJ、装置64から、プロy’/I3(]、+ 
1.)l・・・・1’3 (+1.■)のそれ−どれに
つし)で、」二1己(24) 、 (25)式にしたか
っCあらかしめめられた係数It i 。 β4 (i =O,’I 、2.3)がメモリ602に
入力され、記1意さ、1シる9名画ヌそ位置に対するa
、hは、次のようにして:1(められる。 アドレス光生回路68から送らJlてきた画素座標(Q
++・うから、アドレス変換回路601は、前記(20
)式により1咳当ブロック計号(i、j)を、また前記
(2+):iK:に上りブロック内相刺座漂(11,V
)をa)い゛する5、ブロック番号(i+ j’)にツ
、1応するrni記(:!2)式のli V a ; 
(i 二〇 、 ] 、 2 、3)と前記(2:l 
)式の係数β1(1−0,1,2,3)かメモリ602
か1″、読み出さ、IL、ぞれぞJし演算回踏603゜
60.1に送らALる、 1寅算回路603,604は、ア(−1ノス変換回路6
01の出力である相対座標(u、■)ど前記係数Cχ1
.β1を用いて、それぞれ(22)、 (23)式にし
たがい係数8および係数すを算出する。 この算出された結果が、変換回路66に送られる。 このように、第2の装置によお、は、あらかしめめられ
た陪直線近似の係数αユ、β1を記憶しておくことによ
り、ビジコン・カメラ等で撮影された画像の画素ごとの
譜淡歪補11.とディジタル化 −が、同時に高速に行
われろ、。 〔発明の効果〕 以」;説明したように、本発明によ、ILは、濃&J2
むらの少ない画像を利用して、あるいばあらかじめ陰面
線通((jのための係数ti−割算して記憶、すること
に−より、!Ij!j淡むらを含む画像データの濃淡歪
を高精度に油止でさ、また、変換表等を記1.αしない
ので、メモリ容量を増加することがなく、テイシタル信
号に変換して高速にディジタルil[7,機に入力する
ことにより、高速度の画像f’1M’ ”Jj ’J:
;よび処理が可能とン多〕る。
[Detailed Description of the Invention] [Field of Application of the Invention] The present invention relates to a digital image distortion correction device. This invention relates to a device that corrects d (light distortion) of an input image by It is necessary to convert the analog signal shown into a digital signal and input it to the d1 calculator. Fig. 1 is a schematic diagram of a conventional image input device. is a sensor, 3 is an analog-to-digital converter, and 4 is a storage device.The radiation intensity of the object L is formed by the sensor 2 as a two-dimensional intensity distribution as shown in the image 8, and the intensity at each point is The analog signal shown in FIG. As shown in FIG. 2, a light ray 11 from an object focuses an image as an electric charge on a photosensitive surface 13 of a vidicon or the like in the focal plane by an optical system 12 such as a lens. 444
B is read out, converted to voltage, and sent out via signal line 1G. By the way, the image captured by the input device shown in FIG. 1 contains uneven shading due to the sensitivity of each point on the photosensitive surface 13 of the sensor 2. In other words, even if a light source with a uniform brightness distribution is used for the subject (1), the shading intensity distribution of the sensor output image 8 or the shading intensity distribution of the image 9 read into the It is not uniform on one surface, but causes unevenness. Therefore, in order to use the image data for accurate analysis and processing, it is necessary to remove this unevenness in shading. In order to do this, as shown in Figure 1, the conversion device 6, storage device 7 and conversion table are
A correction device consisting of an input device 5 is connected downstream of the input device. The conversion device 6 converts each image data obtained in the storage device 4 (l('1), that is, the intensity of the sampled image unit, into an intensity with unevenness removed according to the conversion table.2
The conversion table is used to make the a-dark intensity distribution of the image uniform for the subject 1 of a light source with uniform brightness, and the pixel data converted according to this conversion table is stored in the storage device 7. be damaged. This conversion table is created in the order shown in Fig. 3 (a>(b) (c). For example, the unit is cd), and the input radiation intensity is
Let y be the output cuff 6 pressure appearing on the sensor output signal line 16 that slopes to
The relationship (x, 'y') is obtained at the position on the sensor surface corresponding to b. That is, at each point 17 in FIG. 3('a), the radiation intensity X is varied from 0 to x3. This is S-constant data corresponding to each pixel of one light source when measurement is performed using a light source (<-). The 31st "I (b) is a model of the characteristics of an analog-to-digital converter. The range of sensor output voltage 1/Jl (from 0) is the digital data after conversion when it changes from 0 to If the number of levels is 1-, then the analog-to-digital conversion for the sensor 11 voltage y!I
:+'s output Z is given by the following equation. Z = (1,, , the variable table that can be obtained is as follows: For each point 17 in FIG. -3...(2)
Point +9 indicating the relationship between digital converter output 2 and correction strength U. and a curve 20 connecting points 19. Function 2 obtained by interpolating each point 19 in Figure 3(c)
0 can be calculated using the following equation. u = f(7,) (:s) Also, the conversion table is the level of variable Z, O, I...l, -
1 and is expressed by the linear equation ``C''. 10 minutes
′: Can it be approximated by a linear function with accuracy?1. In this way, for each pixel of the image, the following 1.
By doing °6Ii"r"A+, the desired conversion table is obtained 4
I will. Conversion table? In the conventional method of correcting the 11l pale unevenness of a human image using r, the dull unevenness characteristic of the sensor 2 can be completely determined by U+, and )) This method is based on the premise that the function 1 can be completely modeled. However, in the image data that is actually obtained,
〒The function f that expresses the density cannot be measured or is unknown, or even if the function f is measured once, due to the deterioration of the sensor 2, the characteristics of the uneven shading of the image actually obtained and the uneven shading measured in advance In many cases, the characteristics of the two do not match. For such images, conventional shading unevenness correction methods have the disadvantage that it is difficult to remove shading unevenness. Furthermore, in the conventional method for correcting unevenness in agriculture, it is necessary to prepare a conversion table for each pixel on the screen, and there is also the drawback that the memory capacity becomes large because the conversion table is stored. (Object of the invention) The object of the invention is to solve these drawbacks by
It is possible to correct the filtration distortion of image data where the positive leap number f is not clear (it is possible to correct the filtration distortion of uneven image data, and to reduce the storage space of the conversion table to be memorized, it is possible to correct the filtration distortion of uneven image data at high speed and with high precision. An object of the present invention is to provide a digital image distortion correction device capable of correcting image distortion.
Alternatively, the shading is corrected by using image data that has already been corrected for shading, or by using the point where the shading distortion changes smoothly on a regular videocon camera. It is said that [Embodiments of the Invention] First, a device (first device C'i) that uses image data with little 1J"J light unevenness will be explained, and then next, when such image data is not available, By manually calculating the coefficients of the cotton approximation, 1 agricultural distortion can be effectively reduced by 11
11. The device for correcting (second device) will be explained. First: (At position 11, a reference point P (Q, p) is set on the uneven image to correct uneven shading,
For 1) on an image without uneven shading, p and (a
By matching the histograms of images around
A method is used to remove 'a' and light unevenness in the image. . Now, if the pixel intensity of an image with uneven shading is X, and the pixel intensity of an image without uneven shading is y, then in general, the uneven shading correction function [can be approximated by a linear function with an accuracy of -1 minute, so the following equation holds true. y − a :< + l+ ・ (5) Here, a, 1
] is a picture! The function of horse field (Q, p) = 1 Co,
+ IIL b, (p, , p) can be expressed as −(t
! Then, take the hiss [-gram around 13, set the average as l' At points P and l2, a(1), , p), lJ(+1.p) can be obtained by lj for each of the following people. h (A I 1+) = t, t -t+ (R,
+ p)μ Pixels X and 4 of an image with uneven shading! For pixel y of an image without light unevenness, the average intensity μ7. Figure 4 (
, ) and FIG. 4 0]) 1; Fourth
By subtracting Jξ from the average value ■Σ(y) of the frequency shown in figure (b), the above formula (5) is obtained: y = a x + b...
(5) IE (y) =li (ax+b) IE ((y py )2=E (((ax+b) (
aμx +b)) 2)' -E (F12(X tL
×)2)=a2E (CX p×)2 ) −18) Here, σy 2−(y−pJ)2 ゜, 2 = (, Substituting equation (9) into equation (8), σy 2-n
2σX2"' (10), and the above equation (6) is derived. Also, from q)
+ux has the following relationship according to 1) Expression (5) of II. 71, y: F1μ, + b, −(11) Therefore, the above equation (7) is derived from the equation (11). Next, the process given by the above formula (f;) (7) is performed on the reference point P (Q;, pJ) set in the image with grayscale distortion, and a , (0, □ , PJ) + b
Calculate (Q□, ■・J) and use these values to display the screen,
]− at any point (Q, p)”(QI P)
+b(Q'+p) is calculated by spatial interpolation. Using this a(Q, pL b(Q, p), from the pixel intensity x(A, p) of an image with uneven shading, The uneven pixel intensity "(Q + P) is determined by the following HIE. "(Q, p)=a(D,,p)・x(11,p)+b(
A, p) (12) FIG. 5 is an explanatory diagram of the principle of the present invention, and FIG. 6 is a block diagram of a digital image distortion correction apparatus showing a first embodiment of the present invention. An image 23 with uneven shading and an image 21 without uneven shading in FIG. 6 correspond to (a) and (b) in FIG. 5, respectively. Figure 5 (a) is an image with uneven shading divided into nXm blocks, and the shading distortion is corrected by using the surrounding histodarum for each block and making it correspond to the image in Figure 5 (b). do. aThe image 21 with no unevenness and the image 23 with unevenness are as follows.
A control signal from the processing device 25 is inputted from the digital image human power devices 22 and 24 . At this time, as shown in FIGS. 5(a) and 5(b), j, 4;
The peripheral histogram of point f1 (point 1) (Qi, p, >) and the corresponding peripheral histogram of Q with uneven shading are input into i:lli image IRI using an image input device. Perform one autopsy using 22,2/I. The image with shading (4) is transferred to the rljij image recording device 23 after the Hiss I-gram is calculated by nI (11,1). After inputting the image data, the processing unit 25 extracts the histogram from the image input unit l''t +! (i, j), the variance of the His I-gram σX”NIJ), and the average μ of Q corresponding to 7.1
:V(++JL variance '7:V2N+J) is calculated by 81, and a(i, jL b(i, j) is multiplied by 11t. Next, regarding the histogram around the reference point, FIG.
This will be explained with reference to FIG. As shown in Fig. 7, if we take base iiQ point l] in an image with uneven pixel intensity , in order to take an approximation, when the change is slow, P
Set the region to be wide and narrow when there is a sudden change 1. Set the function y(y+~y,4) at the grid point of the region to a! and other points perform spatial interpolation. ! / l: Il l X+" l + / 2"
, 12 χ→−b2y :l = :r 3:< +)
+:, I, y I= a, 1x + b, 1
(+3) a(i, j), 1) (i, j) is the coefficient % D on the grid point
H. C2+ 'l :J + '] -1+ b 1 r
Using b 2 + ”:4 + bd, stone quality is determined as a hidden line function. Each 8, 1) can be described three-dimensionally at a density level of Q< as shown in Figure 8 (+1) (1). ) has the shape shown in a(1),,p)=ag +al 1-1- C2P+
a3 L・Pl, + (12, P) = βυ + β, L + β2
P+β3L·, P (14) Here, (L, P) are relative positions within the block. In Equation (I4), the coefficient a(A+ r
')+h('l+I)), transfer these to the unevenness correction device Vj27, and transfer the image (et) in the image storage device 2G (11!!) to the unevenness correction device 27. By doing so, the corrected image 28 can be obtained as 1').
+ I'): Niga 1 (
Danmura Fl'l! There are 4 pieces in 1 position on the front side. Processing device Jli (z25 or t) sent I) Incoming message r4II
O is 1.1 (Qr P) l l) C1), , p
The determination circuit 29 determines whether the coefficient of ) is the belief that the opening ri of agricultural 1bi5mu1rauraiL<, 1, and izu';h. . The coefficient of take or a<ay pL b(QI P) is 5
1 (for example, l1lJ′i in the hidden surface line model table 33
Transfer the next data and set the contents of table 33 %, j, ff
1j To be able to stop oil from light unevenness 7) 1-, 1a start light unevenness correction (rJ No. 4(ge, counter: 3)
A signal is sent to 0. When the counter 3 is received, the contents of the counter are cleared to 0 and 1, and the contents of the counter are sequentially increased by 1-. Stop-souga performs the following operations. First, the content K of the counter is Hidden Line Motel Table '3
:3 address 81t is sent to the circuit 32, so the address 81 arithmetic circuit 32 calculates the pixel position (jl, p) on the picture and the corresponding block (i
. j) to stone. ■(=(ρ-)) )
From this, calculate β1 in which block (i, j) that point is included. Then, the address cutting 9 circuit 32 calculates H1 from 71 to address AD of the corresponding table and the relative line 1 in the plotter using the following equation. Δu=(i −1)Xm+j −116) Here, m is the number of blocks in the yellow direction. L = Q, -(i-1) X L B-1...(
17) Here, t and B are the number of lines in the block. As a result of H1 calculation, table address Δ, 1. i
Address register 31 of hidden line model table 33
, it is also sent to the phase + 4 line l, ba; tari calculation) 1 (3・] call and : 35. If you want to make the content of the a1-res register 31, the hidden side model (y) coefficient (t1 + ( '21CL3 + (C4
+ ('L +('2+ β・2.θ, read from the hidden surface line model table 33 and transferred to; Performance 5 in the figure
There are 1ift e configurations of vessels 34 and 35. At 18, I 11, cl+a2L or (1-1-β2L), and C11J0α, 17 or 03+βIL.
, is performed.The value calculated by 31 in the arithmetic circuit 44.45 is α14−+z21.+α3+ regardless of the coefficient +1.
α1■-7 and β1 (-〇, +-+
β3+β, 1■-, and each 4 is 't'r'(W)
) is sent to straddle 36 + 37. Detector/l 3 It, address AD of the hidden surface line model table 33 which is the output of the address word 1 address AD
And the opposite number 53. Detect whether L is the same as the previous one [11. If either one is different from the previous one, operation i! 1; 3 Send the signal -'control' to Fl 1 37, and if it differs, send (No. 0 "o"). The purpose of the detection 8343 is to detect the change of blocks within one line. Performance 21 Book 3 C; increments a(Q+p) by ;:1, and arithmetic unit 371 :):a(Q,p) is multiplied by +tl. Figure 1J is No.! 11! , operation on jl! iB3
37 is an internal configuration diagram. detection! If signal 4G from 43 is 1", select circuit/+91J, operation+iH3'I (or 35
), which is the output of αl−1−ri21, (or β1+
β2[-) is selected, and the selection circuit 50 selects O. At this time, the addition H'rf 52 is (11
1, which outputs the corresponding coefficient when the value of the pixel P is F)=0, and stores (χ3+tr 41- (j5 or β3+β11.) in the register 53 corresponding to the selection circuit 51 (J, gain). 1st phrase, '1 blasphemy (Mate 4: Signal 4 from ball 4 [1 is '' (if 1'', selection V (circuits 49 and 50 are addition)) heel 52 l IiF, I i) The contents of register 5:1, which stores the output contents of the mouth, and the contents of register 5:3, which represents the gain, are selected and added and inputted to 52.The output of adder 52 is This is the claimable coefficient a (or b). Since 111/lN] is an integer, by repeating the addition for that number of times, a
(Q, ,p)=(+z'o 4-α1L)-1-1
)(ζλ2→−αr+L) is 1. The 91st! The first one is so talented? The outputs of the input calculators 3G and 37 are transferred to the input calculator 38 and adder 39, and the image storage μ! 1 place 26
Karaashisu, the image with uneven shading (()
Using the strength
By performing rtij 't3' of x -1-b, the following equation can be subtracted by -1. = a x 1b - (12) 'As mentioned above, is the pixel intensity of i11τ1- and b. The uneven density can be corrected by simply superimposing the intensity of the image on the intensity of the image of the base 7j9, which has no unevenness in l farming method. Next, the second device will be explained. Unlike the first device, the second device can store coefficients based on hidden surface lines in advance without increasing the memory capacity when images without uneven shading cannot be obtained. , with high precision! A: It corrects light unevenness. FIG. 12 is a block diagram of a digital 2r image distortion correction device showing a second embodiment of the present invention. The light from the subject 1 is converted into voltage by the sensor 2 and stored in the storage device 4 via the analog-to-digital converter 63. Correction analog/digital converter 6
3 has a conversion characteristic 57 as shown in FIG.
Performs density unevenness correction at the same time as digital conversion. The conversion characteristics 57 in FIG. 13 are as follows: (1,), (2), (3), (
Using equation 4), it can be calculated as: u=(f(Z)) l111x = (g(y):]...(18) Here, as mentioned above, in a normal vidicon sensor,
[Based on the fact that (Z) can be linearly approximated, rr
(y) can also be linearly approximated. Now, suppose that the linear approximation of g(y) in the above equation (18) is expressed by the following equation. g (y) = a y 1 b ・(19) For each pixel of the human image 55 with uneven shading, the upper 1! g(18)
After converting the shading intensity using the formula, the shading of the image read into the storage device 4 becomes uniform as shown in 5G. Next, a method of sampling the sensor output signal will be explained. The image data that is cursed as an electric shock on the sensor photosensitive surface is the first
As shown in FIG. 4, a total of N lines are read out sequentially from the left side of the screen by l lines by electronic beer l. The read sensor output (i t'359 is a voltage that fluctuates with time as shown in Fig. 15. This is sampled in one cycle of Δt1 to obtain M pixel data per line. .Here, the time during which there is no image data between lines is assumed to be Δt2. FIG. 16 is a block diagram of the correction (Tsanalog-to-digital converter) in FIG. Regarding the above (
The sampler 65 generates the coefficients a and b of formula I9), and the sampler 65 receives the timing signal (
Timing signal every Δt1, but ΔL2 between lines
(there is no signal at ), the output signal of sensor 2 is sampled every 8 t1. 71-Response generation port'J868 is the timing generation circuit 6
7 sends an address (line number Q etc., pixel number l on the line) corresponding to the pixel being sampled to the conversion characteristic generator rti, 60;
The conversion characteristic generator 60 corresponds to the above address (
19) The coefficient a, 1) of equation 1 is calculated and sent to the conversion circuit 6. FIG. 17 is an internal configuration diagram of the conversion circuit 6G that is similar to that shown in FIG.
9) Equation 1] is digital-to-analog conversion 'A'+'
r is converted into an analog signal by 661 and sent to an operational amplifier 663, and d in equation (19) is converted to an analog signal by a digital-to-analog converter 662 and sent to an operational amplifier? ! B 664
``That's it'' and it's sent to you. After that, ll! was sampled with sampler 65. II image data is 8 at the data amplifier 664.
From the multiplied -C, an operational amplifier 663 adds b. After that, Ana 1 with the same characteristics as Figure 3 (I))
i? i
It undergoes a conversion corresponding to equation (18) and is sent to the storage device 4. In this way, the timing signal I'ik
1f (By I+, the gradation distortion h1i positive and the sampling/-5tf- conversion are sequentially performed for each pixel at high speed. Figures 18 and 19 are similar to Figure 16)
I illustrating the principle of creating the conversion characteristics of the i% 45 property generator 60
S〈1, raru. As shown in Figure 133, the screen is divided into r] pieces vertically and horizontally. Divided into 111 parts. It is assumed that the coefficient d, 1] of equation (19) is determined by the method described above at each point of the divided plotter 1. At this time, the pixel 70 (coordinates of (Q
ap)), a and b are determined as follows. . As shown in Figure 15, the size of one block is
3 lines XPT 3 pixels. At this time, the block number (i, j) to which the pixel 70 corresponds is determined by the following equation. Also, the relative coordinates (u, ■) in the blog are expressed by the following formula % Formula % At this time, a and b falling at (+, 1.V) are expressed by the following formula. a (u, v) = tto+tz 1 u+a2 v+
a3 uv (22)1) (IIIV) =βo-1
-βI IJ+β2v-1-β: (uv (23) This parallel line approximation increases the number of divided blocks m, +1 by 1・(
If you set JO number or more), you can use j@kawa function for 16 vidicon cameras. . Here, the coefficients α□, β1 of equations (22) and (23) above are
(1=O, i, '2.3) can be calculated from a, l) measured at the grid points of block division as follows. 1.1', Waino, a□, b4 of each grid point in Figure 19
(i=I), I, two! , 3), it can be obtained by the following equation. +zg Hikitori 0 1 Figure 20 shows the conversion characteristic generator 60 shown in Figure 16.
FIG. Human power; iJ, from device 64, proy'/I3(], +
1. ) l...1'3 (+1. ■), (24) and (25) are the roughly specified coefficients It i . β4 (i = O, 'I, 2.3) is input to the memory 602, and a for the position of the 9th famous painting is
, h are determined as follows: 1 (set). The pixel coordinates (Q
++・Ukara, the address conversion circuit 601 has the above (20
) formula, the 1 cough block measurement number (i, j) is obtained, and the above (2+):iK:
) is a) 5, li V a of the rni notation (:!2) formula corresponding to the block number (i + j');
(i20, ], 2, 3) and the above (2:l
) coefficient β1 (1-0, 1, 2, 3) or memory 602
or 1'', read out, IL, and send to the arithmetic circuits 603 and 60.1, respectively.
The relative coordinate (u, ■) which is the output of 01 is the coefficient Cχ1
.. Using β1, coefficient 8 and coefficient S are calculated according to equations (22) and (23), respectively. This calculated result is sent to the conversion circuit 66. In this way, by storing the predetermined bilinear approximation coefficients α and β1 in the second device, it is possible to record the information for each pixel of an image taken with a vidicon camera, etc. Light distortion compensation 11. and digitization at high speed. [Effects of the Invention] As explained, according to the present invention, IL is
By using an image with less unevenness, or by dividing and storing the coefficient ti for ((j) in advance, !Ij! 1.α is not recorded in the conversion table, etc., so there is no need to increase the memory capacity, and by converting it into a digital signal and inputting it into the digital signal at high speed. , high-speed image f'1M' ``Jj 'J:
; many can be read and processed.

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

第1図は従来の画像入力装置の概m3図、第2図は第1
図のセンサの原理図、第3図は第1図における変換表の
作成順序説明図、第4図は本発明で用いZIiijij
像の内索強度の分散曲線図、第5図は本発明の原理説明
図、第6図は本発明の第1の実施例を示すディジタル画
像の歪補正装置のプロッタ(ν1、第7図t、1. J
iい′(Q点周辺のヒストグラ11の説明図、第8図(
:]係数B+、bの陰面線モデルの説明(シ1、第9図
は第+5図における濃淡むl)補正装置のブロック図、
第101°4(J第]昭にお()るπ(算器(34,:
う、5ンの(1守成図、第11.Ji4は第9図におけ
る演−、K(冊(3G。 37)の4vl成図、第]2図は本発明の第2実施例を
ボノーデrシクル両像ゆMi ME装置のプロッタ図5
第1;(図は第12図におりる変換1−Ii性を示す図
、第1・11シ1.L)″き15図はそjしぞれセンザ
出力信シ3の院本化カ一式0) 、iil 明I J、
第161”l ハ第1.2 r31 t:= t+ l
’J ル111i jl−イ:Iアナロク・ティジタル
変換器の構成図、第17し1は第1〔5図にお(づる変
換回路の内部構成図、第1814:”I 、fsよび第
19図は第1iさIにおける変換1・1゛性発生装置の
変換特性作成原理を示す図、第20図は第16図にJハ
づる変換特性発生装置の内部構成図である。 21:、農淡むらのない画像、23:八淡むらのある1
iIii像、22.2/I:画fg!人力装置、25:
処理装置、:g に :画像記イ、α装置、27:譜萩
むらn11正1!’j、 Vfj、28 ’ tili
iIIE画像、29 : 判5L Fffl m、:(
0:カウンタ、3]ニアドレス・レジスタ、32ニアト
レスaf′t”):回路、3:3=陪直線モデル・チー
フル、34へ37:演算H::、38:乗算fli; 
、 39 :加t−1器、44 、 lI 5 :演算
回路、’l (]、50.51 :選択回路、52:加
算器、53.5/l:レジスタ、1; :1 : 11
1i正イ・Iアチーロク・ディジタル変換)1K、64
:人力装置r’i、57:変換特性、58:走査、59
:センリ出力信号、らO:変換特性発生装置、65:す
°ンブラ、06:変換回路、67:タイミング亮生回シ
1()、+; )(ニア1−1ノス発生回路、6 (3
1、662:デCシタル・アナログ変換)(:+、66
3.[36/1:演λつ増幅)1:4.665:アナロ
グ・ディジタル変換)):(,7〔J:画素、71〜7
4:格子点、601:アト1フ六変換回路、602:メ
モリ、603 。 6 (l lI ニル(tつ回路。 第1図 100 101 第2図 第3図 (a)(b)(c) 第4図 メ1xXllyX 第5図 mブロック 第6図 ゐ)7図 り 第8 図 (a) (L、 P、 ) (L2P、) (b) 第9図 aa:<+b 第 1(]χ1 α1+α2L αけαlI。 ゐr 1 1 ::!1 第1214 第 13:χ; セ/す出力市IE(yf 第14 ’>j へ′!■5;4 八5Z1t 第161:4 牛 第171′(1 az 19図 第201χ1 第1頁の続き 0発 明 者 坪 井 晃 @発明者 井原 広− [相]発明者 久 保 裕 川崎市麻生区王禅寺109幡地 株式会社日立製作所シ
ステ日立重大みか町五丁目2番1号 株式会社日立製作
所大みか工場内
Figure 1 is a schematic diagram of a conventional image input device, and Figure 2 is a diagram of a conventional image input device.
Figure 3 is a diagram explaining the creation order of the conversion table in Figure 1, Figure 4 is a diagram showing the principle of the sensor used in the present invention.
A dispersion curve diagram of the inner cable strength of the image, FIG. 5 is a diagram explaining the principle of the present invention, and FIG. 6 is a plotter (ν1, , 1.J
i' (Explanatory diagram of histogram 11 around point Q, Figure 8 (
:] Explanation of the hidden surface line model of the coefficient B+, b (C1, Figure 9 is the shade of Figure +5) Block diagram of the correction device,
101° 4th (Jth) Showa () π (Calculator (34,:
Figure 2 shows the second embodiment of the present invention. Plotter diagram 5 of Shikuru Ryozo Yumi ME device
1; (The figure shows the conversion 1-Ii characteristic in Fig. 12, and Figs. Set 0), il Ming I J,
161”l C 1.2 r31 t:= t+ l
'J Le 111i jl-I: Block diagram of I analog-to-digital converter, No. 17-1 is shown in Fig. 21 is a diagram showing the principle of creating the conversion characteristics of the conversion 1.1 characteristic generator in the 1st I, and FIG. 20 is an internal configuration diagram of the conversion characteristic generator that is added to FIG. 16. Even image, 23: 1 with unevenness
iIiii image, 22.2/I: image fg! Human power equipment, 25:
Processing device, :g to :Image record, α device, 27:Fuhagimura n11 positive 1! 'j, Vfj, 28' tili
iIIE image, 29: Size 5L Fffl m, :(
0: Counter, 3] Near address register, 32 Near address af't"): Circuit, 3: 3 = Paralinear model chifur, to 34 37: Operation H::, 38: Multiplication fli;
, 39: Adder t-1 unit, 44, lI 5: Arithmetic circuit, 'l (], 50.51: Selection circuit, 52: Adder, 53.5/l: Register, 1; : 1: 11
1i positive, I, digital conversion) 1K, 64
: Human power device r'i, 57: Conversion characteristics, 58: Scanning, 59
:Sensor output signal, raO:Conversion characteristic generator, 65:Summerer, 06:Conversion circuit, 67:Timing light regeneration 1(), +; )(Near 1-1 noise generation circuit, 6(3
1, 662: Digital to analog conversion) (: +, 66
3. [36/1: λ amplification) 1:4.665: Analog-to-digital conversion)): (,7 [J: Pixel, 71 to 7
4: Lattice point, 601: AT1F6 conversion circuit, 602: Memory, 603. 6 (l lI nil (t circuit. Fig. 1 100 101 Fig. 2 Fig. 3 (a) (b) (c) Fig. 4 M1xXllyX Fig. 5 m block Fig. 6) 7 diagram Fig. 8 (a) (L, P, ) (L2P,) (b) Figure 9aa:<+b 1st (]χ1 α1+α2L αke αlI. ir 1 1 ::!1 1214th 13th:χ; Se/ Output City IE (yf No. 14 '>j to'!■5;4 85Z1t No. 161:4 Ushi No. 171' (1 az Figure 19 No. 201χ1 Continuation of page 1 0 Inventor Akira Tsuboi @ Inventor Hiroshi Ihara - [Phase] Inventor Yutaka Kubo 109 Hata, Ozenji, Asao-ku, Kawasaki-shi Hitachi, Ltd. System Hitachi University, Mika-cho 5-2-1 Hitachi, Ltd. Omika Factory

Claims (3)

【特許請求の範囲】[Claims] (1)a淡むらを含む画像を入力する手段、および入力
された画像の画像信号をアナログ・テイジタル変換した
後、あるいは変換すると同時に、入力画像を複数ブロッ
クを分割し、該ブロックの各格子点で陰面線近似のため
の係数をめ、該係数とブロック内相対座標により各点の
係数を算出して、a淡むらの補正を行う手段を有するこ
とを特徴とするディジタル画像歪補正装置。
(1) A means for inputting an image containing light unevenness, and after or at the same time converting the image signal of the input image from analog to digital, dividing the input image into a plurality of blocks, and dividing the input image into a plurality of blocks, and each grid point of the block. 1. A digital image distortion correction device, comprising means for determining a coefficient for approximating a hidden surface line, calculating a coefficient for each point based on the coefficient and relative coordinates within the block, and correcting light unevenness.
(2)前記ブロックの各格子点で陰面線近似のための係
数をめる場合、濃淡むらの少ない傅i像のデータを利用
することを特徴とする特a′(請求の範囲第11rrf
記載のティジタル画像歪tili正装置。
(2) Feature a' (Claim No. 11 rrf
A digital image distortion correction apparatus as described.
(3)前記ブロックの各格子点で陰面線近似のための係
数をめる場合、あらかじめ該係数を算出してアナログ・
ディジタル変換装竹内のメモリに記憶しておくことを特
徴とする特許請求の範囲第1項記載のディジタル画像歪
補正装置。
(3) When calculating coefficients for hidden surface line approximation at each grid point of the block, calculate the coefficients in advance and use analog
2. The digital image distortion correction device according to claim 1, wherein the digital image distortion correction device is stored in a memory of the digital conversion device Takeuchi.
JP59264268A 1984-12-14 1984-12-14 Corrector for digital picture distortion Granted JPS60150183A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP59264268A JPS60150183A (en) 1984-12-14 1984-12-14 Corrector for digital picture distortion

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP59264268A JPS60150183A (en) 1984-12-14 1984-12-14 Corrector for digital picture distortion

Publications (2)

Publication Number Publication Date
JPS60150183A true JPS60150183A (en) 1985-08-07
JPH0127467B2 JPH0127467B2 (en) 1989-05-29

Family

ID=17400806

Family Applications (1)

Application Number Title Priority Date Filing Date
JP59264268A Granted JPS60150183A (en) 1984-12-14 1984-12-14 Corrector for digital picture distortion

Country Status (1)

Country Link
JP (1) JPS60150183A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5490222A (en) * 1993-02-12 1996-02-06 Mitsubishi Denki Kabushiki Kaisha Video signal processing apparatus for processing correlative video signals obtained from a plurality of imaging devices

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5517428A (en) * 1978-07-24 1980-02-06 Citizen Watch Co Ltd Construction of alarm watch

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5517428A (en) * 1978-07-24 1980-02-06 Citizen Watch Co Ltd Construction of alarm watch

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5490222A (en) * 1993-02-12 1996-02-06 Mitsubishi Denki Kabushiki Kaisha Video signal processing apparatus for processing correlative video signals obtained from a plurality of imaging devices

Also Published As

Publication number Publication date
JPH0127467B2 (en) 1989-05-29

Similar Documents

Publication Publication Date Title
Debevec et al. Recovering high dynamic range radiance maps from photographs
Battiato et al. High dynamic range imaging for digital still camera: an overview
Zhuo et al. Defocus map estimation from a single image
Robertson et al. Estimation-theoretic approach to dynamic range enhancement using multiple exposures
KR100850729B1 (en) Method and apparatus for enhancing data resolution
CN107516319A (en) A kind of high accuracy simple interactive stingy drawing method, storage device and terminal
Mann et al. Quantigraphic imaging: Estimating the camera response and exposures from differently exposed images
US6804393B2 (en) Method of calculating noise from a digital image utilizing color cross correlation statistics
JPH0273477A (en) Method and device for intensifying digitalized image
CN110276831B (en) Method and device for constructing three-dimensional model, equipment and computer-readable storage medium
Easton Fundamentals of digital image processing
CN114140481A (en) Edge detection method and device based on infrared image
Manders et al. Robust hand tracking using a skin tone and depth joint probability model
JPS60150183A (en) Corrector for digital picture distortion
CN106851141B (en) A kind of asymmetric correction method of remote sensing images
JPS6119077B2 (en)
Heasley Numerical restoration of astronomical images
Johansen et al. Branch points in one-dimensional Gaussian scale space
CN113421292A (en) Three-dimensional modeling detail enhancement method and device
CN112017142A (en) Image color homogenizing method and device
Elhefnawy et al. Piecewise linear model for haze level estimation and an efficient image restoration technique
Lin et al. Photo-consistent synthesis of motion blur and depth-of-field effects with a real camera model
JPH03156306A (en) Direct measurement of volume of three- dimensional characteristic part in binocular stereo image
Sejzei et al. Evaluation of various digital image processing techniques for detecting critical crescent moon and introducing CMD–A tool for critical crescent moon detection
CN113192146B (en) Method for eliminating defocusing error in grid line projection phase shift technology based on deep learning