Changeset 27524
- Timestamp:
- 01/19/23 04:20:45 (2 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/trunk-jpl/src/m/coordsystems/ll2xy.py ΒΆ
r26401 r27524 13 13 x, y = ll2xy(lat, lon, sgn, central_meridian, standard_parallel) 14 14 15 - sgn = Sign of latitude 16 15 - sgn = Sign of latitude 1 : north latitude (default is mer = 45 lat = 70) 16 -1 : south latitude (default is mer = 0 lat = 71) 17 17 """ 18 18 assert sgn == 1 or sgn == -1, 'error: sgn should be either 1 or -1' 19 20 19 20 # Get central_meridian and standard_parallel depending on hemisphere 21 21 if len(args) == 2: 22 delta = args[0]23 slat = args[1]22 delta = args[0] 23 slat = args[1] 24 24 elif len(args) == 0: 25 if sgn == 1:26 delta = 45.27 slat = 70.28 print(' ll2xy: creating coordinates in north polar stereographic (Std Latitude: 70degN Meridian: 45deg)')29 elif sgn == -1:30 delta = 0.31 slat = 71.32 print(' ll2xy: creating coordinates in south polar stereographic (Std Latitude: 71degS Meridian: 0deg)')33 else:34 raise ValueError('sgn should be either 1 or -1')25 if sgn == 1: 26 delta = 45. 27 slat = 70. 28 print(' ll2xy: creating coordinates in north polar stereographic (Std Latitude: 70degN Meridian: 45deg)') 29 elif sgn == -1: 30 delta = 0. 31 slat = 71. 32 print(' ll2xy: creating coordinates in south polar stereographic (Std Latitude: 71degS Meridian: 0deg)') 33 else: 34 raise ValueError('sgn should be either 1 or -1') 35 35 else: 36 raise Exception('bad usage: type "help(ll2xy)" for details')36 raise Exception('bad usage: type "help(ll2xy)" for details') 37 37 38 38 # if lat, lon passed as lists, convert to np.arrays 39 39 if type(lat) != "np.ndarray": 40 lat = np.array(lat)40 lat = np.array(lat) 41 41 if type(lon) != "np.ndarray": 42 lon = np.array(lon)42 lon = np.array(lon) 43 43 44 44 # Conversion constant from degrees to radians
Note:
See TracChangeset
for help on using the changeset viewer.