[python] array 값을 dicom 값으로 바꾸는 코드!

import pydicom
import numpy as np
from pydicom.dataset import Dataset
import glob

img=np.fromfile(data_dir, dtype='int16', sep="")


# DICOM 파일 생성
dicom_dataset = Dataset()
dicom_dataset.PatientName = "John Doe"
dicom_dataset.PatientID = "123456"
dicom_dataset.Modality = "CT"

# 이미지 데이터 설정
dicom_dataset.PixelData = img.tobytes()
dicom_dataset.Rows, dicom_dataset.Columns = img.shape
dicom_dataset.BitsAllocated = 16
dicom_dataset.BitsStored = 16
dicom_dataset.HighBit = 15
dicom_dataset.PixelRepresentation = 0  # unsigned integer
dicom_dataset.SOPInstanceUID = pydicom.uid.generate_uid()
dicom_dataset.PhotometricInterpretation = "MONOCHROME1"


# 필수 속성 추가
dicom_dataset.is_little_endian = True
dicom_dataset.is_implicit_VR = False

# DICOM 파일 저장
dicom_dataset.save_as("output_dicom.dcm")

 

array를 dicom파일로 바꾸는 코드입니다 :)