Changeset 1887 for TI02-CSML/branches/CSML2/csmllibs/csmlfeaturetypes.py

Timestamp:

18/12/06 12:15:46 (14 years ago)

Author:

domlowe

Message:

Initial changes to scanner to use csml 2 featuretypes. Very basic outline of GridSeriesFeature? working.

File:

• TI02-CSML/branches/CSML2/csmllibs/csmlfeaturetypes.py (modified) (4 diffs)

TI02-CSML/branches/CSML2/csmllibs/csmlfeaturetypes.py

@@ -23 +23 @@
         #empty list to hold featureMembers
         self.fms =[]
-        
+                
         #at the moment, only one featuretype per CSML Dataset is supported.
         #get the featuretype of the first representative file in the ffmap object
@@ -45 +45 @@
         try:
             descName=DI.getVariableAttribute('long_name')
+            descName 
         except AttributeError:
             descName = "missing name"
@@ -85 +86 @@
         return fe
         
+    
     def createCSMLGridSeriesFeatures(self):
         #This method assumes that the variables (features) are shared across identically structured files
         #should be supplied with a featurefilemap object (see csmlfiles for FileMapMaker)
+        print 'is this working at all?'
         representativeFiles=self.ffmap.getRepresentativeFiles()
         for repfile in representativeFiles:
@@ -105 +108 @@
                     #it's an axis or bounds not a feature, try next variable
                     continue
-                GridSeriesFeature_element=csml.parser.GridSeriesFeature()
-                GridSeriesFeature_element.id=str(allVarNames[i])
+                gsFeature=csml.parser.GridSeriesFeature()
+                gsFeature.id=str(allVarNames[i])
                 desc = self._getDescriptiveName(DI)
-                GridSeriesFeature_element.description=csml.parser.Description(desc)
+                #GridSeriesFeature_element.description=csml.parser.Description(desc)
+                gsFeature.description=desc
                 #DOMAIN
+                gsCoverage=csml.parser.GridSeriesCoverage()
                 gsDomain=csml.parser.GridSeriesDomain()
-                #DOMAIN REFERENCE
-                tpl=csml.parser.TimePositionList()
-                if self.timestorage =='inline':
-                    tpl.timePositions=self.timeString
-                    tpl.frame='%s:%s'%(self.caltype,self.units)
-                else:
-                    # do something to create a single extract for the times (from the representative file).
-                    tpl.timePositions = csml.csmllibs.csmlfileextracts.createSingleExtract(self.extractType, repfilename, self.timedim, len(self.timeString.split())) 
-                    tpl.frame='%s:%s'%(self.caltype,self.units)
-                gsDomain.domainReference=tpl
-                grid=csml.parser.Grid()
-                #COVERAGE FUNCTION
-                mr =csml.parser.MappingRule(csml.csmllibs.csmlextra.getMappingRule(len(dimNames)))
-                GridSeriesFeature_element.coverageFunction=mr 
-                #RANGESET
-                arrSz = DI.getArraySizeOfVar()
-                try:
-                    strUom = DI.getVariableAttribute('units')
-                except AttributeError:
-                    # if units attribute doesn't exist:
-                    strUom ="dimensionless or units not determined"
-                rs=csml.parser.RangeSet()
-                if self.valuestorage=='inline':
-                    #TO DO, store the rangeset inline - use Datablock class???
-                    pass 
-                else:
-                    #store the rangeSet as an aggregatedArray
-                    aa=csml.parser.AggregatedArray()
-                    aa.arraySize=[]
-                    aa.arraySize.append(arrSz)
-                    aa.uom=strUom
-                    aa.aggType='new' #can it be anything else?
-                    aa.aggIndex='1'
-                    #FileExtract (fe) element will be NetCDF/GRIB/PPExtract element (As defined earlier in ExtractType)
-                    self.extractType= DI.extractType
-                    fe = self._getCorrectExtractType()
-                    varSize=DI.getShapeOfVar()
-                    fe.arraySize=varSize
-                    fe.fileName=self.filesinDir
-                    fe.variableName=allVarNames[i]
-                    aa.component=[fe]
-                    rs.aggregatedArray=aa
-                GridSeriesFeature_element.rangeSet=rs
-                #DOMAIN COMPLEMENT
-                grid.srsName='urn:EPSG:GeographicCRS:4326'
-                numSpDims=len(varSize) -1 
-                grid.srsDimension=str(numSpDims)
-                grid.dimension=str(numSpDims)
-                ge =csml.parser.GridEnvelope(low=DI.getLowLimits(), high=DI.getHighLimits())
-                grid.limits=ge
-                #add an axisName element(s) for  each spatial dimension.
-                # and an ordinate element
-                axes=[]
-                for i in range (1, len(dimNames)):
-                    #axisNames
-                    axisname ='dim'+str(i)
-                    axes.append(axisname)
-                    #ordinates
-                    grid.ordinates=[]
-                    for i in range (1, len(dimNames)):
-                        ord=csml.parser.GridOrdinateDescription()
-                        ord.gridAxesSpanned='dim' + str(i)
-                        ord.sequenceRule=csml.csmllibs.csmlextra.getSeqRule(len(dimNames))
-                        dimName=dimNames[len(dimNames)-i]
-                        ord.definesAxis=dimName
-                        #look up file extract name in dictionary
-                        #(axisid stored in dictionary = current filename + variable name)
-                        axisid=self.repfilename+dimName
-                        if self.spatialstorage=='fileextract':
-                            #refer to extract
-                            ord.axisValues='#'+self.fileExtractDictionary[axisid]
-                        else:
-                            #store inline
-                             DI.setAxis(dimName)
-                             ord.axisValues=csml.csmllibs.csmlextra.cleanString(str(DI.getDataForAxis()))
-                        grid.ordinates.append(ord)
-                    grid.axisNames=axes
-                gsDomain.domainComplement=grid
-                GridSeriesFeature_element.domain=gsDomain
-                GridSeriesFeature_element.parameter=csml.parser.Phenomenon(href='http://badc.rl.ac.uk/localparams#%s'%allVarNames[i])
-                self.fms.append(GridSeriesFeature_element)
+                gsCoverage.gridSeriesDomain=gsDomain
+                gsFeature.value=gsCoverage
+                #gsDomain=csml.parser.GridSeriesDomain()
+                ##DOMAIN REFERENCE
+                #tpl=csml.parser.TimePositionList()
+                #if self.timestorage =='inline':
+                    #tpl.timePositions=self.timeString
+                    #tpl.frame='%s:%s'%(self.caltype,self.units)
+                #else:
+                    ## do something to create a single extract for the times (from the representative file).
+                    #tpl.timePositions = csml.csmllibs.csmlfileextracts.createSingleExtract(self.extractType, repfilename, self.timedim, len(self.timeString.split())) 
+                    #tpl.frame='%s:%s'%(self.caltype,self.units)
+                #gsDomain.domainReference=tpl
+                #grid=csml.parser.Grid()
+                ##COVERAGE FUNCTION
+                #mr =csml.parser.MappingRule(csml.csmllibs.csmlextra.getMappingRule(len(dimNames)))
+                #GridSeriesFeature_element.coverageFunction=mr 
+                ##RANGESET
+                #arrSz = DI.getArraySizeOfVar()
+                #try:
+                    #strUom = DI.getVariableAttribute('units')
+                #except AttributeError:
+                    ## if units attribute doesn't exist:
+                    #strUom ="dimensionless or units not determined"
+                #rs=csml.parser.RangeSet()
+                #if self.valuestorage=='inline':
+                    ##TO DO, store the rangeset inline - use Datablock class???
+                    #pass 
+                #else:
+                    ##store the rangeSet as an aggregatedArray
+                    #aa=csml.parser.AggregatedArray()
+                    #aa.arraySize=[]
+                    #aa.arraySize.append(arrSz)
+                    #aa.uom=strUom
+                    #aa.aggType='new' #can it be anything else?
+                    #aa.aggIndex='1'
+                    ##FileExtract (fe) element will be NetCDF/GRIB/PPExtract element (As defined earlier in ExtractType)
+                    #self.extractType= DI.extractType
+                    #fe = self._getCorrectExtractType()
+                    #varSize=DI.getShapeOfVar()
+                    #fe.arraySize=varSize
+                    #fe.fileName=self.filesinDir
+                    #fe.variableName=allVarNames[i]
+                    #aa.component=[fe]
+                    #rs.aggregatedArray=aa
+                #GridSeriesFeature_element.rangeSet=rs
+                ##DOMAIN COMPLEMENT
+                #grid.srsName='urn:EPSG:GeographicCRS:4326'
+                #numSpDims=len(varSize) -1 
+                #grid.srsDimension=str(numSpDims)
+                #grid.dimension=str(numSpDims)
+                #ge =csml.parser.GridEnvelope(low=DI.getLowLimits(), high=DI.getHighLimits())
+                #grid.limits=ge
+                ##add an axisName element(s) for  each spatial dimension.
+                ## and an ordinate element
+                #axes=[]
+                #for i in range (1, len(dimNames)):
+                    ##axisNames
+                    #axisname ='dim'+str(i)
+                    #axes.append(axisname)
+                    ##ordinates
+                    #grid.ordinates=[]
+                    #for i in range (1, len(dimNames)):
+                        #ord=csml.parser.GridOrdinateDescription()
+                        #ord.gridAxesSpanned='dim' + str(i)
+                        #ord.sequenceRule=csml.csmllibs.csmlextra.getSeqRule(len(dimNames))
+                        #dimName=dimNames[len(dimNames)-i]
+                        #ord.definesAxis=dimName
+                        ##look up file extract name in dictionary
+                        ##(axisid stored in dictionary = current filename + variable name)
+                        #axisid=self.repfilename+dimName
+                        #if self.spatialstorage=='fileextract':
+                            ##refer to extract
+                            #ord.axisValues='#'+self.fileExtractDictionary[axisid]
+                        #else:
+                            ##store inline
+                             #DI.setAxis(dimName)
+                             #ord.axisValues=csml.csmllibs.csmlextra.cleanString(str(DI.getDataForAxis()))
+                        #grid.ordinates.append(ord)
+                    #grid.axisNames=axes
+                #gsDomain.domainComplement=grid
+                #GridSeriesFeature_element.domain=gsDomain
+                #GridSeriesFeature_element.parameter=csml.parser.Phenomenon(href='http://badc.rl.ac.uk/localparams#%s'%allVarNames[i])
+                self.fms.append(gsFeature)
             DI.closeFile()
     ###End of createCSMLGridSeriesFeatures###