Changeset 1966 for TI02-CSML/trunk/csml

Timestamp:

08/01/07 13:55:20 (14 years ago)

Author:

domlowe

Message:

CRS classes refined

File:

• TI02-CSML/trunk/csml/csmllibs/csmlcrs.py (modified) (5 diffs)

TI02-CSML/trunk/csml/csmllibs/csmlcrs.py

@@ +1 @@
+import string
 
 class CRSystem(object):
@@ -7 +8 @@
         self.lonAxis=None
         self.latAxis=None
-        self.srsDimension=len[axes]
-        
+        self.srsDimension=len(self.axes)
+        self.axisLabels=''
+        for axis in self.axes:
+            self.axisLabels=self.axisLabels + axis + ' '
 
 class CRSCatalogue(object):
@@ -16 +19 @@
         
         # define lon lat pressure time CRS:        
-        crs=CRSystem(srsName='ndg:crs:xypt', srsDimension='4', axes =['Long', 'Lat','Pressure','Time'])
+        crs=CRSystem(srsName='ndg:crs:xypt', axes =['Lon', 'Lat','Pressure','Time'])
         crs.lonAxis=0
         crs.latAxis=1
@@ -23 +26 @@
         
         # define lon lat height time CRS:        
-        crs=CRSystem(srsName='ndg:crs:xyht', srsDimension='4', axes =['Long', 'Lat','Height','Time'])
+        crs=CRSystem(srsName='ndg:crs:xyht', axes =['Lon', 'Lat','Height','Time'])
         crs.lonAxis=0
         crs.latAxis=1
@@ -29 +32 @@
         self.systems['ndg:crs:xyht']=crs
         
-        # define test CRS:   
-        crs=CRSystem(srsName='ndg:crs:abcde', srsDimension='5', axes =['Apple', 'Banana','Cat','Dog','Elephant'])
-        self.systems.['ndg:crs:abcde']=crs
+        # define lon lat time CRS:
+        crs=CRSystem(srsName='ndg:crs:xyt', axes =['Lon', 'Lat','Time'])
+        self.systems['ndg:crs:xyt']=crs
         
+        #define unknown CRS:
+        crs=CRSystem(srsName='ndg:crs:unknown', axes=['unknown'])
+        self.systems['ndg:crs:unknown']=crs
         
     def getCRS(self, axes, units):
-        #given any list of axis names and a list of units for these axes return the name of the CRS
-        pass
+        #given any list of axis names and a list of units for these axes attempt to determine the CRS
+        crs=self.systems['ndg:crs:unknown']
+        if len(axes)==3:
+            #it's a 3d crs
+            for axis in axes:
+                for unit in units:
+                    if string.lower(unit) in ['second', 'seconds', 's', 'mins','minute','minutes','hour','hours','h','hr','hrs','day','days']:
+                        #unit is time
+                        return self.systems['ndg:crs:xyt']
+        elif len(axes)==4:
+            #it's a 4d crs
+            pass       
+        return crs
         
 def main():
     cat=CRSCatalogue()
-    #test getCRS
-    print getCRS(axes=['longitude', 'latitude', 'pressure', 't'], units=['','','Pa','s'])
+    
+    #test getting various crs
+    print 'TEST 1:'
+    print '*************************************'
+    axs=['longitude', 'latitude', 'pressure', 't']
+    uns=['','','Pa','s']
+    print 'AXES: %s'%axs
+    print 'UNITS: %s'%uns
+    crs= cat.getCRS(axes=axs, units=uns)
+    print 'srsName = %s' %crs.srsName
+    print 'srsDimension = %s' %crs.srsDimension
+    print 'axisLabels =%s' %crs.axisLabels
+    print '**************************************'
+    
+    print 'TEST 2:'
+    print '*************************************'
+    axs=['longitude', 'latitude', 't']
+    uns=['','','s']  
+    print 'AXES: %s'%axs
+    print 'UNITS: %s'%uns
+    crs=cat.getCRS(axes=axs, units=uns)
+    print 'srsName = %s' %crs.srsName
+    print 'srsDimension = %s' %crs.srsDimension
+    print 'axisLabels = %s' %crs.axisLabels
 
 if __name__=="__main__":