20 import SimpleITK
as sitk
30 +
" <input_directory_with_DICOM_series> <output_directory>"
36 data_directory = sys.argv[1]
37 series_IDs = sitk.ImageSeriesReader.GetGDCMSeriesIDs(data_directory)
40 'ERROR: given directory "'
42 +
'" does not contain a DICOM series.'
45 series_file_names = sitk.ImageSeriesReader.GetGDCMSeriesFileNames(
46 data_directory, series_IDs[0]
50 series_reader.SetFileNames(series_file_names)
57 series_reader.MetaDataDictionaryArrayUpdateOn()
58 series_reader.LoadPrivateTagsOn()
59 image3D = series_reader.Execute()
74 writer.KeepOriginalImageUIDOn()
90 modification_time = time.strftime(
"%H%M%S")
91 modification_date = time.strftime(
"%Y%m%d")
104 direction = filtered_image.GetDirection()
105 series_tag_values = [
106 (k, series_reader.GetMetaData(0, k))
107 for k
in tags_to_copy
108 if series_reader.HasMetaDataKey(0, k)
110 (
"0008|0031", modification_time),
111 (
"0008|0021", modification_date),
112 (
"0008|0008",
"DERIVED\\SECONDARY"),
115 "1.2.826.0.1.3680043.2.1125." + modification_date +
".1" + modification_time,
137 series_reader.GetMetaData(0,
"0008|103e")
138 if series_reader.HasMetaDataKey(0,
"0008|103e")
139 else "" +
" Processed-SimpleITK"
144 for i
in range(filtered_image.GetDepth()):
145 image_slice = filtered_image[:, :, i]
147 for tag, value
in series_tag_values:
148 image_slice.SetMetaData(tag, value)
151 image_slice.SetMetaData(
"0008|0012", time.strftime(
"%Y%m%d"))
153 image_slice.SetMetaData(
"0008|0013", time.strftime(
"%H%M%S"))
155 image_slice.SetMetaData(
157 "\\".join(map(str, filtered_image.TransformIndexToPhysicalPoint((0, 0, i)))),
160 image_slice.SetMetaData(
"0020|0013", str(i))
164 writer.SetFileName(os.path.join(sys.argv[2], str(i) +
".dcm"))
165 writer.Execute(image_slice)