OWSLib 0.7-dev documentation¶

Tom Kralidis · 2013-04-12

$ easy_install OWSLib

$ git clone git://github.com/geopython/OWSLib.git

$ zypper ar http://download.opensuse.org/repositories/Application:/Geo/openSUSE_12.1/ GEO
$ zypper refresh
$ zypper install owslib

$ wget -O /etc/yum.repos.d/CentOS-CBS.repo http://download.opensuse.org/repositories/Application:/Geo/CentOS_6/Application:Geo.repo
$ yum install owslib

$ wget -O /etc/yum.repos.d/RHEL-CBS.repo http://download.opensuse.org/repositories/Application:/Geo/RHEL_6/Application:Geo.repo
$ yum install owslib

>>> from owslib.wms import WebMapService
>>> wms = WebMapService('http://wms.jpl.nasa.gov/wms.cgi', version='1.1.1')
>>> wms.identification.type
'OGC:WMS'
>>> wms.identification.title
'JPL Global Imagery Service'

>>> list(wms.contents)
['global_mosaic', 'global_mosaic_base', 'us_landsat_wgs84', 'srtm_mag', 'daily_terra_721', 'daily_aqua_721', 'daily_terra_ndvi', 'daily_aqua_ndvi', 'daily_terra', 'daily_aqua', 'BMNG', 'modis', 'huemapped_srtm', 'srtmplus', 'worldwind_dem', 'us_ned', 'us_elevation', 'us_colordem']

>>> wms['global_mosaic'].title
'WMS Global Mosaic, pan sharpened'
>>> wms['global_mosaic'].queryable
0
>>> wms['global_mosaic'].opaque
0
>>> wms['global_mosaic'].boundingBox
>>> wms['global_mosaic'].boundingBoxWGS84
(-180.0, -60.0, 180.0, 84.0)
>>> wms['global_mosaic'].crsOptions
['EPSG:4326', 'AUTO:42003']
>>> wms['global_mosaic'].styles
{'pseudo_bright': {'title': 'Pseudo-color image (Uses IR and Visual bands, 542 mapping), gamma 1.5'}, 'pseudo': {'title': '(default) Pseudo-color image, pan sharpened (Uses IR and Visual bands, 542 mapping), gamma 1.5'}, 'visual': {'title': 'Real-color image, pan sharpened (Uses the visual bands, 321 mapping), gamma 1.5'}, 'pseudo_low': {'title': 'Pseudo-color image, pan sharpened (Uses IR and Visual bands, 542 mapping)'}, 'visual_low': {'title': 'Real-color image, pan sharpened (Uses the visual bands, 321 mapping)'}, 'visual_bright': {'title': 'Real-color image (Uses the visual bands, 321 mapping), gamma 1.5'}}

>>> [op.name for op in wms.operations]
['GetCapabilities', 'GetMap']
>>> wms.getOperationByName('GetMap').methods
{'Get': {'url': 'http://wms.jpl.nasa.gov/wms.cgi?'}}
>>> wms.getOperationByName('GetMap').formatOptions
['image/jpeg', 'image/png', 'image/geotiff', 'image/tiff']

>>> img = wms.getmap(   layers=['global_mosaic'],
...                     styles=['visual_bright'],
...                     srs='EPSG:4326',
...                     bbox=(-112, 36, -106, 41),
...                     size=(300, 250),
...                     format='image/jpeg',
...                     transparent=True
...                     )
>>> out = open('jpl_mosaic_visb.jpg', 'wb')
>>> out.write(img.read())
>>> out.close()

>>> from owslib.csw import CatalogueServiceWeb
>>> csw = CatalogueServiceWeb('http://geodiscover.cgdi.ca/wes/serviceManagerCSW/csw')
>>> csw.identification.type
'CSW'
>>> [op.name for op in csw.operations]
['GetCapabilities', 'GetRecords', 'GetRecordById', 'DescribeRecord', 'GetDomain']

>>> csw.getdomain('GetRecords.resultType')
>>> csw.results
{'values': ['results', 'validate', 'hits'], 'parameter': 'GetRecords.resultType', 'type': 'csw:DomainValuesType'}
>>>

>>> csw.getrecords(keywords=['birds'], maxrecords=20)
>>> csw.results
{'matches': 101, 'nextrecord': 21, 'returned': 20}
>>> for rec in csw.records:
...     print csw.records[rec].title
...
ALLSPECIES
NatureServe Canada References
Bird Studies Canada - BirdMap WMS
Parks Canada Geomatics Metadata Repository
Bird Studies Canada - BirdMap WFS
eBird Canada - Survey Locations
WHC CitizenScience WMS
Project FeederWatch - Survey Locations
North American Bird Banding and Encounter Database
Wildlife Habitat Canada CitizenScience WFS
Parks Canada Geomatics Metadata Repository
Parks Canada Geomatics Metadata Repository
Wildlife Habitat Canada CitizenScience WMS
Canadian IBA Polygon layer
Land
Wildlife Habitat Canada CitizenScience WMS
WATER
Parks Canada Geomatics Metadata Repository
Breeding Bird Survey
SCALE
>>>

>>> csw.getrecords(keywords=['birds'],bbox=[-141,42,-52,84])
>>> csw.results
{'matches': 3, 'nextrecord': 0, 'returned': 3}
>>>

>>> csw.getrecords(keywords=['birds', 'fowl'])
>>> csw.results
{'matches': 107, 'nextrecord': 11, 'returned': 10}
>>>

>>> csw.getrecordbyid(id=['9250AA67-F3AC-6C12-0CB9-0662231AA181'])
>>> c.records['9250AA67-F3AC-6C12-0CB9-0662231AA181'].title
'ALLSPECIES'

>>> csw.getrecords(cql='csw:AnyText like "%birds%"')

>>> csw.transaction(ttype='insert', typename='gmd:MD_Metadata', record=open(file.xml).read())

>>> # update ALL records
>>> csw.transaction(ttype='update', typename='csw:Record', propertyname='dc:title', propertyvalue='New Title')
>>> # update records satisfying keywords filter
>>> csw.transaction(ttype='update', typename='csw:Record', propertyname='dc:title', propertyvalue='New Title', keywords=['birds','fowl'])
>>> # update records satisfying BBOX filter
>>> csw.transaction(ttype='update', typename='csw:Record', propertyname='dc:title', propertyvalue='New Title', bbox=[-141,42,-52,84])

>>> # delete ALL records
>>> csw.transaction(ttype='delete', typename='gmd:MD_Metadata')
>>> # delete records satisfying keywords filter
>>> csw.transaction(ttype='delete', typename='gmd:MD_Metadata', keywords=['birds','fowl'])
>>> # delete records satisfying BBOX filter
>>> csw.transaction(ttype='delete', typename='gmd:MD_Metadata', bbox=[-141,42,-52,84])

>>> csw.harvest('http://host/url.xml', 'http://www.isotc211.org/2005/gmd')

>>> from owslib.wps import WebProcessingService
>>> wps = WebProcessingService('http://cida.usgs.gov/climate/gdp/process/WebProcessingService', verbose=False, skip_caps=True)
>>> wps.getcapabilities()
>>> wps.identification.type
'WPS'
>>> wps.identification.title
'Geo Data Portal WPS Processing'
>>> wps.identification.abstract
'Geo Data Portal WPS Processing'
>>> for operation in wps.operations:
...     operation.name
...
'GetCapabilities'
'DescribeProcess'
'Execute'
>>> for process in wps.processes:
...     process.identifier, process.title
...
('gov.usgs.cida.gdp.wps.algorithm.FeatureCoverageIntersectionAlgorithm', 'Feature Coverage WCS Intersection')
('gov.usgs.cida.gdp.wps.algorithm.FeatureCoverageOPeNDAPIntersectionAlgorithm', 'Feature Coverage OPeNDAP Intersection')
('gov.usgs.cida.gdp.wps.algorithm.FeatureCategoricalGridCoverageAlgorithm', 'Feature Categorical Grid Coverage')
('gov.usgs.cida.gdp.wps.algorithm.FeatureWeightedGridStatisticsAlgorithm', 'Feature Weighted Grid Statistics')
('gov.usgs.cida.gdp.wps.algorithm.FeatureGridStatisticsAlgorithm', 'Feature Grid Statistics')
('gov.usgs.cida.gdp.wps.algorithm.PRMSParameterGeneratorAlgorithm', 'PRMS Parameter Generator')
>>>

>>> from owslib.wps import printInputOutput
>>> process = wps.describeprocess('gov.usgs.cida.gdp.wps.algorithm.FeatureWeightedGridStatisticsAlgorithm')
>>> process.identifier
'gov.usgs.cida.gdp.wps.algorithm.FeatureWeightedGridStatisticsAlgorithm'
>>> process.title
'Feature Weighted Grid Statistics'
>>> process.abstract
'This algorithm generates area weighted statistics of a gridded dataset for a set of vector polygon features. Using the bounding-box that encloses ...
>>> for input in process.dataInputs:
...     printInputOutput(input)
...
 identifier=FEATURE_COLLECTION, title=Feature Collection, abstract=A feature collection encoded as a WFS request or one of the supported GML profiles.,...
 Supported Value: mimeType=text/xml, encoding=UTF-8, schema=http://schemas.opengis.net/gml/2.0.0/feature.xsd
 Supported Value: mimeType=text/xml, encoding=UTF-8, schema=http://schemas.opengis.net/gml/2.1.1/feature.xsd
 Supported Value: mimeType=text/xml, encoding=UTF-8, schema=http://schemas.opengis.net/gml/2.1.2/feature.xsd
 Supported Value: mimeType=text/xml, encoding=UTF-8, schema=http://schemas.opengis.net/gml/2.1.2.1/feature.xsd
 Supported Value: mimeType=text/xml, encoding=UTF-8, schema=http://schemas.opengis.net/gml/3.0.0/base/feature.xsd
 Supported Value: mimeType=text/xml, encoding=UTF-8, schema=http://schemas.opengis.net/gml/3.0.1/base/feature.xsd
 Supported Value: mimeType=text/xml, encoding=UTF-8, schema=http://schemas.opengis.net/gml/3.1.0/base/feature.xsd
 Supported Value: mimeType=text/xml, encoding=UTF-8, schema=http://schemas.opengis.net/gml/3.1.1/base/feature.xsd
 Supported Value: mimeType=text/xml, encoding=UTF-8, schema=http://schemas.opengis.net/gml/3.2.1/base/feature.xsd
 Default Value: mimeType=text/xml, encoding=UTF-8, schema=http://schemas.opengis.net/gml/2.0.0/feature.xsd
 minOccurs=1, maxOccurs=1
 identifier=DATASET_URI, title=Dataset URI, abstract=The base data web service URI for the dataset of interest., data type=anyURI
 Allowed Value: AnyValue
 Default Value: None
 minOccurs=1, maxOccurs=1
 identifier=DATASET_ID, title=Dataset Identifier, abstract=The unique identifier for the data type or variable of interest., data type=string
 Allowed Value: AnyValue
 Default Value: None
 minOccurs=1, maxOccurs=2147483647
 identifier=REQUIRE_FULL_COVERAGE, title=Require Full Coverage, abstract=If turned on, the service will require that the dataset of interest ....
 Allowed Value: True
 Default Value: True
 minOccurs=1, maxOccurs=1
 identifier=TIME_START, title=Time Start, abstract=The date to begin analysis., data type=dateTime
 Allowed Value: AnyValue
 Default Value: None
 minOccurs=0, maxOccurs=1
 identifier=TIME_END, title=Time End, abstract=The date to end analysis., data type=dateTime
 Allowed Value: AnyValue
 Default Value: None
 minOccurs=0, maxOccurs=1
 identifier=FEATURE_ATTRIBUTE_NAME, title=Feature Attribute Name, abstract=The attribute that will be used to label column headers in processing output., ...
 Allowed Value: AnyValue
 Default Value: None
 minOccurs=1, maxOccurs=1
 identifier=DELIMITER, title=Delimiter, abstract=The delimiter that will be used to separate columns in the processing output., data type=string
 Allowed Value: COMMA
 Allowed Value: TAB
 Allowed Value: SPACE
 Default Value: COMMA
 minOccurs=1, maxOccurs=1
 identifier=STATISTICS, title=Statistics, abstract=Statistics that will be returned for each feature in the processing output., data type=string
 Allowed Value: MEAN
 Allowed Value: MINIMUM
 Allowed Value: MAXIMUM
 Allowed Value: VARIANCE
 Allowed Value: STD_DEV
 Allowed Value: SUM
 Allowed Value: COUNT
 Default Value: None
 minOccurs=1, maxOccurs=7
 identifier=GROUP_BY, title=Group By, abstract=If multiple features and statistics are selected, this will change whether the processing output ...
 Allowed Value: STATISTIC
 Allowed Value: FEATURE_ATTRIBUTE
 Default Value: None
 minOccurs=1, maxOccurs=1
 identifier=SUMMARIZE_TIMESTEP, title=Summarize Timestep, abstract=If selected, processing output will include columns with summarized statistics ...
 Allowed Value: True
 Default Value: True
 minOccurs=0, maxOccurs=1
 identifier=SUMMARIZE_FEATURE_ATTRIBUTE, title=Summarize Feature Attribute, abstract=If selected, processing output will include a final row of ...
 Allowed Value: True
 Default Value: True
 minOccurs=0, maxOccurs=1
>>> for output in process.processOutputs:
...     printInputOutput(output)
...
 identifier=OUTPUT, title=Output File, abstract=A delimited text file containing requested process output., data type=ComplexData
 Supported Value: mimeType=text/csv, encoding=UTF-8, schema=None
 Default Value: mimeType=text/csv, encoding=UTF-8, schema=None
 reference=None, mimeType=None
>>>

>>> from owslib.wps import GMLMultiPolygonFeatureCollection
>>> polygon = [(-102.8184, 39.5273), (-102.8184, 37.418), (-101.2363, 37.418), (-101.2363, 39.5273), (-102.8184, 39.5273)]
>>> featureCollection = GMLMultiPolygonFeatureCollection( [polygon] )
>>> processid = 'gov.usgs.cida.gdp.wps.algorithm.FeatureWeightedGridStatisticsAlgorithm'
>>> inputs = [ ("FEATURE_ATTRIBUTE_NAME","the_geom"),
...            ("DATASET_URI", "dods://cida.usgs.gov/qa/thredds/dodsC/derivatives/derivative-days_above_threshold.pr.ncml"),
...            ("DATASET_ID", "ensemble_b1_pr-days_above_threshold"),
...            ("TIME_START","2010-01-01T00:00:00.000Z"),
...            ("TIME_END","2011-01-01T00:00:00.000Z"),
...            ("REQUIRE_FULL_COVERAGE","false"),
...            ("DELIMITER","COMMA"),
...            ("STATISTICS","MEAN"),
...            ("GROUP_BY","STATISTIC"),
...            ("SUMMARIZE_TIMESTEP","false"),
...            ("SUMMARIZE_FEATURE_ATTRIBUTE","false"),
...            ("FEATURE_COLLECTION", featureCollection)
...           ]
>>> output = "OUTPUT"
>>> execution = wps.execute(processid, inputs, output = "OUTPUT")
Executing WPS request...
Execution status=ProcessStarted
>>> from owslib.wps import monitorExecution
>>> monitorExecution(execution)

Checking execution status... (location=http://cida.usgs.gov/climate/gdp/process/RetrieveResultServlet?id=6809217153012787208)
Execution status=ProcessSucceeded
Execution status: ProcessSucceeded
Output URL=http://cida.usgs.gov/climate/gdp/process/RetrieveResultServlet?id=6809217153012787208OUTPUT.3cbcd666-a912-456f-84a3-6ede450aca95
>>>

>>> from owslib.wps import WFSQuery, WFSFeatureCollection
>>> wfsUrl = "http://cida.usgs.gov/climate/gdp/proxy/http://igsarm-cida-gdp2.er.usgs.gov:8082/geoserver/wfs"
>>> query = WFSQuery("sample:CONUS_States", propertyNames=['the_geom',"STATE"], filters=["CONUS_States.508","CONUS_States.469"])
>>> featureCollection = WFSFeatureCollection(wfsUrl, query)
>>> # same process submission as above
...

>>> request = open('/Users/cinquini/Documents/workspace-cog/wps/tests/wps_USGSExecuteRequest1.xml','r').read()
>>> execution = wps.execute(None, [], request=request)
Executing WPS request...
Execution status=ProcessStarted
>>> monitorExecution(execution)

Checking execution status... (location=http://cida.usgs.gov/climate/gdp/process/RetrieveResultServlet?id=5103866488472745994)
Execution status=ProcessSucceeded
Execution status: ProcessSucceeded
Output URL=http://cida.usgs.gov/climate/gdp/process/RetrieveResultServlet?id=5103866488472745994OUTPUT.f80e2a78-96a9-4343-9777-be60fac5b256

>>> from owslib.sos import SensorObservationService
>>> from owslib.fes import FilterCapabilities
>>> from owslib.ows import OperationsMetadata
>>> from owslib.crs import Crs
>>> from datetime import datetime

>>> xml = open('sos_ndbc_getcapabilities.xml').read()
>>> ndbc = SensorObservationService(None, xml=xml)

>>> id = ndbc.identification

>>> id.service
'OGC:SOS'

>>> id.version
'1.0.0'

>>> id.title
'National Data Buoy Center SOS'

>>> id.abstract
'National Data Buoy Center SOS'

>>> id.keywords
['Weather', 'Ocean Currents', 'Air Temperature', 'Water Temperature', 'Conductivity', 'Salinity', 'Barometric Pressure', 'Water Level', 'Waves', 'Winds', 'NDBC']

>>> id.fees
'NONE'

>>> id.accessconstraints
'NONE'

>>> p = ndbc.provider

>>> p.name
'National Data Buoy Center'

>>> p.url
'http://sdf.ndbc.noaa.gov/'

>>> sc = p.contact

>>> sc.role
>>> sc.position
>>> sc.instructions
>>> sc.organization
>>> sc.fax
>>> sc.hours

>>> sc.name
'Webmaster'

>>> sc.phone
'228-688-2805'

>>> sc.address
'Bldg. 3205'

>>> sc.city
'Stennis Space Center'

>>> sc.region
'MS'

>>> sc.postcode
'39529'

>>> sc.country
'USA'

>>> sc.email
'webmaster.ndbc@noaa.gov'

>>> o = ndbc.operations
>>> len(o)
3

>>> len(o)
3

>>> getcap = ndbc.get_operation_by_name('GetCapabilities')
>>> isinstance(getcap, OperationsMetadata)
True

>>> getcap = ndbc.get_operation_by_name('getcapabilities')
>>> isinstance(getcap, OperationsMetadata)
True

>>> getcap.constraints
{}

>>> getcap.parameters
{'Sections': {'values': ['ServiceIdentification', 'ServiceProvider', 'OperationsMetadata', 'Contents', 'All']}}

>>> getcap.methods
{'Post': {'url': 'http://sdf.ndbc.noaa.gov/sos/server.php'}, 'Get': {'url': 'http://sdf.ndbc.noaa.gov/sos/server.php'}}

>>> descsen = ndbc.get_operation_by_name('describesensor')

>>> descsen.constraints
{}

>>> descsen.parameters
{'outputFormat': {'values': ['text/xml;subtype="sensorML/1.0.1"']}}

>>> descsen.methods
{'Post': {'url': 'http://sdf.ndbc.noaa.gov/sos/server.php'}, 'Get': {'url': 'http://sdf.ndbc.noaa.gov/sos/server.php'}}

>>> getob = ndbc.get_operation_by_name('getobservation')

>>> getob.constraints
{}

>>> getob.parameters
{'observedProperty': {'values': ['air_temperature', 'air_pressure_at_sea_level', 'sea_water_electrical_conductivity', 'currents', 'sea_water_salinity', 'sea_floor_depth_below_sea_surface', 'sea_water_temperature', 'waves', 'winds']}}

>>> getob.methods
{'Post': {'url': 'http://sdf.ndbc.noaa.gov/sos/server.php'}, 'Get': {'url': 'http://sdf.ndbc.noaa.gov/sos/server.php'}}

>>> filter = ndbc.filters
>>> isinstance(filter, FilterCapabilities)
False

>>> contents = ndbc.contents
>>> len(contents)
848

>>> network = contents['network-all']
>>> network.id
'network-all'

>>> network.name
'urn:ioos:network:noaa.nws.ndbc:all'

>>> network.description
'All stations on the NDBC SOS server'

>>> srs = network.srs
>>> isinstance(srs, Crs)
True
>>> srs.getcodeurn()
'urn:ogc:def:crs:EPSG::4326'
>>> srs.getcode()
'EPSG:4326'

>>> network.bbox
(-179.995, -77.465999999999994, 180.0, 80.810000000000002)
>>> bbsrs = network.bbox_srs
>>> isinstance(bbsrs, Crs)
True
>>> bbsrs.getcodeurn()
'urn:ogc:def:crs:EPSG::4326'
>>> bbsrs.getcode()
'EPSG:4326'

>>> bp = network.begin_position
>>> isinstance(bp, datetime)
True
>>> ep = network.end_position
>>> isinstance(ep, datetime)
True

>>> network.result_model
'om:Observation'

>>> procs = network.procedures
>>> len(procs)
847

>>> network.observed_properties
['http://mmisw.org/ont/cf/parameter/air_temperature', 'http://mmisw.org/ont/cf/parameter/air_pressure_at_sea_level', 'http://mmisw.org/ont/cf/parameter/sea_water_electrical_conductivity', 'http://mmisw.org/ont/cf/parameter/currents', 'http://mmisw.org/ont/cf/parameter/sea_water_salinity', 'http://mmisw.org/ont/cf/parameter/sea_floor_depth_below_sea_surface', 'http://mmisw.org/ont/cf/parameter/sea_water_temperature', 'http://mmisw.org/ont/cf/parameter/waves', 'http://mmisw.org/ont/cf/parameter/winds']

>>> foi = network.features_of_interest
>>> len(foi)
1082

>>> rfs = network.response_formats
>>> len(rfs)
5
>>> rfs[-1]
'application/vnd.google-earth.kml+xml'

>>> rms = network.response_modes
>>> len(rms)
1
>>> rms[0]
'inline'

>>> station = contents['station-zbqn7']
>>> station.id
'station-zbqn7'

>>> station.name
'urn:ioos:station:wmo:zbqn7'

>>> station.description
"Zeke's Basin, North Carolina"

>>> srs = station.srs
>>> isinstance(srs, Crs)
True
>>> srs.getcodeurn()
'urn:ogc:def:crs:EPSG::4326'
>>> srs.getcode()
'EPSG:4326'
>>> station.bbox
(-77.935000000000002, 33.954999999999998, -77.935000000000002, 33.954999999999998)

>>> bbsrs = station.bbox_srs
>>> isinstance(bbsrs, Crs)
True
>>> bbsrs.getcodeurn()
'urn:ogc:def:crs:EPSG::4326'
>>> bbsrs.getcode()
'EPSG:4326'

>>> bp = station.begin_position
>>> isinstance(bp, datetime)
True
>>> ep = station.end_position
>>> isinstance(ep, datetime)
True

>>> station.result_model
'om:Observation'

>>> procs = station.procedures
>>> len(procs)
1
>>> procs[0]
'urn:ioos:station:wmo:zbqn7'

>>> ops = station.observed_properties
>>> len(ops)
3
>>> ops[0]
'http://mmisw.org/ont/cf/parameter/sea_water_electrical_conductivity'

>>> foi = station.features_of_interest
>>> len(foi)
1
>>> foi[0]
'urn:cgi:Feature:CGI:EarthOcean'

>>> rfs = station.response_formats
>>> len(rfs)
5
>>> rfs[0]
'text/xml;schema="ioos/0.6.1"'

>>> rm = station.response_modes
>>> len(rm)
1
>>> rm[0]
'inline'

>>> from owslib.sos import SensorObservationService
>>> from owslib.ows import OperationsMetadata
>>> from owslib.fes import FilterCapabilities
>>> from owslib.crs import Crs
>>> from datetime import datetime

>>> xml = open('sos_ndbc_getcapabilities.xml').read()
>>> ndbc = SensorObservationService(None, xml=xml)

>>> off = ndbc.offerings[1]
>>> offerings = [off.name]
>>> responseFormat = off.response_formats[0]
>>> observedProperties = [off.observed_properties[0]]
>>> #observedProperties = [ndbc.get_operation_by_name('GetObservation').parameters['observedProperty']['values'][0]]
>>> eventTime = "This is not a valid eventTime!"
>>> response = ndbc.get_observation(offerings=offerings, responseFormat=responseFormat, observedProperties=observedProperties, eventTime=eventTime)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
ExceptionReport: 'This is not a valid eventTime!'

>>> off = ndbc.offerings[1]
>>> offerings = [off.name]
>>> responseFormat = off.response_formats[0]
>>> observedProperties = [off.observed_properties[0]]
>>> #observedProperties = [ndbc.get_operation_by_name('GetObservation').parameters['observedProperty']['values'][0]]
>>> eventTime = None

>>> response = ndbc.get_observation(offerings=offerings, responseFormat=responseFormat, observedProperties=observedProperties, eventTime=eventTime)

>>> procedure = "foobar"
>>> outputFormat = ndbc.get_operation_by_name('DescribeSensor').parameters['outputFormat']['values'][0]
>>> response = ndbc.describe_sensor(procedure=procedure, outputFormat=outputFormat)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
ExceptionReport: 'foobar'

>>> procedure = ndbc.offerings[1].procedures[0]
>>> outputFormat = ndbc.get_operation_by_name('DescribeSensor').parameters['outputFormat']['values'][0]
>>> response = ndbc.describe_sensor(procedure=procedure, outputFormat=outputFormat)

>>> from owslib.swe.sensor.sml import SensorML, EventMetadata
>>> from dateutil import parser
>>> import pytz

>>> xml = open('sml_ndbc_station.xml').read()
>>> root = SensorML(xml)
>>> system = root.systems[0]

>>> system.id
'station-41012'

>>> system.description
'Station metadata for 41012 - 40NM ENE of St Augustine, FL'

>>> system.contacts.keys()
['urn:ogc:def:classifiers:OGC:contactType:operator', 'urn:ogc:def:classifiers:OGC:contactType:publisher']

>>> operator = system.get_contact_by_role('urn:ogc:def:classifiers:OGC:contactType:operator')
>>> operator.role
'urn:ogc:def:classifiers:OGC:contactType:operator'
>>> operator.organization
'National Data Buoy Center'
>>> operator.country
'US'

>>> publisher = system.get_contact_by_role('urn:ogc:def:classifiers:OGC:contactType:publisher')
>>> publisher.role
'urn:ogc:def:classifiers:OGC:contactType:publisher'
>>> publisher.organization
'National Data Buoy Center'
>>> publisher.country
'USA'
>>> publisher.phone
'228-688-2805'
>>> publisher.address
'Bldg. 3205'
>>> publisher.city
'Stennis Space Center'
>>> publisher.postcode
'39529'
>>> publisher.email
'webmaster.ndbc@noaa.gov'
>>> publisher.region
'MS'
>>> publisher.url
'http://sdf.ndbc.noaa.gov/'

>>> system.identifiers.keys()
['StationId', 'Long Name', 'Short Name']

>>> sid = system.get_identifier_by_name('StationId')
>>> sid.name
'StationId'
>>> sid.definition
'urn:ioos:def:identifier:NOAA:stationID'
>>> sid.codespace
'http://sdf.ndbc.noaa.gov'
>>> sid.value
'urn:ioos:station:wmo:41012'

>>> system.classifiers.keys()
['Platform Type']
>>> classi = system.get_classifier_by_name('Platform type')
>>> classi.name
'Platform Type'
>>> classi.definition
'urn:ioos:def:classifier:NOAA:platformType'
>>> classi.codespace
'http://sdf.ndbc.noaa.gov'
>>> classi.value
'MOORED BUOY'

>>> system.documents 
[<owslib.swe.sensor.sml.DocumentationMetadata ...>]

>>> doc = system.documents[0]
>>> doc.description
'Handbook of Automated Data Quality Control Checks and Procedures, National Data Buoy Center, August 2009'
>>> doc.format
'pdf'
>>> doc.url
'http://www.ndbc.noaa.gov/NDBCHandbookofAutomatedDataQualityControl2009.pdf'

>>> system.history.keys()
['deployment_start', 'deployment_stop']

>>> his = system.get_history_by_name('deployment_start')
>>> his 
[<owslib.swe.sensor.sml.EventMetadata ...>]
>>> len(his)
2

>>> event = his[0]
>>> event.name
'deployment_start'
>>> parser.parse(event.date).replace(tzinfo=pytz.utc).isoformat()
'2010-01-12T00:00:00+00:00'
>>> event.description
'Deployment start event'
>>> event.documentation_url
'http://sdftest.ndbc.noaa.gov/sos/server.php?service=SOS&request=DescribeSensor&version=1.0.0&outputformat=text/xml;subtype="sensorML/1.0.1"&procedure=urn:ioos:station:wmo:41012:20100112'

>>> from owslib.iso import *
>>> m=MD_Metadata(etree.parse('tests/data/9250AA67-F3AC-6C12-0CB9-0662231AA181_iso.xml')
>>> m.identification.topiccategory
'farming'
>>>

>>> import urllib2
>>> from owslib.etree import etree
>>> from owslib.iso import CodelistCatalogue
>>> e=etree.parse('gmxCodelists.xml')
>>> c=CodelistCatalogue(e)
>>> print c.getcodelistdictionaries()
['MD_RestrictionCode', 'CI_OnLineFunctionCode', 'DQ_EvaluationMethodTypeCode', 'DS_InitiativeTypeCode', 'MD_MaintenanceFrequencyCode', 'MD_SpatialRepresentationTypeCode', 'MD_ImagingConditionCode', 'MD_ObligationCode', 'MD_GeometricObjectTypeCode', 'MD_ScopeCode', 'MD_KeywordTypeCode', 'MD_CharacterSetCode', 'CI_PresentationFormCode', 'MX_ScopeCode', 'CI_DateTypeCode', 'MD_TopicCategoryCode', 'MD_DatatypeCode', 'MD_ProgressCode', 'MD_MediumNameCode', 'MD_ClassificationCode', 'DS_AssociationTypeCode', 'CI_RoleCode', 'MD_TopologyLevelCode', 'MD_DimensionNameTypeCode', 'MD_CoverageContentTypeCode', 'MD_MediumFormatCode', 'MD_PixelOrientationCode', 'MD_CellGeometryCode']
>>> print c.getcodedefinitionidentifiers('CI_RoleCode')
['author', 'processor', 'publisher', 'custodian', 'pointOfContact', 'distributor', 'user', 'resourceProvider', 'originator', 'owner', 'principalInvestigator']

>>> from owslib import crs
>>> c=crs.Crs('EPSG:4326')
>>> c.code
4326
>>> c=crs.Crs('urn:ogc:def:crs:EPSG::4326')
>>> c.authority
'EPSG'
>>> c.axisorder
'yx'
>>> c=crs.Crs('http://www.opengis.net/gml/epsg.xml#4326')
>>> c.code
4326
>>> c.axisorder
'xy'
>>> c=crs.Crs('urn:x-ogc:def:crs:EPSG:6.11:2192')
>>> c.axisorder
'xy'
>>> c.code
2192
>>> c.version
'6.11'
>>> c=crs.Crs('http://www.opengis.net/def/crs/EPSG/0/4326')
>>> c.authority
'EPSG'
>>> c.code
4326
>>> c.axisorder
'yx'

>>> from owslib.wmts import WebMapTileService
>>> wmts = WebMapTileService("http://map1c.vis.earthdata.nasa.gov/wmts-geo/wmts.cgi")
>>> wmts.identification.type
'OGC WMTS'
>>> wmts.identification.version
'1.0.0'
>>> wmts.identification.title
'NASA Global Image Browse Services for EOSDIS'
>>> str.strip(wmts.identification.abstract)
'Near real time imagery from multiple NASA instruments'
>>> wmts.identification.keywords
['World', 'Global']

>>> wmts.provider.name
'National Aeronautics and Space Administration'

>>> wmts.provider.url
'http://earthdata.nasa.gov/'

>>> list(wmts.contents)
['MODIS_Aqua_SurfaceReflectance_Bands121', 'MODIS_Aqua_Cloud_Top_Temp_Night', 'MODIS_Aqua_Sea_Ice', 'MODIS_Aqua_SurfaceReflectance_Bands721', 'MODIS_Terra_SurfaceReflectance_Bands143', 'MODIS_Aqua_Snow_Cover', 'MODIS_Terra_CorrectedReflectance_Bands721', 'MODIS_Terra_Cloud_Top_Pressure_Day', 'AIRS_Dust_Score', 'MODIS_Aqua_Cloud_Top_Pressure_Night', 'MODIS_Aqua_Brightness_Temp_Band31_Day', 'MODIS_Land_Water_Mask', 'OMI_Aerosol_Index', 'MODIS_Aqua_Land_Surface_Temp_Day', 'MODIS_Aqua_Aerosol', 'OMI_Absorbing_Aerosol_Optical_Depth', 'MODIS_Terra_Cloud_Top_Pressure_Night', 'AIRS_Prata_SO2_Index_Day', 'MODIS_Terra_Sea_Ice', 'OMI_Cloud_Pressure', 'MODIS_Aqua_CorrectedReflectance_TrueColor', 'MODIS_Terra_Water_Vapor_5km_Day', 'OMI_Aerosol_Optical_Depth', 'MODIS_Aqua_Water_Vapor_5km_Day', 'MODIS_Terra_SurfaceReflectance_Bands121', 'MODIS_Terra_CorrectedReflectance_Bands367', 'OMI_SO2_Planetary_Boundary_Layer', 'AIRS_CO_Total_Column_Night', 'MODIS_Aqua_Cloud_Top_Pressure_Day', 'AIRS_Precipitation_Day', 'MODIS_Terra_Data_No_Data', 'MODIS_Terra_Brightness_Temp_Band31_Day', 'AIRS_Precipitation_Night', 'MODIS_Terra_Land_Surface_Temp_Day', 'MODIS_Terra_Cloud_Top_Temp_Night', 'MODIS_Terra_Aerosol', 'MODIS_Terra_Brightness_Temp_Band31_Night', 'MODIS_Aqua_Brightness_Temp_Band31_Night', 'AIRS_CO_Total_Column_Day', 'MODIS_Terra_SurfaceReflectance_Bands721', 'MODIS_Terra_Cloud_Top_Temp_Day', 'MODIS_Terra_CorrectedReflectance_TrueColor', 'MODIS_Terra_Land_Surface_Temp_Night', 'MODIS_Aqua_Cloud_Top_Temp_Day', 'OMI_SO2_Middle_Troposphere', 'MODIS_Terra_Water_Vapor_5km_Night', 'OMI_SO2_Lower_Troposphere', 'AIRS_Prata_SO2_Index_Night', 'MODIS_Aqua_Land_Surface_Temp_Night', 'MODIS_Aqua_SurfaceReflectance_Bands143', 'MODIS_Terra_Snow_Cover', 'MODIS_Aqua_CorrectedReflectance_Bands721', 'MODIS_Aqua_Data_No_Data', 'OMI_SO2_Upper_Troposphere_and_Stratosphere', 'MODIS_Aqua_Water_Vapor_5km_Night']

>>> tile = wmts.gettile(layer='MODIS_Terra_CorrectedReflectance_TrueColor', tilematrixset='EPSG4326_250m', tilematrix='0', row=0, column=0, format="image/jpeg")
>>> out = open('nasa_modis_terra_truecolour.jpg', 'wb')
>>> out.write(tile.read())
>>> out.close()

>>> from owslib.waterml.wml11 import WaterML_1_1 as wml

>>> f = open('tests/cuahsi_example_all_sites.xml', 'r')
>>> sites = wml(f).response

>>> sites.query_info.creation_time
datetime.datetime(2009, 6, 12, 10, 47, 54, 531250, tzinfo=tzoffset(None, -25200))
>>> sites.query_info.notes
['ALL Sites(empty request)']
>>> sites.query_info.criteria.method_called
'GetSites'

>>> codes = sites.site_codes
>>> codes
[['USU-LBR-Mendon'], ['USU-LBR-Paradise'], ['USU-LBR-ExpFarm'], ['USU-LBR-SFLower'], ['USU-LBR-EFLower'], ['USU-LBR-EFWeather'], ['USU-LBR-SFUpper'], ['USU-LBR-ParadiseRepeater'], ['USU-LBR-EFRepeater'], ['USU-LBR-Wellsville'], ['USU-LBR-Confluence'], ['10105900']]

>>> sites.site_names
['Little Bear River at Mendon Road near Mendon, Utah', 'Little Bear River at McMurdy Hollow near Paradise, Utah', 'Utah State University Experimental Farm near Wellsville, Utah', 'South Fork Little Bear River below Davenport Creek near Avon, Utah', 'East Fork Little Bear River at Paradise Canal Diversion near Avon, Utah', 'Little Bear River Upper Weather Station near Avon, Utah', 'South Fork Little Bear River above Davenport Creek near Avon, Utah', 'Radio Repeater near Paradise, Utah', 'East Fork Little Bear River Radio Repeater near Avon, Utah', 'Little Bear River near Wellsville, Utah', 'Little Bear River below Confluence of South and East Forks near Avon, Utah', 'Little Bear River at Paradise, Utah']

>>> site = sites[codes[0][0]]
>>> site.geo_coords
[('-111.946402', '41.718473')]
>>> site.latitudes
['41.718473']
>>> site.longitudes
['-111.946402']
>>> info = site.site_info
>>> info.notes
[]
>>> info.site_properties
{'County': 'Cache', 'PosAccuracy_m': '1', 'State': 'Utah', 'Site Comments': 'Located below county road bridge at Mendon Road crossing'}
>>> info.altname
>>> info.elevation
'1345'

>>> f.close()
>>> f = open('tests/cuahsi_example_siteinfo_multiple.xml')
>>> sites = wml(f).response

>>> sites.query_info.criteria.method_called
'GetSiteInfo'

>>> codes = sites.site_codes
>>> codes
[['USU-LBR-Mendon'], ['USU-LBR-Wellsville']]
>>> sites.site_names
['Little Bear River at Mendon Road near Mendon, Utah', 'Little Bear River near Wellsville, Utah']

>>> site = sites[codes[1][0]]
>>> site.geo_coords
[('-111.917649', '41.643457')]

>>> catalog = site[0]

>>> series = catalog[3]
>>> series.properties
{}
>>> series.begin_date_time
datetime.datetime(2007, 11, 5, 14, 30)
>>> series.end_date_time
datetime.datetime(2008, 4, 5, 20, 30)
>>> series.method_id
'23'
>>> series.source_id
'2'
>>> series.value_count
'7309'
>>> series.value_type
'Field Observation'
>>> series.name
'Turbidity'
>>> series.organization
'Utah State University Utah Water Research Laboratory'

>>> site.variable_names
['Battery voltage', 'Oxygen, dissolved percent of saturation', 'Oxygen, dissolved', 'Specific conductance, unfiltered', 'Phosphorus, total as P', 'pH, unfiltered', 'Phosphorus, total as P, filtered', 'Turbidity', 'Solids, total Suspended', 'Gage height', 'Temperature']
>>> site.variable_codes
['USU7', 'USU6', 'USU5', 'USU4', 'USU3', 'USU41', 'USU9', 'USU8', 'USU13', 'USU34', 'USU36', 'USU10', 'USU33', 'USU32', 'USU35', 'USU40', 'USU39']

>>> variable = site['USU7']
>>> variable.variable_name
'Turbidity'
>>> variable.properties
{}
>>> variable.speciation
'Not Applicable'
>>> variable.unit.name
'nephelometric turbidity units'
>>> variable.time_scale.unit.name
'second'

>>> f.close()
>>> f = open('tests/cuahsi_example_get_values.xml')
>>> series = wml(f).response

>>> series.query_info.criteria.method_called
'GetValuesForASite'

>>> series.series_names
[None, None, None, None, None, None, None, None, None, None, None, None]

>>> series.variable_names
['Battery voltage', 'Turbidity', 'Temperature', 'Gage height', 'Battery Voltage']
>>> codes = series.variable_codes
>>> codes
['USU7', 'USU6', 'USU5', 'USU4', 'USU3', 'USU9', 'USU8', 'USU13', 'USU12', 'USU11', 'USU10', 'SDSC45']

>>> var = series.get_series_by_variable(var_code=codes[3])

>>> vals = var[0].values[0]

>>> vals.get_date_values()
[(datetime.datetime(2005, 8, 5, 0, 0), '34.53'), (datetime.datetime(2005, 8, 5, 0, 30), '37.12'), (datetime.datetime(2005, 8, 5, 1, 0), '35.97'), (datetime.datetime(2005, 8, 5, 1, 30), '35.78'), (datetime.datetime(2005, 8, 5, 2, 0), '35.68'), (datetime.datetime(2005, 8, 5, 2, 30), '36.08'), (datetime.datetime(2005, 8, 5, 3, 0), '37.8'), (datetime.datetime(2005, 8, 5, 3, 30), '37.93'), (datetime.datetime(2005, 8, 5, 4, 0), '38.88'), (datetime.datetime(2005, 8, 5, 4, 30), '37.34'), (datetime.datetime(2005, 8, 5, 5, 0), '35.15'), (datetime.datetime(2005, 8, 5, 5, 30), '35.96'), (datetime.datetime(2005, 8, 5, 6, 0), '35.62'), (datetime.datetime(2005, 8, 5, 6, 30), '34.72'), (datetime.datetime(2005, 8, 5, 7, 0), '34.7'), (datetime.datetime(2005, 8, 5, 7, 30), '33.54'), (datetime.datetime(2005, 8, 5, 8, 0), '34.98'), (datetime.datetime(2005, 8, 5, 8, 30), '31.65'), (datetime.datetime(2005, 8, 5, 9, 0), '32.49'), (datetime.datetime(2005, 8, 5, 9, 30), '32.78'), (datetime.datetime(2005, 8, 5, 10, 0), '30.58'), (datetime.datetime(2005, 8, 5, 10, 30), '32.8'), (datetime.datetime(2005, 8, 5, 11, 0), '31.83'), (datetime.datetime(2005, 8, 5, 11, 30), '30.71'), (datetime.datetime(2005, 8, 5, 12, 0), '30.82'), (datetime.datetime(2005, 8, 5, 12, 30), '29.72'), (datetime.datetime(2005, 8, 5, 13, 0), '27.05'), (datetime.datetime(2005, 8, 5, 13, 30), '25.5'), (datetime.datetime(2005, 8, 5, 14, 0), '24.69'), (datetime.datetime(2005, 8, 5, 14, 30), '26.03'), (datetime.datetime(2005, 8, 5, 15, 0), '25.55'), (datetime.datetime(2005, 8, 5, 15, 30), '25.96'), (datetime.datetime(2005, 8, 5, 16, 0), '24.72'), (datetime.datetime(2005, 8, 5, 16, 30), '23.36'), (datetime.datetime(2005, 8, 5, 17, 0), '24.21'), (datetime.datetime(2005, 8, 5, 17, 30), '25.61'), (datetime.datetime(2005, 8, 5, 18, 0), '24.73'), (datetime.datetime(2005, 8, 5, 18, 30), '25.73'), (datetime.datetime(2005, 8, 5, 19, 0), '24.76'), (datetime.datetime(2005, 8, 5, 19, 30), '24.96'), (datetime.datetime(2005, 8, 5, 20, 0), '25.69'), (datetime.datetime(2005, 8, 5, 20, 30), '27.34'), (datetime.datetime(2005, 8, 5, 21, 0), '27.14'), (datetime.datetime(2005, 8, 5, 21, 30), '27.7'), (datetime.datetime(2005, 8, 5, 22, 0), '28.88'), (datetime.datetime(2005, 8, 5, 22, 30), '30.44'), (datetime.datetime(2005, 8, 5, 23, 0), '32.14'), (datetime.datetime(2005, 8, 5, 23, 30), '34.02'), (datetime.datetime(2005, 8, 6, 0, 0), '33.61')]

>>> f.close()
>>> f = open('tests/cuahsi_example_get_variables.xml')
>>> varis = wml(f).response

>>> codes = varis.variable_codes
>>> codes
['USU3', 'USU4', 'USU5', 'USU6', 'USU7', 'USU8', 'USU9', 'USU10', 'USU11', 'USU12', 'USU13', 'USU14', 'USU15', 'USU16', 'USU17', 'USU18', 'USU19', 'USU20', 'USU21', 'USU22', 'USU23', 'USU24', 'USU25', 'USU26', 'USU27', 'USU28', 'USU29', 'USU30', 'USU31', 'USU32', 'USU33', 'USU34', 'USU35', 'USU36', 'USU37', 'USU38', 'USU39', 'USU40', 'USU41', 'USU42', 'USU43', 'SDSC45']
>>> varis.variable_names
['Battery voltage', 'Temperature', 'Gage height', 'Solids, total Suspended', 'Specific conductance, unfiltered', 'Oxygen, dissolved', 'Wind direction', 'Oxygen, dissolved percent of saturation', 'Radiation, incoming shortwave', 'Relative humidity', 'Phosphorus, total as P', 'pH, unfiltered', 'Battery Voltage', 'Phosphorus, total as P, filtered', 'Barometric pressure', 'Discharge', 'Turbidity', 'Precipitation', 'Wind speed']
>>> var = varis[codes[10]]
>>> var.variable_name
'Gage height'
>>> var.no_data_value
'-9999'
>>> var.properties
{}
>>> var.unit.name
'international foot'

#geopython

Author:	Tom Kralidis
Contact:	tomkralidis at hotmail.com
Release:	0.7-dev
Date:	2013-04-12

Standard	Version(s)
OGC WMS	1.1.1
OGC WFS	1.0.0, 1.1.0, 2.0.0
OGC WCS	1.0.0, 1.1.0
OGC WMC	1.1.0
OGC SOS	1.0.0, 2.0.0
OGC SensorML	1.0.1
OGC CSW	2.0.2
OGC WPS	1.0.0
OGC Filter	1.1.0
OGC OWS Common	1.0.0, 1.1.0, 2.0
NASA DIF	9.7
FGDC CSDGM	1998
ISO 19139	2007
Dublin Core	1.1
WMTS	1.0.0
WaterML	1.0, 1.1

OWSLib 0.7-dev documentation¶

Introduction¶

Features¶

Standards Support¶

Installation¶

Requirements¶

Install¶

Usage¶

WMS¶

WFS¶

WCS¶

CSW¶

WMC¶

WPS¶

SOS¶

SensorML¶

ISO¶

CRS Handling¶

Dublin Core¶

NASA DIF¶

FGDC¶

WMTS¶

WaterML¶

Development¶

Testing¶

Support¶

Mailing Lists¶

Submitting Questions to Community¶

Searching the Archives¶

Metrics¶

IRC¶

License¶

Credits¶

Table Of Contents

This Page

Navigation

OWSLib 0.7-dev documentation¶

Introduction¶

Features¶

Standards Support¶

Installation¶

Requirements¶

Install¶

Usage¶

WMS¶

WFS¶

WCS¶

CSW¶

WMC¶

WPS¶

SOS¶

SensorML¶

ISO¶

CRS Handling¶

Dublin Core¶

NASA DIF¶

FGDC¶

WMTS¶

WaterML¶

Development¶

Testing¶

Support¶

Mailing Lists¶

Submitting Questions to Community¶

Searching the Archives¶

Metrics¶

IRC¶

License¶

Credits¶

Table Of Contents

This Page

Quick search

Navigation