JPS60150183A - ディジタル画像歪補正方法 - Google Patents
ディジタル画像歪補正方法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
- digital
- conversion
- coefficient
- distortion correction
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/40—Image enhancement or restoration using histogram techniques
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
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:
;よび処理が可能とン多〕る。
、濃淡むらの少ないセンサで撮った画像データを利用し
、あるいはあらかじめ陰面線近似の係数を記憶すること
により、入力画像の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:
;よび処理が可能とン多〕る。
第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号 株式会社日立製作
所大みか工場内
図のセンサの原理図、第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号 株式会社日立製作
所大みか工場内
Claims (3)
- (1)a淡むらを含む画像を入力する手段、および入力
された画像の画像信号をアナログ・テイジタル変換した
後、あるいは変換すると同時に、入力画像を複数ブロッ
クを分割し、該ブロックの各格子点で陰面線近似のため
の係数をめ、該係数とブロック内相対座標により各点の
係数を算出して、a淡むらの補正を行う手段を有するこ
とを特徴とするディジタル画像歪補正装置。 - (2)前記ブロックの各格子点で陰面線近似のための係
数をめる場合、濃淡むらの少ない傅i像のデータを利用
することを特徴とする特a′(請求の範囲第11rrf
記載のティジタル画像歪tili正装置。 - (3)前記ブロックの各格子点で陰面線近似のための係
数をめる場合、あらかじめ該係数を算出してアナログ・
ディジタル変換装竹内のメモリに記憶しておくことを特
徴とする特許請求の範囲第1項記載のディジタル画像歪
補正装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59264268A JPS60150183A (ja) | 1984-12-14 | 1984-12-14 | ディジタル画像歪補正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59264268A JPS60150183A (ja) | 1984-12-14 | 1984-12-14 | ディジタル画像歪補正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS60150183A true JPS60150183A (ja) | 1985-08-07 |
| JPH0127467B2 JPH0127467B2 (ja) | 1989-05-29 |
Family
ID=17400806
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP59264268A Granted JPS60150183A (ja) | 1984-12-14 | 1984-12-14 | ディジタル画像歪補正方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS60150183A (ja) |
Cited By (1)
| 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)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5517428A (en) * | 1978-07-24 | 1980-02-06 | Citizen Watch Co Ltd | Construction of alarm watch |
-
1984
- 1984-12-14 JP JP59264268A patent/JPS60150183A/ja active Granted
Patent Citations (1)
| 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)
| 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 (ja) | 1989-05-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Battiato et al. | High dynamic range imaging for digital still camera: an overview | |
| Zhuo et al. | Defocus map estimation from a single image | |
| De Villiers et al. | Centi-pixel accurate real-time inverse distortion correction | |
| Flusser et al. | Recognition of blurred images by the method of moments | |
| JP3469582B2 (ja) | ピラミッド画像表示におけるウィーナ変形フィルタを用いる画像ノイズリダクションシステム | |
| Mann et al. | Quantigraphic imaging: Estimating the camera response and exposures from differently exposed images | |
| US4590608A (en) | Topographic feature extraction using sensor array system | |
| Russ | Image processing | |
| CN103106670B (zh) | 一种计算机输入设备的色彩模型自动建立及转换方法 | |
| US6804393B2 (en) | Method of calculating noise from a digital image utilizing color cross correlation statistics | |
| CN103218812B (zh) | 基于摄影测量的树木形态模型参数快速获取方法 | |
| KR20030020357A (ko) | 데이터 해상도를 향상시키는 방법 및 장치 | |
| JPH0273477A (ja) | 画素値修正方法及び装置 | |
| GB1512891A (en) | Method and apparatus for enhancing data | |
| CN110276831B (zh) | 三维模型的建构方法和装置、设备、计算机可读存储介质 | |
| US4653909A (en) | Sinusoidal determination of limiting optical resolution | |
| CN119941982A (zh) | 一种基于结构光图像增强的三维重建方法 | |
| CN117097997A (zh) | 一种反向图像信号处理的噪声图像合成方法 | |
| Manders et al. | Robust hand tracking using a skin tone and depth joint probability model | |
| CN113643446B (zh) | 一种下颌神经管自动标记方法、装置和电子设备 | |
| Easton | Fundamentals of digital image processing | |
| CN114140481A (zh) | 一种基于红外图像的边缘检测方法、装置 | |
| CN113034553B (zh) | 图像配准算法的评估方法、终端及存储介质 | |
| JP2023003764A (ja) | 画像処理装置、画像処理方法、及びプログラム | |
| CN117853356A (zh) | 人体热红外图像数据增强方法、装置、设备及存储介质 |