Todos los datos provienen del paciente clínico, cuya operación de SDR se llevó a cabo en el Hospital BJ Dongzhimen. El protocolo sigue las directrices y fue aprobado por el comité de ética de investigación del Hospital Dongzhimen.
NOTA: El mapa completo del protocolo de reconstrucción del modelo se muestra en la Figura 1. Los datos de tomografía computarizada de alta resolución (HRCT) y los datos de Dixon son materias primas para el modelado; luego, la creación del modelo 3D consiste en el registro y la fusión de imágenes. El modelo digital 3D final se imprime con tecnología PolyJet, que es un proceso de impresión 3D de alta precisión que produce piezas suaves y precisas utilizando una amplia gama de materiales. Para describir exactamente la relación espacial entre la vértebra y el nervio espinal, se utilizan datos de HRCT y series de imágenes de Dixon. El escaneo de Dixon puede identificar imágenes de separación de agua y grasa, en las que la serie de imágenes de fase de agua de Dixon se puede usar para extraer la estructura de los nervios espinales, y la serie de imágenes de fase de Dixon-in se puede usar para verificar el registro de la estructura ósea.

Figura 1: El mapa completo del protocolo. La metodología de investigación de este estudio implica la fusión de secuencias de Dixon de TC y resonancia magnética. Específicamente, la estructura de las vértebras de TC se registra con la estructura de vértebras idéntica contenida en la secuencia de Dixon-in, seguida de la fusión con la secuencia de Dixon-w para el nervio espinal. Haga clic aquí para ver una versión más grande de esta figura.
1. Recopilación y preparación de datos
- TC de alta resolución para vértebras
NOTA: La diferencia de parámetros no es sensible al método de investigación.- Establezca los recursos de datos de la estación de la máquina de TC.
NOTA: Aquí se utiliza la máquina CT SIEMENS-CTAWP73396.
- Abra el software Syngo CT 2012B para recibir datos del protocolo de escaneo SpineRoutine_1. Seleccione el tamaño de píxel y el grosor de corte (ST) del conjunto de datos para adaptarlo al tamaño de las vértebras que se pretende representar en el modelo digital 3D.
- Utilice un ST de 1 mm con un tamaño de matriz de 512 píxeles x 512 píxeles, en el que el espaciado de píxeles sea de 0,3320 mm. El tamaño real del volumen 3D alcanzado es de 512 x 512 x 204 vóxeles.
- Secuencia de Dixon para el nervio espinal
NOTA: En este estudio se utiliza una máquina de resonancia magnética de 1,5 T.- Establezca la resolución de imagen de Dixon como 290 píxeles x 320 píxeles, el espaciado de píxeles como 0,9375 mm y el grosor de corte como 3 mm para obtener datos precisos.
- Establezca el tiempo de repetición en 5.160 ms y el tiempo de eco en 94 ms.
- Asegúrese de que cada capa escaneada conste de imágenes de cuatro fases, que son Dixon-in, Dixon-opp, Dixon-F y Dixon-w.
- Preparar archivos de almacenamiento de datos para la reconstrucción del modelo.
NOTA: Una estructura de almacenamiento de datos bien definida es más conveniente para el trabajo de seguimiento.- Haga una carpeta de proyecto que contenga todos los datos pertenecientes al paciente.
- Prepare diferentes rutas de archivo para los datos de HRCT y MRI-Dixon creando diferentes carpetas para los datos de imágenes digitales y comunicaciones en medicina (DICOM).
- Cree una carpeta separada debajo del proyecto para todos los resultados del análisis.
2. El modelo digital 3D de vértebras
NOTA: Todas las funciones de subproceso provienen de herramientas de software, cuya propiedad pertenece a Beijing Intelligent Entropy Science & Technology Co Ltd.
- Llame al subproceso Dicom2Mat en el lugar de trabajo de MATLAB para obtener el volumen 3D de los archivos DICOM almacenados en la carpeta de datos HRCT.
- Después de someterse al subproceso Dicom2Mat , vea cada segmento dentro del volumen 3D a través de la interfaz gráfica de usuario (GUI), como se muestra en la Figura 2.
- Luego, visualice la distribución de intensidad de los datos de la TCM de las vértebras por la función hist (Figura 3).
- Llame al subproceso NoiseClean para eliminar el ruido de señal formado por el dispositivo en las rutas del archivo de datos HRCT.
- Utilice el subproceso de la función Vértebras bajo la misma ruta para obtener el modelo de vértebras , que también es un volumen 3D pero solo con la estructura ósea (Figura 4). Los parámetros del filtro de paso alto, la intensidad oscila entre 190 y 1.656.
3. El modelo digital 3D del nervio espinal
NOTA: Dixon-in contiene estructura ósea, mientras que Dixon-w describe estructura neural.
- Utilice el subproceso Dicom2Mat en ambas rutas de las secuencias Dixon-in y Dixon-w y obtenga su volumen 3D.
- Además, visualice cada sector individual que constituye un volumen 3D utilizando la GUI presentada en la Figura 5. Acceda a esta visualización una vez que se haya completado el subproceso Dicom2Mat .
- Utilice la función Spinal_Nerve para reconstruir el modelo del nervio espinal con parámetros de filtro de paso alto, la intensidad oscila entre 180 y 643. Debido a que las señales del nervio en la secuencia de Dixon-w son muy altas, extraiga el volumen 3D del nervio espinal filtrando los puntos con baja intensidad.
- Cuando finalice el subproceso Spinal_Nerve , compruebe el modelo generado en la GUI que se muestra en la figura 6.
4. Registro y fusión
NOTA: La idea clave es que la arquitectura ósea está presente tanto en la secuencia de imágenes HRCT como en la secuencia de imágenes Dixon-in.
- Copie los tres volúmenes 3D obtenidos hasta ahora en la ruta del archivo del proyecto realizado en el paso 3.1. Los modelos de HRCT y Dixon-in incluyen la misma estructura de vértebras, y los modelos de Dixon-in y Dixon-w tienen las mismas coordenadas.
- A continuación, coloque los nombres de archivo de los tres modelos en el subproceso vertebra_fusion como entrada para generar el modelo de fusión. Esto se puede visualizar en la Figura 7.
- La fusión suele estar bien hecha. Si es necesario realizar un ajuste fino desde la perspectiva del médico, agregue parámetros de coordenadas en todas las direcciones a la misma función para corregir el modelo de fusión. Si se observan errores leves en la fusión desde una perspectiva clínica, utilice la función vertebra_fusion para ajustar las coordenadas de fusión. Este proceso implica ajustes de parámetros a las seis dimensiones de la dirección de coordenadas (coordenadas XYZ y su rotación).
- Cree una carpeta independiente en el directorio del proyecto para generar el resultado del modelo de fusión.
5. Archivos de modelo digital para impresión 3D
NOTA: Para la fabricación del modelo digital mencionado se utiliza un aparato de impresión 3D completamente desarrollado, con la implementación de triangulaciones de Delaunay. Aquí se utilizó la impresora 3D Stratasys J55 Prime.
- Exporte los modelos de fusión que se utilizarán para la impresión 3D en las secuencias de formato DICOM en la ruta de archivo del directorio de fusión. Utilice el algoritmo Mat2Dicom para ejecutar la operación de exportación introduciendo el modelo de fusión.
- Abra la secuencia de archivos DICOM exportada anteriormente con Materialise Mimics V20. Para realizar la operación de exportación, vaya al menú Exportar en la pestaña Archivo y seleccione el formato VRML. La ruta del archivo para la exportación se puede personalizar libremente de acuerdo con los requisitos del usuario.
- Como la impresión 3D colorida transparente es un servicio profesional, comprima y empaque los archivos VRML y envíelos al proveedor de servicios. El resultado de la impresión 3D se muestra en la Figura 8.