JPH03162685A - 合成画像を形成する方法及び装置 - Google Patents

合成画像を形成する方法及び装置

Info

Publication number
JPH03162685A
JPH03162685A JP1303486A JP30348689A JPH03162685A JP H03162685 A JPH03162685 A JP H03162685A JP 1303486 A JP1303486 A JP 1303486A JP 30348689 A JP30348689 A JP 30348689A JP H03162685 A JPH03162685 A JP H03162685A
Authority
JP
Japan
Prior art keywords
image
value
pixel
pixels
values
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
JP1303486A
Other languages
English (en)
Other versions
JP2581815B2 (ja
Inventor
Xavier A Flinois
ザビア エイ.フリノワ
Stefano E Concina
ステファノ イー.コンシーナ
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.)
Schlumberger Technologies Inc
Original Assignee
Schlumberger Technologies Inc
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 Schlumberger Technologies Inc filed Critical Schlumberger Technologies Inc
Publication of JPH03162685A publication Critical patent/JPH03162685A/ja
Application granted granted Critical
Publication of JP2581815B2 publication Critical patent/JP2581815B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/14Transformations for image registration, e.g. adjusting or mapping for alignment of images
    • G06T3/153Transformations for image registration, e.g. adjusting or mapping for alignment of images using elastic snapping
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/80Geometric correction

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)
  • Closed-Circuit Television Systems (AREA)
  • Testing Of Individual Semiconductor Devices (AREA)
  • Tests Of Electronic Circuits (AREA)
  • Testing Or Measuring Of Semiconductors Or The Like (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 皮血立1 本発明は、高解像度(分解能)画素画像の整合方法及び
装置に関するものであって、更に詳細には、集積回路装
置の動作状態のストロボ電圧コントラスト画像を整合す
る為の方法及び装置に関するちのである.この整合され
た画像は、例えば、集積回路の動的機能障害を解析する
上で有用である. 蓋3Ji奇 超大規模(VLSI)集積回路の内部特徴寸法が減少す
るに従い、より高速で且つより信頼性のある設計及びテ
ストを行なう要求が高まっている.従来の集積回路(I
 C)テスタは,テスト中の装置(DUT)の外部ビン
から情報を検索することが可能であるに過ぎず、従って
爾後の診断を行なうことが制限されていた.従来のテス
タで検知される欠陥は,部品内部の何れかの点における
ズレによって発生される場合がある.DUTが数千個の
ゲートを有するものである場合、欠陥を識別することは
極めて複雑且つ困難な問題となる.典型的なテスト動作
において、従来のICテスタ(ICの入力ビンへ励起を
与え且つIC出力ビンにおいてその結果を測定するもの
)は、IC装置をテストする場合にテストベクトルシー
ケンスにおいて何れかのベクトル(例えば,ベクトルV
)における欠陥を検知するちのである.このテストシー
ケンスは、多数のテストベクトルを有する場合があり、
例えば、100個又はそれ以上のテストベクトルを有す
ることがあり、各ベクトルはDUTのビンへ印加される
例えば入力電圧などのような一組の励起を表わしている
.検知された欠陥の発生源はチップ内のどこかにおいて
、そのシーケンスの最初のベクトル(ベクトル1)とそ
の欠陥が検知されたベクトル(ベクトルV)との間の何
れかのベクトル(例えば、ベクトルa)において発生す
る.従って,この欠陥は,前方へ伝搬し且つベクトルV
においてDUTの外部ビン乃至ボンドパッドに表われる
.従って、ベクトルaにおいて発生する欠陥の性質、位
置及び発生時間を同定即ち識別することが望ましい. 従って、チップの内部探査を行なうことが必要である.
長年の間、好適な解決方法は,低エネルギ(約1key
)の走査型電子顕微鏡(SEM)を使用してチップを探
査することであった.例えば、E.  Menzej2
及びE.Kubaj2ek共著の「集積回路の電子ビー
ムテストの基礎(FunadamentaI2s  o
f  EJ2ectron  Beam  Testi
ng  of  Intergarated  Cir
cuits)J.5スキャニング103−122 (1
983).及びE.  Pgies及びG.  Ott
o共著の「機械的及び電子プロープを使用する集積回路
内部の電圧測定(VoI2tage  Measure
ment  Inside  Integreated
  Circuit  Using  Mechani
cal2  and  Effectron  Pro
bes)J.IVスキャニングエレクトロンマイクロス
コビー1491−1500 (1985)の文献を参照
するとよい. 最近まで、SEMは高度の研究所の装置であって、熟練
した研究者によってのみ使用されていた.1987年に
おいて、rICS  5000Jワークステーションを
ベースとした電子ビームテストスステムがシュルンベル
ジェ社によって市販された.S.   Concina
.G.   Liu、L.  Lattanzi.S.
  Reyfman.N.   Richardson
、「ワークステーションをベースとした電子ビームテス
タにおけるソフトウェア統合(Software  I
ntegration  in  a  Workst
ation  Based  E−Beam  Tes
ter)」、インターナショナル・テスト・コンフェラ
ンス・プロシーリングズ(1 986).N.Ri c
hardson、rVLsI回路デバッグ用の電子ビー
ム探査(E−Beam  Prob ing  for
  VLSI  Circuit  Debug)J 
.BLSIシステムズデザイン(1987).S.  
Concina.N.  Richardson、rI
DS  5000:VLSI用統合診断システム(ID
S  5000:an  Integrated  D
iagnosis  System  for  VL
SI).7vイク口エレクトロニツク・エンジニアリン
グ(1987)などの文献を参照するとよい.更に、N
.   Richardsonに対して発行された米国
特許第4,706,019号及び第4,721,909
号を参照するとよい. 従って、測定及び情報の収集は、IC欠陥を回折する場
合の主要な問題ではない.しかしながら、この様なテス
タでICに関して得られた広範且つ詳細な情報を組織化
することは、迅速且つ効果的な欠陥の診断を行なうため
に重要である.絶対的な問題に対する解答を探すのでは
な<(「なぜこの装置は機能障害を起こすのか?」)、
相対的な問題として処理することが望ましい(「この装
置が既知の良好な装置と異なって動作するのはどこであ
るか、いつであるか且つなぜか?」).この後者の診断
手法によれば、IC内の動作欠陥をトレース即ち追跡す
ることが可能である. この様な診断を行なうために.0(興味のあるテストシ
ーケンスの開始)と最初の欠陥検知の間の興味のある時
間期間を選択する.この期間を、各間隔において、IC
が一定の動作を行なうように複数個の間隔へ分割する.
従って,各間隔は,ICの動作状態に対応する.興味の
ある状態が同定即ち識別された後に、各状態において、
既知の「良好」なIC装置とCUTOとの間での比較を
行なうことが可能である. この様な比較は、ストロボ電圧コントラスト画像で行な
うことが可能である.テスト中の装置(DUT)を、例
えばシュルンベルジエ「IDS5000Jなとのような
従来の電子ビームテストシステムにおける一連のテスト
「ベクトル」で励起させる.各テスト「ベクトル」は、
例えばICのビンへ印加される入力電圧などのような特
定した一組の励起を表わしている.該シーケンスの各ベ
クトルを印加する期間中、CUTは「ストロボウインド
」と呼ばれる時間期間に対する状態にある.電子ビーム
を、該シーケンスの開始に対してあるフェーズで繰返し
パルス動作させることにより、任意の所望のストロボウ
インド内のDUTの状態を表わすストロボ画像を得るこ
とが可能である.このストロボプロセスは、点火時期を
調節するために自動車のエンジンの動作を「凍結状態」
とさせるためにストロボライトを使用するものと類似し
ている.一連の状態画像が1個のスタックを形成する. 既知の良好なIC装置(ゴールデンダイとも呼ばれる)
を使用してこのプロセスを繰返し行ない、同一の一連の
テストベクトルに応答するその状態を表わすストロボ画
像からなるスタックを得る.例えば、各画像は、強度が
0及び255の間の値(8ビット値)によって表わされ
る可変強度の画素からなる5 1 2X5 1 2マト
リクスの形態でのデジタル形式でのグラフ表示とするこ
とが可能である.そのようにして得られた画像は、DU
Tの動作状態を表わしており、DUTにおける動作欠陥
を診断するために格納し且つ使用することが可能である
. 該画像を得た後に、そのスタックを比較することが可能
である.この比較は,画素毎に、DUTの画像からゴー
ルデンダイの画像を減算する形態を取ることが可能であ
り、その場合、比較された画像は同一の組の励起に応答
して発生されるちのである.この比較を画像毎に繰返し
行なって、「差分」画像のスタックを発生させる.一方
,これらの二つのスタックは,二つの異なった組の励起
に応答する単一のDUTの状態を表わすことが可能であ
る.例えば、DUTは、与えられた温度において又は与
えられた入力電圧に対し設計された通りに動作するかも
しれないが、より高い温度又は入力電圧に対しては機能
障害を起こす場合がある.一組の条件下でDUTの正し
い動作を表わすスタックの画像を得ることが可能である
、且つ第二組の条件下において機能障害を表わす第二ス
タックを得ることが可能である.これらのスタックを画
像毎に比較して、「差分」画像のスタックを発生させる
ことも可能である.良好な装置と欠陥性の装置とが同一
の対応で動作する場合には、それらの画像は同一である
.同様に、異なった組の励起が印加される装置は該異な
った組の励起に応答して同一の対応で動作する場合には
、それらの画像は同一である.これら二つの組の画像の
間の偏差は、テスト(即ち欠陥性の)装置におけるズレ
であると考えることが可能である.T.  May.G
.  Scott.E.Meieran.P.  Wi
ener.V.Ra○,rVLSIランダムロジック装
置の動的欠陥画像形成(Dynamic  Fauj2
t  Imaging  of  VLSI  Ran
domLogic  Devices)J.インターナ
ショナル・フィジックス・シンポジウム・プロシーリン
グズ(1984)の文献を参照するとよい.1つの状態
毎に1つの画像を基礎として、この比較処理は、それが
発生した位置から最初に検出された位置への欠陥の伝搬
を表わしている.このような比較処理は、動的欠陥画像
形成(DPI)として知られている. ちちろん、この比較が正確に作用する為には、画像の2
つのスタックが同一の条件を表わすものでなければなら
ない.即ち、チップの区域が同一であり、同一のSEM
差動条件下で同一の倍率で画像形成されたものであるこ
とが必要である.しかしながら、一つのスタックの画像
が、画像を収集した場合にSEMに関してチップの配向
状態が多少異なっていた為に、別のスタックの6のに関
して並進又は回転している場合がある.完全な整合は必
ずしも可能なものではない. 従って、1つのスタックの各画像を他方のスタックの対
応する画像と整合させる為に、一方のスタックの各画像
に関して空間的「ワービング(撓曲)」動作を必要とす
る場合がある.このワービング動作は、2つの画像の比
較(差異)が画素毎になされるので、正確な整合状態と
なる.不整合がある場合には,合成(差分)画像内に表
われる.このような精度を与えるワービング手順は、そ
のワービング動作を実旅する為の専用ハードウエアがな
い場合には、迅速且つ対話的なn F I 隷断におい
て使用する場合に許容できないような長時間を必要とす
る(数十秒). 画像が収集され、フィルターされ、整合され且つ画素毎
に比較されて合成「差分」画像のスタックを形成すると
、一連の画像を介して欠陥伝搬をトラッキングすなわち
追跡することを可能とする差分画像を表示することが可
能である.この「差分」画像を処理して、解析の為に重
要な情報(欠陥伝搬)のみを強調させることが可能であ
る.DPIプロセスを効率的に実行することは、比較す
べき画像の迅速且つ正確な整合を必要とし、従ってこれ
らの画像(又はそれらの差分)が表示されると、DUT
内の欠陥の伝搬を、その装置の機能性について殆ど知識
を有することのないオペレータであってち、容易にトレ
ースすることが可能である. DPIセッションにおける幾つかの作業は繰返し的であ
る.(例えば、画像収集及び処理).例えば、制御した
悲様で得られた画像に対して正しいSEM動作条件が見
付かると、スタック内の画像の残りのちのに対してバッ
チ処理で同一の動作を繰返し行うことが可能である.エ
ンジニアが診断作業に集中することを可能とする為に、
このような繰返し作業を自動化することが望ましい.同
時に、エンジニアが作業に対するパラメーターを設定し
且つ顕著な遅延なく結果を観察することが可能であるよ
うに、対話的作業が可能であることが望ましい. 画像整合用の従来方法は公知である.例えば、従来使用
されているシュルンベルジェN DS5000Jシステ
ムにおいて使用されている方法は、第二画像の3つの選
択した位置とそれぞれ整合されるべき第一画像からの3
つの位置の選択等、浮動小数点値からなるマトリクスの
形態で幾何学的に平坦な変換を計算し、且つそのマトリ
クス値を使用して第一画像の位置を浮動小数点演算を使
用して第二画像の画素と整合した画素を有する合成画像
の位置へ変換させる.この整合の品質は許容可能なちの
であるが、その処理は望ましい程度に迅速なちのではな
い.「対話的」動作を達成する為には、例えば画像整合
等の任意の選択した処理を実施し且つその結果を、例え
ば、2秒を超えることのない期間以内で表示することが
望ましい. 幾何学的演算、フィルター動作及びパターン認識等を実
行する為にコンポリューション即ち畳込みを使用するそ
の他の画像形成システムは従来公知である.このような
システムは、時折、大きなコンポリューション即ち畳込
みの計算を高速化する為に高速フーリエ変換を使用して
いる.データに関してどれ程の演算を実施するかという
ことは予め解っていないので、浮動小数点演算がしばし
ば使用される. 画像処理をスピードアップする別の方法は、高速で特定
の作業を実施すべく設計された専用ハードウエアプロセ
サーを使用することである.しかしながら、このような
プロセサーのコストは、テストシステムのコストを著し
く増加させることがある.従って、その代わりに、多大
な処理時間を失うことなしにこのことを実施することが
可能である場合には、画像処理を実行する為に標準的な
エンジニアリングワークステーションの汎用プロセサー
を使用することが望ましい. 狂一迫 本発明は、以上の点に鑑みなされたものであって、上述
した如き従来技術の欠点を解消し、対話的な「感じ」を
有しているが処理を実施する為の専用ハードウエアを必
要とすることのない動的欠陥画像形成(DPI)を許容
するのに充分な速度を有する例えばシュルンベルジエ「
■DS 5000」等のようなシステム内に設けられる
クイブのワークステーションにおいてストロボ電圧コン
トラスト画像の整合を行う方法を提供することを目的と
する. 本発明の別の目的とするところは、1つの画像の特徴を
別の画像の特徴と画素毎に整合させる方法及び装置を提
供することである. 本発明の更に別の目的とするところは、集積回路のDP
I解析における手助けとして集積回銘動作状態を表わす
ストロボ電圧コントラスト画像を整合する方法及び装置
を提供することである.本発明の更に別の目的とすると
ころは、整合を実施する為に専用の計算を行うハードウ
エアを必要とすることなしに、電子ビームテストシステ
ムと関連する汎用ワークステーションにおける電圧コン
トラスト画像の迅速且つ対話的な整合を提供することで
ある. 透一羞 本発明の一形態によれば、合成画像が第一画像から形成
され、その合成画像が第二画像と画素毎に整合される.
これらの向像の各々は、1組の装置によって表示可能な
画素を有しており且つ表示された場合に画素によって画
定される視覚的特徴を有している.各画素は,表示座標
位置及び強度値によって画定される.その整合は一連の
ステップによって実施することが可能であり、例えば、
(a)第一画像の座標位置を合成画像のそれぞれの座標
位置ヘマップすることが可能な浮動小数点値からなるマ
トリクスを計算し、(b)フローティングポイント値の
マトリクスをそれぞれの固定小数点値表示へ変換し、(
c)該固定小数点表示を使用して、合成画像の各表示座
標位置に対して第一画像の対応する座標位置を固定小数
点演算で計算し、且つ(d)第一画像の対応する座標位
置を取巻く画素の強度値から合成画像の各表示座標位置
に対する画素強度値を外挿し、その際に合成画像を画定
する1組の画素を発生する.合成画像は、第二画像の特
徴と整合された特徴を有する第−画像の変形を有してお
り、合成画像の画素は第二画像の画素のそれぞれの表示
座標位置と整合されている.第一及び第二画像は、電子
ビームテストシステムによって発生され且つ集積回路動
作状態を表わすストロボ電圧コントラスト画像とするこ
とが可能である. 該マトリクスの計算は、表示スクリーン上に第一及び第
二画像を表示し、第一画像の3つの選択した特徴の座標
位置を記録し、合成画像において第一画像の選択した特
徴と整合されるべき第二画像の3つの選択した特徴の座
標位置を記録し、且つ記録した座標位置から該マトリク
スを計算する各ステップを包含することが可能である.
選択した位置を表わす値は、好適には、副画素精度で記
録される.このマトリクスは、好適には、第一画像と第
二画像との間の回転、傾斜及び拡大歪を補正する為の係
数を表わす値からなる2X2ブロックを有すると共に、
第一画像と第二画像との間の点の並進を表わす値からな
る1×2ブロックを有している.このマトリクス値を表
わす固定小数点は、レジスター内に格納することが可能
である.第一画像の対応する座標位置は、好適にはマイ
クロプロセサー装置において,計算する.画素強度値を
外挿する為には、周囲の画素に対する対応する座標位置
の距離を表わす浮動小数点値を計算し、その距離値を固
定小数点値に変換し、且つ固定小数点距離値及び固定小
数点演算を使用するプロセサーにおける周囲の画素強度
値から合成画像画素値を計算することが可能である.こ
の外挿は、合成画像画素値が周囲の画素強度値の平均を
有しており、その平均が距離値によって重み付けされて
いるようなものである.好適実施例においては,この外
挿は、各々の値が合成画像画素強度値に対し1つの周囲
の画素の重み付けした寄与度を表わしている固定小数点
値からなるテーブルを計算し且つそのテーブルの値を格
納する各ステップを有している.各周囲の画素に対して
該テーブルから値を検索し且つ結合して合成画像画素強
度値を形成する.該テーブルの各格納した値は格納した
インデックス値と関連している.該テーブルから値を検
索する為に、少なくとも1個の周囲の画素の値から及び
対応する座標位置に対するこのような周囲の画素の距離
からインデックス値を計算する.その計算したインデッ
クス値を格納したインデックス値とマッチングさせて、
該テーブルから検索されるべき値を選択する.合成画像
の調製即ち形成は、好適には、電子ビームテストブロー
プシステムと関連する適宜にプログラムされた汎用プロ
セサーを使用して実施する. 笈立男 以下、添付の図面を参考に、本発明の具体的実施の対応
について詳細に説明する. 集積回路を回折するのに有用な電子ビームテストプロー
ブシステム10を第1図に概略ブロック図で示してある
.この電子ビームテストブローブシステムは3つの機能
的要素を有している.即ち、電子ビームテストブローブ
12と、回路励起器l4と、ディスプレイターミナルl
8を具備するデータ処理システム16とである.解析す
べき集積回路は、電子ビームテストブローブl2内に配
置し、従ってICの種々の点において電位測定を行うこ
とが可能である.前記洞定を行う点は、バス22を介し
て、データ処理システム16によって電子ビームテスト
ブローブl2へ送られる.回路解析の期間中、テスト信
号パターンが、バス24によってテスト中のICへ接続
されている回路励起器l4によって集積回路へ印加され
る.データ処理システム16は、更に、テスト信号パタ
ーンに対しての電位測定のタイミング及び使用されるテ
スト信号バクーンを特定する為に使用することが可能で
ある.電子ビームテストプローブシステムは、ディスプ
レイターミナル18を介して命令を入力するオペレータ
によって制御される.データ処理システム16をブロッ
ク図で第2図に示してある.データ処理システムl6は
,二つの大きな機能的グループに分割することが可能で
あり、即ち計算即ちコンピュータシステム50と電子ビ
ームテストプローブ制御システム52とである.計算シ
ステム50の新規な側面に関してのみここでは説明する
.尚この基本的なシステムは米国特許第4,706,0
19号に記載されている.計算(コンピュータ)システ
ム50は、テスト中のICの高分解能電圧コントラスト
画像を格納し、処理し且つ表示することが可能である.
コンビュークシステム50は、マイクロプロセサ−54
、RAMメモリー56、ディスクドライブを具備するこ
とが可能なデータ格納モジュール58を有している.更
に、コンピュータシステム50は、例えば少なくと(1
1 000X 1 000画素の解像度を持った高解像
度カラーディスプレイターミナル6lを駆動する為のデ
ィスプレイインターフェース60を有するディスプレイ
ターミナル18を有している.ディスプレイターミナル
l8は、更に、命令を入力する為のキーボード及びディ
スプレイスクリーン上にある点を指示し且つ指定する為
のポインター手段(例えば「マウス」)を有している.
ディスプレイターミナルl8は、更に、グラフィックデ
ィスプレイを発生させるのに必要な時間を改善する為の
それ自身のグラフィック発生器62を有することが可能
である.コンピュータシステム50は、更に、通信バス
66を介してその他のコンピュータシステムと通信する
為のインターフエイス64を有することが可能である. 第3図は、シュルンベルジエr■DS  5000」電
子ビームテストブローブシステムで得られる画像からな
る2つのスタックから選択される一対のストロボ電圧コ
ントラスト画像を示している.各画像は、5 1 2X
5 1 2画素のアレイから形成されており、各画素は
0と255との間のデジタルグレースケール(中間調)
強度値を持っている.各スタックは、最大で100又は
それ以上のこのような画像から構成されている.DPI
解析の為に表示することが可能な差分画像のスタックを
用意即ち作成する為に、このような2つのスタックを画
像毎、画素毎、に比較することが望ましい.第3図の左
側に示した画像は、第3図の右側に示した画像の視覚的
特徴(1、2、3として示してある)と類似してみえる
が整合されていない視覚的特徴(例えば、l、2、3と
して示してある)を有している.この不整合は,多数の
要因のいずれかによる6のであり、例えばテストプロー
プ内の集積回路の回転、傾斜又はその他の異なる物理的
な配白状態、並進、表示倍率等である.第3図の右側に
おける画像内には円形領域20を示してあり、画像の拡
大した部分を示している.第4図は、3つの電圧コント
ラスト画像を示しており、左上の画像(バックグラウン
ド内)は第3図の右上の画像と同じであり、右上の画像
(バ8クグラウンド内)は第3図の左上の画像と同じで
あり,且つ真ん中下側の画像(フォアグラウンド内)は
バックグラウンド内の2つの画像の「差分」画像である
.不整合の為に、差分画像は何等診断をする価値がない
6のである. 第5図は、更に、3つの電圧コントラスト画像を示して
おり、左上の画像(バックグラウンド内)は第4図の左
上の画像と同じであり、右上の画像(フォアグラウンド
内)は第5図の左上の画像と整合させる為に本説明に従
って処理された第4図の右上の画像の変形であり、且つ
真ん中下側の画像(中間グラウンド内)は第5図の左上
及び右上の画像の「差分」である.第5図の中央下側の
画像は、画像が収集されるIC装置内の位置における異
なった電圧レベルを可視的な6のとしている.(例えば
、異なった論理レベルでの表示。)本発明の整合処理は
、大略、一方の画像を「ワーブ(撓曲)コさせるべく作
用し、画像の相対的な幾何学的歪を補正する.従って、
「ワーブ」した画像は、それが整合された画像と集積回
路の同一の区域を表わしている. 画像の各スタックが収集される集積回路が画像収集の期
間中物理的に移動されず,且つ該スタックの収集した画
像が同一の態様で処理される場合には、2つのスタック
の態様する画像の間の不整合は画像毎に同一であるとい
うことに注意すべきである.従って、1つのスタックの
画像を他方のスタックからのそれに対応するちのと正確
に整合する整合処理は、いくつかのスタックを介して画
像毎に繰返し行い整合した画像のスタックを発生するこ
とが可能である. 小数声、表示を使用する従来の整合方法は、ワープされ
るべき画像の512X512画素を、同一のタイプのワ
ークステーション上の専用浮動小数点演算コプロセサー
によっても約80秒で処理することが判明した。本発明
は、固定小数点表示を使用するワープ動作を実施する方
法を提供するものであり、UNIXをベースとしたサン
マイクロシステムズ3/160ワークステーションでC
言語でプログラムした場合に約1.3秒かかることが判
明した.本発明技術による速度は、専用コプロセサーを
必要とすることなしに、動的欠陥画像形成註断において
使用する為に迅速且つ対話的に画像を整合させることを
可能とする。
このような整合の原理について、第6図及び第7図を参
照して以下に説明する.2つのスタックの対応する画像
の間の相対的幾何学的歪の補正は、最初に、第二画像の
座標系(x.y)を第一画像の座標系(U.V)ヘマッ
ピングすることによって行うことが可能である。
第6図は、ボックス100内において、第一画像を画定
する特徴130、132、134、l36、138から
なるグループを示している.ボックス100の特徴は、
グラフィックフレーム内の画素によって形成されており
、それらの(U、■)座標系内に位置されたドットの格
子によって表わされている。ボックス200は、第二画
像の(x.y)座標系に関連して第二画像のそれぞれの
対応する特徴230、232,234、236,238
の位置を示している.ボックス200の特徴230、2
38と相対的なボックスlOOの特徴130−138の
寸法及び間隔の比較から、第二画像の倍率は第一画像の
2倍(2×)であり、即ち第一画像の特徴は第二画像の
特徴の間隔の半分の間隔であることが理解される.ボッ
クス100と200とを比較すると、ボックス100の
特徴はボックス200の対応する特徴の多少下側及び多
少右側へ位置していることが分る。更に、特徴130、
132、134、136、l38のそれらを画定する画
素に対する空間的関係は特徴230,232,234,
236.238のそれらを画定する画素に対する空間的
関係と同一ではないので、第一画像(ボックス300)
の結果的に得られる「ワーブ」した変形を発生する為に
、第二画像(ボックス200)の画素位置において特徴
130,132、134、136、l38に対しての画
素強度値を外挿することが必要である, 第6図のボックス300は、このような結果的に得られ
る画像の特徴330を示しており、それは、約1−1/
2増分上側へ且つ約172増分左へ拡大且つ並進されて
おり、特徴330は第二画像の特徴230の位置と整合
されている.特徴330を画定する画素に対して割当て
られた強度値は、以下の如く、第一画像からの画素値か
ら決定されねばならない. このようなマッピングを実施する為の単調伝達関数を一
般的に以下の如く表現することが可龍である. x=X  (u,  v) y = Y  ( u +  v ) 尚、 X.YcC (即ち、 X及びYは連続的に 派生可能な関数である.) 一般的に説明すると、X及びYのグラフは、使用可能で
はなく、且つ例えば多項式近似等のような経験的モデル
を使用せねばならない.このようなマッピングを実施す
る1つの可能な方法は,本発明の好適実施例を参照して
以下に更に説明する. 整合されるべき画像は離散的空間内にあるので(画像は
画素のアレイである)、x又はyが整数値でない場合に
は、合成即ち結果的に得られるグラフィックフレーム内
の与えられた画素p(x、y)の補正された表示強度値
が、コンポリューション即ち畳込み操作を使用して、ソ
ースグラフィックフレーム内のその画素の対応する位置
を取巻く画素から外挿されねばならない.本発明におい
ては、4個の周囲の画素を使用し、且つコンポリューシ
ョン(畳込み)は双線形外挿となる.このような外挿技
術は、従来公知であり、例えば、W.   Pratt
、「レジクル画像処理(.DIGITAL  IMAG
E  PROSESSING)J.1978,115頁
の文献を参照すると良い.第7図は、第一画像lOOの
(u,v)座標系への第二画像200の(x.y)座標
系内の画素の並進位置である位置440とこの並進した
位置を取巻く画素pi.p2、p3.p4との間の関係
を詳細に示しており、第7図の3番目の字源pは画像強
度を表わしている.第7図中疑問符(?)は、外挿され
るべき画素440の値を表わしている. 0不等記号で右側が大きいと等記号と組合わせた記号δ
とX不等記号で右側が大きいと等記号で組合わせた記号
1及び0不等記号で右側が大きいと等記号の組合わせた
記号δy不等記号で右側が大きいと等記号の組合わせた
記号lで位置440(X+δx.y+δy)における外
挿値は次式で与えられる. p(x+δ x,y十δ y)= [p(x.ylδy+pfx,y+ll ( 1−δy
l ] [1−δX]+[pfx+1,yl  δy+
pfx+1.y+11  (1−δy))   δ X
(11 尚、X及びyの各々は512の値を有している.結果的
に得られる「ワープ」した画像を構築する為にはかなり
多数の演算を行なわねばならないことが明らかである. 標準的なコンポリューション(畳込み)アルゴリズムは
、計算の基本的目的として、浮動小数点演算又は十進数
変数を使用する.しかしながら、浮動少数点変数によっ
て与えられる多数の小数は、DPI解析等の適用におい
て使用される画素画像の処理に対しては冗長であること
が判明した.浮動小数点表示の精度は、小数点の後の桁
数である.0と255との間の強度値を持った画素の5
12X512マトリクスの場合には、3個の桁の精度が
適当であることが判明した.従って、浮動小数点変数を
有する比較的遅い演算は、それらを整数又は「短い」変
数で置換することによって回避した.固定小数点演算を
使用することにより、処理時間が著しく節約された. 固定小数点演算は、整数のビットを2つの部分に分割す
る.下側のビットは、変数の少数部分に取っておかれる
.上側のビットは、変数の仮数に取っておく。32ビッ
トの整数の場合、各部分が16ビットを有していると,
エミュレートした浮動小数7屯変数Zは、5の精度と等
価である.なぜならば、最小変数は以下の如くであるか
らである. 符号付き又は符号無しの変数を使用することが可能であ
り、且つ全ての基本演算を実施することが可能である.
ある演算は極めて高速となり、例えば、2の草による乗
算又は割算は瞬間的である.全く、それらの演算はマイ
クロプロセサーにとって実際には単一シフト演算にすぎ
ず、2ないしは3クロックサイクルを必要とするに過ぎ
ない. それに6拘わらず、浮動小数点表示での乗算は、マイク
ロプロセサーの「内部的」機能ではないので,注意が必
要である. 浮動小数点演算の代わりに固定小数点演算を使用する場
合を説明する為に、ある整数を2つの等しい部分に分割
する例を取って説明する.「1」の表示は、16個の二
進位置、即ち216又は65、536位置だけシフトし
たlである.乗算が内部演算である場合、IX1=2”
で表示され、21 @ X 2 + 8 == 2°で
はない.この結果の2つの部分は別々に計算されねばな
らない.この切断及び接合動作のコストは、仮数に対し
て2つのシフトとを、結果を再構成する為の加算命令と
を加えたものである.しかし、畳込みにおいて発生する
殆どの乗算は整数による外部的乗算であるから(画素座
標及び強度値は整数として表わされている)、このよう
な演算は直接的に行なうことが可能であり、且つ基本的
且つ高速な演算である.本発明の好適実施例においては
、16ビットに短くした値(「ショートJ)を8桁固定
小数点演算と共に使用する.画像の画素は8ビットデー
クである. 上述した式(1)を1例としてとると、実施すべき基本
的演算は次式の如く書き表わされる.p=[  ・ α
 (1−fl)  ・ β(2) 尚、α、β、tN (a及びβは整数値)及びO不等記
号で右側が大きいと等記号との組合わせl不等記号で右
側が大きいと等記号との組合わせの記号1である. 固定小数点演算表示において実施せねばならない、1つ
の減算,整数値による2つの外部乗算及び1つの加算が
ある.10進数表示を使用することを回避することが望
ましいので、本発明は次のごとき置換を行う. 2をkで置換し、k=256・β 最終的な結果は、8だけシフトし(256による割算と
等価)、結果的に得られる画素の値を与える. 例えば、以下の如き条件が成立するちのとする. β=0.  34 α=209、 β=1 98、 k=87l−βは25
6−87=169で表わされる. p′はpの固定小数点表示であり,式(2)は以下の如
くなる. p’=k ・α (256−k)  ・β      
(3)p=p′>s                
    f41尚、「》」はマイクロプロセサーの右シ
フトである.結果を比較すると、 (2)は201.74を与え且つpεN−jp=201
, (3)はp’ =51645を与え、 (4)はp=51645→>8=20 1を与える, 本発明は、好適には、以下の如く、一連のステップで実
施される. (l)2つのスタックの各々から1つずつ、2つの画像
を選択し且つディスプレイターミナルl8上に表示させ
る (2)キーボード又はプリンター手段(それは,マウス
又はその他の適宜の装置とすることが可能である)を使
用して、カーソルを第一画像の任意の選択した特徴の位
置へ移動させ、マーカーを表示された画像上のその位置
に配置させ且つその特徴の座標位置を好適にはRAM5
6内に記録させる.該カーソルを、第二画像の任意の選
択した特徴(好適には、オベレークが第一画像の選択し
た特徴に対応する6のと考える特徴)の位置へ移動させ
、マーカーを表示した画像のその位置上に配置させ且つ
この特徴の座標位置を記録する。
ついで、選択及び記録プロセスを各画像の2つの付加的
特徴に対して繰返し行う.ワーブ動作(伸長十並進)の
振幅をユーザーによって対話的に設定する. 本発明の好適実施形熊においては、コンピュータシステ
ム50は、画像の個々の画素を識別することが可能であ
るように表示されるべき画像の−部が選択された倍率の
ウインド内に入ることを可能とする倍率特徴と、選択し
た特徴の位置を副画素精度で記録する能力を有している
.例えば、コンピュータシステムは、好適には、画像の
一部の4×倍率及びカーソルの配置及び1/4ないし1
/2画素の精度でのカーソル位置の記録を可能としてい
る.画像の拡大部分の一例を第3図の右側の図の円で囲
った領域内に示してある.このような拡大は、画像の適
切の整合を助け、その場合、例えばIC装置内のメタル
導体を表わすラインの幅等のような特徴を最小で2個の
画素によって画定することが可能である. 拡大は本発明の整合方法によって実施することが可能で
あり、拡大すべき画像部分は「ワーブ」されるべき画像
として取扱われる.拡大が4×即ち4倍であると、その
拡大された変形は、もとの画像の選択した部分と同一の
特徴を表わす画素の16倍を有している.拡大された変
形内の各画素の値は、拡大されたちのの画素に対応する
ちとの画像の位置を取巻く画素の値から外挿される。
(3)ついで、浮動小数点値の3×3マトリクスの形態
における幾何学的平坦変換を計算し、それは第一画像の
選択した特徴の位置を第二画像の特徴の中にマップする
. 例えば、第一画像上の3つのマーカーが座標b[i] 
 ・X及びb[i]  −yを有しており且つ第二画像
上の3つのマーカーが座標a[11 ・X及びa[i]
  ・yを有しており、lが1ないし3の間の値を有す
るちのと仮定する.2つの3×3マトリクス「B」及び
「A」が構成され、以下の如く表わされる. B[iH(lI・b(il−x   と   A[i]
[1]J=a[iJ−xB[i]  [tl=b[l]
  ・y   と   A[i]  [1]・a[i]
 ・yB[il  [2] =1.0      と 
  A[i]  [2]=1.0A=M−Bに対する固
有マトリクスMは、所望の幾何学的変換マトリクスであ
る。2×2副マトノクスM [i]  [j]内の値(
i.jは各々○、1の値を有することが可能である)は
、画像の回転、傾斜及び拡大歪を補正する為の係数であ
る.M[O]  [2]の値は、水平並進係数であり、
且つM [tl  [2]の値は垂直並進係数である。
(4)ついで、変換マトリクスM内の値を固定小数点表
示(FPR)へ変換する.本発明に基づいて、FPRは
、好適には、浮動小数点変数の255倍の大きさの整数
である.例えば,0−はOとなり、0.25は64とな
り、0.5は127となる.FPR変換の利点は、実施
されるべき爾後の演算は、整数演算となり、それは例え
ば、マイクロプロセサ−50等の汎用プロセサーにおけ
る浮動小数点演算を実施するより6何倍も高速である. (5)ついで、その結果得られる画像を、ライン毎に走
査して、第一画像の位置を、変換マトリクスのFPR値
を使用して、合成画像の画素位置へ変換する.この変換
を行う場合、且つ合成画像の画素値を計算する為に必要
とされる時間を最小とする為に、テストを実施して、合
成画像の座標内にマップされた場合に、第一画像の与え
られた位置が合成画像の所望の表示フィールド内にある
か否かを決定する.意図するところのものは、512ラ
イン及び512列の画素からなる表示フィルドを有する
第二画像と画素毎にマッチする合成画像を有ることであ
る.変換した画素位置が表示フィールドの外側である場
合には、例えば表示された場合に黒画素を発生する為の
固定強度値が割当てられる. 更に詳細には、合成画像(画像l整合済)が各画素が第
一画像(画像1)の画素値から計算される値を持ってい
る画素から形成されている.画像lにおける位置Ll(
例えば、第7図における位置440)は、画像位置整合
済においては位置Lla(第6図における位置330)
へ変換される.画像l整合済における位置Llaにおけ
る画素は、画像1内の位置Llに対して外挿した値を維
持する.画像1整合済における512X512位置Ll
aの各々に対して、テストを行って、画像l内の対応す
る位置が、画像1整合済の座標系へ変換された場合に、
画像l整合済の表示フィールド内にあるか否かを決定す
る.換言すると、画素位置Llaの水平及び垂直座標の
いずれらが512よりも大きいが又は0よりも小さくな
ることはない.第8図は、この原理を概略示している。
画像lの位置aは、画像l整合済の表示フィールド外側
の位置Aに対応しており、一方画像lの位置bは画像l
整合済の表示フィールド内の位置に対応している. 通常、位置Llは10進数値である.しがし、FPR値
が使用されるので、4つのテスト(表示フィールドの水
平境界に対して2つ及び垂直境界に対して2つ)を、整
数論理演算で実施することが可能である.従って、この
テストは、浮動小数点演算でなされるよりも何倍ら高速
である。このテストが表示されたフィールドのフレーム
の外側の画素を診断する場合、その結果は黒へ設定され
る。
(6)位置Llは小数とすることが可能であるので(即
ち、L1の座標は第一画像の画素間に位置する場合があ
る.画素値を周囲の画素がら外挿即ち補外することが必
要である.この外挿即ち補外は、位置L1を取巻く4個
の画素の重み付けをした平均値として考えることが可能
である.2つの画素の間の距離が1単位であると、各画
素に対して割当てられる重みはl−その位置への割合的
単位距離である. 好適には、この外挿ないしは補外を実施するのに必要な
時間を最小とする為に、予め計算した「ルックアップ」
テーブルを使用する.このようなテーブルを使用するこ
とは、FPRを使用することによって可能とされ、整合
処理のこの段階において、予め計算した値のアレイの要
素番号となることのある整数を取扱うことが必要である
に過ぎない.各結果を個別的に計算する代わりに、イン
デックスが付けられた静的なテーブルの値から単に選択
することが可能である.このようなテーブルは、好適に
は、コンピュータシステム5oの始動と共に自動的に計
算され、且つRAM56の活性メモリー内又はデータ格
納装置58内の「スワップ空間」内に格納することが可
能である.以下に与える例においては、テーブルは2X
65536個の値を有しており、約1 28Kバイトの
格納空間を必要とする. ルックアップテーブルを使用して画素値を外挿する好適
な対応について説明する.興味のある位置はLLは上側
が画素p1及びp2及び下側がp3及びp4の画素によ
って取囲まれている.これら4個の画素の値は、4ビッ
トに圧縮される.δX及びδyの固定小数点表示を使用
する.(δX及びδyは、上述した如く、0と1の間の
小数値である). ついで、3つのステップで平均化を行う。即ち、pt及
びp2の平均を計算し、p3及びp4の平均を計算し、
且つこれら2つの平均を加算する. 最初の平均化はルックアップテーブルを使用し、それは
16ビットのエントリーを有している(各々が画素強度
値pi及びp2に対する4ビット及び各々が位置値δX
及びδyに対する4ビットを表わしている).上の2つ
の画素の寄与はルックアップテーブルにおいて読取られ
る.ルックアップテーブルから値を選択する為に使用さ
れるインデックスは、δX及びδyに対して0と1との
間の小数値を計算し、(興味のある位置はL1及びそれ
を取囲む画素の位置座標から)、δX及びδyに対しF
PRを計算し、画素強度値p1及びp2を抽出し且つこ
れら4つの値を単一のショート(それは、シュルンベル
ジエrIDS  5000」テストシステムのUNIX
をペースとした計算システムに対しては16ビットであ
る.)内に入れることによって、i44fNされる.2
番目の平均化は同一のルックアップテーブルを使用する
.p3、p4、δX、δyから形成されるインデックス
は、テーブル内において、下側の画素の寄与である値へ
ポイントする.その結果得られる画素の値(0と255
の間の値)は、2つの寄与度の和である. シュルンベルジエrlDs  5000Jシステムの場
合、コンピュータシステム(計算システム)50は、適
切にインデックスしたルックアッブテーブルが予め計算
されている場合には、0ないし255の画.素強度値を
有する2つの512×512画素画像を2秒未満で整合
させることが可能である. このようなルックアップテーブルは以下の如くに準備す
ることが可能である.δXのδyは、位置L1の重心位
置を画定し、且つPIおよびP2は上側2つ又は下側2
つの画素のいずれかである。PiU及びPiLが上側及
び下側の画素であるとすると、変換は、P=(1−δy
)*(〔1−δx)*PIU+δX*P2U)+δy*
((1−δx)*PIL+δx*P2L)である。
例えば、δx=0.25、δ,y=0.8.px=40
、p2=200、p3=60及びp4=220であると
、δXを表わす4ビットは8(10進数において)であ
り、δyを表わす4ビットは12であり、p1を表わす
4ビットは2であり、p2を表わす4ビットはl2であ
り、p3を表わす4ビットは3であり、且つp4を表わ
す4ビットはl3である.従って、最初のルックアップ
に対するインデックス値は、3  8  12  12
であり、且つルックアップテーブルからの対応する値は
42であり,結果として得られる画素に対するpi及び
p2の寄与を表わしている。2番目のルックアップに対
するインデックス値は、3813 12であり、且つこ
のルックアップテーブルからの対応する値は54であり
、p3及びp4の寄与を表わしている.これらの値の和
は、96であり、それは結果として得られる画素の値で
ある6本発明の好適実施形態においては、32個のグレ
ーレベル(画素強度値)が使用されているに過ぎないが
、それは許容可能な画像表示品質に対して適切なもので
あることが判明した.これら32個のレベルは、33な
いし64の範囲の値によって表わされる.尚、本発明の
技術的範囲を逸脱することなしに、この好適実施例とは
異なったその他の変形例を使用することも可能である.
結果として得られる画像(第二画像と整合されるように
「ワーブコされた第一画像の変形画像)を、画素毎に強
度値を減算することによって第二画像と比較することが
可能である.各画素の強度は、8ビットによって表わさ
れており、従って0と255と間の値を有する.画像を
減算すると、差分値は−255と+255との間となる
.差分画像は尚且つ画素値でなければならないので、2
55か6255の範囲は通常の0か6255のグレー範
囲(中間調)レベルにマップされねばならない.この差
分画像は、興味のある3つの強度範囲ないしはバンド内
に入る情報から構成されている.最初のバンド(0から
レベルLl)は、第一画像内にはないが第二画像内の特
徴に対応している.2番目のバンド(Llからレベノレ
L2)は、両方の画像における同一の特徴に対応してい
る.3番目のバンド(L2から255)は、第二画像に
はないが第一画像内の特徴に対応している. これらのバンドを検知するために、第9図に示した如く
、3レベルのスレッシュホールドフィルターが使用され
る。従って、初期的に8ビットの深さのグレーレベル(
256個の異なったグレーレベル)を有するヒストグラ
ムは,黒(0からLl)、グレー即ち灰色(LlからL
2)及び白(L2から255)の3つのレベルのスケー
ルヘマップされる. 当業者にとって明らかな如く,集積回路装置のトボグラ
フィ即ち地形的特徴の幾つかの特徴を強調することが可
能である.例えばグレーレベルでの画像を使用する代わ
りに.W色モードを使用することが可能であり、明るい
色(即ち、青、黄色及び赤)を表示させることが可能で
ある.情報を最適化するために、スタックの選択した画
像(「スクー画像」)に関してレベルを対話的に設定し
、且つスタックのすべての又は任意の所望の組の画像に
適用させることが可能である.1つのスタックの画像(
例えば、「ゴールデン装置」画像)と他のスタックの画
像(例えば「テスト装置」画像)の間のコントラスト及
び輝度における際は、スレッシュホー゜ルド値を単に変
化させることによって調節される. 以上、本発明の具体的実施の熊様について詳細に説明し
たが、本発明はこれら具体例にのみ限定されるべきもの
ではなく,本発明の技術的範囲を逸脱することなしに種
々の変形が可能であることはもちろんである.
【図面の簡単な説明】
第1図は電子ビームテストシステムを示した概略ブロッ
ク図、第2図は、第1図に示した電子ビームテストシス
テムにおいて使用されるデータ処理システムを示した概
略ブロック図、第3図ないし第5図はストロボ電圧コン
トラスト画像のスクリーン表示を示した各説明図,第6
図は本発明に基づき画像を整合する場合に行われるステ
ップを概略示した説明図、第7図は周囲の画素の強度値
から画素強度値を外挿することを示した説明図、第8図
は画素位置が合成画像の意図した表示フレム内にあるか
否かを決定する為に画素位置をテストする状熊を示した
説明図、第9図は画像間の際の視覚性を向上させる為の
3レベルスレッシュホールド操作の原理を示した説明図
、である.(符号の説明) lO:電子ビームテストプローブシステムl2:電子ビ
ーブテストプロープ 14:回路励起器 l6:データ処理システム l8:ディスプレーターミナル 22、24:バス 50:コンピュータシステム 60:ディスプレーインターフェース 61:高解像度カラーディスプレーターミナル 62:グラフィック発生器 64:インターフェース 66、通信バス FIG.イ

Claims (18)

    【特許請求の範囲】
  1. 1.第一画像から第二画像と画素毎に整合されている合
    成画像を形成する方法において、前記各画像が1組の装
    置に表示することが可能な画素を有しており且つ表示さ
    れた場合に前記画素によって画定される視覚的特徴を有
    しており、各画素が表示座標位置及び強度値によって画
    定され、 (a)前記第一画像の座標位置を前記合成画像のそれぞ
    れの座標位置へマップさせることが可能な浮動小数点値
    のマトリクスを計算し、 (b)前記浮動小数点値のマトリクスをそれぞれの固定
    小数点表示へ変換し、 (c)前記固定小数点表示を使用して、前記合成画像の
    各表示座標位置に対して前記第一画像の対応する座標位
    置を固定小数点演算で計算し、 (d)前記第一画像の前記対応する座標位置を取巻く画
    素の強度値から前記合成画像の各表示座標位置に対する
    画素強度値を外挿し、その際に前記合成画像を画定する
    1組の画素を発生する、上記各ステップを有しており、
    前記合成画像が前記第二画像の特徴と整合された特徴を
    有する前記第一画像の変形を有しており、前記合成画像
    の画素は前記第二画像の画素のそれぞれの表示座標位置
    と整合していることを特徴とする方法。
  2. 2.特許請求の範囲第一項において、前記第一及び第二
    画像は、電子ビームテストシステムによって発生された
    ストロボ電圧コントラスト画像であり、且つ集積回路動
    作状態を表わすものであることを特徴とする方法。
  3. 3.特許請求の範囲第1項において、前記ステップ(a
    )は、更に、前記第一及び第二画像をディスプレイスク
    リーン上に表示し、前記第一画像の3つの選択した特徴
    の座標位置を記録し、且つ前記第一画像の前記選択した
    特徴が前記合成画像と整合される前記第二画像の3つの
    選択した特徴の座標位置を記録する、上記各ステップを
    有しており、前記マトリクスが前記記録された座標位置
    から計算されることを特徴とする方法。
  4. 4.特許請求の範囲第3項において、前記選択した位置
    を表わす前記値を副画素精度で記録することを特徴とす
    る方法。
  5. 5.特許請求の範囲第1項において、前記マトリクスが
    、前記第一画像と第二画像との間の回転、傾斜及び拡大
    歪の補正の為の係数を表わす2×2にマトリクスの値を
    有していることを特徴とする方法。
  6. 6.特許請求の範囲第1項において、前記ステップ(b
    )が、前記固定小数点表示をレジスター内に格納するス
    テップを有しており、且つ前記ステップ(c)がマイク
    ロプロセサー装置内の前記第一画像の前記対応する座標
    位置を計算するステップを有することを特徴とする方法
  7. 7.特許請求の範囲第1項において、前記ステップ(d
    )が、前記周囲の画素に対する前記対応する座標位置の
    距離を表わす浮動小数点値を計算し、前記距離を固定小
    数点値へ変換し、且つ前記合成画像画素値が前記周囲の
    画素強度値の平均であって前記距離値によって重みずけ
    されている平均を有するように固定小数点演算を使用す
    るプロセサー内において前記固定小数点距離値及び前記
    周囲の画素強度値から前記合成画像画素値を計算する、
    上記各ステップを有することを特徴とする方法。
  8. 8.特許請求の範囲第1項において、更に、各々の値が
    前記合成画像画素強度値に対する1個の前記周囲の画素
    の重み付けした寄与度を表わす固定小数点値のテーブル
    を計算し、且つ前記テーブルの値を格納する、上記各ス
    テップを有しており、且つ前記ステップ(d)が、更に
    、前記各周囲の画素に対する前記テーブルから値を検索
    し且つ前記検索した値を結合して前記合成画像画素強度
    値を形成する、上記各ステップを有することを特徴とす
    る方法。
  9. 9.特許請求の範囲第8項において、前記テーブルの各
    格納した値が格納したインデックス値と関連しており、
    且つ前記テーブルから値を検索する前記ステップが、少
    なくとも一個の前記周囲の画素の値から及び前記対応す
    る座標位置に対する前記少なくとも一個の周囲の画素の
    距離からインデックス値を計算し、且つ前記計算したイ
    ンデックス値を前記格納したインデックス値とマッチン
    グさせて前記テーブルから検索されるべき値を選択する
    、上記各ステップを有することを特徴とする方法。
  10. 10.第一画像から第二画像と整合される合成画像を形
    成する装置において、前記各画像が格納装置内に格納さ
    れる一組の画素を有しており、且つ視覚的ディスプレイ
    ターミナル上に表示することが可能であり且つ表示され
    た場合に前記画素によって画定される視覚的特徴を有し
    ており、各画素は表示座標位置及び強度値によって画定
    され、 (a)前記第一画像の座標位置を前記合成画像のそれぞ
    れの座標位置へマップすることが可能な浮動小数点値か
    らなるマトリクスを計算する手段、 (b)前記浮動小数点値からなるマトリクスをそれぞれ
    の固定小数点表示へ変換させる手段、 (c)前記固定小数点表示を使用して、前記合成画像の
    各表示座標位置に対して前記第一画像の対応する座標位
    置を固定小数点演算で計算する手段、 (d)前記合成画像を画定する一組の画素を発生させる
    為に前記第一画像の前記対応する座標位置を取巻く画素
    の強度値から前記合成画像の各表示座標位置に対し画素
    強度値を外挿する手段、を有しており、前記合成画像を
    画定する前記1組の画素が前記第二画像の特徴と整合し
    た特徴を有する前記第一画像の変形を有しており、前記
    合成画像の画素が前記第二画像の画素のそれぞれの座標
    表示位置と整合されていることを特徴とする装置。
  11. 11.特許請求の範囲第10項において、更に、集積回
    路動作状態を表わすストロボ電圧コントラスト画像とし
    て前記第一及び第二画像を得る為の電子ビームプローブ
    手段を有していることを特徴とする装置。
  12. 12.特許請求の範囲第10項において、更に、前記第
    一及び第二画像をディスプレイスクリーン上に表示し、
    前記第一画像の3つの選択した特徴の座標位置を記録し
    、前記合成画像において前記第一画像の前記選択した特
    徴と整合される前記第二画像の3つの選択した特徴の座
    標位置を記録し、且つ前記記録した座標位置から前記浮
    動小数点値のマトリクスを計算することが可能であるこ
    とを特徴とする装置。
  13. 13.特許請求の範囲第12項において、前記記録手段
    が副画素精度で前記選択した位置を表わす値を記録する
    ことが可能であることを特徴とする装置。
  14. 14.特許請求の範囲第10項において、前記マトリク
    スが、前記第一及び第二画像の間の回転、傾斜及び拡大
    歪の補正の為の係数を表わす2×2マトリクスの値を有
    することを特徴とする装置。
  15. 15.特許請求の範囲第10項において、更に、前記固
    定小数点表示を格納する手段を有することを特徴とする
    装置。
  16. 16.特許請求の範囲第10項において、前記周囲の画
    素に対する前記対応する座標位置の距離を表わす浮動小
    数点値を計算し、前記距離を固定小数点値へ変換し、且
    つ前記合成画像画素値が前記周囲の画素強度値の平均で
    あって前記距離値によって重み付けされている平均を有
    するように固定小数点演算を使用して前記固定小数点距
    離値及び前記周囲の画素強度値から前記合成画像画素値
    を計算することが可能であることを特徴とする装置。
  17. 17.特許請求の範囲第10項において、各々の値が前
    記合成画像画素強度値に対する1個の前記周囲の画素の
    重み付けした寄与度を表わしている固定小数点値のテー
    ブルを計算し、前記テーブルの値を格納し、且つ各前記
    周囲の画素に対する前記テーブルから値を検索し且つ前
    記検索した値を結合して前記合成画像画素強度値を形成
    することが可能であることを特徴とする装置。
  18. 18.特許請求の範囲第17項において、前記テーブル
    の各格納された値が格納されたインデックス値と関連し
    ており、前記テーブルから値を検索し、少なくとも1つ
    の前記周囲の画素の値から及び前記対応する座標位置に
    対しての前記少なくとも1個の周囲の画素の距離からイ
    ンデックス値を計算し、且つ前記計算したインデックス
    値を前記格納したインデックス値とマッチングさせて前
    記テーブルから検索されるべき値を選択することが可能
    であることを特徴とする装置。
JP1303486A 1988-11-23 1989-11-24 合成画像を形成する方法及び装置 Expired - Fee Related JP2581815B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US276,359 1988-11-23
US07/276,359 US5054097A (en) 1988-11-23 1988-11-23 Methods and apparatus for alignment of images

Publications (2)

Publication Number Publication Date
JPH03162685A true JPH03162685A (ja) 1991-07-12
JP2581815B2 JP2581815B2 (ja) 1997-02-12

Family

ID=23056339

Family Applications (1)

Application Number Title Priority Date Filing Date
JP1303486A Expired - Fee Related JP2581815B2 (ja) 1988-11-23 1989-11-24 合成画像を形成する方法及び装置

Country Status (4)

Country Link
US (1) US5054097A (ja)
EP (1) EP0370322B1 (ja)
JP (1) JP2581815B2 (ja)
DE (1) DE68928726T2 (ja)

Families Citing this family (78)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0504944B1 (en) * 1991-03-22 1998-09-23 Nec Corporation Method of analyzing fault using electron beam
WO1993012501A1 (en) * 1991-12-10 1993-06-24 Logitech, Inc. Apparatus and methods for automerging images
JP2940317B2 (ja) * 1992-06-24 1999-08-25 三菱電機株式会社 画像処理装置
US5748195A (en) * 1992-10-29 1998-05-05 International Business Machines Corporation Method and means for evaluating a tetrahedral linear interpolation function
US5432892A (en) * 1992-11-25 1995-07-11 International Business Machines Corporation Volummetric linear interpolation
US5751926A (en) * 1992-12-23 1998-05-12 International Business Machines Corporation Function approximation using a centered cubic packing with tetragonal disphenoid extraction
US5390035A (en) * 1992-12-23 1995-02-14 International Business Machines Corporation Method and means for tetrahedron/octahedron packing and tetrahedron extraction for function approximation
US5401972A (en) * 1993-09-02 1995-03-28 Schlumberger Technologies, Inc. Layout overlay for FIB operations
US6061467A (en) * 1994-05-02 2000-05-09 Cognex Corporation Automated optical inspection apparatus using nearest neighbor interpolation
IL110771A (en) * 1994-08-25 1998-02-08 Holon Holon Hanni A priori and adaptive filtering for detection of signals corrupted by noise
JP3515199B2 (ja) * 1995-01-06 2004-04-05 大日本スクリーン製造株式会社 欠陥検査装置
US5541411A (en) * 1995-07-06 1996-07-30 Fei Company Image-to-image registration focused ion beam system
US6137494A (en) * 1995-08-18 2000-10-24 International Business Machines Corporation Method and means for evaluating a tetrahedral linear interpolation function
US5818951A (en) * 1995-10-13 1998-10-06 Infrared Service Corporation Methods and related apparatus for generating thermographic survey images
US5801719A (en) * 1995-11-27 1998-09-01 Sun Microsystems, Inc. Microprocessor with graphics capability for masking, aligning and expanding pixel bands
US5740278A (en) * 1996-02-16 1998-04-14 Cornell Research Foundation, Inc. Facsimile-based video compression method and system
US5917961A (en) * 1996-10-30 1999-06-29 Hewlett-Packard Co. Image convolution using pre-calculated lookup tables
US6330354B1 (en) * 1997-05-01 2001-12-11 International Business Machines Corporation Method of analyzing visual inspection image data to find defects on a device
US6236222B1 (en) * 1997-11-19 2001-05-22 Philips Electronics North America Corp. Method and apparatus for detecting misalignments in interconnect structures
US6445399B1 (en) * 1998-01-20 2002-09-03 Rafael - Armament Development Authority Ltd. System and method of visual orientation
US6115643A (en) * 1998-02-03 2000-09-05 Mcms Real-time manufacturing process control monitoring method
US6690840B1 (en) 2000-02-08 2004-02-10 Tektronix, Inc. Image alignment with global translation and linear stretch
WO2002040977A2 (en) * 2000-11-17 2002-05-23 Molecular Diagnostics, Inc. Evaluation of microscope slides
US7392287B2 (en) 2001-03-27 2008-06-24 Hemisphere Ii Investment Lp Method and apparatus for sharing information using a handheld device
EP1495446B1 (en) 2002-04-05 2007-10-10 The General Hospital Corporation Systems and method for generating an image
US6943569B1 (en) * 2002-04-12 2005-09-13 Advanced Micro Devices, Inc. Method, system and apparatus to detect defects in semiconductor devices
JP2004021495A (ja) * 2002-06-14 2004-01-22 Mitsubishi Electric Corp 監視システムおよび監視方法
EP1496681A1 (en) * 2003-07-11 2005-01-12 Axalto S.A. A method for printing an image on a support
GB2405764A (en) * 2003-09-04 2005-03-09 Sharp Kk Guided capture or selection of stereoscopic image pairs.
US7160741B2 (en) * 2003-11-06 2007-01-09 Chartered Semiconductor Manufacturing Ltd. Planar voltage contrast test structure and method
JP4904034B2 (ja) 2004-09-14 2012-03-28 ケーエルエー−テンカー コーポレイション レチクル・レイアウト・データを評価するための方法、システム及び搬送媒体
US7664315B2 (en) * 2004-11-03 2010-02-16 Tyzx, Inc. Integrated image processor
US20060098897A1 (en) * 2004-11-10 2006-05-11 Agfa-Gevaert Method of superimposing images
US7383994B2 (en) * 2005-05-03 2008-06-10 Datalogic Scanning, Inc. Methods and systems for forming images of moving optical codes
JP4591211B2 (ja) * 2005-05-31 2010-12-01 富士ゼロックス株式会社 画像処理装置、画像処理方法、媒体、符号読取装置、及びプログラム
US7676077B2 (en) 2005-11-18 2010-03-09 Kla-Tencor Technologies Corp. Methods and systems for utilizing design data in combination with inspection data
US8041103B2 (en) * 2005-11-18 2011-10-18 Kla-Tencor Technologies Corp. Methods and systems for determining a position of inspection data in design data space
WO2007120280A2 (en) * 2005-11-18 2007-10-25 Kla-Tencor Technologies Corporation Methods and systems for utilizing design data in combination with inspection data
US7570796B2 (en) 2005-11-18 2009-08-04 Kla-Tencor Technologies Corp. Methods and systems for utilizing design data in combination with inspection data
US7873233B2 (en) * 2006-10-17 2011-01-18 Seiko Epson Corporation Method and apparatus for rendering an image impinging upon a non-planar surface
US7877722B2 (en) 2006-12-19 2011-01-25 Kla-Tencor Corp. Systems and methods for creating inspection recipes
US8194968B2 (en) 2007-01-05 2012-06-05 Kla-Tencor Corp. Methods and systems for using electrical information for a device being fabricated on a wafer to perform one or more defect-related functions
US8990681B2 (en) * 2007-01-24 2015-03-24 Bluebeam Software, Inc. Method for aligning a modified document and an original document for comparison and difference highlighting
US7962863B2 (en) 2007-05-07 2011-06-14 Kla-Tencor Corp. Computer-implemented methods, systems, and computer-readable media for determining a model for predicting printability of reticle features on a wafer
US7738093B2 (en) 2007-05-07 2010-06-15 Kla-Tencor Corp. Methods for detecting and classifying defects on a reticle
US8213704B2 (en) 2007-05-09 2012-07-03 Kla-Tencor Corp. Methods and systems for detecting defects in a reticle design pattern
US7796804B2 (en) 2007-07-20 2010-09-14 Kla-Tencor Corp. Methods for generating a standard reference die for use in a die to standard reference die inspection and methods for inspecting a wafer
US7711514B2 (en) 2007-08-10 2010-05-04 Kla-Tencor Technologies Corp. Computer-implemented methods, carrier media, and systems for generating a metrology sampling plan
US7975245B2 (en) 2007-08-20 2011-07-05 Kla-Tencor Corp. Computer-implemented methods for determining if actual defects are potentially systematic defects or potentially random defects
US8139844B2 (en) 2008-04-14 2012-03-20 Kla-Tencor Corp. Methods and systems for determining a defect criticality index for defects on wafers
KR101623747B1 (ko) 2008-07-28 2016-05-26 케이엘에이-텐코어 코오포레이션 웨이퍼 상의 메모리 디바이스 영역에서 검출된 결함들을 분류하기 위한 컴퓨터-구현 방법들, 컴퓨터-판독 가능 매체, 및 시스템들
US8781219B2 (en) 2008-10-12 2014-07-15 Fei Company High accuracy beam placement for local area navigation
EP2332162A4 (en) 2008-10-12 2013-03-20 Fei Co HIGHLY PRECISE BAR PLACEMENT FOR NAVIGATION IN LOCAL AREAS
US8775101B2 (en) 2009-02-13 2014-07-08 Kla-Tencor Corp. Detecting defects on a wafer
US8204297B1 (en) 2009-02-27 2012-06-19 Kla-Tencor Corp. Methods and systems for classifying defects detected on a reticle
US8112241B2 (en) 2009-03-13 2012-02-07 Kla-Tencor Corp. Methods and systems for generating an inspection process for a wafer
JP5683839B2 (ja) 2010-05-17 2015-03-11 セミコンダクター・コンポーネンツ・インダストリーズ・リミテッド・ライアビリティ・カンパニー 撮像装置の制御回路
US8781781B2 (en) 2010-07-30 2014-07-15 Kla-Tencor Corp. Dynamic care areas
US9602697B2 (en) * 2010-08-27 2017-03-21 Ricoh Company, Ltd. Color substitution mechanism
US8455822B2 (en) 2010-08-31 2013-06-04 Fei Company Navigation and sample processing using an ion source containing both low-mass and high-mass species
US9170211B2 (en) 2011-03-25 2015-10-27 Kla-Tencor Corp. Design-based inspection using repeating structures
US9087367B2 (en) 2011-09-13 2015-07-21 Kla-Tencor Corp. Determining design coordinates for wafer defects
US8831334B2 (en) 2012-01-20 2014-09-09 Kla-Tencor Corp. Segmentation for wafer inspection
US8826200B2 (en) 2012-05-25 2014-09-02 Kla-Tencor Corp. Alteration for wafer inspection
US9189844B2 (en) 2012-10-15 2015-11-17 Kla-Tencor Corp. Detecting defects on a wafer using defect-specific information
US9053527B2 (en) 2013-01-02 2015-06-09 Kla-Tencor Corp. Detecting defects on a wafer
US9134254B2 (en) 2013-01-07 2015-09-15 Kla-Tencor Corp. Determining a position of inspection system output in design data space
US9311698B2 (en) 2013-01-09 2016-04-12 Kla-Tencor Corp. Detecting defects on a wafer using template image matching
WO2014149197A1 (en) 2013-02-01 2014-09-25 Kla-Tencor Corporation Detecting defects on a wafer using defect-specific and multi-channel information
US9865512B2 (en) 2013-04-08 2018-01-09 Kla-Tencor Corp. Dynamic design attributes for wafer inspection
US9310320B2 (en) 2013-04-15 2016-04-12 Kla-Tencor Corp. Based sampling and binning for yield critical defects
US9265463B1 (en) 2014-12-22 2016-02-23 Medical Metrics, Inc. Methods for determining spine instability and for eliminating the impact of patient effort on stability determinations
JP7294806B2 (ja) * 2015-03-23 2023-06-20 テックインサイツ インコーポレイテッド イメージング装置における歪み補正に関する方法、システム及び装置
US9542732B2 (en) * 2015-04-03 2017-01-10 Cognex Corporation Efficient image transformation
US10275863B2 (en) 2015-04-03 2019-04-30 Cognex Corporation Homography rectification
WO2018170416A1 (en) * 2017-03-16 2018-09-20 Vid Scale, Inc. Floating point to integer conversion for 360-degree video projection format conversion and spherical metrics calculation
US20220084223A1 (en) * 2020-09-14 2022-03-17 The Regents Of The University Of Michigan Focal Stack Camera As Secure Imaging Device And Image Manipulation Detection Method
WO2022187464A1 (en) * 2021-03-03 2022-09-09 Battelle Memorial Institute Systems and methods for inspection of ic devices

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS57138685A (en) * 1981-02-23 1982-08-27 Hitachi Ltd Graphic conversion for graphic indicator
US4610026A (en) * 1982-04-30 1986-09-02 Hitachi, Ltd. Method of and apparatus for enlarging/reducing two-dimensional images
US4578812A (en) * 1982-12-01 1986-03-25 Nec Corporation Digital image processing by hardware using cubic convolution interpolation
JPS59137942A (ja) * 1983-01-28 1984-08-08 Hitachi Ltd 画像位置合わせ方式
JPS6240146A (ja) * 1985-08-14 1987-02-21 Mitsubishi Electric Corp 荷電ビ−ムパタ−ン欠陥検査装置
US4721909A (en) * 1985-08-16 1988-01-26 Schlumberger Technology Corporation Apparatus for pulsing electron beams
US4706019A (en) * 1985-11-15 1987-11-10 Fairchild Camera And Instrument Corporation Electron beam test probe system for analyzing integrated circuits
US4760548A (en) * 1986-06-13 1988-07-26 International Business Machines Corporation Method and apparatus for producing a curve image
US4805123B1 (en) * 1986-07-14 1998-10-13 Kla Instr Corp Automatic photomask and reticle inspection method and apparatus including improved defect detector and alignment sub-systems

Also Published As

Publication number Publication date
EP0370322A2 (en) 1990-05-30
DE68928726D1 (de) 1998-08-13
DE68928726T2 (de) 1999-04-08
JP2581815B2 (ja) 1997-02-12
EP0370322A3 (en) 1991-09-25
EP0370322B1 (en) 1998-07-08
US5054097A (en) 1991-10-01

Similar Documents

Publication Publication Date Title
JPH03162685A (ja) 合成画像を形成する方法及び装置
US5127064A (en) High resolution image compression methods and apparatus
US5604819A (en) Determining offset between images of an IC
Franke et al. Measurement of microdisplacements by machine vision photogrammetry (DISMAP)
EP0504944B1 (en) Method of analyzing fault using electron beam
JP3798842B2 (ja) ビデオ・データにシグニチャ解析を実施する装置および方法
JPH0776974B2 (ja) 多次元視覚解析方法及び装置
JPH06325181A (ja) パターン認識方法
US5265040A (en) Method and device for displaying information on simulation result in a numerical simulation system
US6442733B1 (en) Failure analyzing method and apparatus using two-dimensional wavelet transforms
JP3924796B2 (ja) パターン位置の計測方法および計測装置
Collin et al. Device testing and SEM testing tools
US6606394B1 (en) Surface inspecting apparatus and method therefor
US5150461A (en) Fixed point method of video display scaling
JP3642904B2 (ja) フェイルカウント方法及びその装置
JP2003510568A (ja) パターン比較によるlcd検査方法およびlcd検査装置
EP0391587A2 (en) Method and apparatus for fast generation of warping coefficients
JP2645864B2 (ja) 集積回路試験装置
JPH06187800A (ja) 記憶素子の不良モード解析装置
CN113870255B (zh) Mini LED产品缺陷检测方法及相关设备
JP3237601B2 (ja) メモリlsi検査装置及びメモリlsi検査方法
Pelissier et al. Fully-integrated Dynamic Fault Imaging system for failure analysis and performance enhancement of VLSI
Burton et al. Digital image storage and processing for scanning microscopy
Ning et al. Automatic analysis of moiré fringe patterns by using an image-processing system
JPH045581A (ja) 回路の不良箇所表示方式

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees