Changeset 6424


Ignore:
Timestamp:
Feb 26, 2009, 10:43:04 PM (15 years ago)
Author:
ole
Message:

Experimented with scale factors in redfearn

Location:
anuga_core/source/anuga/coordinate_transforms
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • anuga_core/source/anuga/coordinate_transforms/redfearn.py

    r6404 r6424  
    3838
    3939def redfearn(lat, lon, false_easting=None, false_northing=None,
    40              zone=None, central_meridian=None):
     40             zone=None, central_meridian=None, scale_factor=None):
    4141    """Compute UTM projection using Redfearn's formula
    4242
    … …  
    6363    a = 6378137.0                       #Semi major axis
    6464    inverse_flattening = 298.257222101  #1/f
    65     K0 = 0.9996                         #Central scale factor   
     65    if scale_factor is None:
     66        K0 = 0.9996                         #Central scale factor   
     67    else:
     68        K0 = scale_factor
     69    #print 'scale', K0
    6670    zone_width = 6                      #Degrees
    6771

  • anuga_core/source/anuga/coordinate_transforms/test_redfearn.py

    r6404 r6424  
    300300
    301301            zone, easting, northing = redfearn(lat, lon,
    302                                                central_meridian=137.5)
     302                                               central_meridian=137.5,
     303                                               scale_factor=0.9998154)
    303304
    304305            print
    … …  
    318319        try:
    319320            zone, easting, northing = redfearn(lat, lon,
    320                                                zone=50, central_meridian=137.5)
     321                                               zone=50,
     322                                               central_meridian=137.5,
     323                                               scale_factor=0.968)
    321324        except:
    322325            pass
Note: See TracChangeset for help on using the changeset viewer.