In [3]:

    

from suds.client import *
wsdlurl = 'http://daac.ornl.gov/cgi-bin/MODIS/GLBVIZ_1_Glb_subset/MODIS_webservice.wsdl'
client = Client(wsdlurl)
result = client.service.getsubset(40.115,-110.025,"MOD11A2","LST_Day_1km","A2001001","A2001025",1,1)
print result


    

    




(ModisData){
   xllcorner = -9357990.22
   yllcorner = 4458921.58
   cellsize = 926.625433056
   nrows = 3.0
   ncols = 3.0
   band = "LST_Day_1km"
   units = " Kelvin"
   scale = 0.02
   latitude = 40.115
   longitude = -110.025
   header = "HDFname,Product,Date,Location,Processed_Date,Subset_Data... 
"
   subset[] = 
      "MOD11A2.A2001001.h09v04.005.2006350190424.LST_Day_1km,MOD11A2,A2001001,Lat40.115Lon-110.025Samp3Line3,2006350190424,0,14039,14039,0,0,0,0,0,0
",
      "MOD11A2.A2001009.h09v04.005.2006352085954.LST_Day_1km,MOD11A2,A2001009,Lat40.115Lon-110.025Samp3Line3,2006352085954,13571,13604,13605,13583,13604,13590,13609,13589,13589
",
      "MOD11A2.A2001017.h09v04.005.2006354035405.LST_Day_1km,MOD11A2,A2001017,Lat40.115Lon-110.025Samp3Line3,2006354035405,13290,13302,13403,13289,13364,13364,0,0,0
",
      "MOD11A2.A2001025.h09v04.005.2006355205007.LST_Day_1km,MOD11A2,A2001025,Lat40.115Lon-110.025Samp3Line3,2006355205007,13419,13223,13227,13420,13185,13212,13174,13203,13214
",
 }

In [35]:

    

client.service.getproducts()


    

    
Out[35]:





[MCD12Q1,
 MCD12Q2,
 MCD43A1,
 MCD43A2,
 MCD43A4,
 MOD09A1,
 MOD11A2,
 MOD13Q1,
 MOD15A2,
 MOD15A2GFS,
 MOD16A2,
 MOD17A2_51,
 MOD17A3,
 MYD09A1,
 MYD11A2,
 MYD13Q1,
 MYD15A2]

In [36]:

    

url='http://acdisc.gsfc.nasa.gov:80/opendap/Aqua_MODIS_Level3/MYD08_D3.051/2014/MYD08_D3.A2014001.051.2014002194328.pscs_000500796355.hdf'


    

In [ ]:

    

import netCDF4


    

In [38]:

    

nc = netCDF4.Dataset(url)


    

    




---------------------------------------------------------------------------
RuntimeError                              Traceback (most recent call last)
<ipython-input-38-9e734cf63ec2> in <module>()
      1 import netCDF4
----> 2 nc = netCDF4.Dataset(url)

/home/usgs/anaconda/lib/python2.7/site-packages/netCDF4.so in netCDF4.Dataset.__init__ (netCDF4.c:19478)()

RuntimeError: NetCDF: Malformed or inaccessible DAP DATADDS

In [ ]:

    

nc.variables


    

In [4]:

    

result = client.service.getsubset(40.115,-110.025,"MOD11A2","LST_Day_1km","A2001001","A2001009",100,100)


    

In [33]:

    

result.nrows


    

    
Out[33]:





201.0

In [34]:

    

result.band.


    

    
Out[34]:





LST_Day_1km

In [30]:

    

xml = etree.fromstring(result[0])


    

    




---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
<ipython-input-30-3f74620c503d> in <module>()
----> 1 xml = etree.fromstring(result[0])

/home/usgs/anaconda/lib/python2.7/site-packages/lxml/etree.so in lxml.etree.fromstring (src/lxml/lxml.etree.c:68121)()

/home/usgs/anaconda/lib/python2.7/site-packages/lxml/etree.so in lxml.etree._parseMemoryDocument (src/lxml/lxml.etree.c:102459)()

ValueError: can only parse strings

In [32]:

    

len(result)


    

    
Out[32]:





12

In [25]:

    

b=a[0]


    

In [26]:

    

type(b)


    

    
Out[26]:





suds.sax.text.Text

In [17]:

    

len(c)


    

    




---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
<ipython-input-17-c6494b964a51> in <module>()
----> 1 len(c)

TypeError: len() of unsized object

In [ ]: