Changeset 5663


Ignore:
Timestamp:
Aug 15, 2008, 3:31:46 PM (15 years ago)
Author:
ole
Message:

Got Northern hemisphere non-standard projection test to work using Google Earth Data

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/coordinate_transforms/test_redfearn.py

    r5662 r5663  
    213213        assert allclose(northing, 4949253.934967)
    214214
    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)
    232243
    233244       
Note: See TracChangeset for help on using the changeset viewer.