Summary
Un système d'ingénierie inverse est utilisé pour enregistrer et obtenir des données géométriques détaillées et complètes des plaques vertébrales. Des modèles paramétriques de la plaque vertébrale sont ensuite développés, qui sont bénéfiques à la conception d'implants rachidiens personnalisés, à la réalisation de diagnostics cliniques et au développement de modèles précis d'éléments finis.
Abstract
Des données géométriques détaillées et complètes des plaques d'extrémité des vertèbres sont importantes et nécessaires pour améliorer la fidélité des modèles d'éléments finis de la colonne vertébrale, concevoir et améliorer les implants rachidiens, et comprendre les changements dégénératifs et la biomécanique. Dans ce protocole, un scanner haute vitesse et très précis est utilisé pour convertir les données de morphologie des surfaces d'extrémité en un nuage de points numériques. Dans le système logiciel, le nuage de points est traité et reconstruit en trois dimensions. Ensuite, un protocole de mesure est effectué, impliquant un système de coordonnées 3D défini pour faire de chaque point une coordonnées 3D, trois courbes de surface sagittale et trois courbes frontales qui sont symétriquement montés sur la surface de la plaque d'extrémité, et 11 points équidistants qui sont sélectionnés dans chaque courbe. Des mesures et des analyses spatiales sont enfin effectuées pour obtenir des données géométriques des plaques d'extrémité. Les équations paramétriques représentant la morphologie des courbes et des surfaces sont ajustées en fonction des points caractéristiques. Le protocole proposé, qui est modulaire, fournit une méthode précise et reproductible pour obtenir des données géométriques des plaques vertébrales et peut aider à des études morphologiques plus sophistiquées à l'avenir. Il contribuera également à la conception d'implants rachidiens personnalisés, à la planification d'actes chirurgicaux, à la réalisation de diagnostics cliniques et à l'élaboration de modèles précis d'éléments finis.
Introduction
Une plaque vertébrale est la coquille supérieure ou inférieure du corps vertébral et sert d'interface mécanique pour transférer le stress entre le disque et le corps vertébral1. Il se compose de la jante épiphysher, qui est un labrum osseux fort et solide entourant le bord externe du corps vertébral, et la plaque d'extrémité centrale, qui est mince et poreuse2.
La colonne vertébrale est sujette à un large éventail de troubles dégénératifs, traumatiques et néoplastiques, qui peuvent justifier une intervention chirurgicale. Récemment, les dispositifs de colonne vertébrale tels que les disques artificiels et les cages ont été largement utilisés. Des paramètres morphométriques précis et détaillés des plaques d'extrémité sont nécessaires pour la conception et l'amélioration des implants spinaux avec le contact efficace de prothèse-vertèbre et le potentiel de croissance d'os3. En outre, l'information sur la forme et la géométrie exactes des plaques vertébrales est importante pour comprendre la biomécanique. Bien que la modélisation des éléments finis permet la simulation des vertèbres réelles et a été largement utilisé pour étudier les réponses physiologiques de la colonne vertébrale à diverses conditions de chargement4, cette technique est spécifique au patient et non généralisable à tous les Vertèbres. Il a été suggéré que la variabilité intrinsèque de la géométrie des vertèbres au sein de la population générale devrait être prise en considération lors du développement du modèle d'élément fini5. Par conséquent, les paramètres géométriques des plaques d'extrémité sont propices à la génération de mailles et à l'amélioration de la fidélité dans la modélisation des éléments finis.
Bien que l'importance de l'appariement de la géométrie de la plaque d'extrémité et de la surface de l'implant ait été discutée dans des études antérieures6,7,8, les données sur la morphologie des plaques vertébrales sont rares. La plupart des études précédentes n'ont pas révélé la nature 3D de la plaque d'extrémité9,10,11. Une analyse spatiale est nécessaire pour mieux et bien représenter la morphologie de la plaque d'extrémité12,13,14. En outre, la plupart des études ont employé des techniques de mesure de précision inférieure10,15,16. En outre, un grossissement significatif a été rapporté lorsque les paramètres de géométrie sont mesurés en utilisant la radiographie ou la tomographie calculée (CT)17,18. Bien que l'imagerie par résonance magnétique (IRM) soit considérée comme non invasive, elle est moins précise dans la définition des marges précises des structures ossées11. En raison de l'absence d'un protocole de mesure normalisé, il existe de grandes différences entre les données géométriques existantes.
Ces dernières années, l'ingénierie inverse, qui peut numériser les pièces physiques existantes en modèles solides informatisés, a été de plus en plus appliquée dans le domaine de la médecine. La technique permet de développer une représentation précise du caractère anatomique des surfaces vertébiennes sophistiquées. Le système d'ingénierie inverse comprend deux sous-systèmes : le système d'instrumentation et le système logiciel. Le système d'instrumentation adopté dans ce protocole est muni d'un scanner à plat optique de portée 3D sans contact, à grande vitesse et très précis (précision de 0,02 mm, 1 628 x 1 236 pixels). Le scanner peut capturer efficacement (temps d'entrée 3 s) des informations de morphologie de surface de l'objet cible et le convertir en nuage de points numériques. Le système logiciel (c.-à-d. logiciel d'ingénierie inversée) est une application informatique pour le traitement des données en nuage de points (voir Tableau des matériaux),la reconstruction du modèle de surface 3D, l'édition de courbes et de surface libres, et le traitement des données (voir tableau de Matériaux).
Les objectifs du présent rapport sont (1) de concevoir un protocole de mesure et un algorithme pour obtenir des paramètres quantitatifs des terminaisons vertébrales basés sur une technique d'ingénierie inversée, (2) développer un modèle mathématique qui permet un modèle réaliste représentation des plaques vertébrales sans numériser trop de repères. Ces méthodes seront bénéfiques à la planification de l'acte chirurgical et à la modélisation des éléments finis.
Subscription Required. Please recommend JoVE to your librarian.
Protocol
Cette étude a été approuvée par le comité d'éthique de la recherche en santé de l'institut des auteurs. Comme les os vertébraux cervicaux ont des formes plus complexes19, le protocole utilise les vertèbres cervicales comme une illustration pour faciliter la recherche pertinente.
1. Préparation des matériaux, numérisation et traitement d'image
- Recueillir une vertèbre cervicale sèche sans déformation pathologique ou pièces cassées.
- Placez la vertèbre verticalement dans la plate-forme du scanner (Figure 1, voir Tableau des matériaux), avec la plaque d'extrémité face à l'objectif de la caméra. Utilisez la source de lumière active du scanner. Ensuite, démarrez le processus de numérisation pour obtenir des données de cloud point (. format ASC).
REMARQUE : Selon les images de pré-scan, ajustez le scanner et la position de la vertèbre pour capturer autant d'informations de morphologie de surface que possible. - Ouvrez le logiciel spécialement utilisé pour le traitement des nuages de points (voir Tableau des matériaux). Cliquez sur Import pour importer les données du cloud point et générer le graphique numérique de la vertèbre. Définir le taux d'échantillonnage à 100%, sélectionnez Garder des données complètes sur l'échantillonnage, sélectionnez l'unité de données en tant que millimètres, et cliquez sur Points d'ombre. Utilisez l'outil de sélection Lasso pour sélectionner les points redondants sur le graphique, puis cliquez sur Supprimer pour les supprimer. Cliquez sur Réduire le bruit et définir le niveau de douceur à son maximum pour réduire le bruit et les pointes ( Figure2A, B).
REMARQUE : Il existe des instructions de fonctionnement logicielles de base au bas de l'interface graphique (interface utilisateur graphique). Les points de bruit avec des éperons pointus évidents latéralement ou verticalement doivent être enlevés pour réduire l'erreur. - Cliquez sur Envelopper pour emballer les données d'imagerie dans un fichier de format .stl pour transformer le nuage de points en maille, qui convertira un objet point en un objet polygone.
REMARQUE : Le logiciel d'ingénierie inversée accepte généralement le format 3D de type .stl. - Ouvrez le logiciel spécialement utilisé pour la reconstruction 3D et le traitement des données (voir Tableau des matériaux). Cliquez sur Fichier puis Nouveau dans le sous-menu. Sélectionnez partie dans la liste des types. Cliquez sur Démarrer, puis Forme zetez dans le sous-menu, puis Digitized Shape Editor. Cliquez sur l'icône Import dans la barre d'outils sur le côté droit de l'interface graphique. Dans la fenêtre Import, sélectionnez le fichier de format .stl, puis cliquez sur Appliquer 'gt; OK. Cliquez sur Fit All dans l'icône dans la barre d'outils en bas pour charger l'image reconstruite à la fenêtre principale du logiciel de présentation.
REMARQUE : Les étapes 1.5 à 2.3.3 sont exécutées avec le même logiciel. - Cliquez sur Activer dans la barre d'outils à droite. Dans la fenêtre Activer, sélectionnez Mode Piège 'gt; Type polygonal 'gt; Inside Trap. Ensuite, sélectionnez la plaque vertébrale sur l'image 3D pour enlever les composants vertébraux inutiles, tels que les éléments postérieurs et les ostéophytes (Figure 2C).
2. Quantification de la morphologie 3D de la plaque d'extrémité
- Définition du système de coordonnées 3D de la plaque d'extrémité
- Cliquez sur Démarrer 'gt; Forme dans le sous-menu, puis Generative Shape Design. Cliquez sur l'icône Point dans la barre d'outils à droite. Marquez trois repères anatomiques sur la jante épiphysher : les deux premiers sont les points d'extrémité gauche et droit du bord de fuite de la plaque d'extrémité, respectivement; le troisième est le point médian antérieur.
- Cliquez sur l'icône De ligne dans la barre d'outils sur le côté droit et sélectionnez les deux points de terminaison de bord de fuite pour définir une ligne frontale postérieure. Cliquez sur l'icône Plan, sélectionnez le type d'avion pour qu'il soit normal de la courbe, puis sélectionnez la ligne frontale postérieure et le point médian antérieur pour définir le plan mi-sagittal.
- Cliquez sur Démarrer 'gt; Forme 'gt; Reconstruction de surface rapide. Cliquez sur l'icône de la section planaire, entrez 1 dans l'option de numéro, puis sélectionnez l'image de la plaque d'extrémité et le plan mi-sagittal pour générer une courbe qui se croise. Cliquez sur Courbe de l'icône Scan et sélectionnez l'intersection de la courbe qui se croise et de la jante épiphysétique postérieure. Définissez l'intersection comme le point médian postérieur.
- Cliquez sur Démarrer 'gt; Forme 'gt; Generative Shape Design. Cliquez sur l'icône Ligne et sélectionnez le point médian antérieur et le point médian postérieur pour définir un diamètre mi-sagittal. Cliquez sur l'icône Point, puis Points et Planes Répétition dans le sous-menu. Ensuite, sélectionnez le diamètre mi-sagittal et entrez 1 dans l'option Instance(s) pour définir le point médian du diamètre mi-sagittal.
- Cliquez sur l'icône Axis System dans la barre d'outils en bas. Ensuite, sélectionnez le point médian du diamètre mi-sagittal comme origine, la ligne parallèle à la ligne frontale postérieure comme l'axe x, le diamètre mi-sagittal comme l'axe y, et la ligne pointant vers l'avant et perpendiculaireau au plan x-y comme l'axe z(Figure 3 ).
REMARQUE : Les deux points de terminaison de bord de fuite sont choisis comme points de référence parce qu'ils sont cohérents et montrent une variation minimale en présence d'ostéophytes10.
- Ajustement des courbes et des points caractéristiques sur la surface de la plaque d'extrémité (Figure 4A-D)
- Cliquez sur l'icône Point, puis Points et Planes Répétition dans le sous-menu. Sélectionnez le diamètre mi-sagittal et entrez 3 dans l'option Instance(s) pour diviser le diamètre mi-sagittal également en quatre parties.
- Cliquez sur Démarrer 'gt; Forme 'gt; Reconstruction de surface rapide. Cliquez sur l'icône de la section planaire, entrez 1 dans l'option Nombre, puis sélectionnez l'image de la plaque d'extrémité et le plan x-z pour générer une courbe qui se croise. Cliquez sur Courbe de l'icône Scan et sélectionnez les deux intersections du plan x-z et de la jante épiphysher.
- Définissez la ligne entre les deux intersections comme étant le diamètre mi-frontal. De la même manière, diviser le diamètre mi-frontal également en quatre parties.
REMARQUE : Lorsque la plaque d'extrémité n'est pas symétrique par rapport au plan médio-sagittal, choisissez l'un des deux points d'extrémité de la courbe mi-frontale qui a une distance verticale plus courte par rapport au plan z-y. Ensuite, définissez le diamètre mi-frontal comme 2 fois la longueur du plus court, et divisez-le également en quatre parties. - Cliquez sur la mesure entre l'icône dans la barre d'outils en bas pour mesurer la longueur d'un quart du diamètre mi-sagittal. Cliquez sur l'icône de la section planaire, entrez 2 dans l'option Nombre, entrez la valeur mesurée dans l'option Étape, puis sélectionnez l'image de la plaque d'extrémité et le plan x-z pour générer deux courbes d'ajustement d'un côté de la partie frontale. Cliquez sur Swap pour générer deux courbes d'ajustement de l'autre côté. De la même manière, obtenir les trois autres courbes d'ajustement dans le plan sagittal.
REMARQUE : Les deux courbes d'ajustement mi-frontales se chevauchent avec les deux courbes d'ajustement mi-sagittal. - Sélectionnez 11 points équidistants dans chaque courbe pour les mesures ultérieures. La méthode spécifique est la suivante:
- Si l'on prend l'exemple de la courbe mi-sagittale, divisez également le diamètre mi-sagittal en 10 parties, ce qui se traduit par une somme de 11 points, dont neuf points intermédiaires et deux points d'évaluation (se référer aux étapes 2.1.3 et 2.2.1).
- Passez par chaque point équidistant, obtenez neuf courbes d'ajustement sur la surface de la plaque d'extrémité (voir l'étape 2.2.2). Cliquez sur Courbe de l'icône Scan et sélectionnez l'intersection des courbes d'ajustement et de la courbe mi-sagittale. Enfin, obtenir un total de 66 points sur chaque plaque d'extrémité (11 points par courbe multipliée par six courbes). Cliquez sur l'icône mesurez l'élément dans la barre d'outils en bas pour mesurer les coordonnées de chaque point.
- Mesure des paramètres morphologiques de la plaque d'extrémité
- Paramètre de ligne :
- Cliquez sur la mesure entre l'icône pour mesurer la longueur du paramètre de ligne qui est la distance entre deux points mesurés.
- Paramètres de concavité :
- Créer un plan parallèle à l'avion x-y (Figure 5A):cliquez sur Démarrer 'gt; Forme 'gt; Generative Shape Design. Cliquez sur l'icône Sketch dans la barre d'outils sur le côté droit, puis cliquez sur le plan x-y. Cliquez sur l'icône Cercle, cliquez sur Origine sur la surface de la plaque d'extrémité, faites glisser le curseur de la souris à une distance appropriée, puis cliquez. Cliquez sur l'icône Du banc de sortie, puis l'icône Remplir, puis cliquez.
- Cliquez sur l'icône Offset, sélectionnez le plan rempli et entrez une valeur appropriée dans l'option offset jusqu'à ce qu'elle soit tangente à la partie la plus concave, et zoomez. Cliquez sur Démarrer 'gt; Forme 'gt; Reconstruction de surface rapide. Ensuite, cliquez sur l'icône de courbe 3D pour trouver et créer le point le plus concave. Cliquez sur l'icône Mesure zona pour mesurer les coordonnées du point le plus concave (figure 5B).
- Cliquez sur la mesure entre l'icône, puis sélectionnez le point le plus concave et le plan x-y pour mesurer toute la profondeur de la concavité de la plaque d'extrémité. De même, trouver et créer la profondeur la plus concave sur un plan particulier et mesurer ses coordonnées.
- Cliquez sur l'icône Projection dans la barre d'outils sur le côté droit, puis sélectionnez le point le plus concave et le plan x-y pour obtenir le point projectif. Cliquez sur l'icône Mesure zona pour mesurer les coordonnées du point de projection et déterminer sa distribution en fonction des coordonnées.
- Paramètres de surface :
- Cliquez sur l'icône d'inertie de mesure dans la barre d'outils en bas et cliquez sur la surface de la plaque d'extrémité pour mesurer sa surface. Cliquez sur l'icône Activer et sélectionnez la plaque d'extrémité centrale le long des marges intérieures de l'anneau épiphysélausé (voir l'étape 1.6), puis cliquez sur l'icône d'inertie de mesure pour mesurer sa zone (Figure 5C). Cliquez sur l'icône Activer, puis la plaque d'extrémité centrale, et enfin l'icône Swap dans la fenêtre Activer pour obtenir une jante épiphysher. Ensuite, mesurez sa superficie.
- Paramètre de ligne :
3. Développement du modèle mathématique de surface de surface de plaque d'extrémité
- Détermination de l'ordre d'ajustement de l'équation paramétrique
- Ouvrez le logiciel d'analyse et de visualisation des données (voir Tableau des matériaux). Entrée x [données correspondantes] dans la fenêtre de commande. Cliquez sur Entrez.
REMARQUE : Les « données correspondantes » se réfèrent aux données x-coordonnées des 11 points caractéristiques d'une courbe qui ont été mesurées dans les étapes précédentes. Cliquez sur Entrez après la saisie de chaque commande, la même application aux opérations ultérieures. Les étapes 3.1 à 5.5 sont exécutées uniformément avec le même logiciel. - De la même manière, l'entrée z '[données correspondantes].
- Entrée du code pour i'1:5 z2'polyfit(x,z,i); Z-polyval (z2,x); si la somme (Z-z). fin; fin.
REMARQUE : Le protocole fixe la somme d'erreur des carrés en dessous de 0,01 pour obtenir une précision plus élevée, dont la valeur peut être réajustée pour satisfaire diverses demandes. - Cliquez sur Entrez pour obtenir une valeur C qui est l'ordre d'ajustement souhaité.
- Ouvrez le logiciel d'analyse et de visualisation des données (voir Tableau des matériaux). Entrée x [données correspondantes] dans la fenêtre de commande. Cliquez sur Entrez.
- Ajustement d'équation de paramètre
- Input cftool et cliquez sur Entrez pour mettre en place l'outil de montage de courbe.
- Entrez les coordonnées d'une courbe dans la fenêtre de commande (se référer aux étapes 3.1.1 et 3.1.2). Dans l'outil curve Fitting, sélectionnez les données x-coordonnées lors de l'ajustement des courbes d'avion frontales et des données y-coordinate lors de l'ajustement des courbes d'avion sagittal dans l'option de données x, sélectionnez les données z-coordonnées dans l'option de données y, sélectionnez polynomial, et entrez l'ordre d'ajustement Obtenu. Ensuite, le logiciel produira automatiquement l'équation paramétrique et la bonté de l'ajustement.
REMARQUE : Comme la courbe est une image 2D, l'option de travail par défaut est les options x et y dans l'outil d'ajustement de courbe lors de l'ajustement d'une courbe. - De la même manière, entrez les coordonnées 3D des 66 points et faites correspondre les données de coordonnées aux options d'axe correspondantes. Sélectionnez polynomial et entrez dans l'ordre d'ajustement pour obtenir l'équation paramétrique de la surface de la plaque d'extrémité (figure 6B).
4. Acquisition de données géométriques basées sur l'équation paramétrique
- Entrée x- et y-coordonner les valeurs de n'importe quel point sur la plaque d'extrémité dans la fenêtre de commande.
- Entrée PX1, PX2, PX3....
REMARQUE: Px est les paramètres de l'équation paramétrique qui ont été installés à l'aide polynomiale dans les étapes ci-dessus. - Entrez l'équation et cliquez Sur Entrez pour obtenir le résultat (c.-à-d., format d'entrée : z ' P00 ' P10'x 'P01'y ' P20'x '2 ' P11'x 'y 'p02'y '2 ' P30'x '3 ' P21'x '2 'y ' P12 P2 p. 3,P40x 4, P31x 3,5 , P22x 2,p.
5. Représentation de la plaque d'extrémité basée sur l'équation paramétrique
- Entrée PX1, PX2, PX3.... dans la fenêtre de commande.
- Entrée du code X-N1:0.01:N2;.
REMARQUE : N1-N2 est la gamme de données de l'axe X (c.-à-d., les valeurs des deux points de terminaison de la courbe du milieu de la couronne). - Entrée du code "Y-N 3:0.01:N4;".
- Entrée de l'équation (c.-à-d., z '(x,y)P00 ' P10.x 'P01. . . ' P20. . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 03. . . . . 3 P40. . x.4 p31. . . . . . . .. P22. . .. .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
- Entrée du code ezmesh(z, [N1,N2,N3,N4]) pour obtenir des graphiques de simulation 3D (Figure 6C).
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
À l'aide du scanner à plate-forme optique 3D très précis, les plaques d'extrémité ont été converties en plus de 45 000 points numériques, qui caractérisent adéquatement la morphologie (Figure 2A,B).
Dans le protocole de mesure, l'analyse spatiale des surfaces d'extrémité a été effectuée. Des courbes représentatives ont été ajustées et quantifiées à la surface pour caractériser la morphologie (figure 4B). Les paramètres linéaires ont été mesurés en calculant la distance entre deux points de terminaison. Les mesures obtenues comprennent la profondeur de la cavité et l'emplacement de l'apex de la cavité dans le plan mi-sagittal, en plus de ceux de la concavité de la plaque d'extrémité entière et de toute section spécifique (figure 5B). Les composants des plaques d'extrémité, de la jante épiphysher et de la plaque d'extrémité centrale ont été séparés (figure 5C), et leurs longueurs et leurs zones ont été obtenues commodément.
Un total de 138 plaques vertébrales cervicales ont été numérisées et analysées, et le modèle mathématique de la plaque d'extrémité a été établi. Le protocole fixe les sommes d'erreur au carré en dessous de 0,01, et il a été conclu que l'utilisation de la fonction polynomiale à quatre commandes pourrait atteindre la satisfaction.
L'équation paramétrique de chaque courbe a été déduite en fonction des coordonnées de 11 points : f(x) ' P1x '4 ' P2'x '3 ' P3'x '2 ' P4x 'P'x 'P 5 . P1, P2, P3, P4 et P5 étaient les paramètres, dont les valeurs exactes sont indiquées dans le tableau 1.
L'équation paramétrique représentant les caractéristiques morphologiques de la surface de la plaque d'extrémité est la suivante :
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 ' P03'y '3 ' P4'x'4' P 31 'x '3 'y 'P22 X '2 'y '2 ' P13'x 'y '3 ' P04'y '4
Où: PXYs sont les paramètres, qui ont été déduits des coordonnées pré-mesurées de 66 points (tableau 2).
Figure 1 : Scanner à plat optique 3D sans contact. Le scanner, basé sur la technologie de mesure optique 3D de décalage de phase 3D hétérodérine, comprend la mesure optique (intégration autour de deux caméras et d'un projecteur) et des dispositifs de contrôle. La précision de cet instrument est de 0,02 mm, et les pixels sont 1628 x 1236. Le scanner peut numériser efficacement (temps d'entrée 3 s) la géométrie de surface d'un objet cible. Veuillez cliquer ici pour voir une version plus grande de ce chiffre.
Figure 2 : Nuage de pointe de surface vertébrale et reconstruction 3D de la plaque d'extrémité. (A) et (B) sont les surfaces inférieures et supérieures d'une vertèbre cervicale générée par le logiciel spécialement utilisé pour le traitement des nuages de points, respectivement. (C) et (D) sont la reconstruction 3D des plaques d'extrémité inférieures et supérieures générées par le logiciel spécialement utilisé pour la reconstruction 3D et le traitement des données, respectivement. Les éléments postérieurs et les ostéophytes sont retirés des vertèbres, ne laissant que la plaque d'extrémité. Le plan le mieux adapté est défini à travers les points antérieurs et postérieurs les plus des processus bilatéraux incinate, et les deux courbes formées par le plan le mieux adapté et la plaque d'extrémité sont les limites de l'articulation et de la plaque d'extrémité caudale. Veuillez cliquer ici pour voir une version plus grande de ce chiffre.
Figure 3 : Définition du système de coordonnées 3D de la plaque d'extrémité. Marquage de trois repères anatomiques sur la jante épiphysher : les deux premiers sont les extrémités gauche et droite du bord de fuite de la plaque d'extrémité, respectivement; le troisième est le point médian antérieur. La ligne frontale postérieure est formée par les deux extrémités de bord de fuite, qui définissent le plan mi-sagittal avec le point médian antérieur. Le point médian postérieur est déterminé par le plan mi-sagittal et la jante épiphyseal postérieure, qui forment le diamètre mi-sagittal avec le point médian antérieur. L'origine est le point médian du diamètre mi-sagittal. L'axe y est déterminé par le diamètre mi-sagittal et pointant vers l'avant. L'axe X est la ligne parallèle à la ligne frontale postérieure. L'axe z est normal pour le plan x-y. Veuillez cliquer ici pour voir une version plus grande de ce chiffre.
Figure 4 : Les étapes de l'ajustement des courbes et des points caractéristiques sur la surface de la plaque d'extrémité. (A)Diviser le diamètre mi-sagittal et le diamètre mi-frontal également en quatre parties. (B) Passez par chaque point équidistant, et choisissez six courbes de surface symétriquement, dont trois sont les courbes d'intersection du plan frontal et de la surface de la plaque d'extrémité, et les trois autres dans le plan sagittal. (C) Diviser le diamètre mi-sagittal également en 10 parties. (D) En passant par chaque point équidistant, les plans frontaux et la courbe mi-sagittale forment neuf intersections, ce qui donne une somme de 11 points, ainsi que les deux points d'arrivée. Veuillez cliquer ici pour voir une version plus grande de ce chiffre.
Figure 5 : Mesure de la profondeur et de la surface de la concavité de la plaque d'extrémité. (A) Créer un plan parallèle au plan x-y. (B) Dépensez l'avion jusqu'à ce qu'il soit tangente jusqu'au point le plus concave, et la profondeur de concavité de la plaque d'extrémité est la distance perpendiculaire entre le point le plus concave et le plan x-y. (C) Tracer une ligne le long des marges intérieures de l'anneau épiphysélet pour diviser la plaque d'extrémité dans la plaque d'extrémité centrale et la jante épiphyslée. Veuillez cliquer ici pour voir une version plus grande de ce chiffre.
Figure 6 : La reconstruction 3D et les représentations d'une plaque d'extrémité inférieure. (A) La reconstruction 3D de la surface inférieure de la plaque d'extrémité générée par le logiciel spécialement utilisé pour la reconstruction 3D et le traitement des données. (B) et (C) sont les représentations de la plaque d'extrémité inférieure générée par le logiciel d'analyse et de visualisation des données. Veuillez cliquer ici pour voir une version plus grande de ce chiffre.
| Niveau Endplate | courbe | Paramètres | ||||
| P1 P1 (en) | P2 (en) | P3 (en) | P4 (en) | P5 P5 (en) | ||
| Supérieur C6 | FAC (fac) | 0 | 0 | -0.0128 | -0.0028 | 0.02523 |
| Fmc | 0 | 0 | -0.0199 | 0.00074 | 0.3693 | |
| FPC (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 | |
| SRC (SRC) | 0 | 0.00179 | -0.0096 | 0.04451 | -0.0394 | |
| C6 inférieur | FAC (fac) | 0 | -0.0001 | -0.0225 | 0.00594 | 1.223 |
| Fmc | 0 | 0 | -0.016 | -0.0082 | 1.729 | |
| FPC (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 | |
| SRC (SRC) | 0 | 0.00079 | -0.0295 | -0.0828 | 1.362 | |
Tableau 1 : Paramètres de l'équation pour représenter la courbe de la surface de la plaque d'extrémité. Seules les données de la sixième plaque vertébrale cervicale sont répertoriées. Px - les paramètres de l'équation. Sur chaque plaque d'extrémité, six courbes de surface ont été choisies symétriquement; trois d'entre eux se trouvaient dans le plan frontal et appelaient la courbe antérieure (FAC), la courbe moyenne (FMC) et la courbe postérieure (FPC); les trois autres dans le plan sagittal ont été appelés la courbe gauche (SLC), la courbe moyenne (SMC), et la courbe droite (SRC). Les paramètres d'une valeur absolue inférieure à 0,0001 sont représentés comme 0 ici.
| Paramètres | C3 inf | Sup C4 | C4 inf | Sup C5 | C5 inf | Sup C6 | C6 inf | Sup C7 |
| p00 (en) | 1.989 | 0.4187 | 2.004 | 0.3383 | 1.913 | 0.4276 | 1.779 | 0.5674 |
| p10 (en) | -0.0022 | -0.0043 | 0.00542 | -0.0208 | -0.0111 | 0.0012 | -0.0043 | -0.0052 |
| p01 (en) | -0.0356 | -0.0868 | -0.0537 | -0.0826 | -0.0257 | -0.098 | -0.0407 | -0.0642 |
| p20 (en) | 0.01286 | -0.0252 | -0.0146 | -0.0299 | -0.0253 | -0.0264 | -0.0175 | -0.0088 |
| p11 (en) | 0.00092 | 0.00071 | -0.0009 | 0.00018 | -0.0002 | -0.0012 | 0.00117 | 0.00021 |
| p02 (en) | -0.0529 | -0.0151 | -0.0525 | -0.012 | -0.0418 | -0.0142 | -0.0396 | -0.0134 |
| p30 (en) | 0 | -0.0001 | 0.00013 | 0.00024 | 0.00017 | 0 | 0 | 0 |
| p21 (en) | -0.0011 | 0.00299 | -0.0012 | 0.00363 | -0.0021 | 0.00306 | -0.0019 | 0.00194 |
| p12 (en) | 0 | 0.00048 | -0.0004 | 0.00033 | 0.00014 | 0 | -0.0001 | 0 |
| p03 (en) | 0.00062 | 0.00204 | 0.00089 | 0.00206 | 0.00046 | 0.00208 | 0.00077 | 0.00115 |
| p40 (en) | 0.0002 | 0 | 0.0002 | 0 | 0.00024 | 0 | 0 | 0 |
| p31 (en) | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| p22 (en) | 0.00017 | 0.00013 | 0 | 0.00015 | 0.00015 | 0.00017 | 0.00032 | 0 |
| p13 (en) | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| p04 (en) | 0.00023 | 0.00013 | 0.00024 | 0 | 0 | 0 | 0 | 0 |
Tableau 2 : Paramètres de l'équation paramétrique représentant la morphologie de la surface de la plaque d'extrémité. Px - les paramètres de l'équation; inf - plaque d'extrémité inférieure; sup - plaque d'extrémité supérieure. Les paramètres d'une valeur absolue inférieure à 0,0001 sont représentés comme 0 ici. Ce tableau a été modifié à partir d'une publication précédente3.
| Mesures | Fiabilité intratest | Mesures | RE vs Caliper | ||
| APD (en) | Première-remeasurement | 15,76 à 1,3 | APD (en) | concernant | 16,47 à 1,31 |
| Remesure | 15,86 à 1,61 | Étrier | 16,26 à 1,27 | ||
| Icc | 0.85 | Cronbach alpha | 0.99 | ||
| Cmd | Première-remeasurement | 19,71 à 2,47 | Cmd | concernant | 20,7 à 3,05 |
| Remesure | 19,41 à 2,43 | Étrier | 20,45 à 3,21 | ||
| Icc | 0.96 | Cronbach alpha | 0.99 | ||
Tableau 3 : Fiabilité des mesures. Les données étaient moyennes et standard écart (mm). ICC - coefficient de corrélation intra-classe; APD - diamètre antéro-postérieur; CMD - diamètre médiolatéral central; RE - le système d'ingénierie inverse. Ce tableau a été modifié par rapport à une publication précédente. 3 (en)
| Valeur des mesures | n | Z valeur de coordonnées | T | page | R |
| Points originaux | 15 | 1,75 à 0,87 | 0.26 | 0.8 | 0.98 |
| Points de comparaison | 15 | 1,74 à 0,91 |
Tableau 4 : Validité du modèle géométrique représentant la morphologie de la plaque d'extrémité. Les données sont représentées comme une déviation moyenne et standard (mm). Les points d'origine sont 15 points choisis au hasard sur l'image de reconstruction 3D originale. Points de comparaison - points correspondants générés auto-générés à partir d'équations paramétriques; Coefficient de corrélation R.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
L'ingénierie inverse a été appliquée de plus en plus et avec succès dans le domaine de la médecine, comme la cranioplastie20, orale21, et les implants maxillo-faciaux21. Les mesures d'ingénierie inversées, à savoir la numérisation de surface du produit, se réfèrent à la conversion des informations de surface en données de cloud de point utilisant des équipements et des méthodes de mesure spécifiques. Sur la base de ces données, la modélisation, l'évaluation, l'évaluation, les améliorations et la fabrication de surface complexes peuvent être effectuées. La mesure numérique et le traitement des données sont une technologie de base et clé utilisée dans l'ingénierie inverse.
Dans ce protocole, des informations morphologiques précises et détaillées sur les plaques vertébrales sont enregistrées à l'aide d'un système de balayage optique 3D sans contact, basé sur la multifréquence hétérodène, le changement de phase et la technologie de mesure optique 3D. Le scanner est principalement composé d'appareils de contrôle et d'une mesure optique intégrant deux caméras et un projecteur. Comparé à d'autres instruments de mesure, le scanner est très précis et efficace et évite la numérisation point par point. Lors de la capture de données point-nuage, la tête de balayage n'est généralement pas en contact avec l'objet, de sorte qu'il n'y a pas d'effets de déformation. La fiabilité, la validité et la précision du scanner pour l'enregistrement de la morphologie de surface ont été bien établies2,3,22. La reproductibilité de ces mesures a été vérifiée.
Pour vérifier l'exactitude des mesures prises par le système d'ingénierie inverse, 20 plaques d'extrémité ont été mesurées à l'aide d'un étrier numérique et évaluées à l'aide de Cronbach alpha. Pour la fiabilité intra-test, 16 plaques d'extrémité ont été choisies au hasard à partir des 138 plaques vertébrales et mesurées deux fois à intervalles de 2 semaines, puis évaluées à l'aide d'un coefficient de corrélation intraclasse. Les résultats ont montré un grand accord et fiabilité (tableau 3). Le logiciel d'ingénierie inversée implique des mesures puissantes, le traitement des données, la détection d'erreurs et les fonctions d'édition de courbes et de surface gratuites. Il peut également construire et ajuster intelligemment et efficacement des courbes et des surfaces, et la reconstruction 3D de modèle de surface contribue aux mesures précises23.
Il existe des applications importantes et considérables pour les données d'anatomie détaillées et complètes des vertèbres, telles que la conception d'implants rachidiens, l'amélioration de la fidélité des modèles d'éléments finis de la colonne vertébrale, et le développement de modèles mathématiques. La plaque vertébrale est essentielle au maintien de l'intégrité et de la fonction du disque intervertébral, et elle sert également d'interface mécanique pour transférer le stress. Par conséquent, la quantification de la géométrie de la plaque d'extrémité est importante. Avec l'aide de l'ingénierie inverse, la morphologie de la plaque d'extrémité peut être quantifiée intelligemment et de manière exhaustive. Dans ce protocole, six courbes caractéristiques sont installées sur la surface de chaque plaque d'extrémité, et un système de coordonnées 3D est mis en place pour quantifier la morphologie spatiale.
En outre, un modèle paramétrique de la plaque d'extrémité est développé pour instituer des évaluations quantitatives précises et reproductibles et pour développer des modèles d'éléments finis biomécaniques personnalisés. Le modèle paramétrique de la surface des plaques d'extrémité peut produire des représentations rapides, réalistes et précises qui peuvent être visualisées et commodément analysées par les chercheurs.
L'inclusion d'un plus grand nombre de repères améliorera la précision, mais elle prend beaucoup de temps et coûte cher. Dans ce protocole, il est proposé que 66 points de six courbes de surface soient suffisants pour décrire les caractéristiques morphologiques. Les tests de fiabilité sont également effectués en comparant les valeurs de coordonnées de 15 points choisis au hasard avec les valeurs correspondantes qui sont générées auto à partir d'équations paramétriques. Le résultat révèle que le modèle paramétrique a une bonne fiabilité et une bonne reproductibilité peut servir de représentation réaliste de la surface de la plaque d'extrémité (tableau 4). Il convient de noter que le modèle paramétrique peut être dérivé en fonction d'autres modalités d'imagerie telles que la tomodensitométrie et l'IRM.
Comme les scanners sans contact sont sensibles à la lumière ambiante, il est essentiel de maintenir la lumière ambiante stable, et des sources de lumière actives sont recommandées. S'il y a de la graisse résiduelle sur la surface de la plaque d'extrémité, la poudre infantile de talc doit être barbouillée doucement pour éviter le risque d'être affectée par les caractéristiques de réflectoflant spatial de la surface de l'objet. Les vertèbres cervicales sous-axiales ont une composante spéciale : l'articulation non-covertebral. Pour le distinguer de la plaque d'extrémité, un plan de meilleure qualité est défini en utilisant la méthode la moins carrée. Ensuite, la courbe d'intersection formée par le plan le mieux adapté, et la surface de la plaque d'extrémité est la frontière entre l'articulation non-covertebral et la plaque d'extrémité supérieure (Figure 2D).
L'opération spécifique est la suivante: cliquez sur Démarrer 'gt; Forme 'gt; Generative Shape Design. Cliquez sur l'icône Point dans la barre d'outils à droite, puis sélectionnez les points les plus antérieurs et les plus postérieurs des processus bilatéraux non cinés sur l'image 3D. Cliquez sur l'icône Plan et sélectionnez Mean Through Points dans le type d'avion pour définir le plan le mieux adapté. Cliquez sur Démarrer 'gt; Forme 'gt; Reconstruction de surface rapide. Cliquez sur l'icône De la Section Planaire, puis sélectionnez l'image 3D et le plan le mieux adapté.
Le marquage précis des trois points anatomiques sur la surface de la plaque d'extrémité lors de l'établissement du système de coordonnées 3D est essentiel. Le logiciel d'ingénierie inverse permet un déplacement flexible de l'image de reconstruction et améliore le contraste qui aide à identifier les points de repère. Alternativement, il est important d'évaluer la pertinence du système de coordonnées en fonction de la question de savoir si la ligne d'intersection des plans mi-sagittal et coronal défini est perpendiculaire à la section de la plaque d'arrivée, puis d'ajuster le système en conséquence. Des essais intra-observateurs ont également été évalués, et le résultat a indiqué une bonne fiabilité (tableau 3).
Ce protocole nécessite de multiples compétences et techniques, y compris l'acquisition et le traitement de données en nuage de points, la reconstruction et l'analyse d'images et le développement de modèles paramétriques. Pour un débutant, il peut prendre du temps pour terminer l'ensemble du processus. Cependant, comme seuls quelques modules du logiciel dans ce protocole sont utilisés et la procédure est modulaire, il faut une courte courbe d'apprentissage pour devenir bien expérimenté.
En conclusion, le protocole décrit fournit une méthode précise et reproductible pour obtenir des données géométriques détaillées et complètes des plaques vertébrales. Un modèle paramétrique est également développé sans numériser trop de repères, ce qui est bénéfique pour concevoir des implants rachidiens personnalisés, planifier des actes chirurgicaux, poser des diagnostics cliniques et développer des modèles d'éléments finis précis.
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
Les auteurs ne déclarent aucun intérêt financier concurrent.
Acknowledgments
Ces travaux ont été financés par Key Discipline Construction Project du Pudong Health Bureau de Shanghai (PWZxk2017-08) et par la National Natural Science Foundation of China (81672199). Les auteurs tiens à remercier Wang Lei pour son aide dans la relecture d'une version antérieure et Li Zhaoyang pour son aide dans le développement du modèle paramétrique.
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).
