FR 3161473 A1 20251024 FR 2404139 A 20240422 Procédé (100) de détermination d’une dimension (D) d’une pièce mécanique (PM1, PM2), comprenant les étapes de mesure (101) d’une position de points sur un bord de la section de la pièce, détermination (102) d’une courbe squelette (CSQ), comprenant les étapes de calcul (1021) d’une triangulation des points de la courbe section (CS), et pour chaque triangle de la triangulation, de détermination (1022) d’un centre (CT), et d’un barycentre (BT) dudit triangle, de filtrage (1023) des centres (CT), sur la base des barycentres (BT), pour obtenir un ensemble de centres filtrés (CTF), la courbe squelette (CSQ) déterminée comprenant l’ensemble des centres filtrés (CTF), et de détermination (103) de la dimension (D) de la pièce mécanique, en un point squelette de la courbe squelette (CSQ), en fonction de la position dudit point squelette sur la courbe squelette (CSQ). Figure 2 Procédé de détermination d’une dimension d’une pièce mécanique La présente invention concerne le domaine de la mesure dimensionnelle de pièce mécanique, en particulier par ordinateur, en particulier pour le contrôle qualité dans les processus de fabrication industrielle. Il est connu de mesurer une dimension d’une pièce mécanique à partir de points positionnés sur une section de la pièce mécanique. L’inconvénient des procédés connus est qu’ils sont complexes et que leur fonctionnement est difficile à maîtriser précisément. Il y a donc un besoin de simplification, notamment pour permettre une meilleure maîtrise voire une amélioration des performances. L’invention a donc pour but de proposer une solution à tout ou partie de ces problèmes. A cet effet, la présente invention concerne un procédé de détermination d’une dimension d’une pièce mécanique, le procédé étant mis en œuvre par ordinateur et comprenant les étapes suivantes : - mesure d’une position, dans un plan d’une section de la pièce, d’une pluralité de points positionnés sur un bord de la section de la pièce, le bord de la section de la pièce formant une courbe section dans ledit plan, - détermination d’une courbe squelette comprenant une pluralité de points squelettes, la détermination de la courbe squelette comprenant les étapes suivantes : - calcul d’une triangulation de la pluralité de points de la courbe section, la triangulation comprenant un ensemble de triangles, chaque triangle de l’ensemble de triangles étant formé par trois points de la pluralité de points de la courbe section, - pour chaque triangle de la triangulation, détermination d’un centre dudit triangle, défini par le centre d’un cercle circonscrit audit triangle, et d’un barycentre dudit triangle, défini par le barycentre des sommets dudit triangle, de manière à obtenir un ensemble de barycentres et un ensemble de centres, chaque triangle de la triangulation correspondant à un barycentre de l’ensemble de barycentres et à un centre de l’ensemble de centres, - filtrage des centres, sur la base des barycentres, pour obtenir un ensemble de centres filtrés, la courbe squelette déterminée comprenant l’ensemble des centres filtrés, - détermination de la dimension de la pièce mécanique, en un point squelette de la courbe squelette, en fonction de la position dudit point squelette sur la courbe squelette. Selon un mode de mise en oeuvre, l’invention comprend une ou plusieurs des caractéristiques suivantes, seules ou en combinaison techniquement acceptable. Selon un mode de mise en œuvre, la position des points de la courbe section est mesurée avec un palpeur, ou un capteur radiographique, ledit palpeur ou capteur radiographique étant contrôlé par l’ordinateur, l’ordinateur étant configuré pour recevoir la mesure réalisée par le palpeur, ou capteur radiographique. Selon un mode de mise en œuvre, la triangulation est une triangulation de Delaunay, dans laquelle le cercle circonscrit à un triangle de la triangulation ne contient à l’intérieur de son périmètre aucun sommet d’un autre triangle de la triangulation. Selon un mode de mise en œuvre, ledit point squelette est positionné à une première distance d’une première partie de la courbe section, et à une deuxième distance d’une deuxième partie de la courbe section, la première distance et la deuxième distance étant mesurées transversalement à la courbe squelette, la dimension de la pièce mécanique audit point squelette de la courbe squelette étant égale à une somme de la première distance et de la deuxième distance. Selon un mode de mise en œuvre, une différence entre la première distance et la deuxième distance rapportée à une demi somme de la première distance et de la deuxième distance est inférieure à 1/100, de préférence inférieure à 1/1000, de préférence égale à 0. Selon ces dispositions, les points de la courbe squelette sont équidistants des points de la première partie et de la deuxième partie de la courbe section. Selon un mode de mise en œuvre, la première partie et la deuxième partie de la courbe section sont respectivement de part et d’autre d’une zone de la courbe section formant une cavité de la courbe section. Selon un mode de mise en œuvre, la première distance et la deuxième distance sont mesurées selon une direction transversale à la courbe squelette en un point squelette déterminé de la courbe squelette situé à une distance de la cavité de la courbe section, la direction transversale à la courbe squelette étant déterminée à partir d’un sous-ensemble de points squelettes positionnés autour dudit point squelette déterminé. Selon un mode de mise en œuvre, un nombre de points du sous-ensemble de points squelettes est prédéterminé. Selon un mode de mise en œuvre, l’étape de filtrage comprend : - la détermination qu’un barycentre d’un triangle est un barycentre incohérent, et que le triangle correspondant au barycentre incohérent est un triangle incohérent, et que le centre du triangle incohérent est un centre incohérent, - la suppression dudit centre incohérent de l’ensemble des centres, pour obtenir l’ensemble de centres filtrés. Selon un mode de mise en œuvre, l’étape de détermination de la courbe squelette comprend une extrapolation de points complémentaires à partir de l’ensemble de centres filtrés, de sorte que la courbe squelette comprend les points complémentaires et l’ensemble de centres filtrés. Selon un mode de mise en œuvre, l’extrapolation comprend l’utilisation d’un autre sous-ensemble de points squelettes, de la courbe squelette, pour estimer une tangente à la courbe squelette, ladite tangente définissant une direction d’extrapolation de la courbe squelette vers une cavité de la courbe section. Selon un mode de mise en œuvre, le nombre de points de l’autre sous-ensemble de points squelettes, utilisés pour estimer la direction d’extrapolation de la courbe squelette, est prédéterminé. Selon un mode de mise en œuvre, les points complémentaires sont extrapolés le long d’un morceau de courbe, le morceau de courbe étant par exemple linéaire, ou quadratique. Selon un mode de mise en œuvre, l’étape d’extrapolation de points complémentaires comprend un ordonnancement des points de l’ensemble de centres filtrés. Selon un mode de mise en œuvre, l’étape de détermination qu’un barycentre d’un triangle est un barycentre incohérent comprend : - la construction, à partir des points de la courbe section, d’un polygone de filtrage; - la détermination qu’un barycentre est incohérent si il est à l’extérieur dudit polygone de filtrage. Selon un mode de mise en œuvre, l’étape de construction du polygone de filtrage, comprend la détermination d’un ensemble de sommets dudit polygone de filtrage, chaque sommet étant positionné dans une direction définie à partir d’un sous-ensemble de points de la courbe section, à une distance prédéterminée dudit sous-ensemble de points de la courbe section, la direction définie à partir du sous-ensemble de points de la courbe section étant transversale à une portion de la courbe section passant par les points dudit sous-ensemble de points de la courbe section. Selon un mode de mise en œuvre, la direction définie à partir du sous-ensemble de points de la courbe section est définie par un vecteur propre d’une matrice de covariance calculée à partir des coordonnées des points du sous-ensemble de points de la courbe section. Selon un mode de mise en œuvre, la matrice de covariance s’écrit M=A T *A, où A est une matrice à N lignes et 2 colonnes, N étant un nombre de points du sous-ensemble de points considéré, la première et la deuxième colonne de chaque ligne comprenant respectivement les coordonnées x et y de chaque points du sous-ensemble de points considéré, A T étant une transposée de A, de sorte que M étant une matrice 2x2. Selon un mode de mise en œuvre, la courbe section est fermée de façon à délimiter un intérieur de la courbe section, et dans lequel le sommet du polygone de filtrage est positionné dans la direction transversale à la portion de la courbe section passant par les points dudit sous-ensemble de points de la courbe section, de façon qu’un segment orienté de ladite portion de la courbe section vers ledit sommet, selon ladite direction transversale, est orienté vers l’intérieur de la courbe section. Selon un mode de mise en œuvre, le procédé comprend en outre une étape d’inversion de l’orientation dudit segment orienté, lorsque ledit segment orienté coupe la courbe section en un point d’intersection et lorsqu’un point intermédiaire se trouve à l’extérieur de la courbe section, ledit point intermédiaire étant positionné entre un point origine du segment orienté et ledit point d’intersection, ledit point origine étant situé sur ladite portion de la courbe section, ladite inversion de l’orientation générant un segment orienté inversé, de sens opposé au segment orienté, ledit segment orienté inversé étant dirigé vers l’intérieur de de la courbe section. Selon un autre aspect, l’invention concerne un produit programme d’ordinateur comprenant des instructions de code mettant en œuvre un procédé de détermination de dimension de pièce tel que défini ci-avant, lorsque ces instructions de code sont exécutées par une unité de traitement. Selon un autre aspect, l’invention concerne également un support lisible par ordinateur sur lequel sont enregistrées des instructions de code mettant en œuvre un procédé de détermination de dimension de pièce tel que défini ci-avant. Pour sa bonne compréhension, un mode de réalisation et/ou de mise en oeuvre de l’invention est décrit en référence aux dessins ci-annexés représentant, à titre d’exemple non limitatif, une forme de réalisation ou de mise en œuvre d’un procédé selon un exemple de l’invention. Les mêmes références sur les dessins désignent des éléments similaires ou des éléments dont les fonctions sont similaires. FIG. 1 est une vue en perspective d’une première pièce mécanique, par exemple un bord d’attaque d’aube de turbine, et d’une deuxième pièce mécanique, par exemple une aube de turbine, la première pièce mécanique étant destinée à être montée en superposition sur une surface externe de la deuxième pièce mécanique. FIG. 2 est une vue d’une courbe section d’une pièce mécanique dans un plan de section, et d’une courbe squelette de la courbe section. FIG. 3 est une illustration de deux triangulations différentes de 4 points, avec une première triangulation sur la gauche de la FIG. 3 qui n’est pas une triangulation de Delaunay, et avec une deuxième triangulation sur la droite de la FIG. 3 qui est une triangulation de Delaunay. FIG. 4 est un exemple de triangulation de Delaunay construite à partir des points d’une courbe section, avec les barycentres des triangles représentés par des points ronds et les centres des triangles représentés par les étoiles. FIG. 5 est un exemple de polygone de filtrage construit à partir des points de la courbe section de la figure précédente, les barycentres des triangles situés à l’extérieur du polygone de filtrage ayant été éliminés, ainsi que les centres des cercles circonscrits à ces triangles, de sorte qu’il ne reste plus que les centres des triangles filtrés selon un exemple de l’invention, qui forment la courbe squelette de la courbe section. FIG. 6 est un exemple de sous-ensemble de points de la courbe section utilisé pour déterminer un sommet du polygone de filtrage dans une direction transversale à la portion de courbe section qui passe par les points du sous-ensemble de points de la courbe section. FIG. 7 illustre un cas dans lequel la détermination d’un sommet du polygone de filtrage dans une direction transversale à la portion de courbe section qui passe par les points du sous-ensemble de points de la courbe section comprend une étape d’inversion du segment orienté SO2 pour que ledit segment soit orienté vers l’intérieur de la courbe section. FIG. 8 représente un diagramme de la séquence des étapes du procédé selon un exemple de mise en œuvre de l’invention. Le procédé 100 selon l’invention vise à déterminer une dimension d’une pièce mécanique PM1, PM2, notamment pour vérifier une bonne correspondance entre les dimensions d’une première pièce mécanique PM1 destinée à être rapportée sur une deuxième mécanique PM2, par exemple dans le domaine aéronautique, pour former une aube de turbine. Un exemple d’une pièce mécanique PM1, i.e. un bord d’attaque, destinée à être rapportée sur une deuxième pièce mécanique PM2, i.e. une aube de turbine, pour former une aube de turbine complète, est représentée sur la FIG. 1 . On notera que le procédé 100 peut également s’appliquer à d’autres types de pièces mécaniques, tels que par exemple des aubes de soufflante de turbomachine, ou encore des aubes de compresseur de turbomachine. Pour déterminer une dimension d’une pièce mécanique, le procédé 100 selon l’invention comprend une première étape de mesure 101 d’une position d’une pluralité de points positionnés, dans un plan d’une section de la pièce, sur un bord de la section de la pièce. Le bord de la section de la pièce forme une courbe section CS dans ledit plan de la section, illustrée à titre d’exemple sur la FIG. 2 . Par exemple, la mesure 101 de la position des points de la courbe section CS est mesurée avec un palpeur, ou un capteur radiographique, ledit palpeur ou capteur radiographique pouvant être commandé par un ordinateur, l’ordinateur étant configuré pour contrôler le palpeur ou capteur radiographique, et pour recevoir les mesures réalisées par le palpeur ou capteur radiographique, avant de mettre en œuvre les étapes suivantes du procédé 100. L’étape suivante du procédé 100 est la détermination 102 d’une courbe squelette CSQ comprenant une pluralité de points squelettes. Lesdits points squelettes de la courbe squelette CSQ sont, par définition, chacun situés dans le plan de la courbe section CS, idéalement équidistants d’une première partie B1 de la courbe section CS et d’une deuxième partie B2 de la courbe section CS. La détermination 102 des points squelettes de la courbe squelette CSQ comprend les étapes suivantes : - calcul 1021 d’une triangulation de la pluralité de points de la courbe section, la triangulation comprenant un ensemble de triangles, chaque triangle de l’ensemble de triangles étant formé par trois points de la pluralité de points de la courbe section, - pour chaque triangle de la triangulation, détermination 1022 d’un centre CT dudit triangle, défini par le centre d’un cercle circonscrit audit triangle, et d’un barycentre BT dudit triangle, défini par le barycentre des sommets dudit triangle, de manière à obtenir un ensemble de barycentres BT et un ensemble de centres CT, chaque triangle de la triangulation correspondant à un barycentre BT de l’ensemble de barycentres et à un centre CT de l’ensemble de centres ; la FIG. 4 illustre à titre d’exemple une triangulation obtenue à partir des points de la courbe section CS, les barycentres des triangles de la triangulation représentés par des ronds, et les centres des triangles représentés par des étoiles. - filtrage 1023 des centres CT, sur la base des barycentres BT, pour obtenir un ensemble de centres filtrés CTF. La courbe squelette CSQ ainsi déterminée par le procédé 100 comprend l’ensemble des centres filtrés CTF. Enfin, le procédé 100 comprend la détermination 103 de la dimension D de la pièce mécanique en un point squelette déterminé de la courbe squelette CSQ. Par exemple, la dimension D de la pièce mécanique audit point squelette de la courbe squelette CSQ est égale à une somme d’une première distance D1 du point squelette à la première partie B1 de la courbe section CS, et d’une deuxième distance D2 du point squelette à la deuxième partie B2 de la courbe section CS, la première distance D1 et la deuxième distance D2 étant mesurées transversalement à la courbe squelette. En pratique, les points squelettes, obtenus à partir des centres filtrés CTF, par le procédé 100, ne sont pas parfaitement équidistants de la première partie B1 de la courbe section CS et de la deuxième partie B2 de la courbe section CS. En pratique, par exemple, une différence entre la première distance D1 et la deuxième distance D2 rapportée à une demi somme de la première distance D1 et la deuxième distance D2 est inférieure à 1/100, de préférence inférieure à 1/1000, de préférence égale à 0 dans le cas idéal où le point squelette est parfaitement équidistant de la première partie B1 et de la deuxième partie B2 de la courbe section CS. Plus particulièrement, la première partie B1 et la deuxième partie B2 de la courbe section CS sont respectivement de part et d’autre d’une zone de la courbe section formant une cavité CA de la courbe section. Plus particulièrement encore, la dimension D est mesurée en un point squelette déterminé de la courbe squelette CSQ situé à une distance de la cavité CA de la courbe section CS, la direction transversale à la courbe squelette CSQ étant déterminée, à partir d’un sous-ensemble de points squelettes positionnés autour dudit point squelette déterminé. Par exemple, un nombre de points du sous-ensemble de points squelettes est prédéterminé à cet effet. Selon un exemple de mise en œuvre, l’étape de détermination 102 de la courbe squelette comprend une extrapolation 1024 de points complémentaires PC à partir de l’ensemble de centres filtrés CTF, de sorte que la courbe squelette CSQ comprend les points complémentaires PC et l’ensemble de centres filtrés CTF, comme cela est illustré sur la FIG. 2 . Selon un exemple particulier de mise en œuvre, l’extrapolation 1024 comprend l’utilisation d’un autre sous-ensemble de points squelettes, de la courbe squelette CSQ, pour estimer une tangente à la courbe squelette, ladite tangente définissant une direction d’extrapolation de la courbe squelette CSQ vers une cavité CA de la courbe section CS. Ainsi, par exemple, le nombre de points de l’autre sous-ensemble de points squelettes, utilisés pour estimer la direction d’extrapolation de la courbe squelette, est prédéterminé. A partir de cet autre sous-ensemble de points, une matrice A’ à N’ lignes et 2 colonnes peut être construite, N’ étant ledit nombre de points de l’autre sous-ensemble de points considérés, la première et la deuxième colonne de chaque ligne comprenant respectivement les coordonnées x et y de chaque point de l’autre sous-ensemble de points considérés ; A’ T étant une matrice transposée de A’, M’ = A’ T *A’ est une matrice 2x2. La matrice M peut être diagonalisée et un vecteur propre associé à la valeur propre la plus élevée est déterminé pour définir une direction tangente à une courbe virtuelle qui passe par les points de l’autre sous-ensemble de points ; la courbe squelette peut alors être extrapolée le long de cette direction tangente jusqu’à la cavité CA de la courbe section CS. Selon un exemple plus particulier de mise en œuvre, les points complémentaires sont extrapolés le long d’un morceau de courbe, le morceau de courbe étant par exemple linéaire, ou quadratique. L’étape d’extrapolation 1024 de points complémentaires PC peut en particulier comprendre un ordonnancement 1024bis des points de l’ensemble de centres filtrés CTF. Selon un exemple de mise en œuvre, la triangulation est une triangulation de Delaunay, dans laquelle le cercle circonscrit à un triangle de la triangulation ne contient à l’intérieur de son périmètre aucun sommet d’un autre triangle de la triangulation. La FIG. 3 illustre sur la droite le cas d’une triangulation de Delaunay, alors que la triangulation de gauche ne répond pas à la contrainte selon laquelle le cercle circonscrit à un triangle de la triangulation ne contient à l’intérieur de son périmètre aucun sommet d’un autre triangle de la triangulation. Selon un exemple de mise en œuvre, l’étape de filtrage 1023 comprend les étapes suivantes : - la détermination 10231 qu’un barycentre d’un triangle est un barycentre incohérent, et que le triangle correspondant au barycentre incohérent est triangle incohérent, et que le centre du triangle incohérent est un centre incohérent - la suppression 10232 dudit centre incohérent de l’ensemble des centres, pour obtenir l’ensemble de centres filtrés. Par exemple, l’étape de détermination 10231 qu’un barycentre d’un triangle est un barycentre incohérent comprend : - la construction 102311, à partir des points de la courbe section CS, d’un polygone de filtrage PF; - la détermination 102312 qu’un barycentre est incohérent si il est à l’extérieur dudit polygone de filtrage PF. Ainsi, à titre d’exemple, sur la FIG. 5 , sont représentés les sommets SPF d’un polygone de filtrage PF construit à partir des points de la courbe section CS, et les barycentres et centres des triangles filtrés, i.e. qui restent après l’élimination des triangles dont les barycentres étaient situés à l’extérieur du polygone de filtrage PF. Optionnellement, l’étape de construction 102311 du polygone de filtrage PF comprend la détermination 1023111 d’un ensemble de sommets SPF dudit polygone de filtrage PF ; comme cela est représenté, à titre d’exemple, sur la FIG. 6 , chaque sommet SPF est positionné dans une direction définie DT à partir d’un sous-ensemble de points SEP de la courbe section CS, à une distance prédéterminée dudit sous-ensemble de points SEP de la courbe section CS ; la direction définie DT à partir du sous-ensemble de points SEP de la courbe section CS est transversale à une portion de la courbe section CS qui passe par les points dudit sous-ensemble de points SEP de la courbe section CS. Selon un exemple de mise en œuvre, la direction définie DT à partir du sous-ensemble de points SEP de la courbe section est définie par un vecteur propre, associé à une valeur propre minimale, d’une matrice de covariance calculée à partir des coordonnées des points du sous-ensemble de points SEP de la courbe section CS. Par exemple, la matrice de covariance s’écrit M=A T *A, où A est une matrice à N lignes et 2 colonnes, N étant un nombre de points du sous-ensemble de points SEP considéré, la première et la deuxième colonne de chaque ligne comprenant respectivement les coordonnées x et y de chaque point du sous-ensemble de points SEP considéré, A T est une matrice transposée de A, de sorte que M est une matrice 2x2. Optionnellement, la courbe section CS est fermée de façon à délimiter un intérieur de la courbe section CS ; et le sommet SPF du polygone de filtrage PF est positionné dans la direction DT transversale à la portion de la courbe section CS passant par les points dudit sous-ensemble de points SEP de la courbe section CS, de façon qu’un segment orienté SO, SO1, de ladite portion de la courbe section CS vers ledit sommet SPF, selon ladite direction transversale DT, est orienté vers l’intérieur de la courbe section CS, comme cela est représenté sur la FIG. 6 et sur la FIG. 7 . Il faut noter que, sur l’exemple du segment orienté SO1 illustré sur la FIG. 7 , le segment orienté SO1 est en effet bien orienté vers l’intérieur de la courbe section CS en partant du point origine PO1, placé sur la portion de la courbe section CS qui passe par les points dudit sous-ensemble de points SEP de la courbe section CS, bien que le sommet SPF1 n’est pas à l’intérieur de la courbe section. En revanche le segment orienté SO2, dont le sommet SPF2 est pourtant à l’intérieur de la courbe section CS, n’est pas orienté vers l’intérieur de la courbe section CS, en partant du point origine PO2, placé sur la portion de la courbe section CS qui passe par les points dudit sous-ensemble de points SEP de la courbe section CS. C’est pourquoi la détermination du sommet SPF du polygone de filtrage PF peut en outre comprendre une étape d’inversion de l’orientation dudit segment orienté SO2, illustrée sur la FIG. 7 , lorsque ledit segment orienté SO2 coupe la courbe section CS en un point d’intersection PI2 et qu’un point intermédiaire PI2’ se trouve à l’extérieur de la courbe section, ledit point intermédiaire PI2’ étant positionné entre un point origine PO2 du segment orienté SO2 et ledit point d’intersection PI2, ledit point origine PO2 étant situé sur ladite portion de la courbe section CS. Ladite inversion de l’orientation génère un segment orienté inversé SO2I, de sens opposé au segment orienté SO2, ledit segment orienté inversé SO2I étant dirigé vers l’intérieur de de la courbe section CS. Selon un aspect, l’invention concerne un produit programme d’ordinateur comprenant des instructions de code mettant en œuvre un procédé de détermination de dimension de pièce tel que défini ci-avant, lorsque ces instructions de code sont exécutées par une unité de traitement. Selon un autre aspect, l’invention concerne également un support lisible par ordinateur sur lequel sont enregistrées des instructions de code mettant en œuvre un procédé de détermination de dimension de pièce tel que défini ci-avant Procédé (100) de détermination d’une dimension (D) d’une pièce mécanique (PM1, PM2), le procédé étant mis en œuvre par ordinateur et comprenant les étapes suivantes : - mesure (101) d’une position, dans un plan d’une section de la pièce, d’une pluralité de points positionnés sur un bord de la section de la pièce, le bord de la section de la pièce formant une courbe section (CS) dans ledit plan, - détermination (102) d’une courbe squelette (CSQ) comprenant une pluralité de points squelettes, la détermination (102) de la courbe squelette (CSQ) comprenant les étapes suivantes : - calcul (1021) d’une triangulation de la pluralité de points de la courbe section (CS), la triangulation comprenant un ensemble de triangles, chaque triangle de l’ensemble de triangles étant formé par trois points de la pluralité de points de la courbe section (CS), - pour chaque triangle de la triangulation, détermination (1022) d’un centre (CT) dudit triangle, défini par le centre d’un cercle circonscrit audit triangle, et d’un barycentre (BT) dudit triangle, défini par le barycentre des sommets dudit triangle, de manière à obtenir un ensemble de barycentres (BT) et un ensemble de centres (CT), chaque triangle de la triangulation correspondant à un barycentre (BT) de l’ensemble de barycentres et à un centre (CT) de l’ensemble de centres, - filtrage (1023) des centres (CT), sur la base des barycentres (BT), pour obtenir un ensemble de centres filtrés (CTF), la courbe squelette (CSQ) déterminée comprenant l’ensemble des centres filtrés (CTF), - détermination (103) de la dimension (D) de la pièce mécanique, en un point squelette de la courbe squelette (CSQ), en fonction de la position dudit point squelette sur la courbe squelette (CSQ). Procédé (100) selon la revendication 1, dans lequel ledit point squelette est positionné à une première distance (D1) d’une première partie (B1) de la courbe section (CS), et à une deuxième distance (D2) d’une deuxième partie (B2) de la courbe section (CS), la première distance (D1) et la deuxième distance (D2) étant mesurées transversalement à la courbe squelette, la dimension (D) de la pièce mécanique audit point squelette de la courbe squelette (CSQ) étant égale à une somme de la première distance (D1) et de la deuxième distance (D2). Procédé (100) selon l’une quelconque des revendications 1 ou 2, dans lequel l’étape de filtrage (1023) comprend : - la détermination (10231) qu’un barycentre d’un triangle est un barycentre incohérent, et que le triangle correspondant au barycentre incohérent est un triangle incohérent, et que le centre du triangle incohérent est un centre incohérent, - la suppression (10232) dudit centre incohérent de l’ensemble des centres, pour obtenir l’ensemble de centres filtrés. Procédé (100) selon l’une des revendications 1 à 3, dans lequel l’étape de détermination (102) de la courbe squelette comprend une extrapolation (1024) de points complémentaires à partir de l’ensemble de centres filtrés (CTF), de sorte que la courbe squelette (CSQ) comprend les points complémentaires (PC) et l’ensemble de centres filtrés (CTF). Procédé (100) selon la revendication 4, dans lequel l’étape d’extrapolation (1024) de points complémentaires (PC) comprend un ordonnancement (1024bis) des points de l’ensemble de centres filtrés (CTF). Procédé (100) selon la revendication 3, dans lequel l’étape de détermination (10231) qu’un barycentre d’un triangle est un barycentre incohérent comprend : - la construction (102311), à partir des points de la courbe section (CS), d’un polygone de filtrage (PF); - la détermination (102312) qu’un barycentre est incohérent si il est à l’extérieur dudit polygone de filtrage (PF). Procédé (100) selon la revendication 6, dans lequel l’étape de construction (102311) du polygone de filtrage (PF), comprend la détermination (1023111) d’un ensemble de sommets (SPF) dudit polygone de filtrage (PF), chaque sommet (SPF) étant positionné dans une direction définie (DT) à partir d’un sous-ensemble de points (SEP) de la courbe section (CS), à une distance prédéterminée dudit sous-ensemble de points (SEP) de la courbe section (CS), la direction définie (DT) à partir du sous-ensemble de points (SEP) de la courbe section (CS) étant transversale à une portion de la courbe section (CS) passant par les points dudit sous-ensemble de points (SEP) de la courbe section (CS). Procédé (100) selon la revendication 7, dans lequel la direction définie (DT) à partir du sous-ensemble de points (SEP) de la courbe section est définie par un vecteur propre d’une matrice de covariance calculée à partir des coordonnées des points du sous-ensemble de points (SEP) de la courbe section (CS). Procédé (100) selon l’une des revendications 7 ou 8, dans lequel la courbe section (CS) est fermée de façon à délimiter un intérieur de la courbe section, et dans lequel le sommet (SPF) du polygone de filtrage (PF) est positionné dans la direction (DT) transversale à la portion de la courbe section (CS) passant par les points dudit sous-ensemble de points (SEP) de la courbe section (CS), de façon qu’un segment orienté (SO, SO1, SO2) de ladite portion de la courbe section (CS) vers ledit sommet (SPF), selon ladite direction transversale (DT), est orienté vers l’intérieur de la courbe section (CS). Procédé (100) selon la revendication précédente, comprenant en outre une étape d’inversion de l’orientation dudit segment orienté (SO2), lorsque ledit segment orienté (SO2) coupe la courbe section (CS) en un point d’intersection (PI2) et lorsqu’un point intermédiaire (PI2’) se trouve à l’extérieur de la courbe section, ledit point intermédiaire (PI2’) étant positionné entre un point origine (PO2) du segment orienté (SO2) et ledit point d’intersection (PI2), ledit point origine (PO2) étant situé sur ladite portion de la courbe section (CS), ladite inversion de l’orientation générant un segment orienté inversé (SO2I), de sens opposé au segment orienté (SO2), ledit segment orienté inversé (SO2I) étant dirigé vers l’intérieur de de la courbe section (CS).