JP3225934B2 - Ground collision prediction method and apparatus in air traffic control, and computer-readable recording medium - Google Patents

Ground collision prediction method and apparatus in air traffic control, and computer-readable recording medium

Info

Publication number
JP3225934B2
JP3225934B2 JP33178898A JP33178898A JP3225934B2 JP 3225934 B2 JP3225934 B2 JP 3225934B2 JP 33178898 A JP33178898 A JP 33178898A JP 33178898 A JP33178898 A JP 33178898A JP 3225934 B2 JP3225934 B2 JP 3225934B2
Authority
JP
Japan
Prior art keywords
aircraft
data
velocity vector
vector
altitude
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.)
Expired - Fee Related
Application number
JP33178898A
Other languages
Japanese (ja)
Other versions
JP2000155900A (en
Inventor
善之 永野
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.)
NEC Corp
Original Assignee
NEC Corp
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 NEC Corp filed Critical NEC Corp
Priority to JP33178898A priority Critical patent/JP3225934B2/en
Publication of JP2000155900A publication Critical patent/JP2000155900A/en
Application granted granted Critical
Publication of JP3225934B2 publication Critical patent/JP3225934B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Description

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

【0001】[0001]

【発明の属する技術分野】本発明は、航空機の運航を管
理する航空管制において、3次元の地形データを用いて
航空機の地上への接近及び衝突を予測する航空管制にお
ける地上衝突予測方法および装置並びにコンピュータ読
み取り可能な記録媒体に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method and apparatus for ground collision prediction in air traffic control for predicting approach and collision of an aircraft to the ground using three-dimensional terrain data in air traffic control for managing the operation of an aircraft, and The present invention relates to a computer-readable recording medium.

【0002】[0002]

【従来の技術】航空管制システムにおいて、地上への接
触事故や墜落事故を防ぐために地上への接近や衝突を予
測し、警告を与えることは非常に重要であり、管制シス
テムの中には従来から衝突予測機能が備えられているも
のも存在する。すなわち、従来の航空管制システムにお
ける地上衝突防止は、管制対象のエリアの最低誘導高度
を規定値として持ち、航空機から送られてくる高度情報
に基づいて、航空機の高度が最低誘導高度を下回ったど
うかを判断することにより衝突予測を行っている。
2. Description of the Related Art In an air traffic control system, it is very important to predict approach and collision to the ground and to give a warning in order to prevent a collision accident or a crash on the ground. Some have a collision prediction function. In other words, ground collision prevention in conventional air traffic control systems has the minimum guidance altitude of the area to be controlled as a specified value, and based on the altitude information sent from the aircraft, whether the altitude of the aircraft falls below the minimum guidance altitude The collision prediction is performed by judging.

【0003】[0003]

【発明が解決しようとする課題】しかしながら、上記最
低誘導高度は地形の起伏を表したものではなく、そのエ
リアの最高点に安全性の余裕度を考慮した均一の高度と
なっているため、山間部の空港など、まわりが標高の高
い山に囲まれているような空港周辺の管制には適さない
という問題があった。つまり、この方法では、エリアご
とに均一な最低誘導高度を設定することになるため、山
間部などの地形が込み入っている空港では最低誘導高度
が高く設定されてしまい、効果的な衝突防止と管制の両
立が難しいという問題があった。そのため、3次元の地
形情報に基づいたより精密な衝突防止システムが求めら
れている。これに対して、国土地理院発行の数値地図デ
ータのような3次元の標高データを最低誘導高度の代わ
りに応用することが考えられる。また、従来の航空管制
においては航空機の予測到達範囲を立方体で表していた
が、より詳細な衝突予測を行うためには確率論的に考
え、航空機の位置を中心とした放射状の立体として表す
ことが考えられる。しかしながらその場合、3次元デー
タの衝突予測計算を行わなければならないため、計算量
が膨大になるという問題があった。
However, the above-mentioned minimum guidance altitude does not indicate the undulation of the terrain, but has a uniform height at the highest point of the area in consideration of the safety margin. There is a problem that it is not suitable for the control around the airport, which is surrounded by high altitude mountains, such as the airport of some parts. In other words, with this method, the minimum guidance altitude is set uniformly for each area.Therefore, the minimum guidance altitude is set high at airports that have complicated terrain, such as mountainous areas, and effective collision prevention and control There was a problem that it was difficult to balance. Therefore, a more precise collision prevention system based on three-dimensional terrain information is required. On the other hand, it is conceivable to apply three-dimensional elevation data such as digital map data issued by the Geographical Survey Institute instead of the minimum guidance altitude. Also, in conventional air traffic control, the predicted reach of the aircraft was represented by a cube, but in order to perform a more detailed collision prediction, it should be considered stochastically and represented as a radial solid centering on the position of the aircraft Can be considered. However, in that case, there is a problem that the calculation amount becomes enormous because the collision prediction calculation of the three-dimensional data must be performed.

【0004】本発明は、上記事情に鑑みてなされたもの
で、その目的とするところは、管制対象エリアにおける
3次元の標高データを保持し、3次元的に表された航空
機の予測到達範囲と各標高データとの相対関係を、航空
機の現在位置および速度ベクトルの2つの諸元から判断
することにより、実際の地形に則した地上衝突予測を行
うものであり、短い処理時間で航空機の地上への接近お
よび衝突を確実に予測することができる航空管制におけ
る地上衝突予測方法および装置並びにコンピュータ読み
取り可能な記録媒体を提供することにある。
SUMMARY OF THE INVENTION The present invention has been made in view of the above circumstances, and has as its object to hold three-dimensional elevation data in an area to be controlled, and to provide a predicted arrival range of an aircraft represented three-dimensionally. Judgment of the ground collision based on the actual terrain by judging the relative relationship with each altitude data from the two parameters of the current position and velocity vector of the aircraft. An object of the present invention is to provide a method and apparatus for predicting a ground collision in an air traffic control and a computer-readable recording medium capable of reliably predicting approach and collision of a vehicle.

【0005】[0005]

【課題を解決するための手段】本発明の請求項1は、航
空機の地上への接近、衝突予測を行う航空管制における
地上衝突予測方法において、管制対象エリアの3次元標
高データを保持し、かつ航空機の位置および3次元速度
ベクトルからこの航空機の3次元での予測到達範囲を設
定し、上記航空機の位置および上記標高データに基づい
て差分ベクトルを求め、さらにこの差分ベクトルと上記
航空機の3次元速度ベクトルに基づいて、上記標高デー
タが上記航空機の3次元での予測到達範囲内にあるかど
うかを判定し、差分ベクトルと速度ベクトルの平面成分
に基づいて標高データが航空機の予測到達範囲内にある
かどうかを判定する第1の処理と、この第1の処理によ
って予測到達範囲内と判定された場合に、上記差分ベク
トルを含む垂直面において上記速度ベクトルの投射速度
ベクトルを求め、この投射速度ベクトルと上記差分ベク
トルに基づいて上記標高データが航空機の予測到達範囲
内にあるかどうかを判定する第2の処理とからなること
を特徴とする。本発明の請求項2は、航空機の地上への
接近、衝突予測を行う航空管制における地上衝突予測装
置において、管制対象エリアの3次元標高データを保持
する標高データ記録手段と、3次元における航空機の位
置、飛行方向および速度からなる諸元データを入力する
諸元データ入力手段と、この諸元データ入力手段からの
航空機の位置に基づいて、その位置を含む地域区画とそ
の地域区画に隣接する地域区画を抽出する地域区画抽出
手段と、この地域区画抽出手段により抽出された地域区
画内の各標高データごとに、その標高データが上記航空
機の3次元での予測到達範囲内に入っているかどうかを
判定する判定手段と、この判定手段によって上記標高デ
ータが上記航空機の3次元での予測到達範囲内に入って
いると判定された場合に警告する警告手段とを備え、こ
の判定手段は、標高データが航空機の予測到達範囲の水
平面内に入っているかどうかを平面演算で判断する平面
成分判定部と、高さ方向を含めた演算を行い、上記標高
データが上記航空機の予測到達範囲内に入っているかど
うかを判断する垂直成分判定部とからなることを特徴と
する。 本発明の請求項3は、航空機の地上への接近、衝
突予測処理をコンピュータに実行させるためのプログラ
ムを記録したコンピュータ読み取り可能な記録媒体にお
いて、管制対象エリアの3次元標高データを保持し、か
つ航空機の位置および3次元速度ベクトルからこの航空
機の3次元での予測到達範囲を設定し、上記航空機の位
置および上記標高データに基づいて差分ベクトルを求
め、さらにこの差分ベクトルと上記航空機の3次元速度
ベクトルに基づいて、上記標高データが上記航空機の3
次元での予測到達範囲内にあるかどうかを判定し、差分
ベクトルと速度ベクトルの平面成分に基づいて標高デー
タが航空機の予測到達範囲内にあるかどうかを判定する
第1の処理と、この第1の処理によって予測到達範囲内
と判定された場合に、上記差分ベクトルを含む垂直面に
おいて上記速度ベクトルの投射速度ベクトルを求め、こ
の投射速度ベクトルと上記差分ベクトルに基づいて上記
標高データが航空機の予測到達範囲内にあるかどうかを
判定する第2の処理とをコンピュータに実行させること
を特徴とする。
According to a first aspect of the present invention, there is provided a ground collision prediction method in an air traffic control for predicting an approach of an aircraft to the ground and a collision, wherein three-dimensional elevation data of an area to be controlled is held; A predicted three-dimensional reach of the aircraft is set from the position of the aircraft and the three-dimensional velocity vector, a difference vector is obtained based on the position of the aircraft and the elevation data, and the difference vector and the three-dimensional velocity of the aircraft are further determined. Based on the vector, it is determined whether or not the altitude data is within the three-dimensional predicted reach of the aircraft, and a plane component of the difference vector and the velocity vector is determined.
Altitude data within aircraft's expected range based on
A first process for determining whether or not
Is determined to be within the predicted reach, the difference vector
Projection velocity of the above velocity vector on a vertical plane including the torque
Find the projection velocity vector and the difference vector
The above altitude data is based on the estimated
And a second processing for determining whether the
It is characterized by. Claim 2 of the present invention claims that
Ground collision prediction system in air traffic control that performs approach and collision prediction
3D elevation data of the area to be controlled
Altitude data recording means and aircraft position in three dimensions
Enter specifications data consisting of position, flight direction and speed
From the specification data input means and from the specification data input means.
Based on the location of the aircraft, the area parcel containing that location and its location
Area extraction to extract area areas adjacent to the area area
Means and the regional ward extracted by this regional division extracting means
For each altitude data in the picture, the altitude data is
Whether the aircraft is within the 3D prediction range
Determining means, and the altitude data is determined by the determining means.
Data is within the three-dimensional predicted reach of the aircraft
Warning means for warning when it is determined that
Means that the altitude data is
Plane to determine whether it is inside the plane by plane operation
Performs a calculation including the height direction with the component determination unit, and calculates
Whether the data is within the expected range of the aircraft
And a vertical component determination unit that determines whether
I do. Claim 3 of the present invention is directed to approaching the aircraft to the ground,
A program for causing a computer to execute collision prediction processing
On a computer-readable recording medium on which
And hold the 3D elevation data of the area to be controlled
From this aircraft position and three-dimensional velocity vector
Set the three-dimensional predicted reach of the aircraft and
Difference vector based on the
And the difference vector and the three-dimensional velocity of the aircraft
Based on the vector, the elevation data is
Judge whether it is within the predicted reach in the dimension,
Elevation day based on vector and velocity vector plane components
Determine if the aircraft is within the expected range of the aircraft
The first processing and the first processing
Is determined, the vertical plane including the difference vector
The projection velocity vector of the above velocity vector.
Based on the projection velocity vector and the difference vector
Whether the altitude data is within the expected range of the aircraft
Causing the computer to execute the second processing to be determined
It is characterized by.

【0006】[0006]

【発明の実施の形態】以下、図面を参照して本発明の一
実施形態について説明する。図1は本発明の一実施形態
を示すブロック図である。この図1において、本実施形
態の地上衝突予測装置は、標高データを記憶した標高デ
ータ記憶部1と、航空機の諸元データを入力する諸元デ
ータ入力部2と、評価対象の地域区画を抽出する地域区
画抽出部3と、水平面成分で標高データが航空機の予測
到達範囲にあるかどうかを判断する水平成分判定部4
と、高度方向で接近を判定する垂直成分判定部5と、接
近を警告する警告表示出力部6とを含む。標高データ記
憶部1は、ある緯・経度における標高データの組み合わ
せを、地域区画ごとに格納している。図2は地域区画の
例を表している。地域区画はある特定の大きさの四角い
エリアであり、4辺がそれぞれ隣接する地域区画の一辺
に接している。また図4は、1つの地域区画の内部を表
している。図4を参照すると、各地域区画は縦横それぞ
れ数個から数十のメッシュに等分され、メッシュの交点
ごとに標高データを持つ。諸元データ入力部2は、管制
対象の航空機の現在位置、飛行方向および速度からなる
諸元データを、航空機やレーダ装置から入力する。地域
区画抽出部3は、図2に示すとおり、諸元データにおい
て入力される航空機の現在位置から、その現在位置を含
む地域区画と、その区画に隣接する地域区画を抽出す
る。水平成分判定部4は、抽出された地域区画内の各標
高データごとに、その標高データが管制対象となる航空
機の予測到達範囲の水平面内に入っているかどうかを2
次元の平面演算で判断し、絞り込みを行う。垂直成分判
定部5は、水平成分判定部4の結果がある閾値を下回っ
た場合に高さ方向を含めた3次元の演算を行い、標高デ
ータが予測到達範囲に入っているかどうかを判断する。
警告表示出力部6は、垂直成分判定部5においてある閾
値を下回った場合に、標高データが航空機の予測到達範
囲に入っていることを警告する。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS An embodiment of the present invention will be described below with reference to the drawings. FIG. 1 is a block diagram showing one embodiment of the present invention. In FIG. 1, the ground collision prediction apparatus according to the present embodiment extracts an elevation data storage unit 1 storing elevation data, an specification data input unit 2 for inputting aircraft specification data, and an area block to be evaluated. And a horizontal component judging unit 4 for judging whether or not the altitude data is within the predicted reach of the aircraft in the horizontal plane component
And a vertical component determining unit 5 for determining the approach in the altitude direction, and a warning display output unit 6 for warning the approach. The elevation data storage unit 1 stores a combination of elevation data at a certain latitude / longitude for each area block. FIG. 2 shows an example of a regional block. The regional section is a rectangular area of a specific size, and four sides are in contact with one side of the adjacent regional section. FIG. 4 shows the inside of one area block. Referring to FIG. 4, each regional section is equally divided into several to several tens of meshes in each of the vertical and horizontal directions, and has elevation data at each intersection of the meshes. The specification data input unit 2 inputs specification data including the current position, flight direction, and speed of the aircraft to be controlled from the aircraft and the radar device. As shown in FIG. 2, the regional section extraction unit 3 extracts, from the current position of the aircraft input in the specification data, a regional section including the current position and a regional section adjacent to the current section. The horizontal component determination unit 4 determines for each of the elevation data in the extracted area block whether or not the elevation data is within the horizontal plane of the predicted reach of the aircraft to be controlled.
Judgment is performed by dimensional plane calculation and narrowing down is performed. The vertical component determination unit 5 performs a three-dimensional calculation including the height direction when the result of the horizontal component determination unit 4 falls below a certain threshold value, and determines whether the altitude data falls within the predicted reach.
When the vertical component determination unit 5 falls below a certain threshold, the warning display output unit 6 warns that the altitude data is within the predicted reach of the aircraft.

【0007】次に、本発明の航空管制における地上接近
・衝突予測方法について、図を参照して説明する。図6
において、諸元データ入力部2から入力される管制対象
の航空機の現在位置に関する情報から、まず評価対象の
地域区画を抽出する(ステップA1)。図2を参照する
と、航空機の位置を含む地域区画と、その地域区画を取
り囲む8区画が評価対象の地域区画として抽出される。
Next, a ground approach / collision prediction method in air traffic control according to the present invention will be described with reference to the drawings. FIG.
In step A1, an area block to be evaluated is first extracted from the information on the current position of the aircraft to be controlled input from the specification data input unit 2 (step A1). Referring to FIG. 2, an area section including the position of the aircraft and eight sections surrounding the area section are extracted as area sections to be evaluated.

【0008】次に、抽出された地域区画の中の標高デー
タを1つ抽出する(ステップA2)。抽出された標高デ
ータに対して、その標高点が航空機の予測到達範囲内に
入るかどうかを判断する。図3を参照すると、航空機の
予測到達範囲32は、航空機の位置31からその航空機
の速度ベクトルを中心に、平面、垂直の各方向ごとに角
度の余裕度を考慮した放射状の立体となる。本発明で
は、水平面の形状は航空機の位置を中心とした扇型とし
(図4の42)、高さ方向は計算量を考慮して、最遠点
を垂直とする三角形(図5の52)として予測到達範囲
を定義する。また予測到達範囲の航空機の位置からの大
きさは、速度ベクトルの大きさに比例した値となる。
Next, one altitude data is extracted from the extracted area block (step A2). For the extracted altitude data, it is determined whether or not the altitude point falls within the predicted reach of the aircraft. Referring to FIG. 3, the predicted arrival range 32 of the aircraft is a radial solid from the position 31 of the aircraft and centered on the velocity vector of the aircraft, taking into account the degree of margin in each of the plane and vertical directions. In the present invention, the shape of the horizontal plane is a sector shape centering on the position of the aircraft (42 in FIG. 4), and the height direction is a triangle with the farthest point being vertical in consideration of the amount of calculation (52 in FIG. 5). Is defined as the predicted reach. The size of the predicted reach from the position of the aircraft is a value proportional to the size of the speed vector.

【0009】まず、抽出された標高データごとに、その
水平成分が予測到達範囲内にあるかどうかを判断する。
図4を参照すると、標高データHij(hx,hy,h
z)の平面成分は(hx,hy)である。まず、航空機
の位置T(tx,ty,tz)からみた標高データHi
jの差分ベクトルDを求める(ステップA3)。差分ベ
クトルDは以下の通りとなる。 D(dx,dy,dz)=(hx−tx,hy−ty,
hz−tz) 次に、差分ベクトルDの水平成分(dx,dy)の大き
さが、速度ベクトルVの水平面成分の大きさの定数倍よ
りも大きいかどうかを判断する(ステップA4)。ここ
で定数(第1の規定値)をαとすると、図7のb1式に
表す条件式となる。b1の条件式を満たせば、標高デー
タは予測到達範囲の円の外にあるため、それ以降の計算
の対象外とする。なお、上記第1の規定値αは、航空機
がたとえば約20分間で到達する距離に応じて設定され
ている。
First, it is determined whether or not the horizontal component of each of the extracted elevation data is within the predicted reach.
Referring to FIG. 4, the elevation data Hij (hx, hy, h
The plane component of z) is (hx, hy). First, the altitude data Hi from the position T (tx, ty, tz) of the aircraft
The difference vector D of j is obtained (step A3). The difference vector D is as follows. D (dx, dy, dz) = (hx-tx, hy-ty,
hz-tz) Next, it is determined whether or not the magnitude of the horizontal component (dx, dy) of the difference vector D is larger than a constant multiple of the magnitude of the horizontal component of the velocity vector V (step A4). Here, assuming that a constant (first specified value) is α, a conditional expression represented by the expression b1 in FIG. 7 is obtained. If the conditional expression of b1 is satisfied, the altitude data is outside the circle of the predicted reach, so that it is excluded from subsequent calculations. The first specified value α is set in accordance with the distance that the aircraft reaches in about 20 minutes, for example.

【0010】もしb1の条件式を満たさない場合には、
航空機の位置を中心として速度ベクトルの定数倍の長さ
を半径とする円内(航空機の予測到達範囲内)にあるこ
とになる。そこで、次に差分ベクトルDと速度ベクトル
Vとの平面成分の内積を計算し(ステップA5)、図4
の42に示される扇型の予測到達範囲の範囲内にあるか
どうかを判定する(ステップA6)。差分ベクトルDと
速度ベクトルVとの内積と、差分ベクトルDと速度ベク
トルVとのなす角θaとの間には図7のb2式で表され
る関係がある。このb2式において、cosθaが、第
2の規定値βを代入したcosβよりも小さい場合に
は、角度θaが第2の規定値βよりも大きいことにな
る。したがって、図4の42に示された扇型の範囲外に
あることになるため、それ以降の計算の対象外とする。
ここで第2の規定値βの値は、扇型の中心の角度が45
度の場合、図7のb3のようになる。
If the condition of b1 is not satisfied,
It is within a circle (within the predicted reach of the aircraft) having a radius equal to a constant multiple of the velocity vector around the position of the aircraft. Accordingly, next, the inner product of the plane component of the difference vector D and the velocity vector V is calculated (step A5), and FIG.
It is determined whether it is within the range of the fan-shaped predicted reach indicated by No. 42 (step A6). The inner product of the difference vector D and the speed vector V and the angle θa formed by the difference vector D and the speed vector V have a relationship represented by the equation b2 in FIG. In the equation b2, when cos θa is smaller than cos β obtained by substituting the second specified value β, the angle θa is larger than the second specified value β. Therefore, since it is out of the fan-shaped range shown by 42 in FIG. 4, it is excluded from the subsequent calculations.
Here, the value of the second specified value β is such that the angle of the center of the sector is 45 degrees.
In the case of the degree, it becomes like b3 of FIG.

【0011】一方、図7のb2式の値cosθaがco
sβよりも大きい場合には、標高データの水平面成分
(hx、hy)は扇型で表される予測到達範囲の水平面
内に存在することが判明する。そこで、次に標高データ
の標高値hzが予測到達範囲内にあるかどうかを判断す
るため、航空機の位置Tと標高データHijを結ぶ線を
垂直に切った断面に対して、同様にベクトルの内積演算
を行う。まず、速度ベクトルVと同じ大きさで、ベクト
ルの方向が差分ベクトルDを含む垂直面に一致する(垂
直面に投射した)投射速度ベクトルVbを求める(ステ
ップA7)。投射速度ベクトルVbの値は図7のb4式
のようになる。次に差分ベクトルDと投射速度ベクトル
Vbとの3次元の内積を計算し(ステップA8)、図5
の予測到達範囲52の範囲内にあるかどうかを判定する
(ステップA9)。差分ベクトルDと投射速度ベクトル
Vbとの内積と、差分ベクトルDと投射速度ベクトルV
bとのなす角θbとの間には、図7のb5式で表される
関係がある。このb5式において、cosθbが第3の
規定値γを代入したcosγよりも小さい場合には、航
空機の予測到達範囲の範囲外にあることが判明し、次の
標高データの判別に移る(ステップA11)。ここで、
第3の規定値γはたとえば15〜20度の範囲の値に設
定される。
On the other hand, the value cosθa of the equation b2 in FIG.
When it is larger than sβ, it is found that the horizontal component (hx, hy) of the elevation data exists in the horizontal plane of the predicted reach represented by the fan shape. Then, in order to determine whether or not the altitude value hz of the altitude data is within the predicted reach, the inner product of the vector is similarly calculated for a cross section obtained by cutting a line connecting the aircraft position T and the altitude data Hij vertically. Perform the operation. First, a projection velocity vector Vb having the same size as the velocity vector V and whose vector direction matches the vertical plane including the difference vector D (projected on the vertical plane) is obtained (step A7). The value of the projection velocity vector Vb is as shown by the b4 equation in FIG. Next, a three-dimensional inner product of the difference vector D and the projection velocity vector Vb is calculated (step A8), and FIG.
(Step A9). The inner product of the difference vector D and the projection speed vector Vb, and the difference vector D and the projection speed vector V
There is a relationship expressed by the equation b5 in FIG. 7 between the angle θb and the angle θb. In the equation b5, when cos θb is smaller than cos γ into which the third specified value γ is substituted, it is determined that the aircraft is out of the predicted reach of the aircraft, and the process proceeds to the next elevation data determination (step A11). ). here,
The third specified value γ is set, for example, to a value in the range of 15 to 20 degrees.

【0012】一方、b5式の値cosθbがcosγよ
りも大きい場合には、標高データは垂直成分すなわち高
さ方向でも予測到達範囲内にあることが判明し、警告表
示を行う(ステップA10)。
On the other hand, if the value cosθb of the expression b5 is larger than cosγ, it is determined that the elevation data is within the predicted reach even in the vertical component, that is, in the height direction, and a warning is displayed (step A10).

【0013】このように、本実施形態においては、航空
機の予測到達範囲と各標高データの相対関係を判断する
ため、まず水平面の2次元演算により範囲が重なるかど
うかを判断し、次にその結果がある閾値を下回った場合
のみ高度方向の演算を行うことにより、計算の処理時間
を大幅に低減することができる。なお、本実施形態は、
航空管制システムにおける地上衝突予測方法について説
明したが、航空機搭載用の地上衝突予測システムにも応
用することができる。また、本実施形態において、図1
に示す制御部の動作あるいは図6に示すフローチャート
を実現するためのプログラムをコンピュータ読み取り可
能な記録媒体に記録して、この記録媒体に記録されたプ
ログラムをコンピュータに読み込ませ、実行することに
より、航空機の地上への接近および衝突の予測を行って
もよい。
As described above, in this embodiment, in order to determine the relative relationship between the predicted reach of an aircraft and each elevation data, it is first determined whether or not the ranges overlap by a two-dimensional operation on a horizontal plane, and then the result is determined. By performing the calculation in the altitude direction only when the value falls below a certain threshold, the processing time of the calculation can be significantly reduced. In this embodiment,
Although the ground collision prediction method in the air traffic control system has been described, the invention can also be applied to a ground collision prediction system mounted on an aircraft. In the present embodiment, FIG.
The program for realizing the operation of the control unit shown in FIG. 6 or the flowchart shown in FIG. 6 is recorded on a computer-readable recording medium, and the program recorded on this recording medium is read and executed by a computer, thereby May be predicted for approaching and colliding with the ground.

【0014】[0014]

【発明の効果】本発明の請求項1は、航空機の地上への
接近、衝突予測を行う航空管制における地上衝突予測方
法において、管制対象エリアの3次元標高データを保持
し、かつ航空機の位置および3次元速度ベクトルからこ
の航空機の3次元での予測到達範囲を設定し、上記航空
機の位置および上記標高データに基づいて差分ベクトル
を求め、さらにこの差分ベクトルと上記航空機の3次元
速度ベクトルに基づいて、上記標高データが上記航空機
の3次元での予測到達範囲内にあるかどうかを判定する
ことを特徴とするものであるから、地形状況に合わせた
詳細な地上への接近、衝突予測を行うことができる。こ
の結果、より効率的かつ安全な航空管制を行うことがで
きる。その理由は、管制対象エリアにおける3次元の標
高データを保持し、3次元的に表された航空機の予測到
達範囲と各標高データとの包含関係を、航空機の現在位
置および速度ベクトルの2つの諸元から判断するためで
ある。更に、差分ベクトルと速度ベクトルの平面成分に
基づいて標高データが航空機の予測到達範囲内にあるか
どうかを判定する第1の処理と、この第1の処理によっ
て予測到達範囲内と判定された場合に、上記差分ベクト
ルを含む垂直面において上記速度ベクトルの投射速度ベ
クトルを求め、この投射速度ベクトルと上記差分ベクト
ルに基づいて上記標高データが航空機の予測到達範囲内
にあるかどうかを判定する第2の処理とからなることを
特徴とするものであるから、航空機の予測到達範囲と標
高データの領域の重なり判定処理を簡略化することがで
きる。その結果、計算処理時間を短縮することができ
る。その理由は、航空機の予測到達範囲と各標高データ
の包含関係を判断するために、まず水平面の2次元演算
により領域が重なるかどうかを判断し、次にその結果が
ある閾値を下回った場合のみ、高度方向の演算を行うた
めである。本発明の請求項2は、航空機の地上への接
近、衝突予測を行う航空管制における地上衝突予測装置
において、管制対象エリアの3次元標高データを保持す
る標高データ記録手段と、3次元における航空機の位
置、飛行方向および速度からなる諸元データを入力する
諸元データ入力手段と、この諸元データ入力手段からの
航空機の位置に基づいて、その位置を含む地域区画とそ
の地域区画に隣接する地 域区画を抽出する地域区画抽出
手段と、この地域区画抽出手段により抽出された地域区
画内の各標高データごとに、その標高データが上記航空
機の3次元での予測到達範囲内に入っているかどうかを
判定する判定手段と、この判定手段によって上記標高デ
ータが上記航空機の3次元での予測到達範囲内に入って
いると判定された場合に警告する警告手段とを備えたこ
とを特徴とするものであるから、判定手段によって、地
域区画抽出手段により抽出された地域区画内の各標高デ
ータごとに、その標高データが上記航空機の予測到達範
囲内に入っているかどうかを判定することにより、実際
の地形に則した地上衝突予測を行うことができ、航空機
の地上への接近および衝突を確実に予測することができ
る。更に、この判定手段は、標高データが航空機の予測
到達範囲の水平面内に入っているかどうかを平面演算で
判断する平面成分判定部と、高さ方向を含めた演算を行
い、上記標高データが上記航空機の予測到達範囲内に入
っているかどうかを判断する垂直成分判定部とからなる
ことを特徴とするものであるから、まず平面成分判定部
によって判定を行い、次いで、この判定によって航空機
の予測到達範囲内とされた標高データに対してのみ、垂
直成分判定部による判定を行うことにより、無駄な判定
処理を省くことができて処理時間を短縮することができ
る。本発明の請求項3は、航空機の地上への接近、衝突
予測処理をコンピュータに実行させるためのプログラム
を記録したコンピュータ読み取り可能な記録媒体におい
て、管制対象エリアの3次元標高データを保持し、かつ
航空機の位置および3次元速度ベクトルからこの航空機
の3次元での予測到達範囲を設定し、上記航空機の位置
および上記標高データに基づいて差分ベクトルを求め、
さらにこの差分ベクトルと上記航空機の3次元速度ベク
トルに基づいて、上記標高データが上記航空機の3次元
での予測到達範囲内にあるかどうかを判定する処理をコ
ンピュータに実行させることを特徴とするものであるか
ら、コンピュータによって航空機の地上への接近および
衝突を迅速にかつ確実に予測することができる。更に、
差分ベクトルと速度ベクトルの平面成分に基づいて標高
データが航空機の予測到達範囲内にあるかどうかを判定
する第1の処理と、この第1の処理によって予測到達範
囲内と判定された場合に、上記差分ベクトルを含む垂直
面において上記速度ベクトルの投射速度ベクトルを求
め、この投射速度ベクトルと上記差分ベクト ルに基づい
て上記標高データが航空機の予測到達範囲内にあるかど
うかを判定する第2の処理とをコンピュータに実行させ
ることを特徴とするものであるから、コンピュータによ
って第1の処理および第2の処理に分けて判定を行うこ
とにより、判定処理の簡略化を図ることができる
According to a first aspect of the present invention, there is provided a ground collision prediction method in an air traffic control for predicting an approach of an aircraft to the ground and a collision, wherein three-dimensional elevation data of an area to be controlled is held, and the position and the position of the aircraft are controlled. A three-dimensional predicted reach of the aircraft is set from the three-dimensional velocity vector, a difference vector is determined based on the position of the aircraft and the elevation data, and further based on the difference vector and the three-dimensional velocity vector of the aircraft. And determining whether or not the altitude data is within the three-dimensional predicted reach of the aircraft. Therefore, it is necessary to perform detailed approach to the ground and collision prediction in accordance with the terrain condition. Can be. As a result, more efficient and safe air traffic control can be performed. The reason for this is that three-dimensional elevation data in the control target area is retained, and the inclusion relationship between the predicted reach of the aircraft represented three-dimensionally and each elevation data is determined by the current position and velocity vector of the aircraft. This is to judge from the beginning. Furthermore, the difference vector and the velocity vector
The elevation data is within the expected range of the aircraft based on
A first process for determining whether or not
The difference vector
Of the above velocity vector on the vertical plane containing the
Vector and the difference vector
The above altitude data is within the estimated range of the aircraft based on the
And a second process for determining whether
Because of the characteristics of the aircraft,
It is possible to simplify the process of determining the overlap of high data areas.
Wear. As a result, the calculation processing time can be reduced
You. The reason is that the estimated range of the aircraft and each elevation data
First, to determine the inclusion relationship of
To determine if the areas overlap, and then
Only when the value falls below a certain threshold, the altitude direction
It is. Claim 2 of the present invention relates to the connection of the aircraft to the ground.
Ground collision prediction system in air traffic control that performs near and collision prediction
Holds 3D elevation data of the control target area
Altitude data recording means and aircraft position in three dimensions
Enter specifications data consisting of position, flight direction and speed
From the specification data input means and from the specification data input means.
Based on the location of the aircraft, the area parcel containing that location and its location
Regional compartment extraction to extract the regional compartment adjacent to the regional section of the
Means and the regional ward extracted by this regional division extracting means
For each altitude data in the picture, the altitude data is
Whether the aircraft is within the 3D prediction range
Determining means, and the altitude data is determined by the determining means.
Data is within the three-dimensional predicted reach of the aircraft
Warning means for warning when it is determined that
And the determination means determines
Each elevation data within the area plot extracted by the area plot extraction means
For each aircraft, its altitude data is
By determining whether you are in the enclosure,
Ground collision prediction based on the terrain
Can reliably predict approach and collisions on the ground
You. Further, the determination means determines whether the altitude data is
Whether or not it is within the horizontal plane of the reach
Performs a calculation including the plane component determination unit and height direction
The altitude data falls within the predicted reach of the aircraft.
And a vertical component determination unit that determines whether
First, the plane component determination unit
And then the aircraft
Only for elevation data within the expected
Useless judgment by performing judgment by the direct component judgment unit
Processing can be omitted and processing time can be reduced
You. According to the third aspect of the present invention, an aircraft approaches or collides with the ground.
A program for causing a computer to execute a prediction process
On a computer-readable recording medium on which
To hold the three-dimensional elevation data of the control target area, and
From the position of the aircraft and the three-dimensional velocity vector, this aircraft
Set the 3D predicted reach, and position the aircraft
And a difference vector is obtained based on the elevation data,
Furthermore, the difference vector and the three-dimensional velocity vector of the aircraft
The altitude data is based on the aircraft
The processing to determine whether or not the
Is it characterized by having a computer execute it?
Access to the ground by computer
Collisions can be quickly and reliably predicted. Furthermore,
Elevation based on the plane component of the difference vector and velocity vector
Determine if data is within expected range of aircraft
And a predicted reachable range by the first process.
If it is determined to be within the box, the vertical
Find the projected velocity vector of the above velocity vector on the surface
Because, based on the projection velocity vector and the difference vector
Whether the altitude data is within the expected range of the aircraft
The computer performs a second process of determining whether
Computer.
Thus, the determination can be made separately in the first process and the second process.
Thus, the determination process can be simplified .

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

【図1】 本発明の一実施形態を示すブロック図であ
る。
FIG. 1 is a block diagram showing one embodiment of the present invention.

【図2】 地域区画の説明図である。FIG. 2 is an explanatory diagram of an area block.

【図3】 航空機の予測到達範囲と位置および速度ベク
トルと標高データを説明する説明図である。
FIG. 3 is an explanatory diagram illustrating a predicted arrival range, a position, a velocity vector, and altitude data of an aircraft.

【図4】 地域区画(標高データの管理単位)内の状態
を説明する説明図である。
FIG. 4 is an explanatory diagram illustrating a state in an area block (a unit of management of elevation data).

【図5】 垂直面内における航空機の予測到達範囲と位
置および速度ベクトルと標高データを説明する説明図で
ある。
FIG. 5 is an explanatory diagram illustrating a predicted arrival range, a position, a velocity vector, and altitude data of an aircraft in a vertical plane.

【図6】 本発明の一実施形態の動作を説明する流れ図
(フローチャート)である。
FIG. 6 is a flowchart (flow chart) for explaining the operation of the embodiment of the present invention.

【図7】 関係式を列記した説明図である。FIG. 7 is an explanatory diagram listing relational expressions.

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

1 標高データ記憶部 2 諸元データ入力部 3 地域区画抽出部 4 水平成分判定部 5 垂直成分判定部 6 警告表示出力部 31、41、51 航空機の位置 32、42、52 航空機の予測到達範囲 DESCRIPTION OF SYMBOLS 1 Elevation data storage part 2 Specification data input part 3 Area division extraction part 4 Horizontal component determination part 5 Vertical component determination part 6 Warning display output part 31,41,51 Aircraft position 32,42,52 Predicted arrival range of aircraft

───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.7,DB名) G08G 5/00 - 5/06 G01S 13/91 - 13/93 G05D 1/10 JICSTファイル(JOIS)────────────────────────────────────────────────── ─── Continued on the front page (58) Field surveyed (Int. Cl. 7 , DB name) G08G 5/00-5/06 G01S 13/91-13/93 G05D 1/10 JICST file (JOIS)

Claims (3)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 航空機の地上への接近、衝突予測を行う
航空管制における地上衝突予測方法において、 管制対象エリアの3次元標高データを保持し、かつ航空
機の位置および3次元速度ベクトルからこの航空機の3
次元での予測到達範囲を設定し、上記航空機の位置およ
び上記標高データに基づいて差分ベクトルを求め、さら
にこの差分ベクトルと上記航空機の3次元速度ベクトル
に基づいて、上記標高データが上記航空機の3次元での
予測到達範囲内にあるかどうかを判定する航空管制にお
ける地上衝突予測方法において、差分ベクトルと速度ベクトルの平面成分に基づいて標高
データが航空機の予測到達範囲内にあるかどうかを判定
する第1の処理と、この第1の処理によって予測到達範
囲内と判定された場合に、上記差分ベクトルを含む垂直
面において上記速度ベクトルの投射速度ベクトルを求
め、この投射速度ベクトルと上記差分ベクトルに基づい
て上記標高データが航空機の予測到達範囲内にあるかど
うかを判定する第2の処理とからなることを特徴とする
航空管制における地上衝突予測方法。
A ground collision prediction method in an air traffic control for predicting an approach of an aircraft to the ground and a collision, wherein three-dimensional elevation data of an area to be controlled is retained, and the aircraft's position and three-dimensional velocity vector are used to calculate the 3
A predicted reach in dimension is set, and a difference vector is obtained based on the position of the aircraft and the altitude data. Further, based on the difference vector and the three-dimensional velocity vector of the aircraft, the altitude data is set to 3 In a ground collision prediction method in air traffic control that determines whether the vehicle is within a predicted reach in three dimensions, an elevation based on a plane component of a difference vector and a velocity vector is used.
Determine if data is within expected range of aircraft
And a predicted reachable range by the first process.
If it is determined to be within the box, the vertical
Find the projected velocity vector of the above velocity vector on the surface
Based on this projection velocity vector and the above difference vector
Whether the altitude data is within the expected range of the aircraft
Ground collision prediction method in air traffic control, characterized by comprising a second process for determining emergence.
【請求項2】 航空機の地上への接近、衝突予測を行う
航空管制における地上衝突予測装置において、 管制対象エリアの3次元標高データを保持する標高デー
タ記録手段と、 3次元における航空機の位置、飛行方向および速度から
なる諸元データを入力する諸元データ入力手段と、 この諸元データ入力手段からの航空機の位置に基づい
て、その位置を含む地域区画とその地域区画に隣接する
地域区画を抽出する地域区画抽出手段と、 この地域区画抽出手段により抽出された地域区画内の各
標高データごとに、その標高データが上記航空機の3次
元での予測到達範囲内に入っているかどうかを判定する
判定手段と、 この判定手段によって上記標高データが上記航空機の3
次元での予測到達範囲内に入っていると判定された場合
に警告する警告手段とを備え、 判定手段は、標高データが航空機の予測到達範囲の水平
面内に入っているかどうかを平面演算で判断する平面成
分判定部と、高さ方向を含めた演算を行い、上記 標高デ
ータが上記航空機の予測到達範囲内に入っているかどう
かを判断する垂直成分判定部とからなることを特徴とす
る航空管制における地上衝突予測装置。
2. Predicting approach and collision of an aircraft to the ground
In a ground collision prediction device in air traffic control, elevation data that holds three-dimensional elevation data of the control target area
Data recording means and the position, flight direction and speed of the aircraft in three dimensions
Data input means for inputting specification data, and a position of the aircraft from the specification data input means.
Area adjacent to the area containing the location
A regional division extracting means for extracting a regional division, and each of the regional divisions extracted by the regional division extracting means.
For each altitude data, the altitude data is the third order of the aircraft.
Judge whether it is within the original expected range
Determining means for determining whether the altitude data of the aircraft is 3
If it is determined that the value is within the predicted reach of the dimension
Warning means for warning the aircraft that the altitude data indicates the horizontal
A plane component that determines whether or not it is
And minutes judging section, the operation, including the height direction is performed, the elevation de
Data is within the expected range of the aircraft.
And a vertical component determining unit for determining whether
Ground collision prediction device in air traffic control.
【請求項3】 航空機の地上への接近、衝突予測処理を
コンピュータに実行させるためのプログラムを記録した
コンピュータ読み取り可能な記録媒体において、管制対
象エリアの3次元標高データを保持し、かつ航空機の位
置および3次元速度ベクトルからこの航空機の3次元で
の予測到達範囲を設定し、上記航空機の位置および上記
標高データに基づいて差分ベクトルを求め、さらにこの
差分ベクトルと上記航空機の3次元速度ベクトルに基づ
いて、上記標高データが上記航空機の3次元での予測到
達範囲内にあるかどうかを判定し、 差分ベクトルと速度ベクトルの平面成分に基づいて標高
データが航空機の予測到達範囲内にあるかどうかを判定
する第1の処理と、この第1の処理によって予測到達範
囲内と判定された場合に、上記差分ベクトルを含む垂直
面において上記速度ベクトルの投射速度ベクトルを求
め、この投射速度ベクトルと上記差分ベクトルに基づい
て上記標高データが航空機の予測到達範囲内にあるかど
うかを判定する第2の処理とをコンピュータに実行させ
るためのプログラムを記録したコンピュータ読み取り可
能な記録媒体。
3. An aircraft approaching the ground and a collision prediction process.
Recorded program to be executed by computer
On computer readable storage media,
Holds three-dimensional elevation data of the elephant area and ranks the aircraft
From the position and the three-dimensional velocity vector
Set the predicted reach of the aircraft, the location of the aircraft and the
A difference vector is obtained based on the elevation data, and
Based on the difference vector and the three-dimensional velocity vector of the aircraft
And the altitude data is a three-dimensional prediction of the aircraft.
Judge whether it is within the reach, and calculate the altitude based on the plane component of the difference vector and the velocity vector.
Determine if data is within expected range of aircraft
And a predicted reachable range by the first process.
If it is determined to be within the box, the vertical
Find the projected velocity vector of the above velocity vector on the surface
Based on this projection velocity vector and the above difference vector
Whether the altitude data is within the expected range of the aircraft
The computer performs a second process of determining whether
Computer-readable recording program for
Recording media.
JP33178898A 1998-11-20 1998-11-20 Ground collision prediction method and apparatus in air traffic control, and computer-readable recording medium Expired - Fee Related JP3225934B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP33178898A JP3225934B2 (en) 1998-11-20 1998-11-20 Ground collision prediction method and apparatus in air traffic control, and computer-readable recording medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP33178898A JP3225934B2 (en) 1998-11-20 1998-11-20 Ground collision prediction method and apparatus in air traffic control, and computer-readable recording medium

Publications (2)

Publication Number Publication Date
JP2000155900A JP2000155900A (en) 2000-06-06
JP3225934B2 true JP3225934B2 (en) 2001-11-05

Family

ID=18247656

Family Applications (1)

Application Number Title Priority Date Filing Date
JP33178898A Expired - Fee Related JP3225934B2 (en) 1998-11-20 1998-11-20 Ground collision prediction method and apparatus in air traffic control, and computer-readable recording medium

Country Status (1)

Country Link
JP (1) JP3225934B2 (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3579685B2 (en) 2001-10-24 2004-10-20 独立行政法人電子航法研究所 Aircraft position display method in display device for air traffic control
EP1911675A4 (en) * 2005-06-08 2010-12-08 Toshio Tsuyuki Navigation system
JP2010033349A (en) * 2008-07-29 2010-02-12 Toshiba Corp Contact point computing method and program
US8731810B2 (en) * 2010-12-10 2014-05-20 The Boeing Company Aircraft path conformance monitoring
RU2645734C1 (en) * 2016-10-04 2018-02-28 Акционерное общество "Научно-исследовательский институт Приборостроения имени В.В. Тихомирова" Method of radar-objective detection of dangerous obstacles at low-altitude mission of the flighting apparatus
CN114202967B (en) * 2021-12-22 2023-02-24 中国商用飞机有限责任公司 TCAS avoidance method and system suitable for crowded airspace and display and alarm mechanism thereof

Also Published As

Publication number Publication date
JP2000155900A (en) 2000-06-06

Similar Documents

Publication Publication Date Title
JP6800899B2 (en) Risk-based driver assistance for approaching intersections with limited visibility
EP1203201B1 (en) Hazard detection for flight plans and the like
JP7346707B2 (en) Polyline contour representation of autonomous vehicles
Rummelhard et al. Probabilistic grid-based collision risk prediction for driving application
KR101621649B1 (en) Method for Generating Position Distribution Information of Moving Obstacle, Method and System for Controlling Vehicle Using the Position Distribution Information
US7937191B2 (en) Termination secured route planning
EP3659002B1 (en) Vehicle interface for autonomous vehicle
US11899462B2 (en) Motion-plan validator for autonomous vehicle
US20120123671A1 (en) Method and system for detecting potential collisions by an aircraft
CN112419792B (en) Aircraft flight conflict detection method, system, device and medium
JP3019071B2 (en) Intrusion / collision prediction apparatus and method, and recording medium recording intrusion / collision prediction program
US11608058B2 (en) Method of and system for predicting future event in self driving car (SDC)
EP3471078A1 (en) Method and system for rendering and displaying a perspective view of aircraft taxi operation
JP3225934B2 (en) Ground collision prediction method and apparatus in air traffic control, and computer-readable recording medium
JP2022129175A (en) Vehicle evaluation method and vehicle evaluation device
US10613815B2 (en) Apparatus for positioning an interactive-display device within an interior of an autonomous vehicle
Abramson et al. The Generic Resolution Advisor and Conflict Evaluator (GRACE) for Unmanned Aircraft Detect-And-Avoid Systems
Feng et al. Collision detection algorithm of successive frames for automated air traffic control
CN115440094B (en) Obstacle detection method, device and storage medium for helicopter ground proximity warning
JP7353239B2 (en) Motion prediction device and motion prediction program
Cruck et al. Sense and avoid system for a MALE UAV
Akgunduz et al. Smart collision information processing sensors for fast moving objects
CN117635721A (en) Target positioning method, related system and storage medium
Akgunduz Collision detection in fast moving objects
Radzi Collision Avoidance Algorithm Design for UAV Base on Parametric Theorem and Circle Overlapping Method

Legal Events

Date Code Title Description
A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20000912

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20010731

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

Free format text: PAYMENT UNTIL: 20070831

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20080831

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20080831

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20090831

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20090831

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20100831

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20110831

Year of fee payment: 10

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

Free format text: PAYMENT UNTIL: 20110831

Year of fee payment: 10

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

Free format text: PAYMENT UNTIL: 20120831

Year of fee payment: 11

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

Free format text: PAYMENT UNTIL: 20130831

Year of fee payment: 12

LAPS Cancellation because of no payment of annual fees