Changeset 1012

Timestamp:

Mar 4, 2005, 6:08:57 PM (20 years ago)

Author:

duncan

Message:

increased functionality

File:

• inundation/ga/storm_surge/pyvolution/coordinate_transforms/geo_reference.py (modified) (1 diff)

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

@@ -8 +8 @@
     """
 
-    def __init__(self, xll_corner, yll_corner, zone):
-      self.xll_corner = xll_corner #Easting of lower left corner
-      self.yll_corner = yll_corner #Northing of lower left corner
-      self.zone = zone
+    def __init__(self,
+                 xllcorner=None,
+                 yllcorner=None,
+                 zone=None,
+                 NetCDFObject = None):
+        if not xllcorner == None:
+            self.xllcorner = xllcorner #Easting of lower left corner
+        if not yllcorner == None:
+            self.yllcorner = yllcorner #Northing of lower left corner
+        if not zone == None:
+            self.zone = zone
+        
+        if not NetCDFObject == None:
+            self.read_NetCDF(NetCDFObject)
+        
+    def write_NetCDF(self, outfile):
+        outfile.xllcorner = self.xllcorner 
+        outfile.yllcorner = self.yllcorner
+        outfile.zone = self.zone
 
-    def write_NetCDF(self, outfile):
-      outfile.xllcorner = self.xllcorner 
-      outfile.yllcorner = self.yllcorner
-      outfile.zone = self.zone
+        # might add this as well?
+        #outfile.false_easting = false_easting
+        #outfile.false_northing =false_northing
+    
+        #outfile.projection = projection
+        #outfile.datum = datum
+        #outfile.units = units
 
-      # might add this as well?
-      #outfile.false_easting = false_easting
-      #outfile.false_northing =false_northing
+    def read_NetCDF(self, infile):
+        self.xllcorner = infile.xllcorner
+        self.yllcorner = infile.yllcorner 
+        self.zone = infile.zone
+      
+        #false_easting = infile.false_easting[0]
+        #false_northing = infile.false_northing[0]
+
+    def __repr__(self):
+        return "(%f,%f, %s)" % (self.xllcorner,self.yllcorner, self.zone)
     
-      #outfile.projection = projection
-      #outfile.datum = datum
-      #outfile.units = units
-
-    def read_NetCDF(self, outfile):
-      self.xllcorner = infile.xllcorner[0] 
-      self.yllcorner = infile.yllcorner[0] 
-      self.zone = infile.zone[0]
-      
-      #false_easting = infile.false_easting[0]
-      #false_northing = infile.false_northing[0]
-
+    def __cmp__(self,other):
+        # FIXME (DSG) add a tolerence
+        cmp = 0
+        if not (self.xllcorner == self.xllcorner):
+            cmp = 1
+        if not (self.yllcorner == self.yllcorner):
+            cmp = 1
+        if not (self.zone == self.zone):
+            cmp = 1
+        return cmp
+        
 #-----------------------------------------------------------------------