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
|
// Verify that odd polygons (huge or "narrow") behave as we expect.
// Note that since 2dsphere is spherical, polygons that seem narrow are actually
// rather wide if their latitude (or longitude) range is large.
// @tags: [
// sbe_incompatible,
// ]
var t = db.geo_s2oddshapes;
t.drop();
t.createIndex({geo: "2dsphere"});
var testPoint = {name: "origin", geo: {type: "Point", coordinates: [0.0, 0.0]}};
var testHorizLine = {
name: "horiz",
geo: {type: "LineString", coordinates: [[-2.0, 10.0], [2.0, 10.0]]}
};
var testVertLine = {
name: "vert",
geo: {type: "LineString", coordinates: [[10.0, -2.0], [10.0, 2.0]]}
};
t.insert(testPoint);
t.insert(testHorizLine);
t.insert(testVertLine);
// Test a poly that runs vertically all the way along the meridian.
var tallPoly = {
type: "Polygon",
coordinates: [[[1.0, 89.0], [-1.0, 89.0], [-1.0, -89.0], [1.0, -89.0], [1.0, 89.0]]]
};
// We expect that the testPoint (at the origin) will be within this poly.
var result = t.find({geo: {$within: {$geometry: tallPoly}}});
assert.eq(result.itcount(), 1);
var result = t.find({geo: {$within: {$geometry: tallPoly}}});
assert.eq(result[0].name, 'origin');
// We expect that the testPoint, and the testHorizLine should geoIntersect
// with this poly.
result = t.find({geo: {$geoIntersects: {$geometry: tallPoly}}});
assert.eq(result.itcount(), 2);
result = t.find({geo: {$geoIntersects: {$geometry: tallPoly}}});
// Test a poly that runs horizontally along the equator.
var longPoly = {
type: "Polygon",
coordinates: [[[89.0, 1.0], [-89.0, 1.0], [-89.0, -1.0], [89.0, -1.0], [89.0, 1.0]]]
};
// Thanks to spherical geometry, this poly contains most of the hemisphere.
result = t.find({geo: {$within: {$geometry: longPoly}}});
assert.eq(result.itcount(), 3);
result = t.find({geo: {$geoIntersects: {$geometry: longPoly}}});
assert.eq(result.itcount(), 3);
// Test a poly that is the size of half the earth.
t.drop();
t.createIndex({geo: "2dsphere"});
var insidePoint = {name: "inside", geo: {type: "Point", name: "inside", coordinates: [100.0, 0.0]}};
var outsidePoint = {
name: "inside",
geo: {type: "Point", name: "inside", coordinates: [-100.0, 0.0]}
};
t.insert(insidePoint);
t.insert(outsidePoint);
var largePoly = {
type: "Polygon",
coordinates: [[[0.0, -90.0], [0.0, 90.0], [180.0, 0], [0.0, -90.0]]]
};
result = t.find({geo: {$within: {$geometry: largePoly}}});
assert.eq(result.itcount(), 1);
result = t.find({geo: {$within: {$geometry: largePoly}}});
var point = result[0];
assert.eq(point.name, 'inside');
// Test a poly that is very small. A couple meters around.
t.drop();
t.createIndex({geo: "2dsphere"});
insidePoint = {
name: "inside",
geo: {type: "Point", name: "inside", coordinates: [0.01, 0.0]}
};
outsidePoint = {
name: "inside",
geo: {type: "Point", name: "inside", coordinates: [0.2, 0.0]}
};
t.insert(insidePoint);
t.insert(outsidePoint);
smallPoly = {
type: "Polygon",
coordinates: [[[0.0, -0.01], [0.015, -0.01], [0.015, 0.01], [0.0, 0.01], [0.0, -0.01]]]
};
result = t.find({geo: {$within: {$geometry: smallPoly}}});
assert.eq(result.itcount(), 1);
result = t.find({geo: {$within: {$geometry: smallPoly}}});
point = result[0];
assert.eq(point.name, 'inside');
|