Changeset 868

Timestamp:

Feb 11, 2005, 3:38:42 PM (20 years ago)

Author:

ole

Message:

lat-lon to UTM Done!

Location:

inundation/ga/storm_surge/pyvolution/coordinate_transforms

Files:

• redfearn.py (modified) (4 diffs)
• redfearn.xls (modified) (previous)
• test_redfearn.py (modified) (2 diffs)

inundation/ga/storm_surge/pyvolution/coordinate_transforms/redfearn.py

@@ -1 +1 @@
 
 
-def degminsecs2decimal_degrees(dd,mm,ss):
+def degminsec2decimal_degrees(dd,mm,ss):
     assert abs(mm) == mm
     assert abs(ss) == ss    
@@ -12 +12 @@
     return sign * (abs(dd) + mm/60. + ss/3600.)      
     
+def decimal_degrees2degminsec(dec):
+
+    if dec < 0:
+        sign = -1
+    else:
+        sign = 1
+
+    dec = abs(dec)    
+    dd = int(dec)
+    f = dec-dd
+
+    mm = int(f*60)
+    ss = (f*60-mm)*60
+    
+    return sign*dd, mm, ss
 
 def redfearn(lat, lon):
@@ -32 +47 @@
     zone_width = 6              #Degrees
     false_easting = 500000
-    false_northing = 10000000  #In the southern hemisphere (0 in the northern)
+
+    if lat < 0:
+        false_northing = 10000000  #Southern hemisphere
+    else:
+        false_northing = 0         #Northern hemisphere)
     longitude_of_central_meridian_zone1 = -177  #Needed???
     longitude_of_central_meridian_zone0 = -183    
@@ -116 +135 @@
     #Northing
     term1 = nu*sinphi*cosphi*omega2/2  
-    term2 = nu*sinphi*cosphi3*omega4/24*(4*psi2+psi-t2)
-    term3 = nu*sinphi*cosphi5*omega6/720*\
+    term2 = nu*sinphi*cosphi3*(4*psi2+psi-t2)*omega4/24
+    term3 = nu*sinphi*cosphi5*\
             (8*psi4*(11-24*t2)-28*psi3*(1-6*t2)+\
-             psi2*(1-32*t2)-psi*2*t2+t4-t2)
-    term4 = nu*sinphi*cosphi7*omega8/40320*(1385-3111*t2+543*t4-t6)
-    sum = m + term1 + term2 + term3 + term4
-    northing = false_northing + K0*sum
+             psi2*(1-32*t2)-psi*2*t2+t4-t2)*omega6/720
+    term4 = nu*sinphi*cosphi7*(1385-3111*t2+543*t4-t6)*omega8/40320
+    northing = false_northing + K0*(m + term1 + term2 + term3 + term4)
 
-    print northing, zone
-
-    import sys; sys.exit() 
+    #Easting
+    term1 = nu*omega*cosphi
+    term2 = nu*cosphi3*(psi-t2)*omega3/6
+    term3 = nu*cosphi5*(4*psi3*(1-6*t2)+psi2*(1+8*t2)-2*psi*t2+t4)*omega5/120
+    term4 = nu*cosphi7*(61-479*t2+179*t4-t6)*omega7/5040
+    sum = term1 + term2 + term3 + term4
+    easting = false_easting + K0*sum
     
     return zone, easting, northing

inundation/ga/storm_surge/pyvolution/coordinate_transforms/test_redfearn.py

@@ -1 @@
-
-
-
-#Latitude:  -37 57 03.7203
-#Longitude: 144 25 29.5244
-#Zone:   55    
-#Easting:  273741.297  Northing: 5796489.777 
-#Latitude:   -37  57 ' 3.72030 ''  Longitude: 144  25 ' 29.52440 '' 
-#Grid Convergence:  -1  35 ' 3.65 ''  Point Scale: 1.00023056 
-
 
 from redfearn import *
 from Numeric import allclose
+
+#http://www.cellspark.com/UTM.html
+#http://www.ga.gov.au/nmd/geodesy/datums/redfearn_geo_to_grid.jsp
+
 
 #TEST 1
@@ -23 +17 @@
 #Grid Convergence:  1  47 ' 19.36 ''  Point Scale: 1.00042107 
 
-lat = degminsecs2decimal_degrees(-37,39,10.15610)
-lon = degminsecs2decimal_degrees(143,55,35.38390) 
+lat = degminsec2decimal_degrees(-37,39,10.15610)
+lon = degminsec2decimal_degrees(143,55,35.38390) 
+assert allclose(lat, -37.65282114)
+assert allclose(lon, 143.9264955)
+
+dd,mm,ss = decimal_degrees2degminsec(-37.65282114)
+assert dd==-37
+assert mm==39
+assert allclose(ss, 10.15610)
+
+dd,mm,ss = decimal_degrees2degminsec(143.9264955)
+assert dd==143
+assert mm==55
+assert allclose(ss, 35.38390) 
+
+#print dd,mm,ss
+
+zone, easting, northing = redfearn(lat,lon)
+
+assert zone == 54
+assert allclose(easting, 758173.797)
+assert allclose(northing, 5828674.340)
 
 
-assert allclose(lat, -37.65282114)
 
+#TEST 2
+
+#Latitude:  -37 57 03.7203
+#Longitude: 144 25 29.5244
+#Zone:   55    
+#Easting:  273741.297  Northing: 5796489.777 
+#Latitude:   -37  57 ' 3.72030 ''  Longitude: 144  25 ' 29.52440 '' 
+#Grid Convergence:  -1  35 ' 3.65 ''  Point Scale: 1.00023056 
+
+lat = degminsec2decimal_degrees(-37,57,03.7203)
+lon = degminsec2decimal_degrees(144,25,29.5244) 
+print lat, lon
+
+zone, easting, northing = redfearn(lat,lon)
+
+assert zone == 55
+assert allclose(easting, 273741.297)
+assert allclose(northing, 5796489.777)
+
+
+
+#Test 3
+lat = degminsec2decimal_degrees(-60,0,0)
+lon = degminsec2decimal_degrees(130,0,0) 
+
+zone, easting, northing = redfearn(lat,lon)
+#print zone, easting, northing
+
+assert zone == 52
+assert allclose(easting, 555776.267)
+assert allclose(northing, 3348167.264)
+
+
+#Test 4 (Kobenhavn)
+lat = 55.70248
+dd,mm,ss = decimal_degrees2degminsec(lat)
+print dd, mm, ss
+
+lon = 12.58364
+dd,mm,ss = decimal_degrees2degminsec(lon)
+print dd, mm, ss
 
 zone, easting, northing = redfearn(lat,lon)
 print zone, easting, northing
 
-assert zone == 54
-assert easting == 758173.797
-assert northing == 5828674.340
+assert zone == 33
+assert allclose(easting, 348157.631)
+#assert allclose(northing, 16175612.993) #With false_norting 10000000
+assert allclose(northing, 6175612.993) 
+
+