Я пытаюсь объединить различные файлы DICOM в одну многослойную серию, чтобы визуализировать их в ITK-Snap. Тем не менее, я не могу получить работающую серию DICOM.
Я отсортировал все файлы в зависимости от расположения срезов, и у меня есть несколько упорядоченных отдельных файлов .dcm с их исходной информацией. Я заменил все их исходные UID экземпляров серии одним идентификатором, а их номер серии - одним пользовательским номером серии, который я установил на «999» (чтобы они принадлежали одной серии). Ориентация изображения установлена на [1;0;0;0;1;0]
для всех файлов и толщина среза 8 мм для всех файлов.
Затем я создал массив информационных структур с исходным фрагментом positionings [info(num)]
.
Я пробовал что-то вроде:
for i=1:num %where num is the number of dicom files
k = num2str(i);
dicomwrite(imm,k,info(i),'CreateMode','Copy'); %where imm is the matrix I obtained with dicomread
end
Я получил новый набор файлов dicom, названных числами от 1 до num, однако, когда я пытаюсь открыть серию на ITK-snap, возникает исключение, указывающее, что вектор слишком длинный. Я могу открыть отдельные файлы dicom на ITK-snap, однако, когда в серию входит более одного изображения, а серия визуализируется как 256x212xnum, где num - количество файлов, я сталкиваюсь с исключением.
Что я делаю не так?