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 ", }
client.service.getproducts()
[MCD12Q1, MCD12Q2, MCD43A1, MCD43A2, MCD43A4, MOD09A1, MOD11A2, MOD13Q1, MOD15A2, MOD15A2GFS, MOD16A2, MOD17A2_51, MOD17A3, MYD09A1, MYD11A2, MYD13Q1, MYD15A2]
url='http://acdisc.gsfc.nasa.gov:80/opendap/Aqua_MODIS_Level3/MYD08_D3.051/2014/MYD08_D3.A2014001.051.2014002194328.pscs_000500796355.hdf'
import netCDF4
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
nc.variables
result = client.service.getsubset(40.115,-110.025,"MOD11A2","LST_Day_1km","A2001001","A2001009",100,100)
result.nrows
201.0
result.band.
LST_Day_1km
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
len(result)
12
b=a[0]
type(b)
suds.sax.text.Text
len(c)
--------------------------------------------------------------------------- TypeError Traceback (most recent call last) <ipython-input-17-c6494b964a51> in <module>() ----> 1 len(c) TypeError: len() of unsized object