Changeset 1609 for TI02-CSML

Timestamp:

19/10/06 14:19:40 (13 years ago)

Author:

domlowe

Message:

coverage aggregation and basic data granule creation working

Location:

TI02-CSML/trunk/csml2MolesStuff

Files:

• csml2moles.py (modified) (4 diffs)
• molesWriter.py (modified) (3 diffs)

TI02-CSML/trunk/csml2MolesStuff/csml2moles.py

@@ -6 +6 @@
 import elementtree.ElementTree as etree
 import datetime
-import molesWriter
+import molesWriter as MW
 #this is a fix to the  ElementTree namespace problem that namespaces are usually represented as ns0, ns1, ns2 etc.
 etree._namespace_map.update({
@@ -84 +84 @@
         return self.envelope
 
-    #new process:
+
 def main(args=None):
     '''Get command line arguments, should be paths to csml files
@@ -99 +99 @@
     aggregator=EnvelopeAggregator(env)
  
+    
+    #set up lists to hold stuff
+    schemeIdentifiers=[]
+    repositoryIdentifiers=[]
+    localIdentifiers=[]
+    
     #for each file:
     for f in csmlfilelist:
@@ -106 +112 @@
         aggregator.compareEnvelope(_getEnvelopefromDataset(dataset))
         #get other things:
-        
+        schemeIdentifiers.append('NDG-A0')
+        repositoryIdentifiers.append('badc.nerc.ac.uk')
+        localIdentifiers.append(f)
     #get aggregated envelope:
     finalEnvelope=aggregator.getAggregatedEnvelope()
+    limitN=finalEnvelope.upperCorner.vals.split()[0]
+    limitS=finalEnvelope.lowerCorner.vals.split()[0]
+    limitE=finalEnvelope.upperCorner.vals.split()[1]
+    limitW=finalEnvelope.lowerCorner.vals.split()[1]
     #get other aggregations:
     
     #write MOLES skeleton document:
+    M=MW.MolesDoc()
+    M._createClasses()
+    
+    ''' The MOLES document:
+    create your moles doc by setting attributes of "M.className"
+    the classNames used must be declared in molesWriter.MolesDoc.classList'''
+    
+    #MetaData ID (need to decide on how to define this)
+    dgMID=M.dgMetadataID(schemeIdentifier='NDG-B0', repositoryIdentifier='badc.nerc.ac.uk', localIdentifier='COAPEC_HadCM3_500Yr')
+    
+    #create data granules
+    datagranules=[]
+    for i, file in enumerate(csmlfilelist):
+        DMid=M.dataModelID(schemeIdentifier=schemeIdentifiers[i], repositoryIdentifier=repositoryIdentifiers[i],localIdentifier= localIdentifiers[i])
+        DG  = M.dgDataGranule(dataModelID=DMid)
+        datagranules.append(DG)
+    
+    #create Data Entity    
+    dgDE= M.dgDataEntity(dgDataGranule=datagranules)
+    
+    #create coverage
+    dgBB=M.dgBoundingBox(LimitNorth = limitN, LimitSouth=limitS,LimitEast=limitE, LimitWest=limitW)
+    dgSc=M.dgSpatialCoverage(dgBoundingBox=dgBB)
+    dgTc=M.dgTemporalCoverage()
+    dgST=M.dgSpatioTemporalCoverage(dgSpatialCoverage = dgSc, dgTemporalCoverage=dgTc)
+    dgCv=M.dgCoverage(dgSpatioTemporalCoverage=dgST)
+    
+    #create metadata record
+    dgMR=M.dgMetadataRecord(dgMetadataID=dgMID, dgDataEntity=dgDE, dgCoverage=dgCv)
+    dgMeta=MW.dgMetadata(dgMetadataRecord=dgMR)
+
+    #produce XML 
+    molestree=dgMeta.toXML()
+    #print it out.
+    moles=csml.parser_extra.PrettyPrint(molestree)
+    print moles
+    
+    
     print finalEnvelope.lowerCorner.vals
     print finalEnvelope.upperCorner.vals
     print finalEnvelope.timePosition
     print finalEnvelope.timePosition2
-    
+   
 if __name__=='__main__':
     main()

TI02-CSML/trunk/csml2MolesStuff/molesWriter.py

@@ -10 +10 @@
     def toXML(self,molesFrag):
         for attr in self.__dict__:
+            print attr
             if isinstance(self.__dict__[attr], molesElement):
                 subelem = self.__dict__[attr]
@@ -19 +20 @@
                         frag.append(subfrag)
                     else:
-                        frag.text=attrib
+                        subfrag=ET.Element(attrib)
+                        subfrag.text=subelem.__dict__[attrib]
+                        frag.append(subfrag)
                 molesFrag.append(frag)
             elif isinstance(self.__dict__[attr], list):
@@ -64 +67 @@
     def _createClasses(self):
         #if you want more classes just add their names to this list
-        classList=[\
-        'dataModelID', \
+        classList= \
+        ['dataModelID', \
         'dgDataGranule', \
         'dgDataEntity', \
         'dgMetadataRecord', \
-        'dgMetadataID', \
-        ]
+        'dgMetadataID', 
+        'dgCoverage', \
+        'dgSpatioTemporalCoverage', \
+        'dgSpatialCoverage', \
+        'dgTemporalCoverage', \
+        'dgBoundingBox', \
+        'DateRange']
 
         for className in classList: