Changeset 5663
- Timestamp:
- Aug 15, 2008, 3:31:46 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/coordinate_transforms/test_redfearn.py
r5662 r5663 213 213 assert allclose(northing, 4949253.934967) 214 214 215 #Note Projecting into the Northern Hemisphere Does not coincide 216 #redfearn or ArcMap conversions 217 ## #Testing outside GDA zone (Northern Hemisphere) 218 ## #First test native projection (zone 57) 219 ## zone, easting, northing = redfearn(44,156) 220 ## 221 ## assert zone == 57 222 ## assert allclose(easting, 259473.678944) 223 ## assert allclose(northing, 14876249.1268) 224 ## 225 ## #Then project to zone 59 226 ## zone, easting, northing = redfearn(44,156,zone=56) 227 ## 228 ## assert zone == 56 229 ## assert allclose(easting, 740526.321055) 230 ## assert allclose(northing, 14876249.1268) 231 215 # Note Projecting into the Northern Hemisphere Does not coincide 216 # redfearn or ArcMap conversions. The difference lies in 217 # False Northing which is 0 for the Northern hemisphere 218 # in the redfearn implementation but 10 million in Arc. 219 # 220 # But the redfearn implementation does coincide with 221 # Google Earth (he he) 222 223 #Testing outside GDA zone (Northern Hemisphere) 224 #First test native projection (zone 57) 225 zone, easting, northing = redfearn(44,156) 226 227 # Google Earth interpretation 228 assert zone == 57 229 assert allclose(easting, 259473.69) 230 assert allclose(northing, 4876249.13) 231 232 # ArcMap's interpretation 233 #assert zone == 57 234 #assert allclose(easting, 259473.678944) 235 #assert allclose(northing, 14876249.1268) 236 237 #Then project to zone 56 238 zone, easting, northing = redfearn(44,156,zone=56) 239 240 assert zone == 56 241 assert allclose(easting, 740526.321055) 242 assert allclose(northing, 4876249.13) 232 243 233 244
Note: See TracChangeset
for help on using the changeset viewer.