Summary
Se emplea un sistema de ingeniería inversa para registrar y obtener datos de geometría detallados y completos de placas finales vertebrales. A continuación, se desarrollan modelos paramétricos de placa sorcterácea vertebral, que son beneficiosos para diseñar implantes espinales personalizados, realizar diagnósticos clínicos y desarrollar modelos de elementos finitos precisos.
Abstract
Los datos geométricos detallados y completos de las placas finales de las vértebras son importantes y necesarios para mejorar la fidelidad de los modelos de elementos finitos de la columna vertebral, diseñar y mejorar los implantes espinales, y comprender los cambios degenerativos y la biomecánica. En este protocolo, se emplea un escáner de alta velocidad y alta precisión para convertir los datos morfológicos de las superficies de la placa base en una nube de puntos digital. En el sistema de software, la nube de puntos se procesa y reconstruye en tres dimensiones. A continuación, se realiza un protocolo de medición, que implica un sistema de coordenadas 3D definido para hacer de cada punto una coordenada 3D, tres curvas de superficie sagital y tres frontales que se ajustan simétricamente en la superficie de la placa base, y 11 puntos equidistantes que son seleccionado en cada curva. Finalmente se realizan análisis de medición y espaciales para obtener datos geométricos de las placas finales. Las ecuaciones paramétricas que representan la morfología de curvas y superficies se ajustan en función de los puntos característicos. El protocolo sugerido, que es modular, proporciona un método preciso y reproducible para obtener datos geométricos de placas finales vertebrales y puede ayudar en estudios morfológicos más sofisticados en el futuro. También contribuirá a diseñar implantes espinales personalizados, planificar actos quirúrgicos, realizar diagnósticos clínicos y desarrollar modelos precisos de elementos finitos.
Introduction
Una placa final vertebral es la capa superior o inferior del cuerpo vertebral y sirve como interfaz mecánica para transferir la tensión entre el disco y el cuerpo vertebral1. Consiste en el borde epifisario, que es un labrum óseo fuerte y sólido que rodea el borde exterior del cuerpo vertebral, y la placa central, que es delgada y porosa2.
La columna vertebral está sujeta a una amplia gama de trastornos degenerativos, traumáticos y neoplásicos, que pueden justificar la intervención quirúrgica. Recientemente, los dispositivos espinales como discos artificiales y jaulas han sido ampliamente utilizados. Los parámetros morfométricos precisos y detallados de las placas finales son necesarios para el diseño y la mejora de los implantes espinales con contacto eficaz entre prótesis y vértebras y potencial de crecimiento óseo3. Además, la información sobre la forma exacta y la geometría de las placas finales vertebrales es importante para entender la biomecánica. Aunque el modelado de elementos finitos permite la simulación de las vértebras reales y ha sido ampliamente utilizado para estudiar las respuestas fisiológicas de la columna vertebral a diversas condiciones de carga4,esta técnica es específica del paciente y no generalizable para todos Vértebras. Se ha sugerido que la variabilidad intrínseca de la geometría de las vértebras entre la población general debe tenerse en cuenta al desarrollar el elemento finito modelo5. Por lo tanto, los parámetros geométricos de las placas finales son propicios para la generación de malla y la mejora de la fidelidad en el modelado de elementos finitos.
Aunque la importancia de la coincidencia de la geometría de la placa final y la superficie del implante se ha discutido en estudios anteriores6,7,8, los datos sobre la morfología de las placas finales vertebrales es escasa. La mayoría de los estudios anteriores no han revelado la naturaleza 3D de la placa final9,10,11. Se requiere un análisis espacial para representar mejor y completamente la morfología de la placa final12,13,14. Además, la mayoría de los estudios han empleado técnicas de medición de menor precisión10,15,16. Además, se ha notificado un aumento significativo cuando los parámetros de geometría se miden mediante el empleo de radiografías o tomografías computarizadas (CT)17,18. Aunque la resonancia magnética (RM) se considera no invasiva, es menos precisa a la hora de definir los márgenes precisos de las estructuras óseas11. Debido a la falta de un protocolo de medición estandarizado, existen grandes diferencias entre los datos geométricos existentes.
En los últimos años, la ingeniería inversa, que puede digitalizar las partes físicas existentes en modelos sólidos computarizados, se ha aplicado cada vez más al campo de la medicina. La técnica hace factible desarrollar una representación precisa del carácter anatómico de superficies vértebras sofisticadas. El sistema de ingeniería inversa incluye dos subsistemas: el sistema de instrumentación y el sistema de software. El sistema de instrumentación adoptado en este protocolo tiene un escáner de base plana de rango 3D óptico sin contacto, que es de alta velocidad y muy preciso (precisión 0,02 mm, 1.628 x 1.236 píxeles). El escáner puede capturar de manera eficiente (tiempo de entrada 3 s) información morfológica de la superficie del objeto de destino y convertirla en nube de puntos digital. El sistema de software (es decir, software de ingeniería inversa) es una aplicación informática para el procesamiento de datos de nube sortica (ver Tabla de materiales),reconstrucción del modelo de superficie 3D, edición de curvas y superficies libres y procesamiento de datos (consulte Tabla de Materiales).
Los propósitos del presente informe son (1) diseñar un protocolo de medición y algoritmo para obtener parámetros cuantitativos de placas finales vertebrales basados en una técnica de ingeniería inversa, (2) desarrollar un modelo matemático que permita un modelo realista representación de placas finales vertebrales sin digitalizar demasiados puntos de referencia. Estos métodos serán beneficiosos para la planificación de actos quirúrgicos y el modelado de elementos finitos.
Subscription Required. Please recommend JoVE to your librarian.
Protocol
Este estudio fue aprobado por la junta de ética de la investigación sanitaria del instituto de autores. Como los huesos vertebrales cervicales tienen formas más intrincadas19,el protocolo utiliza las vértebras cervicales como una ilustración para facilitar la investigación relevante.
1. Preparación de materiales, escaneo y procesamiento de imágenes
- Recoger una vértebra cervical seca sin deformación patológica o partes rotas.
- Coloque la vértebra verticalmente en la plataforma del escáner(Figura 1, ver Tabla de materiales), con la placa final mirando hacia la lente de la cámara. Utilice la fuente de luz activa del escáner. A continuación, inicie el proceso de análisis para obtener datos de nube de puntos (. ASC).
NOTA: De acuerdo con las imágenes previas a la exploración, ajuste el escáner y la posición de la vértebra para capturar la mayor cantidad de información de morfología de la superficie como sea posible. - Abra el software especialmente utilizado para procesar nubes de puntos (consulte Tabla de materiales). Haga clic en Importar para importar los datos de nube de puntos y generar el gráfico digital de la vértebra. Establezca la frecuencia de muestreo en 100%, seleccione Mantener datos completos en muestreo, seleccione la unidad de datos como milímetros y haga clic en Puntosde sombra . Utilice la herramienta Selección de lazo para seleccionar puntos redundantes en el gráfico y, a continuación, haga clic en Eliminar para eliminarlos. Haga clic en Reducir ruido y establezca el nivel de suavizado al máximo para reducir el ruido y los picos(Figura 2A,B).
NOTA: Hay instrucciones básicas de funcionamiento del software en la parte inferior de la GUI (interfaz gráfica de usuario). Los puntos de ruido con espolones afilados obvios lateral o verticalmente deben eliminarse para reducir el error. - Haga clic en Ajustar para empaquetar los datos de imagen en un archivo de formato .stl para transformar la nube de puntos en malla, lo que convertirá un objeto de punto en un objeto poligonal.
NOTA: El software de ingeniería inversa suele aceptar el formato 3D de estilo .stl. - Abra el software especialmente utilizado para la reconstrucción 3D y el procesamiento de datos (ver Tabla de Materiales). Haga clic en Archivo y luego en Nuevo en el submenú. Seleccione Pieza en la lista de tipos. Haga clic en Inicioy, a continuación, en Forma en el submenú y, a continuación, en Editor de formas digitalizadas. Haga clic en el icono Importar en la barra de herramientas en el lado derecho de la GUI. En la ventana Importar, seleccione el archivo de formato .stl y, a continuación, haga clic en Aplicar > Aceptar. Haga clic en Ajustar todo en el icono de la barra de herramientas en la parte inferior para cargar la imagen reconstruida en la ventana principal del software de presentación.
NOTA: Los pasos 1.5–2.3.3 se realizan con el mismo software. - Haga clic en Activar en la barra de herramientas del lado derecho. En la ventana Activar, seleccione Modo de reventado > Tipo poligonal > Trampa interior. A continuación, seleccione la placa final vertebral en la imagen 3D para eliminar los componentes vertebrales innecesarios, como los elementos posteriores y los osteofitos(Figura 2C).
2. Cuantificación de la morfología 3D de la placa final
- Definición del sistema de coordenadas 3D de la placa final
- Haga clic en Inicio > Forma en el submenú y, a continuación, en Diseño de forma generativa. Haga clic en el icono Punto de la barra de herramientas situado a la derecha. Marque tres puntos de referencia anatómicos en el borde epifisel: los dos primeros son los puntos finales izquierdo y derecho del borde final de la placa final, respectivamente; el tercero es el punto mediano anterior.
- Haga clic en el icono Línea de la barra de herramientas en el lado derecho y seleccione los dos puntos finales del borde final para definir una línea frontal posterior. Haga clic en el icono Plano, seleccione el tipo de plano que desea curvar y, a continuación, seleccione la línea frontal posterior y el punto mediano anterior para definir el plano de sagital medio.
- Haga clic en Inicio > Forma > Reconstrucción rápida de superficies. Haga clic en el icono Sección plana, escriba 1 en la opción de número y, a continuación, seleccione la imagen de placa final y el plano de sagital medio para generar una curva intersectante. Haga clic en Curva en el icono Escanear y seleccione la intersección de la curva intersectante y el borde epifisal posterior. Defina la intersección como el punto mediano posterior.
- Haga clic en Inicio > Forma > Diseño de forma generativa. Haga clic en el icono Línea y seleccione el punto mediano anterior y el punto mediano posterior para definir un diámetro de sagital medio. Haga clic en el icono Punto y, a continuación, en Points and Planes Repetition en el submenú. A continuación, seleccione el diámetro de sagital medio e introduzca 1 en la opción Instancia(s) para definir el punto medio del diámetro de la sagital media.
- Haga clic en el icono Sistema de eje en la barra de herramientas de la parte inferior. A continuación, seleccione el punto medio del diámetro de la sagital media como origen, la línea paralela a la línea frontal posterior como el eje X, el diámetro medio-sagital como el eje Y, y la línea que apunta hacia adelante y perpendicular al plano x-y como el eje z(Figura 3 ).
NOTA: Los dos puntos finales del borde final se eligen como puntos de referencia porque son coherentes y muestran una variación mínima en presencia de osteofitos10.
- Montaje de curvas y puntos característicos en la superficie de la placa final(Figura 4A–D)
- Haga clic en el icono Punto y, a continuación, en Points and Planes Repetition en el submenú. Seleccione el diámetro de sagital medio e introduzca 3 en la opción Instancia(s) para dividir el diámetro de la sagital media por igual en cuatro partes.
- Haga clic en Inicio > Forma > Reconstrucción rápida de superficies. Haga clic en el icono Sección plana, escriba 1 en la opción Número y, a continuación, seleccione la imagen de placa final y el plano x-z para generar una curva intersectante. Haga clic en Curva en el icono Escanear y seleccione las dos intersecciones del plano x-z y el borde epifiscaso.
- Defina la línea entre las dos intersecciones como el diámetro frontal medio. De la misma manera, divida el diámetro frontal medio por igual en cuatro partes.
NOTA: Cuando la placa final no sea simétrica en relación con el plano med-sagital, elija uno de los dos puntos finales de la curva frontal media que tenga una distancia vertical más corta al plano z-y. A continuación, defina el diámetro frontal medio como 2 veces la longitud del más corto, y divídalo por igual en cuatro partes. - Haga clic en el icono Medir entre en la barra de herramientas en la parte inferior para medir la longitud de un cuarto del diámetro medio de la sagital. Haga clic en el icono Sección plana, introduzca 2 en la opción Número, introduzca el valor medido en la opción Paso y, a continuación, seleccione la imagen de placa final y el plano x-z para generar dos curvas de ajuste en un lado de la pieza frontal. Haga clic en Intercambiar para generar dos curvas de conexión en el otro lado. Del mismo modo, obtenga las otras tres curvas de ajuste en el plano sagital.
NOTA: Las dos curvas de ajuste frontal es superposición con las dos curvas de ajuste de mediados de sagital. - Seleccione 11 puntos equidistantes en cada curva para mediciones posteriores. El método específico es el siguiente:
- Tomando como ejemplo la curva de la sagital media, divida el diámetro de la sagital media por igual en 10 partes, lo que resulta en una suma de 11 puntos, incluidos nueve puntos intermedios y dos puntos finales (consulte los pasos 2.1.3 y 2.2.1).
- Pase por cada punto equidistante, obtenga nueve curvas de conexión en la superficie de la placa base (consulte el paso 2.2.2). Haga clic en Curva en el icono Escanear y seleccione la intersección de las curvas de conexión y la curva de sagital medio. Por último, obtenga un total de 66 puntos en cada placa final (11 puntos por curva multiplicado por seis curvas). Haga clic en el icono Elemento de medida de la barra de herramientas de la parte inferior para medir las coordenadas de cada punto.
- Medición de parámetros morfológicos de placa final
- Parámetro de línea:
- Haga clic en el icono Medir entre para medir la longitud del parámetro de línea que es la distancia entre dos puntos medidos.
- Parámetros de concavidad:
- Cree un plano paralelo al plano x-y(Figura 5A):haga clic en Inicio > Forma > Diseño de forma generativa. Haga clic en el icono Croquis de la barra de herramientas situado a la derecha y, a continuación, haga clic en el plano x-y. Haga clic en el icono Círculo, haga clic en Origen en la superficie de la placa final, arrastre el cursor del ratón a una distancia adecuada y, a continuación, haga clic. Haga clic en el icono Salir de Workbench, luego en el icono Rellenar y, a continuación, haga clic.
- Haga clic en el icono Desplazamiento, seleccione el plano relleno e introduzca un valor adecuado en la opción de desfase hasta que sea tangente a la pieza más cóncava y haga zoom. Haga clic en Inicio > Forma > Reconstrucción rápida de superficies. A continuación, haga clic en el icono de curva 3D para buscar y crear el punto más cóncavo. Haga clic en el icono Medir elemento para medir las coordenadas del punto más cóncavo(figura 5B).
- Haga clic en el icono Medir entre y, a continuación, seleccione el punto más cóncavo y el plano x-y para medir toda la profundidad de la concavidad de la placa final. Del mismo modo, busque y cree la profundidad más cóncava en un plano determinado y mida sus coordenadas.
- Haga clic en el icono Proyección de la barra de herramientas en el lado derecho y, a continuación, seleccione el punto más cóncavo y el plano x-y para obtener el punto proyectivo. Haga clic en el icono Elemento de medida para medir las coordenadas del punto proyectivo y determinar su distribución en función de las coordenadas.
- Parámetros de superficie:
- Haga clic en el icono Medir inercia de la barra de herramientas en la parte inferior y haga clic en la superficie de la placa final para medir su área. Haga clic en el icono Activar y seleccione la placa central a lo largo de los márgenes internos del anillo epifisel (consulte el paso 1.6) y, a continuación, haga clic en el icono Medir inercia para medir su área(Figura 5C). Haga clic en el icono Activar, luego en la placa central y, finalmente, en el icono Intercambiar de la ventana Activar para obtener un borde epifisisco. Luego, mida su área.
- Parámetro de línea:
3. Desarrollo del modelo matemático de superficie de placa final
- Determinar el orden de ajuste de la ecuación paramétrica
- Abra el software de análisis y visualización de datos (consulte Tabla de materiales). Entrada x á [datos correspondientes] en la ventana de comandos. Haga clic en Intro.
NOTA: Los "datos correspondientes" se refieren a los datos de coordenadas x de los 11 puntos característicos de una curva que se ha medido en los pasos anteriores. Haga clic en Intro después de introducir cada comando, con la misma aplicación a las operaciones posteriores. Los pasos 3.1–5.5 se realizan uniformemente con el mismo software. - Del mismo modo, introduzca z á [datos correspondientes].
- Ingrese el código para i-1:5 z2-polyfit(x,z,i); Z-polyval(z2,x); si sum((Z-z).-2)<0.01I break; final; fin.
NOTA: El protocolo establece la suma de errores de cuadrados por debajo de 0,01 para obtener una mayor precisión, cuyo valor se puede reajustar para satisfacer diversas demandas. - Haga clic en Intro para obtener un valor C que sea el orden de ajuste deseado.
- Abra el software de análisis y visualización de datos (consulte Tabla de materiales). Entrada x á [datos correspondientes] en la ventana de comandos. Haga clic en Intro.
- Ajuste de la ecuación de parámetros
- Introduzca cftool y haga clic en Intro para abrir la herramienta de ajuste de curva.
- Introduzca las coordenadas de una curva en la ventana de comandos (consulte los pasos 3.1.1 y 3.1.2). En la herramienta Ajuste de curvas, seleccione datos de coordenadas x al ajustar curvas de plano frontal y datos de coordenadas Y al ajustar curvas de plano sagital en la opción de datos x, seleccione datos de coordenadas z en la opción de datos y, seleccione polinomioe introduzca el orden de ajuste e introduzca el orden de ajuste Obtenido. A continuación, el software mostrará la ecuación paramétrica y la bondad del ajuste automáticamente.
NOTA: Como la curva es una imagen 2D, la opción de trabajo predeterminada son las opciones x e y de la herramienta Ajuste de curva al ajustar una curva. - De la misma manera, introduzca las coordenadas 3D de los 66 puntos y haga coincidir los datos de coordenadas con las opciones de eje correspondientes. Seleccione polinomio e introduzca el orden de ajuste para obtener la ecuación paramétrica de la superficie de la placa final(Figura 6B).
4. Adquisición de datos geométricos basados en la ecuación paramétrica
- Introduzca los valores de coordenadas x e y de cualquier punto de la placa final en la ventana de comandos.
- Entrada PX1, PX2, PX3....
NOTA: Px son los parámetros de la ecuación paramétrica que se han ajustado utilizando polinomio en los pasos anteriores. - Introduzca la ecuación y haga clic en Entrar para obtener el resultado (es decir, formato de entrada: z - P00 + P10*x + P01* y + P20*x 2 + P11*x*y + P02*y 2 + P30*x 3 + P21*x 2 * y + P12 *x*y2 + P03*y3 + P40*x 4 + P31*x 3 *y + P22*x 2 * 2 + 2 + P13*x * y 3 + P04* y 4).
5. Representación de la placa final basada en la ecuación paramétrica
- Entrada PX1, PX2, PX3.... en la ventana de comandos.
- Introduzca el código X-N1:0.01:N2;.
NOTA: N1–N2 es el rango de datos del eje X (es decir, los valores de los dos puntos finales de la curva de la mitad-coronal). - Ingrese el código "Y-N3:0.01:N4;".
- Ingrese la ecuación (es decir, z -(x,y)P00 + P10.*x + P01.*y + P20.*x.-2 + P11.*x.*y + P02.*y.-2 + P30.*x.-3 + P21.*x.-2.*y+ P12.*x.*y.2 + P 03.*y.-3 + P40.*x.-4 + P31.*x.-3.*y + P22.*x.-2.*y.-2+ P13.*x.*y.-3 + P04.*y.-4;).
- Introduzca el código ezmesh(z, [N1,N2,N3,N4]) para obtener gráficos de simulación 3D(Figura 6C).
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
Utilizando el escáner de base plana de rango 3D óptico de alta precisión, las placas finales se convirtieron en más de 45.000 puntos digitales, que caracterizan adecuadamente la morfología(Figura 2A,B).
En el protocolo de medición, se llevó a cabo el análisis espacial de las superficies de la placa final. Se instalaron y cuantificaron curvas representativas en la superficie para caracterizar la morfología(Figura 4B). Los parámetros lineales se midieron calculando la distancia entre dos puntos finales. Las mediciones obtenidas incluyen la profundidad de la concavidad y la ubicación del ápice de la concavidad en el plano mid-sagital, además de las de toda la concavidad de la placa final y cualquier sección específica(Figura 5B). Los componentes de las placas finales, el borde epifisel y la placa central se separaron(Figura 5C),y sus longitudes y áreas se obtuvieron convenientemente.
Se digitalizaron y analizaron un total de 138 placas finales vertebrales cervicales, y se estableció el modelo matemático de la placa final. El protocolo establece las sumas de error al cuadrado por debajo de 0.01, y se concluyó que el uso de la función polinómica de cuatro óses podría lograr satisfacción.
La ecuación paramétrica de cada curva se dedujo en función de las coordenadas de 11 puntos: f(x) - P1*x 4 + P2*x-3 + P3*x-2 + P4*x + P5. P1, P2, P3, P4 y P5 fueron los parámetros, cuyos valores exactos se muestran en el Cuadro 1.
La ecuación paramétrica que representa las características morfológicas de la superficie de la placa final es:
F(x, y)á P00 + P10*x + P01*y + P20*x 2 + P11*x*y + P02* y 2 + P30*x 3+ P21*x 2 *y + P12*x *y 2 +P *x 2* y2 + P13*x*y-3 + P04* y-4
Dónde: PXYs son los parámetros, que se dedujeron de las coordenadas premedidas de 66 puntos (Tabla 2).
Figura 1: El escáner de base plana de rango 3D óptico sin contacto. El escáner, que se basa en la tecnología de medición óptica 3D de desplazamiento de fase multifrecuencia heterodina, incluye medición óptica (integración alrededor de dos cámaras y un proyector) y dispositivos de control. La precisión de este instrumento es de 0,02 mm, y los píxeles son 1628 x 1236. El escáner puede digitalizar eficientemente (tiempo de entrada 3 s) la geometría de superficie de un objeto de destino. Haga clic aquí para ver una versión más grande de esta figura.
Figura 2: La nube puntual de la superficie vertebral y la reconstrucción 3D de la placa final. (A) y (B) son las superficies inferiores y superiores de una vértebra cervical generada por el software especialmente utilizado para el procesamiento de nubes de puntos, respectivamente. (C) y (D) son la reconstrucción 3D de las placas finales inferiores y superiores generadas por el software especialmente utilizado para la reconstrucción 3D y el procesamiento de datos, respectivamente. Los elementos posteriores y osteofitos se eliminan de las vértebras, dejando sólo la placa base. El plano de mejor ajuste se define a través de los puntos anteriores y posteriores de los procesos de uncinado bilateral, y las dos curvas formadas por el plano y la placa final de mejor ajuste son los límites de la junta uncovertebral y la placa final caudal. Haga clic aquí para ver una versión más grande de esta figura.
Figura 3: Definición del sistema de coordenadas 3D de la placa final. Marcado de tres puntos de referencia anatómicos en el borde epifisel: los dos primeros son los puntos finales izquierdo y derecho del borde final de la placa final, respectivamente; el tercero es el punto mediano anterior. La línea frontal posterior está formada por los dos extremos del borde final, que definen el plano medio-sagital con el punto mediano anterior. El punto mediano posterior está determinado por el plano mid-sagital y el borde epifisele posterior, que forman el diámetro medio-sagital con el punto mediano anterior. El origen es el punto medio del diámetro medio-sagital. El eje Y está determinado por el diámetro medio-sagital y apuntando hacia adelante. El eje X es la línea paralela a la línea frontal posterior. El eje z es normal al plano x-y. Haga clic aquí para ver una versión más grande de esta figura.
Figura 4: Los pasos de las curvas y puntos característicos de ajuste en la superficie de la placa final. (A) Divida el diámetro medio-sagital y el diámetro medio-frontal por igual en cuatro partes. (B) Pasar por cada punto equidistante, y elegir seis curvas de superficie simétricamente, tres de las cuales son las curvas de intersección del plano frontal y la superficie de la placa final, y las otras tres en el plano sagital. (C) Divida el diámetro de la sagital media por igual en 10 partes. (D) Pasando por cada punto equidistante, los planos frontales y la curva media-sagital forman nueve intersecciones, dando como resultado una suma de 11 puntos, junto con los dos puntos finales. Haga clic aquí para ver una versión más grande de esta figura.
Figura 5: Medición de la profundidad de la concavidad de la placa final y el área de la superficie. (A) Cree un plano paralelo al plano x-y. (B) Desfase el plano hasta que sea tangente al punto más cóncavo, y la profundidad de la concavidad de la placa final es la distancia perpendicular entre el punto más cóncavo y el plano x-y. (C) Dibuje una línea a lo largo de los márgenes internos del anillo epifisaparal para dividir la placa final en la placa central y el borde epifisal. Haga clic aquí para ver una versión más grande de esta figura.
Figura 6: La reconstrucción 3D y las representaciones de una placa final inferior. (A) La reconstrucción 3D de la superficie inferior de la placa final generada por el software especialmente utilizado para la reconstrucción 3D y el procesamiento de datos. (B) y (C) son las representaciones de la placa final inferior generada por el software de análisis y visualización de datos. Haga clic aquí para ver una versión más grande de esta figura.
| Nivel de placa final | Curva | Parámetros | ||||
| P1 | P2 | P3 | P4 | P5 | ||
| C6 superior | Fac | 0 | 0 | -0.0128 | -0.0028 | 0.02523 |
| Fmc | 0 | 0 | -0.0199 | 0.00074 | 0.3693 | |
| Fpc | 0 | 0 | -0.0329 | 0.00739 | 0.5323 | |
| Slc | 0 | 0.00176 | -0.0113 | -0.0419 | -0.0419 | |
| Smc | 0.00011 | 0.00232 | -0.016 | -0.0986 | 0.4712 | |
| Fuente | 0 | 0.00179 | -0.0096 | 0.04451 | -0.0394 | |
| C6 inferior | Fac | 0 | -0.0001 | -0.0225 | 0.00594 | 1.223 |
| Fmc | 0 | 0 | -0.016 | -0.0082 | 1.729 | |
| Fpc | 0 | 0 | -0.0033 | -0.0033 | 1.404 | |
| Slc | 0.00012 | 0.00087 | -0.0347 | -0.0962 | 1.448 | |
| Smc | 0.00025 | 0.00064 | -0.0495 | -0.0331 | 1.846 | |
| Fuente | 0 | 0.00079 | -0.0295 | -0.0828 | 1.362 | |
Tabla 1: Los parámetros de la ecuación para representar la curva de la superficie de la placa final. Sólo se enumeran los datos de la sexta placa final vertebral cervical. Px : los parámetros de la ecuación. En cada placa final, se eligieron simétricamente seis curvas de superficie; tres de ellos estaban en el plano frontal y se denominaron la curva anterior (FAC), la curva media (FMC) y la curva posterior (FPC); los otros tres en el plano sagital se denominaron la curva izquierda (SLC), la curva media (SMC) y la curva derecha (SRC). Los parámetros con un valor absoluto inferior a 0.0001 se representan como 0 aquí.
| Parámetros | C3 inf | C4 sup | C4 inf | C5 sup | C5 inf | C6 sup | C6 inf | C7 sup |
| p00 | 1.989 | 0.4187 | 2.004 | 0.3383 | 1.913 | 0.4276 | 1.779 | 0.5674 |
| p10 | -0.0022 | -0.0043 | 0.00542 | -0.0208 | -0.0111 | 0.0012 | -0.0043 | -0.0052 |
| p01 | -0.0356 | -0.0868 | -0.0537 | -0.0826 | -0.0257 | -0.098 | -0.0407 | -0.0642 |
| p20 | 0.01286 | -0.0252 | -0.0146 | -0.0299 | -0.0253 | -0.0264 | -0.0175 | -0.0088 |
| p11 | 0.00092 | 0.00071 | -0.0009 | 0.00018 | -0.0002 | -0.0012 | 0.00117 | 0.00021 |
| P02 | -0.0529 | -0.0151 | -0.0525 | -0.012 | -0.0418 | -0.0142 | -0.0396 | -0.0134 |
| p30 | 0 | -0.0001 | 0.00013 | 0.00024 | 0.00017 | 0 | 0 | 0 |
| p21 | -0.0011 | 0.00299 | -0.0012 | 0.00363 | -0.0021 | 0.00306 | -0.0019 | 0.00194 |
| p12 | 0 | 0.00048 | -0.0004 | 0.00033 | 0.00014 | 0 | -0.0001 | 0 |
| p03 | 0.00062 | 0.00204 | 0.00089 | 0.00206 | 0.00046 | 0.00208 | 0.00077 | 0.00115 |
| p40 | 0.0002 | 0 | 0.0002 | 0 | 0.00024 | 0 | 0 | 0 |
| p31 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| p22 | 0.00017 | 0.00013 | 0 | 0.00015 | 0.00015 | 0.00017 | 0.00032 | 0 |
| p13 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| p04 | 0.00023 | 0.00013 | 0.00024 | 0 | 0 | 0 | 0 | 0 |
Tabla 2: Los parámetros de la ecuación paramétrica que representa la morfología de la superficie de la placa final. Px - los parámetros de la ecuación; inf - placa inferior; sup - placa superior. Los parámetros con un valor absoluto inferior a 0.0001 se representan como 0 aquí. Esta tabla se ha modificado de una publicación anterior3.
| Medidas | Fiabilidad intratest | Medidas | RE vs Caliper | ||
| Apd | Primera medida | 15,76 x 1,3 | Apd | RE | 16,47 x 1,31 |
| Medición | 15,86 x 1,61 | Calibrador | 16,26 x 1,27 | ||
| Icc | 0.85 | Cronbach alfa | 0.99 | ||
| Cmd | Primera medida | 19,71 x 2,47 | Cmd | RE | 20,7 x 3,05 |
| Medición | 19,41 x 2,43 | Calibrador | 20,45 x 3,21 | ||
| Icc | 0.96 | Cronbach alfa | 0.99 | ||
Tabla 3: Fiabilidad de las mediciones. Los datos fueron de media: desviación estándar (mm). ICC - coeficiente de correlación intra-clase; APD - diámetro antero-posterior; CMD - diámetro mediolateral central; RE: el sistema de ingeniería inversa. Esta tabla se ha modificado de una publicación anterior. 3
| Valor de las medidas | N | Valor de coordenadaz Z | T | P | R |
| Puntos originales | 15 | 1,75-0,87 | 0.26 | 0.8 | 0.98 |
| Puntos de comparación | 15 | 1,74 x 0,91 |
Tabla 4: La validez del modelo geométrico que representa la morfología de la placa final. Los datos se representan como media de desviación estándar (mm). Los puntos originales son 15 puntos seleccionados al azar en la imagen de reconstrucción 3D original. Puntos de comparación: puntos correspondientes generados automáticamente a partir de ecuaciones paramétricas; R - coeficiente de correlación.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
La ingeniería inversa se ha aplicado cada vez con éxito al campo de la medicina, como la craneoplastia20,oral21e implantes maxilofaciales21. Las mediciones de ingeniería inversa, a saber, la digitalización de la superficie del producto, se refieren a la conversión de información de superficie en datos de nube de puntos que emplean equipos y métodos de medición específicos. Sobre la base de estos datos, se pueden realizar modelos de superficies complejas, evaluación, mejoras y fabricación. La medición digital y el procesamiento de datos son una tecnología básica y clave utilizada en la ingeniería inversa.
En este protocolo, la información morfológica precisa y detallada de las placas finales vertebrales se registra utilizando un sistema de escaneo óptico 3D sin contacto, que se basa en la tecnología de medición óptica 3D, de cambio de fase y de cambio de fase de heterodina. El escáner está hecho principalmente de dispositivos de control y una medición óptica que integra dos cámaras y un proyector. En comparación con otros instrumentos de medición, el escáner es altamente preciso y eficiente y evita el escaneo punto por punto. Al capturar datos de nube de puntos, el cabezal de escaneado normalmente no está en contacto con el objeto, de forma que no haya efectos de deformación. La fiabilidad, validez y precisión del escáner para la grabación de la morfología de la superficie han sido bien establecidas2,3,22. Se ha verificado la replicabilidad de estas mediciones.
Para verificar la precisión de las mediciones tomadas por el sistema de ingeniería inversa, se midieron 20 placas finales utilizando una pinza digital y se evaluaron utilizando Cronbach alpha. Para la fiabilidad dentro de la prueba, se seleccionaron aleatoriamente 16 placas finales de las 138 placas finales vertebrales y se midieron dos veces a intervalos de 2 semanas, luego se evaluaron utilizando un coeficiente de correlación intra-clase. Los resultados mostraron un gran acuerdo y fiabilidad (Tabla 3). El software de ingeniería inversa implica potentes mediciones, procesamiento de datos, detección de errores y funciones de edición de curvas y superficies libres. También puede construir y ajustar de forma inteligente y eficiente curvas y superficies, y la reconstrucción del modelo de superficie 3D contribuye a mediciones precisas23.
Existen aplicaciones importantes y considerables para datos detallados y completos de la anatomía de las vértebras, como el diseño de implantes espinales, la mejora de la fidelidad de los modelos de elementos finitos de la columna vertebral y el desarrollo de modelos matemáticos. La placa final vertebral es esencial para mantener la integridad y la función del disco intervertebral, y también sirve como una interfaz mecánica para transferir tensión. Por lo tanto, la cuantificación de la geometría de la placa final es importante. Con la ayuda de la ingeniería inversa, la morfología de la placa final se puede cuantificar de forma inteligente y completa. En este protocolo, se instalan seis curvas características en la superficie de cada placa final, y se establece un sistema de coordenadas 3D para cuantificar la morfología espacial.
Además, se desarrolla un modelo paramétrico de la placa final para instituir evaluaciones cuantitativas precisas y reproducibles y para desarrollar modelos de elementos finitos biomecánicos personalizados. El modelo paramétrico de la superficie de las placas finales puede producir representaciones rápidas, realistas y precisas que pueden ser visualizadas y analizadas convenientemente por los investigadores.
La inclusión de más puntos de referencia mejorará la precisión, pero es lenta y costosa. En este protocolo, se propone que 66 puntos de seis curvas superficiales sean adecuados para describir las características morfológicas. Las pruebas de fiabilidad también se llevan a cabo comparando los valores de coordenadas de 15 puntos seleccionados aleatoriamente con los valores correspondientes que se generan automáticamente a partir de ecuaciones paramétricas. El resultado revela que el modelo paramétrico tiene una buena fiabilidad y reproducibilidad puede servir como una representación realista de la superficie de la placa final (Tabla 4). Cabe señalar que el modelo paramétrico se puede derivar sobre la base de otras modalidades de diagnóstico por imágenes, como la TC y la RMN.
Dado que los escáneres sin contacto son susceptibles a la luz ambiental, es fundamental mantener la luz ambiental estable y se recomiendan fuentes de luz activas. Si hay grasa residual en la superficie de la placa base, el polvo de talco infantil debe ser manipulado suavemente para evitar el riesgo de verse afectado por las características de reflectancia espacial de la superficie del objeto. Las vértebras cervicales subaxiales tienen un componente especial: la articulación uncovertebral. Para distinguirla de la placa final, se define un plano de mejor ajuste mediante el método de mínimos cuadrados. A continuación, la curva de intersección formada por el plano de mejor ajuste, y la superficie de la placa final es el límite entre la junta uncovertebral y la placa final superior(Figura 2D).
La operación específica es la siguiente: haga clic en Inicio > Forma > Diseño de forma generativa. Haga clic en el icono Punto de la barra de herramientas en el lado derecho y, a continuación, seleccione los puntos anteriores y posteriores de los procesos de uncinar bilaterales en la imagen 3D. Haga clic en el icono Plano y seleccione Media a través de puntos en el tipo de plano para definir el plano que mejor se ajuste. Haga clic en Inicio > Forma > Reconstrucción rápida de superficies. Haga clic en el icono Sección plana y, a continuación, seleccione la imagen 3D y el plano de ajuste recomendado.
El marcado preciso de los tres puntos anatómicos en la superficie de la placa final al establecer el sistema de coordenadas 3D es fundamental. El software de ingeniería inversa permite un cambio flexible de la imagen de reconstrucción y mejora el contraste que ayuda a identificar los puntos de referencia. Alternativamente, es importante evaluar la idoneidad del sistema de coordenadas en función de si la línea intersectante de los planos semisagitato y coronal definidos es perpendicular a la sección de la placa final y, a continuación, ajustar el sistema en consecuencia. También se evaluaron las pruebas intraobservadores y el resultado indicaba una buena fiabilidad (cuadro 3).
Este protocolo requiere múltiples habilidades y técnicas, incluida la adquisición y el procesamiento de datos de nube de puntos, la reconstrucción y el análisis de imágenes y el desarrollo de modelos paramétricos. Para un principiante, puede tomar tiempo completar todo el proceso. Sin embargo, como sólo se utilizan unos pocos módulos del software en este protocolo y el procedimiento es modular, requiere una curva de aprendizaje corta para tener mucha experiencia.
En conclusión, el protocolo descrito proporciona un método preciso y reproducible para obtener datos de geometría detallados y completos de las placas finales vertebrales. Un modelo paramétrico también se desarrolla sin digitalizar demasiados puntos de referencia, lo que es beneficioso para diseñar implantes espinales personalizados, planificar actos quirúrgicos, hacer diagnósticos clínicos y desarrollar modelos precisos de elementos finitos.
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
Los autores no declaran intereses financieros en competencia.
Acknowledgments
Este trabajo fue financiado por Key Discipline Construction Project de Pudong Health Bureau of Shanghai (PWZxk2017-08) y la National Natural Science Foundation of China (81672199). Los autores quieren agradecer a Wang Lei por su ayuda en la corrección de una versión anterior y a Li Zhaoyang por su ayuda en el desarrollo del modelo paramétrico.
Materials
| Name | Company | Catalog Number | Comments |
| Catia | Dassault Systemes, Paris, France | https://www.3ds.com/products-services/catia/ | 3D surface model reconstruction, free curve and surface editing and data processing |
| Geomagic Studio | Geomagic Inc., Morrisville, NC | https://cn.3dsystems.com/software?utm_source=geomagic.com&utm_medium=301 | point cloud data processing |
| MATLAB | The MathWorks Inc., Natick,USA | https://www.mathworks.com/ | analyze data, develop algorithms, and create models |
| Optical 3D range flatbed scanner | Xi’an XinTuo 3D Optical Measurement Technology Co.Ltd., Xi’an, Shaanxi, China | http://www.xtop3d.com/ | acquire surface geometric parameters and convert into digital points |
References
- Wang, Y., Battie, M. C., Boyd, S. K., Videman, T. The osseous endplates in lumbar vertebrae: Thickness, bone mineral density and their associations with age and disk degeneration. Bone. 48, 804-809 (2011).
- Wang, Y., Battie, M. C., Videman, T. A morphological study of lumbar vertebral endplates: radiographic, visual and digital measurements. European Spine Journal. 21, 2316-2323 (2012).
- Feng, H., et al. Morphometry evaluations of cervical osseous endplates based on three dimensional reconstructions. International Orthopaedics. , (2018).
- Liebschner, M. A., Kopperdahl, D. L., Rosenberg, W. S., Keaveny, T. M. Finite element modeling of the human thoracolumbar spine. Spine (Phila Pa 1976). 28, 559-565 (2003).
- Niemeyer, F., Wilke, H. J., Schmidt, H. Geometry strongly influences the response of numerical models of the lumbar spine--a probabilistic finite element analysis. Journal of Biomechanics. 45, 1414-1423 (2012).
- Lin, C. Y., Kang, H., Rouleau, J. P., Hollister, S. J., Marca, F. L. Stress analysis of the interface between cervical vertebrae end plates and the Bryan, Prestige LP, and ProDisc-C cervical disc prostheses: an in vivo image-based finite element study. Spine (Phila Pa 1976). 34, 1554-1560 (2009).
- Cao, J. M., et al. Clinical and radiological outcomes of modified techniques in Bryan cervical disc arthroplasty. Journal of Clinical Neuroscience. 18, 1308-1312 (2011).
- de Beer, N., Scheffer, C. Reducing subsidence risk by using rapid manufactured patient-specific intervertebral disc implants. The Spine Journal. 12, 1060-1066 (2012).
- Chen, H., Zhong, J., Tan, J., Wu, D., Jiang, D. Sagittal geometry of the middle and lower cervical endplates. European Spine Journal. 22, 1570-1575 (2013).
- Tan, S. H., Teo, E. C., Chua, H. C. Quantitative three-dimensional anatomy of cervical, thoracic and lumbar vertebrae of Chinese Singaporeans. European Spine Journal. 13, 137-146 (2004).
- Zhou, S. H., McCarthy, I. D., McGregor, A. H., Coombs, R. R., Hughes, S. P. Geometrical dimensions of the lower lumbar vertebrae--analysis of data from digitised CT images. European Spine Journal. 9, 242-248 (2000).
- Cukovic, S., Devedzic, G., Ivanovic, L., Lukovic, T. Z., Subburaj, K. Development of 3D Kinematic Model of the Spine for Idiopathic Scoliosis Simulation. Computer-Aided Design and Applications. 7, 153-161 (2010).
- Cukovic, S., Devedzic, G. 3D modeling and simulation of scoliosis: An integrated knowledgeware approach. , 411-415 (2015).
- Ćuković, S., et al. Non-Ionizing Three-Dimensional Estimation of Axial Vertebral Rotations in Adolescents Suffering from Idiopathic Scoliosis. , (2018).
- Panjabi, M. M., Duranceau, J., Goel, V., Oxland, T., Takata, K. Cervical human vertebrae. Quantitative three-dimensional anatomy of the middle and lower regions. Spine (Phila Pa 1976). 16, 861-869 (1991).
- Panjabi, M. M., et al. Thoracic human vertebrae. Quantitative three-dimensional anatomy. Spine (Phila Pa 1976). 16, 888-901 (1991).
- Ravi, B., Rampersaud, R. Clinical magnification error in lateral spinal digital radiographs. Spine (Phila Pa 1976). 33, E311-E316 (2008).
- Silva, M. J., Wang, C., Keaveny, T. M., Hayes, W. C. Direct and computed tomography thickness measurements of the human, lumbar vertebral shell and endplate. Bone. 15, 409-414 (1994).
- Langrana, N. A., Kale, S. P., Edwards, W. T., Lee, C. K., Kopacz, K. J. Measurement and analyses of the effects of adjacent end plate curvatures on vertebral stresses. The Spine Journal. 6, 267-278 (2006).
- Chrzan, R., et al. Cranioplasty prosthesis manufacturing based on reverse engineering technology. Medical Science Monitor. 18, (2012).
- De Santis, R., et al. Reverse engineering of mandible and prosthetic framework: Effect of titanium implants in conjunction with titanium milled full arch bridge prostheses on the biomechanics of the mandible. Journal of Biomechanics. 47, 3825-3829 (2014).
- Keating, A. P., Knox, J., Bibb, R., Zhurov, A. I. A comparison of plaster, digital and reconstructed study model accuracy. Journal of Orthodontics. 35, 191-201 (2008).
- Numajiri, T., et al. Designing CAD/CAM Surgical Guides for Maxillary Reconstruction Using an In-house Approach. Journal of Visualized Experiments. , (2018).
