source: TI02-CSML/trunk/csml/testfiles/gridseries/subset.py @ 3296

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

import csml,sys

f='/home/dom/svn/trunk/csml/testfiles/gridseries/testout.xml'
outputdir='GSout'

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

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

#csmltree=csmldoc.toXML()  
#print csml.parser_extra.PrettyPrint(csmltree)
#sys.exit()

#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])

#get  'gml:id' xml attribute
print feature.id

#get the contents of the feature <description> element
print feature.description.CONTENT

#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)
#subsetDictionary={}
##subsetDictionary['time']=('2794-12-1T0:0:0.0', '2804-12-1T0:0:0.0')
#subsetDictionary['time']='2704-12-01T00:00:00.0'
##subsetDictionary['latitude']=(-90,60)
##subsetDictionary['longitude']=(-31.1265,30.258)
#subsetDictionary['latitude']=(-5,5)
#subsetDictionary['longitude']=(-31.1265,30.258)
#request subsetted data from features (can set paths here) and add to container
#subset a feature

#f.subsetToGridSeries(ncname='tmp1.nc', outputdir='.', latitude=(-90, 90), longitude=(-180,180), time='1975-07-15T00:00:00.0')

    #csmlContainer.add(feature.subsetToGridSeries(outputdir,ncname='mygridseries.nc',**subsetDictionary))
csmlContainer.add(feature.subsetToGridSeries(outputdir,ncname='mygridseries.nc',latitude=(-90,90), longitude=(-180,180),time=('2804-12-01T00:00:00.0','2904-12-01T00:00:00.0')))

#csmlContainer.add(feature.subsetToGridSeries(outputdir,ncname='mygridseries.nc',latitude=(60), longitude=(180),time=('2804-12-01T00:00:00.0','2904-12-01T00:00:00.0')))

##subset another feature (in fact it's the same feature again, but this time subsetting to a profileseries)
#subsetDictionary={}
#subsetDictionary['time']=('2804-12-01T00:00:00.0','2904-12-01T00:00:00.0')
#subsetDictionary['latitude']=(75)
#subsetDictionary['longitude']=(45)

###assigning profileseries feature to variable
##print 'subset grid to profileseries'
#psfeature, file, descriptor =feature.subsetToProfileSeries(outputdir, ncname='myprofileseries.nc',**subsetDictionary)
#csmlContainer.add([psfeature,file, descriptor])

#csmlContainer.add(feature.subsetToProfileSeries(outputdir,ncname='myprofileseries.nc',latitude=(75), longitude=(45),time=('2804-12-01T00:00:00.0','2904-12-01T00:00:00.0')))

##Now subset this New ProfileSeries feature into a Profile feature
#print 'subset  profileseries to profile'
#subsetDictionary={}
#subsetDictionary['time']=('2804-12-01T00:00:00.0')
#csmlContainer.add(psfeature.subsetToProfile(outputdir, ncname='myprofile.nc',**subsetDictionary))


#Now try the subset to profile method of the grid series feature -BROKEN
print 'subsetting gridseries to profile'
subsetDictionary={}
subsetDictionary['time']=('2804-12-01T00:00:00.0')
subsetDictionary['latitude']=(80)
subsetDictionary['longitude']=(55)
csmlContainer.add(feature.subsetToProfile(outputdir, ncname='myprofile.nc',**subsetDictionary))

#Now try the subset to PointSeries method of the grid series feature
print 'subsetting gridseries to pointseries'
subsetDictionary={}
subsetDictionary['time']=('2804-12-01T00:00:00.0','2904-12-01T00:00:00.0')
subsetDictionary['latitude']=80
subsetDictionary['longitude']=55
#subsetDictionary['pressure']=(0)  #BROKEN
csmlContainer.add(feature.subsetToPointSeries(outputdir,ncname='mypointseries.nc',**subsetDictionary))



#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