print metadata in YAML format and fix parsing of attribute names; update README examples with better coverages and formatting for YAML output

Author: misev

README.md

@@ -30,8 +30,10 @@ object. Optionally, a username and a password may be specified, if the endpoint
 ```python
 from wcs.service import WebCoverageService
 
-wcs_endpoint = "https://ows.rasdaman.org/rasdaman/ows"
+wcs_endpoint = "https://fairicube.rasdaman.com/rasdaman/ows"
 service = WebCoverageService(wcs_endpoint)
+
+# or with credentials:
 # service = WebCoverageService(wcs_endpoint,
 #                              username=..., password=...)
 ```
@@ -44,88 +46,94 @@ with basic information such as a WGS 84 bounding box and a native bounding box:
 coverages = service.list_coverages()
 ```
 
-Let's inspect a single coverage with name `AvgLandTemp`:
+Let's print information of a single coverage with name `dominant_leaf_type_20m`:
 
 ```python
-avg_land_temp = coverages['AvgLandTemp']
+cov = coverages['dominant_leaf_type_20m']
 
-# print all information
+print(cov)
+```
+```yaml
+dominant_leaf_type_20m:
+  subtype: ReferenceableGridCoverage
+  native CRS: OGC:AnsiDate+EPSG:3035
+  geo bbox:
+    time:
+      min: "2012-01-01"
+      max: "2015-01-01"
+      crs: OGC:AnsiDate
+    Y:
+      min: 900000
+      max: 5500000
+      crs: EPSG:3035
+    X:
+      min: 900000
+      max: 7400000
+      crs: EPSG:3035
+  lon/lat bbox:
+    Lon:
+      min: -56.50514190170437
+      max: 72.9061049341568
+      crs: EPSG:4326
+    Lat:
+      min: 24.28417068794856
+      max: 72.66326966834436
+      crs: EPSG:4326
+  size in bytes: 113000000001
+  additional params:
+    title: Dominant Leaf Type (2012-2015)
+    sizeInBytesWithPyramidLevels: "113000000001"
+```
 
-print(avg_land_temp)
-
-# AvgLandTemp:
-#   subtype: ReferenceableGridCoverage
-#   native CRS: OGC:AnsiDate+EPSG:4326
-#   geo bbox:
-#     ansi:
-#       min: "2000-02-01"
-#       max: "2015-06-01"
-#       crs: OGC:AnsiDate
-#     Lat:
-#       min: -90
-#       max: 90
-#       crs: EPSG:4326
-#     Lon:
-#       min: -180
-#       max: 180
-#       crs: EPSG:4326
-#   lon/lat bbox:
-#     Lon:
-#       min: -180
-#       max: 180
-#       crs: EPSG:4326
-#     Lat:
-#       min: -90
-#       max: 90
-#       crs: EPSG:4326
-#   size in bytes: 4809618404
+Examples for extracting individual details of the coverage:
 
+```python
 # coverage subtype
 
-print(avg_land_temp.subtype)
+print(cov.subtype)
 
 # ReferenceableGridCoverage
 
 # coverage bounding box, containing the CRS and axes
 
-bbox = avg_land_temp.bbox
+bbox = cov.bbox
 
 # full coverage crs identifier
 
 print(bbox.crs)
 
 # https://www.opengis.net/def/crs-compound?
-# 1=https://www.opengis.net/def/crs/OGC/0/AnsiDate&
-# 2=https://www.opengis.net/def/crs/EPSG/0/4326
+# 1=https://www.opengis.net/def/crs/OGC/0/AnsiDate?axis-label="time"&
+# 2=https://www.opengis.net/def/crs/EPSG/0/3035
 
 # coverage crs identifier in shorthand notation
 
 from wcs.model import Crs
+
 print(Crs.to_short_notation(bbox.crs))
 
-# OGC:AnsiDate+EPSG:4326
+# OGC:AnsiDate+EPSG:3035
 
 # get information for the first axis; as it is a temporal axis,
 # the lower_bound and upper_bound are datetime.datetime objects.
 
-axis = bbox.ansi
+axis = bbox.time
 
 # note that these are all equivalent:
-# axis = bbox['ansi']
-# axis = bbox.0
+# axis = bbox['time']
 # axis = bbox[0]
 
 name = axis.name
 lower_bound = axis.low
 upper_bound = axis.high
 print(f'{name}({lower_bound} - {upper_bound})')
-# ansi(2000-02-01 00:00:00+00:00 - 2015-06-01 00:00:00+00:00)
+# time(2012-01-01 00:00:00+00:00 - 2015-01-01 00:00:00+00:00)
 
 # get size in bytes if available
 
-if avg_land_temp.size_bytes is not None:
-    print(avg_land_temp.size_bytes)
-    # 4809618404
+if cov.size_bytes is not None:
+    print(cov.size_bytes)
+    # 113000000001
 ```
 
 ## Full Coverage Information
@@ -139,64 +147,94 @@ More detailed information can be retrieved with the
 object:
 
 ```python
-full_avg_land_temp = service.list_full_info('AvgLandTemp')
+cov = service.list_full_info('dominant_leaf_type_20m')
 
 # print all information
 
-print(full_avg_land_temp)
-
-# AvgLandTemp:
-#   native CRS: OGC:AnsiDate+EPSG:4326
-#   geo bbox:
-#     ansi:
-#       min: "2000-02-01"
-#       max: "2015-06-01"
-#       crs: OGC:AnsiDate
-#       uom: d
-#       type: irregular
-#       coefficients: ["2000-02-01", "2000-03-01", ...
-#                      "2015-05-01", "2015-06-01"]
-#     Lat:
-#       min: -90
-#       max: 90
-#       crs: EPSG:4326
-#       uom: degree
-#       resolution: -0.1
-#       type: regular
-#     Lon:
-#       min: -180
-#       max: 180
-#       crs: EPSG:4326
-#       uom: degree
-#       resolution: 0.1
-#       type: regular
-#   grid bbox:
-#     i:
-#       min: 0
-#       max: 184
-#       resolution: 1
-#       type: regular
-#     j:
-#       min: 0
-#       max: 1799
-#       resolution: 1
-#       type: regular
-#     k:
-#       min: 0
-#       max: 3599
-#       resolution: 1
-#       type: regular
-#   range type fields:
-#     Gray:
-#       type: Quantity
-#       label: Gray
-#       definition: http://www.opengis.net/def/dataType/OGC/0/float32
-#       nil values: 99999
-#       uom: 10^0
-#   metadata:
-#     {
-#       "covMetadata": null
-#     }
+print(cov)
+```
+```yaml
+dominant_leaf_type_20m:
+  native CRS: OGC:AnsiDate+EPSG:3035
+  geo bbox:
+    time:
+      min: 2012-01-01
+      max: 2015-01-01
+      crs: OGC:AnsiDate
+      uom: d
+      type: irregular
+      coefficients:
+        - 2012-01-01
+        - 2015-01-01
+    Y:
+      min: 900000
+      max: 5500000
+      crs: EPSG:3035
+      uom: metre
+      resolution: -20
+      type: regular
+    X:
+      min: 900000
+      max: 7400000
+      crs: EPSG:3035
+      uom: metre
+      resolution: 20
+      type: regular
+  grid bbox:
+    i:
+      min: 0
+      max: 1
+      resolution: 1
+      type: regular
+    j:
+      min: -125000
+      max: 104999
+      resolution: 1
+      type: regular
+    k:
+      min: 0
+      max: 324999
+      resolution: 1
+      type: regular
+  range type fields:
+    dlt:
+      type: Category
+      label: dominant leaf type map of Europe
+      description: >
+        raster coding (thematic pixel values): 
+        0: all non-tree covered areas; 
+        1: broadleaved trees; 
+        2: coniferous trees; 
+        254: unclassifiable (no satellite image available, or clouds, shadows, or snow); 
+        255: outside area
+      definition: https://land.copernicus.eu/en/technical-library/hrl-forest-2012-2015/@@download/file
+      nil values: 250
+  metadata:
+    covMetadata:
+      axes:
+        time:
+          areasOfValidity:
+            area:
+              - 
+                "@start": 2011-01-01T00:00:00.000Z
+                "@end": 2013-12-31T23:59:59.999Z
+              - 
+                "@start": 2014-01-01T00:00:00.000Z
+                "@end": 2016-12-31T23:59:59.999Z
+    rasdamanCoverageMetadata:
+      catalog:
+        title: Dominant Leaf Type (2012-2015)
+        thumbnail: https://fairicube.rasdaman.com/rasdaman/ows/coverage/thumbnail?COVERAGEID=dominant_leaf_type_20m
+        description: Provides at pan-European level in the spatial resolution of 20 m information on the dominant leaf type (broadleaved or coniferous).
+        provenance:
+          "@sourceUrl": https://land.copernicus.eu/en/products/high-resolution-layer-dominant-leaf-type
+          "@providerName": Copernicus
+          "@termsUrl": https://land.copernicus.eu/en/data-policy
+        ourTerms: https://fairicube.rasdaman.com/#terms
+    fairicubeMetadata:
+      "@role": https://codelists.fairicube.eu/metadata/MetadataCatalogLink
+      "@title": Metadata in the FAIRiCUBE Catalog
+      "@href": https://stacapi.eoxhub.fairicube.eu/collections/index/items/dominant_leaf_type_20m
 ```
 
 In addition to the geo `bbox` in native CRS, the 
@@ -207,7 +245,7 @@ type of
 object, except its `crs` attribute is `None`.
 
 ```python
-print(full_avg_land_temp.grid_bbox)
+print(cov.grid_bbox)
 ```
 
 The `range_type` attribute indicates the structure of the cell values
@@ -220,30 +258,33 @@ coverage. Check the documentation of
 for full details.
 
 ```python
-range_type = full_avg_land_temp.range_type
+range_type = cov.range_type
 all_fields = range_type.fields
-field = range_type.Gray
+dlt = range_type.dlt
 
 # note that these are all equivalent:
-# field = range_type['Gray']
-# field = range_type.0
+# field = range_type['dlt']
 # field = range_type[0]
 
 # get all properties of the field
 
-label = field.label
-description = field.description
-definition = field.definition
-nil_values = field.nil_values
-if field.is_quantity:
-  uom = field.uom
+label = dlt.label
+description = dlt.description
+definition = dlt.definition
+nil_values = dlt.nil_values
+if dlt.is_quantity:
+    uom = dlt.uom
 else:
-  codespace = field.codespace
+    codespace = dlt.codespace
 ```
 
 Finally, any coverage metadata is available from the `metadata` attribute,
 which is a nested dict mirroring the XML structure in the *DescribeCoverage* document.
+E.g. to get the link to the FAIRiCUBE catalog entry for this coverage:
 
+```python
+catalog_link = cov.metadata['fairicubeMetadata']['@href']
+```
 
 # Contributing
 

tests/test_service.py

@@ -47,6 +47,18 @@ def test_list_coverages():
     assert cov.bbox.ansi.crs == 'https://www.opengis.net/def/crs/OGC/0/AnsiDate'
 
 
+def test_list_coverages_only_local():
+    service = WebCoverageService("https://fairicube.rasdaman.com/rasdaman/ows")
+    coverages = service.list_coverages(only_local=True)
+    assert all(cov.is_local() for k, cov in coverages.items())
+
+
+def test_list_coverages_all():
+    service = WebCoverageService("https://fairicube.rasdaman.com/rasdaman/ows")
+    coverages = service.list_coverages()
+    assert not all(cov.is_local() for k, cov in coverages.items())
+
+
 def test_list_full_info():
     service = WebCoverageService("https://ows.rasdaman.org/rasdaman/ows")
     cov = service.list_full_info('AvgLandTemp')
@@ -124,9 +136,15 @@ def test_list_full_info():
       nil values: 99999
       uom: 10^0
   metadata:
-    {
-      "covMetadata": null
-    }'''
+    covMetadata: None
+'''
     assert str(cov) == expected
     subset = cov.bbox.ansi["2006-08-01" : "2007-01-01"]
     assert len(subset) == 6
+    assert cov.is_local()
+
+def test_list_full_info2():
+    service = WebCoverageService("https://fairicube.rasdaman.com/rasdaman/ows")
+    cov = service.list_full_info('dominant_leaf_type_20m')
+    catalog_link = cov.metadata['fairicubeMetadata']['@href']
+    assert catalog_link == "https://stacapi.eoxhub.fairicube.eu/collections/index/items/dominant_leaf_type_20m"

wcs/__init__.py

@@ -5,4 +5,4 @@
 """
 __title__ = 'wcs'
 __author__ = 'rasdaman team'
-__version__ = "0.2.2"
+__version__ = "0.2.3"