aboutsummaryrefslogtreecommitdiff
path: root/bh20sequploader/qc_metadata.py
blob: 87fa0b36be05d96312bfaa349abf5f99fdae8ef1 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
import schema_salad.schema
import schema_salad.ref_resolver
import schema_salad.jsonld_context
import logging
import pkg_resources
import logging
import traceback
from rdflib import Graph, Namespace
from pyshex.evaluate import evaluate

metadata_schema = None

def qc_metadata(metadatafile):
    global metadata_schema
    log = logging.getLogger(__name__ )
    if metadata_schema is None:
        schema_resource = pkg_resources.resource_stream(__name__, "bh20seq-schema.yml")
        cache = {"https://raw.githubusercontent.com/arvados/bh20-seq-resource/master/bh20sequploader/bh20seq-schema.yml": schema_resource.read().decode("utf-8")}
        metadata_schema = schema_salad.schema.load_schema("https://raw.githubusercontent.com/arvados/bh20-seq-resource/master/bh20sequploader/bh20seq-schema.yml", cache=cache)

    (document_loader,
     avsc_names,
     schema_metadata,
     metaschema_loader) = metadata_schema

    shex = pkg_resources.resource_stream(__name__, "bh20seq-shex.rdf").read().decode("utf-8")

    if not isinstance(avsc_names, schema_salad.avro.schema.Names):
        raise Exception(avsc_names)

    doc, metadata = schema_salad.schema.load_and_validate(document_loader, avsc_names, metadatafile, True)
    g = schema_salad.jsonld_context.makerdf("workflow", doc, document_loader.ctx)
    rslt, reason = evaluate(g, shex, doc["id"], "https://raw.githubusercontent.com/arvados/bh20-seq-resource/master/bh20sequploader/bh20seq-shex.rdf#submissionShape")

    # As part of QC make sure serialization works too, this will raise
    # an exception if there are invalid URIs.
    g.serialize(format="ntriples")

    if not rslt:
        raise Exception(reason)

    return metadata['sample']['sample_id']