94t billing

94t-billing/.gitignore

@@ -0,0 +1,3 @@
+.ipyn*
+*.ipynb
+netrc

94t-billing/get-study-info.py

@@ -0,0 +1,65 @@
+#!/usr/bin/env python3
+
+import sys
+import datetime
+import netrc
+import requests
+import pyorthanc
+
+host = "micvna.mclean.harvard.edu"
+port = 8042
+server = f"http://{host}:{port}"
+username, _, password = netrc.netrc("netrc").authenticators(host)
+o = pyorthanc.Orthanc(server, username=username, password=password)
+
+if len(sys.argv) < 2:
+    print("Usage: %s <accession number>" % sys.argv[0])
+    sys.exit(1)
+
+accnum = sys.argv[1]
+query = {"AccessionNumber": accnum}
+studies = pyorthanc.find_studies(client=o, query=query)
+study = studies[0]
+
+if len(studies) == 0:
+    print(f"No studies found with query {repr(query)}")
+    sys.exit(1)
+
+
+def duration(study):
+    # use bulk-content as it's vastly faster than the pyorthanc method
+
+    datetime_format = "%Y%m%d%H%M%S.%f"
+
+    data = requests.post(
+        server + "/tools/bulk-content",
+        json={"Resources": [study.identifier], "Level": "Instance"},
+        auth=(username, password),
+    ).json()
+
+    instance_creation_datetimes = [
+        datetime.datetime.strptime(
+            item["MainDicomTags"]["InstanceCreationDate"]
+            + item["MainDicomTags"]["InstanceCreationTime"],
+            datetime_format,
+        )
+        for item in data
+        if "MainDicomTags" in item
+        and "InstanceCreationDate" in item["MainDicomTags"]
+        and "InstanceCreationTime" in item["MainDicomTags"]
+    ]
+
+    return max(instance_creation_datetimes) - min(instance_creation_datetimes)
+
+
+data = {}
+data["date"] = study.date
+data["AccessionNumber"] = study.main_dicom_tags["AccessionNumber"]
+data["protocol"] = study.main_dicom_tags["StudyDescription"]
+data["duration"] = str(duration(study)).split(".")[0]
+
+# print the data as a CSV line
+import csv
+
+writer = csv.writer(sys.stdout)
+writer.writerow(data.values())