123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149 |
- // Boost.Geometry (aka GGL, Generic Geometry Library)
- // Unit Test
- // Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
- // Copyright (c) 2008-2012 Bruno Lalande, Paris, France.
- // Copyright (c) 2009-2012 Mateusz Loskot, London, UK.
- // Parts of Boost.Geometry are redesigned from Geodan's Geographic Library
- // (geolib/GGL), copyright (c) 1995-2010 Geodan, Amsterdam, the Netherlands.
- // Use, modification and distribution is subject to the Boost Software License,
- // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
- // http://www.boost.org/LICENSE_1_0.txt)
- #include <geometry_test_common.hpp>
- #include <boost/geometry/strategies/strategy_transform.hpp>
- #include <boost/geometry/algorithms/transform.hpp>
- #include <boost/geometry/geometries/point.hpp>
- #include <boost/geometry/geometries/point_xy.hpp>
- template <typename T, typename P>
- inline T check_distance(P const& p)
- {
- T x = bg::get<0>(p);
- T y = bg::get<1>(p);
- T z = bg::get<2>(p);
- return sqrt(x * x + y * y + z * z);
- }
- template <typename T>
- void test_transformations_spherical()
- {
- T const input_long = 15.0;
- T const input_lat = 5.0;
- T const expected_long = 0.26179938779914943653855361527329;
- T const expected_lat = 0.08726646259971647884618453842443;
- // Can be checked using http://www.calc3d.com/ejavascriptcoordcalc.html
- // (for phi use long, in radians, for theta use lat, in radians, they are listed there as "theta, phi")
- T const expected_polar_x = 0.084186;
- T const expected_polar_y = 0.0225576;
- T const expected_polar_z = 0.996195;
- // Can be checked with same URL using 90-theta for lat.
- // So for theta use 85 degrees, in radians: 0.08726646259971647884618453842443
- T const expected_equatorial_x = 0.962250;
- T const expected_equatorial_y = 0.257834;
- T const expected_equatorial_z = 0.0871557;
- // 1: Spherical-polar (lat=5, so it is near the pole - on a unit sphere)
- bg::model::point<T, 2, bg::cs::spherical<bg::degree> > sp(input_long, input_lat);
- // 1a: to radian
- bg::model::point<T, 2, bg::cs::spherical<bg::radian> > spr;
- bg::transform(sp, spr);
- BOOST_CHECK_CLOSE(bg::get<0>(spr), expected_long, 0.001);
- BOOST_CHECK_CLOSE(bg::get<1>(spr), expected_lat, 0.001);
- // 1b: to cartesian-3d
- bg::model::point<T, 3, bg::cs::cartesian> pc3;
- bg::transform(sp, pc3);
- BOOST_CHECK_CLOSE(bg::get<0>(pc3), expected_polar_x, 0.001);
- BOOST_CHECK_CLOSE(bg::get<1>(pc3), expected_polar_y, 0.001);
- BOOST_CHECK_CLOSE(bg::get<2>(pc3), expected_polar_z, 0.001);
- BOOST_CHECK_CLOSE(check_distance<T>(pc3), 1.0, 0.001);
- // 1c: back
- bg::transform(pc3, spr);
- BOOST_CHECK_CLOSE(bg::get<0>(spr), expected_long, 0.001);
- BOOST_CHECK_CLOSE(bg::get<1>(spr), expected_lat, 0.001);
- // 2: Spherical-equatorial (lat=5, so it is near the equator)
- bg::model::point<T, 2, bg::cs::spherical_equatorial<bg::degree> > se(input_long, input_lat);
- // 2a: to radian
- bg::model::point<T, 2, bg::cs::spherical_equatorial<bg::radian> > ser;
- bg::transform(se, ser);
- BOOST_CHECK_CLOSE(bg::get<0>(ser), expected_long, 0.001);
- BOOST_CHECK_CLOSE(bg::get<1>(ser), expected_lat, 0.001);
- bg::transform(se, pc3);
- BOOST_CHECK_CLOSE(bg::get<0>(pc3), expected_equatorial_x, 0.001);
- BOOST_CHECK_CLOSE(bg::get<1>(pc3), expected_equatorial_y, 0.001);
- BOOST_CHECK_CLOSE(bg::get<2>(pc3), expected_equatorial_z, 0.001);
- BOOST_CHECK_CLOSE(check_distance<T>(pc3), 1.0, 0.001);
- // 2c: back
- bg::transform(pc3, ser);
- BOOST_CHECK_CLOSE(bg::get<0>(spr), expected_long, 0.001); // expected_long
- BOOST_CHECK_CLOSE(bg::get<1>(spr), expected_lat, 0.001); // expected_lat
- // 3: Spherical-polar including radius
- bg::model::point<T, 3, bg::cs::spherical<bg::degree> > sp3(input_long, input_lat, 0.5);
- // 3a: to radian
- bg::model::point<T, 3, bg::cs::spherical<bg::radian> > spr3;
- bg::transform(sp3, spr3);
- BOOST_CHECK_CLOSE(bg::get<0>(spr3), expected_long, 0.001);
- BOOST_CHECK_CLOSE(bg::get<1>(spr3), expected_lat, 0.001);
- BOOST_CHECK_CLOSE(bg::get<2>(spr3), 0.5, 0.001);
- // 3b: to cartesian-3d
- bg::transform(sp3, pc3);
- BOOST_CHECK_CLOSE(bg::get<0>(pc3), expected_polar_x / 2.0, 0.001);
- BOOST_CHECK_CLOSE(bg::get<1>(pc3), expected_polar_y / 2.0, 0.001);
- BOOST_CHECK_CLOSE(bg::get<2>(pc3), expected_polar_z / 2.0, 0.001);
- BOOST_CHECK_CLOSE(check_distance<T>(pc3), 0.5, 0.001);
- // 3c: back
- bg::transform(pc3, spr3);
- BOOST_CHECK_CLOSE(bg::get<0>(spr3), expected_long, 0.001);
- BOOST_CHECK_CLOSE(bg::get<1>(spr3), expected_lat, 0.001);
- BOOST_CHECK_CLOSE(bg::get<2>(spr3), 0.5, 0.001);
- // 4: Spherical-equatorial including radius
- bg::model::point<T, 3, bg::cs::spherical_equatorial<bg::degree> > se3(input_long, input_lat, 0.5);
- // 4a: to radian
- bg::model::point<T, 3, bg::cs::spherical_equatorial<bg::radian> > ser3;
- bg::transform(se3, ser3);
- BOOST_CHECK_CLOSE(bg::get<0>(ser3), expected_long, 0.001);
- BOOST_CHECK_CLOSE(bg::get<1>(ser3), expected_lat, 0.001);
- BOOST_CHECK_CLOSE(bg::get<2>(ser3), 0.5, 0.001);
- // 4b: to cartesian-3d
- bg::transform(se3, pc3);
- BOOST_CHECK_CLOSE(bg::get<0>(pc3), expected_equatorial_x / 2.0, 0.001);
- BOOST_CHECK_CLOSE(bg::get<1>(pc3), expected_equatorial_y / 2.0, 0.001);
- BOOST_CHECK_CLOSE(bg::get<2>(pc3), expected_equatorial_z / 2.0, 0.001);
- BOOST_CHECK_CLOSE(check_distance<T>(pc3), 0.5, 0.001);
- // 4c: back
- bg::transform(pc3, ser3);
- BOOST_CHECK_CLOSE(bg::get<0>(ser3), expected_long, 0.001);
- BOOST_CHECK_CLOSE(bg::get<1>(ser3), expected_lat, 0.001);
- BOOST_CHECK_CLOSE(bg::get<2>(ser3), 0.5, 0.001);
- }
- int test_main(int, char* [])
- {
- test_transformations_spherical<double>();
- return 0;
- }
|