Posts by marcelvanherk

##
#  @file       simpleDicomQR.py
#  @author     Gareth Price
#  @date       23/9/2016
#  @version    0.01
#
#  @brief      Classes to perform operations on remote PACS using PyNetDicom libraries
#
#  Copyright (C) 2016 The Christie NHS
#  Foundation Trust
##

#  VERSION HISTORY
#  Version:    0.01
#  Date:       23/8/2016
#  Update:     Origin
#  Author:     Gareth Price
#
#  Version:    0.02
#  Date:       06/01/2021
#  Update:     Extracted from classes to simple single file runtime
#  Author:     Gareth Price

#NOTE - This code was writen in 2016 against older libraries. The notes below show the versions of pydicom and pynetdicom
#it was writen against
# - PyDicom 1.0.0a - There are many code changes required between 0.99 and 1.0 versions of pydicom, comes from github
#latest - not an official release https://github.com/darcymason/pydicom
# - PyNetDicom - 0.9.0 - The official PyPI distribution is 2.7 only. The 3.5 port (with significantly different syntax)
#comes from a fork - https://github.com/scaramallion/pynetdicom3
#Both of these are cloned to Downloads and then installed via setup.py

#Using a standard version of conquest 1.5.0b - naming as in distribution localhost, CONQUESTSRV1, 5678

import pydicom
import pynetdicom3
from pynetdicom3.pdu_primitives import SCP_SCU_RoleSelectionNegotiation

# ===================================================================================================
# Retrieve image data from PACS
# ===================================================================================================
class pacsRetrieveClass():

    """Retrieve data from remote PACS"""

    def __init__(self):
        self.connected = False

    def connectPACS(self):
        """"Connect to remote PACS"""
        self.AE = pynetdicom3.AE('localdicom', 9999, scu_sop_class=[pynetdicom3.VerificationSOPClass]+pynetdicom3.QueryRetrieveSOPClassList+pynetdicom3.StorageSOPClassList,
                                 scp_sop_class=[pynetdicom3.VerificationSOPClass]+pynetdicom3.StorageSOPClassList,
                                 transfer_syntax=[pydicom.uid.ImplicitVRLittleEndian,pydicom.uid.ExplicitVRLittleEndian,
                                        pydicom.uid.ExplicitVRBigEndian])  # netdicom local dicom SCU/SCP
        ext_neg = []
        for context in self.AE.presentation_contexts_scu:
            tmp = SCP_SCU_RoleSelectionNegotiation()
            tmp.sop_class_uid = context.AbstractSyntax
            tmp.scu_role = False
            tmp.scp_role = True
            ext_neg.append(tmp)

        self.AE.OnAssociateResponse = self.onAssociateResponse
        self.AE.OnAssociateRequest = self.onAssociateRequest
        self.association = self.AE.associate('localhost', 5678, 'CONQUESTSRV1',
                                             ext_neg=ext_neg)  # Request association with remote SCP/SCU
        self.connected = True

    def onAssociateResponse(self,association):
        print("Association response received")

    def onAssociateRequest(self,association):
        print("Association resquested")
        return True

    def isConnected(self):
        if(not hasattr(self,'association')):
            self.connected = False
            return False
        if(self.association.is_alive()):
            self.connected = True
            return True
        else:
            self.connected = False
            return False

    def echoRemotePACS(self):
        """Simple Echo SCU implementation"""
        if(not self.connected):
            self.connectPACS()
        if self.association.is_established:
            status = self.association.send_c_echo()
            if('Success' in status.status_type):
                self.connected = True
                print("Echo sucessful")
                return True
            else:
                self.connected = False
                print("Echo failed")
                return False

    def queryRemotePACS(self):
        """CFind to get all patient names on the server (ie just the HEAD EXP2 image included with conquest"""
        self.isConnected()
        if (not self.connected):
            self.connectPACS()

        #Build query
        self.query = pydicom.dataset.Dataset()  # Default query for all patients
        self.query.PatientID = ''  # Add additional data elements as required
        self.query.PatientName = ''
        self.query.SeriesInstanceUID = ''
        self.query.QueryRetrieveLevel = "PATIENT"  # Change level as required

        if self.association.is_established:
            print("DICOM query all patients")
            responses = self.association.send_c_find(self.query,query_model='P')  # Patient Root informmation model
            patientIDs = []
            patientNames = []
            for (status, dataset) in responses:
                if (hasattr(dataset, 'PatientID')):
                    patientIDs.append(dataset.PatientID)
                if (hasattr(dataset, 'PatientName')):
                    patientNames.append(dataset.PatientName)
        print("Find query returned patinet IDs and names")
        print(patientIDs)
        print(patientNames)

    def onStore(self,dataset):
        """Method to call when acting as store SCP """

        meta = pydicom.dataset.Dataset()
        meta.TransferSyntaxUID = pydicom.uid.ImplicitVRLittleEndian #Need to explicitly add file metadata and Transfer Syntax to get pixel uncompression to work in pydicom
        meta.MediaStorageSOPClassUID = dataset.SOPClassUID
        meta.MediaStorageSOPInstanceUID = dataset.SOPInstanceUID
        meta.ImplementationClassUID = pynetdicom3.pynetdicom_uid_prefix

        dataset.file_meta = meta    #Need to explicitly add file metadata and Transfer Syntax to get pixel uncompression to work in pydicom
        dataset.is_little_endian = True
        dataset.is_implicit_VR = True


        print("Saving data retrieved via CGet")
        dataset.save_as('E:\\Temp\\dicomExampleCGetData'+str(self.imgNbr)+'.dcm',write_like_original=False)
        self.imgNbr=self.imgNbr+1

        return 0x0000   #Return success

    def retrieveSeriesCGet(self):
        """Given patient ID retrive images"""
        self.isConnected()
        if (not self.connected):
            self.connectPACS()

        self.imgNbr=0   #Update to allow multiple images to be saved with different filenames
        self.AE.on_c_store = self.onStore   #Assign cStore method

        #Query data structure
        self.retrieveQuery = pydicom.dataset.Dataset()
        self.retrieveQuery.QueryRetrieveLevel = "PATIENT"
        self.retrieveQuery.PatientID = '0009703828'

        if self.association.is_established:
            print("DICOM retrieve patient ")
            responses = self.association.send_c_get(self.retrieveQuery, query_model='P')
            if responses is not None:
                for response in responses:
                    pass

#===================================================================================================
#Main
#===================================================================================================
if __name__=="__main__":    #Test for pacsClasses

    pacs = pacsRetrieveClass()
    pacs.connectPACS()
    pacs.echoRemotePACS()
    pacs.queryRemotePACS()
    pacs.retrieveSeriesCGet()