boost_mirror/geometry
1
0
Fork 0
mirror of https://github.com/boostorg/geometry.git synced 2025-05-12 14:01:47 +00:00
Code Issues Packages Projects Releases Wiki Activity

[test][projections] Add esri/iau2000 tests and change names of aeqd projection types.

Browse Source
This commit is contained in:
Adam Wulkiewicz 2017-10-29 03:22:13 +01:00
parent 595d6736d9
commit fe3ae4a41d
2 changed files with 11 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
test/srs/projection_interface.cpp
View File

@ -15,6 +15,7 @@
#include <boost/geometry/geometries/geometries.hpp> #include <boost/geometry/geometries/geometries.hpp>
#include <boost/geometry/srs/epsg.hpp> #include <boost/geometry/srs/epsg.hpp>
#include <boost/geometry/srs/esri.hpp> #include <boost/geometry/srs/esri.hpp>
#include <boost/geometry/srs/iau2000.hpp>
#include <boost/geometry/srs/projection.hpp> #include <boost/geometry/srs/projection.hpp>
#include "check_geometry.hpp" #include "check_geometry.hpp"
@ -28,7 +29,7 @@ int test_main(int, char*[])
typedef point<double, 2, cs::geographic<degree> > point_ll; typedef point<double, 2, cs::geographic<degree> > point_ll;
typedef point<double, 2, cs::cartesian> point_xy; typedef point<double, 2, cs::cartesian> point_xy;
{ {
point_ll pt_ll(1, 1); point_ll pt_ll(1, 1);
point_ll pt_ll2(0, 0); point_ll pt_ll2(0, 0);
@ -49,13 +50,16 @@ int test_main(int, char*[])
point_xy pt_xy(0, 0); point_xy pt_xy(0, 0);
projection<> prj = epsg(2000); projection<> prj = epsg(2000);
prj.forward(pt_ll, pt_xy); prj.forward(pt_ll, pt_xy);
test::check_geometry(pt_xy, "POINT(9413505.3284665551 237337.74515944949)", 0.001); test::check_geometry(pt_xy, "POINT(9413505.3284665551 237337.74515944949)", 0.001);
prj.inverse(pt_xy, pt_ll2); prj.inverse(pt_xy, pt_ll2);
// TODO: investigate this wierd result // TODO: investigate this wierd result
test::check_geometry(pt_ll2, "POINT(-2.4463131191981073 1.5066638962045082)", 0.001); test::check_geometry(pt_ll2, "POINT(-2.4463131191981073 1.5066638962045082)", 0.001);
projection<> prj2 = esri(37001);
projection<> prj3 = iau2000(19900);
} }
{ {
@ -87,6 +91,9 @@ int test_main(int, char*[])
prj.inverse(pt_xy, pt_ll2); prj.inverse(pt_xy, pt_ll2);
// TODO: investigate this wierd result // TODO: investigate this wierd result
test::check_geometry(pt_ll2, "POINT(-2.4463131191981073 1.5066638962045082)", 0.001); test::check_geometry(pt_ll2, "POINT(-2.4463131191981073 1.5066638962045082)", 0.001);
projection<static_esri<37001> > prj2;
projection<static_iau2000<19900> > prj3;
} }
{ {

4
test/srs/projections_static.cpp
View File

@ -84,8 +84,8 @@ void test_all()
geo_point_type utrecht = bg::make<geo_point_type>(5.1213, 52.0907); geo_point_type utrecht = bg::make<geo_point_type>(5.1213, 52.0907);
test_forward<bg::projections::aea_ellipsoid>(amsterdam, utrecht, "+ellps=WGS84 +units=m +lat_1=55 +lat_2=65"); test_forward<bg::projections::aea_ellipsoid>(amsterdam, utrecht, "+ellps=WGS84 +units=m +lat_1=55 +lat_2=65");
test_forward<bg::projections::aeqd_ellipsoid>(amsterdam, utrecht, "+ellps=WGS84 +units=m"); test_forward<bg::projections::aeqd_e>(amsterdam, utrecht, "+ellps=WGS84 +units=m");
test_forward<bg::projections::aeqd_spheroid>(amsterdam, utrecht, "+ellps=WGS84 +units=m"); test_forward<bg::projections::aeqd_s>(amsterdam, utrecht, "+ellps=WGS84 +units=m");
test_forward<bg::projections::airy_spheroid>(amsterdam, utrecht, "+ellps=WGS84 +units=m", 4); test_forward<bg::projections::airy_spheroid>(amsterdam, utrecht, "+ellps=WGS84 +units=m", 4);
test_forward<bg::projections::aitoff_spheroid>(amsterdam, utrecht, "+ellps=WGS84 +units=m", 2); test_forward<bg::projections::aitoff_spheroid>(amsterdam, utrecht, "+ellps=WGS84 +units=m", 2);