1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
|
//>>built
define("dojox/geo/openlayers/GreatCircle", ["dojo/_base/lang",
"dojox/geo/openlayers/GeometryFeature",
"dojox/geo/openlayers/Point",
"dojox/geo/openlayers/LineString"], function(lang, GeometryFeature, Point, lineString){
lang.getObject("geo.openlayers", true, dojox);
dojox.geo.openlayers.GreatCircle = {
toPointArray : function(p1, p2, increment){
// summary:
// Create a geodetic line as an array of OpenLayers.Point.
// descritpion:
// Create a geodetic line as an array of OpenLayers.Point between the point p1
// and the point p2. Result is a polyline approximation for which a new point is
// calculated every <em>increment</em> degrees.
// p1: Point
// The first point of the geodetic line. x and y fields are longitude and
// latitude in decimal degrees.
// p2: Point
// The second point of the geodetic line. x and y fields are longitude and
// latitude in decimal degrees.
// increment: Float
// The value at which a new point is computed.
var startLon = p1.x;
var endLon = p2.x;
var sl = Math.min(startLon, endLon);
var el = Math.max(startLon, endLon);
var d2r = this.DEG2RAD;
var lat1 = p1.y * d2r;
var lon1 = p1.x * d2r;
var lat2 = p2.y * d2r;
var lon2 = p2.x * d2r;
if (Math.abs(lon1 - lon2) <= this.TOLERANCE) {
var l = Math.min(lon1, lon2);
lon2 = l + Math.PI;
}
if (Math.abs(lon2 - lon1) == Math.PI) {
if (lat1 + lat2 == 0.0) {
lat2 += Math.PI / 180000000;
}
}
var lon = sl * d2r;
var elon = el * d2r;
var incr = increment * d2r;
var wp = [];
var k = 0;
var r2d = this.RAD2DEG;
while (lon <= elon) {
lat = Math.atan((Math.sin(lat1) * Math.cos(lat2) * Math.sin(lon - lon2) - Math.sin(lat2) * Math.cos(lat1)
* Math.sin(lon - lon1))
/ (Math.cos(lat1) * Math.cos(lat2) * Math.sin(lon1 - lon2)));
var p = {
x : lon * r2d,
y : lat * r2d
};
wp[k++] = p;
if (lon < elon && (lon + incr) >= elon)
lon = elon;
else
lon = lon + incr;
}
return wp;
},
toLineString : function(p1, p2, increment){
// summary:
// Create a geodetic line as an array of OpenLayers.Geometry.LineString.
// descritpion:
// Create a geodetic line as a OpenLayers.Geometry.LineString between the point p1
// and the point p2. Result is a polyline approximation for which a new point is
// calculated every <em>increment</em> degrees.
// p1: Point
// The first point of the geodetic line. x and y fields are longitude and
// latitude in decimal degrees.
// p2: Point
// The second point of the geodetic line. x and y fields are longitude and
// latitude in decimal degrees.
// increment: Float
// The value at which a new point is computed.
var wp = this.toPointArray(p1, p2, increment);
var ls = new OpenLayers.Geometry.LineString(wp);
return ls;
},
toGeometryFeature : function(p1, p2, increment){
// summary:
// Create a geodetic line as an array of dojox.geo.openlayers.GeometryFeature.
// description:
// Create a geodetic line as a dojox.geo.openlayers.GeometryFeature between the point p1
// ant the point p2. Result is a polyline approximation for which a new point is
// calculated every <em>increment</em> degrees.
// p1: Point
// The first point of the geodetic line. x and y fields are longitude and
// latitude in decimal degrees.
// p2: Point
// The second point of the geodetic line. x and y fields are longitude and
// latitude in decimal degrees.
// increment: Float
// The value at which a new point is computed.
// returns: GeometryFeature
// The geodetic line as a GeometryFeature
var ls = this.toLineString(p1, p2, increment);
return new GeometryFeature(ls);
},
DEG2RAD : Math.PI / 180,
RAD2DEG : 180 / Math.PI,
TOLERANCE : 0.00001
};
return dojox.geo.openlayers.GreatCircle;
});
|
