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

'''apicalls.py  - test the api'''

import csml

f='example.xml'
outputdir = 'out'

#Initialise the CSML container
csmlContainer=csml.API.csmlContainer.Container(outputdir,csmlname='myfile.xml',docID="mycsmldoc_1", metadata="mymetadatastring")

#Initialise and parse the dataset
csmldoc = csml.parser.Dataset(file=f)


#get list of features in the dataset
flist= csmldoc.getFeatureList() 
print '\n Here are all the features in %s:' %f
print flist

#select a feature by name (gml:id)
print '\n Selecting feature with gml:id = %s' %flist[0]
feature=csmldoc.getFeature(flist[0])

f#eature 'bestname'
print feature.id
print feature.description

#get the domain of the feature
print '\n The feature has domain: '
print feature.getDomain()

#get list of allowed subsettings
print '\n the following feature subsetting operations are allowed:'
print feature.getAllowedSubsettings()


#define a selection (based on the domain)
timeSelection=['2007-08-12T10:45:00', '2007-08-12T11:45:00']

#request subsetted data from features (can set paths here) and add to container
#subset a feature
csmlContainer.add(feature.subsetToPointSeries(times=timeSelection,  ncname='myfile.nc'))

#when you have finished, get the contents of the container
csmlbundle=csmlContainer.getContents()
print 'CSML %s'%csmlbundle[0] # csml document
paths=csmlbundle[1:]
print paths #paths to netcdf files