Все данные получены от клинического пациента, чья операция SDR была проведена в больнице BJ Dongzhimen. Протокол соответствует руководящим принципам и был одобрен комитетом по этике исследований больницы Дунчжимэнь.
ПРИМЕЧАНИЕ: Вся карта протокола реконструкции модели показана на рисунке 1. Данные компьютерной томографии высокого разрешения (HRCT) и данные Диксона являются сырьем для моделирования; затем создание 3D-модели состоит из регистрации изображения и слияния. Окончательная цифровая 3D-модель печатается по технологии PolyJet, которая представляет собой высокоточный процесс 3D-печати, позволяющий получать гладкие и точные детали с использованием широкого спектра материалов. Для точного описания пространственных отношений между позвонком и спинномозговым нервом используются данные HRCT и серия изображений Диксона. Сканирование Диксона может идентифицировать изображения разделения воды и жира, в которых серия изображений фазы воды Диксона может быть использована для извлечения структуры спинномозговых нервов, а серия изображений Диксона в фазе может быть использована для проверки регистрации структуры кости.

Рисунок 1: Вся карта протокола. Методология исследования заключается в слиянии последовательностей КТ и магнитно-резонансной последовательности Диксона. В частности, структура КТ-позвонков регистрируется с идентичной структурой позвонков, содержащейся в последовательности Диксона-в, с последующим слиянием с последовательностью Диксона-w для спинномозгового нерва. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
1. Сбор и подготовка данных
- КТ высокого разрешения для позвонков
ПРИМЕЧАНИЕ: Разница параметров не чувствительна к методу исследования.- Установите ресурсы данных со станции КТ.
ПРИМЕЧАНИЕ: Здесь используется компьютерная томография SIEMENS-CTAWP73396.
- Откройте программное обеспечение Syngo CT 2012B , чтобы получить данные из протокола сканирования SpineRoutine_1. Выберите размер пикселя и толщину среза (ST) набора данных, чтобы адаптироваться к размеру позвонков, предназначенных для представления в цифровой 3D-модели.
- Используйте ST 1 мм с размером матрицы 512 пикселей x 512 пикселей, в котором интервал между пикселями составляет 0,3320 мм. Фактический размер достигнутого 3D-объема составляет 512 x 512 x 204 вокселей.
- Последовательность Диксона для спинномозгового нерва
ПРИМЕЧАНИЕ: В этом исследовании используется аппарат МРТ 1,5 Тл.- Для получения точных данных установите разрешение изображения Диксона 290 x 320 пикселей, расстояние между пикселями 0,9375 мм и толщину среза 3 мм.
- Установите время повторения 5,160 мс и время эха 94 мс.
- Убедитесь, что каждый отсканированный слой состоит из четырехфазных изображений: Dixon-in, Dixon-opp, Dixon-F и Dixon-w.
- Подготовьте файлы хранения данных для реконструкции модели.
ПРИМЕЧАНИЕ: Четко определенная структура хранения данных более удобна для последующей работы.- Создайте папку проекта, содержащую все данные, принадлежащие пациенту.
- Подготовьте различные пути к файлам для данных HRCT и MRI-Dixon, создав разные папки для данных цифровой визуализации и коммуникаций в медицине (DICOM).
- Сделайте отдельную папку под проектом для всех результатов анализа.
2. 3D цифровая модель позвонков
ПРИМЕЧАНИЕ: Все функции подпроцессов взяты из программных средств, собственность которых принадлежит Beijing Intelligent Entropy Science & Technology Co Ltd.
- Вызовите подпроцесс Dicom2Mat на рабочем месте MATLAB, чтобы получить 3D-том из файлов DICOM, хранящихся в папке данных HRCT.
- После выполнения подпроцесса Dicom2Mat просмотрите каждый фрагмент в 3D-томе с помощью графического интерфейса пользователя (GUI), как показано на рисунке 2.
- Затем визуализируйте распределение интенсивности данных HRCT позвонков с помощью функции hist (рис. 3).
- Вызовите подпроцесс NoiseClean , чтобы удалить шум сигнала, образованный устройством по путям к файлу данных HRCT.
- Используйте подпроцесс функции «Позвонки » по тому же пути, чтобы получить модель позвонков, которая также представляет собой 3D-объем, но только со структурой кости (рис. 4). Параметры фильтра верхних частот, интенсивность от 190 до 1,656.
3. 3D-цифровая модель спинномозгового нерва
ПРИМЕЧАНИЕ: Диксон-ин содержит костную структуру, в то время как Диксон-w описывает нейронную структуру.
- Используйте подпроцесс Dicom2Mat в обоих путях последовательностей Dixon-in и Dixon-w и получите их 3D-объем.
- Кроме того, визуализируйте каждый отдельный срез, представляющий собой 3D-объем, используя графический интерфейс, представленный на рисунке 5. Получите доступ к этой визуализации после завершения подпроцесса Dicom2Mat .
- Используйте функцию Spinal_Nerve для реконструкции модели спинномозгового нерва с параметрами фильтра верхних частот интенсивностью от 180 до 643. Поскольку сигналы нерва в последовательности Диксона-w очень высоки, извлеките 3D-объем спинномозгового нерва, отфильтровав точки с низкой интенсивностью.
- Когда подпроцесс Spinal_Nerve будет завершен, проверьте модель, сгенерированную в графическом интерфейсе, показанном на рисунке 6.
4. Регистрация и слияние
ПРИМЕЧАНИЕ: Ключевым моментом является то, что костная архитектура присутствует как в HRCT, так и в последовательности визуализации Dixon-in.
- Скопируйте три полученных на данный момент 3D-тома в путь к файлу проекта, созданного на шаге 3.1. Модели от HRCT и Dixon-in включают в себя одно и то же строение позвонков, а модели от Dixon-in и Dixon-w имеют одинаковые координаты.
- Затем поместите имена файлов трех моделей в подпроцесс vertebra_fusion в качестве входных данных для создания модели слияния. Это можно визуализировать на рисунке 7.
- Слияние обычно хорошо сделано. Если тонкая настройка необходима с точки зрения врача, добавьте параметры координат во всех направлениях к одной и той же функции, чтобы скорректировать модель слияния. Если в слиянии наблюдаются небольшие ошибки с клинической точки зрения, используйте функцию vertebra_fusion для точной настройки координат слияния. Этот процесс включает в себя настройку параметров шести измерений координатного направления (координаты XYZ и их вращение).
- Сделайте отдельную папку в каталоге проекта для вывода результата модели слияния.
5. Файлы цифровых моделей для 3D-печати
ПРИМЕЧАНИЕ: Полностью разработанный аппарат для 3D-печати используется для изготовления вышеупомянутой цифровой модели с реализацией триангуляций Делоне. Здесь использовался 3D-принтер Stratasys J55 Prime.
- Экспортируйте модели слияния, которые будут использоваться для 3D-печати, в последовательностях формата DICOM по пути к файлу каталога слияния. Используйте алгоритм Mat2Dicom для выполнения операции экспорта, введя модель слияния.
- Откройте последовательность файлов DICOM, экспортированную ранее, с помощью Materialise Mimics V20. Чтобы выполнить операцию экспорта, перейдите в меню «Экспорт » на вкладке «Файл» и выберите формат VRML. Путь к файлу для экспорта может быть свободно настроен в соответствии с требованиями пользователя.
- Поскольку прозрачная красочная 3D-печать является профессиональной услугой, сжимайте и упаковывайте файлы VRML и отправляйте их поставщику услуг. Результат 3D-печати показан на рисунке 8.