source: TI02-CSML/trunk/csml/Examples/Parsing/demoCreateV2.py @ 3296

import csml.parser
import csml.parser_extra
import csml.csmllibs.csmlextra

  #########################################
    ##  construction of CSML objects ##
  #########################################


# create an empty list to hold featureMembers
fms =[]


#### create a PointSeriesFeature: #####
ptsf=csml.parser.PointSeriesFeature()

#set the id attribute
ptsf.id='pointseriesf'

#set the description
#note, as description is an xml element not an attribute it is of type csString:
ptsf.description=csml.parser.csString('Station BLUEBIRD')

ptsf.location=csml.parser.csString('-50.5, 40.3')

#create a PointSeriesCoverage
ptscvg=csml.parser.PointSeriesCoverage()
ptscvg.id=csml.csmllibs.csmlextra.getRandomID()

#create a TimeSeriesDomain (appropriate domain for a PointSeriesFeature
ptsd=csml.parser.TimeSeries()
ptsd.id=csml.csmllibs.csmlextra.getRandomID()
#create (and populate) a TimePositionList. Using keyword arguements for conciseness.
ptsd.timePositionList=csml.parser.TimePositionList(frame='#pred20060427001',CONTENT='-18 -17 -16 -15 -14 -13 -12 -11 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60')
    
#create a RangeSet
rs=csml.parser.RangeSet()
#set the quantityList attribute of the RangeSet to be a MeasureOrNullList containing the values for the rangeSet
rs.quantityList=csml.parser.MeasureOrNullList(uom='MLUnits.xml#m',CONTENT='0.27 0.25 0.25 0.23 0.22 0.22 0.21 0.17 0.17 0.14 0.14 0.12 0.10 0.08 0.08 0.08 0.11 0.13 0.11 0.11 0.13 0.10 0.06 0.06 0.06 0.01 -0.03 -0.06 -0.09 -0.11 -0.11 -0.13 -0.16 -0.18 -0.17 -0.16 -0.18 -0.16 -0.12 -0.09 -0.08 -0.07 -0.06 -0.05 -0.04 -0.03 -0.03 -0.04 -0.02 -0.01 -0.02 -0.02 -0.02 -0.03 -0.03 -0.02 -0.01 -0.01 -0.02 -0.02 -0.03 -0.04 -0.04 -0.05 -0.04 -0.05 -0.08 -0.09 -0.11 -0.12 -0.12 -0.12 -0.12 -0.12 -0.13 -0.14 -0.15')


#Add the domain and rangeSet as attributes of the coverage
ptscvg.pointSeriesDomain=ptsd
ptscvg.rangeSet=rs

#set the coverage to be the 'value' of the feature
ptsf.value=ptscvg


#the parameter of the feature is of type Phenomenon, here href creates "xlink:href=..."
param=csml.parser.Phenomenon(href='CFStandardNames.xml#temperature')

#set the parameter of the feature to be the parameter just defined
ptsf.parameter=param

#append the feature to the list of featureMembers
fms.append(ptsf)



#create a profile series feature 
            #commented out, was just an approximate template
#psf=csml.parser.ProfileSeriesFeature()
#psf.id='profileseriesf'
#psf.description=csml.parser.csString('sample profile series feature')
#cvg = csml.parser.ProfileSeriesCoverage()
#cvg.id=csml.csmllibs.csmlextra.getRandomID()           
#psd=csml.parser.ProfileSeriesDomain()
#psd.id=csml.csmllibs.csmlextra.getRandomID()
#rs=csml.parser.RangeSet()
#gcT=csml.parser.GridCoordinatesTable()
#psd.coordTransformTable=gcT
#cvg.rangeSet=rs
#cvg.profileSeriesDomain=psd
#psf.value=cvg
#psf.parameter=param #defined earlier
#fms.append(psf)
           
                       
#create a profile feature
pf=csml.parser.ProfileFeature()            
pf.id='profilef'                        
pf.description=csml.parser.csString('sample profile feature')            
cvg = csml.parser.ProfileCoverage()
cvg.id=csml.csmllibs.csmlextra.getRandomID()
pd =csml.parser.ProfileDomain()
pd.id=csml.csmllibs.csmlextra.getRandomID()
pd.pointMembers=csml.parser.csString('1000, 950, 850, 700, 600, 500, 400, 300, 250, 200, 150, 100, 50, 30, 10')
cvg.profileDomain=pd
rs=csml.parser.RangeSet()
rs.quantityList=csml.parser.MeasureOrNullList(uom='MLUnits.xml#m',CONTENT='0.27 0.25 0.25 0.23 0.22 0.22 0.21 0.17 0.17 0.14 0.14 0.12 0.10 0.08 0.08')
cvg.rangeSet=rs
pf.value=cvg
pf.parameter=param
pf.location=csml.parser.csString('75 35')
pf.time=csml.parser.csString('2003-02-02T00:00:00.0')
fms.append(pf)            
                        
                        
#instantiate a FeatureCollection object (and add all featureMembers)
fc=csml.parser.CSMLFeatureCollection(featureMembers=fms)
fc.id=csml.csmllibs.csmlextra.getRandomID()

########### The Dataset  ##############
#Create an Empty Dataset
ds = csml.parser.Dataset()
#set the xml attribute 'id'
ds.id='Test001'

#set the description element 
ds.description=csml.parser.csString('This is a test Water level.Period 26/04/06 06h -> 29/04/06 12h')

#set the metaDataProperty(s)
#if there is only one metaDataProperty you can just do:
#ds.metaDataProperty=csml.parser.csString('Data (c) 2007 CCLRC // www.cclrc.ac.uk')

#but if there are several you need to put them in a list:
mdplist=[]
mdp1=csml.parser.csString('Data (c) 2007 CCLRC // www.cclrc.ac.uk')
mdp2=csml.parser.csString('Another CCLRC MetaDataProperty')
mdplist.append(mdp1)
mdplist.append(mdp2)
ds.metaDataProperty=mdplist

#now add the featureCollection to the dataset
ds.featureCollection=fc
                    
                                      
#dummy storage descriptor
sd=csml.parser.CSMLStorageDescriptor()
sd.id = csml.csmllibs.csmlextra.getRandomID()
ds.storageDescriptor=sd

                                      

#finished creating all the python objects!
########################################

####### Generating XML #################
            
#call the toXML method of the Dataset object:
csmldoc = ds.toXML()
#parse and pretty print the result
strCSML=csml.parser_extra.PrettyPrint(csmldoc)
strCSML=csml.parser_extra.removeInlineNS(strCSML)
print strCSML

#(you could save strCSML to a file at this point)
fname='ptseries.xml'
f=open(fname,'w')
f.write(strCSML)
print 'CSML doc written to %s'%fname