source: TI02-CSML/trunk/Examples/parsing/apicalls.py @ 1493

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

import csml

f='test.xml'

#Initialise and parse the dataset
csmldoc = csml.parser.Dataset()
csmldoc.parse(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[4]
feature=csmldoc.getFeature(flist[4])

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

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

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

#get combined domain
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 ref/complement)
#timeSelection=['2881-3-16T0:0:0.0', '2881-4-16T0:0:0.0', '2881-5-16T0:0:0.0', '2881-6-16T0:0:0.0', '2881-7-16T0:0:0.0', '2881-8-16T0:0:0.0', '2881-9-16T0:0:0.0']
timeSelection=['2794-12-1T0:0:0.0', '2844-12-1T0:0:0.0']

spatialSubsetDictionary= {}
spatialSubsetDictionary['latitude']=(-90,90)
#spatialSubsetDictionary['longitude']=(330,359)
spatialSubsetDictionary['longitude']=(-31.265,30.258)

#request subsetted data from feature (can set paths here)
subsetCSML, subsetNetCDF, arraySize=feature.subsetToGridSeries(timeSelection,  csmlpath='myfile.xml', ncpath='myfile.nc',**spatialSubsetDictionary)

#parse and pretty print the result
# strCSML=parser_extra.PrettyPrint(subsetCSML)
# strCSML=parser_extra.removeInlineNS(strCSML)
# print strCSML #csml document (string)
print subsetCSML
print subsetNetCDF # netcdf file (file)
print 'arraySize: %s' %arraySize