Brevet FR3156962A1 — publié le 20250620

Procédé (10) de modélisation fonctionnelle d'un élément de relief (3) présent dans une surface, comprenant les étapes suivantes : - acquisition (11) d'une image tridimensionnelle de ladite surface, cette image se présentant sous la forme d'un nuage de points décrivant cette surface ; - sélection (12) d'une pluralité de points du nuage de points marquant un élément de relief cible joignant lesdits points sélectionnés, cet élément de relief cible s'étendant entre un premier point sélectionné et un deuxième point sélectionné ; - extraction (13) du nuage de points d'un sous-nuage de points autour de l'élément de relief cible, ce sous-nuage de points décrivant une sous-surface ; - modélisation (14) de la sous-surface au moyen d'une surface de forme libre ; - détermination (15), à partir d'un modèle de la sous-surface, d'une première ligne topographique de la sous-surface ; - détermination (16) d'une deuxième ligne topographique et d'une troisième ligne topographique de part et d'autre de la première ligne. Figure à publier avec l’abrégé : Fig 2 Procédé de modélisation fonctionnelle d'un élément de relief présent dans une surface tridimensionnelle, dispositif électronique mettant en œuvre ledit procédé La présente invention a trait aux procédés de modélisation fonctionnelle d'un élément de relief compris dans une surface tridimensionnelle, ainsi qu'aux produits programme d'ordinateur pour adapter des dispositifs électroniques de sorte que ces derniers mettent en œuvre de tels procédés. On entend par « élément de relief » un sillon ou une saillie linéaire présente sur une surface tridimensionnelle. Un sillon désigne toute dépression linéaire, enfoncement linéaire, pli ou rainure présents à la surface d'un objet tels qu'un sillon cutané, un sillon nasogénien, un pli, une ride, une ridule, ou une entaille dans une peau, une fissure ou une fente dans un matériau ou une tranchée ou une faille dans un sol. Quant au terme « saillie linéaire », il fait référence à toute élévation, protubérance ou, plus généralement, à toute déformation locale où une partie de la surface s'étend au-dessus du niveau environnant le long d'une ligne telle qu'une ligne de suture, une varice ou une cicatrice dans une peau, une ligne de soudure ou des nervures sur la surface d'une pièce manufacturée ou une vague à la surface de l'eau. Un élément de relief peut être naturel, intentionnellement créé pour diverses raisons fonctionnelles ou esthétiques ou résultant de contraintes, d'usure ou d'autres facteurs. Il est connu de l'état de la technique des méthodes pour la détection et la modélisation d'éléments de relief présents sur la surface d'un objet réel. Cependant, ces méthodes se limitent à une représentation numérique de l'apparence dans l'espace, c'est-à-dire à des caractéristiques géométriques en trois dimensions (3D), des éléments de relief. Une telle représentation statique ne permet pas d’étudier convenablement le comportement dynamique d’un élément de relief. Un objet de la présente invention est de proposer un ou plusieurs modèles fonctionnels d'un élément de relief, autrement dit des modèles qui vont au-delà d'une simple représentation géométrique 3D statique dans le but de simuler le comportement, le fonctionnement et/ou les interactions de cet élément de relief dans un environnement virtuel. Un autre objet de la présente invention est de produire un modèle tridimensionnel fonctionnel, à savoir un modèle géométrique dont l'utilisation ne se limite pas à la visualisation mais permet aussi la mesure, la simulation mécanique et l'évaluation dynamique d'un élément de relief. Il est proposé, en premier lieu, un procédé de modélisation fonctionnelle d'un élément de relief présent dans une surface, ce procédé comprenant les étapes suivantes : acquisition d'une image tridimensionnelle de ladite surface, cette image se présentant sous la forme d'un nuage de points décrivant cette surface ; sélection d'une pluralité de points du nuage de points marquant un élément de relief cible joignant lesdits points sélectionnés, cet élément de relief cible s'étendant entre un premier point sélectionné et un deuxième point sélectionné de ladite pluralité de points sélectionnés ; extraction du nuage de points d'un sous-nuage de points autour de l'élément de relief cible, ce sous-nuage de points décrivant une sous-surface de ladite surface ; modélisation de la sous-surface au moyen d'une surface de forme libre paramétrée à l'aide de points de contrôle ; détermination, à partir d'un modèle de la sous-surface, d'une première ligne topographique de la sous-surface joignant le premier point sélectionné et le deuxième point sélectionné ; détermination, à partir du modèle de la sous-surface, d'une deuxième ligne topographique et d'une troisième ligne topographique de part et d'autre de la première ligne topographique. Diverses caractéristiques supplémentaires peuvent être prévues, seules ou en combinaison : l'élément de relief est un sillon ; la première ligne topographique est une ligne de fond du sillon ; l'élément de relief est une saillie linéaire ; la première ligne topographique est une ligne de crête de la saillie linéaire ; la surface déformable est une surface B-Spline cubique uniforme ; la forme initiale de la surface de forme libre est choisie en fonction de données métier ; la surface est initialisée par un treillis de points de contrôle compris dans un plan intégrant le premier point sélectionné et le deuxième point sélectionné ; l’étape de détermination de la deuxième ligne topographique comprend une recherche d’un changement de pente, par rapport à un axe perpendiculaire à la direction longitudinale de l'élément de relief cible, entre un premier segment sensiblement rectiligne et un deuxième segment sensiblement rectiligne successifs de la sous-surface en partant de la première ligne topographique. Il est proposé, en deuxième lieu, un produit programme d’ordinateur implémenté sur un support de mémorisation et susceptible d’être mis en œuvre au sein d’une unité de traitement d’un dispositif électronique et comprenant des instructions pour la mise en œuvre du procédé présenté ci-dessus. L’invention concerne en outre un support de mémorisation lisible par un ordinateur comportant les instructions d’un tel produit programme d’ordinateur. L’invention concerne un dispositif électronique comprenant : une unité de traitement ; une mémoire de données ; une mémoire de programmes comprenant les instructions de programme d’un tel produit programme d’ordinateur. D'autres caractéristiques et avantages de l'invention apparaîtront plus clairement et de manière concrète à la lecture de la description ci-après de modes de réalisation, laquelle est faite en référence aux dessins annexés dans lesquels : FIG. 1 illustre schématiquement une surface comprenant un élément de relief selon divers modes de réalisation ; FIG. 2 illustre schématiquement des étapes d’un procédé de modélisation fonctionnelle d'un élément de relief présent dans une surface selon divers modes de réalisation ; FIG. 3 illustre schématiquement un sous-nuage de points extrait d’un nuage de points décrivant ladite surface ; FIG. 4 illustre schématiquement une surface de forme libre selon divers modes de réalisation ; FIG. 5 illustre schématiquement des lignes topographiques caractérisant un élément de relief selon divers modes de réalisation ; FIG. 6 illustre un dispositif électronique adapté pour mettre en œuvre un procédé selon l’invention. La FIG. 1 illustre un exemple illustratif d’une image 1 en trois dimensions ou 3D. Une telle image est celle d'une zone corporelle, en l’espèce un visage humain, comprenant un ou plusieurs éléments de relief 3 cutanés tels qu'une ride. En variante, une telle image pourrait décrire une ridule, une cicatrice, une varice, une ligne de suture d’une toute autre partie d’un corps humain. Une telle image 1 pourrait être celle d'une coque comprenant une fissure, une ligne de soudure ou une ligne de collage en lien avec un objet physique. En lien avec la FIG. 1 , la FIG. 2 présente un procédé 10 de modélisation fonctionnelle d'un élément de relief présent dans une surface 3D conforme à l’invention, ledit procédé 10 étant destiné à être mis en œuvre par une unité de traitement d’un dispositif électronique tel qu’un téléphone mobile intelligent, un ordinateur portable ou fixe. Un tel procédé 10 comporte une première étape 11 d’acquisition d'une image tridimensionnelle de ladite surface à l’aide de moyens de capture dudit dispositif électronique ou coopérant avec ce dernier. Cette étape 11 consiste in fine , à produire une image se présentant sous la forme d'un nuage 2 de points décrivant cette surface capturée. La FIG. 1 illustre un tel nuage 2 de points assemblé en un volume de points qui représente la topologie de la surface à étudier. Chaque point de ce nuage 2 de points est représenté par des coordonnées ou par une donnée de distance et un angle solide dans un repère prédéfini (absolu ou relatif) à trois dimensions. Ce repère est le plus souvent associé au dispositif de prise de vue de l'image 3D 1 tel qu'une caméra 3D, une caméra de profondeur, un dispositif de photogrammétrie ou un scanner 3D (scanner laser, Lidar, scanner stéréoscopique, scanner à lumière modulée ou scanner à silhouette par exemple). Bien entendu, des attributs ou des données associées à chaque point du nuage 2 de points peuvent, en outre, comprendre d’autres informations directement issues de l’image capturée, telles que des informations d’intensité lumineuse, de couleur ou de niveaux de gris. L’invention prévoit en outre que ladite étape 11 comporte un traitement consistant à analyser une telle image capturée pour en déduire des informations d’ordre sémantique ou « connaissances métier » et ainsi enrichir les attributs ou données associés aux points du nuage 2. Une telle analyse permet ainsi d’identifier le corps physique ou une zone dudit corps sur lequel ou laquelle porte la capture d’image (que ledit corps soit un corps humain, animal, un végétal ou minéral, voire un objet matériel, pour déterminer des propriétés physiques (par exemple, des propriétés physiques propres à un tissu, un organe, un os ou un muscle pour un matériau ou corps biologique), et/ou des propriétés biomécaniques (par exemple, la rigidité, l'élasticité, la ténacité, la viscoélasticité ou la résistance à la compression, etc.) et ainsi produire des attributs complémentaires associés aux points du nuage 2 de points. Un tel traitement peut s’appuyer sur des techniques de reconnaissance de formes ou d’identification d’objets mettant en œuvre avantageusement mais non limitativement des techniques de « Machine Learning » ou plus largement toute techniques d'apprentissage automatique se fondant sur des approches mathématiques et statistiques pour donner à un ordinateur la capacité d'« apprendre » à partir de données, sans être explicitement programmés à dessein. Un tel traitement peut s’appuyer également sur un ou plusieurs modèles de surfaces (ou peaux) capturables, modèles traduisant de telles propriétés physiques et biomécaniques selon la nature desdites surfaces, pour produire lesdits attributs complémentaires évoqués ci-dessus ou encore sur des entrées ou saisies de l’utilisateur du dispositif électronique mettant en œuvre le procédé 10. Ce dernier comporte une étape de sélection 12, selon laquelle une pluralité de points 21-22 du nuage 2 de points est sélectionnée pour marquer un élément de relief 30 cible joignant lesdits points 21-22 sélectionnés. Cet élément de relief 30 cible s'étend entre un premier point 21 sélectionné et un deuxième point 22 sélectionné de ladite pluralité de points 21-22 sélectionnés. Ces points 21-22 sont, par exemple, sélectionnés dans une vue 3D à l'aide d'un moyen de sélection de type souris d'ordinateur ou équivalent. Il s'agit d'un marquage initial de l'élément de relief 3 à modéliser pour désigner sa position approximative sur la surface décrite par le nuage 2 de points. Dans un mode de réalisation, l'élément de relief 30 cible est marqué seulement par un premier et un deuxième point 21-22 sélectionnés représentant les deux extrémités de cet élément de relief 30 cible ayant la forme d'un segment. En alternative, une pluralité de points 21-22 sélectionnés intermédiaires (notamment, deux ou trois) peut être envisagée pour marquer un élément de relief 30 cible courbe. La sélection des points 21-22 par un opérateur (par exemple, un praticien médical ou un technicien) permet, avantageusement, de récupérer de la connaissance métier. Dans ce cas, les points 21-22 sélectionnés marquant les deux extrémités de l'élément de relief 30 cible sont considérés comme étant précisément acquises. Dans un mode de réalisation, les « connaissance métier » comprennent des caractéristiques de la structure anatomique étudiée. En liaison avec la FIG. 3 , comme l’illustre ladite FIG. 3 , un sous-nuage 20 de points autour de l'élément de relief 30 cible est extrait du nuage 2 de points lors d'une étape d'extraction 13, ce sous-nuage 20 de points décrivant une sous-surface de la surface initiale. Le sous-nuage 20 de points intègre, par exemple, les points du nuage 2 de points situés à moins d'une distance prédéfinie autour de l’élément de relief 30 cible. Le sous-nuage 20 de points est donc associé à une sous-surface locale au sillon ou à la saillie linéaire. Comme l’indique ladite FIG. 3 , un repère cartésien à trois dimensions R(Ox, Oy, Oz) est associé au sous-nuage 20 extrait de sorte que les deux points 21, 22 d'extrémité de l'élément de relief 30 cible soient sur l'axe Oy à égale distance du centre O du repère. Autrement dit, les deux extrémités de l'élément de relief 30 cible sont centrées sur l'origine du repère cartésien sur l'axe Oy (c’est-à-dire situées sur l'axe Oy symétriquement de part et d’autre de l’origine O du repère R). L'origine O du repère est mise au milieu du segment joignant les deux extrémités de l'élément de relief 30 cible, tel qu’il est illustré par la FIG. 4 . Ces deux extrémités ont pour coordonnées (0, -y0, 0) et (0, y0, 0) où -y0 et y0 sont respectivement l'ordonnée du premier et du deuxième points 21, 22 sélectionnés délimitant l'élément de relief 30 cible. En outre, le plan xOy est choisi de sorte que la projection sur ce plan du sous-nuage 20 de points soit maximale. En effet, une rotation autour de l'axe Oy est appliquée au sous-nuage 20 de points ou, de manière équivalente au repère R(Ox, Oy, Oz), jusqu'à ce que la projection sur le plan xOy de la sous-surface décrite par le sous-nuage 20 de points soit maximale. Cette sous-surface est donc au plus près du plan xOy et sa forme est principalement décrite par des élévations et/ou des dépressions suivant l'axe Oz. Une telle disposition du repère R(Ox, Oy, Oz) met, avantageusement, « à plat » sur le plan xOy le sous-nuage 20 de points d'intérêt, l'axe Oz donnant la hauteur/profondeur de l'élément de relief 30 cible. Cette étape vise, en particulier, à focaliser les étapes suivantes sur la sous-surface d'intérêt autour de l'élément de relief 30 cible à modéliser. S'agissant d'isométries appliquées au sous-nuage 20 de points, il est possible de revenir dans le repère initial associé au nuage 2 de points avec des transformations inverses, les calculs de mesure étant identiques dans les deux repères. En lien avec les figures 2 et 4, lors d'une étape de modélisation 14 d’un procédé 10 selon l’invention, la sous-surface décrite par le sous-nuage 20 de points est modélisée au moyen d'une surface 23 de forme libre (continue et dérivable) de type surface B-Spline, une surface NURBS (« Non Uniform Rational Basis Splines » selon une terminologie anglo-saxonne), une surface de Bézier, ou une surface T-Spline. Comme l’indique la FIG. 4 , la surface 23 de forme libre est paramétrée à l'aide de points 24 de contrôle. L'étape de modélisation 14 ajuste de façon itérative la surface 23 de forme libre jusqu'à ce qu'elle représente le sous-nuage 20 de points en passant à travers celui-ci de manière suffisamment proche pour en gommer les artefacts tout en mettant en évidence ses caractéristiques de forme. L'ajustement de la surface 23 de forme libre s'effectue par déplacement de ses points 24 de contrôle. En effet, en se reportant à la FIG. 4 , un algorithme itératif permet d'affiner la mise en correspondance des points 24 de contrôle de la surface 23 de forme libre avec la sous-surface décrite par le sous-nuage 20 de points, en minimisant la distance qui les sépare. Cet algorithme met en œuvre un procédé de recalage itératif pour mettre en coïncidence la surface 23 de forme libre avec la sous-surface décrite par le sous-nuage 20 de points, en déplaçant un ou plusieurs points 24 de contrôle de la surface 23 de forme libre. Dans un mode de réalisation, la surface 23 de forme libre est une surface B-Spline cubique uniforme. Avantageusement, une surface B-Spline cubique uniforme permet de reproduire les variations locales de forme tout en maintenant un niveau de continuité et de dérivabilité en accord avec la surface à modéliser. Dans un mode de réalisation, une surface B-Spline cubique uniforme est associée à un treillis local de seize points 24 de contrôle (4x4), ce treillis local étant une partie d'un treillis global de NxM points 24 de contrôle qui couvre la sous-surface décrite par le sous-nuage 20 de points à modéliser. Le déplacement d’un point 24 de contrôle produit une déformation de la surface 23 de forme libre sur les quatre carreaux voisins d’ordre un et les douze carreaux voisins d’ordre deux (les déformations sur les douze voisins d’ordre deux étant assez faibles par rapport à celles appliquées sur les quatre carreaux voisins d’ordre un). Dans un mode de réalisation, une surface 23 de forme libre est, d'abord, initialisée par un treillis de NxM points 24 de contrôle disposés de manière régulière dans le plan xOy (i.e. leur coordonnée z vaut 0) et comprenant le segment joignant les deux extrémités sur l'axe Oy de l'élément de relief 30 cible. Plus généralement, la surface 23 de forme libre est initialisée par un treillis de points 24 de contrôle compris dans un plan intégrant le premier point 21 sélectionné et le deuxième point 22 sélectionné. Ensuite, pour chaque point 24 de contrôle, un déplacement D selon la direction Oz maximisant un critère de proximité prédéfini des quatre carreaux voisins aux points du sous-nuage 20 de points qu’ils recouvrent est déterminé. Le critère de proximité est, par exemple, la somme des carrés des distances (selon la direction Oz) des points du sous-nuage 20 à la surface 23 de forme libre. Le point 24 de contrôle dont le déplacement D serait le plus grand (c.à.d. le plus distant de la sous-surface décrite par le sous-nuage de points) est déplacé de la distance D. Ces deux dernières étapes peuvent être réitérées tant que le critère de proximité prédéfini est supérieur à une valeur seuil prédéfinie. Il en résulte, avantageusement, un modèle surfacique, par exemple de type B-Spline, ajusté à la sous-surface décrite par le sous-nuage 20 de points. Ce modèle est défini par une ou plusieurs fonctions paramétriques de la surface libre, en fonction de la position finale des points 24 de contrôle. Les points 24 de contrôle forment des carreaux élémentaires de surface (ou «patchs ») permettant de mettre en évidence et d’éliminer le bruit et les artéfacts (points aberrants dans le sous-nuage 20 de points), de représenter localement le sous-nuage 20 de points par une surface lisse et avec une précision maîtrisée, et de recaler des données et/ou intégrer des nouvelles données. Dans un autre mode de réalisation, la forme initiale de la surface 23 de forme libre est choisie en fonction de données métier, notamment la structure géométrique ou anatomique du sous-nuage 20 de points (par exemple, les relations entre les phalanges d’une main ou celles entre les éléments d'un visage, les contraintes physiologiques sur ces composantes, les propriétés mécaniques de l’ensemble des muscles et des tendons, ainsi que de la peau, l’existence de sillons ou de veines apparentes), permettant d'accélérer la convergence des étapes ci-dessus pour la modélisation du sous-nuage 20 de points et améliorer sa précision. Une surface 23 de forme libre dont la géométrie initiale est similaire à celle du sous-nuage 20 de points cible permet, avantageusement, de réduire le nombre d'itérations des étapes de modélisation précitées. Avantageusement, des données métier indiquant des points représentatifs de la sous-surface à modéliser permettent, en outre, de réduire un éventuel bruit compris dans les données de position des points du sous-nuage 20 de points. En liaison avec la FIG. 2 et comme illustré par la FIG. 5 , à partir du modèle obtenu de la sous-surface, le procédé 10 comporte une étape 15 de détermination d’une première ligne topographique 31 de ladite sous-surface joignant le premier point 21 sélectionné et le deuxième point 22 sélectionné. Cette première ligne topographique 31 correspond, en effet, à une ligne de fond (dans le cas d'un sillon) ou une ligne de crête (dans le cas d'une saillie linéaire) de l'élément de relief 30 cible. La première ligne topographique est contrainte au niveau des deux extrémités sélectionnées. Le premier point 21 sélectionné et le deuxième point 21 sélectionné constituent deux points de passage prédéfinis de la première ligne topographique 31. La première ligne topographique 31 joignant les points 21, 22 sélectionnés peut être visualisée comme une ligne brisée (Spline de degré un) ou comme une ligne « souple » (Spline de degré trois). Cette première ligne topographique 31 constitue une première ligne d’accrochage d'un modèle décrivant/reproduisant la ligne de fond ou la ligne de crête de l'élément de relief 30 cible. Partant d'une ligne brisée initiale joignant les points 21, 22 sélectionnés au préalable, ladite étape 15 consiste à échantillonner cette ligne brisée uniformément selon l'axe Oy, lorsque le nombre des points 21-22 sélectionnés est inférieur à une valeur prédéfinie (par exemple, 10). Autrement dit, en disposant d'une ligne brisée définie par les points 21-22 sélectionnés approximant la première ligne topographique 31, une nouvelle ligne brisée comportant plus de points ou, plus généralement, un nombre prédéfini de points en fonction de l’application, est déterminée. Cet échantillonnage permet, avantageusement, d'avoir un nombre suffisant de points intermédiaires entre les deux extrémités de l'élément de relief 30 cible. Dans un mode de réalisation, la première ligne topographique 31 est déterminée à l'aide d'un modèle élastique contraint; soit une ligne brisée joignant deux point fixes (les deux extrémités de l'élément de relief 30) et pouvant se déformer jusqu’à épouser la ligne de fond d'un sillon ou la ligne de crête d'une saillie linéaire dans la mesure d’une contrainte maximale d’élasticité. Pour chaque point intermédiaire (qu'il soit sélectionné ou ajouté par échantillonnage, à l'exception des extrémités), un extremum local sur la sous-surface autour de ce point intermédiaire est recherché, à savoir le point le plus bas (dans le cas d'un sillon) ou le point le plus haut (dans le cas d'une saillie linéaire) dans un voisinage immédiat dudit point intermédiaire (par exemple, à moins d'une distance prédéfinie dudit point intermédiaire). Cet extremum local remplace le point intermédiaire correspondant. Cette recherche d'extremums locaux est itérée tant que la longueur de la première ligne topographique 31 est inférieure à une proportion prédéfinie de la longueur initiale de la première ligne topographique initiale ou tant que l'écart entre deux longueurs successives de la première ligne topographique 31 est inférieur à un seuil prédéfini (par exemple, 2% ou 5%). La ligne ainsi déterminée est utilisée comme support pour une courbe Spline cubique qui caractérise la première ligne topographique 31 de la sous-surface décrite par le sous-nuage 20 de points. Une deuxième ligne topographique 32 et une troisième ligne topographique 33 de part et d'autre de la première ligne topographique 31 sont également déterminées, en une étape 16 du procédé 10, à partir du modèle de la sous-surface décrite par le sous-nuage 20 de points. La deuxième ligne topographique 32 et la troisième ligne topographique 33 correspondent aux bords de l'élément de relief 30 cible de part et d'autre de la première ligne topographique 31. Les deux lignes de bord sont, dans un mode de réalisation, déterminées en faisant évoluer un plan parallèle au plan xOy dans la direction de l'axe Oz. Pour cela, l’étape 16 consiste à faire évoluer verticalement un plan horizontal (c’est-à-dire, pour lequel z est constant) et à conserver la ligne d'intersection de ce plan avec la sous-surface qui est la plus proche ou, plus généralement, qui est située à moins d'une distance prédéfinie de la première ligne topographique 31. Des connaissances métier peuvent, avantageusement, être utilisées pour définir ladite distance prédéfinie en fonction de l'application envisagée (par exemple, la détection de rides/ridules, de fissures/fentes ou la détection de varices/varices réticulaires.) Dans un autre mode de réalisation, les deux lignes de bord de l'élément de relief sont détectées en cherchant « un point de cassure » sur la sous-surface dans une direction perpendiculaire à la direction longitudinale de l'élément de relief. Ce point de cassure marque un changement de pente par rapport à l'axe Oz (ou tout autre axe perpendiculaire à la direction longitudinale de l'élément de relief 30 cible) entre un premier segment sensiblement rectiligne et un deuxième segment sensiblement rectiligne successifs (le premier segment et le deuxième segment étant joints en ce point de cassure) de la sous-surface en partant de la première ligne topographique 31. Le point de cassure est détectable approximativement sur « l’arrondi » du passage de l’intérieur à l’extérieur de l'élément de relief. En prenant un premier point A du fond/crête de l'élément de relief et un deuxième point B à l'extérieur de l'élément de relief dans la direction de l'axe Ox, le point de cassure C est le point joignant deux sections sensiblement rectilignes AC et CB. Par exemple, lorsque la première ligne topographique 31 est composée d'une pluralité de points (par exemple 10 points) intermédiaires entre le premier point 21 sélectionné et le deuxième point 22 sélectionné, l’étape 16 consiste alors à mettre en œuvre le processus suivant pour chacun desdits points intermédiaires : échantillonnage de l'intersection de la sous-surface avec le plan xOy passant par ledit point intermédiaire ; recherche du premier changement de direction (par exemple, à l'aide d'un produit vectoriel des points d’échantillonnage successifs) supérieur à un seuil prédéfini dans les deux sections de part et d'autre dudit point intermédiaire pour détecter, respectivement, un premier point de cassure compris dans la première ligne de bord et un deuxième point de cassure compris dans la deuxième ligne de bord. Avantageusement, l’un quelconque des modes de réalisation décrits ci-dessus pour la détection des lignes de bord peut être adopté suivant la morphologie de l'élément de relief et l’utilisation visée de la modélisation fonctionnelle de cet élément de relief (par exemple, l'application médicale envisagée de la modélisation fonctionnelle d'une ride). Le procédé 10 permet ainsi de disposer d'un modèle de forme libre contraint au niveau des deux extrémités sélectionnées de l'élément de relief 30 cible et de trois lignes (ligne de fond ou de crête et les deux lignes de bord). Une étape d'optimisation supplémentaire du modèle obtenu peut être envisagée en reprenant l'étape d’optimisation précédente, mais en considérant les lignes topographiques 31-33 déterminées comme des contraintes à respecter. Une telle étape permet d'améliorer la fiabilité et la robustesse du modèle tridimensionnel. Il en résulte avantageusement un modèle permettant toutes les simulations et mesures sur une zone d'intérêt de la surface étudiée. Il est, par exemple, possible de récupérer divers descripteurs de l'élément de relief tels que sa longueur, sa courbure, la variation de sa largeur, de sa profondeur ou de sa hauteur sur sa longueur, son volume, ses déformations éventuelles tenant compte des propriétés mécaniques (ou biomécaniques) de la surface, ses éventuelles réactions à des actions prédéfinies (injection de produit, tensions, pressions par exemple). Lorsqu'il s'agit d'une surface cutanée (par exemple, la surface d'un visage), les modes de réalisation décrits ci-dessus permettent, avantageusement, une caractérisation exhaustive du relief cutané humain (par exemple, une caractérisation des rides ou d’une cicatrice selon de nombreux critères géométriques tels que par exemple la longueur, la largeur, la profondeur, le volume, ou la courbure) de sorte à pouvoir rechercher plus facilement un traitement adapté, de suivre dans le temps ou d'évaluer le succès d'une intervention ou l'efficacité d'un traitement (un produit de soin par exemple). Cette modélisation tridimensionnelle fonctionnelle inclut des informations sur la manière dont l'élément de relief se comporte ou interagit avec son environnement ou en réponse à une action prédéfinie. Pour cela, le modèle tridimensionnel fonctionnel obtenu intègre, avantageusement, des paramètres mécaniques en rapport avec l’étude biomécanique de la structure anatomique étudiée. Avantageusement, les modes de réalisation présentés ci-dessus, dont notamment ceux de la modélisation de la sous-surface au moyen d'une surface de forme libre, permettent d'intégrer des données de modalités différentes, de résolutions différentes, et non recalées géométriquement dans l’espace, d'éliminer le bruit et les artefacts et de maîtriser la précision avec laquelle les résultats sont produits. Comme l’indique la FIG. 6 , un dispositif électronique 100 agencé pour mettre en œuvre un procédé conforme à l’invention (tel que décrit en lien avec la FIG. 2 ) comporte une unité centrale de traitement 111 (sous la forme d’un ou plusieurs microprocesseurs ou microcontrôleurs) commandant, par des signaux acheminés par un bus de communication symbolisé sur la FIG. 6 par des doubles flèches en traits simples, des éléments électroniques dont une mémoire. Cette dernière comprend une mémoire de données 112 et une mémoire de programmes 113, lesdites mémoires 112 et 113 pouvant ne former qu’une seule et même entité physique. On entend par « mémoire » toute mémoire informatique qu’elle soit volatile ou non. Une mémoire non volatile est une mémoire informatique dont la technologie conserve ses données en l’absence d’une alimentation en énergie électrique. Elle peut contenir des données résultant de saisies, de calculs, de mesures et/ou des instructions de programmes. Les principales mémoires non volatiles actuellement disponibles sont inscriptibles électriquement telles que la technologie EPROM (« Erasable Programmable Read-Only Memory », selon une terminologie anglo-saxonne) ou encore inscriptibles et effaçables électriquement telles que les technologies EEPROM (« Electrically-Erasable Programmable Read-Only Memory », selon une terminologie anglo-saxonne), flash, SSD (« Solid-State Drive », selon une terminologie anglo-saxonne), etc. Les mémoires non volatiles se distinguent des mémoires dites « volatiles » dont les données sont perdues en l’absence d’une alimentation électrique. Les principales mémoires volatiles actuellement disponibles exploitent les technologies RAM (« Random Access Memory » selon une terminologie anglo-saxonne ou encore nommée « mémoire vive »), DRAM (mémoire vive dynamique, nécessitant une réactualisation régulière), SRAM (mémoire vive statique nécessitant une telle réactualisation lors d’une sous-alimentation électrique), DPRAM ou VRAM (particulièrement adaptées à la vidéo), etc. Un dispositif électronique 100 peut comporter en outre un ou plusieurs capteurs d’images matricielles 114 sous la forme d’une caméra de profondeur (ou une caméra 3D). Un tel dispositif électronique 100 peut comporter des moyens de communication 115 avec le monde extérieur sous la forme d’une unité d'entrée et une unité de sortie. Pour fonctionner, un tel dispositif électronique 100 comporte généralement une source d’énergie électrique 116, externe ou interne sous la forme d’une ou plusieurs batteries par exemple. L’unité de traitement 111 peut également comporter des moyens de pilotage 117 d'une interface homme-machine d’entrée et/ou de sortie. On entend par « interface homme-machine de sortie », tout moyen, employé seul ou en combinaison, permettant de sortir ou délivrer une représentation graphique, haptique, sonore ou, plus généralement, perceptible par l'humain. Une telle interface homme-machine de sortie peut consister de manière non exhaustive en un ou plusieurs écrans, haut-parleurs ou autres moyens alternatifs adaptés. On entend par « interface homme-machine d’entrée », un clavier informatique, un dispositif de pointage, un écran tactile, un microphone ou, plus généralement, toute interface agencée pour traduire une gestuelle ou une consigne émise par un humain en données de commande ou de paramétrage. Avantageusement, les interfaces homme-machine d'entrée et de sortie peuvent ne constituer qu'une seule et même entité physique. Un tel dispositif électronique 100 peut consister en un téléphone portable intelligent (ou « smartphone » selon une terminologie anglo-saxonne), une tablette électronique, un ordinateur de bureau ou portable. L’agencement du procédé de modélisation fonctionnelle conforme à l’invention, de nature à pouvoir être mis en œuvre par un tel dispositif 100 que nous pourrions qualifier de « grand public », permet l’exploitation de moyens de capture non spécifiques ou spécialisés pour la mise en œuvre de l’invention. Un tel dispositif 100, par exemple un smartphone, peut être adapté en chargeant dans sa mémoire de programmes 113 un programme d'ordinateur P comportant des instructions pour provoquer, lors de leurs exécutions, la mise en œuvre d’un procédé conforme à l’invention. Ce programme d'ordinateur P peut utiliser n'importe quel langage de programmation, et se présenter sous la forme de codes source, codes objet, ou de codes intermédiaires entre codes source et codes objet, tel que dans une forme interprétée, partiellement ou entièrement compilée, ou sous n'importe quelle autre forme souhaitable. Procédé (10) de modélisation fonctionnelle d'un élément de relief (3) présent dans une surface (1), ce procédé (10) étant destiné à être mis en œuvre par une unité de traitement (111) d’un dispositif électronique (100) et comprenant les étapes suivantes : acquisition (11) d'une image tridimensionnelle (1) de ladite surface, cette image se présentant sous la forme d'un nuage (2) de points décrivant cette surface ; sélection (12) d'une pluralité de points (21-22) du nuage de points marquant un élément de relief (30) cible joignant lesdits points (21-22) sélectionnés, cet élément de relief (30) cible s'étendant entre un premier point (21) sélectionné et un deuxième point (22) sélectionné de ladite pluralité de points sélectionnés ; extraction (13) du nuage (2) de points d'un sous-nuage (20) de points (21-22) autour de l'élément de relief (30) cible, ce sous-nuage de points décrivant une sous-surface de ladite surface ; modélisation (14) de la sous-surface au moyen d'une surface (23) de forme libre paramétrée à l'aide de points (24) de contrôle ; détermination (15), à partir d'un modèle de la sous-surface, d'une première ligne topographique (31) de la sous-surface joignant le premier point (21) sélectionné et le deuxième point (22) sélectionné ; détermination (16), à partir du modèle de la sous-surface, d'une deuxième ligne topographique (32) et d'une troisième ligne topographique (33) de part et d'autre de la première ligne topographique (31). Procédé (10) selon la revendication précédente, caractérisé en ce que l'élément de relief (3) est un sillon. Procédé (10) selon la revendication précédente, caractérisé en ce que la première ligne topographique (31) est une ligne de fond du sillon. Procédé (10) selon la revendication 1, caractérisé en ce que l'élément de relief (3) est une saillie linéaire. Procédé (10) selon la revendication précédente, caractérisé en ce que la première ligne topographique (31) est une ligne de crête de la saillie linéaire. Procédé (10) selon l’une quelconque des revendications précédentes, caractérisé en ce que la surface (23) de forme libre est une surface B-Spline cubique uniforme. Procédé (10) selon l’une quelconque des revendications précédentes, caractérisé en ce que la forme initiale de la surface (23) de forme libre est choisie en fonction de données métier. Procédé (10) selon l’une quelconque des revendications 1 à 6, caractérisé en ce que la surface (23) est initialisée par un treillis de points (24) de contrôle compris dans un plan intégrant le premier point (21) sélectionné et le deuxième point (22) sélectionné. Procédé (10) selon l’une quelconque des revendications précédentes, caractérisé en ce que l’étape de détermination (16) de la deuxième ligne topographique (32) comprend une recherche d’un changement de pente, par rapport à un axe perpendiculaire à la direction longitudinale de l'élément de relief (30) cible, entre un premier segment sensiblement rectiligne et un deuxième segment sensiblement rectiligne successifs de la sous-surface en partant de la première ligne topographique (31). Dispositif électronique (100) comprenant : une unité traitement (111) ; une mémoire de données (112) ; une mémoire de programmes (113) ; ledit dispositif électronique étant caractérisé en ce que l’unité de traitement (111) est agencée pour mettre en œuvre un procédé (10) de modélisation fonctionnelle d'un élément de relief (3) présent dans une surface (1) selon l’une quelconque des revendications précédentes. Produit programme d’ordinateur (P) implémenté sur un support mémoire, susceptible d’être mis en œuvre au sein d’une unité de traitement informatique et comprenant des instructions pour la mise en œuvre d'un procédé de modélisation fonctionnelle d'un élément de relief présent dans une surface selon l'une quelconque des revendications 1 à 9. Support de mémorisation lisible par un ordinateur comportant les instructions d’un produit programme d’ordinateur selon la précédente revendication.