define REC relationships

bricksrc/relationships.py

@@ -1,5 +1,6 @@
 from rdflib import Literal
-from .namespaces import A, OWL, RDFS, BRICK, VCARD, QUDT, SDO, RDF, BSH, XSD
+from .namespaces import A, OWL, RDFS, BRICK, VCARD, QUDT, SDO, RDF, BSH, XSD, REC
+from .env import env
 
 """
 Defining Brick relationships
@@ -194,3 +195,19 @@
         RDFS.label: Literal("is sub-meter of", lang="en"),
     },
 }
+
+# add REC relationships by mining them from the REC ontology
+rec = env.get_graph("https://w3id.org/rec")
+# for all objects of sh:path, read out the sh:nodeKind and sh:datatype
+query = """SELECT ?path ?nodeKind ?datatype WHERE {
+    ?p sh:path ?path .
+    OPTIONAL { ?p sh:nodeKind ?nodeKind }
+    OPTIONAL { ?p sh:datatype ?datatype }
+}"""
+for row in rec.query(query):
+    if row["path"] not in relationships:
+        relationships[row["path"]] = {}
+    if row["datatype"]:
+        relationships[row["path"]][A] = [OWL.DatatypeProperty]
+    if row["nodeKind"]:
+        relationships[row["path"]][A] = [OWL.ObjectProperty, OWL.AsymmetricProperty, OWL.IrreflexiveProperty]

generate_brick.py

@@ -631,12 +631,13 @@ def define_relationships(definitions, superprop=None, graph=G):
         return
 
     for prop, propdefn in definitions.items():
-        if isinstance(prop, str):
+        if not isinstance(prop, URIRef):
             prop = BRICK[prop]
         if superprop is not None:
             graph.add((prop, RDFS.subPropertyOf, superprop))
 
-        graph.add((prop, RDFS.subPropertyOf, BRICK.Relationship))
+        if prop.startswith(BRICK):
+            graph.add((prop, RDFS.subPropertyOf, BRICK.Relationship))
 
         # define property types
         prop_types = propdefn.get(A, [])
@@ -650,7 +651,8 @@ def define_relationships(definitions, superprop=None, graph=G):
         define_relationships(subproperties_def, prop, graph=graph)
 
         # generate a SHACL Property Shape for this relationship
-        propshape = BSH[f"{prop.split('#')[-1]}Shape"]
+        qname = graph.namespace_manager.qname(prop)
+        propshape = BSH[f"{qname.replace(':','_')}Shape"]
         graph.add((propshape, A, SH.PropertyShape))
         graph.add((propshape, SH.path, prop))
         if "range" in propdefn.keys():