WO2016009465A1 - 空域情報処理装置、空域情報処理方法、空域情報処理プログラムが格納された非一時的なコンピュータ可読媒体 - Google Patents

空域情報処理装置、空域情報処理方法、空域情報処理プログラムが格納された非一時的なコンピュータ可読媒体 Download PDF

Info

Publication number
WO2016009465A1
WO2016009465A1 PCT/JP2014/003783 JP2014003783W WO2016009465A1 WO 2016009465 A1 WO2016009465 A1 WO 2016009465A1 JP 2014003783 W JP2014003783 W JP 2014003783W WO 2016009465 A1 WO2016009465 A1 WO 2016009465A1
Authority
WO
WIPO (PCT)
Prior art keywords
closed curve
line segment
point
airspace
transfer image
Prior art date
Application number
PCT/JP2014/003783
Other languages
English (en)
French (fr)
Inventor
雅彦 石田
Original Assignee
日本電気株式会社
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 日本電気株式会社 filed Critical 日本電気株式会社
Priority to JP2016533997A priority Critical patent/JP6233517B2/ja
Priority to PCT/JP2014/003783 priority patent/WO2016009465A1/ja
Priority to US15/326,695 priority patent/US20170206663A1/en
Publication of WO2016009465A1 publication Critical patent/WO2016009465A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/11—Region-based segmentation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20—Instruments for performing navigational calculations
    • G—PHYSICS
    • G08—SIGNALLING
    • G08G—TRAFFIC CONTROL SYSTEMS
    • G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
    • G08G5/0047—Navigation or guidance aids for a single aircraft
    • G08G5/006—Navigation or guidance aids for a single aircraft in accordance with predefined flight zones, e.g. to avoid prohibited zones

Definitions

  • an air route such as an aircraft can be expressed using a line segment connecting two points on a true sphere.
  • the aircraft ensures safety by flying in the airspace permitted to operate among the airspaces set in the sky.
  • an appropriate airspace design is necessary for ensuring the safety of the aircraft.
  • the examination direction of the sides of the polygon for defining the airspace is taken into consideration, and it is determined which area on the left and right of the circulation direction is the airspace.
  • Japanese Patent Application No. 2013-271712 which is a Japanese patent application, proposes a method of detecting intersections of line segments constituting airspaces and determining inside / outside airspaces for various airspaces.
  • An airspace information processing apparatus includes one or a plurality of line segments on a spherical surface, and all or part of a closed curve representing the contour of the airspace is on the spherical surface so as not to intersect with the closed curve.
  • Transfer means for generating a transfer image transferred to another position, and a determination for having an intersection point with the transfer image without having an intersection with another line segment forming the closed curve from the line segment forming the closed curve
  • Line segment generating means for generating a line segment
  • air space recognition means for recognizing an area where the line segment exists among the two areas on the spherical surface delimited by the closed curve as the air space.
  • a non-transitory computer-readable medium in which an airspace information processing program according to an aspect of the present invention is stored is composed of one or a plurality of line segments on a spherical surface, and all or a part of a closed curve representing an airspace contour is represented by the closed curve.
  • the line segment is present in the two regions on the spherical surface separated by the closed curve in the process of generating a determination line segment having an intersection with the transfer image and the air space recognition unit.
  • a process for recognizing an area as the airspace is executed by a computer.
  • FIG. 1 It is a figure which shows the line segment production
  • Embodiment 1 The airspace information processing apparatus 100 according to the first embodiment will be described.
  • the airspace information processing apparatus 100 is an apparatus that collectively handles information on a plurality of airspaces that are not defined in the circulation direction and are delimited by a closed curve composed of one or more line segments on a spherical surface.
  • the airspace information processing apparatus 100 is configured using hardware resources such as a computer system.
  • FIG. 1 is a diagram showing a line segment L connecting between the point P 1 and point P 2 of the perfect sphere CB.
  • V a is a unit normal vector with respect to the plane PL1 to which the line segment L connecting the point P 1 and the point P 2 belongs.
  • the plane PL1 is a plane including the center of the true sphere CB.
  • EQ is the equator of the true sphere CB.
  • Unit normal vector V a with respect to the plane PL1 is expressed by the following equation (1).
  • a point on the line segment L can be defined as a point that also satisfies the following expression (3).
  • FIG. 2 is a diagram showing a circle CC1 on the true sphere CB.
  • Circle CC1 on sphericity CB can be understood as a collection of point a distance r from a point P 0.
  • Position vector of the point P on the circumference of the circle CC1 satisfy each vector equation of the following equation using the position vector of the point P 0 (4).
  • R represents the radius of the true sphere CB.
  • V d is the unit normal vector of a plane circle CC1 belongs coincides with the position vector of the point P 0.
  • s d is a cosine of an angle formed by the point P 0 and the point P on the true sphere CB, and is expressed by the following equation (5).
  • the arc on the true sphere CB can be understood as a set of points having a distance r from the point P 0 on the true sphere CB.
  • FIG. 3 is a diagram showing an arc CC2 on the true sphere CB whose counterclockwise direction is from the start point to the end point.
  • R represents the radius of the true sphere CB.
  • V e is a unit normal vector of the plane to which the arc CC2 belongs, and coincides with the position vector of the point P 0 .
  • s d is the cosine of the angle formed by the point P 0 and the point P on the true sphere, and is expressed by the following equation (7).
  • FIG. 4 is a diagram showing an arc CC3 on a true sphere CB whose direction from the start point to the end point is clockwise.
  • R represents the radius of the true sphere CB.
  • V e is the unit normal vector of a plane circular arc CC3 belongs, is the position vector in the opposite direction of the point P 0.
  • s e is equal to the cosine of the angle formed by the point P 0 and the arbitrary point P on the arc on the true sphere CB, has a negative sign, and is expressed by the following formula (9).
  • FIG. 5 is a diagram illustrating an example of the airspace provided on the true sphere CB.
  • the airspace A is surrounded by a closed curve composed of line segments L A1 to L A4 so as to be separated from the external region. Since the airspace in FIG. 5 is merely an example, the number of line segments surrounding the airspace A can be one (that is, a circle on the true sphere CB) or a plurality other than four. In the example of FIG.
  • Line segment information Designation of one or more line segments surrounding the airspace.
  • Direction information Designation of the direction (counterclockwise, clockwise) when a closed curve composed of line segments surrounding the airspace is viewed from outside the true sphere.
  • the airspace information processing apparatus 100 handles a considerably large airspace on the true sphere. Therefore, it is necessary to collectively handle airspace information created by different entities such as organizations, corporations, and states.
  • the start point and end point of the line segment (for example, points P 1 and P 2 in FIG. 1) may be given as the line segment information.
  • information for designating the route of the line segment such as the above equation (3) may be added to the line segment information. That is, the line segment information is uniquely determined mathematically. Therefore, even when the airspace definition rules differ among organizations, corporations, nations, etc. that handle airspace, it is sufficient to express the line segments surrounding the airspace in some way, so the difference in line segment information does not matter.
  • the direction of the closed curve is artificially determined. Therefore, the direction of the closed curve may differ among organizations, corporations, nations, etc. that handle airspace. For example, it can be assumed that the direction of the closed curve is defined as counterclockwise in Country A and the direction of the closed curve is defined as clockwise in Country B. In this case, since the direction of the closed curve is defined to be counterclockwise in the system using the airspace information of country A, let's input the line segment information created in country B to the system of country A to recognize the airspace. Then, the airspace that Country B's line segment information is trying to indicate would be recognized as outside the airspace in the country A system. That is, in such a case, there arises a problem that the airspace is erroneously recognized.
  • direction information for each piece of line segment information created by different entities such as organizations, corporations, and countries.
  • entities such as organizations, corporations, and countries.
  • the airspace covering the wide area like the airspace information processing apparatus 100 according to the present embodiment is handled, so the direction in which the direction of the closed curve is designated in the line segment information that designates the airspace It does not have a function to add information. Even if the direction information is added, not only the information input to the system is increased, but the same problem as described above occurs if the direction information is incorrectly specified.
  • the area of the airspace delimited by the closed curve is usually smaller than half the surface area on the earth, as is clear from its use. Therefore, if the area of the airspace is compared with the area of the area outside the airspace, it can be determined that the smaller area is the airspace.
  • the amount of calculation becomes enormous, and it is not suitable as a process for simply recognizing the airspace. In particular, when handling a plurality of airspaces, the system simply performs enormous calculations just for recognizing the airspace, which is not realistic.
  • the airspace information processing apparatus 100 can recognize the airspace accurately with a small amount of calculation from airspace information in which the direction of the closed curve is not unified.
  • the airspace information processing apparatus 100 will be specifically described below.
  • FIG. 6 is a diagram schematically illustrating a basic configuration of the airspace information processing apparatus 100 according to the first embodiment.
  • the airspace information processing apparatus 100 includes a transfer unit 2, a line segment generation unit 3, and an airspace recognition unit 4.
  • FIG. 7 is a diagram schematically illustrating a configuration example of the airspace information processing apparatus 100 to which the configuration of the peripheral device is added.
  • the transfer unit 2 in addition to the transfer unit 2, the line segment generation unit 3, and the airspace recognition unit 4 shown in FIG. 6, a closed curve reading unit 1 and a storage unit 5 are added.
  • the transfer unit 2 includes a transfer processing unit 21 and an intersection detection unit 22.
  • the airspace information processing apparatus 100 is configured to divide a determination target closed curve based on a positional relationship between a closed curve (determination target closed curve and a transfer image) indicating the outlines of two spatially separated airspaces and a line segment drawn therebetween. The inside / outside determination is performed.
  • FIG. 8 is a flowchart showing the airspace information processing operation of the airspace information processing apparatus 100 according to the first embodiment.
  • Step S1 Reading of the determination target closed curve AZ1
  • the closed curve reading unit 1 reads the determination target closed curve AZ1.
  • the determination target closed curve AZ1 is not given a circular direction, and the determination target closed curve AZ1 only shows the outline of the airspace, and any of the two regions on the true sphere CB divided by the determination target closed curve AZ1 It is unknown whether it is an airspace.
  • the closed curve reading unit 1 reads line segment information defining the determination target closed curve AZ1 stored in the storage unit 5 in advance. In the case of the example shown in FIG. 5, the closed curve reading unit 1 reads information indicating the line segments L A1 to L A4 constituting the closed curve indicating the airspace A.
  • the closed curve reading unit 1 can output information indicating the read determination target closed curve AZ1 to the transfer unit 2 and the line segment generation unit 3.
  • Step S2 Generation of Transfer Image (Reverse Transfer Image) AZ2
  • the transfer processing unit 21 of the transfer unit 2 generates a transfer image AZ2 in which the closed curve reading unit 1 is transferred to another position on the true sphere CB.
  • the transfer processing unit 21 generates a reverse transfer image obtained by transferring the determination target closed curve AZ1 at a point-symmetrical position with respect to the center of the true sphere CB as a transfer image AZ2.
  • FIG. 9 is a diagram illustrating a relationship between the determination target closed curve AZ1 and the transfer image AZ2.
  • the transfer image AZ2 (displayed by a broken line) exists on the back side with the center O of the true sphere CB interposed therebetween.
  • Step S3 Intersection detection
  • the airspace information processing apparatus 100 performs inside / outside determination of the determination target closed curve AZ1 based on the positional relationship between two spatially separated airspaces and a line segment drawn between them. It is. Therefore, it is necessary to ensure that the determination target closed curve AZ1 and the transfer image AZ2 are spatially separated. Therefore, here, the transfer processing unit 21 of the transfer unit 2 determines whether or not the determination target closed curve AZ1 and the transfer image AZ2 have an intersection. Note that the intersection here does not include the contact point between the determination target closed curve AZ1 and the transfer image AZ2. In other words, when the determination target closed curve AZ1 and the transfer image AZ2 have an intersection, the determination of the rotation direction is impossible, and the processing is stopped.
  • Step S4 Line Segment Generation
  • the line generation unit 3 determines A line segment passing through the transfer image AZ2 is generated from a point on the line segment closest to the transfer image AZ2 among the line segments L A1 to L A4 of the target closed curve AZ1.
  • FIG. 10 is a flowchart of line segment generation processing in the airspace information processing apparatus 100 according to the first embodiment.
  • Step S41 An arbitrary point P0 (also referred to as a first point) is set on an arbitrary line segment that constitutes the airspace.
  • Step S42 From the point P0, a temporary line segment Lp (also referred to as a first line segment) having an intersection with the line segment constituting the transfer image AZ2 is drawn.
  • Lp also referred to as a first line segment
  • Step S43 The intersection of the line segment Lp and the line segment of the determination target closed curve AZ1 other than the line segment where the point P0 is set is obtained.
  • Step S44 Among the intersection points obtained above, the one closest to the transfer image AZ2 is selected as the point PA.
  • the intersection point includes a point P0 that is an end point of the line segment Lp.
  • Step S45 of the provisional line segment Lp a section between the point PA and any point on the transfer image AZ2 is set as the determination line segment Ld.
  • the definition of any point on the transfer image AZ2 is not limited to this.
  • FIGS. 11 to 17 are diagrams showing line segment generation in a crescent moon-shaped airspace sandwiched between two arcs.
  • 14 and 15 are diagrams showing line segment generation in a circular airspace.
  • 16 and 17 are diagrams showing line segment generation in a rectangular airspace surrounded by four line segments.
  • Step S15 Airspace recognition
  • the two areas divided by the closed curve representing the airspace are defined as A1 on the left side and A2 on the right side when the boundary line is circulated in the defined direction of the airspace. Since it is known that the transfer image AZ2 is outside the determination target closed curve AZ1, it is clear that the determination line segment Ld is emitted outward from the line segment defining the determination target closed curve AZ1.
  • the determination line segment Ld exists in the right side, that is, in the right area A2 when viewed from the line segment having the intersection with the determination line segment Ld (that is, the point PA), it is determined that the left area A1 indicates an air space. be able to.
  • the determination line segment Ld exists in the left side, that is, the left side area A1 when viewed from the line segment having the intersection with the determination line segment Ld (that is, the point PA), it is determined that the right side area A2 indicates an air space. Can do.
  • step S5 by determining whether the determination line segment Ld is in the left or right region of the closed curve (the line segment constituting the airspace), the left or right region of the closed curve is the determination target closed curve. Whether AZ1 is represented can be recognized.
  • the revolving direction of the recognized airspace may be set so as to match the revolving direction of the closed curve set in the airspace information processing apparatus 100.
  • the circulation direction of the airspace is defined as counterclockwise, the direction in which the determination line segment Ld is viewed to the right is the circulation direction.
  • the circulation direction of the airspace is defined as clockwise, the direction in which the determination line segment Ld is viewed to the left is the circulation direction.
  • the transfer image AZ2 is generated.
  • the entire airspace is not necessarily transferred, and only a part of the closed curve constituting the determination target closed curve AZ1 may be transferred.
  • a part of the closed curve to be transferred may be a point instead of a line segment.
  • a line segment Lp passing through the transferred line segment or the transferred point may be generated.
  • the position where the transferred point exists on the line segment Lp is also referred to as an intersection point for convenience.
  • this is limited to the case where it is clear that the transfer point is not included in the determination target closed curve AZ1.
  • the above intersection detection (step S12) can be omitted, which is advantageous from the viewpoint of reducing the amount of calculation.
  • Embodiment 2 An air space information processing apparatus according to the second embodiment will be described.
  • a modified example of the method for generating the transfer image AZ2 will be described.
  • the reverse transfer image of the determination target closed curve AZ1 is used as the transfer image AZ2.
  • the transfer image AZ2 can be used as long as it does not have an intersection with the determination target closed curve AZ1, a modification of the method for generating the transfer image AZ2 can be applied.
  • Modified example 1 Airspace center-of-gravity method
  • the center-of-gravity point G of the determination target closed curve AZ1 is obtained, and a vector OG connecting the center-of-gravity point G and the center O of the true sphere CB is obtained.
  • a copy obtained by rotating and copying the determination target closed curve AZ1 by a predetermined angle (for example, 90 °, 120 °, 180 °, etc.) with a vector perpendicular to the vector OG passing through the center O of the true sphere CB as a rotation axis is transferred.
  • a predetermined angle for example, 90 °, 120 °, 180 °, etc.
  • Modification 2 Vector Average Method For example, a plurality of points (XYZ orthogonal coordinates) are set at equal intervals on the closed curve surrounding the determination target closed curve AZ1, and an average vector of the position vectors of the set plurality of points is obtained. Then, a copy obtained by rotating and copying the determination target closed curve AZ1 by a predetermined angle (for example, 90 °, 120 °, 180 °, etc.) with a vector perpendicular to the average vector passing through the center O of the sphere CB as a rotation axis is transferred. Let it be image AZ2. In this case, the calculation of the average vector is easier than the calculation of the center-of-gravity point G of the determination target closed curve AZ1, so that the amount of calculation can be reduced.
  • a predetermined angle for example, 90 °, 120 °, 180 °, etc.
  • Modification 3 Average Latitude / Longitude For example, a plurality of points are set at equal intervals on a closed curve surrounding the determination target closed curve AZ1, and an average latitude / longitude coordinate composed of average values of the latitude / longitude of the plurality of set points is obtained. In addition, a vector connecting the average latitude and longitude coordinates and the center O of the true sphere CB is obtained. Then, the determination target closed curve AZ1 is rotated and copied by a predetermined angle (for example, 90 °, 120 °, 180 °, etc.) with a vector passing through the center O of the true sphere CB and perpendicular to the obtained vector as a rotation axis. This is designated as a transfer image AZ2. In this case, the calculation of the average latitude and longitude coordinates is easier than the calculation of the barycentric point G of the determination target closed curve AZ1, so that the calculation amount can be reduced.
  • a predetermined angle for example, 90 °, 120 °, 180 °,
  • Step S15 When the determination target closed curve AZ1 and the transfer image AZ2_X do not have an intersection, the transfer image AZ2_X is set as the transfer image AZ2, and the process ends.
  • the rotation amount is arbitrary.
  • the determination target closed curve AZ1 and the transfer image AZ2 may be rotated a plurality of times until there is no intersection.
  • Modifications 1 to 5 may be used in appropriate combination.
  • Expression (19) includes the square root of the discriminant D. Therefore, the solution of the equation (14) representing the intersection point Pc needs to be classified according to the value of the discriminant D.
  • FIG. 34 is a diagram illustrating a case where the reference circles match and two line segments are separated.
  • the four points of the start point PS1 of the line segment L1, the end point PE1 of the line segment L1, the start point PS2 of the line segment L2, and the end point PE2 of the line segment L2 are intersection points Pc .
  • FIG. 35 is a diagram illustrating a case where the reference circles match and the start point of one line segment overlaps the end point of the other line segment.
  • the intersection point Pc is a point that is the start point PS1 of the line segment L1 and is also the end point PE2 of the line segment L2, the end point PE1 of the line segment L1, and the start point PS2 of the line segment L2.
  • the intersection of the reference circle C 1 and the reference circle C 2 is not necessarily an intersection between the line segment L 1 and the line segment L 2. Therefore, in order to distinguish the intersection of the reference circle C 1 and the reference circle C 2 and the intersection of the line segment L 1 and the line segment L 2 , the intersection of the reference circle C 1 and the reference circle C 2 detected above. Are referred to as candidate points.
  • Step SS4 When the discriminant D is larger than 0, the detection unit 35 calculates a candidate point P c1 .
  • Step SS5 The detection unit 35 performs an intersection determination process for the candidate point P c1 .
  • the intersection determination process will be described later.
  • Step SS10 The detection unit 35 performs an intersection determination process for the candidate point P c0 .
  • the intersection determination process will be described later.
  • Step SS13 The detection unit 35 performs an intersection determination process for the start point PS2 of the line segment L2.
  • FIG. 43 is a flowchart showing intersection determination processing.
  • Step SR2 A range verification process is performed to determine whether the determination target point PJ exists on the line segment L1. Details of the range verification process will be described later. If the determination target point PJ does not exist on the line segment L1, the process ends.
  • Step SR4 When the determination target point PJ exists on the line segment L1 and L2, the determination target point PJ is registered as a candidate point.
  • FIG. 44 is a flowchart showing the range verification process.
  • LJ line segment to be verified
  • Step SA1 It is determined whether the determination target line segment LJ is a circle.
  • Step SA2 If the determination target line segment LJ is not a circle, it is determined whether the line segment is a superior arc.
  • Step SA4 When the determination target line segment LJ is inferior arc, it is determined whether or not both Expression (28) and Expression (29) are satisfied. When both Expression (28) and Expression (29) are satisfied, the determination target point PJ is on the determination target line segment LJ (YES determination). When at least one of Expression (28) and Expression (29) is not satisfied, the determination target point PJ does not exist on the determination target line segment LJ (NO determination).
  • the intersection detection unit 22 configures a closed curve that surrounds the determination target closed curve AZ1 by detecting the intersection point of the two line segments. It is specifically and easily determined whether or not the determination target closed curve AZ1 and the transfer image AZ2 have an intersection by applying to each of the line segments to be performed and each of the line segments constituting the closed curve surrounding the transfer image AZ2. It can be understood that it can be detected.

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Image Analysis (AREA)
  • Image Generation (AREA)
  • Traffic Control Systems (AREA)

Abstract

　空域情報処理装置は、転写部（２）、線分生成部（３）及び空域認識部（４）を有する。転写部（２）は、球面ＣＢ上の１又は複数の線分からなり、空域の輪郭を現す判定対象閉曲線（ＡＺ１）の全部又は一部を、判定対象閉曲線（ＡＺ１）と交点を有しないように球面ＣＢ上の別の位置に転写した転写像（ＡＺ２）を生成する。線分生成部（３）は、判定対象閉曲線（ＡＺ１）を構成する線分から、判定対象閉曲線（ＡＺ１）を構成する他の線分と交点を有することなく、かつ、転写像（ＡＺ２）と交点を有する判定用線分（Ｌｄ）を生成する。空域認識部（４）は、判定対象閉曲線（ＡＺ１）によって区切られる球面ＣＢ上の２つの領域のうち、判定用線分（Ｌｄ）が存在する領域を、空域として認識する。

Description

空域情報処理装置、空域情報処理方法、空域情報処理プログラムが格納された非一時的なコンピュータ可読媒体
　本発明は空域情報処理装置、空域情報処理方法、空域情報処理プログラムが格納された非一時的なコンピュータ可読媒体に関する。
　今日、地球上での移動体監視を行うため、様々なナビゲーションシステムが運用されている。輸送機械の中でも移動距離が長大である航空機の運航を管理するには、広範囲での方位、距離の算出を行う必要が有る。航空機のナビゲーションシステムでは、一般に、国家の領土及び領空、又は飛行情報区（ＦＩＲ：Flight Information Region）程度における広域の範囲で、大縮尺の空間情報を正確かつ効率的に処理することが求められる。
　例えば、航空機等の空路は、真球上の２地点間を結ぶ線分を用いて表現することができる。この際、航空機等の安全確保のため、２本の空路が交差するか否かを判定することが極めて重要である。また、航空機は、上空に設定された空域のうち、運航が許可された空域を飛行することで、安全を確保している。この際、隣接する空域が重複してしまうと、重複した領域に複数の航空機が進入してしまい、安全確保の観点から問題である。よって、上記のナビゲーションシステムでは、航空機の安全確保のため、適切な空域設計が必要である。
　一例として、地球上の多角形に対する任意の点の内外判定を行う位置関係判定方法が提案されている。この例では、空域を定義するための多角形の辺の調査方向（換言すれば、閉曲線の周回方向）を考慮し、周回方向に対して左右いずれの領域が空域であるかを判定する。
　また、日本国特許出願である特願２０１３－２７１７１２では、様々な空域について、空域を構成する線分の交点検出や、空域の内外判定を行う手法が提案されている。
特開２０１２－８８９０２号公報
　しかし、発明者は上述の手法には以下に示す問題点が有ることを見出した。飛行方式や空域設計仕様によっては、国家間や大陸間にまたがる大きな空域を取り扱うことが求められる場合が有る。この場合、例えば、国ごと又は空域ごとに、空域を定義している閉曲線の周回方向が異なる場合が想定し得る。これに対し、特許文献１の手法は、閉曲線の周回方向（多角形の辺の調査方向）を考慮しているものの、取り扱う空域を定義する閉曲線の方向が統一されていない場合の取り扱いについては考慮していない。閉曲線の周回方向が異なる空域が混在している複数の空域を、特許文献１の手法で取り扱うと、周回方向の相違により空域の内外を誤って認識するなど、空域設計上許容できないエラーが生じるおそれが有る。
　本発明は、上記の事情に鑑みて成されたものであり、本発明の目的は、周回方向の定められていない複数の空域を統合して取り扱うことである。
　本発明は、上記の事情に鑑みて成されたものであり、本発明の目的は、地球上の任意の形状、大きさの領域での位置関係の判定を正確かつ高精度に行うことである。
　本発明の一態様である空域情報処理装置は、球面上の１又は複数の線分からなり、空域の輪郭を現す閉曲線の全部又は一部を、前記閉曲線と交点を有しないように前記球面上の別の位置に転写した転写像を生成する転写手段と、前記閉曲線を構成する線分から、前記閉曲線を構成する他の線分と交点を有することなく、かつ、前記転写像と交点を有する判定用線分を生成する線分生成手段と、前記閉曲線によって区切られる前記球面上の２つの領域のうち、前記線分が存在する領域を、前記空域として認識する空域認識手段と、を備えるものである。
　本発明の一態様である空域情報処理方法は、転写手段に、球面上の１又は複数の線分からなり、空域の輪郭を現す閉曲線の全部又は一部を、前記閉曲線と交点を有しないように前記球面上の別の位置に転写した転写像を生成させ、線分生成手段に、前記閉曲線を構成する線分から、前記閉曲線を構成する他の線分と交点を有することなく、かつ、前記転写像と交点を有する判定用線分を生成させ、空域認識手段に、前記閉曲線によって区切られる前記球面上の２つの領域のうち、前記線分が存在する領域を、前記空域として認識させるものである。
　本発明の一態様である空域情報処理プログラムが格納された非一時的なコンピュータ可読媒体は、球面上の１又は複数の線分からなり、空域の輪郭を現す閉曲線の全部又は一部を、前記閉曲線と交点を有しないように前記球面上の別の位置に転写した転写像を生成する処理と、線分生成手段に、前記閉曲線を構成する線分から、前記閉曲線を構成する他の線分と交点を有することなく、かつ、前記転写像と交点を有する判定用線分を生成する処理と、空域認識手段に、前記閉曲線によって区切られる前記球面上の２つの領域のうち、前記線分が存在する領域を、前記空域として認識する処理と、をコンピュータに実行させるものである。
　本発明によれば、周回方向の定められていない複数の空域を統合して取り扱うことができる。
真球上の２地点間を結ぶ線分を示す図である。 真球上の円を示す図である。 始点から終点への方向が反時計回りの真球上の円弧を示す図である。 始点から終点への方向が時計回りの真球上の円弧を示す図である。 真球上に設けられた空域の例を示す図である。 実施の形態１にかかる空域情報処理装置の基本構成を模式的に示す図である。 周辺装置の構成を追加した実施の形態１にかかる空域情報処理装置の構成例を模式的に示す図である。 実施の形態１にかかる空域情報処理装置の空域情報処理動作を示すフローチャートである。 判定対象閉曲線と転写像との関係を示す図である。 実施の形態１にかかる空域情報処理装置での線分生成処理を示すフローャートである。 ２本の円弧で挟まれた三日月形の空域における線分生成を示す図である。 ２本の円弧で挟まれた三日月形の空域における線分生成を示す図である。 ２本の円弧で挟まれた三日月形の空域における線分生成を示す図である。 円形の空域における線分生成を示す図である。 円形の空域における線分生成を示す図である。 ４本の線分で囲まれた矩形の空域における線分生成を示す図である。 ４本の線分で囲まれた矩形の空域における線分生成を示す図である。 ３軸回転法の処理を示すフローチャートである。 実施の形態３にかかる交点検出部の構成を模式的に示すブロック図である。 基本形状データベースに含まれる情報を示す図である。 空域情報データベースに含まれる情報を示す図である。 演算部の基本的構成を模式的に示すブロック図である。 交点検出部の交点検出動作を示すフローチャートである。 真球上の始点から終点への方位が東向きである場合を示す図である。 真球上の始点から終点への方位が西向きである場合を示す図である。 真球上の線分を示す図である。 真球上の２本の線分を示す図である。 ２つの基準円が２つの交点を有する（交接する）場合を示す図である。 ２つの基準円が分離の関係にある場合を示す図である。 ２つの基準円が内包の関係にある場合を示す図である。 ２つの基準円が外接の関係にある場合を示す図である。 ２つの基準円が内接の関係にある場合を示す図である。 ２つの基準円が一致する場合を示す図である。 基準円が一致し、かつ、２本の線分が分離している場合を示す図である。 基準円が一致し、かつ、一方の線分の始点と他方の線分の終点とが重なる場合を示す図である。 基準円が一致し、かつ、２本の線分間に１つの重複部分がある場合を示す図である。 基準円が一致し、かつ、一方の線分の始点と他方の線分の終点とが重なるとともに２本の線分間に１つの重複部分が有る場合を示す図である。 基準円が一致し、かつ、２本の線分間に２つの重複部分が有る場合を示す図である。 中心角Ψが２πである場合（Ψ＝２π）の線分を示す図である。 中心角Ψがπ以上かつ２πよりも小さい場合（π≦Ψ＜２π）の線分を示す図である。 中心角Ψがπより小さい場合（０＜Ψ＜π）の線分を示す図である。 交点検出部での線分の交点検出動作を示すフローチャートである。 交点判定処理を示すフローチャートである。 範囲検証処理を示すフローチャートである。
　実施の形態１
　実施の形態１にかかる空域情報処理装置１００について説明する。空域情報処理装置１００は、球面上において１又は複数の線分からなる閉曲線で区切られる周回方向の定められていない複数の空域の情報を統合して取り扱う装置である。空域情報処理装置１００は、例えばコンピュータシステムなどのハードウェア資源を用いて構成される。
　まず、空域を理解するための前提として、閉曲線を構成する線分について説明する。真球上での線分は、以下の３つに大別できる。
真球上における２地点を最短で結ぶ線分
　真球ＣＢ上（地表面上）の地点Ｐ１と地点Ｐ２との間を結ぶ線分について説明する。図１は、真球ＣＢ上の地点Ｐ１と地点Ｐ２との間を結ぶ線分Ｌを示す図である。Ｖａは、地点Ｐ１と地点Ｐ２との間を結ぶ線分Ｌが属する平面ＰＬ１に対する単位法線ベクトルである。平面ＰＬ１は、真球ＣＢの中心を含む面である。ＥＱは真球ＣＢの赤道である。平面ＰＬ１に対する単位法線ベクトルＶａは、以下の式（１）で表される。
Figure JPOXMLDOC01-appb-M000001
　真球ＣＢ上の地点Ｐ１と地点Ｐ２とを結ぶ線分Ｌ上の点をＰ、ｓａを単位法線ベクトルＶａと点Ｐの位置ベクトルとがなす角の余弦とすると、ｓａは以下の式（２）で表される。
Figure JPOXMLDOC01-appb-M000002
　単位法線ベクトルＶａと線分Ｌとが直交することは明らかであるので、余弦Ｓａは０となる。したがって、線分Ｌ上の点をＰは、以下の式（３）も満たす点として定義できる。
Figure JPOXMLDOC01-appb-M000003
真球上の円
　真球ＣＢ上での円について説明する。図２は、真球ＣＢ上の円ＣＣ１を示す図である。真球ＣＢ上の円ＣＣ１は、ある点Ｐ０からの距離ｒの地点の集合として理解することができる。円ＣＣ１の円周上の点Ｐの位置ベクトルは、点Ｐ０の位置ベクトルを用いた以下の式（４）の各ベクトル方程式を満たす。Ｒは、真球ＣＢの半径を示す。Ｖｄは、円ＣＣ１が属する平面の単位法線ベクトルであり、点Ｐ０の位置ベクトルに一致する。
Figure JPOXMLDOC01-appb-M000004
 　ｓｄは、真球ＣＢ上において点Ｐ０と点Ｐとがなす角の余弦であり、以下の式（５）で表される。
Figure JPOXMLDOC01-appb-M000005
真球の２地点を結ぶ円弧
　真球ＣＢ上の円弧について説明する。真球ＣＢ上の円弧は、真球ＣＢ上の点Ｐ０から距離ｒの点の集合として理解することができる。
　円弧の始点から終点への方向が反時計回りの場合について説明する。図３は、始点から終点への方向が反時計回りの真球ＣＢ上の円弧ＣＣ２を示す図である。２点間の方向が反時計回りの場合、円弧ＣＣ２上の点Ｐの位置ベクトルは、以下の式（６）の各ベクトル方程式を満たす。Ｒは、真球ＣＢの半径を示す。Ｖｅは、円弧ＣＣ２が属する平面の単位法線ベクトルであり、点Ｐ０の位置ベクトルに一致する。
Figure JPOXMLDOC01-appb-M000006
 　ｓｄは、真球上において地点Ｐ０と点Ｐとがなす角の余弦であり、以下の式（７）で表される。
Figure JPOXMLDOC01-appb-M000007
　円弧の始点から終点への方向が時計回りの場合について説明する。図４は、始点から終点への方向が時計回りの真球ＣＢ上の円弧ＣＣ３を示す図である。２点間の方向が時計回りの場合、円弧ＣＣ３上の点Ｐの位置ベクトルは、以下の式（８）の各ベクトル方程式を満たす。Ｒは、真球ＣＢの半径を示す。Ｖｅは、円弧ＣＣ３が属する平面の単位法線ベクトルであり、点Ｐ０の位置ベクトルと逆方向である。
Figure JPOXMLDOC01-appb-M000008
 　ｓｅは、真球ＣＢ上において地点Ｐ０と円弧上の任意の点Ｐとがなす角の余弦に等しく、かつ符号が負であり、以下の式（９）で表される。
Figure JPOXMLDOC01-appb-M000009
　次いで、真球上に設定される空域について説明する。図５は、真球ＣＢ上に設けられた空域の例を示す図である。図５では、空域Ａが線分ＬＡ１～ＬＡ４からなる閉曲線で囲まれることで、外部の領域と区切られている。図５の空域はあくまで例示であるので、空域Ａを囲む線分の数は、１本（すなわち、真球ＣＢ上の円）又は４本以外の複数本とすることができる。図５の例では、線分ＬＡ１～ＬＡ４で構成される閉曲線を真球の外部から見て反時計回りに進む場合に、線分から見て左側に見える真球上の領域が、空域Ａとして定義される。よって、この場合、線分から見て右側に見える真球上の領域は、空域Ａの外側の領域として定義される。
　以上、要約すると、空域を定義する場合、以下の２つの情報が必要であることが理解できる。
（１）線分情報
　空域を囲む１又は複数の線分の指定。
（２）方向情報
　空域を囲む線分で構成される閉曲線を真球の外部から見た場合の向き（反時計回り、時計回り）の指定。
　しかし、本実施の形態にかかる空域情報処理装置１００では、真球上の相当程度に大きな空域を取り扱うことを想定している。そのため、団体、法人、国家等の異なる主体が作成した空域情報を一括して取り扱う必要が生じる。
　この場合、空域を囲む線分のそれぞれを指定するには、線分の始点及び終点（例えば、図１の点Ｐ１及びＰ２）を線分情報として与えればよい。また、始点と終点とを結ぶ経路が一意に定義されていない場合には、上述の式（３）のような線分の経路を指定する情報を線分情報に加えればよい。つまり、線分情報については、数学的に一意に定められるものである。よって、空域を取り扱う団体、法人、国家等の間で、空域の定義ルールが異なる場合でも、空域を囲む線分を何らかの方式で表現すれば足りるので、線分情報の相違が問題にはならない。
　これに対し、方向情報については、以下の理由で取扱いに注意を要する。方向情報については、閉曲線の向きは人為的に決められるものである。よって、閉曲線の向きは、空域を取り扱う団体、法人、国家等の間では異なることがある。例えば、Ａ国においては閉曲線の向きが反時計回りと規定され、Ｂ国においては閉曲線の向きが時計回りと規定されている場合が想定し得る。この場合、Ａ国の空域情報を用いるシステムでは閉曲線の向きは反時計回りであると定義されるので、Ｂ国で作成された線分情報をＡ国のシステムに入力して空域を認識させようとすると、Ｂ国の線分情報が示そうとした空域はＡ国のシステムでは空域の外側だと認識されてしまう。つまり、このような場合には、空域が誤って認識されてしまう問題が生じる。
　これを避けるには、団体、法人、国家等の異なる主体で作成された線分情報ごとに、方向情報を指定することも考え得る。しかし、既存のシステムでは、本実施の形態にかかる空域情報処理装置１００のような広域にわたる空域を取り扱うことを想定していないので、空域を指定する線分情報に、閉曲線の向きを指定する方向情報を付加する機能を有していない。また、方向情報を付加するにしても、システムに入力する情報が増えるだけでなく、方向情報の指定を誤れば上述と同様の問題が生じる。
　また、閉曲線で区切られた空域の面積は、その用途から明らかなように、通常は地球上の表面積の半分よりも小さい。よって、空域の面積と空域の外側の領域の面積とを比較すれば、面積の小さい方が空域であることを判別することはできる。しかし、球面上の閉曲線で区切られる領域の面積を求める場合、計算量が膨大になり、単に空域を認識させるための処理としては適当でない。特に、複数の空域を取り扱う場合、システムが単に空域を認識するためだけに膨大な計算を行うことになり、現実的でない。
　これに対し、本実施の形態にかかる空域情報処理装置１００は、閉曲線の向きが統一されていない空域情報から、少ない計算量にて正確に空域を認識できるものである。以下、空域情報処理装置１００について、具体的に説明する。
　図６は、実施の形態１にかかる空域情報処理装置１００の基本構成を模式的に示す図である。空域情報処理装置１００は、転写部２、線分生成部３及び空域認識部４を有する。
　図７は、周辺装置の構成を追加した空域情報処理装置１００の構成例を模式的に示す図である。図７では、図６で示した転写部２、線分生成部３及び空域認識部４の他に、閉曲線読込部１及び記憶部５が追加されている。なお、図２では、転写部２は、転写処理部２１と交点検出部２２とを有する。
　以下、本実施の形態にかかる空域情報処理装置１００の動作について説明する。空域情報処理装置１００は、空間的に離隔した２つの空域の輪郭を示す閉曲線（判定対象閉曲線及び転写像）とその間に引かれた線分との位置関係に基づいて、判定対象閉曲線が区切る空域の内外判定を行うものである。図８は、実施の形態１にかかる空域情報処理装置１００の空域情報処理動作を示すフローチャートである。
ステップＳ１：判定対象閉曲線ＡＺ１の読み込み
　まず、閉曲線読込部１が、判定対象閉曲線ＡＺ１を読み込む。この時点において、判定対象閉曲線ＡＺ１には周回方向は与えられておらず、判定対象閉曲線ＡＺ１は空域の輪郭を現すのみで、判定対象閉曲線ＡＺ１で区切られる真球ＣＢ上の２つの領域のいずれが空域であるかは不明である。具体的には、この際、閉曲線読込部１は、記憶部５に予め記憶されている判定対象閉曲線ＡＺ１を規定する線分情報を読み込む。図５に示す例の場合、閉曲線読込部１は、空域Ａを示す閉曲線を構成する線分ＬＡ１～ＬＡ４を示す情報を読み込む。閉曲線読込部１は、読み込んだ判定対象閉曲線ＡＺ１を示す情報を、転写部２と線分生成部３とに出力することができる。
ステップＳ２：転写像（反転転写像）ＡＺ２の生成
　転写部２の転写処理部２１は、閉曲線読込部１を真球上ＣＢ上の別の位置に転写した転写像ＡＺ２を生成する。本実施の形態では、転写処理部２１は、判定対象閉曲線ＡＺ１を真球上ＣＢの中心に対して点対称の位置に転写した反転転写像を、転写像ＡＺ２として生成する。図９は、判定対象閉曲線ＡＺ１と転写像ＡＺ２との関係を示す図である。図９では、判定対象閉曲線ＡＺ１が真球ＣＢの手前側に存在するので、転写像ＡＺ２（破線で表示）は、真球ＣＢの中心Ｏを挟んで裏側に存在する。
ステップＳ３：交点検出
　上述のとおり、空域情報処理装置１００は、空間的に離隔した２つの空域とその間に引かれた線分との位置関係に基づいて、判定対象閉曲線ＡＺ１の内外判定を行うものである。よって、判定対象閉曲線ＡＺ１と転写像ＡＺ２とが空間的に離隔していることが担保されている必要が有る。よって、ここで、転写部２の転写処理部２１は、判定対象閉曲線ＡＺ１と転写像ＡＺ２とが交点を有するか否かを判定する。なお、ここでいう交点とは、判定対象閉曲線ＡＺ１と転写像ＡＺ２との接点は含まないものとする。換言すれば判定対象閉曲線ＡＺ１と転写像ＡＺ２とが交点を有する場合、周回方向の判定は不能であるので、処理を中止する。
ステップＳ４：線分生成
　判定対象閉曲線ＡＺ１と転写像ＡＺ２とが空間的に離隔している場合（判定対象閉曲線ＡＺ１と転写像ＡＺ２とが交点を有しない場合）、線分生成部３は、判定対象閉曲線ＡＺ１の線分ＬＡ１～ＬＡ４のうちで転写像ＡＺ２に最も近い線分上の点から、転写像ＡＺ２を通過する線分を生成する。
　線分生成（ステップＳ１４）について、より詳しく説明する。図１０は、実施の形態１にかかる空域情報処理装置１００での線分生成処理を示すフローチャートである。
ステップＳ４１
　空域を構成する線分のうちの任意の線分上に、任意の点Ｐ０（第１の点とも称する）を設定する。
ステップＳ４２
　点Ｐ０から、転写像ＡＺ２を構成する線分と交点を有する仮の線分Ｌｐ（第１線分とも称する）を引く。
ステップＳ４３
　点Ｐ０を設定した線分以外の判定対象閉曲線ＡＺ１の線分と、線分Ｌｐとの交点を求める。
ステップＳ４４
　上記で求めた交点のうち、転写像ＡＺ２に最も近いものを、点ＰＡとして選定する。なお、ここでいう交点とは、線分Ｌｐの端点である点Ｐ０を含むものとする。
ステップＳ４５
　仮の線分Ｌｐのうちで、点ＰＡと転写像ＡＺ２上のいずれかの点との間の区間を、判定用線分Ｌｄとして設定する。ここでは、転写像ＡＺ２上のいずれかの点として、例えば、転写像ＡＺ２を構成する線分と仮の線分Ｌｐとの交点のうちで、判定対象閉曲線ＡＺ１に最も近い点（第２の点とも称する）を用いるものとする。但し、転写像ＡＺ２上のいずれかの点の定義は、これに限られない。
　以上のステップＳ４１～Ｓ４５により、上述したステップＳ１４での線分生成を行うことができる。線分生成の例を図１１～図１７に示す。図１１～図１７では、図面の簡略のため、空域を平面上で近似的に表している。図１１～図１３は、２本の円弧で挟まれた三日月形の空域における線分生成を示す図である。図１４及び図１５は、円形の空域における線分生成を示す図である。図１６及び図１７は、４本の線分で囲まれた矩形の空域における線分生成を示す図である。
ステップＳ１５：空域認識 
　図１０に戻り、空域情報処理装置１００の空域情報処理動作を引き続き説明する。
　空域を表す閉曲線で区切られる２領域を、空域の定義された方向に境界線を周回したときに左側にある領域をＡ１、右側にある領域をＡ２とする。転写像ＡＺ２は、判定対象閉曲線ＡＺ１の外部にあることが分かっているので、判定用線分Ｌｄは、判定対象閉曲線ＡＺ１を定義する線分から外側に向けて出射していることは明らかである。
　この場合、判定用線分Ｌｄとの交点（すなわち点ＰＡ）を有する線分から見て右側、すなわち、右側領域Ａ２に判定用線分Ｌｄが存在する場合、左側領域Ａ１が空域を示すと判定することができる。
　また、判定用線分Ｌｄとの交点（すなわち点ＰＡ）を有する線分から見て左側、すなわち、左側領域Ａ１に判定用線分Ｌｄが存在する場合、右側領域Ａ２が空域を示すと判定することができる。
　以上説明したように、ステップＳ５では、判定用線分Ｌｄが閉曲線（空域を構成する線分）の左右いずれの領域に存在するかを判定することで、閉曲線の左右いずれの領域が判定対象閉曲線ＡＺ１を表すかを認識することができる。
　その後、空域情報処理装置１００で設定されている閉曲線の周回方向に合うように、認識した空域の周回方向を設定してもよい。例えば、空域の周回方向を反時計回りとして定義する場合、判定用線分Ｌｄを右に見る方向が周回方向である。空域の周回方向を時計回りとして定義する場合、判定用線分Ｌｄを左に見る方向が周回方向である。
　なお、上述では、転写像ＡＺ２を生成するものとして説明したが、必ずしも空域全体を転写する必要はなく、判定対象閉曲線ＡＺ１を構成する閉曲線上の一部のみを転写してもよい。また、転写する閉曲線上の一部は線分でなく、点でもよい。そして、転写した線分または転写した点を通る線分Ｌｐを生成してもよい。なお、線分Ｌｐが転写した点を通る場合、線分Ｌｐ上で転写した点が存在する位置を、便宜上交点とも称する。ただし、これは、転写点が判定対象閉曲線ＡＺ１に含まれないことが明らかな場合に限られる。この場合、上述の交点検出（ステップＳ１２）を省略することができ、計算量削減の観点からは有利である。
　また、転写点に代えて、判定対象閉曲線ＡＺ１に含まれないことが明らかな別の点を用いてもよい。例えば、実用上、南極点を含む空域が設定されることは考え難いため、南極点を上述の別の点として用いることも考え得る。
　実施の形態２
　実施の形態２にかかる空域情報処理装置について説明する。本実施の形態では、転写像ＡＺ２の生成方法の変形例について説明する。実施の形態１では、判定対象閉曲線ＡＺ１の反転転写像を転写像ＡＺ２として用いた。しかし、転写像ＡＺ２は、判定対象閉曲線ＡＺ１と交点を有しないものであれば使用可能なので、転写像ＡＺ２の生成方法の変形例を適用することができる。
変形例１：空域重心法
　判定対象閉曲線ＡＺ１の重心点Ｇを求め、重心点Ｇと真球ＣＢの中心Ｏとを結ぶベクトルＯＧを求める。そして、真球ＣＢの中心Ｏを通るベクトルＯＧに垂直なベクトルを回転軸として、所定の角度（例えば、９０°、１２０°、１８０°など）だけ判定対象閉曲線ＡＺ１を回転複写したものを、転写像ＡＺ２とする。この場合、判定対象閉曲線ＡＺ１の重心点Ｇの算出に際し、相応の計算量が必要となる。
変形例２：ベクトル平均法
　例えば、判定対象閉曲線ＡＺ１を囲む閉曲線に等間隔に複数の点（ＸＹＺ直交座標）を設定し、設定した複数の点の位置ベクトルの平均ベクトルを求める。そして、真球ＣＢの中心Ｏを通る平均ベクトルに垂直なベクトルを回転軸として、所定の角度（例えば、９０°、１２０°、１８０°など）だけ判定対象閉曲線ＡＺ１を回転複写したものを、転写像ＡＺ２とする。この場合、判定対象閉曲線ＡＺ１の重心点Ｇの算出に比べて、平均ベクトルの算出は容易であるので、計算量の削減ができる。
変形例３：緯度経度平均法
　例えば、判定対象閉曲線ＡＺ１を囲む閉曲線に等間隔に複数の点を設定し、設定した複数の点の緯度経度の平均値で構成される平均緯度経度座標を求め、かつ、平均緯度経度座標と真球ＣＢの中心Ｏとを結ぶベクトルを求める。そして、真球ＣＢの中心Ｏを通り、かつ、求めたベクトルに垂直なベクトルを回転軸として、所定の角度（例えば、９０°、１２０°、１８０°など）だけ判定対象閉曲線ＡＺ１を回転複写したものを、転写像ＡＺ２とする。この場合、判定対象閉曲線ＡＺ１の重心点Ｇの算出に比べて、平均緯度経度座標の算出は容易であるので、計算量の削減ができる。
変形例４：構成点抽出法
　判定対象閉曲線ＡＺ１を囲む閉曲線上の任意の点と真球ＣＢの中心Ｏとを結ぶベクトルを求める。そして、真球ＣＢの中心Ｏを通り、かつ、求めたベクトルに垂直なベクトルを回転軸として、所定の角度（例えば、９０°、１２０°、１８０°など）だけ判定対象閉曲線ＡＺ１を回転複写したものを、転写像ＡＺ２とする。この場合、判定対象閉曲線ＡＺ１を囲む閉曲線上の任意の点を求めるだけで足りるので、計算量の削減ができる。
変形例５：構成点ペア抽出法
　判定対象閉曲線ＡＺ１を囲む閉曲線上に距離が最大となる２点を設定し、設定した２点間の中点を求める。さらに、中点と真球ＣＢの中心Ｏとを結ぶベクトルを求める。そして、真球ＣＢの中心Ｏを通り、かつ、求めたベクトルに垂直なベクトルを回転軸として、所定の角度（例えば、９０°、１２０°、１８０°など）だけ判定対象閉曲線ＡＺ１を回転複写したものを、転写像ＡＺ２とする。この場合、中点を求めるだけで足りるので、計算量の削減ができる。
変形例６：３軸回転法
　上記で説明した変形例における３次元での任意角度の座標回転は、３つの浮動小数点パラメータと三角関数演算を含むため、計算誤差が生じやすい。これに対し、変形例６は、原理的に計算誤差が生じない転写像ＡＺ２の生成方法にかかるものである。
　真球ＣＢに、例えば図５に示すように、ＸＹＺ軸を設定する。この場合、真球ＣＢを地球とした場合、Ｚ軸は地軸に該当する。なお、Ｘ軸を第２回転軸、Ｙ軸を第３回転軸、Ｚ軸を第１回転軸とも称する。第１回転軸（Ｚ軸）回りの回転は、Ｘ－Ｙ平面（Ｘ－Ｙ座標）がＺ軸を中心に回転することを意味する。第２回転軸（Ｘ軸）回りの回転は、Ｙ－Ｚ平面（Ｙ－Ｚ座標）がＸ軸を中心に回転することを意味する。第３回転軸（Ｙ軸）回りの回転は、Ｚ－Ｘ平面（Ｚ－Ｘ座標）がＹ軸を中心に回転することを意味する。
　以下、３軸回転法について具体的に説明する。図１８は、３軸回転法の処理を示すフローチャートである。図１８で説明するステップＳ１０～Ｓ１９は、図８のステップＳ１２及びＳ１３に対応するものである。
ステップＳ１０
　まず、判定対象閉曲線ＡＺ１をＺ軸回りに１８０°回転させ、転写像を生成する。このときの転写像ＡＺ２＿Ｚと表記する。この際、判定対象閉曲線ＡＺ１上の座標（ｘ,ｙ,ｚ）は、座標（－ｘ,－ｙ,ｚ）に転写される。すなわち、回転とは言いつつも、実際には判定対象閉曲線ＡＺ１を囲む閉曲線を定義する座標情報のｘ座標及びｙ座標の符号を反転させる簡易な演算を行うのみでよいことが理解できる。
ステップＳ１１
　判定対象閉曲線ＡＺ１と転写像ＡＺ２＿Ｚとが、交点を有するかを検出する。ここでの交点検出は、既述のステップＳ１２と同様の手法を用いることができる。
ステップＳ１２
　判定対象閉曲線ＡＺ１と転写像ＡＺ２＿Ｚとが交点を有しない場合、転写像ＡＺ２＿Ｚを転写像ＡＺ２として設定し、処理を終了する。
ステップＳ１３
　判定対象閉曲線ＡＺ１と転写像ＡＺ２＿Ｚとが交点を有する場合、判定対象閉曲線ＡＺ１をＸ軸回りに１８０°回転させ、新たな転写像を生成する。このときの転写像ＡＺ２＿Ｘと表記する。この際、判定対象閉曲線ＡＺ１上の座標（ｘ,ｙ,ｚ）は、座標（ｘ,－ｙ,－ｚ）に転写される。すなわち、回転とは言いつつも、実際には判定対象閉曲線ＡＺ１を囲む閉曲線を定義する座標情報のｙ座標及びｚ座標の符号を反転させる簡易な演算を行うのみでよいことが理解できる。
ステップＳ１４
　判定対象閉曲線ＡＺ１と転写像ＡＺ２＿Ｘとが、交点を有するかを検出する。ここでの交点検出は、既述のステップＳ１２と同様の手法を用いることができる。
ステップＳ１５
　判定対象閉曲線ＡＺ１と転写像ＡＺ２＿Ｘとが交点を有しない場合、転写像ＡＺ２＿Ｘを転写像ＡＺ２として設定し、処理を終了する。
ステップＳ１６
　判定対象閉曲線ＡＺ１と転写像ＡＺ２＿Ｘとが交点を有する場合、判定対象閉曲線ＡＺ１をＹ軸回りに１８０°回転させ、新たな転写像を生成する。このときの転写像ＡＺ２＿Ｙと表記する。この際、判定対象閉曲線ＡＺ１上の座標（ｘ,ｙ,ｚ）は、座標（－ｘ,ｙ,－ｚ）に転写される。すなわち、回転とは言いつつも、実際には判定対象閉曲線ＡＺ１を囲む閉曲線を定義する座標情報のｘ座標及びｚ座標の符号を反転させる簡易な演算を行うのみでよいことが理解できる。
ステップＳ１７
　判定対象閉曲線ＡＺ１と転写像ＡＺ２＿Ｙとが、交点を有するかを検出する。ここでの交点検出は、既述のステップＳ１２と同様の手法を用いることができる。
ステップＳ１８
　判定対象閉曲線ＡＺ１と転写像ＡＺ２＿Ｙとが交点を有しない場合、転写像ＡＺ２＿Ｙを転写像ＡＺ２として設定し、処理を終了する。
ステップＳ１９
　判定対象閉曲線ＡＺ１と転写像ＡＺ２＿Ｙとが交点を有する場合、転写像ＡＺ２の作成を中止し、処理を終了する。
　以上、３軸回転法は、上述の通り、回転とは言いつつも、実際には判定対象閉曲線ＡＺ１を囲む閉曲線を定義する座標情報の符号を反転させる簡易な演算を行うのみでよい。よって、上述の変形例１～５に比べて、計算量をより低減することができる。
　なお、実施の形態１で説明した反転転写像は、判定対象閉曲線ＡＺ１を真球ＣＢの中心を挟んで点対称な位置に転写して生成するため、判定対象閉曲線ＡＺ１の周回方向と反転転写像の周回方向とは逆転する。これに対し、３軸回転法では、周回方向をそのまま保存して転写像ＡＺ２を生成することができる。
　他にも、例えば変形例１～５において、回転量は任意である。また、判定対象閉曲線ＡＺ１と転写像ＡＺ２とが交点を持たなくなるまで、複数回回転させてもよい。また、複数回回転させる場合には、変形例１～５を適宜組み合わせて用いてもよい。
　実施の形態３
　実施の形態３にかかる地理情報管理装置について説明する。本実施の形態では、図８のステップＳ１３で説明した交点検出の具体例について説明する。図１９は、実施の形態３にかかる交点検出部２２の構成を模式的に示すブロック図である。交点検出部２２は、記憶装置３１、演算部３２及びバス３３を有する。交点検出部２２は、例えばコンピュータシステムなどのハードウェア資源を用いて構成される。
　記憶装置３１は、データが格納されたデータベースや、演算部３２での処理に供されるプログラムを記憶することが可能である。記憶装置３１は、例えばハードディスクドライブ、フラッシュメモリなどの各種の記憶装置を適用することが可能である。具体的には、記憶装置３１には、基本形状データベースＤ１、空域情報データベースＤ２が格納される。
　基本形状データベースＤ１は、予め与えられる固有の情報である。図２０は、基本形状データベースＤ１に含まれる情報を示す図である。基本形状データベースＤ１は、例えば真球ＣＢ（地球）の半径Ｒを含む。
　空域情報データベースＤ２は、真球ＣＢ上での線分や空域を示す座標情報を含む。図２１は、空域情報データベースＤ２に含まれる情報を示す図である。空域情報データベースＤ２は、真球ＣＢにおける航空機の座標Ｐ（Ｘ，Ｙ，Ｚ）や、２つの地点間を結ぶ線分（空路）、空域名、空域の形状（円、矩形など）や範囲を示す情報が含まれる。空域情報データベースＤ２は、例えばＰ（Ｘ，Ｙ，Ｚ）、線分の始点緯度／経度、線分の終点緯度／経度、空域形状、空域の範囲を表す線分（大圏、緯線、経線）、空域の範囲を表す円や円弧の情報、円を表すための中心緯度／経度及び半径を含む。
　また、記憶装置３１には、後述する線分の交点検出の演算処理を規定するプログラムＰＲＧ１を格納することも可能である。
　演算部３２は、記憶装置３１から、プログラム及びデータベースを読み出し、必要な演算処理を行うことができる。演算部３２は、例えばＣＰＵ（Central Processing Unit）により構成される。
　図２２は、演算部３２の基本的構成を模式的に示すブロック図である。演算部３２は、候補点検出部３４及び検出部３５を有する。候補点検出部３４及び検出部３５の詳細については、後述する。
　続いて、交点検出部２２の交点検出の動作について説明する。図２３は、交点検出部２２の交点検出動作を示すフローチャートである。
ステップＳ２１
　まず、演算部３２は、プログラムＰＲＧ１を読み込む。プログラムＰＲＧ１は、基本形状データベースＤ１及び空域情報データベースＤ２を用いて、真球ＣＢ上の２本の線分が交点を有するかを判定するためのプログラムである。これにより、演算部３２は、候補点検出部３４及び検出部を有する形状判定装置として機能する。プログラムＰＲＧ１は、例えば記憶装置３１から読み出される。
　なお、この例では演算部３２がコンピュータにより構成され、プログラムＰＲＧ１を読み込むものとして説明した。しかし演算部３２は、物理的実体を有する候補点検出部３４及び検出部が内部に形成された装置として構成できることは、いうまでもない。
ステップＳ２２
　次いで、演算部３２は、記憶装置３１から基本形状データベースＤ１及び空域情報データベースＤ２を読み出す。
ステップＳ２３
　演算部３２は、プログラムＰＲＧ１で規定される数式に、基本形状データベースＤ１及び空域情報データベースＤ２に含まれる情報を代入し、交点検出を行う。
　以下、ステップＳ２３での交点検出の詳細について、具体的に説明する。真球ＣＢ（地表面上）の地点を示すにあたり、以下の説明で用いる数式及び図では、ベクトル量には上付き矢印を付している。説明の簡略化のため、ベクトル量は全て規格化されている。具体的には、真球ＣＢ上の点を表す位置ベクトルは、基本形状データベースＤ１に含まれる真球ＣＢの半径Ｒで除することで、規格化位置ベクトルとなっている。以下では、説明の簡略化のため、規格化されたベクトルを単にベクトルと称する。
　真球ＣＢ上では、空域は、１又は複数の互いに交差しない線分で囲まれる領域として定義できる。一般に、真球ＣＢ上における線分は、円弧となる。円弧は、閉曲線である円上において始点と終点とで挟まれた区間として表すことができる。以下では、本実施の形態にかかる交点検出を理解するための前提として、真球ＣＢ上での線分の表現方法について説明する。
　真球上における２地点を最短経路で結ぶ線分、真球上の円、真球の２地点を結ぶ円弧については、実施の形態１で説明したので、省略する。
同一緯度の２地点間を結ぶ緯線
　真球ＣＢ上（地表面上）の同一緯度の地点Ｐ１と地点Ｐ２との間を結ぶ緯線について説明する。真球ＣＢ上（地表面上）の緯線は、真球ＣＢ上の同一緯度における２地点間の航程線として理解することができる。
　地点Ｐ１（始点）から地点Ｐ２（終点）への方位が東向きである場合について説明する。図２４は、真球ＣＢ上の地点Ｐ１から地点Ｐ２への方位が東向きである場合を示す図である。真球ＣＢ上の地点Ｐ１と地点Ｐ２とが存在する緯線上の点をＰとすると、点Ｐの位置ベクトルは、式（１０）に示す各ベクトル方程式を満足する。なお、Ｖｂは、地点Ｐ１と地点Ｐ２とが存在する緯線が属する平面ＰＬ２に対する単位法線ベクトルである。極点Ｎは、真球ＣＢの北極点である。平面ＰＬ２は緯線に対して平行であるので、単位法線ベクトルＶｂと極点Ｎの位置ベクトルは一致する。
Figure JPOXMLDOC01-appb-M000010
 　ｓｂは、地点Ｐ１及び地点Ｐ２の緯度θと赤道面とがなす角の正弦であり、以下の式（１１）で表される。
Figure JPOXMLDOC01-appb-M000011
　地点Ｐ１（始点）から地点Ｐ２（終点）への方位が西向きである場合について説明する。図２５は、真球ＣＢ上の地点Ｐ１から地点Ｐ２への方位が西向きである場合を示す図である。真球ＣＢ上の地点Ｐ１と地点Ｐ２とが存在する緯線上の点をＰとすると、点Ｐの位置ベクトルは、式（１２）に示す各ベクトル方程式を満足する。なお、Ｖｃは、地点Ｐ１と地点Ｐ２とが存在する緯線が属する平面ＰＬ３に対する単位法線ベクトルである。ここで、真球ＣＢの極点Ｓ（地球の南極）を定義する。極点Ｓを示す位置ベクトルは、以下の式（１２）で表される。平面ＰＬ３は緯線に対して平行であるので、単位法線ベクトルＶｃと極点Ｓを示す位置ベクトルは一致する。
Figure JPOXMLDOC01-appb-M000012
 　ｓｃは、地点Ｐ１及び地点Ｐ２の緯度θと赤道面とがなす角の正弦に等しく、かつ符号が逆であり、以下の式（１３）で表される。
Figure JPOXMLDOC01-appb-M000013
　次いで、交点検出における線分の取り扱いについて説明する。以下では、真球ＣＢ上の線分である円弧をその一部に含む円を基準円と称し、この場合、円弧が基準円に属すると称する。
　図２６は、真球ＣＢ上の線分Ｌを示す図である。この例では、真球ＣＢ上の円弧である線分Ｌが属する基準円をＣとする。また、基準円Ｃの周上の点をＰとする。基準円を真球ＣＢの上空から俯瞰した場合に、基準円の周上を始点ＰＳから反時計回りに終点ＰＥに至る経路を、基準円Ｃに属する線分Ｌと定義する。図２６及び以降の図では、北極をＮ、南極をＳと表示している。赤道を、ＥＱと表示している。
　基準円Ｃ上の点Ｐの位置ベクトルは、以下の式（１４）を満たす。式（１４）において、ｓは基準円Ｃの半径（曲率半径）を示すパラメータである。Ｖは基準円Ｃが属する平面に対する単位法線ベクトルである。
Figure JPOXMLDOC01-appb-M000014
　上述の前提のもと、真球ＣＢ上に２本の線分Ｌ１及びＬ２が存在する例について検討する。図２７は、真球ＣＢ上の２本の線分Ｌ１及びＬ２を示す図である。２本の線分Ｌ１及びＬ２を取り扱うにあたり、線分Ｌ１が属する基準円をＣ１、線分Ｌ２が属する基準円をＣ２とする。基準円Ｃ１の半径（曲率半径）を示すパラメータをｓ１、基準円Ｃ２の半径（曲率半径）を示すパラメータをｓ２とする。基準円Ｃ１が属する平面に対する単位法線ベクトルをＶ１とする。基準円Ｃ２が属する平面に対する単位法線ベクトルをＶ２とする。基準円Ｃ１の周上の点をＰ１、基準円Ｃ２の周上の点をＰ２とする。この場合、式（１４）より、以下の式（１５）が得られる。
Figure JPOXMLDOC01-appb-M000015
　演算部３２の候補点検出部３４は、基準円Ｃ１と基準円Ｃ２との交点（候補点）を検出する。この検出においては、以下で説明する判別式Ｄを用いて、交点を検出する。以下、判別式Ｄの導出について説明する。
　基準円Ｃ１と基準円Ｃ２との交点をＰｃとする。交点Ｐｃの位置ベクトルは、以下の式（１６）で定義することができる。式（１６）において、β、γ及びδは、任意の実数である。
Figure JPOXMLDOC01-appb-M000016
　交点Ｐｃは式（１５）のいずれも満たす必要がある。そこで、式（１６）を式（１５）の各式に代入して、以下の式（１７）が得られる。
Figure JPOXMLDOC01-appb-M000017
　式（１７）を、β及びγについて解くと、以下の式（１８）が得られる。
Figure JPOXMLDOC01-appb-M000018
　また、交点Ｐｃにおいては、以下の式（１９）が成立する。
Figure JPOXMLDOC01-appb-M000019
 　　式（１６）を用いて式（１９）を展開すると、以下の式（２０）が得られる。
Figure JPOXMLDOC01-appb-M000020
　式（２０）に式（１８）を代入し、δについて解くと、以下の式（２１）が得られる。
Figure JPOXMLDOC01-appb-M000021
 　式（２１）に示すＤは、交点有無の判別式であり、以下の式（２２）で表される。
Figure JPOXMLDOC01-appb-M000022
　式（１９）は、判別式Ｄの平方根を含む。そのため、交点Ｐｃを表す式（１４）の解は、判別式Ｄの値により場合分けが必要となる。
判別式Ｄが正の値をとる場合（Ｄ＞０）
　判別式Ｄが正の値をとる場合、δは絶対値が等しい正負の２値をとる。よって、交点Ｐｃを表す式（１６）の解は２つ得られる。すなわち、この場合には、基準円Ｃ１と基準円Ｃ２とは、真球ＣＢ上の２つの交点Ｐｃ１及びＰｃ２で交差する。図２８は、基準円Ｃ１と基準円Ｃ２とが２つの交点を有する（交接する）場合を示す図である。
　式（１６）に式（１８）及び式（２１）を代入することで、交点Ｐｃ１及びＰｃ２の位置ベクトルは、以下の式（２３）で表される。
Figure JPOXMLDOC01-appb-M000023
判別式Ｄが負の値をとる場合（Ｄ＜０）
　判別式Ｄが負の値をとる場合、δは虚数解となるので、基準円Ｃ１と基準円Ｃ２とは交点を有しない。基準円Ｃ１と基準円Ｃ２とが交点を有しない場合、基準円Ｃ１と基準円Ｃ２とは、分離又は内包の関係にある。図２９は、基準円Ｃ１と基準円Ｃ２とが分離の関係にある場合を示す図である。この場合、図２９に示すように、基準円Ｃ１と基準円Ｃ２とは空間的に離隔しており、交点を有しない。図３０は、基準円Ｃ１と基準円Ｃ２とが内包の関係にある場合を示す図である。この場合、図３０に示すように、基準円Ｃ１と基準円Ｃ２とは真球ＣＢ上で領域を共有するものの、基準円Ｃ１を構成する線分と基準円Ｃ２を構成する線分とは、交点を有しない。
判別式Ｄが０の場合（Ｄ＝０）
　判別式Ｄが０の場合、δも０となる。この場合、基準円Ｃ１と基準円Ｃ２とは接している状態にある。基準円Ｃ１と基準円Ｃ２とが接している状態は、２つに分けて考えることができる。１つは、基準円Ｃ１と基準円Ｃ２とが、交点Ｐｃを接点として外接又は内接する場合である。もう１つは、基準円Ｃ１と基準円Ｃ２とが一致する場合である。
基準円Ｃ１と基準円Ｃ２とが外接又は内接する場合
　判別式Ｄが０で、かつ、以下の式（２４）を満たす場合には、基準円Ｃ１と基準円Ｃ２とは、１つの交点を有する。
Figure JPOXMLDOC01-appb-M000024
 　この場合の基準円Ｃ１と基準円Ｃ２との交点Ｐｃ０の位置ベクトルは、式（１６）に式（１８）及び式（２１）を代入することで、以下の式（２５）で表される。
Figure JPOXMLDOC01-appb-M000025
 　図３１は、基準円Ｃ１と基準円Ｃ２とが外接の関係にある場合を示す図である。この例では、基準円Ｃ１と基準円Ｃ２とは、交点Ｐｃ０で外接する。図３２は、基準円Ｃ１と基準円Ｃ２とが内接の関係にある場合を示す図である。この例では、基準円Ｃ１は交点Ｐｃ０で基準円Ｃ２と内接する。
基準円Ｃ１と基準円Ｃ２とが一致する場合
　また、判別式Ｄが０で、かつ、以下の式（２６）を満たす場合には、基準円Ｃ１と基準円Ｃ２とは一致する。
Figure JPOXMLDOC01-appb-M000026
 　図３３は、基準円Ｃ１と基準円Ｃ２とが一致する場合を示す図である。この例では、基準円Ｃ２は、基準円Ｃ１と同一の円となる。この場合、基準円Ｃ１及び基準円Ｃ２の全周の任意の位置に交点が存在する。この場合、２本の線分のそれぞれの始点及び終点を交点とする。
　図３４は、基準円が一致し、かつ、２本の線分が分離している場合を示す図である。この例では、線分Ｌ１の始点ＰＳ１、線分Ｌ１の終点ＰＥ１、線分Ｌ２の始点ＰＳ２、線分Ｌ２の終点ＰＥ２の４点が交点Ｐｃとなる。
　図３５は、基準円が一致し、かつ、一方の線分の始点と他方の線分の終点とが重なる場合を示す図である。この例では、線分Ｌ１の始点ＰＳ１であるとともに線分Ｌ２の終点ＰＥ２でもある点、線分Ｌ１の終点ＰＥ１、線分Ｌ２の始点ＰＳ２の３点が交点Ｐｃとなる。
　図３６は、基準円が一致し、かつ、２本の線分間に１つの重複部分がある場合を示す図である。この例では、線分Ｌ１の始点ＰＳ１、線分Ｌ１の終点ＰＥ１、線分Ｌ２の始点ＰＳ２、線分Ｌ２の終点ＰＥ２の４点が交点Ｐｃとなる。
　図３７は、基準円が一致し、かつ、一方の線分の始点と他方の線分の終点とが重なるとともに２本の線分間に１つの重複部分が有る場合を示す図である。この例では、線分Ｌ１の始点ＰＳ１であるとともに線分Ｌ２の終点ＰＥ２でもある点、線分Ｌ１の終点ＰＥ１、線分Ｌ２の始点ＰＳ２の３点が交点Ｐｃとなる。
　図３８は、基準円が一致し、かつ、２本の線分間に２つの重複部分が有る場合を示す図である。この例では、線分Ｌ１の始点ＰＳ１、線分Ｌ１の終点ＰＥ１、線分Ｌ２の始点ＰＳ２、線分Ｌ２の終点ＰＥ２の４点が交点Ｐｃとなる。
　以上、２つの基準円が交点を持つか否か、或いは２つの基準円が一致するか否かについて説明したが、２本の線分が交点を持つか否かは、基準円上における線分の区間を考慮しなければならない。つまり、基準円Ｃ１と基準円Ｃ２との交点が線分Ｌ１及び線分Ｌ２の区間内に存在しない場合には、線分Ｌ１と線分Ｌ２とは交点を有しない。
　従って、本実施の形態では、基準円Ｃ１と基準円Ｃ２との交点は、必ずしも線分Ｌ１と線分Ｌ２との交点とはならない。よって、基準円Ｃ１と基準円Ｃ２との交点と、線分Ｌ１と線分Ｌ２との交点とを区別するため、上述で検出した基準円Ｃ１と基準円Ｃ２との交点を、候補点と称する。
　以下、検出部３５が、基準円Ｃ１上の線分Ｌ１が式（１４）で表される候補点Ｐｃを含むか否かを判定する方法について説明する。判定に当たり、線分Ｌ１の中心角Ψにより、場合分けを行う。
中心角Ψがπ以上２π以下である場合（π≦Ψ≦２π）
　図３９は、中心角Ψが２πである場合（Ψ＝２π）の線分Ｌ１を示す図である。中心角Ψが２πの場合、候補点Ｐｃは線分Ｌ１上に存在する。また、図４０は、中心角Ψがπ以上かつ２πよりも小さい場合（π≦Ψ＜２π）の線分Ｌ１を示す図である。この場合、線分Ｌ１は、半円弧又は優弧となり、以下の式（２７）を満たす。
Figure JPOXMLDOC01-appb-M000027
 　候補点Ｐｃは、以下の式（２８）又は式（２９）を満たす場合、線分Ｌ１上に存在する。
Figure JPOXMLDOC01-appb-M000028
Figure JPOXMLDOC01-appb-M000029
中心角Ψがπよりも小さい場合（０＜Ψ＜π）
　図４１は、中心角Ψがπより小さい場合（０＜Ψ＜π）の線分Ｌ１を示す図である。この場合、円弧は劣弧となり、以下の式（３０）を満たす。
Figure JPOXMLDOC01-appb-M000030
 　候補点Ｐｃは、上述の式（２８）及び式（２９）をともに満たす場合、線分Ｌ１上に存在する。
　以上、線分Ｌ１が交点を有するかを判定する方法について説明したが、線分Ｌ２についても同様に交点を有するか否かを判定することができる。
　以上より、線分Ｌ１及び線分Ｌ２が同じ候補点Ｐｃを含む場合、これが交点Ｐｃであると判定できる。この場合、線分Ｌ１と線分Ｌ２とは、２点で交わり（この場合を交接という）、接し、又は一致していると判定することができる。
　以下、上記の交点検出（図２３のステップＳ２３）の手順について整理する。図４２は、交点検出部２２での線分の交点検出動作を示すフローチャートである。
ステップＳＳ１
　候補点検出部３４は、判別式Ｄを算出する。
ステップＳＳ２
　候補点検出部３４は、判別式Ｄが０よりも小さいか否かを判定する。これにより、候補点が存在するかを判定できる。判別式Ｄが０よりも小さい場合、候補点は存在しない。判別式Ｄが０以上の場合、少なくとも１つ以上の候補点が存在する。
ステップＳＳ３
　判別式Ｄが０以上の場合、検出部３５は、判別式Ｄが０であるかを判定する。
ステップＳＳ４
　判別式Ｄが０よりも大きい場合、検出部３５は、候補点Ｐｃ１を算出する。
ステップＳＳ５
　検出部３５は、候補点Ｐｃ１について、交点判定処理を行う。交点判定処理については後述する。
ステップＳＳ６
　検出部３５は、候補点Ｐｃ２を算出する。
ステップＳＳ７
　検出部３５は、候補点Ｐｃ２について交点判定処理を行う。交点判定処理については後述する。
ステップＳＳ８
　判別式Ｄが０の場合の場合、検出部３５は、式（３１）を満たすか判定する。
Figure JPOXMLDOC01-appb-M000031
ステップＳＳ９
　式（３１）を満たす場合、検出部３５は、候補点Ｐｃ０を算出する。
ステップＳＳ１０
　検出部３５は、候補点Ｐｃ０について交点判定処理を行う。交点判定処理については後述する。
ステップＳＳ１１
　式（３１）を満たさない場合、検出部３５は、線分Ｌ１の始点ＰＳ１について、交点判定処理を行う。
ステップＳＳ１２
　検出部３５は、線分Ｌ１の終点ＰＥ１について、交点判定処理を行う。
ステップＳＳ１３
　検出部３５は、線分Ｌ２の始点ＰＳ２について、交点判定処理を行う。
ステップＳＳ１４
　検出部３５は、線分Ｌ２の終点ＰＥ２について、交点判定処理を行う。
　次いで、交点判定処理について説明する。図４３は、交点判定処理を示すフローチャートである。
ステップＳＲ１
　判定対象点ＰＪとして、直前のステップで算出した候補点を設定する。
ステップＳＲ２
　判定対象点ＰＪが、線分Ｌ１上に存在するかを判定する範囲検証処理を行う。範囲検証処理の詳細は、後述する。判定対象点ＰＪが線分Ｌ１上に存在しない場合には、処理を終了する。
ステップＳＲ３
　判定対象点ＰＪが線分Ｌ１上に存在する場合、判定対象点ＰＪが線分Ｌ２上に存在するかを判定する範囲検証処理を行う。範囲検証処理の詳細は、後述する。判定対象点ＰＪが線分Ｌ２上に存在しない場合には、処理を終了する。
ステップＳＲ４
　判定対象点ＰＪが線分Ｌ１上及びＬ２上に存在する場合、判定対象点ＰＪを候補点として登録する。
　上述のステップＳＲ２及びＳＲ３における範囲検証処理について説明する。図４４は、範囲検証処理を示すフローチャートである。ここでは、検証の対象となる線分をＬＪと称する。
ステップＳＡ１
　判定対象線分ＬＪが円であるかを判定する。
ステップＳＡ２
　判定対象線分ＬＪが円ではない場合、線分が優弧であるかを判定する。
ステップＳＡ３
　判定対象線分ＬＪが優弧または半円弧である場合、式（２８）及び式（２９）の少なくともいずれか一方を満たすか判定する。式（２８）及び式（２９）の少なくともいずれか一方を満たす場合、判定対象点ＰＪは、判定対象線分ＬＪ上にある（ＹＥＳ判定）。式（２８）及び式（２９）のいずれも満さない場合、判定対象点ＰＪは、判定対象線分ＬＪ上に存在しない（ＮＯ判定）。
ステップＳＡ４
　判定対象線分ＬＪが劣弧である場合、式（２８）及び式（２９）の両方を満たすか判定する。式（２８）及び式（２９）の両方を満たす場合、判定対象点ＰＪは、判定対象線分ＬＪ上にある（ＹＥＳ判定）。式（２８）及び式（２９）の少なくとも一方を満たさない場合、判定対象点ＰＪは、判定対象線分ＬＪ上に存在しない（ＮＯ判定）。
　以上より、本実施の形態によれば、真球上に設定された２本の線分が交点を有するか否かを確実に判定することができる。これにより、真球上の円弧で表される２本の空路が交差するか否か、空域を構成する線分が交差するか否かを、確実に判定することが可能である。
　上述では、一般的な２本の線分が交点を有するか否かについて説明したが、交点検出部２２は、上記の２本の線分が交点検出を、判定対象閉曲線ＡＺ１を囲む閉曲線を構成する線分のそれぞれと、転写像ＡＺ２を囲む閉曲線を構成する線分のそれぞれと、に適用することで、判定対象閉曲線ＡＺ１と転写像ＡＺ２とが交点を有するか否かを、具体的かつ容易に検出可能であることが理解できる。
　なお、本発明は上記実施の形態に限られたものではなく、趣旨を逸脱しない範囲で適宜変更することが可能である。
　以上では、空域情報処理装置及びこの装置で行われる空域情報処理方法について説明した。しかし、本発明は、これに限定されるものではない。本発明は、任意の処理を、ＣＰＵ（Central Processing Unit）にコンピュータプログラムを実行させることにより実現することも可能である。
　プログラムは、様々なタイプの非一時的なコンピュータ可読媒体（non-transitory computer readable medium）を用いて格納され、コンピュータに供給することができる。非一時的なコンピュータ可読媒体は、様々なタイプの実体のある記録媒体（tangible storage medium）を含む。非一時的なコンピュータ可読媒体の例は、磁気記録媒体（例えばフレキシブルディスク、磁気テープ、ハードディスクドライブ）、光磁気記録媒体（例えば光磁気ディスク）、ＣＤ－ＲＯＭ（Read Only Memory）、ＣＤ－Ｒ、ＣＤ－Ｒ／Ｗ、半導体メモリ（例えば、マスクＲＯＭ、ＰＲＯＭ（Programmable ROM）、ＥＰＲＯＭ（Erasable PROM）、フラッシュＲＯＭ、ＲＡＭ（random access memory））を含む。また、プログラムは、様々なタイプの一時的なコンピュータ可読媒体（transitory computer readable medium）によってコンピュータに供給されてもよい。一時的なコンピュータ可読媒体の例は、電気信号、光信号、及び電磁波を含む。一時的なコンピュータ可読媒体は、電線及び光ファイバ等の有線通信路、又は無線通信路を介して、プログラムをコンピュータに供給できる。
　以上、実施の形態を参照して本願発明を説明したが、本願発明は上記によって限定されるものではない。本願発明の構成や詳細には、発明のスコープ内で当業者が理解し得る様々な変更をすることができる。
１　閉曲線読込部
２　転写部
３　線分生成部
４　空域認識部
５　記憶部
２１　転写処理部
２２　交点検出部
２２　交点検出部
３１　記憶装置
３２　演算部
３３　バス
３４　候補点検出部
３５　検出部
１００　空域情報処理装置
ＡＺ１　判定対象閉曲線
ＡＺ２　転写像
Ｃ、Ｃ１、Ｃ２、基準円
ＣＢ　真球
ＣＣ１　円
ＣＣ２、ＣＣ３　円弧
Ｄ１　基本形状データベース
Ｄ２　空域情報データベース
Ｌｄ　判定用線分
ＬＪ　判定対象線分
Ｌｐ　仮の線分
Ｏ　真球の中心
ＰＡ、Ｐ０　点

Claims (12)

  1. 　球面上の１又は複数の線分からなり、空域の輪郭を現す閉曲線の全部又は一部を、前記閉曲線と交点を有しないように前記球面上の別の位置に転写した転写像を生成する転写手段と、
    　前記閉曲線を構成する線分から、前記閉曲線を構成する他の線分と交点を有することなく、かつ、前記転写像と交点を有する判定用線分を生成する線分生成手段と、
    　前記閉曲線によって区切られる前記球面上の２つの領域のうち、前記判定用線分が存在する領域を、前記空域として認識する空域認識手段と、を備える、
    　空域情報処理装置。
  2. 　前記転写手段は、前記閉曲線を、球の中心に対して点対称の位置に転写して前記転写像を生成する、
    　請求項１に記載の空域情報処理装置。
  3. 　前記転写手段は、前記閉曲線を、球の中心を通る回転軸まわりに所定の角度だけ回転させた位置に転写して前記転写像を生成する、
    　請求項１に記載の空域情報処理装置。
  4. 　前記転写手段は、
    　前記閉曲線を、球の中心を通る第１回転軸まわりに所定の角度だけ回転させた位置に転写して第１転写像を生成し、
    　前記第１転写像が前記閉曲線と交点を有しない場合、前記第１転写像を前記転写像として設定し、
    　前記第１転写像が前記閉曲線と交点を有する場合、前記閉曲線を、球の中心を通り、かつ、前記第１回転軸に対して垂直な第２回転軸まわりに所定の角度だけ回転させた位置に転写して第２転写像を生成し、
    　前記第２転写像が前記閉曲線と交点を有しない場合、前記第２転写像を前記転写像として設定し、
    　前記第２転写像が前記閉曲線と交点を有する場合、前記閉曲線を、球の中心を通り、かつ、前記第１回転軸及び第２回転軸に対して垂直な第３回転軸まわりに所定の角度だけ回転させた位置に転写して第３転写像を生成し、
    　前記第３転写像が前記閉曲線と交点を有しない場合、前記第３転写像を前記転写像として設定する、
    　請求項３に記載の空域情報処理装置。
  5. 　前記球は地球に対応し、前記第１回転軸は地軸である、
    　請求項４に記載の空域情報処理装置。
  6. 　前記回転軸は、
    　前記閉曲線上の複数の座標の平均座標と前記球の中心とを通る線に対して垂直である、
    　請求項３に記載の空域情報処理装置。
  7. 　前記回転軸は、
    　前記閉曲線上の複数の点の緯度及び経度の平均で表される座標と前記球の中心とを通る線に対して垂直である、
    　請求項３に記載の空域情報処理装置。
  8. 　前記回転軸は、
    　前記閉曲線上の１つの点と前記球の中心とを通る線に対して垂直である、
    　請求項３に記載の空域情報処理装置。
  9. 　前記回転軸は、
    　前記閉曲線上の２つの点を結ぶ線の中点と前記球の中心とを通る線に対して垂直である、
    　請求項３に記載の空域情報処理装置。
  10. 　前記線分生成手段は、
    　前記閉曲線を構成する複数の線分のいずれかに、第１の点を設定し、
    　前記第１の点と前記転写像を構成する線分上の第２の点とを結ぶ第１線分を生成し、
    　前記第１線分と前記閉曲線を構成する前記複数の線分との交点を全て検出し、
    　前記第１線分のうち、検出した前記交点のうちで最も前記第２の点に近いものと前記第２の点との間の区間を、前記判定用線分として設定する、
    　請求項１乃至９のいずれか一項に記載の空域情報処理装置。
  11. 　転写手段に、球面上の１又は複数の線分からなり、空域の輪郭を現す閉曲線の全部又は一部を、前記閉曲線と交点を有しないように前記球面上の別の位置に転写した転写像を生成させ、
    　線分生成手段に、前記閉曲線を構成する線分から、前記閉曲線を構成する他の線分と交点を有することなく、かつ、前記転写像と交点を有する判定用線分を生成させ、
    　空域認識手段に、前記閉曲線によって区切られる前記球面上の２つの領域のうち、前記判定用線分が存在する領域を、前記空域として認識させる。
    　空域情報処理方法。
  12. 　球面上の１又は複数の線分からなり、空域の輪郭を現す閉曲線の全部又は一部を、前記閉曲線と交点を有しないように前記球面上の別の位置に転写した転写像を生成する処理と、
    　線分生成手段に、前記閉曲線を構成する線分から、前記閉曲線を構成する他の線分と交点を有することなく、かつ、前記転写像と交点を有する判定用線分を生成する処理と、
    　空域認識手段に、前記閉曲線によって区切られる前記球面上の２つの領域のうち、前記判定用線分が存在する領域を、前記空域として認識する処理と、をコンピュータに実行させる、
    　空域情報処理プログラムが格納された非一時的なコンピュータ可読媒体。
PCT/JP2014/003783 2014-07-17 2014-07-17 空域情報処理装置、空域情報処理方法、空域情報処理プログラムが格納された非一時的なコンピュータ可読媒体 WO2016009465A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2016533997A JP6233517B2 (ja) 2014-07-17 2014-07-17 空域情報処理装置、空域情報処理方法、空域情報処理プログラム
PCT/JP2014/003783 WO2016009465A1 (ja) 2014-07-17 2014-07-17 空域情報処理装置、空域情報処理方法、空域情報処理プログラムが格納された非一時的なコンピュータ可読媒体
US15/326,695 US20170206663A1 (en) 2014-07-17 2014-07-17 Airspace information processing device, airspace information processing method, and non-transitory computer-readable medium storing airspace information processing program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2014/003783 WO2016009465A1 (ja) 2014-07-17 2014-07-17 空域情報処理装置、空域情報処理方法、空域情報処理プログラムが格納された非一時的なコンピュータ可読媒体

Publications (1)

Publication Number Publication Date
WO2016009465A1 true WO2016009465A1 (ja) 2016-01-21

Family

ID=55077991

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2014/003783 WO2016009465A1 (ja) 2014-07-17 2014-07-17 空域情報処理装置、空域情報処理方法、空域情報処理プログラムが格納された非一時的なコンピュータ可読媒体

Country Status (3)

Country Link
US (1) US20170206663A1 (ja)
JP (1) JP6233517B2 (ja)
WO (1) WO2016009465A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105206114B (zh) * 2015-08-10 2018-08-21 华为技术有限公司 飞行控制、许可、安全维护方法和装置、服务器、飞行器
CN109559567B (zh) * 2018-12-12 2020-11-24 四川九洲空管科技有限责任公司 一种双层网格化空域资源动态状态管理方法及系统
CN114543816B (zh) * 2022-04-25 2022-07-12 深圳市赛特标识牌设计制作有限公司 一种基于物联网的导视方法、装置以及系统
CN115330989B (zh) * 2022-10-13 2023-01-20 南通红运金属科技有限公司 一种用于激光束切割的控制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003085569A (ja) * 2001-09-14 2003-03-20 Canon Inc 内外点判定アルゴリズム
JP2009133798A (ja) * 2007-12-03 2009-06-18 Furuno Electric Co Ltd 位置検出方法および位置検出装置
JP2012088902A (ja) * 2010-10-19 2012-05-10 Fujitsu Advanced Engineering Ltd 位置関係判定装置、位置関係判定方法、および位置関係判定プログラム
WO2014108957A1 (ja) * 2013-01-08 2014-07-17 日本電気株式会社 座標変換装置、座標変換プログラムが格納された非一時的なコンピュータ可読媒体、座標変換方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2549641A (en) * 1947-10-31 1951-04-17 Rca Corp Display of three-dimensional information
US7421340B2 (en) * 2005-02-28 2008-09-02 Vectronix Ag Method, apparatus and computer program for azimuth determination e.g. for autonomous navigation applications
FR2925712B1 (fr) * 2007-12-21 2010-01-01 Thales Sa Procede pour l'aide a l'atterrissage d'aeronef utilisant un gps et un mls dans le cadre d'une approche axiale calculee.
US9879993B2 (en) * 2010-12-23 2018-01-30 Trimble Inc. Enhanced bundle adjustment techniques
US8862395B2 (en) * 2011-01-31 2014-10-14 Raytheon Company Coded marker navigation system and method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003085569A (ja) * 2001-09-14 2003-03-20 Canon Inc 内外点判定アルゴリズム
JP2009133798A (ja) * 2007-12-03 2009-06-18 Furuno Electric Co Ltd 位置検出方法および位置検出装置
JP2012088902A (ja) * 2010-10-19 2012-05-10 Fujitsu Advanced Engineering Ltd 位置関係判定装置、位置関係判定方法、および位置関係判定プログラム
WO2014108957A1 (ja) * 2013-01-08 2014-07-17 日本電気株式会社 座標変換装置、座標変換プログラムが格納された非一時的なコンピュータ可読媒体、座標変換方法

Also Published As

Publication number Publication date
JPWO2016009465A1 (ja) 2017-07-13
JP6233517B2 (ja) 2017-11-22
US20170206663A1 (en) 2017-07-20

Similar Documents

Publication Publication Date Title
EP3345164B1 (en) Mesh generation system and method
JP6233517B2 (ja) 空域情報処理装置、空域情報処理方法、空域情報処理プログラム
EP3828720A2 (en) Method and apparatus for merging data of building blocks, device and storage medium
JP6132033B2 (ja) 形状判定装置、形状判定プログラム、形状判定方法
US8675993B2 (en) Methods and apparatus for patch-based shape from shading
WO2019119187A1 (zh) 一种无人机的航线规划方法、控制设备及存储介质
WO2014108957A1 (ja) 座標変換装置、座標変換プログラムが格納された非一時的なコンピュータ可読媒体、座標変換方法
EP4024369A1 (en) Flight leg termination visualization systems and methods for flight leg termination visualization
US10108172B2 (en) Spiral toolpaths for high-speed machining of polygonal pockets
US11308700B2 (en) Augmented reality visualization of underground pipelines using geospatial databases and KM markers
JP2009133798A (ja) 位置検出方法および位置検出装置
US10146888B2 (en) Systems and methods for criteria analysis prototyping
KR20210107558A (ko) 라벨 생성 방법
JP2006334684A (ja) 面沿い加工工具経路決定方法、面沿い加工工具経路決定プログラム、面沿い加工工具経路決定装置
Mitchell Deterministic Linear Time for Maximal Poisson‐Disk Sampling using Chocks without Rejection or Approximation
JP6260688B2 (ja) 空域情報処理装置、空域情報処理方法、空域情報処理プログラム
US11074749B2 (en) Planar surface detection
US20160035131A1 (en) 3d-cg processing device and 3d-cg processing method
JP2007121242A (ja) タイヤの形状の比較方法及びタイヤの形状の比較方法用コンピュータプログラム
JP2007293391A (ja) 設計支援装置、設計支援方法、およびプログラム
US20140316751A1 (en) System and method of generating geometry of a swept volume of a spun tool
Giang Method of Orienting Curves for Some 2D and 3D Geodesic Problems in Computational Geometry
JP2020035087A (ja) 地理情報処理装置、曲線近似方法及び曲線近似プログラム
KR101434380B1 (ko) 그래픽 가속기 기반의 3차원 폴리곤 모델에서의 특징선 추출 방법 및 장치
JPH11202919A (ja) 板金ｃａｄ図面からの立体図作成方法及びそのプログラムを記憶した記憶媒体

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 14897434

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2016533997

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 15326695

Country of ref document: US

122 Ep: pct application non-entry in european phase

Ref document number: 14897434

Country of ref document: EP

Kind code of ref document: A1