In [ ]:
from IPython.core.display import HTML
HTML('<iframe src=http://www.nodc.noaa.gov/geoportal/ width=900 height=280></iframe>')
In [1]:
from owslib.csw import CatalogueServiceWeb
In [2]:
# connect to CSW, explore it's properties
endpoint = 'http://www.ngdc.noaa.gov/geoportal/csw' # NGDC Geoportal
#endpoint = 'http://data.nodc.noaa.gov/geoportal/csw' # NODC Geoportal: collection level
#endpoint = 'http://geodiscover.cgdi.ca/wes/serviceManagerCSW/csw' # NRCAN CUSTOM
#endpoint = 'http://geoport.whoi.edu/gi-cat/services/cswiso' # USGS Woods Hole GI_CAT
#endpoint = 'http://cida.usgs.gov/gdp/geonetwork/srv/en/csw' # USGS CIDA Geonetwork
#endpoint = 'http://www.nodc.noaa.gov/geoportal/csw' # NODC Geoportal: granule level
csw = CatalogueServiceWeb(endpoint,timeout=30)
csw.version
Out[2]:
In [3]:
[op.name for op in csw.operations]
Out[3]:
In [4]:
#bbox=[-141,42,-52,84]
bbox=[-71.5, 39.5, -63.0, 46]
csw.getrecords(keywords=['temperature'],bbox=bbox,maxrecords=20)
#csw.getrecords(keywords=['sea_water_temperature'],maxrecords=20)
csw.results
Out[4]:
In [6]:
for rec,item in csw.records.iteritems():
print item.title
In [7]:
print(csw.records.keys())
In [8]:
# choose a sample record
a=csw.records['satellite.G1.ssta.1day']
In [9]:
print a.title
In [10]:
# unfortunately the "uris" property is empty
print a.uris
In [11]:
# yet I can see the URIs here:
print a.xml
In [12]:
# lets look at the references
a.references
Out[12]:
In [13]:
# get specific ServiceType URL from records
def service_urls(records,service_string='urn:x-esri:specification:ServiceType:OPeNDAP'):
urls=[]
for key,rec in records.iteritems():
#create a generator object, and iterate through it until the match is found
#if not found, gets the default value (here "none")
url = next((d['url'] for d in rec.references if d['scheme'] == service_string), None)
if url is not None:
urls.append(url)
return urls
In [14]:
dap_urls = service_urls(csw.records,service_string='urn:x-esri:specification:ServiceType:OPeNDAP')
print dap_urls
In [15]:
# find all the WMS ServiceType URLs
wms_urls = service_urls(csw.records,service_string='urn:x-esri:specification:ServiceType:WMS')
print wms_urls
In [18]:
a.uris
Out[18]:
In [19]:
type a.uris
In [20]:
type(a.uris)
Out[20]:
In [ ]: