@@ -715,6 +715,9 @@ def testGeodesicIntersectionAtAntimeridian(self):
715
715
716
716
# calculation should be ellipsoid dependent!
717
717
da .setEllipsoid ("Phobos2000" )
718
+ ####
719
+ #### THIS FAILS UNDER PROJ 6 -- we need to set the sourceCrs to a Phobos based CRS. (which is...?)
720
+ ####
718
721
lat , fract = da .latitudeGeodesicCrossesAntimeridian (QgsPointXY (- 175.15030911497356669 , 8.59851183021221033 ),
719
722
QgsPointXY (175.76717768974583578 ,
720
723
8.93749416467257873 ))
@@ -773,6 +776,9 @@ def testGeodesicLine(self):
773
776
774
777
# different ellipsoid, should be respected
775
778
da .setEllipsoid ("Phobos2000" )
779
+ ####
780
+ #### THIS FAILS UNDER PROJ 6 -- we need to set the sourceCrs to a Phobos based CRS. (which is...?)
781
+ ####
776
782
g = QgsGeometry .fromMultiPolylineXY (da .geodesicLine (QgsPointXY (121.4 , 6.4 ), QgsPointXY (- 121.6 , 2.8 ),
777
783
1000000 , False ))
778
784
self .assertEqual (g .asWkt (5 ),
@@ -850,6 +856,9 @@ def testSplitGeometryAtAntimeridian(self):
850
856
851
857
# different ellipsoid - should change intersection latitude
852
858
da .setEllipsoid ("Phobos2000" )
859
+ ####
860
+ #### THIS FAILS UNDER PROJ 6 -- we need to set the sourceCrs to a Phobos based CRS. (which is...?)
861
+ ####
853
862
g = da .splitGeometryAtAntimeridian (QgsGeometry .fromWkt ('LineString(179 10, -179 -20)' ))
854
863
self .assertEqual (g .asWkt (3 ), 'MultiLineString ((179 10, 180 -5.459),(-180 -5.459, -179 -20))' )
855
864
0 commit comments