The present invention relates to a coordinate transformation device, a non-transitory computer readable medium storing a coordinate transformation program, and a coordinate transformation method.
Today, various navigation systems have been put into practice to monitor vehicles on earth. In order to manage the operation of aircraft whose travel distance is longer than other carriers, it is necessary to calculate the azimuth and distance of the aircraft in a wide area. Aircraft navigation systems are generally required to process large-scale spatial information accurately and effectively in a wide area such as a country's territory and air space, or a flight information region (FIR).
As an example of such navigation systems, a method for calculating a distance using piece-wise linear interpolation has been proposed (Patent Literature 1). In this method, the distance between two points is divided into segments of a certain length, and the distance of each segment is calculated. Further, a correction is made to the connecting part between the segments, which makes it possible to calculate the distance between two points. Accordingly, the calculation can be simplified and performed efficiently.
As another example of the navigation systems, a distance calculation method in which the shape of the earth (spheroid) is taken into consideration at a radio station tower located on land has been proposed (Patent Literature 2). In this method, a highly accurate distance calculation can be achieved within a line-of-sight distance from the radio station tower by performing a correction calculation in consideration of the spheroid shape.
However, the inventor has found that the above-mentioned methods have the following problems. In the method disclosed in Patent Literature 1, a highly accurate distance calculation is achieved by using piece-wise linear interpolation. However, if two points to be compared are far away from each other, the number of segments between the two points increases. As a result, there is a problem that the processing efficiency deteriorates and calculation errors accumulate, leading to a large error in the calculation. In other words, this method is not suitable in principle for calculation of a distance between points which are far away from each other. Moreover, the method is not suitable for calculation of a shortest path.
Furthermore, the method disclosed in Patent Literature 2 has a problem that the method cannot be applied if the two points to be compared are far away and over the horizon. That is, the method cannot be applied to the calculation of a distance of, for example, an aircraft, over a wide area.
The present invention has been made in view of the above-mentioned circumstances, and an object of the present invention is to provide a coordinate transformation device, a non-transitory computer readable medium storing a coordinate transformation program, and a coordinate transformation method which are capable of calculating coordinates on a spheroid at a high speed and with a minimum error.
A coordinate transformation device according to an exemplary aspect of the present invention includes: storage means for storing information specifying a shape of a spheroid; information indicating coordinates on the spheroid described using the information specifying the shape of the spheroid; and information indicating a transformation parameter for transforming the information indicating coordinates on the spheroid into information indicating coordinates on a true sphere; and operation means for generating a formula for performing coordinate transformation processing by substituting the transformation parameter into a predetermined formula, transforming the information indicating coordinates on the spheroid into coordinates on the true sphere by substituting, into the generated formula, the information specifying the shape of the spheroid and the information indicating coordinates on the spheroid, and outputting the transformed information indicating coordinates on the true sphere.
A non-transitory computer readable medium storing a coordinate transformation program according to an exemplary aspect of the present invention causes a computer to execute processing including: reading information specifying a shape of a spheroid, information indicating coordinates on the spheroid described using the information specifying the shape of the spheroid, and information indicating a transformation parameter used for coordinate transformation processing; generating a formula for performing the coordinate transformation processing by substituting the transformation parameter into a predetermined formula; transforming the information indicating coordinates on the spheroid into information indicating coordinates on a true sphere by substituting, into the generated formula, the information specifying the shape of the spheroid and the information indicating coordinates on the spheroid; and outputting the transformed information indicating coordinates on the true sphere.
A coordinate transformation method according to an exemplary aspect of the present invention includes: generating a formula for performing coordinate transformation processing by substituting a transformation parameter into a predetermined formula, the transformation parameter being used for the coordinate transformation processing; substituting, into the generated formula, information specifying a shape of a spheroid and information indicating coordinates on the spheroid described using the information specifying the shape of the spheroid, and transforming the information indicating coordinates on the spheroid into information indicating coordinates on a true sphere; and outputting the transformed information indicating coordinates on the true sphere.
According to the present invention, it is possible to provide a coordinate transformation device, a non-transitory computer readable medium storing a coordinate transformation program, and a coordinate transformation method which are capable of calculating coordinates on a spheroid at a high speed and with a minimum error.
Exemplary embodiments of the present invention will be described below with reference to the drawings. In the drawings, the same elements are denoted by the same reference numerals, and a repeated explanation is omitted as necessary.
First, a geographical information management device 100 according to a first exemplary embodiment of the present invention will be described. The geographical information management device 100 is configured using a hardware resource such as a computer system.
In general, calculations for control and navigation of vehicles that travel above the earth, such as aircraft, require arithmetic processing by digitizing map information. However, the shape of the earth is not a true sphere but rather a spheroid which is compressed in the north-south direction and has a radius that is maximum in the vicinity of the equator. Accordingly, if the arithmetic processing is to be performed using coordinate values on the spheroid, a vast amount of complicated processing is required.
The geographical information management device 100 performs coordinate transformation processing for projecting coordinates on a spheroid as coordinates on a true sphere. The arithmetic processing is performed using the projected coordinates on the true sphere, thereby simplifying the processing, speeding up the processing, and achieving the calculations for control and navigation of vehicles traveling above the earth, with small hardware resources. That is, the geographical information management device 100 is an example of a coordinate transformation device that is configured so as to be able to execute coordinate transformation processing.
The input device 1 is used to input data to the geographical information management device 100 from the outside. As the input device 1, various data input means, such as a keyboard, a mouse, a DVD (Digital Versatile Disc) drive, and network connection, can be used.
The storage device 2 can store a database in which data provided through the input device 1 is stored, and can also store a program provided for processing in the operation unit 3. As the storage device 2, various storage devices, such as a hard disk drive and a flash memory, can be used. Specifically, the storage device 2 stores a basic shape database D1, an airspace information database D2, a parameter information database D3, and a projection information database D4.
The basic shape database D1 is preliminarily given specific information. The basic shape database D1 will be described in detail later.
The airspace information database D2 is, for example, information input through the input device 1, and includes coordinates of an aircraft to be monitored on the spheroid, and coordinate information about an airspace on the spheroid. The airspace information database D2 will be described in detail later.
The parameter information database D3 includes parameters for transforming coordinate information on the spheroid in the airspace information database D2 into coordinate information on the true sphere. In other words, the coordinates on the spheroid are projected onto the true sphere. The transformed coordinates on the true sphere are hereinafter referred to as projected coordinates.
The projection information database D4 includes projected coordinate information on the true sphere. The projection information database D4 includes, for example, a projection space database D41, a projection line segment database D42, and a projection point database D43. The projection space database D41 includes coordinate information specifying a predetermined three-dimensional space on the true sphere. The projection line segment database D42 includes coordinate information specifying a line segment that connects two points on the true sphere. The projection point database D43 includes coordinate information indicating a predetermined point on the true sphere.
The storage device 2 stores a coordinate transformation program PRG1 for specifying arithmetic processing for coordinate transformation.
The operation unit 3 can read out the program and database from the storage device 2, and can perform necessary arithmetic processing. The operation unit 3 consists of, for example, a CPU (Central Processing Unit).
The display device 4 displays the coordinates, flight information, and the like of an aircraft in a visible manner according to the operation result of the operation unit 3. As the display device 4, various display devices such as a liquid crystal display monitor can be used.
Next, an operation of the geographical information management device 100 will be described.
First, the operation unit 3 reads out the coordinate transformation program PRG1. The coordinate transformation program PRG1 is a program for transforming coordinates on the spheroid into projected coordinates on the true sphere by using the basic shape database D1, the airspace information database D2, and the parameter information database D3. The coordinate transformation program PRG1 is read out from, for example, the storage device 2.
Next, the operation unit 3 reads out the basic shape database D1, the airspace information database D2, and the parameter information database D3 from the storage device 2.
The operation unit 3 generates a formula for performing coordinate transformation by substituting the information included in the basic shape database D1 and the parameter information database D3 into a formula specified by the coordinate transformation program PRG1.
The operation unit 3 transforms coordinates on the spheroid into projected coordinates on the true sphere by substituting the coordinate information included in the airspace information database D2 into the generated formula. That is, the operation unit 3 transforms the airspace information database D2 into the projection information database D4. The coordinate transformation in Step S14 will be described in detail later.
The operation unit 3 outputs to the outside the projected coordinates on the true sphere that are included in the projection information database D4. For example, the operation unit 3 outputs the projection information database D4 to the storage device 2.
After that, the operation unit 3 performs arithmetic processing necessary for control or navigation calculation of a vehicle, by using the projection information database D4. In this case, for example, the operation unit 3 reads out the program for specifying arithmetic processing from the storage device 2 and performs the arithmetic processing. For example, the operation unit 3 can calculate the azimuth and distance between points by using the projection information database D4, and can obtain an intersection between different line segments. For example, in the case of calculating the azimuth and distance between points, the operation unit 3 reads out an azimuth/distance calculation program from the storage device 2. In the case of obtaining an intersection between different line segments, the operation unit 3 reads out an intersection calculation program from the storage device 2. In the case of displaying these calculation results, the operation unit 3 reads out an image transformation program and processes the calculation results into a form that can be displayed. The information processed into a form that can be displayed is output to the display device 4.
Next, a method for representing coordinates on the spheroid that are included in the airspace information database D2 will be described.
Coordinates on the spheroid are represented by using, for example, a WGS84 coordinate system.
where f represents flattening. If the spheroid EB represents the earth, f=1/298.257223563 holds. The symbol “e” represents eccentricity.
The basic shape database D1 will now be described in detail. The basic shape database D1 is preliminarily given specific information.
The airspace information database D2 will be described in detail.
The projection information database D4 will be described in detail.
Next, properties of latitude line spacing and longitude line spacing on the spheroid EB will be considered. In general, the latitude line spacing on the true sphere is constant. On the other hand, on the spheroid EB, at the altitude 0 (ground surface), the latitude line spacing is minimum in the vicinity of the equator and is maximum in the vicinity of the north pole or the south pole. A latitude line spacing Dθ at the ground surface on the spheroid EB in the WGS84 coordinate system is represented by the following Formula (2) based on Formula (1).
Next, the longitude line spacing will be considered. The longitude line spacing on the true sphere is in proportion to the cosine of the latitude. At the altitude 0 (ground surface), the longitude line spacing takes a minimum value “a” at the equator and takes a minimum value “0” at the north pole or the south pole. The same holds true for the spheroid EB. A longitude line spacing Dφ at the ground surface on the true sphere and on the spheroid EB in the WGS84 coordinate system is represented by the following Formula (3) based on Formula (1).
A method for representing projected coordinates on the true sphere included in the projection information database D4 will be described.
X=(R+h)cos Θ cos Φ
Y=(R+h)cos Θ sin Φ
Z=(R+h)sin Θ (4)
When Formula (4) is compared with Formula (1), it can be understood that the projected coordinates on the true sphere CB can be expressed in the form of a simpler function than that of the coordinates on the spheroid EB.
A latitude line spacing Δθ at the ground surface on the true sphere is represented by the following Formula (5) based on Formula (4).
A longitude line spacing Δφ at the ground surface on the true sphere is represented by the following Formula (6) based on Formula (4).
In this exemplary embodiment, coordinates on the spheroid EB in the WGS84 coordinate system described above are transformed into projected coordinates that are projected on the true sphere CB. Transformation parameters for transforming coordinates on the spheroid EB in the WGS84 coordinate system described above into projected coordinates on the true sphere CB will be described below.
In order to project the spheroid EB onto the true sphere CB with as small a distortion as possible, it is necessary that the latitude line spacing Dθ on the spheroid EB and the latitude line spacing Δθ on the true sphere CB be equal (Δθ=Dθ) at the altitude h=0 and at a reference latitude θ0. Accordingly, the following Formula (7) is obtained by Formulas (2) and (5).
In order to project the spheroid EB onto the true sphere CB with as small a distortion as possible, it is necessary that the longitude line spacing Dφ on the spheroid EB and the longitude line spacing Δφ on the true sphere CB be equal (Δφ=Dφ) at the altitude h=0 and at the reference latitude θ0. Accordingly, the following Formula (8) is obtained by Formulas (3) and (6).
It is necessary that a primary change rate in the latitude direction of the latitude line spacing on the spheroid EB and a primary change rate in the latitude direction of the latitude line spacing on the true sphere CB be equal at the reference latitude θ0. Accordingly, the following Formula (9) is obtained by Formulas (2) and (5).
It is necessary that a primary change rate in the latitude direction of the longitude line spacing on the spheroid EB and a primary change range in the latitude direction of the longitude line spacing on the true sphere CB be equal at the reference latitude θ0. Accordingly, the following Formula (10) is obtained by Formulas (3) and (6).
Assume that a secondary change rate in the latitude direction of the longitude line spacing on the spheroid EB and a secondary change rate in the latitude direction of the longitude line spacing on the true sphere CB are equal at the reference latitude θ0. Accordingly, the following Formula (11) is obtained by Formulas (3) and (6).
Transformation parameters are calculated as follows by using the above formulas. First, a transformation parameter Pr for transforming the equatorial radius a of the spheroid EB in the WGS84 coordinate system into the radius R of the true sphere CB is calculated. The following Formula (12) is first obtained from Formulas (7), (9), and (10).
Formula (12)=Formula (10)×Formula (9)/Formula (7), that is:
Formula (12) is equal to the first term on the left side of Formula (11).
The following Formula (13) is obtained from Formulas (1) and (2).
Formula (13)=Formula (2)×{Formula (1)}2/R2, that is:
Formula (13) is equal to the second term on the left side of Formula (11).
Accordingly, by substituting Formulas (12) and (13) into Formula (11), the transformation parameter Pr for transforming the equatorial radius a of the spheroid EB in the WGS84 coordinate system into the radius R of the true sphere CB can be calculated as shown in the following Formula (14).
Next, a transformation parameter for transforming the longitude φ on the spheroid EB into the projected longitude Φ on the true sphere CB is calculated. The following Formula (15) is first obtained by Formula (8).
Formula (15)=Formula (8)×{Formula (1)}2/R, that is:
The following Formula (16) is obtained by substituting Formula (14) into Formula (15).
Next, the following Formula (17) is obtained by Formulas (7) and (10).
Formula (17)=Formula (10)/Formula (7), that is:
The following Formula (18) is obtained by Formulas (16) and (17).
Formula (18)={Formula (16)}2+{Formula (17)}2, that is:
Accordingly, a transformation parameter λφ for transforming the longitude φ on the spheroid EB represented by the following Formula (19) into the projected longitude Φ on the true sphere CB is obtained by Formula (18).
Therefore, the projected longitude Φ is represented by the following Formula (20).
Φ=λφ(φ−φ0) (20)
The coordinate transformation in Step S14 will be described in detail below. A method for transforming the reference latitude θ0 on the spheroid EB into a projection reference latitude Θ0 on the true sphere CB will be described. The following Formula (21) is obtained by Formulas (16) and (17).
Formula (21)=Formula (17)/Formula (16), that is:
Therefore, the projection reference latitude Θ0 is represented by the following Formula (22).
A method for transforming the latitude θ on the spheroid EB into the projected latitude Θ on the true sphere will be described. In order to project the spheroid EB onto the true sphere CB with as small a distortion as possible, a condition is set in which the latitude line spacing Dθ on the spheroid EB and the latitude line spacing Δθ on the true sphere CB are equal at the altitude h=0 and at a given latitude θ (equal latitude line spacing condition). In this case, Formula (7) is transformed to obtain the following Formula (23).
Formula (23) is transformed to obtain the following Formula (24).
Note that the above-mentioned set condition is merely an example. It is also possible to set other conditions, such as a condition in which the longitude line spacing Dφ on the spheroid EB and the longitude line spacing Δφ on the true sphere CB are equal (Δφ=Dφ) at a given latitude θ; a condition in which an area on the spheroid EB and an area on the true sphere CB are equal (equivalent) at a given latitude θ; and a condition in which a direction on the spheroid EB and a direction on the true sphere CB coincide with each other (conformal) at a given latitude θ.
where Π represents an elliptic integral of the third kind.
In the actual coordinate transformation, the operation using an ellipse function is complicated. Accordingly, the use of an approximation using Helmert's expansion (for example, Non Patent Literature 1) facilitates the coordinate transformation. A calculation using Helmert's expansion will be described below. A third flattening n is defined by the following Formula (25) using the flattening f.
Approximation (26) of Formula (24) is obtained by using Formula (25).
Here, coefficients in each term on the right side of Formula (26) are respectively defined as the following Formulas (27A) to (27F).
Formula (26) is transformed by using Formulas (27A) to (27F), and the projected latitude Θ on the true sphere CB is represented by the following Formula (28).
Θ(θ)=λθθ+λ2 sin 2θ+λ4 sin 4θ+λ6 sin 6θ+λ8 sin 8θ+λ0 (28)
The parameter information database D3 will be described in detail based on the above-mentioned calculation results of correction parameters.
Next, a change in the latitude line spacing will be described.
The three-dimensional positions of two points in a projection coordinate system on the true sphere CB facilitate the calculation of a distance d between a point P1 and a point P2 on the ground surface at the altitude h=0. However, to facilitate a subsequent numerical calculation, normalization vectors are used in Formulas (29) to (31).
The position of the point P1 is represented by the following Formula (29) based on Formula (4).
P
1=(X1,Y1,Z1)
X
1=cos Θ1 cos Φ1
Y
1=cos Θ1 sin Φ1
Z
1=sin Θ1 (29)
The position of the point P2 is represented by the following Formula (30) based on Formula (4).
P
2=(X2,Y2,Z2)
X
2=cos Θ2 cos Φ2
Y
2=cos Θ2 sin Φ2
Z
2=sin Θ2 (30)
Therefore, the distance d between the point P1 and the point P2 is represented by the following Formula (31).
d=R cos−1((P1·P2)) (31)
This corresponds to the distance calculation formula (following Formula (32)) in spherical trigonometry.
d=R cos−1(sin Θ1 sin Θ2+cos Θ1 cos Θ2 cos(Φ2−Φ1)) (32)
As shown in
Next, a change in the longitude line spacing will be described.
Next, an error rate of each of latitude line spacing and longitude line spacing when the reference latitude θ0 is changed will be described. The term “error rate” herein used refers to a value indicating how far the projected latitude line spacing and projected longitude line spacing obtained after the transformation is from the latitude line spacing and longitude line spacing on the spheroid EB.
As described above, the error rate Errθ of latitude line spacing is sufficiently small at each reference latitude. However, the generation mode of the error rate Errφ of longitude line spacing is changed by changing the reference latitude θ0.
Note that as shown in
If the reference latitude θ0 is 35 degrees north latitude, the error rate Errφ of longitude line spacing is less than 0.01% in a range from 13 degrees north latitude to 54 degrees north latitude, and the error rate Errφ of longitude line spacing is less than 0.04% in a range from the equator (0 degrees north latitude) to 63 degrees north latitude.
As shown in
Furthermore, the use of the coordinate transformation method according to this exemplary embodiment makes it possible to transform coordinates on the spheroid EB into coordinates on the true sphere CB with an error rate of less than 0.04% in a wide area ranging from the equator to the central Siberia in the north-south direction and from the Indian Ocean to the central Pacific Ocean in the east-west direction.
In a case where the geographical information management device 100 according to this exemplary embodiment is mounted in, for example, an airliner, if correction parameters are calculated by determining an appropriate reference latitude and an appropriate reference longitude, a relatively long-distance flight can be dealt with by performing the correction parameter calculation once, in view of a flight leg of a typical airliner.
For example, a flight starting from London will be considered. Flights from London to Tokyo or New York can be achieved with an error rate of 0.01% or less by performing the correction parameter calculation once. Flights from London to Rio de Janeiro can be achieved with an error rate of 0.04% or less by performing the correction parameter calculation once. In order to achieve flights from London to Rio de Janeiro with an error rate of 0.01% or less, it is necessary to calculate the correction parameters twice. The following table shows the latitude range of each operation route (flight leg), the reference latitude, and the error rate.
In the above table, the names of departure and destination airports are shown by ICAO (International Civil Aviation Organization) code. In the ICAO code, EGLL represents Heathrow Airport (London, England); RJAA represents New Tokyo International Airport (Narita Airport, Japan); John F. Kennedy International Airport (New York, U.S.A); and SBGL represents Antonio Carlos Jobim International Airport (Rio de Janeiro, Brazil).
On the other hand, in order to perform navigation computation with an error rate of less than 0.01% by using, for example, a map of orthonormal coordinates, it is necessary to perform a recalculation in about thirty minutes per degree. Thus, frequent calculations are required and a computer with high throughput is required to perform the calculations in vehicles such as airliners, which leads to an increase in size of the calculation system. Therefore, the navigation computation using a map of orthonormal coordinates is not realistic.
On the other hand, the geographical information management device 100 enables highly accurate coordinate transformation as described above by only calculating the correction parameters once every several hours. Therefore, the geographical information management device 100 can be easily mounted in vehicles, such as airlines, in which downsizing of the calculation system is required.
According to the above configuration, coordinates on the earth, which is a spheroid, can be projected onto a true sphere with a minimum error in an area larger than ever before. This enables the arithmetic processing using coordinate information to be performed on the true sphere on which mathematical processing can be easily performed. Consequently, the information about the operation of an aircraft can be accurately processed in a small device at a high speed.
Next, a geographical information management device 200 according to a second exemplary embodiment of the present invention will be described. The geographical information management device 200 has a configuration in which the function of reversely transforming projected coordinates on the true sphere CB into coordinates of the WGS84 coordinate system on the spheroid EB is added to the geographical information management device 100 according to the first exemplary embodiment.
In the geographical information management device 100, the transformation of coordinates of the WGS84 coordinate system on the spheroid EB into projected coordinates on the true sphere CB has been described above so as to facilitate arithmetic processing using the coordinate information. However, in order for a pilot of an airliner in which the geographical information management device is mounted, for example, to use the information obtained by the arithmetic processing at the projected coordinates on the true sphere CB to operate the airliner, it is necessary to inversely transform the projected coordinates into the coordinates of the WGS84 coordinate system based on which the airliner is actually flying. Accordingly, the geographical information management device 200 has the function of inversely transforming projected coordinates on the true sphere CB into coordinates of the WGS84 coordinate system on the spheroid EB.
The storage device 6 has a configuration in which a parameter information database D5 is added to the storage device 2. To simplify the drawing,
Next, the operation of the geographical information management device 200 will be described.
First, the operation unit 3 reads out a coordinate inverse transformation program PRG2. The coordinate inverse transformation program PRG2 is a program that transforms projected coordinates on the true sphere into coordinates on the spheroid by using the basic shape database D1, the projection information database D4, and the parameter information database D5. The coordinate inverse transformation program PRG2 is read out from, for example, the storage device 6.
Next, the operation unit 3 reads out the basic shape database D1, the projection information database D4, and the parameter information database D5.
The operation unit 3 generates a formula for performing coordinate transformation by substituting the information included in the basic shape database D1 and the parameter information database D5 into a formula specified by the coordinate inverse transformation program PRG2.
The operation unit 3 transforms coordinates on the true sphere into coordinates on the spheroid by substituting the projected coordinate information included in the projection information database D4 into the generated formula. In other words, the operation unit 3 inversely transforms the projection information database D4 into the airspace information database D2. The coordinate inverse transformation in Step S24 will be described in detail later.
The operation unit 3 outputs the coordinates on the spheroid EB included in the airspace information database D2 to the outside. For example, the operation unit 3 outputs the airspace information database D2 to the storage device 2. Further, for example, the operation unit 3 processes the coordinate information included in the airspace information database D2 into a form that can be displayed, and outputs the processed coordinate information to the display device 4.
The coordinate inverse transformation in Step S24 will be described in detail below. Specifically, the projected longitude φ on the true CB can be inversely transformed into the longitude φ on the spheroid EB by the following Formula (33) which is obtained by transforming Formula (20).
Specifically, the projected latitude Θ on the true sphere CB can be inversely transformed into the latitude θ on the spheroid EB by Formulas (26) and (27). However, even when it is intended to obtain a function for inversely transforming the projected latitude Θ on the spheroid EB into the latitude on the true sphere CB from Formula (26), the function cannot be described by an elementary function. Accordingly, the latitude θ on the spheroid EB may be obtained by a numerical calculation using Formula (26). As the numerical calculation, for example, the so-called Newton-Raphson method can be used.
In the case of using the Newton-Raphson method, it is necessary that Formula (26) can be differentiated by the latitude θ. A derivative Δ (θ) of the function Θ (θ) with respect to θ is represented by the following Formulas (34) and (35A) to (35E).
From the above formulas, it can be understood that the projected coordinates on the true sphere CB can be transformed into the coordinates on the spheroid EB by adding δ0, δ2, δ4, δ6, and δ8, which are respectively represented by Formulas (35A) to (35E), to the parameter information database. Specifically, the parameter information database D5 of this exemplary embodiment has a configuration in which δ0 to δ8 which are respectively represented by Formulas (35A) to (35E) are added to the parameter information database D3.
Note that the Newton-Raphson method is a general solution to equations by numerical calculation, and thus the detailed description thereof is omitted.
Next, a geographical information management device 300 according to a third exemplary embodiment of the present invention will be described. The geographical information management device 300 has the same configuration as that of the geographical information management device 100 according to the first exemplary embodiment.
In the case of transforming the latitude θ on the spheroid EB into the projected latitude Θ on the true sphere CB, the condition in which the latitude line spacing Dθ on the spheroid EB and the latitude line spacing Δθ on the true sphere CB are equal (Δθ=Dθ) at the altitude h=0 and at a given latitude θ (equal latitude line spacing condition) is set in the geographical information management device 100. On the other hand, in the geographical information management device 300, a condition in which the longitude line spacing Dφ on the spheroid EB and the longitude line spacing Δφ on the true sphere CB are equal (Δφ=Dφ) at the altitude h=0 and at a given latitude θ (equal longitude line spacing condition) is set. In this case, the following Formula (36) is obtained by transforming Formula (8).
Formula (36) is transformed by using Formula (20), to thereby obtain the following Formula (37).
By Formula (37), the projected latitude Θ can be expressed as the following Formula (38).
Other calculations for transforming the latitude θ on the spheroid EB into the projected latitude Θ on the true sphere CB are similar to those of the first exemplary embodiment, and thus a description thereof is omitted.
In view of the above, it can be understood that even in the case where the condition in which the longitude line spacing Dφ on the spheroid EB and the longitude line spacing Δφ on the true sphere CB are equal (Δφ=Dφ) at the altitude h=0 and at a given latitude θ (equal longitude line spacing condition) is set, the spheroid EB can be projected onto the true sphere CB with as small a distortion as possible.
Next, a geographical information management device 400 according to a fourth exemplary embodiment of the present invention will be described. The geographical information management device 400 has the same configuration as that of the geographical information management device 100 according to the first exemplary embodiment.
In the case of transforming the latitude θ on the spheroid EB into the projected latitude Θ on the true sphere CB, the condition in which the latitude line spacing Dθ on the spheroid EB and the latitude line spacing Δθ on the true sphere CB are equal (Δθ=Dθ) at the altitude h=0 and at a given latitude θ (equal latitude line spacing condition) is set in the geographical information management device 100. On the other hand, in the geographical information management device 400, a condition in which an area DθDφ formed by the latitude line spacing and the longitude line spacing on the spheroid EB and an areas ΔθΔφ formed by the latitude line spacing and the longitude line spacing on the true sphere CB are equal (ΔθΔφ=DθDφ) at the altitude h=0 and at a given latitude θ (equivalence condition) is set. In this case, Formulas (7) and (8) are transformed to obtain the following Formula (39).
Formula (39) is transformed by using Formula (20), to thereby obtain the following Formula (40).
By Formula (40), the projected latitude Θ can be expressed as the following Formula (41).
Other calculations for transforming the latitude θ on the spheroid EB into the projected latitude Θ on the true sphere CB are similar to those of the first exemplary embodiment, and thus a description thereof is omitted.
In view of the above, it can be understood that even in the case where the condition in which the area DθDφ formed by the latitude line spacing and the longitude line spacing on the spheroid EB and the area ΔθΔφ formed by the latitude line spacing and the longitude line spacing on the true sphere CB are equal (ΔθΔφ=DθDφ) at the altitude h=0 and at a given latitude θ (equivalence condition) is set, the spheroid EB can be projected onto the true sphere CB with as small a distortion as possible.
Next, a geographical information management device 500 according to a fifth exemplary embodiment of the present invention will be described. The geographical information management device 500 has the same configuration as that of the geographical information management device 100 according to the first exemplary embodiment.
In the case of transforming the latitude θ on the spheroid EB into the projected latitude Θ on the true sphere CB, the condition in which the latitude line spacing Dθ on the spheroid EB and the latitude line spacing Δθ on the true sphere CB are equal (Δθ=Dθ) at the altitude h=0 and at a given latitude θ (equal latitude line spacing condition) is set in the geographical information management device 100. On the other hand, in the geographical information management device 500, a condition in which an angle Dθ/Dφ formed by the latitude line spacing and the longitude line spacing on the spheroid EB and an angle Δθ/Δφ formed by the latitude line spacing and the longitude line spacing on the true sphere CB are equal (Δθ/Δφ=Dθ/Dφ) at the altitude h=0 and at a given latitude θ (conformal condition) is set. In this case, Formulas (7) and (8) are transformed to obtain the following Formula (42).
Formula (42) is transformed by using Formula (20), to thereby obtain the following Formula (43).
By Formula (43), the projected latitude Θ is defined as the following Formula (44).
Other calculations for transforming the latitude θ on the spheroid EB into the projected latitude Θ on the true sphere CB are similar to those of the first exemplary embodiment, and thus a description thereof is omitted.
In view of the above, it can be understood that even in the case where the condition in which the angle Dθ/Dφ formed by the latitude line spacing and the longitude line spacing on the spheroid EB and the angle Δθ/Δφ formed by the latitude line spacing and the longitude line spacing on the true sphere CB are equal (Δθ/Δφ=Dθ/Dφ) at the altitude h=0 and at a given latitude θ is set, the spheroid EB can be projected with as small a distortion as possible.
Next, a geographical information management device 600 according to a sixth exemplary embodiment of the present invention will be described. The geographical information management device 600 has the same configuration as that of the geographical information management device 100 according to the first exemplary embodiment.
In the case of transforming the latitude θ on the spheroid EB into the projected latitude Θ on the true sphere CB, the condition in which the latitude line spacing Dθ on the spheroid EB and the latitude line spacing Δθ on the true sphere CB are equal (Δθ=Dθ) at the altitude h=0 and at a given latitude θ (equal latitude line spacing condition) is set in the geographical information management device 100. On the other hand, in the geographical information management device 600, a condition in which a function Dθf(Dφ) of the latitude line spacing and the longitude line spacing on the spheroid EB and a function Δθf(Δφ) of the latitude line spacing and the longitude line spacing on the true sphere CB are equal at the altitude h=0 and at a given latitude θ is set. In this case, Formulas (7) and (8) are transformed to obtain the following Formula (45).
Assuming that primitive functions are respectively represented by FD and FΔ (especially, when f represents a rational function, the integration can be analytically performed without fail), the relationship between the latitude θ on the spheroid EB and the projected latitude Θ on the true sphere CB is defined as the following Formula (46).
aF
D(θ)−aFD(θ0)=RFΔ(Θ)−RFΔ(Θ0) (46)
Other calculations for transforming the latitude θ on the spheroid EB into the projected latitude Θ on the true sphere CB are similar to those of the first exemplary embodiment, and thus a description thereof is omitted.
In view of the above, it can be understood that even in the case where the condition in which the function Dθf(Dφ) of the latitude line spacing and the longitude line spacing on the spheroid EB and the function Δθf(Δφ) of the latitude line spacing and the longitude line spacing on the true sphere are equal at the altitude h=0 and at a given latitude θ is set, the spheroid EB can be projected onto the true sphere CB with as small a distortion as possible.
Next, a geographical information management device 700 according to a seventh exemplary embodiment of the present invention will be described. The geographical information management device 700 has the same configuration as that of the geographical information management device 100 according to the first exemplary embodiment. In this exemplary embodiment, an example of the operation method of the geographical information management device 700 on the true sphere CB will be described. The operation method described in this exemplary embodiment is carried out by the operation unit 3. In the case of setting an airspace for aircraft training or the like, an airway along which an aircraft is operated, and the like by using coordinates on the true sphere CB that are obtained by transforming coordinates on the spheroid EB, it is necessary to perform an appropriate operation on the true sphere CB. Examples of various operations on the true sphere CB will be described below.
Note that in formulas used in the following description, a vector quantity is represented by an arrow above a symbol. To simplify the explanation, all vector quantities are normalized. In the drawings described below, the true sphere CB indicates the earth and EQ represents the equator of the earth.
For example, in order to obtain a required flight time of an aircraft, it is essential to obtain a distance between two points. The calculation of the distance between the point P1 and the point P2 on the true sphere CB will now be described. An arc-like path with the distance d from the point P1 as a start point to the point P2 as an end point is set on the true sphere CB.
The position vector of the point P2 can be defined by the above Formula (30).
The distance d between the point P1 and the point P2 on the true sphere CB (on the ground surface) can be obtained by the above Formula (31) by using Formulas (29) and (30). The above Formula (32) can be obtained by expanding Formula (31). Formula (32) corresponds to the distance calculation formula in spherical trigonometry. This makes it possible to obtain the distance d between the point P1 and the point P2.
The calculation of an azimuth Ψ between the point P1 and the point P2 on the true sphere CB (on the ground surface) will be described.
{right arrow over (N)}=(0,0,1) (47)
The plane PL2 is a plane that is perpendicular to the position vector of the point P1 and includes the center of the earth. The plane PL2 is a plane that is parallel to the vector of the cross product between the position vector of the pole N and the position vector of the point P1, and is also parallel to the vector of the cross product between this cross product and the position vector of the vector and point P1.
By using the position vectors of the point P1 and the point P2, x is defined by the following Formula (48) and y is defined by the following Formula (49).
x=({right arrow over (P2)}·{right arrow over (N)}−({right arrow over (N)}·{right arrow over (P1)}){right arrow over (P1)}) (48)
y=({right arrow over (P2)}·{right arrow over (N)}×{right arrow over (P1)}) (49)
When y>0 holds, the azimuth Ψ from the point P1 to the point P2 on the true sphere CB (on the ground surface) is represented by the following Formula (50).
When y<0 holds, the azimuth Ψ from the point P1 to the point P2 on the true sphere CB (on the ground surface) is represented by the following Formula (51).
When y=0 and x>0 hold, the azimuth Ψ from the point P1 to the point P2 on the true sphere CB (on the ground surface) is represented by the following Formula (52).
ψ=0 (52)
When y=0 and x<0 hold, the azimuth Ψ from the point P1 to the point P2 on the true sphere CB (on the ground surface) is represented by the following Formula (53).
ψ=π (53)
A minimum path between two points on the true sphere and a method for interpolating the minimum path
A minimum path between the point P1 and the point P2 on the true sphere CB (on the ground surface) and a method for interpolating the minimum path will be described.
The case of interpolating the minimum path between the point P1 and the point P2 will now be considered. x is defined by the following Formula (55) and y is defined by the following Formula (56). In this case, (x, y) represents a two-dimensional vector which is obtained by projecting the point P2 onto a plane including the point P1 and the point P2, assuming that a component parallel to the point P1 is plotted on an x-axis and a component perpendicular to the point P1 is plotted on a y-axis.
x=({right arrow over (P2)}·{right arrow over (P1)}) (55)
y=({right arrow over (P2)}·{right arrow over (Va)}×{right arrow over (P1)}) (56)
By Formulas (55) and (56), a central angle Ψ of an arc that connects the point P1 and the point P2 on the plane including the point P1 and the point P2 is represented by the following Formula (57).
That is, the segment to be interpolated corresponds to the point P that satisfies each vector equation in Formula (54), and the segment is within the central angle Ψ of the arc that connects the point P1 and the point P2. Accordingly, when a parameter t is used, the position vector of an interpolation point P(t) between the point P1 and the point P2 can be defined by the following formula (58).
{right arrow over (P(t))}={right arrow over (P1)} cos t+({right arrow over (Va)}×{right arrow over (P1)})sin t(0≦t≦ψ) (58)
In Formula (58), the minimum path between the point P1 and the point P2 can be interpolated by setting the parameter t as appropriate.
A Latitude Line that Connects Two Points at the Same Latitude and Interpolation the Latitude Line
A latitude line that connects the point P1 and the point P2 at the same latitude on the true sphere CB (on the ground surface) and interpolation of the latitude line will be described. A latitude line on the true sphere CB (on the ground surface) can be recognized as a rhumb line between two points at the same latitude on the true sphere CB.
A case where the azimuth from the point P1 (start point) to the point P2 (end point) is east will be described.
{right arrow over (Vb)}={right arrow over (N)}=(0,0,1)
(Vb·P)=cos α
|Vb·P|=sin α (59)
α represents an angle formed between the pole N and the latitude Θ (θ) of each of the point P1 and the point P2, and is represented by the following Formula (60).
A case where the azimuth from the point P1 (start point) to the point P2 (end point) is west will be described.
V
c
=S=(0,0,−1)
(Vc·P)=cos α
|Vc×P|=sin α (61)
α represents an angle formed between the pole S and the latitude Θ (θ) of each of the point P1 and the point P2, and is represented by the following Formula (62).
By using the position vector of the point P1 and the position vector of the point P2, x is defined by the following Formula (63) and y is defined by the following Formula (64).
x=({right arrow over (P2)}·{right arrow over (P1)}−({right arrow over (V)}·{right arrow over (P1)}){right arrow over (V)}) (63)
y=({right arrow over (P2)}·{right arrow over (V)}×{right arrow over (P1)}) (64)
The case of interpolating the latitude line between the point P1 and the point P2 will be considered. When the latitude line to be interpolated is a minor arc, the central angle Ψ of the arc that connects the point P1 and the point P2 can be expressed by the following Formula (65).
When the latitude line to be interpolated is a major arc, the central angle Ψ of the arc that connects the point P1 and the point P2 can be expressed by the following Formula (66).
When the latitude line to be interpolated is a semi-circumference, the central angle Ψ of the arc that connects the point P1 and the point P2 can be expressed by the following Formula (67).
ψ=π (67)
Specifically, when the azimuth is east, the segment to be interpolated corresponds to the point P that satisfies each vector equation in Formulas (59) and (60), and also corresponds to the segment of the azimuth Ψ between the point P1 and the point P2. When the azimuth Ψ is west, the segment to be interpolated corresponds to the point P that satisfies each vector equation in Formulas (61) and (62), and also corresponds to the segment within the central angle Ψ of the arc that connects the point P1 and the point P2. Accordingly, when the parameter t is used, the position vector of the interpolation point P(t) between the point P1 and the point P2 can be defined by the following Formula (68). In Formula (68), the unit normal vector V is one of the unit normal vector Vb and the unit normal vector Vc.
{right arrow over (P(t))}={right arrow over (V)} cos α+({right arrow over (P1)}−({right arrow over (V)}·{right arrow over (P1)})cos t+({right arrow over (V)}×{right arrow over (P1)})sin t(0≦t≦ψ) (68)
In Formula (68), the rhumb line between the point P1 and the point P2 at the same latitude line can be interpolated by setting the parameter t as appropriate.
A Path Obtained by Projecting a Circle on the Spheroid onto the True Sphere and a Method for Interpolating the Path
A path on the true sphere CB that is obtained when a circle on the spheroid EB is projected onto the true sphere CB, and a method for interpolating the path will be described.
{right arrow over (Vd)}={right arrow over (P0)}
({right arrow over (Vd)}·{right arrow over (P)})=cos α
|Vd×P|=sin α (69)
α represents an angle formed between the point P0 and the point P on the true sphere, and the angle can be represented by the following Formula (70).
In the case of interpolating a circle, a point on a uniform circumference is set as the point P1 which serves as an interpolation start point and an interpolation end point. For example, the point P1 is a point where the circumference intersects with the longitude line that connects the center of the circle and the pole. In this case, the point P1 can be represented by the following Formula (71). Note that as described above, the pole N represents the pole (north pole) of the true sphere CB.
The central angle Ψ of the arc on the path C1 is 2π as represented by Formula (72).
ψ=2π (72)
The case of interpolating the path C1 will now be considered. The segment to be interpolated corresponds to the point P that satisfies each vector equation in Formulas (69) and (70). Accordingly, when the parameter t is used, the position vector of the interpolation point P(t) on the path C1 can be defined by the following Formula (73).
{right arrow over (P(t))}={right arrow over (P0)} cos α+({right arrow over (P1)}−({right arrow over (P0)}·{right arrow over (P1)}){right arrow over (P0)})cos t+({right arrow over (P0)}×{right arrow over (P1)})sin t(0≦t≦ψ) (73)
In Formula (73), the path C1 that is obtained by projecting a circle on the spheroid EB onto the true sphere CB can be interpolated by setting the parameter t as appropriate.
A Path Obtained by Projecting an Arc Connecting Two Points on the Spheroid onto the True Sphere and a Method for Interpolating the Path
A path on the true sphere CB that is obtained when an arc on the spheroid EB is projected onto the true sphere CB, and a method for interpolating the path will be described. An arc which has the radius r, connects two points, and is centered on a certain point P0WGS on the spheroid EB can be interpolated as a set of points with the distance r from the point P0 on the true sphere CB that is obtained by projecting the point P0WGS.
A case where the direction from the start point to the end point of the arc is counterclockwise will be described.
{right arrow over (Vd)}={right arrow over (P0)}
({right arrow over (Vd)}·{right arrow over (P)})=cos α
|{right arrow over (Vd)}×{right arrow over (P)}|=sin α (74)
α represents an angle formed between the point P0 and the interpolation point P on the true sphere, and is represented by the following Formula (75).
A case where the direction from the start point to the end point of the arc is clockwise will be described.
{right arrow over (Ve)}=−{right arrow over (P0)}
(Ve·P)=cos α
|Ve×P|=sin α (76)
α represents an angle formed between the point P0 and the interpolation P on the true sphere, and is represented by the following Formula (77).
A case where the minimum path between the point P1 (start point) and the point P2 (end point) of each of the paths C2 and C3 will now be considered. The following description is common to both cases in which the direction of the path is counterclockwise and the direction of the path is clockwise. By using the position vector of the point P1 and the position vector of the point P2, x is defined by the following Formula (78) and y is defined by the following Formula (79). Note that the point P1 is the start point of the path C2 or the path C3 and the point P2 is the end point of the path C2 or the path C3. In Formulas (78) and (79), the unit normal vector V is one of the unit normal vector Vd and the unit normal vector Ve.
x=({right arrow over (P2)}·{right arrow over (P1)}−({right arrow over (V)}·{right arrow over (P1)}){right arrow over (V)}) (78)
y=({right arrow over (P2)}·{right arrow over (V)}·{right arrow over (P1)}) (79)
Based on Formulas (78) and (79), when the arc that connects the point P1 and the point P2 on the plane including the point P1 and the point P2 is a minor arc, the central angle Ψ of the arc can be represented by the following Formula (80).
When the arc that connects the point P1 and the point P2 on the plane including the point P1 and the point P2 is a major arc, the central angle Ψ of the arc can be represented by the following Formula (81).
When the arc that connects the point P1 and the point P2 on the plane including the point P1 and the point P2 is a semi-circumference, the central angle Ψ of the arc can be represented by the following Formula (82).
ψ=π (82)
Specifically, when the path is counterclockwise (path C2), the segment to be interpolated corresponds to the point P that satisfies Formulas (74) and (75), and the segment is within the central angle Ψ of the arc that connects the point P1 and the point P2. When the path is clockwise (path C3), the segment to be interpolated corresponds to the point P that satisfies Formulas (76) and (77), and the segment is within the central angle Ψ of the arc that connects the point P1 and the point P2. Accordingly, when the parameter t is used, the position vector of the interpolation point P(t) between the point P1 and the point P2 can be defined by the following Formula (83).
{right arrow over (P(t))}={right arrow over (P0)} cos α+({right arrow over (P1)}−({right arrow over (P0)}·{right arrow over (P1)}){right arrow over (P0)})cos t+({right arrow over (P0)}×{right arrow over (P1)})sin t(0≦t≦ψ) (83)
When paths which can be expressed as a great circle, a latitude line, a circle, and an arc on the spheroid EB are projected onto the true sphere CB with a small error, the paths can be described using similar vector equations with a minimum error in a wide range. Accordingly, the interpolation method on the true sphere EB can be generalized regardless of the type of paths. The generalized interpolation method can be expressed as the following Formula (84) using the parameter t. The point P1 is the start point of the path to be interpolated, and is represented as a position vector in the formula. V represents a unit normal vector for a plane to which the path to be interpolated belongs.
{right arrow over (P(t))}=V cos α+({right arrow over (P1)}−({right arrow over (V)}·{right arrow over (P1)}){right arrow over (V)})cos t+({right arrow over (V)}×{right arrow over (P1)})sin t(0≦t≦ψ) (84)
A block distance L of the path can be represented by the following Formula (85).
L=Rψ sin α (85)
On the basis of a minimum path, a latitude line, a circle, or an arc defined on the spheroid EB, a line segment having a constant distance d is defined on the side (one or both of the right side and the left side) of the minimum path, the latitude line, the circle, or the arc in some cases. This line segment is generally referred to as an offset line segment.
By defining the offset line segment on both sides of a certain path, for example, an area having a constant spacing d on both sides of the path can be held as a buffering area. The safety of aircraft, such as avoidance of a collision between aircraft, can be ensured by preventing this buffering area from overlapping other paths. The offset line segment may be set so as to hold a space with a constant width on the basis of a reference line.
{right arrow over (V)}={right arrow over (Vf)}
({right arrow over (Vf)}·{right arrow over (P)})=cos α
|{right arrow over (Vf)}×{right arrow over (P)}|=sin α (86)
αc represents an angle formed between Vf and a given point on the reference path. The positive (negative) sign corresponds to the case where an offset is taken on the right (left) side with respect to the reference line. α in Formula (86) is represented by the following Formula (87) using αc.
The point P is represented by the following Formula (88).
That is, when the parameter t is used, the position vector of the interpolation point P(t) can be defined by the following Formula (89).
{right arrow over (P(t))}={right arrow over (V)} cos α+({right arrow over (P1)}−({right arrow over (Vf)}·{right arrow over (P1)}){right arrow over (Vf)})cos t+({right arrow over (Vf)}×{right arrow over (P1)})sin t(0≦t≦ψ) (89)
In Formula (89), the offset line segment can be interpolated by setting the parameter t as appropriate.
Basically, the interpolation of the offset line segment at inflection points of the continuous reference line segment can be performed in the same manner as in the interpolation of an arc.
{right arrow over (V)}={right arrow over (P0)}
({right arrow over (V)}·{right arrow over (P)})=cos α
|{right arrow over (V)}×{right arrow over (P)}|=sin α (90)
α represents a value obtained by dividing the offset distance d by the radius R of the true sphere, and is represented by the following Formula (91).
The position vector of the point P1 is represented by the following Formula (92).
The position vector of the point P2 is defined as the following Formula (93).
When Formula (92) is substituted into Formula (90), the following Formula (94) is obtained.
When Formula (93) is substituted into Formula (90), the following Formula (95) is obtained.
In Formulas (94) and (95), the double sign in parentheses indicates inflection to the right (+) and inflection to the left (−).
By using the above Formulas (94) and (95), the interpolation between the endpoints (the point P1 and the point P2) of the arc portion of the offset line segment in the vicinity of the inflection point can be performed by performing the line segment interpolation described above. Specifically, the offset line segment in the vicinity of the inflection point of the continuous line segment can be interpolated by using the above Formula (83).
Basically, the interpolation of endpoints of the offset line segment with respect to the starting and terminating ends of the reference line segment can be performed in the same manner as in the interpolation of an arc.
{right arrow over (V)}={right arrow over (P0)}
({right arrow over (V)}·{right arrow over (P)})=cos α
|{right arrow over (V)}×{right arrow over (P)}|=sin α (96)
α represents a value obtained by dividing the offset distance d by the radius R of the true sphere, and is represented by the following Formula (97).
In this case, the position vector of the endpoint P1 is represented by the following Formula (98).
In this case, the position vector of the endpoint P2 is represented by Formula (99).
The interpolation between the endpoints (between the point P1 and the point P2) at the ends of the offset line segment can be performed by interpolating the above-mentioned line segment by using the above Formulas (98) and (99). Specifically, the offset line segment at the endpoints of the reference line segment can be interpolated using the above Formula (83).
Next, a geographical information management device 800 according to an eighth exemplary embodiment of the present invention will be described. The geographical information management device 800 has the same configuration as that of the geographical information management device 100 according to the first exemplary embodiment. This exemplary embodiment illustrates an example in which an airspace is set on the true sphere by using the operation method described in the seventh exemplary embodiment.
Airspaces have various shapes. For example, in the case of a square airspace, the airspace can be defined by four line segments on the true sphere CB. In the case of a circular airspace, the airspace can be defined by a path obtained by projecting a circle on the spheroid EB onto the true sphere CB. In the case of a fan-shaped airspace, the airspace can be defined by a path obtained by projecting an arc on the spheroid EB onto the true sphere CB and two line segments on the true sphere CB. In the case of an airspace having a predetermined width with respect to a certain line segment, the airspace can be defined by two offset line segments with respect to a predetermined line segment on the true sphere CB and two line segments on the true sphere CB that connect the endpoints of the two offset line segments.
Specific examples of airspaces will be described below.
The airspace A1 is formed of 22 inflection points PP1 to PP22 and 21 line segments that connect the inflection points to each other. The line segments are each interpolated by the method described in the sixth exemplary embodiment.
The airspace A2 which is obtained by linear interpolation as a comparative example is formed by linear interpolation between the adjacent inflection points of the 22 inflection points PP1 to PP22 by using the Universal Transverse Mercator (UTM) projection. Note that since the transformation error is large in the UTM projection, the UTM projection can be applied only to a latitude range of about ±3°. Accordingly, the UTM projection is not suitable for a wide airspace such as the flight information region in which Japan is responsible for the air traffic control. Therefore, in this exemplary embodiment, the UTM projection, which is generally used, is employed as a comparative example so as to understand the superiority of the projection onto the true sphere EB according to the above exemplary embodiments.
As shown in
If such a gap is generated, various problems occur in the management of the airspace. First, for example, in the case of using the UTM projection, a difference in normal coordinates for transformation may result in a variation in error. For example, in the case of setting the flight information region under the control of each country, if each country employs different normal coordinates using coordinate systems that are convenient for their own purposes, a malfunction, such as overlapping of the airspaces that are not supposed to overlap each other, occurs. As a result, the airspace information varies from country to country, which makes it difficult to share the airspace information between different air traffic control authorities. As a result, the airspaces and routes cannot be optimized in a wide area.
Second, there can be a case where it comes to a conclusion, due to an error, that routes that originally overlap each other are not overlapping. In this case, there is a possibility that a plurality of aircraft enter the same airspace at the same time, which poses a problem in terms of ensuring the safety of aircraft.
Third, even in the case of setting a route that has a minimum path, a route deviating from the minimum path may be set due to an error. In this case, the aircraft is caused to fly an extra distance, which results in an increase in fuel consumption. This is disadvantageous from the viewpoint of economic operation of aircraft.
These problems can be solved by geodesic interpolation using coordinates on the spheroid that is approximate to the actual shape of the earth. However, as described above, the geodesic interpolation using coordinates on the spheroid requires complicated calculations such as elliptic integral. Therefore, it is necessary to use techniques such as sectional measurement and repeated calculation. This results in the necessity of a large number of calculation resources.
Accordingly, in this exemplary embodiment, the design of an airspace and setting of a route can be easily performed by only projecting coordinates on the spheroid onto the true sphere by calculations of, for example, the interpolation on the true sphere as described in the sixth exemplary embodiment, to be specific, elementary functions and arithmetic operations. Further, as described in the above exemplary embodiments, the error in the interpolation on the true sphere is smaller than that in a general method, such as the UTM projection, and thus is suitable for an application to a wide area. Furthermore, unlike in the UTM projection, a variation in error due to a difference in normal coordinates is small. This facilitates sharing of information among, for example, flight information regions.
Note that the present invention is not limited to the above exemplary embodiments and can be modified as appropriate without departing from the scope of the invention. For example, the above exemplary embodiments have illustrated the case where the geographical information management device, which is an example of the coordinate transformation device, is mounted on a vehicle such as an aircraft that travels above the earth. However, this is merely an example. The coordinate transformation device (geographical information management device) according to exemplary embodiments described above can be mounted on other vehicles that travel on the earth, such as vehicles that fly through the air, other than aircraft, vehicles that travel on the land, ships that travel on the sea, and submersible vehicles that travel in the sea. The coordinate transformation device (geographical information management device) according to exemplary embodiments described above can be incorporated not only in vehicles, but also in an operation management system and the like of vehicles, such as a control system for controlling aircraft.
In the above exemplary embodiments, the present invention is described as a hardware configuration, but the present invention is not limited to this. According to the present invention, any processing can be implemented by causing a CPU (Central Processing Unit) to execute a computer program. The program can be stored and provided to a computer using any type of non-transitory computer readable media. Non-transitory computer readable media include any type of tangible storage media. Examples of non-transitory computer readable media include magnetic storage media (such as floppy disks, magnetic tapes, hard disk drives, etc.), optical magnetic storage media (e.g. magneto-optical disks), CD-ROM (Read Only Memory), CD-R, CD-R/W, and semiconductor memories (such as mask ROM, PROM (Programmable ROM), EPROM (Erasable PROM), flash ROM, RAM (Random Access Memory), etc.). The program may be provided to a computer using any type of transitory computer readable media. Examples of transitory computer readable media include electric signals, optical signals, and electromagnetic waves. Transitory computer readable media can provide the program to a computer via a wired communication line, such as electric wires and optical fibers, or a wireless communication line.
While the present invention has been described above with reference to exemplary embodiments, the present invention is not limited to the above exemplary embodiments. The configuration and details of the present invention can be modified in various ways which can be understood by those skilled in the art within the scope of the invention.
This application is based upon and claims the benefit of priority from Japanese patent application No. 2013-001110, filed on Jan. 8, 2013, and Japanese patent application No. 2013-178833, filed on Aug. 30, 2013, the disclosure of which is incorporated herein in its entirety by reference.
Number | Date | Country | Kind |
---|---|---|---|
2013-001110 | Jan 2013 | JP | national |
2013-178833 | Aug 2013 | JP | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/JP2013/007221 | 12/9/2013 | WO | 00 |