From b6749b81e09abad7944a4af21417e23cd8ad3460 Mon Sep 17 00:00:00 2001 From: mdroe Date: Mon, 4 Jun 2012 18:58:11 +0000 Subject: Use an entirely different method for area computation. It first transforms to a 2D plane by interpolating the great circle arcs through a Lambert azimuthal equal area projection. Then a standard 2D polygon area method is used. Make sure the polygons always go clockwise to aid in area computation. Fix union calculation -- it should be removing interior edges, not interior nodes. Fix some numerical out-of-range problems in great_circle_arc.py Remove the serial method for multi_union -- it no longer generates polygons that are compatible with calculating the area. Make debugging images more explanatory by putting a meaningful title at the top. git-svn-id: http://svn.stsci.edu/svn/ssb/stsci_python/stsci_python/branches/sphere@17200 fe389314-cf27-0410-b35b-8c050e845b92 Former-commit-id: fb841a481025150f631be6d64cf995dda592ecbc --- lib/test/test_intersection.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'lib/test/test_intersection.py') diff --git a/lib/test/test_intersection.py b/lib/test/test_intersection.py index adb79f0..12b3ce7 100644 --- a/lib/test/test_intersection.py +++ b/lib/test/test_intersection.py @@ -70,7 +70,7 @@ class intersection_test: lengths = np.array([len(x._points) for x in intersections]) assert np.all(lengths == [lengths[0]]) areas = np.array([x.area() for x in intersections]) - assert_array_almost_equal(areas, areas[0]) + assert_array_almost_equal(areas, areas[0], decimal=3) return run -- cgit