def read_dicom(dicom_path): # Function that reads a dicom file and writes to a text file dataset = ReadData.read_dicom(dicom_path) ds = dataset.DeformableRegistrationSequence[1].DeformableRegistrationGridSequence[0].VectorGridData ds = np.array(ds).astype(np.float64) print(ds) with open("dvf.raw", "wb") as f: f.write(ds)
def read_dicom(dicom_path): # Function that reads a dicom file and writes to a text file dataset = ReadData.read_dicom(dicom_path) # vector_grid = dataset.DeformableRegistrationSequence[1].DeformableRegistrationGridSequence[0].VectorGridData # vector_grid = np.array(vector_grid).astype(np.float64) # # with open("dvf.raw", "wb") as f: # f.write(vector_grid) ReadData.write_dicom(dataset, "image_test")
def get_rigid_transforms(reg_path): dataset = ReadData.read_dicom(reg_path) # Pre-Deformation pre_def = dataset.DeformableRegistrationSequence[ 1].PreDeformationMatrixRegistrationSequence[ 0].FrameOfReferenceTransformationMatrix pre_def = np.array(pre_def).astype(np.float64) pre_def = np.reshape(pre_def, (4, 4)) # Post-Deformation post_def = dataset.DeformableRegistrationSequence[ 1].PostDeformationMatrixRegistrationSequence[ 0].FrameOfReferenceTransformationMatrix post_def = np.array(post_def).astype(np.float64) post_def = np.reshape(post_def, (4, 4)) return pre_def, post_def
def read_dicom(dicom_path): ds = ReadData.read_dicom(dicom_path) #petct_data = # dvf_data = ReadData.load_dvf_data(ds) ReadData.write_dicom(ds, "pct_dicom")