I’m using the SimpleITK API for python to load medical images in DICOM format (PET and CT). However, it appears that SimpleITK doesn’t apply the rescale slope and rescale intercept. My code right now is:
import SimpleITK as sitk root_pet = 'path_dicom_files' reader = sitk.ImageSeriesReader() dicom_names = reader.GetGDCMSeriesFileNames(root_pet) reader.SetFileNames(dicom_names) reader.MetaDataDictionaryArrayUpdateOn() reader.LoadPrivateTagsOn() imagePET = reader.Execute() pet_array = sitk.GetArrayFromImage(imagePET)
However, the range is: [-31796.78, 31074.89], but it should be [0, 4883997.44]. Can anyone help as to why the PET image is not rescaled? It seems like it’s still in the int32 format right now.