Changeset 6400

Timestamp:

Feb 24, 2009, 4:49:08 PM (15 years ago)

Author:

ole

Message:

Added data from Leharne for use with arbitrary meridian.

File:

• anuga_core/source/anuga/coordinate_transforms/test_redfearn.py (modified) (15 diffs)

anuga_core/source/anuga/coordinate_transforms/test_redfearn.py

@@ -122 +122 @@
 
     def test_UTM_6_nonstandard_projection(self):
-        #Test 6 (Geraldton, WA)
-
-        #First test native projection (zone 50)
+        """test_UTM_6_nonstandard_projection
+
+        Test that projections can be forced to
+        use other than native zone.
+
+        Data is from Geraldton, WA
+        """
+        
+
+        # First test native projection (zone 50)
         zone, easting, northing = redfearn(-29.233299999,114.05)
 
@@ -131 +138 @@
         assert num.allclose(northing, 6762559.15978)
 
-        #Testing using the native zone
+        # Testing using the native zone
         zone, easting, northing = redfearn(-29.233299999,114.05, zone=50)
 
@@ -138 +145 @@
         assert num.allclose(northing, 6762559.15978)
 
-        #Then project to zone 49
+        # Then project to zone 49
         zone, easting, northing = redfearn(-29.233299999,114.05,zone=49)
 
@@ -147 +154 @@
        
 
-        
-
-        #First test native projection (zone 49)
+        # First test native projection (zone 49)
         zone, easting, northing = redfearn(-29.1333,113.9667)
 
@@ -156 +161 @@
         assert num.allclose(northing, 6773605.46384)
 
-        #Then project to zone 50
+        # Then project to zone 50
         zone, easting, northing = redfearn(-29.1333,113.9667,zone=50)
 
@@ -163 +168 @@
         assert num.allclose(northing, 6773440.04726)
 
-        #Testing point on zone boundary
-        #First test native projection (zone 50)
+        # Testing point on zone boundary
+        # First test native projection (zone 50)
         zone, easting, northing = redfearn(-29.1667,114)
 
@@ -171 +176 @@
         assert num.allclose(northing, 6769820.01453)
 
-        #Then project to zone 49
+        # Then project to zone 49
         zone, easting, northing = redfearn(-29.1667,114,zone=49)
 
@@ -178 +183 @@
         assert num.allclose(northing, 6769820.01453)
 
-        #Testing furthest point in Geraldton scenario)
-        #First test native projection (zone 49)
+        # Testing furthest point in Geraldton scenario)
+        # First test native projection (zone 49)
         zone, easting, northing = redfearn(-28.2167,113.4167)
 
@@ -186 +191 @@
         assert num.allclose(northing, 6876426.38578)
 
-        #Then project to zone 50
+        # Then project to zone 50
         zone, easting, northing = redfearn(-28.2167,113.4167,zone=50)
 
@@ -193 +198 @@
         assert num.allclose(northing, 6873587.50926)
 
-        #Testing outside GDA zone (New Zeland)
-        #First test native projection (zone 60)
+        # Testing outside GDA zone (New Zeland)
+        # First test native projection (zone 60)
         zone, easting, northing = redfearn(-44,178)
 
@@ -201 +206 @@
         assert num.allclose(northing, 5127641.114461)
 
-        #Then project to zone 59
+        # Then project to zone 59
         zone, easting, northing = redfearn(-44,178,zone=59)
 
@@ -208 +213 @@
         assert num.allclose(northing, 5104249.395469)
 
-        #Then skip three zones 57 (native 60)
+        # Then skip three zones 57 (native 60)
         zone, easting, northing = redfearn(-44,178,zone=57)
 
@@ -223 +228 @@
         # Google Earth (he he)
         
-        #Testing outside GDA zone (Northern Hemisphere)
-        #First test native projection (zone 57)
+        # Testing outside GDA zone (Northern Hemisphere)
+        # First test native projection (zone 57)
         zone, easting, northing = redfearn(44,156)
 
@@ -237 +242 @@
         #assert num.allclose(northing, 14876249.1268)
         
-        #Then project to zone 56
+        # Then project to zone 56
         zone, easting, northing = redfearn(44,156,zone=56)
 
@@ -267 +272 @@
     #    #assert allclose(northing, 6181725.1724276)
 
+
+    def test_nonstandard_meridian(self):
+        """test_nonstandard_meridian
+
+        This test will verify that redfearn can be used to project
+        points using an arbitrary central meridian.
+        """
+
+
+        # FIXME: To do using csv file
+        pass
+
+
+        # The file projection_test_points.csv contains 10 points
+        # which straddle the boundary between UTM zones 53 and 54.
+        # They have been projected using a central meridian of 137.5
+        # degrees (the boundary is 138 so it is pretty much right
+        # in the middle of zones 53 and 54).
+
+        
+    
+    
+
     def test_convert_lats_longs(self):