JP2000155051A - Phase state analyzing method for equiphase stripes - Google Patents

Phase state analyzing method for equiphase stripes

Info

Publication number
JP2000155051A
JP2000155051A JP32863398A JP32863398A JP2000155051A JP 2000155051 A JP2000155051 A JP 2000155051A JP 32863398 A JP32863398 A JP 32863398A JP 32863398 A JP32863398 A JP 32863398A JP 2000155051 A JP2000155051 A JP 2000155051A
Authority
JP
Japan
Prior art keywords
phase
fringe
fringes
phase distribution
data
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
JP32863398A
Other languages
Japanese (ja)
Other versions
JP3497393B2 (en
Inventor
Nobuaki Ueki
伸明 植木
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.)
Fujinon Corp
Original Assignee
Fuji Photo Optical Co 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 Fuji Photo Optical Co Ltd filed Critical Fuji Photo Optical Co Ltd
Priority to JP32863398A priority Critical patent/JP3497393B2/en
Publication of JP2000155051A publication Critical patent/JP2000155051A/en
Application granted granted Critical
Publication of JP3497393B2 publication Critical patent/JP3497393B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Instruments For Measurement Of Length By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

PROBLEM TO BE SOLVED: To improve the precision of a phase state analysis in an observation region by applying a masking process to portions corresponding to regions where no stripe information exists in the phase distribution applied with a phase unlapping process in the analyzing method for the phase state of a checked body having an observation region split into multiple portions by separation bands where no stripe information exists. SOLUTION: In this phase state analyzing method, stripes scanning is applied in four steps to a checked face split into multiple observation regions by separation bands (S1). The phase distribution of the checked face is obtained from four interference stripes intensities obtained by the stripes scanning (S2). The phase distribution is folded between -π and π, thus phase unlapping is applied by the phase unlapping method, and the continuous phase distribution is obtained (S3). A masking process is applied to portions having 0 modulation corresponding to the separation bands (S4). Finally, an inclination correction is made to the continuous phase distribution applied with the masking process by the method of least squares (S5).

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、光反射体表面の形
状あるいは光透過体の厚みムラや屈折率分布を干渉縞あ
るいはモアレ縞を用いて解析測定する際の等位相縞の位
相状態解析方法に関し、特に、被検体表面が、互いに略
面一の状態に設定された複数の観察領域からなる場合、
あるいは被検体内部の屈折率分布が略均一に設定された
複数の観察領域からなる場合において、領域分割のため
の、縞情報が得られない分離帯が存在しても、この複数
の観察領域全体に亘る位相状態を解析測定し得る等位相
縞の位相状態解析方法に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method for analyzing the phase state of equiphase fringes when analyzing and measuring the shape of the surface of a light reflector or the thickness unevenness or refractive index distribution of a light transmitting body using interference fringes or moiré fringes. Regarding, particularly, when the subject surface is composed of a plurality of observation areas set to be substantially flush with each other,
Alternatively, in the case where a plurality of observation regions in which the refractive index distribution inside the subject is set to be substantially uniform, even if there is a separation band for which stripe information cannot be obtained for region division, the plurality of observation regions are The present invention relates to a phase state analysis method for equiphase fringes, which can analyze and measure a phase state over a range.

【0002】[0002]

【従来の技術】従来より、干渉縞等の等高線縞を解析し
て被検体の表面形状を高精度に計測する手法が知られて
いる。
2. Description of the Related Art Conventionally, there has been known a method of analyzing contour fringes such as interference fringes and measuring the surface shape of a subject with high accuracy.

【0003】例えば、IC基板として利用されるシリコ
ンウエハを吸着支持して、ウエハの表面チェックや表面
加工を行なうための吸着チャック装置においては試料を
吸着支持する基準面を光波長オーダーの面精度で加工し
なければならず、この基準面の面精度を測定するために
反射型の干渉縞計測による手法が用いられる。
For example, in a suction chuck apparatus for suction-supporting a silicon wafer used as an IC substrate and performing surface checking and surface processing of the wafer, a reference surface for suction-supporting a sample is provided with surface accuracy on the order of light wavelength. In order to measure the surface accuracy of the reference surface, a technique based on reflection-type interference fringe measurement is used.

【0004】ところで、この吸着チャック装置は、図7
で示される如き形状とされており、台座本体1の上面に
は面精度の高い基準面2が形成され、その面内には複数
の吸引孔3が平均的に分布され、この吸引孔3における
吸引作用により、例えば、薄板状のシリコンウエハ4が
この基準面2上に吸着保持されるようになっている。
By the way, this suction chuck device is shown in FIG.
A reference surface 2 with high surface accuracy is formed on the upper surface of the pedestal main body 1, and a plurality of suction holes 3 are evenly distributed in the surface. By the suction action, for example, a thin silicon wafer 4 is sucked and held on the reference surface 2.

【0005】このような基準面2の面精度を干渉縞計測
による手法を用いて計測する場合には、基準面2の全領
域における干渉縞情報を取り込んで(必要であれば縞走
査を行なって)、この取り込んだ情報に基づき所定の演
算を行なって位相分布(高さデータ)を得ることにな
る。
When the surface accuracy of the reference plane 2 is measured using the technique based on interference fringe measurement, interference fringe information in the entire area of the reference plane 2 is acquired (if necessary, fringe scanning is performed. ), A predetermined operation is performed based on the acquired information to obtain a phase distribution (height data).

【0006】ところが、上記した所定の演算によって
は、得られる位相値が(-π,π)の主値の範囲に折り
畳まれるため(位相ラッピング)、ダイナミックレンジ
の大きな位相に対しては位相値が不連続となり2πの整
数倍の不確定値を有するものとなる。したがって、実際
の表面形状に沿した位相分布を得るためには、このよう
に(-π,π)の主値の範囲に折り畳まれた位相分布Φ
(x,y)から元の連続な位相分布を求める位相アンラ
ッピング処理を施す必要がある。
However, since the obtained phase value is folded into the range of the main value of (-π, π) (phase wrapping) according to the above-mentioned predetermined calculation, the phase value is large for a phase having a large dynamic range. It becomes discontinuous and has an uncertain value that is an integral multiple of 2π. Therefore, in order to obtain the phase distribution along the actual surface shape, the phase distribution Φ folded in the range of the main value of (−π, π) as described above is obtained.
It is necessary to perform a phase unwrapping process for obtaining an original continuous phase distribution from (x, y).

【0007】このような位相アンラッピング処理を行な
う場合、ノイズの少ない滑らかな位相分布については単
純なアルゴリズムを用いて高精度な連続位相分布Φ
(x,y)を求めることが可能であるが、ノイズが多
く、モジュレーション(干渉縞振幅に依存する量)が低
い領域を含む場合には極めて難解な問題となり、連続位
相分布Φ(x,y)を求めることが困難となる。
When such a phase unwrapping process is performed, a highly accurate continuous phase distribution Φ is obtained using a simple algorithm for a smooth phase distribution with little noise.
Although it is possible to obtain (x, y), it becomes a very difficult problem when a region including a large amount of noise and low modulation (an amount depending on the interference fringe amplitude) is included, and the continuous phase distribution Φ (x, y) is obtained. ) Is difficult to obtain.

【0008】このような事情の下、近年になって、ノイ
ズが多かったりモジュレーションが低い場合において
も、位相アンラッピング処理を良好に行なうことができ
る手法が提案されている。
[0008] Under such circumstances, in recent years, there has been proposed a technique capable of performing a phase unwrapping process satisfactorily even when noise is large or modulation is low.

【0009】この手法は干渉縞のコントラストがより良
い(モジュレーションがより高い)部分に着目し、この
ような部分から順に位相アンラッピングしていく方法
で、計算機科学の分野で最小木問題(Minimum Spanning
Tree Problem)と称されるグラフ問題を応用したもの
であり、特にそれを用いた振幅最大木法が有効なものと
して知られている(第55回(1994年)応用物理学会学
術講演会予稿集P803参照)。
This method focuses on a portion where interference fringes have higher contrast (higher modulation) and sequentially unwraps the phase from such a portion. In the field of computer science, the minimum tree problem (Minimum Spanning) is used.
Tree Problem) is an application of a graph problem, and the maximum-amplitude tree method using it is known to be particularly effective (The 55th (1994) Proceedings of the Japan Society of Applied Physics) P803).

【0010】[0010]

【発明が解決しようとする課題】しかしながら、図7に
示す吸着チャック装置等においては、所定の処理が終了
した後において、過大な力を要することなくウエハ4を
基準面2から除去することができるように基準面2が格
子状の溝5によって複数の領域に分割されている。
However, in the suction chuck device and the like shown in FIG. 7, the wafer 4 can be removed from the reference surface 2 without requiring excessive force after the predetermined processing is completed. As described above, the reference surface 2 is divided into a plurality of regions by the lattice-shaped grooves 5.

【0011】このように溝5により観察領域が完全に複
数の領域に分割されている場合には、この溝5から有効
な干渉縞情報が得られないため、有効な干渉縞情報が得
られる観察領域が、干渉縞データ上において各々完全に
分割された領域となるから、基準面2に対する干渉縞情
報に基づき、この干渉縞情報が無い領域を含んだ状態で
所定の演算を行なって得られた位相分布に位相アンラッ
ピング処理を施すと解析データに誤差が生じるという問
題があった。
In the case where the observation area is completely divided into a plurality of areas by the groove 5, effective interference fringe information cannot be obtained from the groove 5, so that observation interference information can be obtained. Since the area is a completely divided area on the interference fringe data, the area is obtained by performing a predetermined calculation based on the interference fringe information with respect to the reference plane 2 while including the area without the interference fringe information. When the phase unwrapping process is performed on the phase distribution, there is a problem that an error occurs in the analysis data.

【0012】すなわち、上述した位相アンラッピング手
法を用いた場合においては、上記基準面2における、干
渉縞が生じない溝5に相当する部分がどうしても誤差要
因となり、基準面2の全体形状を正確に解析することは
難しい。
That is, in the case where the above-described phase unwrapping method is used, a portion corresponding to the groove 5 where the interference fringes do not occur on the reference surface 2 is an error factor, and the entire shape of the reference surface 2 is accurately determined. Difficult to analyze.

【0013】また、干渉縞計測の手法を用いて被検面の
表面形状を得る場合、被検体設定部に設定された被検面
はどうしても若干の傾斜がついてしまうため、この被検
面の形状データを得る場合には、最小二乗法等を用いて
その傾きを補正する処理が最終的に必要となるが、位相
アンラッピングを行なって得られたデータに対してこの
ような傾き補正処理を加えると、干渉縞が発生していな
い前述した段差部分がデータ全体に対して重みとなり正
確な傾き補正ができないという問題もある。
When the surface shape of the test surface is obtained by using the interference fringe measurement technique, the test surface set in the test object setting section is necessarily slightly inclined. When obtaining data, a process of correcting the inclination using the least squares method or the like is finally required. However, such a tilt correction process is added to data obtained by performing phase unwrapping. In addition, there is also a problem that the above-described step portion where no interference fringe occurs becomes a weight for the entire data and accurate tilt correction cannot be performed.

【0014】このような問題は、被検体が透過物体から
なり、その観察領域が、干渉縞情報を得ることが困難な
分離帯により複数個に分割されている場合に、この観察
領域内部の屈折率分布、さらにはこの観察領域の厚みム
ラを透過型の干渉縞計測による手法を用いて測定する場
合においても同様に生じることになる。
[0014] Such a problem is caused by the fact that a subject is composed of a transmissive object and its observation area is divided into a plurality of parts by a separation band in which it is difficult to obtain interference fringe information. This also occurs when the rate distribution and the thickness unevenness of the observation area are measured using a transmission interference fringe measurement technique.

【0015】本発明はこのような事情に鑑みなされたも
ので、等位相縞を解析して被検体の表面形状、屈折率分
布あるいは厚みムラを測定する場合において、該被検体
が分離帯により複数の観察領域に分割され、かつ該分離
帯の縞情報を得ることができない場合においても、その
分離帯におけるデータが観察領域の位相解析に影響を与
えることがなく、高精度で観察領域の位相状態測定を行
なうことができる等位相縞の位相状態解析方法を提供す
ることを目的とするものである。
The present invention has been made in view of such circumstances, and when measuring the surface shape, refractive index distribution, or thickness unevenness of an object by analyzing equiphase fringes, the object is separated by a plurality of separation bands. Is divided into observation regions, and even when the stripe information of the separation band cannot be obtained, the data in the separation band does not affect the phase analysis of the observation region. It is an object of the present invention to provide a phase state analysis method of equiphase fringes that can perform measurement.

【0016】[0016]

【課題を解決するための手段】本発明の等位相縞の位相
状態解析方法は、縞情報が存在しない分離帯と、この分
離帯により複数個に分割され、該分離帯を挟んで互いに
隣接する領域端部の位相差が、観察される等位相縞で1
/2縞よりも小さい縞情報値を有する観察領域とからな
る被検体の縞情報を取り込み、前記観察領域に前記分離
帯を含めて前記被検体の縞情報に基づく位相計算を行な
い、次いで、該位相計算により位相ラッピングされた位
相分布の位相アンラッピング処理を行ない、該位相アン
ラッピング処理が施された位相分布において、前記縞情
報が存在しない分離帯に対応する部分にマスキング処理
を施すことを特徴とするものである。
According to the method for analyzing the phase state of equiphase fringes according to the present invention, a separating band having no fringe information and a plurality of divided bands separated by the separating band are adjacent to each other with the separating band interposed therebetween. The phase difference at the end of the region is 1
The fringe information of the subject including the observation region having a fringe information value smaller than / 2 fringes is taken in, the phase is calculated based on the fringe information of the subject including the separation band in the observation region, and then, A phase unwrapping process is performed on a phase distribution that has been phase-wrapped by phase calculation, and a masking process is performed on a portion of the phase distribution that has been subjected to the phase unwrapping process, the portion corresponding to the separation band where the fringe information does not exist. It is assumed that.

【0017】また、マスキング処理が施された位相分布
に所定の関数補正処理を施すことを特徴とするものであ
る。
Further, the present invention is characterized in that a predetermined function correction process is performed on the phase distribution subjected to the masking process.

【0018】上記「分離帯を挟んで互いに隣接する領域
端部」とは、分離帯を挟んで互いに隣接する2つの観察
領域における、対向する分離帯境界部分を意味するもの
とする。
The above-mentioned "ends of the regions adjacent to each other with the separator in between" mean the boundary portions of the opposing separators in the two observation regions adjacent to each other with the separator in between.

【0019】[0019]

【発明の実施の形態】以下、本発明に係る実施形態につ
いて図面を参照しつつ具体的に説明する。
Embodiments of the present invention will be specifically described below with reference to the drawings.

【0020】図1は、本実施形態に係る等位相縞の位相
状態解析方法を説明するためのフローチャートである。
本実施形態方法においては、フィゾー型の反射型レーザ
干渉計を用いており、まず、ステップ1(S1)に示す
如く、分離帯によって複数の観察領域に分割された、光
反射体である被検面に対し、4ステップの縞走査(フリ
ンジスキャン)を行う。なお、分離帯は観察領域にピン
トを合わせた場合、干渉縞が生じないものとされてい
る。また、分離帯を挟んで隣接する2つの観察領域にお
ける、対向する分離体境界部分の位相差は1/2縞(垂
直入射の光波干渉計で入射光波長の1/4)よりも良好
であることが必要であり、さらに観察領域においてモジ
ュレーションの高い縞情報を得られることが必要であ
る。
FIG. 1 is a flowchart for explaining a phase state analysis method for equiphase fringes according to the present embodiment.
In the method of the present embodiment, a Fizeau-type reflection laser interferometer is used. First, as shown in step 1 (S1), a test object which is a light reflector divided into a plurality of observation regions by a separation band. A four-step fringe scan is performed on the surface. It should be noted that the separation band does not generate interference fringes when focused on the observation area. Further, the phase difference at the boundary between the opposing separators in the two observation regions adjacent to each other with the separation band interposed therebetween is better than 縞 fringe (の of the incident light wavelength by a vertically incident light wave interferometer). It is necessary to obtain fringe information with high modulation in the observation region.

【0021】次に、この縞走査により得られた4つの干
渉縞強度I、I、I、Iから被検面の位相分布
Φ(x,y)を求める(S2)。なお、この場合の位相
分布Φ(x,y)は平面内の一直線上(図2の白矢印
上)における縞情報に基づいて求められる。
Next, the phase distribution Φ (x, y) of the test surface is determined from the four interference fringe intensities I 1 , I 2 , I 3 and I 4 obtained by the fringe scanning (S2). In this case, the phase distribution Φ (x, y) is obtained based on stripe information on a straight line in the plane (on the white arrow in FIG. 2).

【0022】この位相分布Φ(x,y)は−πからπの
間に畳み込まれた位相ラッピングデータとされているの
で、前述した振幅最大木法等の位相アンラッピング手法
を用いて位相アンラッピングを行い、連続する位相分布
Φ(x,y)を求める(S3)。
Since this phase distribution Φ (x, y) is assumed to be phase wrapping data convolved between −π and π, the phase unwrapping method such as the above-mentioned maximum amplitude tree method is used. Lapping is performed to obtain a continuous phase distribution Φ (x, y) (S3).

【0023】次に、連続位相分布Φ(x,y)におい
て、上記分離帯に相当する、モジュレーションが0の部
分についてマスキング処理を行う(S4)。ここでマス
キング処理とは、例えば除去したいデータに対して0
を、そのまま残したいデータに対して1を各々乗じる演
算処理をいうが、不要なデータのみを除去することがで
きれば、その他の態様のマスキング処理によることも可
能である。
Next, in the continuous phase distribution Φ (x, y), a masking process is performed on a portion where the modulation is 0, which corresponds to the above-mentioned separation band (S4). Here, the masking process means that, for example, 0
Is multiplied by 1 with respect to the data to be left as it is. However, if only unnecessary data can be removed, masking processing in another mode can be used.

【0024】最後に、このマスキング処理を施された連
続位相分布Φ(x,y)に対して最小二乗法を用いて傾
き補正を行う(S5)。このときパワー補正、ザイデル
収差の関数での補正、ツェルニケ多項式での補正等の他
の関数補正を同時に行ってもよい。
Finally, inclination correction is performed on the masked continuous phase distribution Φ (x, y) using the least squares method (S5). At this time, other function corrections such as power correction, correction using a function of Seidel aberration, and correction using a Zernike polynomial may be performed at the same time.

【0025】また、上述したステップ1(S1)とステ
ップ2(S2)との間において、所望の観察領域のみを
抽出するマスキング処理を位相アンラッピング計算前に
行なうようにすれば、その後の計算処理時間の短縮化を
図ることができ、形状解析処理の高速化を図ることがで
きる。
If the masking process for extracting only a desired observation region is performed before the phase unwrapping calculation between the above-mentioned steps 1 (S1) and 2 (S2), the subsequent calculation process The time can be reduced, and the shape analysis processing can be speeded up.

【0026】以下、本実施形態について、より詳細に説
明する。
Hereinafter, the present embodiment will be described in more detail.

【0027】レーザ干渉計で縞走査法を用いた場合、例
えば4ステップ法だと、4回のフリンジスキャンにおけ
る干渉縞強度I、I、I,Iは以下の各式
(1)、(2)、(3)、(4)により表される。 I(x,y)=I(x,y)[1+γ(x,y)cos{Φ(x,y)}] ・・・(1) I(x,y)=I(x,y)[1+γ(x,y)cos{Φ(x,y)+π/2}] ・・・(2) I(x,y)=I(x,y)[1+γ(x,y)cos{Φ(x,y)+π}] ・・・(3) I(x,y)=I(x,y)[1+γ(x,y)cos{Φ(x,y)+3π/2}] ・・・(4)
When a fringe scanning method is used in a laser interferometer, for example, in the case of a four-step method, the interference fringe intensities I 1 , I 2 , I 3 , and I 4 in four fringe scans are represented by the following equations (1). , (2), (3), and (4). I 1 (x, y) = I 0 (x, y) [1 + γ (x, y) cos {Φ (x, y)}] (1) I 2 (x, y) = I 0 (x , Y) [1 + γ (x, y) cos {Φ (x, y) + π / 2}] (2) I 3 (x, y) = I 0 (x, y) [1 + γ (x, y ) cos {Φ (x, y) + π}] (3) I 4 (x, y) = I 0 (x, y) [1 + γ (x, y) cos {Φ (x, y) + 3π / 2}] ... (4)

【0028】ここで、x、yは座標、Φ(x,y)は位
相、I(x,y)は各点での平均光強度、γ(x,
y)は干渉縞のモジュレーションを各々表す。なお、モ
ジュレーションとは、干渉縞振幅と略相関関係を有する
ものと考えて良く、具体的には下記式(5)により表さ
れる。
Here, x and y are coordinates, Φ (x, y) is a phase, I 0 (x, y) is an average light intensity at each point, and γ (x, y)
y) represents the interference fringe modulation, respectively. The modulation may be considered to have a substantially correlation with the interference fringe amplitude, and is specifically expressed by the following equation (5).

【0029】[0029]

【数1】 (Equation 1)

【0030】上述した4つの式(1)、(2)、
(3)、(4)から、位相Φ(x,y)を求めると下記
式(6)となる。
The above four equations (1), (2),
When the phase Φ (x, y) is obtained from (3) and (4), the following equation (6) is obtained.

【0031】[0031]

【数2】 (Equation 2)

【0032】このとき、上記式(6)が逆正接関数とな
っていることからも明らかなように、位相Φ(x,y)
は−πからπの間に畳み込まれている。なお、このよう
な状態のデータを通常位相ラッピングデータと称してい
る。この式で分母が0の場合(すなわちI(x,y)
=I(x,y)の場合)は、分子が正であればπ/2
となり、負であれば−π/2となる。
At this time, as apparent from the fact that the above equation (6) is an arctangent function, the phase Φ (x, y)
Is convolved between -π and π. The data in such a state is usually called phase wrapping data. When the denominator is 0 in this formula (that is, I 1 (x, y)
= I 3 (x, y)) is π / 2 if the molecule is positive
And if negative, -π / 2.

【0033】また、干渉縞が存在しない部分(図2の十
字線部分)は、I(x,y)=I (x,y)=I
(x,y)=I(x,y)=0となるため、位相ラッ
ピングデータにおいてπ/2もしくは−π/2となる。
すなわち、干渉縞の存在しない領域でのデータは誤差と
なってしまう。
Further, a portion where no interference fringes exist (10 in FIG. 2)
Character part) is I1(X, y) = I 2(X, y) = I3
(X, y) = I4Since (x, y) = 0, the phase
It becomes π / 2 or -π / 2 in the ping data.
In other words, data in an area where no interference fringes exist is considered to be an error.
turn into.

【0034】しかし、このような干渉縞の存在しない領
域(分離帯)で分けられた両側の観察領域の、分離帯境
界部分の高さ(位相)が、略同じ高さ(少なくとも位相
差がπ未満、すなわち、He−Neレーザを用いたフィ
ゾー型干渉計の場合には、入射光の波長を考慮すると高
低差が0.16μm以下である。したがって、被検面の
傾き調整を十分に行う必要がある。)であれば、前述し
た振幅最大木法等の位相アンラッピング方法を用いて解
析を行った場合、干渉縞の存在する観察領域のみを解析
する場合には計算誤差は生じていない。換言すれば、干
渉縞の存在しない分離帯部分について計算誤差が生じて
いる。
However, the heights (phases) of the separation zone boundary portions of the observation areas on both sides divided by such an area where no interference fringes exist (separation zones) are substantially the same height (at least the phase difference is π). In other words, in the case of a Fizeau interferometer using a He-Ne laser, the height difference is 0.16 μm or less in consideration of the wavelength of incident light. If the analysis is performed using the above-described phase unwrapping method such as the maximum amplitude tree method, no calculation error occurs when only the observation region where the interference fringe exists is analyzed. In other words, a calculation error occurs in the separation band where no interference fringes exist.

【0035】そこで、このような事実を考慮に入れ、本
実施形態においては、干渉縞が存在しない分離帯と、干
渉縞の存在する複数の観察領域とを含んだ全領域につい
て上記位相アンラッピング方法による計算を行なった
後、モジュレーションが0の領域(この段階ではモジュ
レーション値の判断が容易である)についてマスキング
処理を行って上記干渉縞が存在しない領域のデータを排
除し、干渉縞の存在しない分離帯部分についての計算誤
差が上記観察領域の解析値に影響を与えないようにして
いる。
Taking this fact into consideration, in the present embodiment, the phase unwrapping method is applied to the entire region including the separation band where no interference fringes exist and the plurality of observation regions where the interference fringes exist. , The masking process is performed on the region where the modulation is 0 (the modulation value is easy to judge at this stage) to eliminate the data of the region where the interference fringes do not exist, and the separation where the interference fringes do not exist is performed. The calculation error for the band portion does not affect the analysis value of the observation region.

【0036】また、干渉縞計測の手法を用いて被検面の
表面形状を得る場合、被検体設定部に設定された被検面
はどうしても若干の傾斜がついてしまうため、この被検
面の実際の形状データを得る場合には、最小二乗法等を
用いてその傾きを補正する処理が最終的に必要となる。
しかし、位相アンラッピングを行なって得られた連続位
相分布に対してこのような傾き補正処理を加えると、干
渉縞が発生していない前述した分離帯に相当する部分が
データ全体に対して重みとなり正確な傾き補正ができな
いという問題があるが、本実施形態においては、モジュ
レーションが0となる領域についてマスキング処理を行
って上記干渉縞が存在しない領域のデータを排除してい
るので、最小二乗法により正確な傾き補正を行うことが
でき、正しい形状測定を実施することができる。なお、
干渉縞情報の無い分離帯を挟んで隣接する領域の位相差
が1/2縞より小さい傾きであることが条件となる。
Further, when the surface shape of the test surface is obtained by using the interference fringe measurement technique, the test surface set in the test object setting section is inevitably slightly inclined. When the shape data is obtained, a process of correcting the inclination using the least square method or the like is finally required.
However, when such a slope correction process is applied to the continuous phase distribution obtained by performing the phase unwrapping, a portion corresponding to the above-described separation band where no interference fringe is generated becomes a weight for the entire data. Although there is a problem that accurate tilt correction cannot be performed, in the present embodiment, masking processing is performed on an area where modulation is 0 to eliminate data in an area where the interference fringes do not exist. Accurate inclination correction can be performed, and correct shape measurement can be performed. In addition,
The condition is that the phase difference between the regions adjacent to each other across the separation band having no interference fringe information has a slope smaller than 1/2 fringe.

【0037】また、上記縞走査法においては4ステップ
法の代わりに、5ステップ法などの、3ステップ以上の
他のステップ法を用いて行うことも可能である。
Further, in the above fringe scanning method, it is possible to use another step method of three or more steps such as a five step method instead of the four step method.

【0038】〈実施例〉以下、具体的な実施例について
シミュレーションデータを、従来技術と対比しつつ解析
することにより説明する。
<Embodiments> Specific embodiments will be described below by analyzing simulation data in comparison with the prior art.

【0039】まず、被検面内に干渉縞の生じない2本の
十字ライン(分離帯)が存在して観察領域が分割され、
観察領域は完全に平面な状態とされている場合について
想定し、このような被検面を有する被検体を少し傾ける
(分離帯を挟んだ領域の高低差は1/2縞未満)ように
してフィゾー型のレーザ干渉計にセッティングした状態
について想定する。
First, the observation area is divided by two cross lines (separation zones) in which no interference fringes occur in the surface to be inspected.
It is assumed that the observation region is completely flat, and the subject having such a test surface is slightly tilted (the height difference in the region sandwiching the separation band is less than 1/2 fringe). It is assumed that the state is set in a Fizeau laser interferometer.

【0040】このように想定された被検面についての干
渉縞シミュレーションデータを図2に示す。なお、図2
の(A)は前述した式(1)に対応する干渉縞データ
を、(B)は前述した式(3)に対応する干渉縞データ
を示すものである。
FIG. 2 shows simulation data of interference fringes for the test surface assumed as described above. Note that FIG.
(A) shows the interference fringe data corresponding to the above equation (1), and (B) shows the interference fringe data corresponding to the above equation (3).

【0041】この干渉縞データに対し、位相ラッピング
データの中心横断面(白矢印線上)の表面形状は、図3
に示すような鋸歯状に形成される。なお、十字線に相当
する部分は干渉縞が存在せず、π/2ラジアンに設定さ
れるため、グラフ上に歪が形成され、これが解析誤差の
発生原因となる。なお、解析処理を行なうソフトプログ
ラムの設定によっては、上記π/2ラジアンとならずに
−π/2ラジアンとなる場合もある。
With respect to this interference fringe data, the surface shape of the center cross section (on the white arrow line) of the phase wrapping data is shown in FIG.
Is formed in a sawtooth shape as shown in FIG. Note that the portion corresponding to the crosshair has no interference fringes and is set to π / 2 radians, so that a distortion is formed on the graph, which causes an analysis error. Note that, depending on the setting of the software program for performing the analysis process, the value may be -π / 2 radians instead of the above-mentioned π / 2 radians.

【0042】次に、この図3に示す位相ラッピングデー
タを、前述した振幅最大木法により位相アンラッピング
すると、図4に表わされたグラフの如くなる。十字線に
相当する、干渉縞が存在しない部分は直線グラフ上で歪
となっていることが分かる。
Next, when the phase wrapping data shown in FIG. 3 is subjected to phase unwrapping by the above-described maximum amplitude tree method, a graph shown in FIG. 4 is obtained. It can be seen that the portion corresponding to the crosshair and having no interference fringes is distorted on the straight line graph.

【0043】次に、この図4で表わされたデータについ
て、最小二乗法により傾き補正を行うと、全領域に亘っ
て、誤差成分が重みとなり、図5に示す如く、正確な傾
き補正が困難となる。すなわち、図5においてグラフ直
線部分が所定の傾きをもってしまう。
Next, when the data shown in FIG. 4 is subjected to the inclination correction by the least square method, the error component becomes a weight over the entire area, and as shown in FIG. 5, accurate inclination correction is performed. It will be difficult. That is, the straight line portion of the graph has a predetermined slope in FIG.

【0044】そこで、最小二乗法により傾き補正を行な
う前に、モジュレーションが0となる部分(十字線に相
当する部分)のデータにマスキング処理を施すようにす
る。ここで、モジュレーションの値の変化を図6に示
す。マスキング処理は、モジュレーションの値が1であ
れば該データに1を、0であれば該データに0を各々乗
算する。
Therefore, before performing the inclination correction by the least-squares method, masking processing is performed on the data of the portion where the modulation is 0 (the portion corresponding to the cross hair). Here, a change in the value of the modulation is shown in FIG. The masking process multiplies the data by 1 if the modulation value is 1, and multiplies the data by 0 if the modulation value is 0.

【0045】このように、傾き補正を行う前にモジュレ
ーションが0となる領域にマスキング処理を施すことに
よって、傾き補正を行なった場合にもグラフ直線部分の
傾きを確実に補正することができる。
As described above, by performing the masking process on the area where the modulation becomes 0 before performing the inclination correction, the inclination of the straight line portion of the graph can be surely corrected even when the inclination is corrected.

【0046】以下、本実施形態のものにおいて、より良
い解析結果を得る上で注意すべき点について列挙する。
Hereinafter, points to be noted in obtaining a better analysis result in the embodiment will be enumerated.

【0047】(1)モジュレーションの良好な観察領
域、すなわち、基準面との反射率がほぼ等しく干渉縞の
コントラストが良好な観察領域であることが好ましい。
(1) It is preferable that the observation region has good modulation, that is, the observation region has almost the same reflectance with the reference surface and has good contrast of interference fringes.

【0048】(2)干渉計の照射光量を、干渉縞コント
ラストが常に最良となるような状態に調整しておくこと
が好ましい。
(2) It is preferable to adjust the irradiation light amount of the interferometer so that the interference fringe contrast always becomes the best.

【0049】(3)被検面に干渉計のピントが合ってい
ることが好ましい。ピントが合っていないと干渉縞が、
回折現象の影響からエッジ部分で流れたようになってし
まい誤差要因となる。
(3) It is preferable that the interferometer is focused on the surface to be measured. If not focused, interference fringes
Due to the effect of the diffraction phenomenon, it flows as if at the edge, which causes an error.

【0050】(4)空気の揺らぎや振動が極力小さい測
定環境で測定することが好ましい。このような要件が満
たされないと、フリンジスキャンが正確に行われず、測
定誤差が発生する他、干渉縞の存在する領域におけるモ
ジュレーションが低くなってしまう。
(4) It is preferable to perform measurement in a measurement environment in which air fluctuation and vibration are as small as possible. If such requirements are not satisfied, fringe scanning will not be performed accurately, causing measurement errors and lowering the modulation in the region where the interference fringes exist.

【0051】(5)被検面の平面度が良く、かつ、アラ
イメントは縞がもっとも少なくなる状態にして測定する
ことが好ましい。
(5) It is preferable that the measurement is performed in a state where the flatness of the surface to be inspected is good and the alignment is minimized.

【0052】なお、本発明の等位相縞の位相状態解析方
法としては上述した実施形態のものに限られるものでは
なく、その他の種々の態様の変更が可能である。例え
ば、上述した実施形態のものでは被検体が光反射体であ
る場合に、その表面形状(凹凸形状や段差形状)の測定
に適用するものであるが、本発明方法は被検体が光透過
体である場合に、その厚みムラや屈折率分布の測定にも
適用することができる。
The method for analyzing the phase state of equiphase fringes of the present invention is not limited to the above-described embodiment, and various other modifications can be made. For example, in the above-described embodiment, when the object is a light reflector, the method is applied to measurement of the surface shape (irregular shape or step shape) of the light reflector. In this case, the method can be applied to the measurement of the thickness unevenness and the refractive index distribution.

【0053】ただし、被検体が光透過体である場合に、
本発明方法を適用する際にも、上記実施形態と同様の前
提条件を満足することが必要である。
However, when the subject is a light transmitting body,
When applying the method of the present invention, it is necessary to satisfy the same prerequisites as in the above embodiment.

【0054】すなわち、分離帯によって分割された両観
察領域の対向する分離帯境界部分の光入射方向の位相差
(厚みムラや屈折率分布に伴なう位相差)が等位相縞で
1/2縞未満であることが必要であり、したがって、上
記対向する分離帯の高低差は、例えばフィゾー型やマイ
ケルソン型の干渉計を用いる場合にはλ/4未満、マッ
ハツェンダ型の干渉計ではλ/2未満であることが必要
である(ただし、λは光源からの光波長)。また、分離
帯には縞情報が存在しないことが必要であり、測定領域
の縞情報はコントラストのよいものであることが必要で
あり、さらに、傾きは分離帯を挟んでいる干渉縞の存在
する領域の位相差が1/2縞未満となることが必要であ
る。
That is, the phase difference in the light incident direction (the phase difference accompanying the thickness unevenness and the refractive index distribution) at the boundary between the opposite separation zones of the two observation regions divided by the separation zone is equal to 1/2 of the equal phase fringe. It is necessary to be less than the fringes. Therefore, the height difference between the opposing separation zones is, for example, less than λ / 4 when a Fizeau type or Michelson type interferometer is used, and λ / It is necessary to be less than 2 (where λ is the light wavelength from the light source). In addition, it is necessary that no stripe information exists in the separation band, the stripe information in the measurement area needs to have good contrast, and the inclination has the presence of interference fringes sandwiching the separation band. It is necessary that the phase difference of the region is less than 1/2 fringe.

【0055】さらに、本実施形態方法では、干渉計装置
として光反射型フィゾーのレーザ干渉計を用いている
が、これ以外の種々の干渉計を使用し得る。また、照射
光源としては、レーザ光源に代えて白色光源とすること
も可能である。
Furthermore, in the method of the present embodiment, a light reflection type Fizeau laser interferometer is used as the interferometer device, but various other interferometers may be used. Further, as the irradiation light source, a white light source may be used instead of the laser light source.

【0056】また、縞情報の種類としては上述した干渉
縞情報に限られるものではなく、モアレ縞情報をも対象
とし得る。ただし、モアレ縞の縞情報においては、1縞
に相当する量が干渉縞の縞情報とは異なるものであり、
例えば平行光照明を用いた垂直観察の格子照射法におい
ては、格子ピッチをP、照明角度をθとしたとき、1縞
はP/tanθの量に相当する。
The type of fringe information is not limited to the above-described interference fringe information, but may be moire fringe information. However, in the fringe information of the moire fringe, the amount corresponding to one fringe is different from the fringe information of the interference fringe.
For example, in a grid irradiation method for vertical observation using parallel light illumination, when a grid pitch is P and an illumination angle is θ, one fringe corresponds to the amount of P / tan θ.

【0057】さらに、本発明の等位相縞の位相状態解析
方法としては、縞走査法により得られた複数の縞情報に
基いて被検体の位相分布Φ(x,y)を得る場合に限ら
れず、1回の取込みによって得られた縞情報に基づき上
記位相分布Φ(x,y)を得る場合のものに適用するこ
とも可能である。
Further, the phase state analysis method of equiphase fringes of the present invention is not limited to the case where the phase distribution Φ (x, y) of the subject is obtained based on a plurality of pieces of fringe information obtained by the fringe scanning method. The present invention can also be applied to a case where the phase distribution Φ (x, y) is obtained based on the stripe information obtained by one capture.

【0058】[0058]

【発明の効果】以上説明したように、本発明の等位相縞
の位相状態解析方法によれば、位相アンラッピング処理
が施された位相分布において、縞情報が存在しない領域
に対応する部分にマスキング処理を施し、最終段階にお
いて分離帯に関するデータは除去されるようにしてお
り、この分離帯のデータが解析値全体に及ぼす影響を排
除できるから、被検体の観察領域における位相状態解析
の精度の向上を図ることができる。また、この後におい
て位相分布に傾き補正等の関数補正処理を施す場合に
は、その補正の精度を向上させることができる。
As described above, according to the phase state analysis method for equiphase fringes of the present invention, masking is performed on a portion corresponding to a region where no fringe information exists in the phase distribution subjected to the phase unwrapping process. In the final stage, the data on the separation zone is removed, and the influence of the data on the separation zone on the entire analysis value can be eliminated, thus improving the accuracy of the phase state analysis in the observation region of the subject. Can be achieved. Further, if a function correction process such as a tilt correction is performed on the phase distribution thereafter, the accuracy of the correction can be improved.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明の等位相縞の位相状態解析方法に係る実
施形態を説明するためのフローチャート
FIG. 1 is a flowchart for explaining an embodiment according to a phase state analysis method of equiphase fringes of the present invention.

【図2】本発明方法の実施形態に係るシミュレーション
において用いる干渉縞データを示す図
FIG. 2 is a diagram showing interference fringe data used in a simulation according to an embodiment of the method of the present invention.

【図3】図2の干渉縞データに対応する位相ラッピング
データを示すグラフ
FIG. 3 is a graph showing phase wrapping data corresponding to the interference fringe data of FIG. 2;

【図4】図3の位相ラッピングデータに対して位相アン
ラッピング処理を施した後のデータを示すグラフ
FIG. 4 is a graph showing data obtained by performing a phase unwrapping process on the phase wrapping data of FIG. 3;

【図5】図4のデータに対して、解析誤差を含んだ状態
で最小二乗法により傾き補正したデータを示すグラフ
5 is a graph showing data obtained by correcting the inclination of the data of FIG. 4 by a least squares method in a state including an analysis error.

【図6】モジュレーションの値の変化を示すグラフFIG. 6 is a graph showing a change in modulation value.

【図7】干渉縞計測による測定手法が適用される被検面
形状を示す斜視図
FIG. 7 is a perspective view showing a shape of a test surface to which a measurement technique based on interference fringe measurement is applied.

【符号の説明】[Explanation of symbols]

1 台座本体 2 基準面 3 吸引孔 4 シリコンウエハ 5 溝 DESCRIPTION OF SYMBOLS 1 Pedestal main body 2 Reference surface 3 Suction hole 4 Silicon wafer 5 Groove

Claims (2)

【特許請求の範囲】[Claims] 【請求項1】 縞情報が存在しない分離帯と、この分離
帯により複数個に分割され、該分離帯を挟んで互いに隣
接する領域端部の位相差が、観察される等位相縞で1/
2縞よりも小さい縞情報値を有する観察領域とからなる
被検体の縞情報を取り込み、 前記観察領域に前記分離帯を含めて前記被検体の縞情報
に基づく位相計算を行ない、 次いで、該位相計算により位相ラッピングされた位相分
布の位相アンラッピング処理を行ない、 該位相アンラッピング処理が施された位相分布におい
て、前記縞情報が存在しない分離帯に対応する部分にマ
スキング処理を施すことを特徴とする等位相縞の位相状
態解析方法。
1. A separation band having no fringe information and a plurality of divisions by the separation band, and a phase difference between end portions of regions adjacent to each other with the separation band therebetween is 1 /
Taking in fringe information of the subject consisting of an observation region having a fringe information value smaller than two fringes, performing a phase calculation based on the fringe information of the subject including the separation band in the observation region; Performing a phase unwrapping process on the phase distribution that has been phase-wrapped by calculation, and performing a masking process on a portion corresponding to a separation band where the fringe information does not exist in the phase distribution on which the phase unwrapping process has been performed. To analyze the phase state of equiphase fringes.
【請求項2】 前記マスキング処理が施された位相分布
に所定の関数補正処理を施すことを特徴とする請求項1
記載の等位相縞の位相状態解析方法。
2. The method according to claim 1, wherein a predetermined function correction process is performed on the phase distribution subjected to the masking process.
The phase state analysis method for the described equiphase fringes.
JP32863398A 1998-11-18 1998-11-18 Phase state analysis method for equiphase fringes Expired - Fee Related JP3497393B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP32863398A JP3497393B2 (en) 1998-11-18 1998-11-18 Phase state analysis method for equiphase fringes

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP32863398A JP3497393B2 (en) 1998-11-18 1998-11-18 Phase state analysis method for equiphase fringes

Publications (2)

Publication Number Publication Date
JP2000155051A true JP2000155051A (en) 2000-06-06
JP3497393B2 JP3497393B2 (en) 2004-02-16

Family

ID=18212453

Family Applications (1)

Application Number Title Priority Date Filing Date
JP32863398A Expired - Fee Related JP3497393B2 (en) 1998-11-18 1998-11-18 Phase state analysis method for equiphase fringes

Country Status (1)

Country Link
JP (1) JP3497393B2 (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002183738A (en) * 2000-12-18 2002-06-28 Fuji Photo Optical Co Ltd Analysis method for stripe image with separated area
JP2002195806A (en) * 2000-12-26 2002-07-10 Mitsutoyo Corp Method and device for unwrapping two-dimensional phase data in interferometer
JP2003042728A (en) * 2001-07-27 2003-02-13 Fuji Photo Optical Co Ltd Phase unwrapping method for fringe image analysis
JP2005274317A (en) * 2004-03-24 2005-10-06 Fujinon Corp Interferometer device for virtual contact surface measurement
JP2006079275A (en) * 2004-09-08 2006-03-23 Olympus Corp Phase connection method, its device, and phase connection program
JP2009236706A (en) * 2008-03-27 2009-10-15 Kobe Steel Ltd Shape calculator, shape calculation program, shape calculating method, and shape-measuring device
JP2009250875A (en) * 2008-04-09 2009-10-29 Canon Inc Wavefront measuring method and wavefront measuring apparatus using the wavefront measuring method
JP2010008150A (en) * 2008-06-25 2010-01-14 Kobe Steel Ltd Shape measuring apparatus
JP2011085411A (en) * 2009-10-13 2011-04-28 Nippon Steel Corp Device and method for measuring shape, and program
JP2017083448A (en) * 2015-10-30 2017-05-18 ケーエルエー−テンカー コーポレイション Method and system for reginal phase unwrapping with pattern-assisted correction
JPWO2021117198A1 (en) * 2019-12-12 2021-06-17

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4578674B2 (en) * 2000-12-18 2010-11-10 富士フイルム株式会社 Analysis method of fringe image with separated region
JP2002183738A (en) * 2000-12-18 2002-06-28 Fuji Photo Optical Co Ltd Analysis method for stripe image with separated area
JP2002195806A (en) * 2000-12-26 2002-07-10 Mitsutoyo Corp Method and device for unwrapping two-dimensional phase data in interferometer
JP2003042728A (en) * 2001-07-27 2003-02-13 Fuji Photo Optical Co Ltd Phase unwrapping method for fringe image analysis
JP2005274317A (en) * 2004-03-24 2005-10-06 Fujinon Corp Interferometer device for virtual contact surface measurement
JP2006079275A (en) * 2004-09-08 2006-03-23 Olympus Corp Phase connection method, its device, and phase connection program
JP2009236706A (en) * 2008-03-27 2009-10-15 Kobe Steel Ltd Shape calculator, shape calculation program, shape calculating method, and shape-measuring device
JP2009250875A (en) * 2008-04-09 2009-10-29 Canon Inc Wavefront measuring method and wavefront measuring apparatus using the wavefront measuring method
JP2010008150A (en) * 2008-06-25 2010-01-14 Kobe Steel Ltd Shape measuring apparatus
JP2011085411A (en) * 2009-10-13 2011-04-28 Nippon Steel Corp Device and method for measuring shape, and program
JP2017083448A (en) * 2015-10-30 2017-05-18 ケーエルエー−テンカー コーポレイション Method and system for reginal phase unwrapping with pattern-assisted correction
JPWO2021117198A1 (en) * 2019-12-12 2021-06-17
WO2021117198A1 (en) * 2019-12-12 2021-06-17 オリンパス株式会社 Sample structure measuring device and sample structure measuring method
JP7277610B2 (en) 2019-12-12 2023-05-19 株式会社エビデント SAMPLE STRUCTURE MEASURING DEVICE AND SAMPLE STRUCTURE MEASURING METHOD

Also Published As

Publication number Publication date
JP3497393B2 (en) 2004-02-16

Similar Documents

Publication Publication Date Title
JP5340539B2 (en) Method and system for surface interference analysis and related applications
TWI394930B (en) Method and apparatus for analyzing low-coherence interferometry signals for obtaining information about thin film structures
JP5290322B2 (en) Analysis of surface structure using scanning interferometry
US7468799B2 (en) Scanning interferometry for thin film thickness and surface measurements
US7494892B2 (en) Method of measuring warpage of rear surface of substrate
JP2000155051A (en) Phase state analyzing method for equiphase stripes
US8553231B2 (en) Method and apparatus for determining the height of a number of spatial positions on a sample defining a profile of a surface through white light interferometry
Gödecke et al. Model-based characterisation of complex periodic nanostructures by white-light Mueller-matrix Fourier scatterometry
US6768554B2 (en) Fringe analysis method using fourier transform
JP6758736B1 (en) Optical measurement system and optical measurement method
JP4610117B2 (en) Fourier transform fringe analysis method and apparatus
US8976337B2 (en) Method of measuring mark position and measuring apparatus
US7956630B1 (en) Real-time effective-wavelength error correction for HDVSI
JP4132365B2 (en) Analytical fringe correction method in interference fringe measurement using multiple wavelengths
CN112945152B (en) Wafer flatness detection device based on double-side grazing incidence common path self-interference technology
JP2950579B2 (en) Method and apparatus for measuring toroidal surface
JP2000337836A (en) Shape measuring method with multi-wavelength interference fringe
US7898672B1 (en) Real-time scanner-nonlinearity error correction for HDVSI
TWI290617B (en) A device and a method for detecting object profile
CN116164668A (en) Method for obtaining surface morphology of wafer
CN117422674A (en) Laser measurement and surface shape characteristic data management method
KR20010059009A (en) A method for solving of a expose process

Legal Events

Date Code Title Description
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20031113

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R360 Written notification for declining of transfer of rights

Free format text: JAPANESE INTERMEDIATE CODE: R360

R371 Transfer withdrawn

Free format text: JAPANESE INTERMEDIATE CODE: R371

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20071128

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20081128

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091128

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091128

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101128

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101128

Year of fee payment: 7

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313113

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101128

Year of fee payment: 7

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111128

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121128

Year of fee payment: 9

LAPS Cancellation because of no payment of annual fees