COMPARE(sg_bucket_span(-89.9), 12.0);
}
+void testBasic()
+{
+ SGBucket b1(5.1, 55.05);
+ COMPARE(b1.get_chunk_lon(), 5);
+ COMPARE(b1.get_chunk_lat(), 55);
+ COMPARE(b1.get_x(), 0);
+ COMPARE(b1.get_y(), 0);
+ COMPARE(b1.gen_index(), 3040320);
+ COMPARE(b1.gen_base_path(), "e000n50/e005n55");
+ VERIFY(b1.isValid());
+
+ SGBucket b2(-10.1, -43.8);
+ COMPARE(b2.get_chunk_lon(), -11);
+ COMPARE(b2.get_chunk_lat(), -44);
+ COMPARE(b2.get_x(), 3);
+ COMPARE(b2.get_y(), 1); // latitude chunks numbered bottom to top, it seems
+ COMPARE(b2.gen_base_path(), "w020s50/w011s44");
+ VERIFY(b2.isValid());
+
+ SGBucket b3(123.48, 9.01);
+ COMPARE(b3.get_chunk_lon(), 123);
+ COMPARE(b3.get_chunk_lat(), 9);
+ COMPARE(b3.get_x(), 3);
+ COMPARE(b3.get_y(), 0);
+ COMPARE(b3.gen_base_path(), "e120n00/e123n09");
+ VERIFY(b3.isValid());
+
+ SGBucket defBuck;
+ VERIFY(!defBuck.isValid());
+
+ b3.make_bad();
+ VERIFY(!b3.isValid());
+
+ SGBucket atAntiMeridian(180.0, 12.3);
+ VERIFY(atAntiMeridian.isValid());
+ COMPARE(atAntiMeridian.get_chunk_lon(), -180);
+ COMPARE(atAntiMeridian.get_x(), 0);
+
+ SGBucket atAntiMeridian2(-180.0, -78.1);
+ VERIFY(atAntiMeridian2.isValid());
+ COMPARE(atAntiMeridian2.get_chunk_lon(), -180);
+ COMPARE(atAntiMeridian2.get_x(), 0);
+
+// check comparisom operator overload
+ SGBucket b4(5.11, 55.1);
+ VERIFY(b1 == b4); // should be equal
+ VERIFY(b1 == b1);
+ VERIFY(b1 != defBuck);
+ VERIFY(b1 != b2);
+
+// check wrapping/clipping of inputs
+ SGBucket wrapMeridian(-200.0, 45.0);
+ COMPARE(wrapMeridian.get_chunk_lon(), 160);
+
+ SGBucket clipPole(48.9, 91);
+ COMPARE(clipPole.get_chunk_lat(), 89);
+}
+
void testPolar()
{
SGBucket b1(0.0, 89.92);
COMPARE(b1.get_x(), 0);
COMPARE(b1.get_y(), 7);
+ COMPARE_EP(b1.get_highest_lat(), 90.0);
+ COMPARE_EP(b1.get_width_m(), 10.0);
+
+ COMPARE(b2.get_chunk_lat(), 89);
+ COMPARE(b2.get_chunk_lon(), 0);
+ COMPARE(b2.get_x(), 0);
+ COMPARE(b2.get_y(), 7);
+
COMPARE(b1.gen_index(), b2.gen_index());
+ SGGeod actualNorthPole1 = b1.get_corner(2);
+ SGGeod actualNorthPole2 = b1.get_corner(3);
+ COMPARE_EP(actualNorthPole1.getLatitudeDeg(), 90.0);
+ COMPARE_EP(actualNorthPole1.getLongitudeDeg(), 12.0);
+ COMPARE_EP(actualNorthPole2.getLatitudeDeg(), 90.0);
+ COMPARE_EP(actualNorthPole2.getLongitudeDeg(), 0.0);
+
SGBucket b3(-2, 89.88);
SGBucket b4(-7, 89.88);
COMPARE(b3.gen_index(), b4.gen_index());
COMPARE(b5.get_x(), 0);
COMPARE(b5.get_y(), 0);
COMPARE(b5.gen_index(), b6.gen_index());
+ COMPARE_EP(b5.get_highest_lat(), -90.0);
+ COMPARE_EP(b5.get_width_m(), 10.0);
+
+ SGGeod actualSouthPole1 = b5.get_corner(0);
+ SGGeod actualSouthPole2 = b5.get_corner(1);
+ COMPARE_EP(actualSouthPole1.getLatitudeDeg(), -90.0);
+ COMPARE_EP(actualSouthPole1.getLongitudeDeg(), -180);
+ COMPARE_EP(actualSouthPole2.getLatitudeDeg(), -90.0);
+ COMPARE_EP(actualSouthPole2.getLongitudeDeg(), -168);
- // no automatic wrapping of these values occurs
SGBucket b7(200, 89.88);
- COMPARE(b7.get_chunk_lon(), 192);
+ COMPARE(b7.get_chunk_lon(), -168);
}
+// test the tiles just below the pole (between 86 & 89 degrees N/S)
void testNearPolar()
{
SGBucket b1(1, 88.5);
COMPARE(b1.get_chunk_lon(), 0);
COMPARE(b1.get_chunk_lat(), 88);
VERIFY(b1.gen_index() != b2.gen_index());
-
+
SGBucket b3(176.1, 88.5);
COMPARE(b3.get_chunk_lon(), 176);
SGBucket b4(-178, 88.5);
- // the fix for the polar cap issue, causes this to report -180,
- // not -184
COMPARE(b4.get_chunk_lon(), -180);
}
void testOffset()
{
+ // bucket just below the 22 degree cutoff, so the next tile north
+ // is twice the width
+ SGBucket b1(-59.8, 21.9);
+ COMPARE(b1.get_chunk_lat(), 21);
+ COMPARE(b1.get_chunk_lon(), -60);
+ COMPARE(b1.get_x(), 1);
+ COMPARE(b1.get_y(), 7);
+ // offset vertically
+ SGBucket b2(b1.sibling(0, 1));
+ COMPARE(b2.get_chunk_lat(), 22);
+ COMPARE(b2.get_chunk_lon(), -60);
+ COMPARE(b2.get_x(), 0);
+ COMPARE(b2.get_y(), 0);
+
+ COMPARE(b2.gen_index(), sgBucketOffset(-59.8, 21.9, 0, 1));
+
+ // offset vertically and horizontally. We compute horizontal (x)
+ // movement at the target latitude, so this should move 0.25 * -3 degrees,
+ // NOT 0.125 * -3 degrees.
+ SGBucket b3(b1.sibling(-3, 1));
+ COMPARE(b3.get_chunk_lat(), 22);
+ COMPARE(b3.get_chunk_lon(), -61);
+ COMPARE(b3.get_x(), 1);
+ COMPARE(b3.get_y(), 0);
+
+ COMPARE(b3.gen_index(), sgBucketOffset(-59.8, 21.9, -3, 1));
}
-int main(int argc, char* argv[])
+void testPolarOffset()
{
- testBucketSpans();
-
- SGBucket b1(5.1, 55.05);
- COMPARE(b1.get_chunk_lon(), 5);
- COMPARE(b1.get_chunk_lat(), 55);
+ SGBucket b1(-11.7, -89.6);
+ COMPARE(b1.get_chunk_lat(), -90);
+ COMPARE(b1.get_chunk_lon(), -12);
COMPARE(b1.get_x(), 0);
- COMPARE(b1.get_y(), 0);
- COMPARE(b1.gen_index(), 3040320);
- COMPARE(b1.gen_base_path(), "e000n50/e005n55");
+ COMPARE(b1.get_y(), 3);
- SGBucket b2(-10.1, -43.8);
- COMPARE(b2.get_chunk_lon(), -11);
- COMPARE(b2.get_chunk_lat(), -44);
- COMPARE(b2.get_x(), 3);
- COMPARE(b2.get_y(), 1); // latitude chunks numbered bottom to top, it seems
- COMPARE(b2.gen_base_path(), "w020s50/w011s44");
+ // offset horizontally
+ SGBucket b2(b1.sibling(-2, 0));
+ COMPARE(b2.get_chunk_lat(), -90);
+ COMPARE(b2.get_chunk_lon(), -36);
+ COMPARE(b2.get_x(), 0);
+ COMPARE(b2.get_y(), 3);
- SGBucket b3(123.48, 9.01);
- COMPARE(b3.get_chunk_lon(), 123);
- COMPARE(b3.get_chunk_lat(), 9);
- COMPARE(b3.get_x(), 3);
- COMPARE(b3.get_y(), 0);
- COMPARE(b3.gen_base_path(), "e120n00/e123n09");
+ COMPARE(b2.gen_index(), sgBucketOffset(-11.7, -89.6, -2, 0));
+
+// offset and wrap
+ SGBucket b3(-170, 89.1);
+ SGBucket b4(b3.sibling(-1, 0));
+ COMPARE(b4.get_chunk_lat(), 89);
+ COMPARE(b4.get_chunk_lon(), 168);
+ COMPARE(b4.get_x(), 0);
+ COMPARE(b4.get_y(), 0);
+
+ COMPARE(b4.gen_index(), sgBucketOffset(-170, 89.1, -1, 0));
+
+
+ SGBucket b5(177, 87.3);
+ SGBucket b6(b5.sibling(1, 1));
+ COMPARE(b6.get_chunk_lat(), 87);
+ COMPARE(b6.get_chunk_lon(), -180);
+ COMPARE(b6.get_x(), 0);
+ COMPARE(b6.get_y(), 3);
+
+ COMPARE(b6.gen_index(), sgBucketOffset(177, 87.3, 1, 1));
+
+ // offset vertically towards the pole
+ SGBucket b7(b1.sibling(0, -5));
+ VERIFY(!b7.isValid());
+
+ VERIFY(!SGBucket(0, 90).sibling(0, 1).isValid());
+}
+
+// test behaviour of bucket-offset near the anti-meridian (180-meridian)
+void testOffsetWrap()
+{
+ // near the equator
+ SGBucket b1(-179.8, 16.8);
+ COMPARE(b1.get_chunk_lat(), 16);
+ COMPARE(b1.get_chunk_lon(), -180);
+ COMPARE(b1.get_x(), 1);
+ COMPARE(b1.get_y(), 6);
+
+ SGBucket b2(b1.sibling(-2, 0));
+ COMPARE(b2.get_chunk_lat(), 16);
+ COMPARE(b2.get_chunk_lon(), 179);
+ COMPARE(b2.get_x(), 7);
+ COMPARE(b2.get_y(), 6);
+ COMPARE(b2.gen_index(), sgBucketOffset(-179.8, 16.8, -2, 0));
+
+
+}
+
+int main(int argc, char* argv[])
+{
+ testBucketSpans();
+ testBasic();
testPolar();
testNearPolar();
testOffset();
+ testOffsetWrap();
+ testPolarOffset();
cout << "all tests passed OK" << endl;
return 0; // passed