In [12]:
from owslib.csw import CatalogueServiceWeb
from owslib import fes
In [13]:
endpoint = 'http://comt.sura.org:8000/pycsw' # NODC/UAF Geoportal: granule level
csw = CatalogueServiceWeb(endpoint,timeout=60)
print csw.version
In [14]:
csw.get_operation_by_name('GetRecords').constraints
Out[14]:
In [15]:
val = 'ROMS'
filter1 = fes.PropertyIsLike(propertyname='apiso:AnyText',literal=('*%s*' % val),
escapeChar='\\',wildCard='*',singleChar='?')
filter_list = [ filter1 ]
In [16]:
csw.getrecords2(constraints=filter_list,maxrecords=100,esn='full')
print len(csw.records.keys())
for rec in list(csw.records.keys()):
print csw.records[rec].title
In [17]:
choice=np.random.choice(list(csw.records.keys()))
print(csw.records[choice].title)
csw.records[choice].references
Out[17]:
Add bounding box constraint. To specify lon,lat order for bbox (which we want to do so that we can use the same bbox with either geoportal server or pycsw requests), we need to request the bounding box specifying the CRS84 coordinate reference system. The CRS84 option is available in pycsw 1.1.10+. The ability to specify the crs in the bounding box request is available in owslib 0.8.12+. For more info on the bounding box problem and how it was solved, see this pycsw issue, this geoportal server issue, and this owslib issue
In [18]:
bbox = [-87.40, 24.25, -74.70, 36.70]
bbox_filter = fes.BBox(bbox,crs='urn:ogc:def:crs:OGC:1.3:CRS84')
filter_list = [filter1, bbox_filter]
csw.getrecords2(constraints=filter_list, maxrecords=1000)
In [19]:
print len(csw.records.keys())
for rec in list(csw.records.keys()):
print csw.records[rec].title
In [19]:
In [19]:
In [ ]: