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
|
t = db.geo_s2multi;
t.drop();
t.ensureIndex({geo: "2dsphere"});
// Let's try the examples in the GeoJSON spec.
multiPointA = {
"type": "MultiPoint",
"coordinates": [[100.0, 0.0], [101.0, 1.0]]
};
assert.writeOK(t.insert({geo: multiPointA}));
multiLineStringA = {
"type": "MultiLineString",
"coordinates": [[[100.0, 0.0], [101.0, 1.0]], [[102.0, 2.0], [103.0, 3.0]]]
};
assert.writeOK(t.insert({geo: multiLineStringA}));
multiPolygonA = {
"type": "MultiPolygon",
"coordinates": [
[[[102.0, 2.0], [103.0, 2.0], [103.0, 3.0], [102.0, 3.0], [102.0, 2.0]]],
[
[[100.0, 0.0], [101.0, 0.0], [101.0, 1.0], [100.0, 1.0], [100.0, 0.0]],
[[100.2, 0.2], [100.8, 0.2], [100.8, 0.8], [100.2, 0.8], [100.2, 0.2]]
]
]
};
assert.writeOK(t.insert({geo: multiPolygonA}));
assert.eq(3, t.find({
geo: {$geoIntersects: {$geometry: {"type": "Point", "coordinates": [100, 0]}}}
}).itcount());
assert.eq(3,
t.find({
geo: {$geoIntersects: {$geometry: {"type": "Point", "coordinates": [101.0, 1.0]}}}
}).itcount());
// Inside the hole in multiPolygonA
assert.eq(
0, t.find({
geo: {$geoIntersects: {$geometry: {"type": "Point", "coordinates": [100.21, 0.21]}}}
}).itcount());
// One point inside the hole, one out.
assert.eq(
3,
t.find({
geo: {
$geoIntersects:
{$geometry: {"type": "MultiPoint", "coordinates": [[100, 0], [100.21, 0.21]]}}
}
}).itcount());
assert.eq(
3,
t.find({
geo: {
$geoIntersects: {
$geometry:
{"type": "MultiPoint", "coordinates": [[100, 0], [100.21, 0.21], [101, 1]]}
}
}
}).itcount());
// Polygon contains itself and the multipoint.
assert.eq(2, t.find({geo: {$geoWithin: {$geometry: multiPolygonA}}}).itcount());
partialPolygonA = {
"type": "Polygon",
"coordinates": [[[102.0, 2.0], [103.0, 2.0], [103.0, 3.0], [102.0, 3.0], [102.0, 2.0]]]
};
assert.writeOK(t.insert({geo: partialPolygonA}));
// Polygon contains itself, the partial poly, and the multipoint
assert.eq(3, t.find({geo: {$geoWithin: {$geometry: multiPolygonA}}}).itcount());
assert.eq(1, t.find({geo: {$geoWithin: {$geometry: partialPolygonA}}}).itcount());
// Itself, the multi poly, the multipoint...
assert.eq(3, t.find({geo: {$geoIntersects: {$geometry: partialPolygonA}}}).itcount());
|
